Skip to content

Commit c6f21ab

Browse files
tutripatTulika Tripathyalibuild
authored
[PWGLF] New control histograms for normalization in MFT analysis (#11150)
Co-authored-by: Tulika Tripathy <tripathy@clralicemac08.local> Co-authored-by: ALICE Action Bot <alibuild@cern.ch>
1 parent ea60ce6 commit c6f21ab

File tree

1 file changed

+74
-28
lines changed

1 file changed

+74
-28
lines changed

‎PWGLF/Tasks/GlobalEventProperties/dndeta-mft-pp.cxx‎

Lines changed: 74 additions & 28 deletions
Original file line numberDiff line numberDiff line change
@@ -162,6 +162,18 @@ struct PseudorapidityDensityMFT {
162162
registry.add({"EventsNtrkZvtx_gt0",
163163
"; N_{trk}; #it{z}_{vtx} (cm); events",
164164
{HistType::kTH2F, {MultAxis, ZAxis}}});
165+
registry.add({"Tracks/2Danalysis/EventsNtrkZvtx_all",
166+
"; N_{trk}; #it{z}_{vtx} (cm); events",
167+
{HistType::kTH2F, {MultAxis, ZAxis}}});
168+
registry.add({"Tracks/2Danalysis/EventsNtrkZvtx_sel8",
169+
"; N_{trk}; #it{z}_{vtx} (cm); events",
170+
{HistType::kTH2F, {MultAxis, ZAxis}}});
171+
registry.add({"Tracks/2Danalysis/EventsNtrkZvtx_sel8_inelgt0",
172+
"; N_{trk}; #it{z}_{vtx} (cm); events",
173+
{HistType::kTH2F, {MultAxis, ZAxis}}});
174+
registry.add({"Tracks/2Danalysis/EventsNtrkZvtx_sel8_inelfwdgt0",
175+
"; N_{trk}; #it{z}_{vtx} (cm); events",
176+
{HistType::kTH2F, {MultAxis, ZAxis}}});
165177
registry.add({"Tracks/Control/DCAXY",
166178
" ; DCA_{XY} (cm)",
167179
{HistType::kTH1F, {DCAxyAxis}}});
@@ -319,7 +331,18 @@ struct PseudorapidityDensityMFT {
319331
registry.add({"Tracks/Control/woOrp/woOrpEtaZvtx_gt0",
320332
"; #eta; #it{z}_{vtx} (cm); tracks",
321333
{HistType::kTH2F, {EtaAxis, ZAxis}}}); //
322-
334+
registry.add({"Tracks/2Danalysis/EtaZvtx",
335+
"; #eta; #it{z}_{vtx} (cm); tracks",
336+
{HistType::kTH2F, {EtaAxis, ZAxis}}}); //
337+
registry.add({"Tracks/2Danalysis/EtaZvtx_sel8",
338+
"; #eta; #it{z}_{vtx} (cm); tracks",
339+
{HistType::kTH2F, {EtaAxis, ZAxis}}}); //
340+
registry.add({"Tracks/2Danalysis/EtaZvtx_sel8_inelgt0",
341+
"; #eta; #it{z}_{vtx} (cm); tracks",
342+
{HistType::kTH2F, {EtaAxis, ZAxis}}}); //
343+
registry.add({"Tracks/2Danalysis/EtaZvtx_sel8_inelfwdgt0",
344+
"; #eta; #it{z}_{vtx} (cm); tracks",
345+
{HistType::kTH2F, {EtaAxis, ZAxis}}}); //
323346
registry.add({"Tracks/Control/woOrp/woOrpTracksPhiEta",
324347
"; #varphi; #eta; tracks",
325348
{HistType::kTH2F, {PhiAxis, EtaAxis}}}); //
@@ -539,25 +562,34 @@ struct PseudorapidityDensityMFT {
539562
FiCentralTracks const& midtracks, aod::Tracks const&)
540563
{
541564
registry.fill(HIST("EventSelection"), 1.);
542-
if (!useEvSel || (useEvSel && collision.sel8())) {
543-
registry.fill(HIST("EventSelection"), 2.);
544-
auto z = collision.posZ();
545-
if ((z >= cfgVzCut1) && (z <= cfgVzCut2)) {
565+
auto perCollisionSample = sampleCentral->sliceByCached(
566+
o2::aod::track::collisionId, collision.globalIndex(), cache);
567+
auto Ntrk = perCollisionSample.size();
568+
auto z = collision.posZ();
569+
registry.fill(HIST("EventsNtrkZvtx"), Ntrk, z);
570+
if ((z >= cfgVzCut1) && (z <= cfgVzCut2)) {
571+
registry.fill(HIST("Tracks/2Danalysis/EventsNtrkZvtx_all"), Ntrk, z);
572+
for (auto& retrack : retracks) {
573+
auto track = retrack.mfttrack();
574+
if ((cfgnEta1 < track.eta()) && (track.eta() < cfgnEta2) && track.nClusters() >= cfgnCluster && retrack.ambDegree() > 0) {
575+
registry.fill(HIST("Tracks/2Danalysis/EtaZvtx"), track.eta(), z);
576+
}
577+
}
578+
if (!useEvSel || (useEvSel && collision.sel8())) {
579+
registry.fill(HIST("EventSelection"), 2.);
580+
registry.fill(HIST("Tracks/2Danalysis/EventsNtrkZvtx_sel8"), Ntrk, z);
546581
registry.fill(HIST("EventSelection"), 3.);
547-
auto perCollisionSample = sampleCentral->sliceByCached(
548-
o2::aod::track::collisionId, collision.globalIndex(), cache);
549-
auto Ntrk = perCollisionSample.size();
550582
std::unordered_set<int> uniqueEvents;
551583
std::unordered_set<int> uniqueEventsAmb;
552584
std::unordered_set<int> uniqueCollisions;
553585
std::unordered_set<int> uniqueCollisionsAmb;
554586
std::unordered_set<int> eventsInelMFT;
555587
std::unordered_set<int> eventsInel;
556-
registry.fill(HIST("EventsNtrkZvtx"), Ntrk, z);
557588
if (midtracks.size() > 0) {
558589
registry.fill(HIST("EventSelection"), 4.);
559590
registry.fill(HIST("EventSelection"), 10.);
560591
registry.fill(HIST("EventsNtrkZvtx_gt0"), Ntrk, z);
592+
registry.fill(HIST("Tracks/2Danalysis/EventsNtrkZvtx_sel8_inelgt0"), Ntrk, z);
561593
eventsInel.insert(collision.globalIndex());
562594
}
563595
if (perCollisionSample.size() > 0) {
@@ -567,15 +599,28 @@ struct PseudorapidityDensityMFT {
567599
registry.fill(HIST("EventSelection"), 12.);
568600
}
569601
int64_t i = 0.0, j = 0.0, k = 0.0;
602+
for (auto& retrack : retracks) {
603+
auto track = retrack.mfttrack();
604+
if ((cfgnEta1 < track.eta()) && (track.eta() < cfgnEta2) && track.nClusters() >= cfgnCluster && retrack.ambDegree() > 0) {
605+
registry.fill(HIST("Tracks/2Danalysis/EtaZvtx_sel8"), track.eta(), z);
606+
if (midtracks.size() > 0 && retrack.ambDegree() > 0) {
607+
registry.fill(HIST("Tracks/2Danalysis/EtaZvtx_sel8_inelgt0"), track.eta(), z);
608+
}
609+
}
610+
}
570611
if (retracks.size() > 0) {
571612
registry.fill(HIST("EventSelection"), 5.);
613+
if (midtracks.size() > 0) {
614+
registry.fill(HIST("Tracks/2Danalysis/EventsNtrkZvtx_sel8_inelfwdgt0"), Ntrk, z);
615+
}
572616
for (auto& retrack : retracks) {
573617
auto track = retrack.mfttrack();
574618

575619
if ((cfgnEta1 < track.eta()) && (track.eta() < cfgnEta2) && track.nClusters() >= cfgnCluster) {
576620
registry.fill(HIST("TracksEtaZvtx"), track.eta(), z);
577621
if (midtracks.size() > 0 && retrack.ambDegree() > 0) {
578622
registry.fill(HIST("Tracks/EtaZvtx_gt0"), track.eta(), z);
623+
registry.fill(HIST("Tracks/2Danalysis/EtaZvtx_sel8_inelfwdgt0"), track.eta(), z);
579624
eventsInelMFT.insert(retrack.bestCollisionId());
580625
}
581626
if (retrack.ambDegree() != 0) {
@@ -852,29 +897,30 @@ struct PseudorapidityDensityMFT {
852897
if (std::abs(charge) < 3.) {
853898
continue;
854899
}
855-
856-
registry.fill(HIST("TracksEtaZvtxGen_t"), particle.eta(),
857-
mcCollision.posZ());
858-
if (perCollisionMCSampleCentral.size() > 0) {
859-
registry.fill(HIST("TracksEtaZvtxGen_gt0t"), particle.eta(),
900+
if (cfgnEta1 < particle.eta() && particle.eta() < cfgnEta2) {
901+
registry.fill(HIST("TracksEtaZvtxGen_t"), particle.eta(),
860902
mcCollision.posZ());
861-
registry.fill(HIST("TracksPhiEtaGen_gt0t"), particle.phi(), particle.eta());
862-
}
863-
if (atLeastOne) {
864-
registry.fill(HIST("TracksEtaZvtxGen"), particle.eta(),
865-
mcCollision.posZ());
866-
registry.fill(HIST("TracksPtEtaGen"), particle.pt(), particle.eta());
867-
if (atLeastOne_gt0) {
868-
registry.fill(HIST("TracksEtaZvtxGen_gt0"), particle.eta(),
903+
if (perCollisionMCSampleCentral.size() > 0) {
904+
registry.fill(HIST("TracksEtaZvtxGen_gt0t"), particle.eta(),
869905
mcCollision.posZ());
870-
registry.fill(HIST("TracksPhiEtaGen_gt0"), particle.phi(), particle.eta());
906+
registry.fill(HIST("TracksPhiEtaGen_gt0t"), particle.phi(), particle.eta());
907+
}
908+
if (atLeastOne) {
909+
registry.fill(HIST("TracksEtaZvtxGen"), particle.eta(),
910+
mcCollision.posZ());
911+
registry.fill(HIST("TracksPtEtaGen"), particle.pt(), particle.eta());
912+
if (atLeastOne_gt0) {
913+
registry.fill(HIST("TracksEtaZvtxGen_gt0"), particle.eta(),
914+
mcCollision.posZ());
915+
registry.fill(HIST("TracksPhiEtaGen_gt0"), particle.phi(), particle.eta());
916+
}
871917
}
872-
}
873918

874-
registry.fill(HIST("TracksPhiEtaGen"), particle.phi(), particle.eta());
875-
registry.fill(HIST("TracksPhiZvtxGen"), particle.phi(),
876-
mcCollision.posZ());
877-
registry.fill(HIST("TracksPtEtaGen_t"), particle.pt(), particle.eta());
919+
registry.fill(HIST("TracksPhiEtaGen"), particle.phi(), particle.eta());
920+
registry.fill(HIST("TracksPhiZvtxGen"), particle.phi(),
921+
mcCollision.posZ());
922+
registry.fill(HIST("TracksPtEtaGen_t"), particle.pt(), particle.eta());
923+
}
878924
}
879925
}
880926
}

0 commit comments

Comments
 (0)