fLHC(kFALSE),
fRandomPz(kFALSE),
fNoHeavyQuarks(kFALSE),
- fEventTime(0.),
fIsoft(4),
fNtMax(150),
fIpop(1),
fLHC(kFALSE),
fRandomPz(kFALSE),
fNoHeavyQuarks(kFALSE),
- fEventTime(0.),
fIsoft(1),
fNtMax(150),
fIpop(1),
Float_t polar[3] = {0,0,0};
Float_t origin[3] = {0,0,0};
Float_t origin0[3] = {0,0,0};
+ Float_t time0 = 0.;
Float_t p[3];
Float_t tof;
fTrials = 0;
for (j = 0;j < 3; j++)
origin0[j] = fOrigin[j];
+ time0 = fTimeOrigin;
if(fVertexSmear == kPerEvent) {
Vertex();
for (j=0; j < 3; j++)
origin0[j] = fVertex[j];
+ time0 = fTime;
}
Float_t sign = (fRandomPz && (Rndm() < 0.5))? -1. : 1.;
fVertex[0] = origin0[0];
fVertex[1] = origin0[1];
fVertex[2] = origin0[2];
+ fTime = time0;
// First select parent particles
for (Int_t i = 0; i < np; i++) {
} // selected
} // particle loop final state
- //Time of the interactions
- Float_t tInt = 0.;
- if (fPileUpTimeWindow > 0.)
- tInt = fPileUpTimeWindow * (2. * gRandom->Rndm() - 1.);
-
// Write particles to stack
for (Int_t i = 0; i<np; i++) {
TParticle *iparticle = (TParticle *) fParticles.At(i);
origin[0] = origin0[0]+iparticle->Vx()/10;
origin[1] = origin0[1]+iparticle->Vy()/10;
origin[2] = origin0[2]+iparticle->Vz()/10;
- fEventTime = 0.;
-
- if (TestBit(kVertexRange)) {
- fEventTime = sign * origin0[2] / 2.99792458e10;
- tof = kconv * iparticle->T() + fEventTime;
- } else {
- tof = kconv * iparticle->T();
- fEventTime = tInt;
- if (fPileUpTimeWindow > 0.) tof += tInt;
- }
+ tof = time0+kconv * iparticle->T();
+
imo = -1;
TParticle* mother = 0;
if (hasMother) {
fHeader->SetTrials(fTrials);
// Event Vertex
fHeader->SetPrimaryVertex(fVertex);
- fHeader->SetInteractionTime(fEventTime);
+ fHeader->SetInteractionTime(fTime);
fCollisionGeometry = fHeader;
AddHeader(fHeader);
Int_t fLHC; // Assume LHC as lab frame
Bool_t fRandomPz; // Randomise sign of pz event by event
Bool_t fNoHeavyQuarks; // If true no heavy quarks are produced
- Float_t fEventTime; // The event time
Int_t fIsoft; // ISOFT (D=1): select Default AMPT or String Melting
Int_t fNtMax; // NTMAX: number of timesteps (D=150)
Int_t fIpop; // (D=1,yes;0,no) flag for popcorn mechanism(netbaryon stopping)
// check if stable
Bool_t Stable(TParticle* particle) const;
- ClassDef(AliGenAmpt, 3) // AliGenerator interface to Ampt
+ ClassDef(AliGenAmpt, 4) // AliGenerator interface to Ampt
};
#endif