- //
- // The cocktail
- AliGenCocktail *gener = new AliGenCocktail();
-
- //
- // Charm production by Pythia
- AliGenPythia * genpyc = new AliGenPythia(230);
- genpyc->SetProcess(kPyCharmPbPbMNR);
- genpyc->SetStrucFunc(kCTEQ4L);
- genpyc->SetPtHard(2.1,-1.0);
- genpyc->SetEnergyCMS(5500.);
- genpyc->SetNuclei(208,208);
- genpyc->SetYRange(-999,999);
- genpyc->SetForceDecay(kAll);
- genpyc->SetFeedDownHigherFamily(kFALSE);
- genpyc->SetCountMode(AliGenPythia::kCountParents);
- //
- // Beauty production by Pythia
- AliGenPythia * genpyb = new AliGenPythia(9);
- genpyb->SetProcess(kPyBeautyPbPbMNR);
- genpyb->SetStrucFunc(kCTEQ4L);
- genpyb->SetPtHard(2.75,-1.0);
- genpyb->SetEnergyCMS(5500.);
- genpyb->SetNuclei(208,208);
- genpyb->SetYRange(-999,999);
- genpyb->SetForceDecay(kAll);
- genpyb->SetFeedDownHigherFamily(kFALSE);
- genpyb->SetCountMode(AliGenPythia::kCountParents);
- //
- // Hyperons
- //
- AliGenSTRANGElib *lib = new AliGenSTRANGElib();
- Int_t particle;
- // Xi
- particle = kXiMinus;
- AliGenParam *genXi = new AliGenParam(16,particle,lib->GetPt(particle),lib->GetY(particle),lib->GetIp(particle));
- genXi->SetPtRange(0., 12.);
- genXi->SetYRange(-1.1, 1.1);
- genXi->SetForceDecay(kNoDecay);
-
- //
- // Omega
- particle = kOmegaMinus;
- AliGenParam *genOmega = new AliGenParam(10,particle,lib->GetPt(particle),lib->GetY(particle),lib->GetIp(particle));
- genOmega->SetPtRange(0, 12.);
- genOmega->SetYRange(-1.1, 1.1);
- genOmega->SetForceDecay(kNoDecay);
-
- //
- // Central Hijing
- AliGenHijing *genHi = HijingStandard();
- genHi->SwitchOffHeavyQuarks(kTRUE);
- genHi->SetImpactParameterRange(0.,5.);
- //
- // Add everything to the cocktail and shake ...
- gener->AddGenerator(genHi, "Hijing cent1", 1);
- gener->AddGenerator(genpyc, "Extra charm", 1);
- gener->AddGenerator(genpyb, "Extra beauty", 1);
- gener->AddGenerator(genXi, "Xi" , 1);
- gener->AddGenerator(genOmega, "Omega", 1);
- gGener = gener;
- }
- break;
- default: break;
+ //
+ // The cocktail
+ AliGenCocktail *gener = new AliGenCocktail();
+
+ //
+ // Charm production by Pythia
+ AliGenPythia * genpyc = new AliGenPythia(230);
+ genpyc->SetProcess(kPyCharmPbPbMNR);
+ genpyc->SetStrucFunc(kCTEQ4L);
+ genpyc->SetPtHard(2.1,-1.0);
+ genpyc->SetEnergyCMS(5500.);
+ genpyc->SetNuclei(208,208);
+ genpyc->SetYRange(-999,999);
+ genpyc->SetForceDecay(kAll);
+ genpyc->SetFeedDownHigherFamily(kFALSE);
+ genpyc->SetCountMode(AliGenPythia::kCountParents);
+ //
+ // Beauty production by Pythia
+ AliGenPythia * genpyb = new AliGenPythia(9);
+ genpyb->SetProcess(kPyBeautyPbPbMNR);
+ genpyb->SetStrucFunc(kCTEQ4L);
+ genpyb->SetPtHard(2.75,-1.0);
+ genpyb->SetEnergyCMS(5500.);
+ genpyb->SetNuclei(208,208);
+ genpyb->SetYRange(-999,999);
+ genpyb->SetForceDecay(kAll);
+ genpyb->SetFeedDownHigherFamily(kFALSE);
+ genpyb->SetCountMode(AliGenPythia::kCountParents);
+ //
+ // Hyperons
+ //
+ AliGenSTRANGElib *lib = new AliGenSTRANGElib();
+ Int_t particle;
+ // Xi
+ particle = kXiMinus;
+ AliGenParam *genXi = new AliGenParam(16,particle,lib->GetPt(particle),lib->GetY(particle),lib->GetIp(particle));
+ genXi->SetPtRange(0., 12.);
+ genXi->SetYRange(-1.1, 1.1);
+ genXi->SetForceDecay(kNoDecay);
+
+ //
+ // Omega
+ particle = kOmegaMinus;
+ AliGenParam *genOmega = new AliGenParam(10,particle,lib->GetPt(particle),lib->GetY(particle),lib->GetIp(particle));
+ genOmega->SetPtRange(0, 12.);
+ genOmega->SetYRange(-1.1, 1.1);
+ genOmega->SetForceDecay(kNoDecay);
+
+ //
+ // Central Hijing
+ AliGenHijing *genHi = HijingStandard();
+ genHi->SwitchOffHeavyQuarks(kTRUE);
+ genHi->SetImpactParameterRange(0.,5.);
+ //
+ // Add everything to the cocktail and shake ...
+ gener->AddGenerator(genHi, "Hijing cent1", 1);
+ gener->AddGenerator(genpyc, "Extra charm", 1);
+ gener->AddGenerator(genpyb, "Extra beauty", 1);
+ gener->AddGenerator(genXi, "Xi" , 1);
+ gener->AddGenerator(genOmega, "Omega", 1);
+ gener->SetProjectile("A", 208, 82) ;
+ gener->SetTarget("A", 208, 82) ;
+ gGener = gener;