PWGEM/PhotonMeson: update createPCM and histograms (#3255)

Author: dsekihat

PWGEM/PhotonMeson/Core/HistogramsLibrary.cxx

@@ -49,6 +49,8 @@ void o2::aod::emphotonhistograms::DefineHistograms(THashList* list, const char*
     list->Add(new TH1F("hQoverPt", "q/pT;q/p_{T} (GeV/c)^{-1}", 1000, -50, 50));
     list->Add(new TH2F("hEtaPhi", "#eta vs. #varphi;#varphi (rad.);#eta", 180, 0, TMath::TwoPi(), 40, -2.0f, 2.0f));
     list->Add(new TH2F("hDCAxyz", "DCA xy vs. z;DCA_{xy} (cm);DCA_{z} (cm)", 1000, -50.0f, 50.0f, 1000, -50.0f, 50.0f));
+    list->Add(new TH2F("hNclsTPC_Pt", "number of TPC clusters vs. pT", 1000, 0, 10, 161, -0.5, 160.5));
+    list->Add(new TH2F("hNcrTPC_Pt", "number of TPC crossed rows vs. pT", 1000, 0, 10, 161, -0.5, 160.5));
     list->Add(new TH1F("hNclsTPC", "number of TPC clusters", 161, -0.5, 160.5));
     list->Add(new TH1F("hNcrTPC", "number of TPC crossed rows", 161, -0.5, 160.5));
     list->Add(new TH1F("hChi2TPC", "chi2/number of TPC clusters", 100, 0, 10));

PWGEM/PhotonMeson/Core/HistogramsLibrary.h

@@ -100,8 +100,10 @@ void FillHistClass(THashList* list, const char* subGroup, T const& obj)
     reinterpret_cast<TH1F*>(list->FindObject("hQoverPt"))->Fill(obj.sign() / obj.pt());
     reinterpret_cast<TH2F*>(list->FindObject("hEtaPhi"))->Fill(obj.phi(), obj.eta());
     reinterpret_cast<TH2F*>(list->FindObject("hDCAxyz"))->Fill(obj.dcaXY(), obj.dcaZ());
-    reinterpret_cast<TH1F*>(list->FindObject("hNclsTPC"))->Fill(obj.tpcNClsFound());
+    reinterpret_cast<TH2F*>(list->FindObject("hNclsTPC_Pt"))->Fill(obj.pt(), obj.tpcNClsFound());
+    reinterpret_cast<TH2F*>(list->FindObject("hNcrTPC_Pt"))->Fill(obj.pt(), obj.tpcNClsCrossedRows());
     reinterpret_cast<TH1F*>(list->FindObject("hNclsITS"))->Fill(obj.itsNCls());
+    reinterpret_cast<TH1F*>(list->FindObject("hNclsTPC"))->Fill(obj.tpcNClsFound());
     reinterpret_cast<TH1F*>(list->FindObject("hNcrTPC"))->Fill(obj.tpcNClsCrossedRows());
     reinterpret_cast<TH1F*>(list->FindObject("hTPCNcr2Nf"))->Fill(obj.tpcCrossedRowsOverFindableCls());
     reinterpret_cast<TH1F*>(list->FindObject("hTPCNcls2Nf"))->Fill(obj.tpcFoundOverFindableCls());

PWGEM/PhotonMeson/TableProducer/createPCM.cxx

@@ -83,8 +83,8 @@ struct createPCM {
   Configurable<float> maxeta{"maxeta", 0.9, "eta acceptance for single track"};
   Configurable<int> mincrossedrows{"mincrossedrows", 10, "min crossed rows"};
   Configurable<float> maxchi2tpc{"maxchi2tpc", 4.0, "max chi2/NclsTPC"};
-  Configurable<float> maxchi2its{"maxchi2its", 36.0, "max chi2/NclsITS"};
-  Configurable<float> maxpt_itsonly{"maxpt_itsonly", 1.0, "max pT for ITSonly tracks"};
+  Configurable<float> maxchi2its{"maxchi2its", 5.0, "max chi2/NclsITS"};
+  Configurable<float> maxpt_itsonly{"maxpt_itsonly", 0.5, "max pT for ITSonly tracks"};
   Configurable<float> min_tpcdEdx{"min_tpcdEdx", 30.0, "min TPC dE/dx"};
   Configurable<float> max_tpcdEdx{"max_tpcdEdx", 110.0, "max TPC dE/dx"};
   Configurable<float> margin_r{"margin_r", 7.0, "margin for r cut"};
@@ -330,15 +330,18 @@ struct createPCM {
       }
     }
 
-    bool isITSonly = track.hasITS() & !track.hasTPC() & !track.hasTRD() & !track.hasTOF();
-    if (isITSonly) {
+    if (track.hasITS()) {
       if (track.itsChi2NCl() > maxchi2its) {
         return false;
       }
-      if (track.pt() > maxpt_itsonly) {
-        return false;
+      bool isITSonly = track.hasITS() & !track.hasTPC() & !track.hasTRD() & !track.hasTOF();
+      if (isITSonly) {
+        if (track.pt() > maxpt_itsonly) {
+          return false;
+        }
       }
     }
+
     return true;
   }