Skip to content

Commit d325b69

Browse files
alibuildalibuild
committed
Please consider the following formatting changes
1 parent 25abc2a commit d325b69

1 file changed

Lines changed: 126 additions & 121 deletions

File tree

‎PWGCF/EbyEFluctuations/Tasks/netchargeFluctuations.cxx‎

Lines changed: 126 additions & 121 deletions
Original file line numberDiff line numberDiff line change
@@ -45,7 +45,7 @@
4545
#include "TRandom3.h"
4646

4747
#include <string>
48-
#include <vector>
48+
#include <vector>
4949

5050
using namespace o2;
5151
using namespace o2::framework;
@@ -88,10 +88,10 @@ struct NetchargeFluctuations {
8888
Configurable<double> threshold{"threshold", 1e-6, "Delta eta bin count"};
8989

9090
// Event selections
91-
Configurable<bool> cSel8Trig{"cSel8Trig", true, "Sel8 (T0A + T0C) Selection Run3"}; // sel8
92-
Configurable<bool> cInt7Trig{"cInt7Trig", true, "kINT7 MB Trigger"}; // kINT7
93-
Configurable<bool> cSel7Trig{"cSel7Trig", true, "Sel7 (V0A + V0C) Selection Run2"}; // sel7
94-
Configurable<bool> cDcaXy{"cDcaXy", true, "Dca XY cut"};
91+
Configurable<bool> cSel8Trig{"cSel8Trig", true, "Sel8 (T0A + T0C) Selection Run3"}; // sel8
92+
Configurable<bool> cInt7Trig{"cInt7Trig", true, "kINT7 MB Trigger"}; // kINT7
93+
Configurable<bool> cSel7Trig{"cSel7Trig", true, "Sel7 (V0A + V0C) Selection Run2"}; // sel7
94+
Configurable<bool> cDcaXy{"cDcaXy", true, "Dca XY cut"};
9595
Configurable<bool> cDcaZ{"cDcaZ", true, "Dca Z cut"};
9696
Configurable<bool> cTpcCr{"cTpcCr", true, "tpc crossrows"};
9797
Configurable<bool> cItsChi{"cItsChi", true, "ITS chi"};
@@ -101,67 +101,67 @@ struct NetchargeFluctuations {
101101
ConfigurableAxis centBining{"centBining", {0, 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100}, "Centrality/Multiplicity percentile bining"};
102102
Configurable<bool> cPileupReject{"cPileupReject", true, "Pileup rejection"}; // pileup
103103
Configurable<bool> cfgUseGoodItsLayerAllCut{"cfgUseGoodItsLayerAllCut", true, "Good ITS Layers All"}; // pileup
104-
Configurable<bool> cTFBorder{"cTFBorder", false, "Timeframe Border Selection"}; // pileup
105-
Configurable<bool> cNoItsROBorder{"cNoItsROBorder", false, "No ITSRO Border Cut"}; // pileup
106-
Configurable<bool> cItsTpcVtx{"cItsTpcVtx", true, "ITS+TPC Vertex Selection"}; // pileup
107-
Configurable<bool> cZVtxTimeDiff{"cZVtxTimeDiff", false, "z-vtx time diff selection"}; // pileup
108-
Configurable<bool> cPVcont{"cPVcont", true, "primary vertex contributor"};
104+
Configurable<bool> cTFBorder{"cTFBorder", false, "Timeframe Border Selection"}; // pileup
105+
Configurable<bool> cNoItsROBorder{"cNoItsROBorder", false, "No ITSRO Border Cut"}; // pileup
106+
Configurable<bool> cItsTpcVtx{"cItsTpcVtx", true, "ITS+TPC Vertex Selection"}; // pileup
107+
Configurable<bool> cZVtxTimeDiff{"cZVtxTimeDiff", false, "z-vtx time diff selection"}; // pileup
108+
Configurable<bool> cPVcont{"cPVcont", true, "primary vertex contributor"};
109109

110110
O2_DEFINE_CONFIGURABLE(cfgEvSelMultCorrelation, bool, true, "Multiplicity correlation cut")
111-
struct : ConfigurableGroup {
111+
struct : ConfigurableGroup {
112112

113-
O2_DEFINE_CONFIGURABLE(cfgMultPVT0CCutEnabled, bool, true, "Enable PV multiplicity vs T0C centrality cut")
114-
O2_DEFINE_CONFIGURABLE(cfgMultGlobalFT0CCutEnabled, bool, true, "Enable globalTracks vs FT0C multiplicity cut")
115-
O2_DEFINE_CONFIGURABLE(cfgMultGlobalPVCutEnabled, bool, false, "Enable globalTracks vs PV multiplicity cut")
116-
117-
Configurable<std::vector<double>> cfgMultPVT0CCutPars{ "cfgMultPVT0CCutPars",
118-
std::vector<double>{187.621, -5.14575, 0.0716601, -0.000586642, 2.02818e-06, 51.2929, -1.66644, 0.0354762, -0.000389809, 1.55365e-06},
119-
"PV multiplicity vs T0C centrality cut parameter values"};
113+
O2_DEFINE_CONFIGURABLE(cfgMultPVT0CCutEnabled, bool, true, "Enable PV multiplicity vs T0C centrality cut")
114+
O2_DEFINE_CONFIGURABLE(cfgMultGlobalFT0CCutEnabled, bool, true, "Enable globalTracks vs FT0C multiplicity cut")
115+
O2_DEFINE_CONFIGURABLE(cfgMultGlobalPVCutEnabled, bool, false, "Enable globalTracks vs PV multiplicity cut")
120116

121-
Configurable<std::vector<double>> cfgMultGlobalFT0CCutPars{ "cfgMultGlobalFT0CCutPars",
122-
std::vector<double>{135.561, -3.7818, 0.0536562, -0.000445155, 1.55429e-06, 38.2336, -1.2568, 0.0270932, -0.000301034, 1.21234e-06},
123-
"globalTracks vs FT0C cut parameter values"};
117+
Configurable<std::vector<double>> cfgMultPVT0CCutPars{"cfgMultPVT0CCutPars",
118+
std::vector<double>{187.621, -5.14575, 0.0716601, -0.000586642, 2.02818e-06, 51.2929, -1.66644, 0.0354762, -0.000389809, 1.55365e-06},
119+
"PV multiplicity vs T0C centrality cut parameter values"};
124120

125-
Configurable<std::vector<double>> cfgMultGlobalPVCutPars{ "cfgMultGlobalPVCutPars",
126-
std::vector<double>{100., -2., 0.05, -0.0003, 1e-06, 30., -1.0, 0.02, -0.0002, 8e-07},
127-
"globalTracks vs PV cut parameter values"};
121+
Configurable<std::vector<double>> cfgMultGlobalFT0CCutPars{"cfgMultGlobalFT0CCutPars",
122+
std::vector<double>{135.561, -3.7818, 0.0536562, -0.000445155, 1.55429e-06, 38.2336, -1.2568, 0.0270932, -0.000301034, 1.21234e-06},
123+
"globalTracks vs FT0C cut parameter values"};
128124

129-
std::vector<double> multPVT0CCutPars;
130-
std::vector<double> multGlobalFT0CPars;
131-
std::vector<double> multGlobalPVCutPars;
125+
Configurable<std::vector<double>> cfgMultGlobalPVCutPars{"cfgMultGlobalPVCutPars",
126+
std::vector<double>{100., -2., 0.05, -0.0003, 1e-06, 30., -1.0, 0.02, -0.0002, 8e-07},
127+
"globalTracks vs PV cut parameter values"};
132128

133-
TF1* fMultPVT0CCutLow = nullptr;
134-
TF1* fMultPVT0CCutHigh = nullptr;
135-
TF1* fMultGlobalFT0CCutLow = nullptr;
136-
TF1* fMultGlobalFT0CCutHigh = nullptr;
137-
TF1* fMultGlobalPVCutLow = nullptr;
138-
TF1* fMultGlobalPVCutHigh = nullptr;
129+
std::vector<double> multPVT0CCutPars;
130+
std::vector<double> multGlobalFT0CPars;
131+
std::vector<double> multGlobalPVCutPars;
139132

140-
} cfgFunCoeff;
133+
TF1* fMultPVT0CCutLow = nullptr;
134+
TF1* fMultPVT0CCutHigh = nullptr;
135+
TF1* fMultGlobalFT0CCutLow = nullptr;
136+
TF1* fMultGlobalFT0CCutHigh = nullptr;
137+
TF1* fMultGlobalPVCutLow = nullptr;
138+
TF1* fMultGlobalPVCutHigh = nullptr;
139+
140+
} cfgFunCoeff;
141141

142142
// CCDB efficiency histograms
143143
TH1D* efficiency = nullptr;
144-
144+
145145
Filter collisionFilter = nabs(aod::collision::posZ) <= vertexZcut;
146-
Filter trackFilter = (nabs(aod::track::eta) < etaCut) && (aod::track::pt > ptMinCut ) && (aod::track::pt < ptMaxCut) && (requireGlobalTrackInFilter());
146+
Filter trackFilter = (nabs(aod::track::eta) < etaCut) && (aod::track::pt > ptMinCut) && (aod::track::pt < ptMaxCut) && (requireGlobalTrackInFilter());
147147

148-
using MyCollisionsRun2 = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::CentRun2V0Ms, aod::Mults>>;
149-
using MyCollisionRun2 = MyCollisionsRun2::iterator;
148+
using MyCollisionsRun2 = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::CentRun2V0Ms, aod::Mults>>;
149+
using MyCollisionRun2 = MyCollisionsRun2::iterator;
150150

151-
using MyCollisionsRun3 = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms, aod::CentFT0Cs, aod::Mults>>;
152-
using MyCollisionRun3 = MyCollisionsRun3::iterator;
151+
using MyCollisionsRun3 = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms, aod::CentFT0Cs, aod::Mults>>;
152+
using MyCollisionRun3 = MyCollisionsRun3::iterator;
153153

154-
using MyTracks = soa::Filtered<soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::TrackSelection>>;
155-
using MyTrack = MyTracks::iterator;
154+
using MyTracks = soa::Filtered<soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::TrackSelection>>;
155+
using MyTrack = MyTracks::iterator;
156156

157-
using MyMCCollisionsRun2 =soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::CentRun2V0Ms, aod::Mults, aod::McCollisionLabels>>;
158-
using MyMCCollisionRun2 = MyMCCollisionsRun2::iterator;
157+
using MyMCCollisionsRun2 = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::CentRun2V0Ms, aod::Mults, aod::McCollisionLabels>>;
158+
using MyMCCollisionRun2 = MyMCCollisionsRun2::iterator;
159159

160-
using MyMCCollisionsRun3 = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms, aod::CentFT0Cs, aod::Mults, aod::McCollisionLabels>>;
161-
using MyMCCollisionRun3 = MyMCCollisionsRun3::iterator;
160+
using MyMCCollisionsRun3 = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms, aod::CentFT0Cs, aod::Mults, aod::McCollisionLabels>>;
161+
using MyMCCollisionRun3 = MyMCCollisionsRun3::iterator;
162162

163-
using MyMCTracks = soa::Filtered<soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::StoredTracks, aod::TrackSelection, aod::McTrackLabels>>;
164-
using MyMCTrack = MyMCTracks::iterator;
163+
using MyMCTracks = soa::Filtered<soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::StoredTracks, aod::TrackSelection, aod::McTrackLabels>>;
164+
using MyMCTrack = MyMCTracks::iterator;
165165

166166
void init(o2::framework::InitContext&)
167167
{
@@ -325,44 +325,44 @@ Configurable<std::vector<double>> cfgMultGlobalPVCutPars{ "cfgMultGlobalPVCutPa
325325

326326
cfgFunCoeff.multPVT0CCutPars = cfgFunCoeff.cfgMultPVT0CCutPars;
327327
cfgFunCoeff.multGlobalFT0CPars = cfgFunCoeff.cfgMultGlobalFT0CCutPars;
328-
cfgFunCoeff.multGlobalPVCutPars = cfgFunCoeff.cfgMultGlobalPVCutPars;
328+
cfgFunCoeff.multGlobalPVCutPars = cfgFunCoeff.cfgMultGlobalPVCutPars;
329329

330-
cfgFunCoeff.fMultPVT0CCutLow =
330+
cfgFunCoeff.fMultPVT0CCutLow =
331331
new TF1("fMultPVT0CCutLow",
332332
"[0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x - 2.0*([5]+[6]*x+[7]*x*x+[8]*x*x*x+[9]*x*x*x*x)",
333333
0, 100);
334-
cfgFunCoeff.fMultPVT0CCutLow->SetParameters(&(cfgFunCoeff.multPVT0CCutPars[0]));
334+
cfgFunCoeff.fMultPVT0CCutLow->SetParameters(&(cfgFunCoeff.multPVT0CCutPars[0]));
335335

336-
cfgFunCoeff.fMultPVT0CCutHigh =
336+
cfgFunCoeff.fMultPVT0CCutHigh =
337337
new TF1("fMultPVT0CCutHigh",
338338
"[0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x + 2.*([5]+[6]*x+[7]*x*x+[8]*x*x*x+[9]*x*x*x*x)",
339339
0, 100);
340-
cfgFunCoeff.fMultPVT0CCutHigh->SetParameters(&(cfgFunCoeff.multPVT0CCutPars[0]));
340+
cfgFunCoeff.fMultPVT0CCutHigh->SetParameters(&(cfgFunCoeff.multPVT0CCutPars[0]));
341341

342-
cfgFunCoeff.fMultGlobalFT0CCutLow =
343-
new TF1("fMultGlobalFT0CCutLow",
344-
"[0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x - 2.*([5]+[6]*x+[7]*x*x+[8]*x*x*x+[9]*x*x*x*x)",
345-
0, 100);
346-
cfgFunCoeff.fMultGlobalFT0CCutLow->SetParameters(&(cfgFunCoeff.multGlobalFT0CPars[0]));
342+
cfgFunCoeff.fMultGlobalFT0CCutLow =
343+
new TF1("fMultGlobalFT0CCutLow",
344+
"[0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x - 2.*([5]+[6]*x+[7]*x*x+[8]*x*x*x+[9]*x*x*x*x)",
345+
0, 100);
346+
cfgFunCoeff.fMultGlobalFT0CCutLow->SetParameters(&(cfgFunCoeff.multGlobalFT0CPars[0]));
347347

348-
cfgFunCoeff.fMultGlobalFT0CCutHigh =
349-
new TF1("fMultGlobalFT0CCutHigh",
350-
"[0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x + 2.*([5]+[6]*x+[7]*x*x+[8]*x*x*x+[9]*x*x*x*x)",
351-
0, 100);
352-
cfgFunCoeff.fMultGlobalFT0CCutHigh->SetParameters(&(cfgFunCoeff.multGlobalFT0CPars[0]));
348+
cfgFunCoeff.fMultGlobalFT0CCutHigh =
349+
new TF1("fMultGlobalFT0CCutHigh",
350+
"[0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x + 2.*([5]+[6]*x+[7]*x*x+[8]*x*x*x+[9]*x*x*x*x)",
351+
0, 100);
352+
cfgFunCoeff.fMultGlobalFT0CCutHigh->SetParameters(&(cfgFunCoeff.multGlobalFT0CPars[0]));
353353

354-
cfgFunCoeff.fMultGlobalPVCutLow =
354+
cfgFunCoeff.fMultGlobalPVCutLow =
355355
new TF1("fMultGlobalPVCutLow",
356-
"[0] + [1]*x - 5.*([2] + [3]*x)",
356+
"[0] + [1]*x - 5.*([2] + [3]*x)",
357357
0, 100);
358-
cfgFunCoeff.fMultGlobalPVCutLow->SetParameters(&(cfgFunCoeff.multGlobalPVCutPars[0]));
358+
cfgFunCoeff.fMultGlobalPVCutLow->SetParameters(&(cfgFunCoeff.multGlobalPVCutPars[0]));
359359

360-
cfgFunCoeff.fMultGlobalPVCutHigh =
360+
cfgFunCoeff.fMultGlobalPVCutHigh =
361361
new TF1("fMultGlobalPVCutHigh",
362-
"[0] + [1]*x + 5.*([2] + [3]*x)",
362+
"[0] + [1]*x + 5.*([2] + [3]*x)",
363363
0, 100);
364-
cfgFunCoeff.fMultGlobalPVCutHigh->SetParameters(&(cfgFunCoeff.multGlobalPVCutPars[0]));
365-
364+
cfgFunCoeff.fMultGlobalPVCutHigh->SetParameters(&(cfgFunCoeff.multGlobalPVCutPars[0]));
365+
366366
if (cfgLoadEff) {
367367
ccdb->setURL(cfgUrlCCDB.value);
368368
ccdb->setCaching(true);
@@ -376,34 +376,34 @@ cfgFunCoeff.fMultGlobalFT0CCutHigh->SetParameters(&(cfgFunCoeff.multGlobalFT0CPa
376376
}
377377
}
378378

379-
bool eventSelected(const float& globalNch, const float& pvTrack, const float& centrality)
380-
{
381-
if (cfgFunCoeff.cfgMultPVT0CCutEnabled) {
379+
bool eventSelected(const float& globalNch, const float& pvTrack, const float& centrality)
380+
{
381+
if (cfgFunCoeff.cfgMultPVT0CCutEnabled) {
382382

383-
if (pvTrack < cfgFunCoeff.fMultPVT0CCutLow->Eval(centrality))
384-
return false;
385-
if (pvTrack > cfgFunCoeff.fMultPVT0CCutHigh->Eval(centrality))
386-
return false;
387-
}
383+
if (pvTrack < cfgFunCoeff.fMultPVT0CCutLow->Eval(centrality))
384+
return false;
385+
if (pvTrack > cfgFunCoeff.fMultPVT0CCutHigh->Eval(centrality))
386+
return false;
387+
}
388388

389-
if (cfgFunCoeff.cfgMultGlobalFT0CCutEnabled) {
389+
if (cfgFunCoeff.cfgMultGlobalFT0CCutEnabled) {
390390

391-
if (globalNch < cfgFunCoeff.fMultGlobalFT0CCutLow->Eval(centrality))
392-
return false;
393-
if (globalNch > cfgFunCoeff.fMultGlobalFT0CCutHigh->Eval(centrality))
394-
return false;
395-
}
391+
if (globalNch < cfgFunCoeff.fMultGlobalFT0CCutLow->Eval(centrality))
392+
return false;
393+
if (globalNch > cfgFunCoeff.fMultGlobalFT0CCutHigh->Eval(centrality))
394+
return false;
395+
}
396396

397-
if (cfgFunCoeff.cfgMultGlobalPVCutEnabled) {
398-
399-
if (globalNch < cfgFunCoeff.fMultGlobalPVCutLow->Eval(pvTrack))
400-
return false;
401-
if (globalNch > cfgFunCoeff.fMultGlobalPVCutHigh->Eval(pvTrack))
402-
return false;
403-
}
397+
if (cfgFunCoeff.cfgMultGlobalPVCutEnabled) {
404398

405-
return true;
406-
}
399+
if (globalNch < cfgFunCoeff.fMultGlobalPVCutLow->Eval(pvTrack))
400+
return false;
401+
if (globalNch > cfgFunCoeff.fMultGlobalPVCutHigh->Eval(pvTrack))
402+
return false;
403+
}
404+
405+
return true;
406+
}
407407

408408
template <RunType run, typename C>
409409
bool selCollision(C const& coll, float& cent, float& mult)
@@ -414,7 +414,7 @@ bool eventSelected(const float& globalNch, const float& pvTrack, const float& ce
414414
if constexpr (run == kRun3) {
415415
if (cSel8Trig && !coll.sel8()) {
416416
return false;
417-
}
417+
}
418418
if (cFT0M) {
419419
cent = coll.centFT0M(); // centrality for run3 using FT0M
420420
mult = coll.multFT0M();
@@ -483,7 +483,7 @@ bool eventSelected(const float& globalNch, const float& pvTrack, const float& ce
483483
{
484484
if (!track.isGlobalTrack())
485485
return false;
486-
if (cPVcont && !track.isPVContributor())
486+
if (cPVcont && !track.isPVContributor())
487487
return false;
488488
if (std::fabs(track.eta()) >= etaCut)
489489
return false;
@@ -549,28 +549,29 @@ bool eventSelected(const float& globalNch, const float& pvTrack, const float& ce
549549
histogramRegistry.fill(HIST("subsample/posneg"), cent, sampleIndex, posneg);
550550
}
551551

552-
553552
template <RunType run, typename C, typename T>
554553
void calculationData(C const& coll, T const& tracks)
555554
{
556-
float cent = -1, mult = -1;
555+
float cent = -1, mult = -1;
557556
histogramRegistry.fill(HIST("QA/hVtxZ_before"), coll.posZ());
558557
if (!selCollision<run>(coll, cent, mult)) {
559558
return;
560559
}
561-
562-
float globalNch = tracks.size();
563-
float pvTrack = coll.multNTracksPV();
560+
561+
float globalNch = tracks.size();
562+
float pvTrack = coll.multNTracksPV();
564563

565564
histogramRegistry.fill(HIST("QA/hCentFT0C"), cent);
566-
histogramRegistry.fill(HIST("QA/hNchGlobal"), globalNch);
567-
histogramRegistry.fill(HIST("QA/hNchPV"), pvTrack);
568-
565+
histogramRegistry.fill(HIST("QA/hNchGlobal"), globalNch);
566+
histogramRegistry.fill(HIST("QA/hNchPV"), pvTrack);
567+
569568
histogramRegistry.fill(HIST("MultCorrelationPlots/globalTracks_PV_bef"), pvTrack, globalNch);
570569
histogramRegistry.fill(HIST("MultCorrelationPlots/globalTracks_FT0C_bef"), cent, globalNch);
571-
histogramRegistry.fill(HIST("MultCorrelationPlots/PV_FT0C_bef"), cent, pvTrack);
572-
573-
if(cfgEvSelMultCorrelation && !eventSelected(globalNch, pvTrack, cent)) { return;}
570+
histogramRegistry.fill(HIST("MultCorrelationPlots/PV_FT0C_bef"), cent, pvTrack);
571+
572+
if (cfgEvSelMultCorrelation && !eventSelected(globalNch, pvTrack, cent)) {
573+
return;
574+
}
574575
histogramRegistry.fill(HIST("MultCorrelationPlots/globalTracks_PV_aft"), pvTrack, globalNch);
575576
histogramRegistry.fill(HIST("MultCorrelationPlots/globalTracks_FT0C_aft"), cent, globalNch);
576577
histogramRegistry.fill(HIST("MultCorrelationPlots/PV_FT0C_aft"), cent, pvTrack);
@@ -613,7 +614,7 @@ bool eventSelected(const float& globalNch, const float& pvTrack, const float& ce
613614
termp = fpos * (fpos - 1);
614615
termn = fneg * (fneg - 1);
615616
posneg = fpos * fneg;
616-
617+
617618
histogramRegistry.fill(HIST("data/cent_nchTotal"), cent, nchTotal);
618619
histogramRegistry.fill(HIST("data/cent_nchTotalCor"), cent, nchTotalCor);
619620
histogramRegistry.fill(HIST("data/nch_nchCor"), nch, nchCor);
@@ -638,18 +639,20 @@ bool eventSelected(const float& globalNch, const float& pvTrack, const float& ce
638639
return;
639640
}
640641

641-
int globalNch = inputTracks.size();
642-
int pvTrack = coll.multNTracksPV();
642+
int globalNch = inputTracks.size();
643+
int pvTrack = coll.multNTracksPV();
643644

644645
histogramRegistry.fill(HIST("QA/hCentFT0C"), cent);
645-
histogramRegistry.fill(HIST("QA/hNchGlobal"), globalNch);
646-
histogramRegistry.fill(HIST("QA/hNchPV"), pvTrack);
647-
646+
histogramRegistry.fill(HIST("QA/hNchGlobal"), globalNch);
647+
histogramRegistry.fill(HIST("QA/hNchPV"), pvTrack);
648+
648649
histogramRegistry.fill(HIST("MultCorrelationPlots/globalTracks_PV_bef"), pvTrack, globalNch);
649650
histogramRegistry.fill(HIST("MultCorrelationPlots/globalTracks_FT0C_bef"), cent, globalNch);
650-
histogramRegistry.fill(HIST("MultCorrelationPlots/PV_FT0C_bef"), cent, pvTrack);
651-
652-
if(cfgEvSelMultCorrelation && !eventSelected(globalNch, pvTrack, cent)) { return;}
651+
histogramRegistry.fill(HIST("MultCorrelationPlots/PV_FT0C_bef"), cent, pvTrack);
652+
653+
if (cfgEvSelMultCorrelation && !eventSelected(globalNch, pvTrack, cent)) {
654+
return;
655+
}
653656
histogramRegistry.fill(HIST("MultCorrelationPlots/globalTracks_PV_aft"), pvTrack, globalNch);
654657
histogramRegistry.fill(HIST("MultCorrelationPlots/globalTracks_FT0C_aft"), cent, globalNch);
655658
histogramRegistry.fill(HIST("MultCorrelationPlots/PV_FT0C_aft"), cent, pvTrack);
@@ -775,9 +778,10 @@ bool eventSelected(const float& globalNch, const float& pvTrack, const float& ce
775778
if (!selCollision<run>(coll, cent, mult))
776779
return;
777780

778-
int globalNch = tracks.size();
779-
int pvTrack = coll.multNTracksPV();
780-
if(cfgEvSelMultCorrelation && !eventSelected(globalNch, pvTrack, cent)) return;
781+
int globalNch = tracks.size();
782+
int pvTrack = coll.multNTracksPV();
783+
if (cfgEvSelMultCorrelation && !eventSelected(globalNch, pvTrack, cent))
784+
return;
781785

782786
if (!(cent >= centMin && cent < centMax))
783787
return;
@@ -852,9 +856,10 @@ bool eventSelected(const float& globalNch, const float& pvTrack, const float& ce
852856
if (!selCollision<run>(coll, cent, mult))
853857
return;
854858

855-
int globalNch = inputTracks.size();
856-
int pvTrack = coll.multNTracksPV();
857-
if(cfgEvSelMultCorrelation && !eventSelected(globalNch, pvTrack, cent)) return;
859+
int globalNch = inputTracks.size();
860+
int pvTrack = coll.multNTracksPV();
861+
if (cfgEvSelMultCorrelation && !eventSelected(globalNch, pvTrack, cent))
862+
return;
858863

859864
if (!(cent >= centMin && cent < centMax))
860865
return;

0 commit comments

Comments
 (0)