new changes

Author: JosueMtzGar

PWGMM/Lumi/Tasks/lumiStability.cxx

@@ -82,8 +82,9 @@ struct lumiStabilityTask {
     return false;
   }
 
-  void processFDDFT0(aod::FT0s const& ft0s, aod::FDDs const& fdds, BCsRun3 const& bcs)
+  void processMain(aod::FT0s const& ft0s, aod::FDDs const& fdds, aod::FV0As const& fv0s, BCsRun3 const& bcs)
   {
+    int CountNormal(0), CountPastProtec(0);
     for (auto const& fdd : fdds) {
       auto bc = fdd.bc_as<BCsRun3>();
       if (bc.timestamp() == 0) {
@@ -98,21 +99,41 @@ struct lumiStabilityTask {
       int maxDeltaBC = 5; // maximum difference
       bool pastActivityFDD = false;
       while (deltaBC < maxDeltaBC) {
+        deltaIndex++;
         if (bc.globalIndex() - deltaIndex < 0) {
           break;
         }
-        deltaIndex++;
         const auto& bc_past = bcs.iteratorAt(bc.globalIndex() - deltaIndex);
-        deltaBC = globalBC = bc_past.globalBC();
+        deltaBC = globalBC - bc_past.globalBC();
         if (deltaBC < maxDeltaBC) {
           pastActivityFDD |= bc_past.has_fdd();
         }
       }
+      deltaIndex = 0;
+      deltaBC = 0;
 
-      if (pastActivityFDD == false) {
+      bool futureActivityFDD = false;
+      while (deltaBC < maxDeltaBC) {
+        deltaIndex++;
+        if (bc.globalIndex() + deltaIndex >= bcs.size()) {
+          break;
+        }
+        const auto& bc_future = bcs.iteratorAt(bc.globalIndex() + deltaIndex);
+        deltaBC = bc_future.globalBC() - globalBC;
+        if (deltaBC < maxDeltaBC) {
+          futureActivityFDD |= bc_future.has_fdd();
+        }
+      }
+
+      CountNormal++;
+      if (pastActivityFDD == true || futureActivityFDD == true) {
+        CountPastProtec++;
         continue;
-        std::cout << "BC skipped!" << std::endl;
       }
+      /*if (pastActivityFDD == true) {
+        CountPastProtec++;
+        continue;
+      }*/
 
       std::bitset<8> fddTriggers = fdd.triggerMask();
       bool vertex = fddTriggers[o2::fdd::Triggers::bitVertex];
@@ -173,13 +194,51 @@ struct lumiStabilityTask {
 
     for (auto const& ft0 : ft0s) {
       auto bc = ft0.bc_as<BCsRun3>();
-      if (bc.timestamp() == false) {
+      if (bc.timestamp() == 0) {
         continue;
       }
 
       Long64_t globalBC = bc.globalBC();
       int localBC = globalBC % nBCsPerOrbit;
 
+      int deltaIndex = 0; // backward move counts
+      int deltaBC = 0;    // current difference wrt globalBC
+      int maxDeltaBC = 5; // maximum difference
+      bool pastActivityFT0 = false;
+      while (deltaBC < maxDeltaBC) {
+        deltaIndex++;
+        if (bc.globalIndex() - deltaIndex < 0) {
+          break;
+        }
+        const auto& bc_past = bcs.iteratorAt(bc.globalIndex() - deltaIndex);
+        deltaBC = globalBC - bc_past.globalBC();
+        if (deltaBC < maxDeltaBC) {
+          pastActivityFT0 |= bc_past.has_ft0();
+        }
+      }
+      deltaIndex = 0;
+      deltaBC = 0;
+
+      bool futureActivityFT0 = false;
+      while (deltaBC < maxDeltaBC) {
+        deltaIndex++;
+        if (bc.globalIndex() + deltaIndex >= bcs.size()) {
+          break;
+        }
+        const auto& bc_future = bcs.iteratorAt(bc.globalIndex() + deltaIndex);
+        deltaBC = bc_future.globalBC() - globalBC;
+        if (deltaBC < maxDeltaBC) {
+          futureActivityFT0 |= bc_future.has_ft0();
+        }
+      }
+
+      if (pastActivityFT0 == true || futureActivityFT0 == true) {
+        continue;
+      }
+      /*if (pastActivityFT0 == true) {
+        continue;
+      }*/
+
       std::bitset<8> fT0Triggers = ft0.triggerMask();
       bool vertex = fT0Triggers[o2::fdd::Triggers::bitVertex];
 
@@ -207,21 +266,55 @@ struct lumiStabilityTask {
         }
       }
     } // loop over FT0 events
-  }   // end processFDDFT0
-
-  PROCESS_SWITCH(lumiStabilityTask, processFDDFT0, "Process FDD and FT0 to lumi stability analysis", true);
 
-  void processV0(aod::FV0As const& fv0s, aod::BCsWithTimestamps const&)
-  {
     for (auto const& fv0 : fv0s) {
-      auto bc = fv0.bc_as<BCsWithTimestamps>();
-      if (bc.timestamp() == false) {
+      auto bc = fv0.bc_as<BCsRun3>();
+      if (bc.timestamp() == 0) {
         continue;
       }
 
       Long64_t globalBC = bc.globalBC();
       int localBC = globalBC % nBCsPerOrbit;
 
+      int deltaIndex = 0; // backward move counts
+      int deltaBC = 0;    // current difference wrt globalBC
+      int maxDeltaBC = 5; // maximum difference
+      bool pastActivityV0A = false;
+      while (deltaBC < maxDeltaBC) {
+        deltaIndex++;
+        if (bc.globalIndex() - deltaIndex < 0) {
+          break;
+        }
+        const auto& bc_past = bcs.iteratorAt(bc.globalIndex() - deltaIndex);
+        deltaBC = globalBC - bc_past.globalBC();
+        if (deltaBC < maxDeltaBC) {
+          pastActivityV0A |= bc_past.has_fv0a();
+        }
+      }
+      deltaIndex = 0;
+      deltaBC = 0;
+
+      bool futureActivityV0A = false;
+      while (deltaBC < maxDeltaBC) {
+        deltaIndex++;
+        if (bc.globalIndex() + deltaIndex >= bcs.size()) {
+          break;
+        }
+        const auto& bc_future = bcs.iteratorAt(bc.globalIndex() + deltaIndex);
+        deltaBC = bc_future.globalBC() - globalBC;
+        if (deltaBC < maxDeltaBC) {
+          futureActivityV0A |= bc_future.has_fv0a();
+        }
+      }
+
+      if (pastActivityV0A == true || futureActivityV0A == true) {
+        continue;
+      }
+
+      /*if (pastActivityV0A == true) {
+        continue;
+      }*/
+
       std::bitset<8> fv0Triggers = fv0.triggerMask();
       bool aOut = fv0Triggers[o2::fdd::Triggers::bitAOut];
       bool aIn = fv0Triggers[o2::fdd::Triggers::bitAIn];
@@ -247,9 +340,12 @@ struct lumiStabilityTask {
         }
       }
     } // loop over V0 events
-  }   // end processV0
+    std::cout << "************ >>>>>>>>>>>>>> "
+              << "Whithout Past Protection: " << CountNormal << " "
+              << "Avoided Whith Past Protection: " << CountPastProtec << "<<<<<<<<<<<< ********************" << std::endl;
+  } // end processMain
 
-  PROCESS_SWITCH(lumiStabilityTask, processV0, "Process V0 to lumi stability analysis", true);
+  PROCESS_SWITCH(lumiStabilityTask, processMain, "Process FDD and FT0 to lumi stability analysis", true);
 };
 
 WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)