AliGenPythia::AliGenPythia():
AliGenMC(),
fProcess(kPyCharm),
+ fItune(-1),
fStrucFunc(kCTEQ5L),
fKineBias(0.),
fTrials(0),
AliGenPythia::AliGenPythia(Int_t npart)
:AliGenMC(npart),
fProcess(kPyCharm),
+ fItune(-1),
fStrucFunc(kCTEQ5L),
fKineBias(0.),
fTrials(0),
fRL = 0x0;
}
//
- fPythia->ProcInit(fProcess,fEnergyCMS,fStrucFunc);
+ fPythia->ProcInit(fProcess,fEnergyCMS,fStrucFunc, fItune);
// Forward Paramters to the AliPythia object
fDecayer->SetForceDecay(fForceDecay);
// Switch off Heavy Flavors on request
void MakeHeader();
void GeneratePileup();
Process_t fProcess; //Process type
+ Int_t fItune; // Pythia tune > 6.4
StrucFunc_t fStrucFunc; //Structure Function
Float_t fKineBias; //!Bias from kinematic selection
Int_t fTrials; //!Number of trials for current event
pythia.Copy(*this);
}
-void AliPythia::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfunc)
+void AliPythia::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfunc, Int_t itune)
{
// Initialise the process to generate
if (!AliPythiaRndm::GetPythiaRandom())
}
//
// Initialize PYTHIA
+//
+// Select the tune
+ if (itune > -1) Pytune(itune);
+
+//
SetMSTP(41,1); // all resonance decays switched on
Initialize("CMS","p","p",fEcms);
virtual Int_t CheckedLuComp(Int_t kf);
// Pythia initialisation for selected processes
virtual void ProcInit
- (Process_t process, Float_t energy, StrucFunc_t strucfunc);
+ (Process_t process, Float_t energy, StrucFunc_t strucfunc, Int_t itune = -1);
// treat protons as inside nuclei
virtual void SetNuclei(Int_t a1, Int_t a2, Int_t pdf);
// Print particle properties