KVSimReader_MMM_asym.cpp

//Created by KVClassFactory on Tue Jul  6 17:30:28 2010
//Author: Eric Bonnet
 
#include "KVSimReader_MMM_asym.h"
 
ClassImp(KVSimReader_MMM_asym)
 
 
 
 
KVSimReader_MMM_asym::KVSimReader_MMM_asym()
{
   // Default constructor
   init();
}
 
 
 
 
KVSimReader_MMM_asym::KVSimReader_MMM_asym(KVString filename)
{
   init();
   ConvertEventsInFile(filename);
}
 
 
 
 
void KVSimReader_MMM_asym::SetBoost(TVector3& vQP, TVector3& vQC)
{
   fApplyBoost = kTRUE;
   fBoostQP = vQP;
   fBoostQC = vQC;
}
 
 
 
 
Bool_t KVSimReader_MMM_asym::ReadEventHeader(KVSimEvent* evt)
{
   Int_t natt = 0;
   mult = 0;
 
   ReadLine(" ");
 
   Int_t res = GetNparRead();
   switch (res) {
      case 0:
         return kFALSE;
      default:
         //nlus = toks->GetEntries();
         idx = 0;
         mult = GetIntReadPar(idx++);
         natt = 5 * mult + 1;
         if (natt != res) {
            Info("ReadEvent", "Nombre de parametres (%d) different de celui attendu (%d)", res, natt);
            return kFALSE;
         }
         return kTRUE;
   }
   return kFALSE;
}
 
 
 
 
Bool_t KVSimReader_MMM_asym::ReadNucleus(KVSimNucleus* nuc)
{
 
   Int_t aa = GetIntReadPar(idx++);
   Int_t zz = GetIntReadPar(idx++);
 
   nuc->SetZ(zz);
   nuc->SetA(aa);
 
   Double_t px = GetDoubleReadPar(idx++);
   Double_t py = GetDoubleReadPar(idx++);
   Double_t pz = GetDoubleReadPar(idx++);
 
   nuc->SetMomentum(px, py, pz);
 
   if (fApplyBoost) {
      // add 2 kinematical frames to the particle:
      // the first one (labelled "QP") with the particle's momentum
      // in the QP Frame
      // the second one (labelled "QC") with the particle's momentum
      // in the QC Frame
      nuc->SetFrame("QP", {fBoostQP, kTRUE});
      nuc->SetFrame("QC", {fBoostQC, kTRUE});
   }
 
   return kTRUE;
 
}