KVRDP_INDRAFAZIA_DropUncorrEventsInFirstRunFile.cpp

#include "KVRDP_INDRAFAZIA_DropUncorrEventsInFirstRunFile.h"
#include "KVMultiDetArray.h"
 
 
 
 
Bool_t KVRDP_INDRAFAZIA_DropUncorrEventsInFirstRunFile::IsRequired(TString dataset, Int_t runnumber)
{
   fRunlist = GetDataSetEnv(dataset,"SkipUncorEventsAtRunBegin","");
   int index_multiplier = GetDataSetEnv(dataset,"DataSet.RunFileIndexMultiplier.raw", -1.0);
   auto first_file = (runnumber%index_multiplier == 0);
   runnumber/=index_multiplier;
   return (first_file && fRunlist.Contains(runnumber));
}
 
 
 
 
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 (i.e. without any" << std::endl;
    std::cout << "'run-index multiplier' being applied). 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)