AliGenerator* MbCocktail()
{
comment = comment.Append(" pp at 14 TeV: Pythia low-pt, no heavy quarks + J/Psi from parameterisation");
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(14000.);
pythia->SwitchHFOff();
// J/Psi parameterisation
AliGenParam* jpsi = new AliGenParam(1, AliGenMUONlib::kJpsi, "CDF scaled", "Jpsi");
jpsi->SetPtRange(0.,100.);
jpsi->SetYRange(-8., 8.);
jpsi->SetPhiRange(0., 360.);
jpsi->SetForceDecay(kAll);
gener->AddGenerator(pythia, "Pythia", 1.);
gener->AddGenerator(jpsi, "J/Psi", 8.e-4);
return gener;
}
AliGenerator* JPsiPbPb()
{
AliGenParam *jpsiPbPb = new AliGenParam(1, AliGenMUONlib::kJpsi,"CDF PbPb 3.94");
jpsiPbPb->SetMomentumRange(0,999);
jpsiPbPb->SetPtRange(0,50.);
jpsiPbPb->SetYRange(-4.5,-2.);
jpsiPbPb->SetPhiRange(0., 360.);
jpsiPbPb->SetCutOnChild(1);
jpsiPbPb->SetChildPhiRange(0.,360.);
jpsiPbPb->SetChildThetaRange(0.,180.);
jpsiPbPb->SetForceDecay(kDiMuon);
jpsiPbPb->SetTrackingFlag(1);
return jpsiPbPb;
}
AliGenerator* FlatJPsi()
{
AliGenParam *flatjpsi = new AliGenParam(1, AliGenMUONlib::kJpsi,"Flat");
flatjpsi->SetMomentumRange(0,999);
flatjpsi->SetPtRange(0,20.);
flatjpsi->SetYRange(-4.2,-2.3);
flatjpsi->SetPhiRange(0., 360.);
flatjpsi->SetCutOnChild(1);
flatjpsi->SetChildPhiRange(0.,360.);
flatjpsi->SetChildThetaRange(0.,180.);
flatjpsi->SetForceDecay(kDiMuon);
flatjpsi->SetTrackingFlag(1);
return flatjpsi;
}
AliGenerator* GenParamJpsi() {
AliGenParam *gener = new AliGenParam(5, AliGenMUONlib::kJpsi);
gener->SetMomentumRange(0,999);
gener->SetPtRange(0,100.);
gener->SetYRange(-4.0, -2.5);
gener->SetPhiRange(0., 360.);
gener->SetChildThetaRange(171.0,177.0);
gener->SetChildMomentumRange(4.0, 999.);
gener->SetForceDecay(kDiMuon);
gener->SetTrackingFlag(1);
gener->SetCutOnChild(1);
return gener;
}
AliGenerator* JPsi7TeV()
/*{
AliGenParam *jpsi7TeV = new AliGenParam(1, AliGenMUONlib::kJpsi,"CDF pp 7");
jpsi7TeV->SetMomentumRange(0,999);
jpsi7TeV->SetPtRange(0,20.);
jpsi7TeV->SetYRange(-4.2,-2.3);
jpsi7TeV->SetPhiRange(0., 360.);
jpsi7TeV->SetCutOnChild(1);
jpsi7TeV->SetChildPhiRange(0.,360.);
jpsi7TeV->SetChildThetaRange(0.,180.);
jpsi7TeV->SetForceDecay(kDiMuon);
jpsi7TeV->SetTrackingFlag(1);
return jpsi7TeV;
}*/
{
comment = comment.Append(" pp at 7 TeV: Pythia low-pt, no heavy quarks + J/Psi from parameterisation");
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);
pythia->SwitchHFOff();
// J/Psi parameterisation
// AliGenParam* jpsi = new AliGenParam(1, AliGenMUONlib::kJpsi, "CDF scaled", "Jpsi");
AliGenParam* jpsi = new AliGenParam(1, AliGenMUONlib::kJpsi,"CDF pp 7");
jpsi->SetPtRange(0.,100.);
// jpsi->SetYRange(-8., 8.);
jpsi->SetYRange(-1.5, 1.5); //this to speed up !!!
jpsi->SetPhiRange(0., 360.);
jpsi->SetForceDecay(kDiElectron); // this has to be changed to force J/psi->e+e-
//gener->AddGenerator(jpsi, "J/Psi", 8.e-4); // expected yield
gener->AddGenerator(jpsi, "J/Psi", 1.); // 1 J/psi per event
gener->AddGenerator(pythia, "Pythia", 1.);
return gener;
}