ITS: cut for real primary particles for efficiency calculation (#2125)

* cut for real primary particles for efficiency calculation * add more explicative plot titles
AliceO2Group · Feb 8, 2024 · 3ed2699 · 3ed2699
1 parent 02dc412
commit 3ed2699
Showing 1 changed file with 19 additions and 16 deletions.
diff --git a/Modules/ITS/src/ITSTrackSimTask.cxx b/Modules/ITS/src/ITSTrackSimTask.cxx
@@ -205,13 +205,14 @@ void ITSTrackSimTask::monitorData(o2::framework::ProcessingContext& ctx)
       info[i][mc].isPrimary = mcTrack.isPrimary();
       if (mcTrack.isPrimary()) {
         hPrimaryGen_pt->Fill(mcTrack.GetPt());
-      }
 
-      hDenTrue_r->Fill(distance);
-      hDenTrue_pt->Fill(mcTrack.GetPt());
-      hDenTrue_eta->Fill(mcTrack.GetEta());
-      hDenTrue_phi->Fill(mcTrack.GetPhi());
-      hDenTrue_z->Fill(mcTrack.Vz());
+        // True Generated primaries: denominator of the efficiency plots
+        hDenTrue_r->Fill(distance);
+        hDenTrue_pt->Fill(mcTrack.GetPt());
+        hDenTrue_eta->Fill(mcTrack.GetEta());
+        hDenTrue_phi->Fill(mcTrack.GetPhi());
+        hDenTrue_z->Fill(mcTrack.Vz());
+      }
     }
   }
 
@@ -253,12 +254,14 @@ void ITSTrackSimTask::monitorData(o2::framework::ProcessingContext& ctx)
           hPrimaryReco_pt->Fill(info[MCinfo.getEventID()][MCinfo.getTrackID()].pt);
         }
       } else {
-        hNumRecoValid_pt->Fill(info[MCinfo.getEventID()][MCinfo.getTrackID()].pt);
-        hNumRecoValid_phi->Fill(info[MCinfo.getEventID()][MCinfo.getTrackID()].phi);
-        hNumRecoValid_eta->Fill(info[MCinfo.getEventID()][MCinfo.getTrackID()].eta);
-        hNumRecoValid_z->Fill(info[MCinfo.getEventID()][MCinfo.getTrackID()].z);
-        hNumRecoValid_r->Fill(info[MCinfo.getEventID()][MCinfo.getTrackID()].r);
         if (info[MCinfo.getEventID()][MCinfo.getTrackID()].isPrimary) {
+          // True primaries reconstructed: numerator of efficiency plots
+          hNumRecoValid_pt->Fill(info[MCinfo.getEventID()][MCinfo.getTrackID()].pt);
+          hNumRecoValid_phi->Fill(info[MCinfo.getEventID()][MCinfo.getTrackID()].phi);
+          hNumRecoValid_eta->Fill(info[MCinfo.getEventID()][MCinfo.getTrackID()].eta);
+          hNumRecoValid_z->Fill(info[MCinfo.getEventID()][MCinfo.getTrackID()].z);
+          hNumRecoValid_r->Fill(info[MCinfo.getEventID()][MCinfo.getTrackID()].r);
+
           hTrackImpactTransvValid->Fill(ip[0]);
           hPrimaryReco_pt->Fill(info[MCinfo.getEventID()][MCinfo.getTrackID()].pt);
         }
@@ -357,7 +360,7 @@ void ITSTrackSimTask::createAllHistos()
 
   hDuplicate_pt = new TEfficiency("Duplicate_pt", "#it{p}_{T} fraction of mother duplicate track;#it{p}_{T} (GeV/#it{c});Fraction", nb, xbins);
   addObject(hDuplicate_pt);
-  hEfficiency_pt = new TEfficiency("efficiency_pt", "#it{p}_{T} efficiency of tracking; #it{p}_{T} (GeV/#it{c}); {p}_{T}Efficiency", nb, xbins);
+  hEfficiency_pt = new TEfficiency("efficiency_pt", "Primary pions with 7cls - tracking efficiency vs #it{p}_{T}; #it{p}_{T} (GeV/#it{c}); {p}_{T}Efficiency", nb, xbins);
   addObject(hEfficiency_pt);
   hFakeTrack_pt = new TEfficiency("faketrack_pt", "#it{p}_{T} fake-track rate;#it{p}_{T} (GeV/#it{c});Fake-track rate", nb, xbins);
   addObject(hFakeTrack_pt);
@@ -366,7 +369,7 @@ void ITSTrackSimTask::createAllHistos()
   hNumRecoFake_pt = new TH1D("NumRecoFake_pt", "", nb, xbins);
   hDenTrue_pt = new TH1D("DenTrueMC_pt", "", nb, xbins);
 
-  hEfficiency_phi = new TEfficiency("efficiency_phi", "#phi efficiency of tracking;#phi;Efficiency", 60, 0, TMath::TwoPi());
+  hEfficiency_phi = new TEfficiency("efficiency_phi", "Primary pions with 7cls - tracking efficiency vs #phi;#phi;Efficiency", 60, 0, TMath::TwoPi());
   addObject(hEfficiency_phi);
   hFakeTrack_phi = new TEfficiency("faketrack_phi", "#phi fake-track rate;#phi;Fake-track rate", 60, 0, TMath::TwoPi());
   addObject(hFakeTrack_phi);
@@ -377,7 +380,7 @@ void ITSTrackSimTask::createAllHistos()
   hNumRecoFake_phi = new TH1D("NumRecoFake_phi", "", 60, 0, TMath::TwoPi());
   hDenTrue_phi = new TH1D("DenTrueMC_phi", "", 60, 0, TMath::TwoPi());
 
-  hEfficiency_eta = new TEfficiency("efficiency_eta", "#eta efficiency of tracking;#eta;Efficiency", 30, -1.5, 1.5);
+  hEfficiency_eta = new TEfficiency("efficiency_eta", "Primary pions with 7 cls - tracking efficiency vs #eta;#eta;Efficiency", 30, -1.5, 1.5);
   addObject(hEfficiency_eta);
   hFakeTrack_eta = new TEfficiency("faketrack_eta", "#eta fake-track rate;#eta;Fake-track rate", 30, -1.5, 1.5);
   addObject(hFakeTrack_eta);
@@ -388,7 +391,7 @@ void ITSTrackSimTask::createAllHistos()
   hNumRecoFake_eta = new TH1D("NumRecoFake_eta", "", 30, -1.5, 1.5);
   hDenTrue_eta = new TH1D("DenTrueMC_eta", "", 30, -1.5, 1.5);
 
-  hEfficiency_r = new TEfficiency("efficiency_r", "r efficiency of tracking;r (cm);Efficiency", 100, 0, 5);
+  hEfficiency_r = new TEfficiency("efficiency_r", "Primary pions with 7 cls - tracking efficiency vs r;r (cm);Efficiency", 100, 0, 5);
   addObject(hEfficiency_r);
   hFakeTrack_r = new TEfficiency("faketrack_r", "r fake-track rate;r (cm);Fake-track rate", 100, 0, 5);
   addObject(hFakeTrack_r);
@@ -399,7 +402,7 @@ void ITSTrackSimTask::createAllHistos()
   hNumDuplicate_r = new TH1D("NumDuplicate_r", "", 100, 0, 5);
   hDenTrue_r = new TH1D("DenTrueMC_r", "", 100, 0, 5);
 
-  hEfficiency_z = new TEfficiency("efficiency_z", "z efficiency of tracking;z (cm);Efficiency", 101, -5, 5);
+  hEfficiency_z = new TEfficiency("efficiency_z", "Primary pions with 7 cls - tracking efficiency vs z;z (cm);Efficiency", 101, -5, 5);
   addObject(hEfficiency_z);
   hFakeTrack_z = new TEfficiency("faketrack_z", "z fake-track rate;z (cm);Fake-track rate", 101, -5, 5);
   addObject(hFakeTrack_z);