@@ -84,22 +84,22 @@ struct StrangeCascTrack {
8484 Configurable<bool > doApplyTPCPIDOmega{" doApplyTPCPIDOmega" , true , " apply tpc pid to dau tracks (Omega)" };
8585 Configurable<bool > doApplyTOFPIDOmega{" doApplyTOFPIDOmega" , true , " apply tof pid to dau tracks (Omega)" };
8686 Configurable<bool > doCompetingMassRej{" doCompetingMassRej" , true , " competing mass rejection (Omega)" };
87- // efficiency and purity corrections on the fly (warning: to be avoided because interpolation causes errors):
88- // only correct by pt
89- Configurable<bool > doApplyEfficiency1D{" doApplyEfficiency1D" , false , " apply efficiency correction" };
90- Configurable<bool > doPropagateEfficiency1D{" doPropagateEfficiency1D" , false , " apply efficiency propagation" };
91- Configurable<bool > doApplyPurity1D{" doApplyPurity1D" , false , " apply purity correction" };
92- Configurable<bool > doPropagatePurity1D{" doPropagatePurity1D" , false , " apply purity propagation" };
93- // correct by both pt and mult
94- Configurable<bool > doApplyEfficiency2D{" doApplyEfficiency2D" , false , " apply efficiency correction" };
95- Configurable<bool > doPropagateEfficiency2D{" doPropagateEfficiency2D" , false , " apply efficiency propagation" };
96- Configurable<bool > doApplyPurity2D{" doApplyPurity2D" , false , " apply purity correction" };
97- Configurable<bool > doPropagatePurity2D{" doPropagatePurity2D" , false , " apply purity propagation" };
98- Configurable<std::string> ccdburl{" ccdburl" , " http://alice-ccdb.cern.ch" , " url of the ccdb repository to use" };
99- Configurable<std::string> efficiencyCCDBPathpp{" efficiencyCCDBPathpp" , " Users/y/yparovia/LHC24f4d" , " Path of the efficiency and purity corrections (pp)" };
100- Configurable<std::string> efficiencyCCDBPathPbPb{" efficiencyCCDBPathPbPb" , " Users/y/yparovia/LHC25f3" , " Path of the efficiency and purity corrections (PbPb)" };
101- Configurable<std::string> efficiencyCCDBPathOO{" efficiencyCCDBPathOO" , " Users/y/yparovia/LHC25h3" , " Path of the efficiency and purity corrections (OO)" };
102- Configurable<std::string> efficiencyCCDBPathpO{" efficiencyCCDBPathpO" , " Users/y/yparovia/LHC25h2" , " Path of the efficiency and purity corrections (pO)" };
87+ // * // efficiency and purity corrections on the fly (warning: to be avoided because interpolation causes errors):
88+ // * // only correct by pt
89+ // * Configurable<bool> doApplyEfficiency1D{"doApplyEfficiency1D", false, "apply efficiency correction"};
90+ // * Configurable<bool> doPropagateEfficiency1D{"doPropagateEfficiency1D", false, "apply efficiency propagation"};
91+ // * Configurable<bool> doApplyPurity1D{"doApplyPurity1D", false, "apply purity correction"};
92+ // * Configurable<bool> doPropagatePurity1D{"doPropagatePurity1D", false, "apply purity propagation"};
93+ // * // correct by both pt and mult
94+ // * Configurable<bool> doApplyEfficiency2D{"doApplyEfficiency2D", false, "apply efficiency correction"};
95+ // * Configurable<bool> doPropagateEfficiency2D{"doPropagateEfficiency2D", false, "apply efficiency propagation"};
96+ // * Configurable<bool> doApplyPurity2D{"doApplyPurity2D", false, "apply purity correction"};
97+ // * Configurable<bool> doPropagatePurity2D{"doPropagatePurity2D", false, "apply purity propagation"};
98+ // * Configurable<std::string> ccdburl{"ccdburl", "http://alice-ccdb.cern.ch", "url of the ccdb repository to use"};
99+ // * Configurable<std::string> efficiencyCCDBPathpp{"efficiencyCCDBPathpp", "Users/y/yparovia/LHC24f4d", "Path of the efficiency and purity corrections (pp)"};
100+ // * Configurable<std::string> efficiencyCCDBPathPbPb{"efficiencyCCDBPathPbPb", "Users/y/yparovia/LHC25f3", "Path of the efficiency and purity corrections (PbPb)"};
101+ // * Configurable<std::string> efficiencyCCDBPathOO{"efficiencyCCDBPathOO", "Users/y/yparovia/LHC25h3", "Path of the efficiency and purity corrections (OO)"};
102+ // * Configurable<std::string> efficiencyCCDBPathpO{"efficiencyCCDBPathpO", "Users/y/yparovia/LHC25h2", "Path of the efficiency and purity corrections (pO)"};
103103
104104 // axes
105105 struct : ConfigurableGroup {
@@ -119,7 +119,7 @@ struct StrangeCascTrack {
119119 Configurable<float > cutZVertex{" cutZVertex" , 10 .0f , " max Z-vertex position" };
120120 Configurable<bool > cutDoSel8{" cutDoSel8" , true , " choose events with sel8" };
121121 Configurable<bool > cutDoNoPileup{" cutDoNoPileup" , true , " choose events with no bunch pileup" };
122- Configurable<bool > cutDoGoodFT0{" cutDoGoodFT0" , true , " choose events with IsGoodZvtxFT0VsPV" };
122+ Configurable<bool > cutDoGoodFT0 {" cutDoGoodFT0" , true , " choose events with IsGoodZvtxFT0VsPV" };
123123 // cascade cuts
124124 Configurable<bool > cutDoPropagateDCA{" cutDoPropagateDCA" , false , " choose events with sel8" };
125125 Configurable<float > cutPropDCAtoPVxy{" cutPropDCAtoPVxy" , 0 .02f , " max cascade dca to PV in xy - propagated" };
@@ -132,7 +132,8 @@ struct StrangeCascTrack {
132132 Configurable<std::vector<float >> cutMinCascCosPaVsPt{
133133 " cutMinCascCosPaVsPt" ,
134134 {0.993 , 0.994 , 0.995 , 0.996 , 0.997 , 0.997 , 0.998 , 0.998 , 0.999 , 0.999 , 0.999 },
135- " Min Casc CosPA per pT bin (same binning as axisPt)" };
135+ " Min Casc CosPA per pT bin (same binning as axisPt)"
136+ };
136137 Configurable<float > cutRapidity{" cutRapidity" , 0 .5f , " max rapidity" };
137138 Configurable<float > cutDauEta{" cutDauEta" , 1 .0f , " max eta of dau tracks" };
138139 Configurable<float > cutCompMassRej{" cutCompMassRej" , 0 .008f , " Competing mass rejection" };
@@ -404,10 +405,11 @@ struct StrangeCascTrack {
404405 bool passedCascCosPA = true ;
405406 double casccospa = cascade.casccosPA (collision.posX (), collision.posY (), collision.posZ ());
406407 const auto & edges = axesConfig.axisPt .value ;
407- int ptBin = std::upper_bound (edges.begin (), edges.end (), cascade.pt ()) - edges.begin () - 1 ;
408+ int ptBin = std::upper_bound (edges.begin (), edges.end (), cascade.pt ())
409+ - edges.begin () - 1 ;
408410 if (ptBin < 0 || ptBin >= (int )selCuts.cutMinCascCosPaVsPt ->size ()) {
409411 ptBin = 1 ;
410- } // safety check - if pt bin not determined, default to loosest cut
412+ } // safety check - if pt bin not determined, default to loosest cut
411413 if (casccospa < selCuts.cutMinCascCosPaVsPt ->at (ptBin)) {
412414 passedCascCosPA = false ;
413415 passedAllSels = false ;
@@ -974,22 +976,22 @@ struct StrangeCascTrack {
974976 histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/GenRecPt" ), genPt, pt);
975977 if (passedAllSelsOmega || passedAllSelsXi) {
976978 histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec/GenRecPt" ), genPt, pt);
977- if (fillTruthXi || fillTruthOmega)
978- histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/GenRecPt" ), genPt, pt);
979+ if (fillTruthXi || fillTruthOmega)
980+ histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/GenRecPt" ), genPt, pt);
979981 }
980982 if (fillTruthXi)
981983 histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/GenRecRapidityXi" ), genYXi, cascade.yXi ());
982984 if (passedAllSelsXi) {
983985 histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec/GenRecRapidityXi" ), genYXi, cascade.yXi ());
984- if (fillTruthXi)
985- histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/GenRecRapidityXi" ), genYXi, cascade.yXi ());
986+ if (fillTruthXi)
987+ histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/GenRecRapidityXi" ), genYXi, cascade.yXi ());
986988 }
987989 if (fillTruthOmega)
988990 histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/GenRecRapidityOmega" ), genYOmega, cascade.yOmega ());
989991 if (passedAllSelsOmega) {
990992 histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec/GenRecRapidityOmega" ), genYOmega, cascade.yOmega ());
991- if (fillTruthOmega)
992- histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/GenRecRapidityOmega" ), genYOmega, cascade.yOmega ());
993+ if (fillTruthOmega)
994+ histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/GenRecRapidityOmega" ), genYOmega, cascade.yOmega ());
993995 }
994996 }
995997 }
@@ -1000,7 +1002,7 @@ struct StrangeCascTrack {
10001002 {
10011003 // check if cut configuration is valid
10021004 if (selCuts.cutMinCascCosPaVsPt ->size () != axesConfig.axisPt .value .size () - 1 ) {
1003- LOGF (fatal, " cutMinCascCosPaVsPt size does not match axisPt binning" );
1005+ LOGF (fatal, " cutMinCascCosPaVsPt size does not match axisPt binning" );
10041006 }
10051007 // for all events processing
10061008 histos.add (" NoSel-Events/EvCounter" , " Event Counter" , kTH1D , {{1 , 0 , 1 }});
@@ -1119,17 +1121,17 @@ struct StrangeCascTrack {
11191121 histos.add (" MC/Gen/EvCounter" , " Event Counter" , kTH1D , {{1 , 0 , 1 }});
11201122 histos.add (" MC/Gen/Mult" , " Event charged multiplicty (all gen events)" , kTH1D , {{3200 , 0 , 3200 }});
11211123 histos.add (" MC/Gen/RapidityXi" , " Generated y (Xi)" , kTH1D , {{200 , -1.0 , 1.0 }});
1122- histos.add (" MC/Gen/Xi" , " Xi" , kTH2D , {{3200 , 0 , 3200 }, {100 , 0.0 , 10.0 }}); // generated Xis
1123- histos.add (" MC/Gen/PrimaryXi" , " Xi primaries" , kTH2D , {{3200 , 0 , 3200 }, {100 , 0.0 , 10.0 }}); // generated primary Xis
1124- histos.add (" MC/Gen/PrimaryXiRapidity" , " Xi primaries in |y|" , kTH2D , {{3200 , 0 , 3200 }, {100 , 0.0 , 10.0 }}); // generated primary Xis in selected rapidity range
1125- histos.add (" MC/Gen/PrimaryXiMinusRapidity" , " Xi- primaries in |y| " , kTH2D , {{3200 , 0 , 3200 }, {100 , 0.0 , 10.0 }}); // generated primary Xi-minus in selected rapidity range
1126- histos.add (" MC/Gen/PrimaryXiPlusRapidity" , " Xi+ primaries in |y|" , kTH2D , {{3200 , 0 , 3200 }, {100 , 0.0 , 10.0 }}); // generated primary Xi-plus in selected rapidity range
1124+ histos.add (" MC/Gen/Xi" , " Xi" , kTH2D , {{3200 , 0 , 3200 }, {100 , 0.0 , 10.0 }}); // generated Xis
1125+ histos.add (" MC/Gen/PrimaryXi" , " Xi primaries" , kTH2D , {{3200 , 0 , 3200 }, {100 , 0.0 , 10.0 }}); // generated primary Xis
1126+ histos.add (" MC/Gen/PrimaryXiRapidity" , " Xi primaries in |y|" , kTH2D , {{3200 , 0 , 3200 }, {100 , 0.0 , 10.0 }}); // generated primary Xis in selected rapidity range
1127+ histos.add (" MC/Gen/PrimaryXiMinusRapidity" , " Xi- primaries in |y| " , kTH2D , {{3200 , 0 , 3200 }, {100 , 0.0 , 10.0 }}); // generated primary Xi-minus in selected rapidity range
1128+ histos.add (" MC/Gen/PrimaryXiPlusRapidity" , " Xi+ primaries in |y|" , kTH2D , {{3200 , 0 , 3200 }, {100 , 0.0 , 10.0 }}); // generated primary Xi-plus in selected rapidity range
11271129 histos.add (" MC/Gen/RapidityOmega" , " Generated y (Omega)" , kTH1D , {{200 , -1.0 , 1.0 }});
1128- histos.add (" MC/Gen/Omega" , " Omega" , kTH2D , {{3200 , 0 , 3200 }, {100 , 0.0 , 10.0 }}); // generated Omegas
1129- histos.add (" MC/Gen/PrimaryOmega" , " Omega primaries" , kTH2D , {{3200 , 0 , 3200 }, {100 , 0.0 , 10.0 }}); // generated primary Omegas
1130- histos.add (" MC/Gen/PrimaryOmegaRapidity" , " Omega primaries in |y|" , kTH2D , {{3200 , 0 , 3200 }, {100 , 0.0 , 10.0 }}); // generated primary Omegas in selected rapidity range
1131- histos.add (" MC/Gen/PrimaryOmegaMinusRapidity" , " Omega- primaries in |y| " , kTH2D , {{3200 , 0 , 3200 }, {100 , 0.0 , 10.0 }}); // generated primary Omega-minus in selected rapidity range
1132- histos.add (" MC/Gen/PrimaryOmegaPlusRapidity" , " Omega+ primaries in |y|" , kTH2D , {{3200 , 0 , 3200 }, {100 , 0.0 , 10.0 }}); // generated primary Omega-plus in selected rapidity range
1130+ histos.add (" MC/Gen/Omega" , " Omega" , kTH2D , {{3200 , 0 , 3200 }, {100 , 0.0 , 10.0 }}); // generated Omegas
1131+ histos.add (" MC/Gen/PrimaryOmega" , " Omega primaries" , kTH2D , {{3200 , 0 , 3200 }, {100 , 0.0 , 10.0 }}); // generated primary Omegas
1132+ histos.add (" MC/Gen/PrimaryOmegaRapidity" , " Omega primaries in |y|" , kTH2D , {{3200 , 0 , 3200 }, {100 , 0.0 , 10.0 }}); // generated primary Omegas in selected rapidity range
1133+ histos.add (" MC/Gen/PrimaryOmegaMinusRapidity" , " Omega- primaries in |y| " , kTH2D , {{3200 , 0 , 3200 }, {100 , 0.0 , 10.0 }}); // generated primary Omega-minus in selected rapidity range
1134+ histos.add (" MC/Gen/PrimaryOmegaPlusRapidity" , " Omega+ primaries in |y|" , kTH2D , {{3200 , 0 , 3200 }, {100 , 0.0 , 10.0 }}); // generated primary Omega-plus in selected rapidity range
11331135 // generated events with reco >= 1:
11341136 histos.add (" MC/EvRec/EvCounter" , " Event Counter" , kTH1D , {{1 , 0 , 1 }});
11351137 histos.add (" MC/EvRec/Mult" , " Event charged multiplicty (gen events with reco>=1)" , kTH1D , {{3200 , 0 , 3200 }});
@@ -1139,11 +1141,11 @@ struct StrangeCascTrack {
11391141 histos.add (" MC/EvRec/PrimaryXiRapidity" , " Xi primaries in |y|" , kTH2D , {axesConfig.axisPt , axesConfig.axisMult }); // generated primary Xis in selected rapidity range in reco>=1
11401142 histos.add (" MC/EvRec/PrimaryXiMinusRapidity" , " Xi- primaries in |y| " , kTH2D , {axesConfig.axisPt , axesConfig.axisMult }); // generated primary Xi-minus in selected rapidity range in reco>=1
11411143 histos.add (" MC/EvRec/PrimaryXiPlusRapidity" , " Xi+ primaries in |y|" , kTH2D , {axesConfig.axisPt , axesConfig.axisMult }); // generated primary Xi-plus in selected rapidity range in reco>=1
1142- histos.add (" MC/EvRec/Omega" , " Omega" , kTH2D , {axesConfig.axisPt , axesConfig.axisMult }); // generated Omegas in reco>=1
1143- histos.add (" MC/EvRec/PrimaryOmega" , " Omega primaries" , kTH2D , {axesConfig.axisPt , axesConfig.axisMult }); // generated primary Omegas in reco>=1
1144- histos.add (" MC/EvRec/PrimaryOmegaRapidity" , " Omega primaries in |y|" , kTH2D , {axesConfig.axisPt , axesConfig.axisMult }); // generated primary Omegas in selected rapidity range in reco>=1
1145- histos.add (" MC/EvRec/PrimaryOmegaMinusRapidity" , " Omega- primaries in |y|" , kTH2D , {axesConfig.axisPt , axesConfig.axisMult }); // generated primary Omega-minus in selected rapidity range in reco>=1
1146- histos.add (" MC/EvRec/PrimaryOmegaPlusRapidity" , " Omega+ primaries in |y|" , kTH2D , {axesConfig.axisPt , axesConfig.axisMult }); // generated primary Omega-plus in selected rapidity range in reco>=1
1144+ histos.add (" MC/EvRec/Omega" , " Omega" , kTH2D , {axesConfig.axisPt , axesConfig.axisMult }); // generated Omegas in reco>=1
1145+ histos.add (" MC/EvRec/PrimaryOmega" , " Omega primaries" , kTH2D , {axesConfig.axisPt , axesConfig.axisMult }); // generated primary Omegas in reco>=1
1146+ histos.add (" MC/EvRec/PrimaryOmegaRapidity" , " Omega primaries in |y|" , kTH2D , {axesConfig.axisPt , axesConfig.axisMult }); // generated primary Omegas in selected rapidity range in reco>=1
1147+ histos.add (" MC/EvRec/PrimaryOmegaMinusRapidity" , " Omega- primaries in |y|" , kTH2D , {axesConfig.axisPt , axesConfig.axisMult }); // generated primary Omega-minus in selected rapidity range in reco>=1
1148+ histos.add (" MC/EvRec/PrimaryOmegaPlusRapidity" , " Omega+ primaries in |y|" , kTH2D , {axesConfig.axisPt , axesConfig.axisMult }); // generated primary Omega-plus in selected rapidity range in reco>=1
11471149 // label filter statistic bins for standard cascs
11481150 histos.get <TH1>(HIST (" Standard/Rec/FiltersXi" ))->GetXaxis ()->SetBinLabel (1 , " p_{T}" );
11491151 histos.get <TH1>(HIST (" Standard/Rec/FiltersXi" ))->GetXaxis ()->SetBinLabel (2 , " gen" );
@@ -1291,35 +1293,35 @@ struct StrangeCascTrack {
12911293 auto slicedGenCascs = genCascs.sliceBy (cascsPerMcCollision, genColl.globalIndex ());
12921294 for (auto const & casc : slicedGenCascs) {
12931295 if (casc.straMCCollisionId () != genColl.index ())
1294- continue ; // safety check for correct slicing
1296+ continue ; // safety check for correct slicing
12951297 // fill generated Xi and Omega (all) for statistics
12961298 double cascPt = casc.ptMC ();
12971299 if (isValidPDG (casc, " Xi" ))
1298- histos.fill (HIST (" MC/Gen/Xi" ), genMult, cascPt);
1300+ histos.fill (HIST (" MC/Gen/Xi" ), genMult, cascPt);
12991301 if (isValidPDG (casc, " Omega" ))
1300- histos.fill (HIST (" MC/Gen/Omega" ), genMult, cascPt);
1302+ histos.fill (HIST (" MC/Gen/Omega" ), genMult, cascPt);
13011303 if (casc.isPhysicalPrimary ()) {
13021304 // fill generated primary Xi and Omega for statistics
1303- if (isValidPDG (casc, " Xi" ))
1304- histos.fill (HIST (" MC/Gen/PrimaryXi" ), genMult, cascPt);
1305- if (isValidPDG (casc, " Omega" ))
1306- histos.fill (HIST (" MC/Gen/PrimaryOmega" ), genMult, cascPt);
1305+ if (isValidPDG (casc, " Xi" ))
1306+ histos.fill (HIST (" MC/Gen/PrimaryXi" ), genMult, cascPt);
1307+ if (isValidPDG (casc, " Omega" ))
1308+ histos.fill (HIST (" MC/Gen/PrimaryOmega" ), genMult, cascPt);
13071309 // fill generated primary Xi within rapidity for corrections
1308- if (isValidPDG (casc, " Xi" ) && std::abs (casc.rapidityMC (0 )) < selCuts.cutRapidity ) {
1309- histos.fill (HIST (" MC/Gen/PrimaryXiRapidity" ), genMult, cascPt);
1310- if (casc.pdgCode () == PDG_t::kXiMinus )
1311- histos.fill (HIST (" MC/Gen/PrimaryXiMinusRapidity" ), genMult, cascPt);
1312- if (casc.pdgCode () == PDG_t::kXiPlusBar )
1313- histos.fill (HIST (" MC/Gen/PrimaryXiPlusRapidity" ), genMult, cascPt);
1314- }
1310+ if (isValidPDG (casc, " Xi" ) && std::abs (casc.rapidityMC (0 )) < selCuts.cutRapidity ) {
1311+ histos.fill (HIST (" MC/Gen/PrimaryXiRapidity" ), genMult, cascPt);
1312+ if (casc.pdgCode () == PDG_t::kXiMinus )
1313+ histos.fill (HIST (" MC/Gen/PrimaryXiMinusRapidity" ), genMult, cascPt);
1314+ if (casc.pdgCode () == PDG_t::kXiPlusBar )
1315+ histos.fill (HIST (" MC/Gen/PrimaryXiPlusRapidity" ), genMult, cascPt);
1316+ }
13151317 // fill generated primary Omega within rapidity for corrections
1316- if (isValidPDG (casc, " Omega" ) && std::abs (casc.rapidityMC (2 )) < selCuts.cutRapidity ) {
1317- histos.fill (HIST (" MC/Gen/PrimaryOmegaRapidity" ), genMult, cascPt);
1318- if (casc.pdgCode () == PDG_t::kOmegaMinus )
1319- histos.fill (HIST (" MC/Gen/PrimaryOmegaMinusRapidity" ), genMult, cascPt);
1320- if (casc.pdgCode () == PDG_t::kOmegaPlusBar )
1321- histos.fill (HIST (" MC/Gen/PrimaryOmegaPlusRapidity" ), genMult, cascPt);
1322- }
1318+ if (isValidPDG (casc, " Omega" ) && std::abs (casc.rapidityMC (2 )) < selCuts.cutRapidity ) {
1319+ histos.fill (HIST (" MC/Gen/PrimaryOmegaRapidity" ), genMult, cascPt);
1320+ if (casc.pdgCode () == PDG_t::kOmegaMinus )
1321+ histos.fill (HIST (" MC/Gen/PrimaryOmegaMinusRapidity" ), genMult, cascPt);
1322+ if (casc.pdgCode () == PDG_t::kOmegaPlusBar )
1323+ histos.fill (HIST (" MC/Gen/PrimaryOmegaPlusRapidity" ), genMult, cascPt);
1324+ }
13231325 }
13241326 // fill gen rapidity for statistics
13251327 if (isValidPDG (casc, " Xi" ))
@@ -1330,13 +1332,13 @@ struct StrangeCascTrack {
13301332 // look at rec. ev. >= 1 for corrections
13311333 auto slicedRecColls = recColls.sliceBy (perMcCollision, genColl.globalIndex ());
13321334 if (slicedRecColls.size () > 0 ) {
1333- histos.fill (HIST (" MC/EvRec/EvCounter" ), 0.5 );
1334- histos.fill (HIST (" MC/EvRec/Mult" ), genMult);
1335+ histos.fill (HIST (" MC/EvRec/EvCounter" ), 0.5 );
1336+ histos.fill (HIST (" MC/EvRec/Mult" ), genMult);
13351337 }
13361338 for (auto const & recColl : slicedRecColls) {
13371339 int64_t genCollId = recColl.straMCCollisionId ();
13381340 if (genColl.index () != genCollId)
1339- continue ; // safety check for correct slicing
1341+ continue ; // safety check for correct slicing
13401342 // fill multiplicity-centrality distribution
13411343 double recoMult = (doProcessIons) ? recColl.centFT0C () : recColl.centFT0M ();
13421344 histos.fill (HIST (" MC/EvRec/MultCent" ), recoMult, genMult);
@@ -1403,4 +1405,4 @@ WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)
14031405 return WorkflowSpec{
14041406 adaptAnalysisTask<StrangeCascTrack>(cfgc),
14051407 };
1406- }
1408+ }
0 commit comments