KVSimReader_HaC.cpp

#include "KVSimReader_HaC.h"
 
 
 
KVSimReader_HaC::KVSimReader_HaC()
{
   tree_name = "HaC";
   tree_title = "Hybrid alpha-Cluster model";
   root_file_name = "HaC_events.root";
}
 
 
 
 
KVSimReader_HaC::KVSimReader_HaC(const TString& filename)
   : KVSimReader_HaC()
{
   ConvertAndSaveEventsInFile(filename);
}
 
 
 
 
void KVSimReader_HaC::ReadFile()
{
   while (IsOK()) {
      while (ReadEvent()) {
         if (nevt % 1000 == 0) Info("ReadFile", "%d evts lus", nevt);
         if (HasToFill()) FillTree();
      }
   }
}
 
 
 
 
Bool_t KVSimReader_HaC::ReadEvent()
{
   // Format of events:
   //   1  1001.0000       0       7    -32.98796701      1.90742430 => event number,simulation time,l-value (entrance channel),number of fragments,?,?
   //1    -21.80360935     37.17158811     77.88473322     -0.46023991      0.78299889      1.66679455    -28.30000000     18.65061322      0.00000000
   //nalpha,xcm,ycm,zcm,kxcm,kycm,kzcm,total internal energy,kinetic cm energy,internal kinetic energy
   //1  => id of each alpha in cluster
   //2    -32.98385616     18.06693460     32.40992559     -0.71499466      0.36699478      0.70946383    -58.66470598     11.91719999      0.20506607
   //2    7  => id of each alpha in cluster
   //3      8.23280728    -52.40935312     51.83598017      0.17541902     -1.11296131      1.12388664    -86.46620999     39.34742042      1.90408960
   //3    4    6 => id of each alpha in cluster
   //1      0.08819134    105.14792154     49.71340857     -0.01218202      2.29301067      1.10528160    -28.30000000     33.47455814      0.00000000
   //5
   //1     82.12723574    -37.59507061   -129.88880017      1.75078647     -0.83889006     -2.81229792    -28.30000000     60.11693332      0.00000000
   //8
   //1    -84.28224914     -0.48406407   -130.11314158     -1.80061962     -0.00914512     -2.82881943    -28.30000000     57.90240833      0.00000000
   //9
   //1     65.13972189     16.85381518    -87.92399172      1.42598734      0.37691999     -1.92154637    -28.30000000     30.32639124      0.00000000
   //10
   //
   //    2  1001.0000       0       8    -32.98597627      2.21586178 ==> next event
 
   evt->Clear();
   auto res = ReadLineAndCheck(6," ");
   int nfra = 0;
   switch (res) {
      case KVFileReader::ReadStatus::EmptyLine:
      case KVFileReader::ReadStatus::CommentLine:
         return kTRUE;
 
      case KVFileReader::ReadStatus::OK:
         nevt = GetIntReadPar(0);
         evt->GetParameters()->SetValue("SimulationTime",GetDoubleReadPar(1));
         evt->GetParameters()->SetValue("AngularMomentum", GetIntReadPar(2));
         nfra = GetIntReadPar(3);
         evt->SetNumber(nevt);
         evt->GetParameters()->SetValue("nfra", nfra);
         break;
 
      case KVFileReader::ReadStatus::EndOfFile:
         KVError::Info(this,"ReadEvent","End of file reached after %d events",nevt);
         return kFALSE;
 
      default:
         KVError::Error(this,"ReadEvent","ReadLineAndCheck status = %s",GetReadStatus(res).Data());
         return kFALSE;
   }
   if(nfra)
   {
      for(int i=0; i<nfra; ++i)
      {
         nuc = (KVSimNucleus*)evt->AddParticle();
         if(!ReadNucleus())
            return kFALSE;
      }
   }
   return kTRUE;
}
 
 
 
 
Bool_t KVSimReader_HaC::ReadNucleus()
{
   // Format of events:
   //   1  1001.0000       0       7    -32.98796701      1.90742430 => event number,simulation time,l-value (entrance channel),number of fragments,?,?
   //1    -21.80360935     37.17158811     77.88473322     -0.46023991      0.78299889      1.66679455    -28.30000000     18.65061322      0.00000000
   //nalpha,xcm,ycm,zcm,kxcm,kycm,kzcm,total internal energy,kinetic cm energy,internal kinetic energy
   //1  => id of each alpha in cluster
   //2    -32.98385616     18.06693460     32.40992559     -0.71499466      0.36699478      0.70946383    -58.66470598     11.91719999      0.20506607
   //2    7  => id of each alpha in cluster
   //3      8.23280728    -52.40935312     51.83598017      0.17541902     -1.11296131      1.12388664    -86.46620999     39.34742042      1.90408960
   //3    4    6 => id of each alpha in cluster
   //1      0.08819134    105.14792154     49.71340857     -0.01218202      2.29301067      1.10528160    -28.30000000     33.47455814      0.00000000
   //5
   //1     82.12723574    -37.59507061   -129.88880017      1.75078647     -0.83889006     -2.81229792    -28.30000000     60.11693332      0.00000000
   //8
   //1    -84.28224914     -0.48406407   -130.11314158     -1.80061962     -0.00914512     -2.82881943    -28.30000000     57.90240833      0.00000000
   //9
   //1     65.13972189     16.85381518    -87.92399172      1.42598734      0.37691999     -1.92154637    -28.30000000     30.32639124      0.00000000
   //10
   //
   //    2  1001.0000       0       8    -32.98597627      2.21586178 ==> next event
   auto res = ReadLineAndCheck(10, " ");
   switch (res) {
      case KVFileReader::ReadStatus::OK:
         {
            int nalpha = GetIntReadPar(0);
            nuc->SetZ(nalpha*2);
            nuc->SetA(nalpha*4);
            nuc->SetPosition(GetDoubleReadPar(1),GetDoubleReadPar(2),GetDoubleReadPar(3));
            nuc->SetMomentum(nalpha*KVNucleus::hbar*TVector3(GetDoubleReadPar(4),GetDoubleReadPar(5),GetDoubleReadPar(6)));
            nuc->GetParameters()->SetValue("Total internal energy [MeV]", GetDoubleReadPar(7));
            nuc->GetParameters()->SetValue("Kinetic CM energy [MeV]", GetDoubleReadPar(8));
            nuc->GetParameters()->SetValue("Internal kinetic energy [MeV/alpha]", GetDoubleReadPar(9));
            auto reres = ReadLineAndCheck(1,"|");//not used as delimiter: KVString::Begin() just returns the full string
            if(reres!=KVFileReader::ReadStatus::OK)
            {
               KVError::Error(this,"ReadNucleus","ReadLineAndCheck(1) status = %s",GetReadStatus(reres).Data());
               return kFALSE;
            }
            KVNumberList idlist(GetReadPar(0));
            nuc->GetParameters()->SetValue("Alpha id list", idlist.AsString());
            break;
         }
      default:
         KVError::Error(this,"ReadNucleus","ReadLineAndCheck(10) status = %s",GetReadStatus(res).Data());
         return kFALSE;
   }
   return kTRUE;
}
 
 
ClassImp(KVSimReader_HaC)