KVSimReader_HIPSE.cpp

//Created by KVClassFactory on Wed Jul  7 12:15:56 2010
//Author: bonnet
 
#include "KVSimReader_HIPSE.h"
 
ClassImp(KVSimReader_HIPSE)
 
 
 
 
void KVSimReader_HIPSE::define_output_filename()
{
   SetROOTFileName(Form("HIPSE_%s_%s_%.1fAMeV",
                        proj.GetSymbol(), targ.GetSymbol(), ebeam));
   SetTreeTitle(Form("HIPSE primary events %s + %s %.1f MeV/nuc.",
                     proj.GetSymbol(), targ.GetSymbol(), ebeam));
}
 
 
 
 
KVSimReader_HIPSE::KVSimReader_HIPSE()
{
   // Default constructor
   init();
}
 
 
 
 
KVSimReader_HIPSE::KVSimReader_HIPSE(KVString filename)
{
   init();
   ConvertEventsInFile(filename);
}
 
 
 
 
Bool_t KVSimReader_HIPSE::ReadFileHeader()
{
 
   // create histogram for impact parameter distribution
   h1 = new TH1F("impact_parameter", "distri", 200, 0, 20);
   AddObject(h1);
   h1->SetDirectory(nullptr);//avoid double deletion on ROOT file close
 
   auto res = ReadLineAndCheck(2, " ");
   switch (res) {
      case KVFileReader::ReadStatus::EmptyLine:
         return kFALSE;
      case KVFileReader::ReadStatus::OK:
         AddInfo({"Aproj", GetIntReadPar(0)});
         AddInfo({"Zproj", GetIntReadPar(1)});
         proj.SetZandA(GetReadPar(1).Atoi(), GetReadPar(0).Atoi());
         break;
      default:
         return kFALSE;
   }
 
   res = ReadLineAndCheck(2, " ");
   switch (res) {
      case KVFileReader::ReadStatus::EmptyLine:
         return kFALSE;
      case KVFileReader::ReadStatus::OK:
         AddInfo({"Atarg", GetIntReadPar(0)});
         AddInfo({"Ztarg", GetIntReadPar(1)});
         targ.SetZandA(GetReadPar(1).Atoi(), GetReadPar(0).Atoi());
         break;
      default:
         return kFALSE;
   }
 
   res = ReadLineAndCheck(1, " ");
   switch (res) {
      case KVFileReader::ReadStatus::EmptyLine:
         return kFALSE;
      case KVFileReader::ReadStatus::OK:
         AddInfo({"Ebeam", GetDoubleReadPar(0)});
         ebeam = GetDoubleReadPar(0);
         break;
      default:
         return kFALSE;
   }
 
   define_output_filename();
 
   KVError::Info(this, "ReadFileHeader", "Reading %s", GetTreeTitle().Data());
 
   return kTRUE;
}
 
 
 
 
Bool_t KVSimReader_HIPSE::ReadEventHeader(KVSimEvent* evt)
{
 
   mult = 0, mtotal = 0;
 
   /*---------------------------------------------
     mul_vrai = multiplicity of charged particles
     mult     = total multiplicity (i.e. including neutrons)
   //---------------------------------------------
   */
   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);   //mul_vrai
         mtotal = GetIntReadPar(1);
         evt->GetParameters()->SetValue("mult", mtotal);
         break;
      default:
         return kFALSE;
   }
 
   /*---------------------------------------------
     Esa       = excitation per nucleon
     vcm       = center of mass energy
     Bparstore = impact parameter
   //---------------------------------------------
   */
   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"));
         break;
      default:
         return kFALSE;
   }
 
   /*---------------------------------------------
     energetic information
     excitat        : total excitation energy
     xmassav        : Q-value
     ekinav         : total kinetic energy at freeze-out
     epotav         : total potential energy at freeze-out
     erotav         : total rotational energy at freeze-out
   //---------------------------------------------
   */
   res = ReadLineAndCheck(2, " ");
   switch (res) {
      case KVFileReader::ReadStatus::EmptyLine:
         return kFALSE;
      case KVFileReader::ReadStatus::OK:
         evt->GetParameters()->SetValue("excitat", GetDoubleReadPar(0));
         evt->GetParameters()->SetValue("xmassav", GetDoubleReadPar(1));
         break;
      default:
         return kFALSE;
   }
 
   res = ReadLineAndCheck(3, " ");
   switch (res) {
      case KVFileReader::ReadStatus::EmptyLine:
         return kFALSE;
      case KVFileReader::ReadStatus::OK:
         evt->GetParameters()->SetValue("ekinav", GetDoubleReadPar(0));
         evt->GetParameters()->SetValue("epotav", GetDoubleReadPar(1));
         evt->GetParameters()->SetValue("erotav", GetDoubleReadPar(2));
         break;
 
      default:
         return kFALSE;
   }
   return kTRUE;
}
 
 
 
 
Bool_t KVSimReader_HIPSE::ReadEvent(KVSimEvent* evt)
{
   evt->SetNumber(GetEventNumber());
   for (Int_t mm = 0; mm < mult; mm += 1) {
      auto nuc = evt->AddParticle();
      if (!ReadNucleus(nuc)) return kFALSE;
   }
 
   return kTRUE;
}
 
 
 
 
Bool_t KVSimReader_HIPSE::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));
         break;
 
      default:
         return kFALSE;
   }
 
   int AA = nuc->GetA();
 
   /*---------------------------------------------------------
     exci(I)  : excitation energy
     ether(I) : not used
     spinx(i),spiny(i),spinz(i) : angular momentum (hbar units)
   //---------------------------------------------
   */
   res = ReadLineAndCheck(2, " ");
   switch (res) {
      case KVFileReader::ReadStatus::EmptyLine:
         Info("ReadNucleus", "case 0 line est vide");
         return kFALSE;
 
      case KVFileReader::ReadStatus::OK:
         nuc->SetExcitEnergy(GetDoubleReadPar(0)*AA);
         nuc->GetParameters()->SetValue("exci", GetDoubleReadPar(0));
         nuc->GetParameters()->SetValue("ether", GetDoubleReadPar(1));
         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:
         //On effectue la meme rotation que les impulsions ... à vérifier
         nuc->SetAngMom(GetDoubleReadPar(1), GetDoubleReadPar(2), GetDoubleReadPar(0));
         break;
 
      default:
         return kFALSE;
   }
 
   return kTRUE;
}
 
 
 
 
void KVSimReader_HIPSE::NotifyChangeOfOutputFile()
{
   // Clear the impact parameter histogram when starting a new output file.
   //
   // This way the distribution in the histogram only concerns the events in the current file.
 
   h1->Reset();
}