@@ -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,21 +546,22 @@ 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;
551551
552552 // number of parameters required in cfgEvSelsMultPv and cfgEvSelsMult. (5 Low + 5 High)
553- int nParams = 10 ;
554-
555- if (paramsMultPVCut.size () < nParams) LOGF (fatal, " cfgEvSelsMultPv not set properly.. size = %d (should be 10) --> Check your config files!" , paramsMultPVCut.size ());
556- else if (paramsMultCut.size () < nParams) LOGF (fatal, " cfgEvSelsMult not set properly.. size = %d (should be 10) --> Check your config files!" , paramsMultCut.size ());
553+ int nParams = 10 ;
554+
555+ if (paramsMultPVCut.size () < nParams)
556+ LOGF (fatal, " cfgEvSelsMultPv not set properly.. size = %d (should be 10) --> Check your config files!" , paramsMultPVCut.size ());
557+ else if (paramsMultCut.size () < nParams)
558+ LOGF (fatal, " cfgEvSelsMult not set properly.. size = %d (should be 10) --> Check your config files!" , paramsMultCut.size ());
557559 else {
558560 fMultPVCutLow ->SetParameters (paramsMultPVCut[0 ], paramsMultPVCut[1 ], paramsMultPVCut[2 ], paramsMultPVCut[3 ], paramsMultPVCut[4 ]);
559561 fMultPVCutHigh ->SetParameters (paramsMultPVCut[5 ], paramsMultPVCut[6 ], paramsMultPVCut[7 ], paramsMultPVCut[8 ], paramsMultPVCut[9 ]);
560562 fMultCutLow ->SetParameters (paramsMultCut[0 ], paramsMultCut[1 ], paramsMultCut[2 ], paramsMultCut[3 ], paramsMultCut[4 ]);
561563 fMultCutHigh ->SetParameters (paramsMultCut[5 ], paramsMultCut[6 ], paramsMultCut[7 ], paramsMultCut[8 ], paramsMultCut[9 ]);
562564 }
563-
564565 }
565566
566567 if (cfgTrackSelsUseAdditionalTrackCut) {
@@ -686,10 +687,9 @@ struct FlowSP {
686687 return false ;
687688 weight_nue = 1 . / eff;
688689 int sizeAcc = cfg.mAcceptance .size ();
689- if (sizeAcc > pID){
690+ if (sizeAcc > pID) {
690691 weight_nua = cfg.mAcceptance [pID]->getNUA (phi, eta, vtxz);
691- }
692- else
692+ } else
693693 weight_nua = 1 ;
694694 return true ;
695695 }
@@ -1010,17 +1010,17 @@ struct FlowSP {
10101010
10111011 static constexpr std::string_view Time[] = {" before/" , " after/" };
10121012 // NOTE: species[kUnidentified] = "" (when no PID)
1013- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hPt" ), track.pt (), wacc * weff);
1014- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hPhi" ), track.phi (), wacc);
1015- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hPhi_uncorrected" ), track.phi ());
1016- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hEta" ), track.eta (), wacc);
1017- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hPhi_Eta_vz" ), track.phi (), track.eta (), vz);
1018- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hPhi_Eta_vz_corrected" ), track.phi (), track.eta (), vz, wacc);
1019- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hDCAxy_pt" ), track.pt (), track.dcaXY (), wacc * weff);
1020- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hDCAz_pt" ), track.pt (), track.dcaZ (), wacc * weff);
1021- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hSharedClusters_pt" ), track.pt (), track.tpcFractionSharedCls (), wacc * weff);
1022- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hCrossedRows_pt" ), track.pt (), track.tpcNClsFound (), wacc * weff);
1023- registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hCrossedRows_vs_SharedClusters" ), track.tpcNClsFound (), track.tpcFractionSharedCls (), wacc * weff);
1013+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hPt" ), track.pt (), wacc * weff);
1014+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hPhi" ), track.phi (), wacc);
1015+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hPhi_uncorrected" ), track.phi ());
1016+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hEta" ), track.eta (), wacc);
1017+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hPhi_Eta_vz" ), track.phi (), track.eta (), vz);
1018+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hPhi_Eta_vz_corrected" ), track.phi (), track.eta (), vz, wacc);
1019+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hDCAxy_pt" ), track.pt (), track.dcaXY (), wacc * weff);
1020+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hDCAz_pt" ), track.pt (), track.dcaZ (), wacc * weff);
1021+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hSharedClusters_pt" ), track.pt (), track.tpcFractionSharedCls (), wacc * weff);
1022+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hCrossedRows_pt" ), track.pt (), track.tpcNClsFound (), wacc * weff);
1023+ registry.fill (HIST (Charge[ct]) + HIST (Species[pt]) + HIST (" QA/" ) + HIST (Time[ft]) + HIST (" hCrossedRows_vs_SharedClusters" ), track.tpcNClsFound (), track.tpcFractionSharedCls (), wacc * weff);
10241024 }
10251025
10261026 template <FillType ft, ChargeType ct, typename TrackObject>
0 commit comments