KVRDP_INDRAFAZIA_DropUncorrEventsInFirstRunFile.cpp

#include "KVRDP_INDRAFAZIA_DropUncorrEventsInFirstRunFile.h"
#include "KVMultiDetArray.h"
 
 
 
 
Bool_t KVRDP_INDRAFAZIA_DropUncorrEventsInFirstRunFile::IsRequired(TString dataset, const run_index_t& runnumber)
{
   fRunlist.Set(GetDataSetEnv(dataset, "SkipUncorEventsAtRunBegin", ""));
   auto first_file = !runnumber.has_index();
   return (first_file && fRunlist.Contains(runnumber.run()));
}
 
 
 
 
void KVRDP_INDRAFAZIA_DropUncorrEventsInFirstRunFile::PrintPatchInfo() const
{
   std::cout << "Avoid using uncorrelated INDRA/FAZIA events at beginning of run" << std::endl;
 
   std::cout << "\nDue to the way Narval handles data buffers, at the beginning of a new run" << std::endl;
   std::cout << "('run' in the sense of Narval: a run begins when 'Start' is pressed, it can" << std::endl;
   std::cout << "be written in several files) there is typically a buffer of only INDRA or" << std::endl;
   std::cout << "only FAZIA events (a few hundreds), uncorrelated with anything else." << std::endl;
 
   std::cout << "\nIn order to be sure to only treat correlated data which could have been" << std::endl;
   std::cout << "merged if required (if the INDRA-FAZIA time coincidence is in the window)," << std::endl;
   std::cout << "this patch will skip all events until the first merged INDRA-FAZIA event" << std::endl;
   std::cout << "is read, thereafter all events are accepted." << std::endl;
 
   std::cout << "\nThis is only applied to the runs found in the following dataset variable:" << std::endl;
 
   std::cout << "\n~~~" << std::endl;
   std::cout << " [dataset].SkipUncorEventsAtRunBegin: [runlist]" << std::endl;
   std::cout << "~~~" << std::endl;
 
   std::cout << "\nwhere the run numbers to be given are those of Narval. Obviously, when no correlation is" << std::endl;
   std::cout << "expected (alpha source or other calibration runs, or running independently" << std::endl;
   std::cout << "without merger), the condition should not be applied (it could reject all" << std::endl;
   std::cout << "events)." << std::endl;
}
 
 
 
 
Bool_t KVRDP_INDRAFAZIA_DropUncorrEventsInFirstRunFile::SkipEvent(KVMultiDetArray* mda)
{
   if (!fReject) return kFALSE;
   fReject = !(mda->GetArray("INDRA")->HandledRawData() && mda->GetArray("FAZIA")->HandledRawData());
   if (fReject) {
      ++events_skipped;
   }
   else
      Info("SkipEvent", "Beginning analysis after skipping %d events...", events_skipped);
 
   return fReject;
}
 
 
ClassImp(KVRDP_INDRAFAZIA_DropUncorrEventsInFirstRunFile)