@@ -115,7 +115,7 @@ struct FlowSP {
115115 ConfigurableAxis axisCentrality{" axisCentrality" , {10 , 0 , 100 }, " Centrality bins for vn " };
116116 ConfigurableAxis axisNch = {" axisNch" , {400 , 0 , 4000 }, " Global N_{ch}" };
117117 ConfigurableAxis axisMultpv = {" axisMultpv" , {400 , 0 , 4000 }, " N_{ch} (PV)" };
118- // Configurables containing vector
118+ // Configurables containing vector
119119 Configurable<std::vector<double >> cfgEvSelsMultPv{" cfgEvSelsMultPv" , std::vector<double >{2389.99 , -83.8483 , 1.11062 , -0.00672263 , 1.54725e-05 , 4067.4 , -145.485 , 2.27273 , -0.0186308 , 6.5501e-05 }, " Multiplicity cuts (PV) first 5 parameters cutLOW last 5 cutHIGH (Default is +-3sigma pass4) " };
120120 Configurable<std::vector<double >> cfgEvSelsMult{" cfgEvSelsMult" , std::vector<double >{1048.48 , -31.4568 , 0.287794 , -0.00046847 , -3.5909e-06 , 2610.98 , -83.3983 , 1.0893 , -0.00735094 , 2.26929e-05 }, " Multiplicity cuts (Global) first 5 parameters cutLOW last 5 cutHIGH (Default is +-3sigma pass4) " };
121121
@@ -371,7 +371,7 @@ struct FlowSP {
371371 registry.add (" incl/QA/after/hCrossedRows_vs_SharedClusters" , " " , {HistType::kTH2D , {axisCl, axisShCl}});
372372
373373 if (cfgFillQABefore)
374- registry.addClone (" incl/QA/after/" , " incl/QA/before/" );
374+ registry.addClone (" incl/QA/after/" , " incl/QA/before/" );
375375 }
376376
377377 if (doprocessMCReco) {
@@ -546,18 +546,19 @@ struct FlowSP {
546546 fMultCutLow = new TF1 (" fMultCutLow" , " [0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x" , 0 , 100 );
547547 fMultCutHigh = new TF1 (" fMultCutHigh" , " [0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x" , 0 , 100 );
548548
549- std::vector<double > paramsMultPVCut = cfgEvSelsMultPv;
549+ std::vector<double > paramsMultPVCut = cfgEvSelsMultPv;
550550 std::vector<double > paramsMultCut = cfgEvSelsMult;
551-
552- if (paramsMultPVCut.size () < 10 ) LOGF (fatal, " cfgEvSelsMultPv not set properly.. size = %d (should be 10) --> Check your config files!" , paramsMultPVCut.size ());
553- else if (paramsMultCut.size () < 10 ) LOGF (fatal, " cfgEvSelsMult not set properly.. size = %d (should be 10) --> Check your config files!" , paramsMultCut.size ());
551+
552+ if (paramsMultPVCut.size () < 10 )
553+ LOGF (fatal, " cfgEvSelsMultPv not set properly.. size = %d (should be 10) --> Check your config files!" , paramsMultPVCut.size ());
554+ else if (paramsMultCut.size () < 10 )
555+ LOGF (fatal, " cfgEvSelsMult not set properly.. size = %d (should be 10) --> Check your config files!" , paramsMultCut.size ());
554556 else {
555557 fMultPVCutLow ->SetParameters (paramsMultPVCut[0 ], paramsMultPVCut[1 ], paramsMultPVCut[2 ], paramsMultPVCut[3 ], paramsMultPVCut[4 ]);
556558 fMultPVCutHigh ->SetParameters (paramsMultPVCut[5 ], paramsMultPVCut[6 ], paramsMultPVCut[7 ], paramsMultPVCut[8 ], paramsMultPVCut[9 ]);
557559 fMultCutLow ->SetParameters (paramsMultCut[0 ], paramsMultCut[1 ], paramsMultCut[2 ], paramsMultCut[3 ], paramsMultCut[4 ]);
558560 fMultCutHigh ->SetParameters (paramsMultCut[5 ], paramsMultCut[6 ], paramsMultCut[7 ], paramsMultCut[8 ], paramsMultCut[9 ]);
559561 }
560-
561562 }
562563
563564 if (cfgTrackSelsUseAdditionalTrackCut) {
@@ -683,10 +684,9 @@ struct FlowSP {
683684 return false ;
684685 weight_nue = 1 . / eff;
685686 int sizeAcc = cfg.mAcceptance .size ();
686- if (sizeAcc > pID){
687+ if (sizeAcc > pID) {
687688 weight_nua = cfg.mAcceptance [pID]->getNUA (phi, eta, vtxz);
688- }
689- else
689+ } else
690690 weight_nua = 1 ;
691691 return true ;
692692 }
@@ -1007,17 +1007,17 @@ struct FlowSP {
10071007
10081008 static constexpr std::string_view Time[] = {" before/" , " after/" };
10091009 // NOTE: species[kUnidentified] = "" (when no PID)
1010- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hPt" ), track.pt (), wacc * weff);
1011- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hPhi" ), track.phi (), wacc);
1012- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hPhi_uncorrected" ), track.phi ());
1013- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hEta" ), track.eta (), wacc);
1014- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hPhi_Eta_vz" ), track.phi (), track.eta (), vz);
1015- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hPhi_Eta_vz_corrected" ), track.phi (), track.eta (), vz, wacc);
1016- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hDCAxy_pt" ), track.pt (), track.dcaXY (), wacc * weff);
1017- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hDCAz_pt" ), track.pt (), track.dcaZ (), wacc * weff);
1018- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hSharedClusters_pt" ), track.pt (), track.tpcFractionSharedCls (), wacc * weff);
1019- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hCrossedRows_pt" ), track.pt (), track.tpcNClsFound (), wacc * weff);
1020- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hCrossedRows_vs_SharedClusters" ), track.tpcNClsFound (), track.tpcFractionSharedCls (), wacc * weff);
1010+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hPt" ), track.pt (), wacc * weff);
1011+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hPhi" ), track.phi (), wacc);
1012+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hPhi_uncorrected" ), track.phi ());
1013+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hEta" ), track.eta (), wacc);
1014+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hPhi_Eta_vz" ), track.phi (), track.eta (), vz);
1015+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hPhi_Eta_vz_corrected" ), track.phi (), track.eta (), vz, wacc);
1016+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hDCAxy_pt" ), track.pt (), track.dcaXY (), wacc * weff);
1017+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hDCAz_pt" ), track.pt (), track.dcaZ (), wacc * weff);
1018+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hSharedClusters_pt" ), track.pt (), track.tpcFractionSharedCls (), wacc * weff);
1019+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hCrossedRows_pt" ), track.pt (), track.tpcNClsFound (), wacc * weff);
1020+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hCrossedRows_vs_SharedClusters" ), track.tpcNClsFound (), track.tpcFractionSharedCls (), wacc * weff);
10211021 }
10221022
10231023 template <FillType ft, ChargeType ct, typename TrackObject>
0 commit comments