adding stat histos for lumi computation

Author: rolavick

PWGUD/TableProducer/twoTracksEventTableProducer.cxx

@@ -116,7 +116,7 @@ struct TwoTracksEventTableProducer {
     Configurable<bool> preselBothAreTOFtracks{"preselBothAreTOFtracks", false, {"Both tracks are required to hit TOF."}};
     Configurable<bool> preselUseMinMomentumOnBothTracks{"preselUseMinMomentumOnBothTracks", false, {"Both tracks are required to fill requirement on minimum momentum."}};
     Configurable<float> preselMinTrackMomentum{"preselMinTrackMomentum", 0.1, {"Requirement on minimum momentum of the track."}};
-    Configurable<float> preselSystemPtCut{"preselSystemPtCut", 0.0, {"By default, cut on maximum system pT."}};
+    Configurable<float> preselSystemPtCut{"preselSystemPtCut", 2.0, {"By default, cut on maximum system pT."}};
     Configurable<bool> preselUseOppositeSystemPtCut{"preselUseOppositeSystemPtCut", false, {"Negates the system pT cut (cut on minimum system pT)."}};
     Configurable<float> preselMinInvariantMass{"preselMinInvariantMass", 2.0, {"Requirement on minimum system invariant mass."}};
     Configurable<float> preselMaxInvariantMass{"preselMaxInvariantMass", 5.0, {"Requirement on maximum system invariant mass."}};
@@ -137,8 +137,10 @@ struct TwoTracksEventTableProducer {
 
     mySetITShitsRule(cutGlobalTrack.cutITShitsRule);
 
-    histos.add("Reco/hSelections", "Effect of selections;;Number of  events (-)", HistType::kTH1D, {{50, 0.5, 50.5}});
-    histos.add("Truth/hTroubles", "Counter of unwanted issues;;Number of  troubles (-)", HistType::kTH1D, {{15, 0.5, 15.5}});
+    histos.add("Reco/hSelections", "Effect of selections;;Number of events (-)", HistType::kTH1D, {{50, 0.5, 50.5}});
+	  histos.add("Reco/hNanalyzedPerRun", "N analyzed events per run;Run number (-);Number of analyzed events (-)", HistType::kTH1D, {{1, 0., 1.}});
+	  histos.add("Reco/hNselectedPerRun", "N selected events per run;Run number (-);Number of selected events (-)", HistType::kTH1D, {{1, 0., 1.}});
+    histos.add("Truth/hTroubles", "Counter of unwanted issues;;Number of troubles (-)", HistType::kTH1D, {{15, 0.5, 15.5}});
 
   } // end init
 
@@ -358,6 +360,9 @@ struct TwoTracksEventTableProducer {
   void processDataSG(FullSGUDCollision const& collision,
                      FullUDTracks const& tracks)
   {
+	  int run = collision.runNumber();
+	  const char* srun = Form("%d", run);
+	  histos.get<TH1>(HIST("Reco/hNanalyzedPerRun"))->Fill(srun, 1);
 
     nSelection = 0;
     histos.get<TH1>(HIST("Reco/hSelections"))->Fill(nSelection);
@@ -490,6 +495,8 @@ struct TwoTracksEventTableProducer {
     float tofEP[2] = {trk1.tofExpMom(), trk2.tofExpMom()};
     float infoZDC[4] = {collision.energyCommonZNA(), collision.energyCommonZNC(), collision.timeZNA(), collision.timeZNC()};
 
+	  histos.get<TH1>(HIST("Reco/hNselectedPerRun"))->Fill(srun, 1);
+
     twoTracks(collision.runNumber(), collision.globalBC(), countTracksPerCollision, collision.numContrib(), countGoodNonPVtracks, collision.posX(), collision.posY(), collision.posZ(),
               collision.flags(), collision.occupancyInTime(), collision.hadronicRate(), collision.trs(), collision.trofs(), collision.hmpr(),
               collision.tfb(), collision.itsROFb(), collision.sbp(), collision.zVtxFT0vPV(), collision.vtxITSTPC(),
@@ -500,7 +507,7 @@ struct TwoTracksEventTableProducer {
               tpcSignal, tpcEl, tpcMu, tpcPi, tpcKa, tpcPr, tpcIP,
               tofSignal, tofEl, tofMu, tofPi, tofKa, tofPr, tofEP);
   }
-  PROCESS_SWITCH(TwoTracksEventTableProducer, processDataSG, "Iterate UD tables with measured data created by SG-Candidate-Producer.", false);
+  PROCESS_SWITCH(TwoTracksEventTableProducer, processDataSG, "Iterate UD tables with measured data created by SG-Candidate-Producer.", true);
 
   PresliceUnsorted<aod::UDMcParticles> partPerMcCollision = aod::udmcparticle::udMcCollisionId;
   PresliceUnsorted<FullMCSGUDCollisions> colPerMcCollision = aod::udcollision::udMcCollisionId;