// Is a generalisation of AliGenParam class for correlated pairs of hadrons.
// In this version quark pairs and fragmentation functions are obtained from
// ~2.10^6 Pythia6.214 events generated with kCharmppMNRwmi & kBeautyppMNRwmi,
-// CTEQ5L PDF and Pt_hard = 2.76 GeV/c for p-p collisions at 7, 10 and 14 TeV,
+// CTEQ5L PDF and Pt_hard = 2.76 GeV/c for p-p collisions at 2.76, 7, 8, 10 and 14 TeV,
// and with kCharmppMNR (Pt_hard = 2.10 GeV/c) & kBeautyppMNR (Pt_hard = 2.75 GeV/c),
-// CTEQ4L PDF for Pb-Pb at 3.94 TeV, for p-Pb & Pb-p at 8.8 TeV.
+// CTEQ4L PDF for Pb-Pb at 2.76 and 3.94 TeV, for p-Pb & Pb-p at 5 and 8.8 TeV.
// Decays are performed by Pythia.
// Author: S. Grigoryan, LPC Clermont-Fd & YerPhI, Smbat.Grigoryan@cern.ch
// July 07: added quarks in the stack (B. Vulpescu)
// of HF-hadrons, when the cuts on their children are satisfied (L. Manceau)
// Oct 11: added Pb-Pb at 2.76 TeV (S. Grigoryan)
// June 12: added p-Pb & Pb-p at 5 TeV (S. Grigoryan)
+// April 13: added p-p at 2.76 and 8 TeV (S. Grigoryan)
//
//-------------------------------------------------------------------------
// How it works (for the given flavor and p-p energy):
//
// add the following typical lines in Config.C
/*
- if (!strcmp(option,"corr")) {
// An example for correlated charm or beauty hadron pair production at 14 TeV
// AliGenCorrHF *gener = new AliGenCorrHF(1, 4, 14); // for charm, 1 pair per event
gener->SetOrigin(0,0,0); //vertex position
gener->SetSigma(0,0,0); //Sigma in (X,Y,Z) (cm) on IP position
gener->SetForceDecay(kSemiMuonic);
- gener->SetTrackingFlag(0);
+ gener->SetSelectAll(kTRUE); //Force the transport of all particles.
+ //Necessary while using a different
+ //option than kAll for SetForceDecay
+ gener->SetTrackingFlag(1); //1: Decay during transport,
+ //0: No Decay during transport
gener->Init();
-}
*/
-// and in aliroot do e.g. gAlice->Run(10,"Config.C") to produce 10 events.
// One can include AliGenCorrHF in an AliGenCocktail generator.
//--------------------------------------------------------------------------
if (fQuark == 5) {
if (fEnergy == 7)
fFileName = "$ALICE_ROOT/EVGEN/dataCorrHF/BeautyPP7PythiaMNRwmi.root";
+ else if (fEnergy == 8)
+ fFileName = "$ALICE_ROOT/EVGEN/dataCorrHF/BeautyPP8PythiaMNRwmi.root";
else if (fEnergy == 10)
fFileName = "$ALICE_ROOT/EVGEN/dataCorrHF/BeautyPP10PythiaMNRwmi.root";
else if (fEnergy == 14)
fFileName = "$ALICE_ROOT/EVGEN/dataCorrHF/BeautyPP14PythiaMNRwmi.root";
+ else if (fEnergy == 2)
+ fFileName = "$ALICE_ROOT/EVGEN/dataCorrHF/BeautyPP276PythiaMNRwmi.root";
else if (fEnergy == 3)
fFileName = "$ALICE_ROOT/EVGEN/dataCorrHF/BeautyPbPb276PythiaMNR.root";
else if (fEnergy == 4)
fQuark = 4;
if (fEnergy == 7)
fFileName = "$ALICE_ROOT/EVGEN/dataCorrHF/CharmPP7PythiaMNRwmi.root";
+ else if (fEnergy == 8)
+ fFileName = "$ALICE_ROOT/EVGEN/dataCorrHF/CharmPP8PythiaMNRwmi.root";
else if (fEnergy == 10)
fFileName = "$ALICE_ROOT/EVGEN/dataCorrHF/CharmPP10PythiaMNRwmi.root";
else if (fEnergy == 14)
fFileName = "$ALICE_ROOT/EVGEN/dataCorrHF/CharmPP14PythiaMNRwmi.root";
+ else if (fEnergy == 2)
+ fFileName = "$ALICE_ROOT/EVGEN/dataCorrHF/CharmPP276PythiaMNRwmi.root";
else if (fEnergy == 3)
fFileName = "$ALICE_ROOT/EVGEN/dataCorrHF/CharmPbPb276PythiaMNR.root";
else if (fEnergy == 4)
// particle decay related initialization
- if (gMC) fDecayer = gMC->GetDecayer();
+ if (TVirtualMC::GetMC()) fDecayer = TVirtualMC::GetMC()->GetDecayer();
fDecayer->SetForceDecay(fForceDecay);
fDecayer->Init();
TParticle* iparticle = (TParticle *) particles->At(i);
Int_t kf = iparticle->GetPdgCode();
Int_t jpa = iparticle->GetFirstMother()-1;
+ Int_t ksc = iparticle->GetStatusCode();
// RS: note, the conversion mm->cm is done now in the decayer. The time is ignored here!
och[0] = origin0[0]+iparticle->Vx();
och[1] = origin0[1]+iparticle->Vy();
PushTrack(fTrackIt*trackIt[i], iparent, kf,
pc, och, polar,
- 0, kPDecay, nt, wgtch);
+ 0, kPDecay, nt, wgtch,ksc);
pParent[i] = nt;
KeepTrack(nt);
fNprimaries++;
return;
}
-