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);
}
 
 
 
 
Bool_t KVSimReader_HIPSE_asym::ReadEventHeader(KVSimEvent* evt)
{
 
   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:
         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));
         h1->Fill(evt->GetParameters()->GetDoubleValue("Bparstore"));
         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;
   }
   return kTRUE;
}
 
 
 
 
Bool_t KVSimReader_HIPSE_asym::ReadNucleus(KVSimNucleus* nuc)
{
 
   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);
         }
         break;
      default:
         return kFALSE;
   }
 
   return true;
}
 
 
 
 
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",
                        proj.GetSymbol(), targ.GetSymbol(), ebeam));
   SetTreeTitle(Form("HIPSE secondary events %s + %s %.1f MeV/nuc.",
                   proj.GetSymbol(), targ.GetSymbol(), ebeam));
}