Bool_t KVSimReader_SMF_asym::ReadEvent(){
evt->Clear();
Int_t mult=0;
Int_t res = ReadLineAndCheck(1," ");
switch (res){
case 0:
return kFALSE;
case 1:
mult = GetIntReadPar(0);
break;
default:
return kFALSE;
}
evt->SetNumber(nv->GetIntValue("event_number"));
evt->GetParameters()->SetValue("sub_number",ndes);
for (Int_t mm=0; mm<mult; mm+=1){
nuc = (KVSimNucleus* )evt->AddParticle();
if (!ReadNucleus()) return kFALSE;
}
nevt+=1;
ndes+=1;
return kTRUE;
}
Bool_t KVSimReader_MMM_asym::ReadEvent(){
evt->Clear();
Int_t mult=0,natt=0;
ReadLine(" ");
Int_t res = GetNparRead();
switch (res){
case 0:
return kFALSE;
default:
//nlus = toks->GetEntries();
idx = 0;
mult = GetIntReadPar(idx++);
natt = 5*mult+1;
if (natt!=res){
Info("ReadEvent","Nombre de parametres (%d) different de celui attendu (%d)",res,natt);
return kFALSE;
}
for (Int_t mm=0; mm<mult; mm+=1){
nuc = (KVSimNucleus* )evt->AddParticle();
ReadNucleus();
}
evt->SetNumber(nevt);
nevt+=1;
return kTRUE;
}
}
Bool_t KVSimReader_ELIE::ReadEvent()
{
// Events begin with event number, multiplicity and (?) impact parameter,
// followed by list of particles:
// 0 61 0
// 0 11 24.0 57.4156777729 306.493830835 107.248038513
// 1 2 3.0 44.7736647511 50.6065561553 127.526000603
// 2 7 15.0 57.7807658304 121.672537111 55.9601768553
// 3 9 19.0 35.5421741719 70.5034488937 236.567079783
// etc.
evt->Clear();
Int_t res = ReadLineAndCheck(3," ");
Int_t mult=0;
switch (res){
case 0:
Info("ReadEvent","case 0 line est vide");
return kFALSE;
case 1:
evt->SetNumber(GetIntReadPar(0));
mult=GetIntReadPar(1);
evt->GetParameters()->SetValue("b",GetDoubleReadPar(2));
evt->GetParameters()->SetValue("mult",mult);
break;
default:
return kFALSE;
}
for (Int_t mm=0; mm<mult; mm++){
nuc = (KVSimNucleus* )evt->AddParticle();
if (!ReadNucleus()) return kFALSE;
}
nevt++;
return kTRUE;
}
Bool_t KVSimReader_MMM::ReadEvent()
{
evt->Clear();
Int_t mult = 0, natt = 0;
ReadLine(" ");
Int_t res = GetNparRead();
switch (res) {
case 0:
Info("ReadEvent", "case 0 line est vide");
return kFALSE;
default:
idx = 0;
mult = GetIntReadPar(idx++);
natt = 6 * mult + 4 + 1;
if (natt != res) {
Info("ReadEvent", "Nombre de parametres (%d) different de celui attendu (%d)", res, natt);
return kFALSE;
}
evt->GetParameters()->SetValue("Eth", GetDoubleReadPar(idx++));
evt->GetParameters()->SetValue("Eclb", GetDoubleReadPar(idx++));
evt->GetParameters()->SetValue("Eexci", GetDoubleReadPar(idx++));
evt->GetParameters()->SetValue("Ebind", GetDoubleReadPar(idx++));
for (Int_t mm = 0; mm < mult; mm += 1) {
nuc = (KVSimNucleus*)evt->AddParticle();
ReadNucleus();
}
evt->SetNumber(nevt);
nevt += 1;
return kTRUE;
}
}
