break;
case kPyMb:
case kPyMbNonDiffr:
+ case kPyMbMSEL1:
case kPyJets:
case kPyDirectGamma:
break;
//
// Configure detector (EMCAL like)
//
- fPythia->SetPARU(51, fPycellEtaMax);
- fPythia->SetMSTU(51, fPycellNEta);
- fPythia->SetMSTU(52, fPycellNPhi);
+ fPythia->SetPARU(51, fPycellEtaMax);
+ fPythia->SetMSTU(51, fPycellNEta);
+ fPythia->SetMSTU(52, fPycellNPhi);
//
// Configure Jet Finder
//
- fPythia->SetPARU(58, fPycellThreshold);
- fPythia->SetPARU(52, fPycellEtSeed);
- fPythia->SetPARU(53, fPycellMinEtJet);
- fPythia->SetPARU(54, fPycellMaxRadius);
- fPythia->SetMSTU(54, 2);
+ fPythia->SetPARU(58, fPycellThreshold);
+ fPythia->SetPARU(52, fPycellEtSeed);
+ fPythia->SetPARU(53, fPycellMinEtJet);
+ fPythia->SetPARU(54, fPycellMaxRadius);
+ fPythia->SetMSTU(54, 2);
//
// This counts the total number of calls to Pyevnt() per run.
fTrialsRun = 0;
fDyBoost = 0;
Warning("Init","SetNuclei used. Use SetProjectile + SetTarget instead. fDyBoost has been reset to 0\n");
}
-
+
if (fQuench) {
fPythia->InitQuenching(0., 0.1, 0.6e6, 0);
}
-
}
void AliGenPythia::Generate()
if (fProcess != kPyMb && fProcess != kPyJets &&
fProcess != kPyDirectGamma &&
fProcess != kPyMbNonDiffr &&
+ fProcess != kPyMbMSEL1 &&
fProcess != kPyW && fProcess != kPyZ &&
fProcess != kPyCharmppMNRwmi && fProcess != kPyBeautyppMNRwmi) {
SetMSEL(0);
SetMSUB(95,1); // low pt production
+ AtlasTuning();
+ break;
+ case kPyMbMSEL1:
+ ConfigHeavyFlavor();
+// Intrinsic <kT^2>
+ SetMSTP(91,1);// Width (1=gaussian) primordial kT dist. inside hadrons
+ SetPARP(91,1.); // <kT^2> = PARP(91,1.)^2
+ SetPARP(93,5.); // Upper cut-off
+// Set Q-quark mass
+ SetPMAS(4,1,1.2); // Charm quark mass
+ SetPMAS(5,1,4.78); // Beauty quark mass
+ SetPARP(71,4.); // Defaut value
+// Atlas Tuning
AtlasTuning();
break;
case kPyJets:
kPyCharmPbPbMNR, kPyD0PbPbMNR, kPyDPlusPbPbMNR, kPyDPlusStrangePbPbMNR, kPyBeautyPbPbMNR,
kPyCharmpPbMNR, kPyD0pPbMNR, kPyDPluspPbMNR, kPyDPlusStrangepPbMNR, kPyBeautypPbMNR,
kPyCharmppMNR, kPyCharmppMNRwmi, kPyD0ppMNR, kPyDPlusppMNR, kPyDPlusStrangeppMNR,
- kPyBeautyppMNR, kPyBeautyppMNRwmi, kPyW, kPyZ,
+ kPyBeautyppMNR, kPyBeautyppMNRwmi, kPyW, kPyZ, kPyMbMSEL1,
kPyOldUEQ2ordered, kPyOldUEQ2ordered2, kPyOldPopcorn}
Process_t;
/*