AliGenerator* MbPythiaTunePerugia0bele()
{
comment = comment.Append(" pp: Pythia (Perugia0) bele (1 bbbar per event, 1 b-hadron in |y|<1, 1 electron in |y|<2");
//
// Pythia
AliGenPythia* pythia = new AliGenPythia(-1);
pythia->SetMomentumRange(0, 999999.);
pythia->SetThetaRange(0., 180.);
pythia->SetYRange(-2.,2.);
pythia->SetPtRange(0,1000.);
pythia->SetProcess(kPyBeautyppMNRwmi);
pythia->SetEnergyCMS(energy);
// Tune
// 320 Perugia 0
pythia->SetTune(320);
pythia->UseNewMultipleInteractionsScenario();
//
// decays
pythia->SetCutOnChild(1);
pythia->SetPdgCodeParticleforAcceptanceCut(11);
pythia->SetChildYRange(-2,2);
pythia->SetChildPtRange(0,10000.);
return pythia;
}
AliGenerator* MbPythiaTunePerugia0BtoJpsi2mu() {
comment = comment.Append(" pp: Pythia (Perugia0) BtoJpsi (1 bbbar per event, 1 b-hadron, 1 J/psi");
// Pythia
AliGenPythia* pythia = new AliGenPythia(-1);
pythia->SetMomentumRange(0, 999999.);
pythia->SetYRange(-4.5, -2.0);
pythia->SetPtRange(0,1000.);
pythia->SetProcess(kPyBeautyppMNRwmi);
pythia->SetEnergyCMS(energy);
// Tune
// 320 Perugia 0
pythia->SetTune(320);
pythia->UseNewMultipleInteractionsScenario();
//
// decays
pythia->SetCutOnChild(1);
pythia->SetPdgCodeParticleforAcceptanceCut(443);
pythia->SetChildYRange(-4.5, -2.0);
pythia->SetChildPtRange(0,10000.);
//
// decays
pythia->SetForceDecay(kBJpsiDiMuon);
return pythia;
}
AliGenerator* CreateGenerator(Tune_t tune, Float_t energy)
{
if (tune == -1) {
// phojet
AliGenDPMjet* gener = new AliGenDPMjet(1);
gener->SetProcess(kDpmMb);
gener->SetProjectile("P", 1, 1);
gener->SetTarget("P", 1, 1);
gener->SetEnergyCMS(energy);
return gener;
}
if (tune != kPyTuneAtlasCSC && tune != kPyTuneCDFA && tune != kPyTuneCMS6D6T && tune != kPyTunePerugia0) {
Printf("Unknown pythia tune, quitting");
exit(1);
} else {
AliGenPythia * gener = new AliGenPythia(1);
//
// set pythia tune
gener->SetTune(tune);
// structure function
if(tune == kPyTuneAtlasCSC) {
cout << "Setting structure function" << endl;
gener->SetStrucFunc(kCTEQ61);
}
if(tune == kPyTuneCMS6D6T) {
cout << "Setting structure function" << endl;
gener->SetStrucFunc(kCTEQ6l);
}
if(tune == kPyTunePerugia0) {
cout << "Setting new parton shower" << endl;
gener->UseNewMultipleInteractionsScenario();
}
// charm, beauty, charm_unforced, beauty_unforced, jpsi, jpsi_chi, mb
gener->SetProcess(kPyMb);
// Centre of mass energy
gener->SetEnergyCMS(energy);
// Set target/projectile // Is this needded?
gener->SetProjectile("P", 1, 1);
gener->SetTarget("P", 1, 1);
return gener;
}
}
AliGenerator* MbPythiaTunePerugia0()
{
comment = comment.Append(" pp: Pythia low-pt (Perugia0)");
//
// Pythia
AliGenPythia* pythia = new AliGenPythia(-1);
pythia->SetMomentumRange(0, 999999.);
pythia->SetThetaRange(0., 180.);
pythia->SetYRange(-12.,12.);
pythia->SetPtRange(0,1000.);
pythia->SetProcess(kPyMb);
pythia->SetEnergyCMS(energy);
// Tune
// 320 Perugia 0
pythia->SetTune(320);
pythia->UseNewMultipleInteractionsScenario();
//
return pythia;
}
AliGenerator* MbPythiaTunePerugia0Jpsi()
{
comment = comment.Append("pp Cocktail: Pythia (Perugia0) + 1 Jpsi forced to dielectrons");
//Generating a cocktail
AliGenCocktail *gener = new AliGenCocktail();
gener->UsePerEventRates();
//
// Pythia
AliGenPythia* pythia = new AliGenPythia(-1);
pythia->SetMomentumRange(0, 999999.);
pythia->SetThetaRange(0., 180.);
pythia->SetYRange(-12.,12.);
pythia->SetPtRange(0,1000.);
pythia->SetProcess(kPyMb);
pythia->SetEnergyCMS(energy);
// Tune
// 320 Perugia 0
pythia->SetTune(320);
pythia->UseNewMultipleInteractionsScenario();
//
// JPsi is the second ingredient of the cocktail
//AliGenParam* jpsi = new AliGenParam(1, AliGenMUONlib::kJpsi, "CDF scaled", "Jpsi"); // 14 TeV
AliGenParam* jpsi=0x0;
if(JpsiHarderPt) jpsi = new AliGenParam(1, AliGenMUONlib::kJpsi, "CDF pp 8.8", "Jpsi"); // 8.8 TeV
else jpsi = new AliGenParam(1, AliGenMUONlib::kJpsi, "CDF pp 7", "Jpsi"); // 7 TeV
//AliGenParam* jpsi = new AliGenParam(1, AliGenMUONlib::kJpsi, "CDF pp 7", "Jpsi"); // 7 TeV
//AliGenParam* jpsi = new AliGenParam(1, AliGenMUONlib::kJpsi, "CDF pp 8.8", "Jpsi"); // 8.8 TeV
jpsi->SetPtRange(0.,999.);
jpsi->SetYRange(-1.0, 1.0);
jpsi->SetPhiRange(0.,360.);
//jpsi->SetChildYRange(-0.9,0.9); // decaying electrons should be in TRD acceptance
jpsi->SetForceDecay(kDiElectron);
jpsi->Init();
gener->AddGenerator(pythia, "Pythia", 1.); // 1 is weighting factor : Rate per event
gener->AddGenerator(jpsi,"JPsi",1);
// gener->Init();
return gener;
}
AliGenerator* MbPythiaTunePerugia0bchadr()
{
comment = comment.Append(" pp: Pythia (Perugia0) bchadr (1 bbbar per event, 1 c-hadron in |y|<1, chadrons decay to hadrons");
//
// Pythia
AliGenPythia* pythia = new AliGenPythia(-1);
pythia->SetMomentumRange(0, 999999.);
pythia->SetThetaRange(0., 180.);
pythia->SetYRange(-1.,1.);
pythia->SetPtRange(0,1000.);
pythia->SetProcess(kPyBeautyppMNRwmi);
pythia->SetEnergyCMS(energy);
// Tune
// 320 Perugia 0
pythia->SetTune(320);
pythia->UseNewMultipleInteractionsScenario();
//
// decays
pythia->SetForceDecay(kHadronicD);
return pythia;
}