Float_t polar[3] = {0,0,0};
Float_t origin[3] = {0,0,0};
+ Float_t time = 0.;
Float_t origin0[3] = {0,0,0};
+ Float_t time0 = 0.;
Float_t p[3];
Float_t mass, energy;
Vertex();
for (j=0; j < 3; j++) origin0[j] = fVertex[j];
+ time0 = fTime;
// Generate one event
origin[0] = vrho*TMath::Cos(vphi + evrot);
origin[1] = vrho*TMath::Sin(vphi + evrot);
origin[2] = iparticle->Vz();
-
+ time = iparticle->T();
+
imo = -1;
- TParticle* mother = 0;
+ // TParticle* mother = 0;
if (hasMother) {
imo = iparticle->GetFirstMother();
- mother = (TParticle *) fParticles.At(imo);
+ // mother = (TParticle *) fParticles.At(imo);
} // if has mother
Bool_t tFlag = (!hasDaughter);
// printf("Pushing Track %d with status %d mother %d\n", kf, tFlag, imo>=0?idsOnStack[imo]:imo);
PushTrack(tFlag,imo>=0?idsOnStack[imo]:imo,kf,
p[0],p[1],p[2],energy,
- origin[0],origin[1],origin[2],iparticle->T(),
+ origin[0],origin[1],origin[2],time,
polar[0],polar[1],polar[2],
hasMother ? kPDecay:kPNoProcess,nt);
idsOnStack[i] = nt;
origin[0] = vrho*TMath::Cos(vphi + evrot);
origin[1] = vrho*TMath::Sin(vphi + evrot);
origin[2] = iparticle->Vz();
-
+ time = iparticle->T();
+
imo = -1;
- TParticle* mother = 0;
+ // TParticle* mother = 0;
if (hasMother) {
imo = iparticle->GetFirstMother();
- mother = (TParticle *) fParticles.At(imo);
+ // mother = (TParticle *) fParticles.At(imo);
} // if has mother
Bool_t tFlag = (hasDaughter);
// printf("Pushing Track %d with status %d mother %d\n", kf, tFlag, imo>=0?idsOnStack[imo]:imo);
PushTrack(tFlag,imo>=0?idsOnStack[imo]:imo,kf,
p[0],p[1],p[2],energy,
- origin[0],origin[1],origin[2],iparticle->T(),
+ origin[0],origin[1],origin[2],time,
polar[0],polar[1],polar[2],
hasMother ? kPDecay:kPNoProcess,nt);
idsOnStack[i] = nt;
origin[0] = origin0[0]+vrho*TMath::Cos(vphi + evrot);
origin[1] = origin0[1]+vrho*TMath::Sin(vphi + evrot);
origin[2] = origin0[2]+iparticle->Vz();
+ time = time0+iparticle->T();
imo = nt;
// mother = (TParticle *) fParticles.At(nt);
// printf("Pushing Track %d with status %d mother %d\n", kf, tFlag, imo);
PushTrack(tFlag,imo,kf,
p[0],p[1],p[2],energy,
- origin[0],origin[1],origin[2],iparticle->T(),
+ origin[0],origin[1],origin[2],time,
polar[0],polar[1],polar[2],
hasMother ? kPDecay:kPNoProcess,nt);
fNprimaries++;
eventVertex[0] = origin0[0];
eventVertex[1] = origin0[1];
eventVertex[2] = origin0[2];
+ Float_t eventTime = time0;
// Builds the event header, to be called after each event
AliGenEventHeader* header = new AliGenHijingEventHeader("Therminator");
((AliGenHijingEventHeader*) header)->SetNProduced(fNprimaries);
((AliGenHijingEventHeader*) header)->SetPrimaryVertex(eventVertex);
+ ((AliGenHijingEventHeader*) header)->SetInteractionTime(eventTime);
((AliGenHijingEventHeader*) header)->SetImpactParameter(0.0);
((AliGenHijingEventHeader*) header)->SetTotalEnergy(0.0);
((AliGenHijingEventHeader*) header)->SetHardScatters(0);
// ((AliGenHijingEventHeader*) header)->SetTrials(fTrials);
// Event Vertex
header->SetPrimaryVertex(fVertex);
+ header->SetInteractionTime(fTime);
AddHeader(header);
fCollisionGeometry = (AliGenHijingEventHeader*) header;
- delete idsOnStack;
+ delete [] idsOnStack;
// gAlice->SetGenEventHeader(header);
}
fFOHSlocation.Append(aroot);
fFOHSlocation.Append("/TTherminator/data/LHC500C4050");
}
+ else if (strstr(set, "LHC276TC0005")) {
+ AliWarning(Form("AliGenTherminator: Selected default Lhyquid hypersurface"));
+ AliWarning(Form(" Pb-Pb collisions, 2.76TeV, centrality 00-05 percent"));
+ AliWarning(Form(" freeze-out criteria Tf=145 MeV"));
+ AliWarning(Form(" for details see $(ALICE_ROOT)/TTherminator/data/LHC276TC0005/FO.txt"));
+ fFOHSlocation.Append(aroot);
+ fFOHSlocation.Append("/TTherminator/data/LHC276TC0005");
+ }
+ else if (strstr(set, "LHC276TC0510")) {
+ AliWarning(Form("AliGenTherminator: Selected default Lhyquid hypersurface"));
+ AliWarning(Form(" Pb-Pb collisions, 2.76TeV, centrality 05-10 percent"));
+ AliWarning(Form(" freeze-out criteria Tf=145 MeV"));
+ AliWarning(Form(" for details see $(ALICE_ROOT)/TTherminator/data/LHC276TC0510/FO.txt"));
+ fFOHSlocation.Append(aroot);
+ fFOHSlocation.Append("/TTherminator/data/LHC276TC0510");
+ }
+ else if (strstr(set, "LHC276TC1020")) {
+ AliWarning(Form("AliGenTherminator: Selected default Lhyquid hypersurface"));
+ AliWarning(Form(" Pb-Pb collisions, 2.76TeV, centrality 10-20 percent"));
+ AliWarning(Form(" freeze-out criteria Tf=145 MeV"));
+ AliWarning(Form(" for details see $(ALICE_ROOT)/TTherminator/data/LHC276TC1020/FO.txt"));
+ fFOHSlocation.Append(aroot);
+ fFOHSlocation.Append("/TTherminator/data/LHC276TC1020");
+ }
+ else if (strstr(set, "LHC276TC2030")) {
+ AliWarning(Form("AliGenTherminator: Selected default Lhyquid hypersurface"));
+ AliWarning(Form(" Pb-Pb collisions, 2.76TeV, centrality 20-30 percent"));
+ AliWarning(Form(" freeze-out criteria Tf=145 MeV"));
+ AliWarning(Form(" for details see $(ALICE_ROOT)/TTherminator/data/LHC276TC2030/FO.txt"));
+ fFOHSlocation.Append(aroot);
+ fFOHSlocation.Append("/TTherminator/data/LHC276TC2030");
+ }
+ else if (strstr(set, "LHC276TC3040")) {
+ AliWarning(Form("AliGenTherminator: Selected default Lhyquid hypersurface"));
+ AliWarning(Form(" Pb-Pb collisions, 2.76TeV, centrality 30-40 percent"));
+ AliWarning(Form(" freeze-out criteria Tf=145 MeV"));
+ AliWarning(Form(" for details see $(ALICE_ROOT)/TTherminator/data/LHC276TC3040/FO.txt"));
+ fFOHSlocation.Append(aroot);
+ fFOHSlocation.Append("/TTherminator/data/LHC276TC3040");
+ }
+ else if (strstr(set, "LHC276TC4050")) {
+ AliWarning(Form("AliGenTherminator: Selected default Lhyquid hypersurface"));
+ AliWarning(Form(" Pb-Pb collisions, 2.76TeV, centrality 40-50 percent"));
+ AliWarning(Form(" freeze-out criteria Tf=145 MeV"));
+ AliWarning(Form(" for details see $(ALICE_ROOT)/TTherminator/data/LHC276TC4050/FO.txt"));
+ fFOHSlocation.Append(aroot);
+ fFOHSlocation.Append("/TTherminator/data/LHC276TC4050");
+ }
else {
AliWarning(Form("Did not find Lhyquid set %s", set));
AliWarning(Form("Reverting to default: current directory"));