KVSimReader_SMF.cpp

//Created by KVClassFactory on Mon Jul  5 16:21:41 2010
//Author: bonnet
 
#include "KVSimReader_SMF.h"
 
ClassImp(KVSimReader_SMF)
 
 
 
KVSimReader_SMF::KVSimReader_SMF()
{
   // Default constructor
   init();
}
 
 
 
 
KVSimReader_SMF::KVSimReader_SMF(KVString filename)
{
   init();
   ConvertEventsInFile(filename);
}
 
 
 
 
 
Bool_t KVSimReader_SMF::ReadEventHeader(KVSimEvent* evt)
{
 
   KVString snom;
   auto res = ReadLineAndCheck(1, " ");
   switch (res) {
      case KVFileReader::ReadStatus::EmptyLine:
         return kFALSE;
      case KVFileReader::ReadStatus::OK:
         snom.Form("%s", GetReadPar(0).Data());
         snom.ReplaceAll("evt_", "");
         //Info("ReadHeader","lecture %d",snom.Atoi());
         event_number = snom.Atoi();
         break;
      default:
         return kFALSE;
   }
   return kTRUE;
}
 
 
 
 
Bool_t KVSimReader_SMF::ReadEvent(KVSimEvent* evt)
{
   mult = 0;
   auto res = ReadLineAndCheck(1, " ");
   switch (res) {
      case KVFileReader::ReadStatus::EmptyLine:
         return kFALSE;
      case KVFileReader::ReadStatus::OK:
         mult = GetIntReadPar(0);
         break;
      default:
         return kFALSE;
   }
   evt->SetNumber(event_number);
   for (Int_t mm = 0; mm < mult; mm += 1) {
      auto nuc = evt->AddParticle();
      if (!ReadNucleus(nuc)) return kFALSE;
   }
   return kTRUE;
}
 
 
 
 
Bool_t KVSimReader_SMF::ReadNucleus(KVSimNucleus* nuc)
{
 
   ReadLine(" ");
   Int_t res = GetNparRead();
   Int_t npar = 0;
   switch (res) {
      case 0:
         Info("ReadNucleus", "case 0 line est vide");
         return kFALSE;
 
      default:
 
         ReadLineAndAdd(" ");
         npar = GetNparRead();
         if (npar != 12) {
            Info("ReadNucleus", "Nombre de parametres (%d) different de celui attendu(%d)", npar, 12);
            return kFALSE;
         }
         nuc->SetZ(TMath::Nint(GetDoubleReadPar(1)));
         nuc->SetA(TMath::Nint(GetDoubleReadPar(0)));
         nuc->GetParameters()->SetValue("density", GetDoubleReadPar(2));
         nuc->SetExcitEnergy(GetDoubleReadPar(3));
         nuc->GetParameters()->SetValue("spin", GetDoubleReadPar(4));
         //Axe "faisceau dans SMF z"
         nuc->SetMomentum(GetDoubleReadPar(5), GetDoubleReadPar(6), GetDoubleReadPar(7));
         nuc->SetPosition(GetDoubleReadPar(8), GetDoubleReadPar(9), GetDoubleReadPar(10));
         nuc->GetParameters()->SetValue("avoir", GetDoubleReadPar(11));
 
         return kTRUE;
   }
 
}