KVDirectivity.cpp

#include "KVDirectivity.h"
 
ClassImp(KVDirectivity)
 
 
 
 
void KVDirectivity::Copy(TObject& a) const
{
 
   KVVarGlob1::Copy(a);// copy attributes of KVVarGlob1 base object
 
   KVDirectivity& aglob = (KVDirectivity&)a;
   aglob.total_mod_trans = total_mod_trans;
   aglob.total_trans = total_trans;
}
 
 
 
 
 
void KVDirectivity::Calculate(void)
{
   double mod_tot_trans = total_trans.Mag();
   SetValue(mod_tot_trans / total_mod_trans);
}
 
 
 
 
 
void KVDirectivity::init()
{
   fType = KVVarGlob::kOneBody; // this is a 1-body variable
#ifdef USING_ROOT6
   SetSelection({"V>Vcm", [](const KVNucleus * n)
   {
      return n->GetVpar() > 0;
   }
                });
#else
   SetSelection("_NUC_->GetVpar()>0");
#endif
}
 
 
 
 
 
void KVDirectivity::fill(const KVNucleus* n)
{
   TVector3 trans_mom = n->GetTransverseMomentum();
   total_trans += trans_mom;
   total_mod_trans += trans_mom.Mag();
}