KVSimReader_MMM.cpp

//Created by KVClassFactory on Tue Jul  6 16:48:05 2010
//Author: Eric Bonnet
 
#include "KVSimReader_MMM.h"
 
ClassImp(KVSimReader_MMM)
 
 
 
 
KVSimReader_MMM::KVSimReader_MMM()
{
   // Default constructor
   init();
}
 
 
 
 
KVSimReader_MMM::KVSimReader_MMM(KVString filename)
{
   init();
   ConvertEventsInFile(filename);
}
 
 
 
 
Bool_t KVSimReader_MMM::ReadEventHeader(KVSimEvent* evt)
{
   Int_t natt = 0;
   mult = 0;
 
   ReadLine(" ");
 
   Int_t res = GetNparRead();
   switch (res) {
      case 0:
         return kFALSE;
      default:
         idx = 0;
         mult = GetIntReadPar(idx++);
         natt = 6 * mult + 4 + 1;
         if (natt != res) {
            Info("ReadEvent", "Nombre de parametres (%d) different de celui attendu (%d)", res, natt);
            return kFALSE;
         }
         evt->GetParameters()->SetValue("Eth", GetDoubleReadPar(idx++));
         evt->GetParameters()->SetValue("Eclb", GetDoubleReadPar(idx++));
         evt->GetParameters()->SetValue("Eexci", GetDoubleReadPar(idx++));
         evt->GetParameters()->SetValue("Ebind", GetDoubleReadPar(idx++));
         return kTRUE;
   }
   return kFALSE;
}
 
 
 
 
Bool_t KVSimReader_MMM::ReadEvent(KVSimEvent* evt)
{
   for (Int_t mm = 0; mm < mult; mm += 1) {
      auto nuc = evt->AddParticle();
      ReadNucleus(nuc);
   }
 
   evt->SetNumber(GetEventNumber());
 
   return kTRUE;
}
 
 
 
 
Bool_t KVSimReader_MMM::ReadNucleus(KVSimNucleus* nuc)
{
 
   Int_t aa = GetIntReadPar(idx++);
   nuc->SetZ(GetIntReadPar(idx++));
   nuc->SetA(aa);
   Double_t exci = GetDoubleReadPar(idx++);
   Double_t px = GetDoubleReadPar(idx++);
   Double_t py = GetDoubleReadPar(idx++);
   Double_t pz = GetDoubleReadPar(idx++);
 
   nuc->SetExcitEnergy(exci);
   nuc->SetMomentum(px, py, pz);
 
   return kTRUE;
 
}