fPHOSEta(0.13),
fEMCALMinPhi(79.),
fEMCALMaxPhi(191.),
- fEMCALEta(0.71)
-
+ fEMCALEta(0.71),
+ fItune(-1)
{
// Default Constructor
fEnergyCMS = 5500.;
fPHOSEta(0.13),
fEMCALMinPhi(79.),
fEMCALMaxPhi(191.),
- fEMCALEta(0.71)
+ fEMCALEta(0.71),
+ fItune(-1)
{
// default charm production at 5. 5 TeV
// semimuonic decay
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
//Used in some processes to selected child properties
Bool_t CheckKinematicsOnChild();
void GetSubEventTime();
-
+ virtual void SetTune(Int_t itune) {fItune = itune;}
protected:
// adjust the weight from kinematic cuts
void AdjustWeights() const;
Float_t fEMCALMinPhi; // Minimum phi EMCAL
Float_t fEMCALMaxPhi; // Maximum phi EMCAL
Float_t fEMCALEta; // Maximum eta EMCAL
-
+ Int_t fItune; // Pythia tune
private:
AliGenPythiaPlus(const AliGenPythiaPlus &Pythia);
AliGenPythiaPlus & operator=(const AliGenPythiaPlus & rhs);
pythia.Copy(*this);
}
-void AliPythia6::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfunc)
+void AliPythia6::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfunc, Int_t /*tune*/)
{
// Initialise the process to generate
if (!AliPythiaRndm::GetPythiaRandom())
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 tune);
virtual void GenerateEvent() {Pyevnt();}
virtual void GenerateMIEvent() {Pyevnw();}
virtual void HadronizeEvent() {Pyexec();}
// Convert to compressed code and print result (for debugging only)
virtual Int_t CheckedLuComp(Int_t kf) = 0;
// Pythia initialisation for selected processes
- virtual void ProcInit (Process_t process, Float_t energy, StrucFunc_t strucfunc) = 0;
+ virtual void ProcInit (Process_t process, Float_t energy, StrucFunc_t strucfunc, Int_t tune) = 0;
virtual void GenerateEvent() = 0;
virtual void GenerateMIEvent() = 0;
virtual Int_t GetNumberOfParticles() = 0;
pythia.Copy(*this);
}
-void AliPythia8::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfunc)
+void AliPythia8::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfunc, Int_t tune)
{
// Initialise the process to generate
if (!AliPythiaRndm::GetPythiaRandom())
// Multiple interactions on.
ReadString("PartonLevel:MI = on");
// Double Gaussian matter distribution.
- ReadString("MultipleInteractions:bProfile = 2");
- ReadString("MultipleInteractions:coreFraction = 0.5");
- ReadString("MultipleInteractions:coreRadius = 0.4");
+ ReadString("MultipartonInteractions:bProfile = 2");
+ ReadString("MultipartonInteractions:coreFraction = 0.5");
+ ReadString("MultipartonInteractions:coreRadius = 0.4");
// pT0.
- ReadString("MultipleInteractions:pTmin = 2.0");
+ ReadString("MultipartonInteractions:pTmin = 2.0");
// Reference energy for pT0 and energy rescaling pace.
- ReadString("MultipleInteractions:ecmRef = 1800.");
- ReadString("MultipleInteractions:ecmPow = 0.25");
+ ReadString("MultipartonInteractions:ecmRef = 1800.");
+ ReadString("MultipartonInteractions:ecmPow = 0.25");
// String drawing almost completely minimizes string length.
// SetPARP(85,0.9);
// SetPARP(86,0.95);
ReadString("SoftQCD:minBias = on");
ReadString("SoftQCD:singleDiffractive = on");
ReadString("SoftQCD:doubleDiffractive = on");
- AtlasTuning();
+ if (tune == -1) AtlasTuning();
break;
case kPyMbDefault:
// Minimum Bias pp-Collisions
ReadString("SoftQCD:minBias = on");
ReadString("SoftQCD:singleDiffractive = on");
ReadString("SoftQCD:doubleDiffractive = on");
+ ReadString("SoftQCD:doubleDiffractive = on");
+ if (tune > -1) ReadString(Form("Tune:pp = %3d", tune));
break;
case kPyLhwgMb:
// Les Houches Working Group 05 Minimum Bias pp-Collisions: hep-ph/0604120
// SetPARP(85,0.33); // Regulates gluon prod. mechanism
// SetPARP(86,0.66); // Regulates gluon prod. mechanism
ReadString("SigmaProcess:factorMultFac = 1.");
+
}
void AliPythia8::SetPtHardRange(Float_t ptmin, Float_t ptmax)
// convert to compressed code and print result (for debugging only)
virtual Int_t CheckedLuComp(Int_t /*kf*/) {return -1;}
// Pythia initialisation for selected processes
- virtual void ProcInit (Process_t process, Float_t energy, StrucFunc_t strucfunc);
+ virtual void ProcInit (Process_t process, Float_t energy, StrucFunc_t strucfunc, Int_t tune);
virtual void GenerateEvent();
virtual void GenerateMIEvent();
virtual void HadronizeEvent();