fProcess(kPyMb),
fEcms(0.),
fStrucFunc(kCTEQ5L),
+ fCellJet(),
fEtSeed(0.),
fMinEtJet(0.),
fRJet(0.),
+ fClusterJet(),
fYScale(0.),
fPtScale(0.),
fNJetMin(0),
fProcess(kPyMb),
fEcms(0.),
fStrucFunc(kCTEQ5L),
+ fCellJet(),
fEtSeed(0.),
fMinEtJet(0.),
fRJet(0.),
+ fClusterJet(),
fYScale(0.),
fPtScale(0.),
fNJetMin(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
// With parton showers on we are generating "Z inclusive process"
ReadString("PartonLevel:ISR = on");
ReadString("PartonLevel:FSR = on");
- break;
+ break;
+ case kPyZgamma:
+ //Inclusive production of Z/gamma*
+ //f fbar -> Z/gamma
+ ReadString("WeakSingleBoson:ffbar2gmZ = on");
+ // Initial/final parton shower on (Pythia default)
+ // With parton showers on we are generating "Z inclusive process"
+ ReadString("PartonLevel:ISR = on");
+ ReadString("PartonLevel:FSR = on");
+ break;
+ case kPyMBRSingleDiffraction:
+ ReadString("Diffraction:PomFlux = 5");
+ ReadString("SoftQCD:singleDiffractive = on");
+ break;
+ case kPyMBRDoubleDiffraction:
+ ReadString("Diffraction:PomFlux = 5");
+ ReadString("SoftQCD:doubleDiffractive = on");
+ break;
+ case kPyMBRCentralDiffraction:
+ ReadString("Diffraction:PomFlux = 5");
+ ReadString("SoftQCD:centralDiffractive = on");
+ break;
case kPyMbWithDirectPhoton:
case kPyBeautyJets:
case kPyMbAtlasTuneMC09:
// 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)
void AliPythia8::SetIntrinsicKt(Float_t kt)
{
+// Set the intrinsic kt
ReadString("BeamRemnants:primordialKT = on");
ReadString("BeamRemnants:primordialKTsoft = 0.");
ReadString(Form("BeamRemnants:primordialKThard = %13.3f", kt));