SetMC(AliPythia::Instance());
fPythia=(AliPythia*) fMCEvGen;
+ printf("Pythia Init %p", fPythia);
+
//
fParentWeight=1./Float_t(fNpart);
//
if (fQuench) {
fPythia->SetMSTJ(1, 0);
}
-
+ printf("Pythia pyevnt %p", fPythia);
fPythia->Pyevnt();
if (fQuench) {
fPythia->Pyexec();
}
+ if (gAlice) {
+ if (gAlice->GetEvNumber()>=fDebugEventFirst &&
+ gAlice->GetEvNumber()<=fDebugEventLast) fPythia->Pylist(1);
+ }
- if (gAlice->GetEvNumber()>=fDebugEventFirst &&
- gAlice->GetEvNumber()<=fDebugEventLast) fPythia->Pylist(1);
fTrials++;
-
+ printf("Pythia import t %p %p ", fPythia, fParticles);
fPythia->ImportParticles(fParticles,"All");
Boost();
{
// Adjust the weights after generation of all events
//
- TParticle *part;
- Int_t ntrack=gAlice->GetMCApp()->GetNtrack();
- for (Int_t i=0; i<ntrack; i++) {
- part= gAlice->GetMCApp()->Particle(i);
- part->SetWeight(part->GetWeight()*fKineBias);
+ if (gAlice) {
+ TParticle *part;
+ Int_t ntrack=gAlice->GetMCApp()->GetNtrack();
+ for (Int_t i=0; i<ntrack; i++) {
+ part= gAlice->GetMCApp()->Particle(i);
+ part->SetWeight(part->GetWeight()*fKineBias);
+ }
}
}
jets[3][i]);
}
}
- gAlice->SetGenEventHeader(fHeader);
+ if (gAlice) gAlice->SetGenEventHeader(fHeader);
}