Some extra debug functionality

Author: ddobrigk

Common/DataModel/Multiplicity.h

@@ -75,6 +75,12 @@ DECLARE_SOA_COLUMN(MultNTracksGlobal, multNTracksGlobal, int);       //!
 
 DECLARE_SOA_COLUMN(BCNumber, bcNumber, int); //!
 
+// even further QA: timing information for neighboring events
+DECLARE_SOA_COLUMN(TimeToPrePrevious, timeToPrePrevious, float); //!
+DECLARE_SOA_COLUMN(TimeToPrevious, timeToPrevious, float); //!
+DECLARE_SOA_COLUMN(TimeToNext, timeToNext, float); //!
+DECLARE_SOA_COLUMN(TimeToNeNext, timeToNeNext, float); //!
+
 } // namespace mult
 DECLARE_SOA_TABLE(FV0Mults, "AOD", "FV0MULT", //! Multiplicity with the FV0 detector
                   mult::MultFV0A, mult::MultFV0C,
@@ -109,6 +115,10 @@ DECLARE_SOA_TABLE(MultsExtra, "AOD", "MULTEXTRA", //!
                   mult::MultAllTracksTPCOnly, mult::MultAllTracksITSTPC,
                   mult::BCNumber);
 
+DECLARE_SOA_TABLE(MultNeighs, "AOD", "MULTNEIGH", //!
+                  mult::TimeToPrePrevious, mult::TimeToPrevious, 
+                  mult::TimeToNext, mult::TimeToNeNext);
+
 // for QA purposes
 DECLARE_SOA_TABLE(MultsGlobal, "AOD", "MULTGLOBAL", //! counters that use Track Selection (optional)
                   mult::MultNTracksGlobal);

Common/TableProducer/multiplicityExtraTable.cxx

@@ -32,6 +32,7 @@ const int nBCsPerOrbit = o2::constants::lhc::LHCMaxBunches;
 
 struct MultiplicityExtraTable {
   Produces<aod::MultsBC> multBC;
+  Produces<aod::MultNeighs> multNeigh;
 
   o2::ccdb::CcdbApi ccdbApi;
   Service<o2::ccdb::BasicCCDBManager> ccdb;
@@ -50,7 +51,7 @@ struct MultiplicityExtraTable {
 
   using BCsWithRun3Matchings = soa::Join<aod::BCs, aod::Timestamps, aod::Run3MatchedToBCSparse>;
 
-  void process(BCsWithRun3Matchings::iterator const& bc, aod::FV0As const&, aod::FT0s const&, aod::FDDs const&, aod::Zdcs const&)
+  void processBCs(BCsWithRun3Matchings::iterator const& bc, aod::FV0As const&, aod::FT0s const&, aod::FDDs const&, aod::Zdcs const&)
   {
     bool Tvx = false;
     bool isFV0OrA = false;
@@ -153,6 +154,33 @@ struct MultiplicityExtraTable {
 
     multBC(multFT0A, multFT0C, multFV0A, multFDDA, multFDDC, multZNA, multZNC, multZEM1, multZEM2, multZPA, multZPC, Tvx, isFV0OrA, multFV0TriggerBits, multFT0TriggerBits, multFDDTriggerBits, multBCTriggerMask, collidingBC);
   }
+
+  void processCollisionNeighbors(aod::Collisions const& collisions)
+  {
+    std::vector<float> timeArray;
+    timeArray.resize(collisions.size(), 1e+3);
+
+    for (const auto& collision : collisions) {
+      timeArray[collision.globalIndex()] = collision.collisionTime(); 
+    }
+
+    float deltaPrevious = 1e+6, deltaPrePrevious = 1e+6; 
+    float deltaNext = 1e+6, deltaNeNext = 1e+6; 
+    for (const auto& collision : collisions) {
+      int ii = collision.globalIndex();
+
+      if(ii-1>=0) deltaPrevious = timeArray[ii]-timeArray[ii-1];
+      if(ii-2>=0) deltaPrePrevious = timeArray[ii]-timeArray[ii-2];
+      if(ii+1<collisions.size()) deltaNext = timeArray[ii+1]-timeArray[ii];
+      if(ii+2<collisions.size()) deltaNeNext = timeArray[ii+2]-timeArray[ii];
+
+      multNeigh(deltaPrePrevious, deltaPrevious, deltaNext, deltaNeNext);
+    }
+  }
+
+  // Process switches
+  PROCESS_SWITCH(MultiplicityExtraTable, processBCs, "Produce BC tables", true);
+  PROCESS_SWITCH(MultiplicityExtraTable, processCollisionNeighbors, "Produce neighbor timing tables", true);
 };
 
 WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)

Common/Tasks/centralityStudy.cxx

@@ -46,6 +46,7 @@ struct centralityStudy {
   Configurable<bool> requireIsVertexTOFmatched{"requireIsVertexTOFmatched", true, "require events with at least one of vertex contributors matched to TOF"};
   Configurable<bool> requireIsVertexTRDmatched{"requireIsVertexTRDmatched", true, "require events with at least one of vertex contributors matched to TRD"};
   Configurable<bool> rejectSameBunchPileup{"rejectSameBunchPileup", true, "reject collisions in case of pileup with another collision in the same foundBC"};
+  Configurable<float> minTimeDelta{"minTimeDelta", -1.0f, "reject collision if another collision is this close or less in time"};
 
   // Configurable Axes
   ConfigurableAxis axisMultFT0C{"axisMultFT0C", {2000, 0, 100000}, "FT0C amplitude"};
@@ -61,12 +62,13 @@ struct centralityStudy {
   // For centrality studies if requested
   ConfigurableAxis axisCentrality{"axisCentrality", {100, 0, 100}, "FT0C percentile"};
   ConfigurableAxis axisPVChi2{"axisPVChi2", {300, 0, 30}, "FT0C percentile"};
+  ConfigurableAxis axisDeltaTime{"axisDeltaTime", {300, 0, 300}, "#Delta time"};
 
   void init(InitContext&)
   {
     if (doprocessCollisions) {
       const AxisSpec axisCollisions{100, -0.5f, 99.5f, "Number of collisions"};
-      histos.add("hCollisionSelection", "hCollisionSelection", kTH1D, {{10, -0.5f, +9.5f}});
+      histos.add("hCollisionSelection", "hCollisionSelection", kTH1D, {{20, -0.5f, +19.5f}});
       histos.get<TH1>(HIST("hCollisionSelection"))->GetXaxis()->SetBinLabel(1, "All collisions");
       histos.get<TH1>(HIST("hCollisionSelection"))->GetXaxis()->SetBinLabel(2, "sel8 cut");
       histos.get<TH1>(HIST("hCollisionSelection"))->GetXaxis()->SetBinLabel(3, "posZ cut");
@@ -77,6 +79,7 @@ struct centralityStudy {
       histos.get<TH1>(HIST("hCollisionSelection"))->GetXaxis()->SetBinLabel(8, "kIsVertexTOFmatched");
       histos.get<TH1>(HIST("hCollisionSelection"))->GetXaxis()->SetBinLabel(9, "kIsVertexTRDmatched");
       histos.get<TH1>(HIST("hCollisionSelection"))->GetXaxis()->SetBinLabel(10, "kNoSameBunchPileup");
+      histos.get<TH1>(HIST("hCollisionSelection"))->GetXaxis()->SetBinLabel(11, "Neighbour rejection");
 
       histos.add("hFT0C_Collisions", "hFT0C_Collisions", kTH1D, {axisMultUltraFineFT0C});
       histos.add("hNPVContributors", "hNPVContributors", kTH1D, {axisMultUltraFinePVContributors});
@@ -99,6 +102,7 @@ struct centralityStudy {
       histos.add("hNITSOnlyTracksVsCentrality", "hNITSOnlyTracksVsCentrality", kTH2F, {axisCentrality, axisMultPVContributors});
       histos.add("hNGlobalTracksVsCentrality", "hNGlobalTracksVsCentrality", kTH2F, {axisCentrality, axisMultPVContributors});
       histos.add("hPVChi2VsCentrality", "hPVChi2VsCentrality", kTH2F, {axisCentrality, axisPVChi2});
+      histos.add("hDeltaTimeVsCentrality", "hDeltaTimeVsCentrality", kTH2F, {axisCentrality, axisDeltaTime});
     }
   }
 
@@ -152,6 +156,19 @@ struct centralityStudy {
     }
     histos.fill(HIST("hCollisionSelection"), 9 /* Not at same bunch pile-up */);
 
+    //do this only if information is available
+    if constexpr (requires { collision.timeToNext(); }) {
+      float timeToNeighbour = TMath::Min(
+        std::abs(collision.timeToNext()),
+        std::abs(collision.timeToPrevious())
+        );
+      histos.fill(HIST("hDeltaTimeVsCentrality"), collision.centFT0C(), timeToNeighbour);
+      if(timeToNeighbour<minTimeDelta){
+        return;
+      }
+      histos.fill(HIST("hCollisionSelection"), 10 /* has suspicious neighbour */);
+    }
+
     // if we got here, we also finally fill the FT0C histogram, please
     histos.fill(HIST("hNPVContributors"), collision.multPVTotalContributors());
     histos.fill(HIST("hFT0C_Collisions"), collision.multFT0C());
@@ -177,7 +194,7 @@ struct centralityStudy {
     genericProcessCollision(collision);
   }
 
-  void processCollisionsWithCentrality(soa::Join<aod::Mults, aod::MultsExtra, aod::MultSelections, aod::CentFT0Cs, aod::MultsGlobal>::iterator const& collision)
+  void processCollisionsWithCentrality(soa::Join<aod::Mults, aod::MultsExtra, aod::MultSelections, aod::CentFT0Cs, aod::MultsGlobal, aod::MultNeighs>::iterator const& collision)
   {
     genericProcessCollision(collision);
   }