#include "AliGenHijing.h"
#include "AliGenHijingEventHeader.h"
-#include "AliRun.h"
#include "AliHijingRndm.h"
+#include "AliLog.h"
+#include "AliRun.h"
ClassImp(AliGenHijing)
fNoHeavyQuarks(kFALSE)
{
// Constructor
- SetEnergyCMS(5500.);
+ fEnergyCMS = 5500.;
AliHijingRndm::SetHijingRandom(GetRandom());
}
{
// Default PbPb collisions at 5. 5 TeV
//
- SetEnergyCMS(5500.);
+ fEnergyCMS = 5500.;
fName = "Hijing";
fTitle= "Particle Generator using HIJING";
//
fTrials = 0;
for (j = 0;j < 3; j++) origin0[j] = fOrigin[j];
+
if(fVertexSmear == kPerEvent) {
Vertex();
for (j=0; j < 3; j++) origin0[j] = fVertex[j];
Float_t sign = (fRandomPz && (Rndm() < 0.5))? -1. : 1.;
+
while(1)
{
// Generate one event
Int_t np = fParticles.GetEntriesFast();
- printf("\n **************************************************%d\n",np);
Int_t nc = 0;
if (np == 0 ) continue;
Int_t i;
origin[0] = origin0[0]+iparticle->Vx()/10;
origin[1] = origin0[1]+iparticle->Vy()/10;
origin[2] = origin0[2]+iparticle->Vz()/10;
- tof = kconv * iparticle->T() + sign * origin0[2] / 3.e10;
- if (fPileUpTimeWindow > 0.) tof += tInt;
-
+ if (TestBit(kVertexRange)) {
+ tof = kconv * iparticle->T() + sign * origin0[2] / 2.99792458e10;
+ } else {
+ tof = kconv * iparticle->T();
+ if (fPileUpTimeWindow > 0.) tof += tInt;
+ }
imo = -1;
TParticle* mother = 0;
if (hasMother) {
delete[] newPos;
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) {
jev += nc;
if (jev >= fNpart || fNpart == -1) {
fKineBias = Float_t(fNpart)/Float_t(fTrials);
- printf("\n Trials: %i %i %i\n",fTrials, fNpart, jev);
+ AliInfo(Form("\n Trials: %i %i %i\n",fTrials, fNpart, jev));
break;
}
}
((AliGenHijingEventHeader*) header)->SetSpectators(fProjectileSpecn, fProjectileSpecp,
fTargetSpecn,fTargetSpecp);
((AliGenHijingEventHeader*) header)->SetReactionPlaneAngle(fHijing->GetHINT1(20));
+// printf("Impact Parameter %13.3f \n", fHijing->GetHINT1(19));