#include <TClonesArray.h>
#include "AliGenHydjet.h"
+#include "AliLog.h"
#include "AliGenHydjetEventHeader.h"
#include "AliRun.h"
#include "AliPythiaRndm.h"
AliGenHydjet::AliGenHydjet(Int_t npart) :
AliGenMC(npart),
//initial parameters
- fEnergyCMS(5500.), //Energy cms
fFrame("CMS"), // Reference frame
fAtomicWeigth(207), // Projectile-Target atomic weight
fIfbtype(0), // centrality type
{
// Default PbPb collisions at 5. 5 TeV
//
+ fEnergyCMS = 5500.; //Energy cms
fName = "Hydjet";
fTitle = "Particle Generator using Hydjet";
// Set random number generator
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;
kf = 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;
}
while(1)
fHydjet->ImportParticles(&fParticles,"All");
Int_t np = fParticles.GetEntriesFast();
- printf("\n **************************************************%d\n",np);
Int_t nc = 0;
if (np == 0 ) continue;
Int_t i;
fVertex[0] = origin0[0];
fVertex[1] = origin0[1];
fVertex[2] = origin0[2];
+ fTime = time0;
//
// Now select the final state particles
origin[0] = fVertex[0]+iparticle->Vx()/10;
origin[1] = fVertex[1]+iparticle->Vy()/10;
origin[2] = fVertex[2]+iparticle->Vz()/10;
- tof = kconv*iparticle->T();
+ tof = fTime + kconv*iparticle->T();
imo = -1;
//imo = iparticle->GetFirstMother();
} // particle loop
delete[] pSelected;
- printf("\n I've put %i particles on the stack \n",nc);
+ AliInfo(Form("\n I've put %i particles on the stack \n",nc));
if (nc > 0) break;
} // event loop
MakeHeader();
// Event Vertex
header->SetPrimaryVertex(fVertex);
+ header->SetInteractionTime(fTime);
AddHeader(header);
fCollisionGeometry = (AliGenHydjetEventHeader*) header;
}