Replaced first time interval with job average

Author: jackal1-66

MC/config/common/external/generator/TPCLoopers.C

@@ -218,6 +218,11 @@ class GenTPCLoopers : public Generator
                 LOG(info) << "Interaction Time records has " << mInteractionTimeRecords.size() << " entries.";
                 mCollisionContext->printCollisionSummary();
             }
+            for (int c = 0; c < mInteractionTimeRecords.size() - 1; c++)
+            {
+                mIntTimeRecMean += mInteractionTimeRecords[c + 1].bc2ns() - mInteractionTimeRecords[c].bc2ns();
+            }
+            mIntTimeRecMean /= (mInteractionTimeRecords.size() - 1); // Average interaction time record used for the reference
         }
 
         Bool_t generateEvent() override
@@ -230,11 +235,11 @@ class GenTPCLoopers : public Generator
             {
                 unsigned int nLoopers, nLoopersPairs, nLoopersCompton;
                 LOG(info) << "mCurrentEvent is " << mCurrentEvent;
-                LOG(info) << "Current event time: " << ((mCurrentEvent < mInteractionTimeRecords.size() - 1) ? std::to_string(mInteractionTimeRecords[mCurrentEvent + 1].bc2ns() - mInteractionTimeRecords[mCurrentEvent].bc2ns()) : "Final Event till the end") << " ns";
+                LOG(info) << "Current event time: " << ((mCurrentEvent < mInteractionTimeRecords.size() - 1) ? std::to_string(mInteractionTimeRecords[mCurrentEvent + 1].bc2ns() - mInteractionTimeRecords[mCurrentEvent].bc2ns()) : std::to_string(mIntTimeRecMean)) << " ns";
                 LOG(info) << "Current time offset wrt BC: " << mInteractionTimeRecords[mCurrentEvent].getTimeOffsetWrtBC() << " ns";
-                mTimeLimit = (mCurrentEvent < mInteractionTimeRecords.size() - 1) ? mInteractionTimeRecords[mCurrentEvent + 1].bc2ns() - mInteractionTimeRecords[mCurrentEvent].bc2ns() : mInteractionTimeRecords[1].bc2ns() - mInteractionTimeRecords[0].bc2ns();
+                mTimeLimit = (mCurrentEvent < mInteractionTimeRecords.size() - 1) ? mInteractionTimeRecords[mCurrentEvent + 1].bc2ns() - mInteractionTimeRecords[mCurrentEvent].bc2ns() : mIntTimeRecMean;
                 // With Flat Gas number of loopers are adapted based on time interval widths
-                nLoopers = mFlatGasNumber * (mTimeLimit / mStartingTime);
+                nLoopers = mFlatGasNumber * (mTimeLimit / mIntTimeRecMean);
                 nLoopersPairs = static_cast<unsigned int>(std::round(nLoopers * mLoopsFractionPairs));
                 nLoopersCompton = nLoopers - nLoopersPairs;
                 SetNLoopers(nLoopersPairs, nLoopersCompton);
@@ -421,14 +426,13 @@ class GenTPCLoopers : public Generator
                 else
                 {
                     mFlatGasNumber = number;
-                    mStartingTime = mInteractionTimeRecords[1].bc2ns() - mInteractionTimeRecords[0].bc2ns(); // Time width of first event
                 }
             }
             else
             {
                 mFlatGasNumber = -1;
             }
-            LOG(info) << "Flat gas loopers: " << (mFlatGas ? "ON" : "OFF") << ", Number of loopers per event: " << mFlatGasNumber;
+            LOG(info) << "Flat gas loopers: " << (mFlatGas ? "ON" : "OFF") << ", Reference loopers number per event: " << mFlatGasNumber;
         }
 
         void setFractionPairs(float &fractionPairs)
@@ -468,7 +472,7 @@ class GenTPCLoopers : public Generator
         std::vector<o2::InteractionTimeRecord> mInteractionTimeRecords; // Interaction time records from collision context
         Bool_t mFlatGas = false;                                        // Flag to indicate if flat gas loopers are used
         Int_t mFlatGasNumber = -1;                                      // Number of flat gas loopers per event
-        double mStartingTime = 1.0;                                     // First BC time interval in ns used for reference in adaptive loopers
+        double mIntTimeRecMean = 1.0;                                   // Average interaction time record used for the reference
         double mTimeLimit = 0.0;                                        // Time limit for the current event
         float mLoopsFractionPairs = 0.08;                               // Fraction of loopers from Pairs
 };