@@ -117,7 +117,7 @@ struct FlowSP {
117117 // Connect to ccdb
118118 Service<ccdb::BasicCCDBManager> ccdb;
119119
120- int counter = 0 ;
120+ int counter = 0 ;
121121
122122 // struct to hold the correction histos/
123123 struct Config {
@@ -370,30 +370,30 @@ struct FlowSP {
370370 // Fitted for LHC23zzh_pass4
371371 fMultPVCutLow = new TF1 (" fMultPVCutLow" , " [0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x" , 0 , 100 );
372372 fMultPVCutLow ->SetParameters (2942.55 , -103.111 , 1.4397 , -0.00974862 , 2.71433e-05 );
373- if (cfgnSigmaMultCuts == 2 )
373+ if (cfgnSigmaMultCuts == 2 )
374374 fMultPVCutLow ->SetParameters (2665.68 , -93.3784 , 1.27137 , -0.00818936 , 2.115e-05 );
375- if (cfgnSigmaMultCuts == 3 )
375+ if (cfgnSigmaMultCuts == 3 )
376376 fMultPVCutLow ->SetParameters (2389.99 , -83.8483 , 1.11062 , -0.00672263 , 1.54725e-05 );
377377
378378 fMultPVCutHigh = new TF1 (" fMultPVCutHigh" , " [0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x" , 0 , 100 );
379379 fMultPVCutHigh ->SetParameters (3508.13 , -124.831 , 1.87871 , -0.0145343 , 4.80688e-05 );
380- if (cfgnSigmaMultCuts == 2 )
380+ if (cfgnSigmaMultCuts == 2 )
381381 fMultPVCutHigh ->SetParameters (3787.93 , -135.184 , 2.07683 , -0.0165997 , 5.68725e-05 );
382- if (cfgnSigmaMultCuts == 3 )
382+ if (cfgnSigmaMultCuts == 3 )
383383 fMultPVCutHigh ->SetParameters (4067.4 , -145.485 , 2.27273 , -0.0186308 , 6.5501e-05 );
384384
385385 fMultCutLow = new TF1 (" fMultCutLow" , " [0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x" , 0 , 100 );
386386 fMultCutLow ->SetParameters (1566.5 , -48.2114 , 0.529522 , -0.00235284 , 3.01132e-06 );
387- if (cfgnSigmaMultCuts == 2 )
387+ if (cfgnSigmaMultCuts == 2 )
388388 fMultCutLow ->SetParameters (1307.92 , -39.9168 , 0.412675 , -0.00148081 , 1.10868e-07 );
389- if (cfgnSigmaMultCuts == 3 )
389+ if (cfgnSigmaMultCuts == 3 )
390390 fMultCutLow ->SetParameters (1048.48 , -31.4568 , 0.287794 , -0.00046847 , -3.5909e-06 );
391391
392392 fMultCutHigh = new TF1 (" fMultCutHigh" , " [0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x" , 0 , 100 );
393393 fMultCutHigh ->SetParameters (2089.73 , -65.9772 , 0.816781 , -0.00496563 , 1.34314e-05 );
394- if (cfgnSigmaMultCuts == 2 )
394+ if (cfgnSigmaMultCuts == 2 )
395395 fMultCutHigh ->SetParameters (2350.39 , -74.6939 , 0.953287 , -0.006162 , 1.80808e-05 );
396- if (cfgnSigmaMultCuts == 3 )
396+ if (cfgnSigmaMultCuts == 3 )
397397 fMultCutHigh ->SetParameters (2610.98 , -83.3983 , 1.0893 , -0.00735094 , 2.26929e-05 );
398398 }
399399
@@ -596,13 +596,15 @@ struct FlowSP {
596596
597597 phimodn += o2::constants::math::PI / 18.0 ; // to center gap in the middle
598598 phimodn = fmod (phimodn, o2::constants::math::PI / 9.0 );
599- if (cfgFillQAHistos) registry.fill (HIST (" QA/before/pt_phi" ), track.pt (), phimodn);
599+ if (cfgFillQAHistos)
600+ registry.fill (HIST (" QA/before/pt_phi" ), track.pt (), phimodn);
600601
601602 if (cfgUseAdditionalTrackCut) {
602603 if (phimodn < fPhiCutHigh ->Eval (track.pt ()) && phimodn > fPhiCutLow ->Eval (track.pt ()))
603604 return false ; // reject track
604605 }
605- if (cfgFillQAHistos) registry.fill (HIST (" QA/after/pt_phi" ), track.pt (), phimodn);
606+ if (cfgFillQAHistos)
607+ registry.fill (HIST (" QA/after/pt_phi" ), track.pt (), phimodn);
606608 registry.fill (HIST (" hTrackCount" ), trackSel_TPCBoundary);
607609 return true ;
608610 }
@@ -833,7 +835,7 @@ struct FlowSP {
833835 }
834836
835837 for (const auto & track : tracks) {
836- if (cfgFillQAHistos)
838+ if (cfgFillQAHistos)
837839 registry.fill (HIST (" QA/before/hPt_inclusive" ), track.pt ());
838840
839841 registry.fill (HIST (" hTrackCount" ), trackSel_FilteredTracks);
@@ -847,7 +849,7 @@ struct FlowSP {
847849
848850 if (track.sign () == 0.0 )
849851 continue ;
850-
852+
851853 registry.fill (HIST (" hTrackCount" ), trackSel_ZeroCharge);
852854 bool pos = (track.sign () > 0 ) ? true : false ;
853855
@@ -872,7 +874,7 @@ struct FlowSP {
872874
873875 registry.fill (HIST (" hTrackCount" ), trackSel_ParticleWeights);
874876
875- if (cfgFillQAHistos)
877+ if (cfgFillQAHistos)
876878 registry.fill (HIST (" QA/after/hPt_inclusive" ), track.pt (), wacc * weff);
877879
878880 // // constrain angle to 0 -> [0,0+2pi]
@@ -930,8 +932,8 @@ struct FlowSP {
930932 fillEventQA<kAfter >(collision, tracks);
931933
932934 for (const auto & track : tracks) {
933- double wacc = 1.0 ;
934- double weff = 1.0 ;
935+ double wacc = 1.0 ;
936+ double weff = 1.0 ;
935937
936938 auto mcParticle = track.mcParticle ();
937939 if (!mcParticle.isPhysicalPrimary ())
@@ -944,7 +946,7 @@ struct FlowSP {
944946 continue ;
945947 bool pos = (track.sign () > 0 ) ? true : false ;
946948
947- if (cfgFillQAHistos){
949+ if (cfgFillQAHistos) {
948950 registry.fill (HIST (" QA/before/hPt_inclusive" ), track.pt ());
949951 if (pos) {
950952 registry.fill (HIST (" QA/before/hPt_positive" ), track.pt ());
@@ -974,7 +976,7 @@ struct FlowSP {
974976 Filter mcCollFilter = nabs(aod::mccollision::posZ) < cfgVtxZ;
975977 void processMCGen (soa::Filtered<aod::McCollisions>::iterator const & mcCollision, soa::SmallGroups<soa::Join<aod::McCollisionLabels, aod::Collisions, aod::CentFT0Cs, aod::CentFV0As, aod::CentFT0CVariant1s, aod::CentFT0Ms, aod::CentNGlobals>> const & collisions, aod::McParticles const & particles)
976978 {
977- if (collisions.size () != 1 ){ // check if MC collision is only reconstructed once! (https://indico.cern.ch/event/1425820/contributions/6170879/attachments/2947721/5180548/DDChinellato-O2AT4-HandsOn-03a.pdf)
979+ if (collisions.size () != 1 ) { // check if MC collision is only reconstructed once! (https://indico.cern.ch/event/1425820/contributions/6170879/attachments/2947721/5180548/DDChinellato-O2AT4-HandsOn-03a.pdf)
978980 return ;
979981 }
980982 float centrality = -1 ;
@@ -995,7 +997,6 @@ struct FlowSP {
995997 if (centrality < cfgCentMin || centrality > cfgCentMax)
996998 return ;
997999
998-
9991000 float vtxz = mcCollision.posZ ();
10001001
10011002 for (const auto & particle : particles) {
@@ -1011,7 +1012,7 @@ struct FlowSP {
10111012 registry.fill (HIST (" trackMCGen/after/pt_gen_incl" ), particle.pt ());
10121013 registry.fill (HIST (" trackMCGen/after/phi_eta_vtxZ_gen" ), particle.phi (), particle.eta (), vtxz);
10131014 }
1014- counter++;
1015+ counter++;
10151016 }
10161017 PROCESS_SWITCH (FlowSP, processMCGen, " Process analysis for MC generated events" , false );
10171018};
0 commit comments