new plot: n_clusters per track vs pt

Author: iravasen

Diff for: Modules/ITS/include/ITS/ITSTrackTask.h

@@ -71,6 +71,7 @@ class ITSTrackTask : public TaskInterface
 
   std::vector<TObject*> mPublishedObjects;
   TH1D* hNClusters;
+  TH2D* hNClustersPt;
   std::unique_ptr<TH1DRatio> hTrackEta;
   std::unique_ptr<TH1DRatio> hTrackPhi;
   TH1D* hVerticesRof;
@@ -114,6 +115,7 @@ class ITSTrackTask : public TaskInterface
   int nVertices = 0;
   double mChipBins[2125]; // x bins for cos(lambda) plot
   double mCoslBins[25];   // y bins for cos(lambda) plot
+  double ptBins[141];     // pt bins
 
   o2::itsmft::TopologyDictionary* mDict;
 };

Diff for: Modules/ITS/src/ITSTrackTask.cxx

@@ -44,6 +44,7 @@ ITSTrackTask::ITSTrackTask() : TaskInterface()
 ITSTrackTask::~ITSTrackTask() // make_shared objects will be delete automatically
 {
   delete hNClusters;
+  delete hNClustersPt;
   delete hVertexCoordinates;
   delete hVertexRvsZ;
   delete hVertexZ;
@@ -67,6 +68,12 @@ void ITSTrackTask::initialize(o2::framework::InitContext& /*ctx*/)
   mDoNorm = o2::quality_control_modules::common::getFromConfig<int>(mCustomParameters, "doNorm", mDoNorm);
   mInvMasses = o2::quality_control_modules::common::getFromConfig<int>(mCustomParameters, "InvMasses", mInvMasses);
 
+  // pt bins definition: 20 MeV/c width up to 1 GeV/c, 100 MeV/c afterwards
+  ptBins[0] = 0.;
+  for (int i = 1; i < 141; i++) {
+    ptBins[i] = i <= 50 ? ptBins[i - 1] + 0.02 : ptBins[i - 1] + 0.1;
+  }
+
   createAllHistos();
   publishHistos();
 
@@ -200,6 +207,7 @@ void ITSTrackTask::monitorData(o2::framework::ProcessingContext& ctx)
       hTrackPhi->getNum()->Fill(out.getPhi());
       hAngularDistribution->getNum()->Fill(Eta, out.getPhi());
       hNClusters->Fill(track.getNumberOfClusters());
+      hNClustersPt->Fill(track.getPt(), track.getNumberOfClusters());
 
       hTrackPtVsEta->Fill(out.getPt(), Eta);
       hTrackPtVsPhi->Fill(out.getPt(), out.getPhi());
@@ -402,6 +410,7 @@ void ITSTrackTask::reset()
   ILOG(Debug, Devel) << "Resetting the histograms" << ENDM;
   hAngularDistribution->Reset();
   hNClusters->Reset();
+  hNClustersPt->Reset();
   hTrackPhi->Reset();
   hTrackEta->Reset();
   hVerticesRof->Reset();
@@ -450,6 +459,11 @@ void ITSTrackTask::createAllHistos()
   formatAxes(hNClusters, "Number of clusters per Track", "Counts", 1, 1.10);
   hNClusters->SetStats(0);
 
+  hNClustersPt = new TH2D("NClustersPt", "NClustersPt;#it{p}_{T} (GeV/#it{c}); Number of clusters per Track", 140, ptBins, 4, 3.5, 7.5);
+  addObject(hNClustersPt);
+  formatAxes(hNClustersPt, "#it{p}_{T} (GeV/#it{c})", "Number of clusters per Track", 1, 1.10);
+  hNClustersPt->SetStats(0);
+
   hTrackEta = std::make_unique<TH1DRatio>("EtaDistribution", "EtaDistribution", 40, -2.0, 2.0, true);
   if (mDoNorm) {
     hTrackEta->SetBit(TH1::kIsAverage);