@@ -121,10 +121,10 @@ struct NPCascCandidate {
121121 float bachPionTOFNSigma;
122122 bool sel8;
123123 float multFT0C;
124- float multFT0A ;
124+ float multFV0A ;
125125 float multFT0M;
126126 float centFT0C;
127- float centFT0A ;
127+ float centFV0A ;
128128 float centFT0M;
129129 int multNTracksGlobal;
130130 uint32_t toiMask;
@@ -177,8 +177,8 @@ struct NonPromptCascadeTask {
177177
178178 using TracksExtData = soa::Join<aod::TracksIU, aod::TracksCovIU, aod::TracksExtra, aod::pidTPCFullKa, aod::pidTPCFullPi, aod::pidTPCFullPr, aod::pidTOFFullKa, aod::pidTOFFullPi, aod::pidTOFFullPr>;
179179 using TracksExtMC = soa::Join<aod::TracksIU, aod::TracksCovIU, aod::TracksExtra, aod::McTrackLabels, aod::pidTPCFullKa, aod::pidTPCFullPi, aod::pidTPCFullPr, aod::pidTOFFullKa, aod::pidTOFFullPi, aod::pidTOFFullPr>;
180- using CollisionCandidatesRun3 = soa::Join<aod::Collisions, aod::EvSels, aod::FT0Mults, aod::CentFT0Cs, aod::CentFT0As , aod::CentFT0Ms, aod::MultsGlobal>;
181- using CollisionCandidatesRun3MC = soa::Join<aod::Collisions, aod::McCollisionLabels, aod::EvSels, aod::FT0Mults, aod::CentFT0Cs, aod::CentFT0As , aod::CentFT0Ms, aod::MultsGlobal>;
180+ using CollisionCandidatesRun3 = soa::Join<aod::Collisions, aod::EvSels, aod::FT0Mults, aod::FV0Mults, aod:: CentFT0Cs, aod::CentFV0As , aod::CentFT0Ms, aod::MultsGlobal>;
181+ using CollisionCandidatesRun3MC = soa::Join<aod::Collisions, aod::McCollisionLabels, aod::EvSels, aod::FT0Mults, aod::FV0Mults, aod:: CentFT0Cs, aod::CentFV0As , aod::CentFT0Ms, aod::MultsGlobal>;
182182
183183 Preslice<TracksExtData> perCollision = aod::track::collisionId;
184184 Preslice<TracksExtMC> perCollisionMC = aod::track::collisionId;
@@ -205,6 +205,7 @@ struct NonPromptCascadeTask {
205205 Configurable<std::string> cfgTriggersOfInterest{" cfgTriggersOfInterest" , " fTrackedOmega,fOmegaHighMult" , " Triggers of interest, comma separated for Zorro" };
206206
207207 Configurable<float > cfgMaxMult{" cfgMaxMult" , 8000 .f , " Upper range of multiplicty histo" };
208+ Configurable<float > cfgMaxMultFV0{" cfgMaxMultFV0" , 10000 .f , " Upper range of multiplicty FV0 histo" };
208209 Configurable<float > cfgMinMult{" cfgMinMult" , 3000 .f , " Lower range of FT0M histo in zoomed histo" };
209210 Configurable<float > cfgMaxCent{" cfgMaxCent" , 8 .0025f , " Upper range of FT0M histo" };
210211
@@ -224,9 +225,11 @@ struct NonPromptCascadeTask {
224225 std::unordered_map<std::string, std::shared_ptr<TH2>> mHistsPerRunNtracktVsCentZoom ;
225226
226227 int nBinsMult = cfgMaxMult;
228+ int nBinsMultFV0 = cfgMaxMultFV0;
227229 int nBinsMultZoom = cfgMaxMult - cfgMinMult;
228230 int nBinsCentZoom = (cfgMaxCent + 0.0025 ) / 0.005 ;
229231 AxisSpec multAxis = {nBinsMult, 0 , cfgMaxMult, " Multiplicity FT0M" };
232+ AxisSpec multAxisFV0 = {nBinsMultFV0, 0 , cfgMaxMultFV0, " Multiplicity FT0M" };
230233 AxisSpec centAxis = {101 , -0.025 , 101.025 , " Centrality" };
231234 AxisSpec centAxisZoom = {nBinsCentZoom, -0.0025 , cfgMaxCent, " Centrality" };
232235 AxisSpec multAxisZoom = {nBinsMultZoom, cfgMinMult, cfgMaxMult, " Multiplicity FT0M" };
@@ -274,6 +277,12 @@ struct NonPromptCascadeTask {
274277 auto cutsXi{std::get<std::shared_ptr<TH2>>(mRegistry .add (" h_PIDcutsXi" , " ;;Invariant mass (GeV/#it{c}^{2})" , HistType::kTH2D , {{6 , -0.5 , 5.5 }, {125 , 1.296 , 1.346 }}))};
275278 mRegistry .add (" hMultVsCent" , " hMultVsCent" , HistType::kTH2F , {centAxis, multAxis});
276279 mRegistry .add (" hMultVsCentZoom" , " hMultVsCentZoom" , HistType::kTH2F , {centAxisZoom, multAxisZoom});
280+ mRegistry .add (" hMultFV0VshNTracks" , " hMultFV0VshNTracks" , HistType::kTH2F , {nTracksAxis, multAxisFV0});
281+ mRegistry .add (" hNTracksVsCentFV0A" , " hNTracksVsCentFV0A" , HistType::kTH2F , {nTracksAxis, centAxis});
282+ mRegistry .add (" hMultFT0CVshNTracks" , " hMultFT0CVshNTracks" , HistType::kTH2F , {nTracksAxis, multAxis});
283+ mRegistry .add (" hNTracksVsCentFT0C" , " hNTracksVsCentFT0C" , HistType::kTH2F , {nTracksAxis, centAxis});
284+
285+
277286 mRegistry .add (" hNTracksVsCent" , " hNTracksVsCent" , HistType::kTH2F , {centAxis, nTracksAxis});
278287 mRegistry .add (" hNTracksVsCentZoom" , " hNTracksVsCentZoom" , HistType::kTH2F , {centAxisZoom, nTracksAxis});
279288
@@ -362,9 +371,13 @@ struct NonPromptCascadeTask {
362371 mHistsPerRunNtracktVsCentZoom [histNameTvCZ]->Fill (coll.centFT0M (), coll.multNTracksGlobal ());
363372 // run integrated histos
364373 mRegistry .fill (HIST (" hMultVsCent" ), coll.centFT0M (), coll.multFT0M ());
374+ mRegistry .fill (HIST (" hMultFV0VshNTracks" ), coll.multNTracksGlobal (), coll.multFV0A ());
375+ mRegistry .fill (HIST (" hNTracksVsCentFV0A" ), coll.multNTracksGlobal (), coll.centFV0A ());
376+ mRegistry .fill (HIST (" hMultFT0CVshNTracks" ), coll.multNTracksGlobal (), coll.multFT0C ());
377+ mRegistry .fill (HIST (" hNTracksVsCentFT0C" ), coll.multNTracksGlobal (), coll.centFT0C ());
365378 mRegistry .fill (HIST (" hMultVsCentZoom" ), coll.centFT0M (), coll.multFT0M ());
366379 mRegistry .fill (HIST (" hNTracksVsCent" ), coll.centFT0M (), (float )coll.multNTracksGlobal ());
367- mRegistry .fill (HIST (" hNTracksVsCentZoom" ), coll.centFT0M (), coll.multNTracksGlobal ());
380+ mRegistry .fill (HIST (" hNTracksVsCentZoom" ), coll.centFT0M (), ( float ) coll.multNTracksGlobal ());
368381 }
369382 };
370383
@@ -587,7 +600,7 @@ struct NonPromptCascadeTask {
587600 cascITSclusters, protonTrack.itsNCls (), pionTrack.itsNCls (), bachelor.itsNCls (), protonTrack.tpcNClsFound (), pionTrack.tpcNClsFound (), bachelor.tpcNClsFound (),
588601 protonTrack.tpcNSigmaPr (), pionTrack.tpcNSigmaPi (), bachelor.tpcNSigmaKa (), bachelor.tpcNSigmaPi (),
589602 protonTrack.hasTOF (), pionTrack.hasTOF (), bachelor.hasTOF (),
590- protonTrack.tofNSigmaPr (), pionTrack.tofNSigmaPi (), bachelor.tofNSigmaKa (), bachelor.tofNSigmaPi (), collision.sel8 (), collision.multFT0C (), collision.multFT0A (), collision.multFT0M (), collision.centFT0C (), collision.centFT0A (), collision.centFT0M (), collision.multNTracksGlobal (), toiMask, collision.selection_bit (aod::evsel::kNoSameBunchPileup )});
603+ protonTrack.tofNSigmaPr (), pionTrack.tofNSigmaPi (), bachelor.tofNSigmaKa (), bachelor.tofNSigmaPi (), collision.sel8 (), collision.multFT0C (), collision.multFV0A (), collision.multFT0M (), collision.centFT0C (), collision.centFV0A (), collision.centFT0M (), collision.multNTracksGlobal (), toiMask, collision.selection_bit (aod::evsel::kNoSameBunchPileup )});
591604 }
592605 }
593606
@@ -607,7 +620,7 @@ struct NonPromptCascadeTask {
607620 c.protonTPCNSigma , c.pionTPCNSigma , c.bachKaonTPCNSigma , c.bachPionTPCNSigma ,
608621 c.protonHasTOF , c.pionHasTOF , c.bachHasTOF ,
609622 c.protonTOFNSigma , c.pionTOFNSigma , c.bachKaonTOFNSigma , c.bachPionTOFNSigma ,
610- c.sel8 , c.multFT0C , c.multFT0A , c.multFT0M , c.centFT0C , c.centFT0A , c.centFT0M , c.multNTracksGlobal , c.toiMask , c.noSameBunchPileup );
623+ c.sel8 , c.multFT0C , c.multFV0A , c.multFT0M , c.centFT0C , c.centFV0A , c.centFT0M , c.multNTracksGlobal , c.toiMask , c.noSameBunchPileup );
611624 }
612625 }
613626
@@ -643,7 +656,7 @@ struct NonPromptCascadeTask {
643656 c.cascNClusITS , c.protonNClusITS , c.pionNClusITS , c.bachNClusITS , c.protonNClusTPC , c.pionNClusTPC , c.bachNClusTPC , c.protonTPCNSigma ,
644657 c.pionTPCNSigma , c.bachKaonTPCNSigma , c.bachPionTPCNSigma , c.protonHasTOF , c.pionHasTOF , c.bachHasTOF ,
645658 c.protonTOFNSigma , c.pionTOFNSigma , c.bachKaonTOFNSigma , c.bachPionTOFNSigma ,
646- c.sel8 , c.multFT0C , c.multFT0A , c.multFT0M , c.centFT0C , c.centFT0A , c.centFT0M ,
659+ c.sel8 , c.multFT0C , c.multFV0A , c.multFT0M , c.centFT0C , c.centFV0A , c.centFT0M ,
647660 particle.pt (), particle.eta (), particle.phi (), mcCollision.posX (), mcCollision.posY (), mcCollision.posZ (),
648661 particle.pdgCode (), mcCollision.posX () - particle.vx (), mcCollision.posY () - particle.vy (),
649662 mcCollision.posZ () - particle.vz (), mcCollision.globalIndex () == recCollision.mcCollisionId (), c.hasFakeReassociation , motherDecayDaughters, c.multNTracksGlobal , c.toiMask , c.noSameBunchPileup );
0 commit comments