KVSubEventMaker.cpp

#include "KVSubEventMaker.h"
#include <random>
#include <algorithm>
ClassImp(KVSubEventMaker)
 
 
 
 
void KVSubEventMaker::Copy(TObject& a) const
{
   KVVarGlob::Copy(a);// copy attributes of KVVarGlob base object
 
   // Now copy any additional attributes specific to this class:
   // To copy a specific field, do as follows:
   //
   //     aglob.field=field;
   //
   // If setters and getters are available, you can also proceed as follows
   //
   //    aglob.SetField(GetField());
   //
}
 
 
 
 
 
void KVSubEventMaker::Init(void)
{
}
 
 
 
 
 
void KVSubEventMaker::Reset(void)
{
   particles.clear();
}
 
 
 
 
 
void KVSubEventMaker::Calculate(void)
{
   // Set group for required fraction of particles (starting from the first one)
 
   auto mult = particles.size();
   // perform a random shuffle of the particles in case there is a correlation
   // between position in the event and type of particle, where it was detected,
   // etc. etc.
   std::random_device rd;
   std::mt19937 g(rd());
   std::shuffle(particles.begin(), particles.end(), g);
 
   auto required = mult * GetParameter("FRACTION");
   // nothing guarantees that the required number of particles is a whole number
   auto req_int = std::floor(required);
   auto req_frac = required - req_int;
   // if so, we use req_int particles with probability 1-req_frac and req_int+1
   // particles with probability req_frac. The mean number of particles used
   // is exactly the value 'required'.
   if (req_frac > 1.e-10) {
      if (gRandom->Uniform() < req_frac) ++req_int;
   }
   int i = 0;
   for (auto p : particles) {
      if (i < req_int) p->AddGroup(GetName());
      ++i;
   }
}
 
 
 
 
 
Double_t KVSubEventMaker::getvalue_int(Int_t index) const
{
 
   return 0;
}
 
 
 
 
void KVSubEventMaker::fill(const KVNucleus* n)
{
   // add to internal list
   particles.push_back(const_cast<KVNucleus*>(n));
}
 
 
 
 
 
void KVSubEventMaker::init()
{
   // PRIVATE method
   fType = KVVarGlob::kOneBody; // this is a N-body variable
}