KVSimReader_HIPSE_asym.cpp

//Created by KVClassFactory on Fri Jul  2 15:16:15 2010
//Author: bonnet
 
#include "KVSimReader_HIPSE_asym.h"
#include "KV2Body.h"
#include "TRandom.h"
 
ClassImp(KVSimReader_HIPSE_asym)
 
 
 
 
 
KVSimReader_HIPSE_asym::KVSimReader_HIPSE_asym()
{
   // Default constructor
   init();
}
 
 
 
 
KVSimReader_HIPSE_asym::KVSimReader_HIPSE_asym(KVString filename) : KVSimReader_HIPSE()
{
   init();
   ConvertEventsInFile(filename);
}
 
 
 
 
KVSimReader_HIPSE_asym::~KVSimReader_HIPSE_asym()
{
   // Destructor
}
 
 
 
 
Bool_t KVSimReader_HIPSE_asym::ReadEvent()
{
 
   evt->Clear();
   Int_t mult = 0, mtotal = 0;
   auto res = ReadLineAndCheck(2, " ");
   switch (res) {
      case KVFileReader::ReadStatus::EmptyLine:
         Info("ReadEvent", "case 0 line est vide");
         return kFALSE;
      case KVFileReader::ReadStatus::OK:
         evt->SetNumber(nevt);
         mult = GetIntReadPar(0);
         mtotal = GetIntReadPar(1);
         evt->GetParameters()->SetValue("mult", mtotal);
 
         break;
      default:
 
         return kFALSE;
   }
 
   /*---------------------------------------------
      Esa       = excitation per nucleon
      vcm       = center of mass energy
      Bparstore = impact parameter
      PhiPlan   = angle of the reaction plane
   //---------------------------------------------
   */
 
   fPhiPlan = gRandom->Rndm() * 2.*TMath::Pi();
   rr.SetXEulerAngles(0., 0., fPhiPlan);
 
   res = ReadLineAndCheck(3, " ");
   switch (res) {
      case KVFileReader::ReadStatus::EmptyLine:
         return kFALSE;
      case KVFileReader::ReadStatus::OK:
         evt->GetParameters()->SetValue("Esa", GetDoubleReadPar(0));
         evt->GetParameters()->SetValue("vcm", GetDoubleReadPar(1));
         evt->GetParameters()->SetValue("Bparstore", GetDoubleReadPar(2));
         evt->GetParameters()->SetValue("PhiPlan", fPhiPlan);
         break;
      default:
 
         return kFALSE;
   }
 
 
   res = ReadLineAndCheck(3, " ");
   switch (res) {
      case KVFileReader::ReadStatus::EmptyLine:
         return kFALSE;
      case KVFileReader::ReadStatus::OK:
 
         break;
      default:
 
         return kFALSE;
   }
 
   evt->SetNumber(nevt);
   for (Int_t mm = 0; mm < mult; mm += 1) {
      nuc = (KVSimNucleus*)evt->AddParticle();
      if (!ReadNucleus()) return kFALSE;
   }
 
   nevt += 1;
   return kTRUE;
 
}
 
 
 
 
Bool_t KVSimReader_HIPSE_asym::ReadNucleus()
{
 
   auto res = ReadLineAndCheck(3, " ");
   switch (res) {
      case KVFileReader::ReadStatus::EmptyLine:
         Info("ReadNucleus", "case 0 line est vide");
         return kFALSE;
 
      case KVFileReader::ReadStatus::OK:
         /*
         proven = 0 -> fusion of the QP and QT
         proven = 1 -> QP
         proven = 2 -> QT
         proven > 2 -> other
         */
 
         nuc->SetZ(GetIntReadPar(1));
         nuc->SetA(GetIntReadPar(0));
         nuc->GetParameters()->SetValue("proven", GetDoubleReadPar(2));
 
         break;
 
      default:
 
 
         return kFALSE;
   }
 
   res = ReadLineAndCheck(3, " ");
   switch (res) {
      case KVFileReader::ReadStatus::EmptyLine:
         Info("ReadNucleus", "case 0 line est vide");
         return kFALSE;
 
      case KVFileReader::ReadStatus::OK:
         //Axe "faisceau dans HIPSE x -> on effectue une rotation X,Y,Z -> Y,Z,X"
         nuc->SetMomentum(GetDoubleReadPar(1), GetDoubleReadPar(2), GetDoubleReadPar(0));
         {
            TVector3 vv = nuc->GetVelocity();
            vv *= rr;
            nuc->SetVelocity(vv);
            return kTRUE;
         }
      default:
 
 
         return kFALSE;
   }
 
 
}
 
 
 
 
void KVSimReader_HIPSE_asym::define_output_filename()
{
   // ROOT file called: HIPSE_[PROJ]_[TARG]_[EBEAM]AMeV_ASYM.root
   // Call after reading file header
 
   SetROOTFileName(Form("HIPSE_%s_%s_%.1fAMeV_ASYM.root",
                        proj.GetSymbol(), targ.GetSymbol(), ebeam));
}
 
 
 
 
 
void KVSimReader_HIPSE_asym::ConvertEventsInFile(KVString filename)
{
   if (!OpenFileToRead(filename)) return;
   if (!ReadHeader()) return;
   define_output_filename();
   tree_title.Form("HIPSE secondary events %s + %s %.1f MeV/nuc.",
                   proj.GetSymbol(), targ.GetSymbol(), ebeam);
   Run();
   CloseFile();
}