class AliGenLightNuclei: public AliGenCocktail
Afterburner to generate light nuclei for event generators
such as PYTHIA and PHOJET
Light nuclei are generated whenever a cluster of nucleons is found
within a sphere of radius p0 (coalescence momentum), i.e. have the
same momentum.
By default it starts with He4 nuclei which are the most stable,
then He3 nuclei and tritons and finally deuterons. It can also generate
a single nucleus species by disabling the others.
Sample code for PYTHIA:
AliGenLightNuclei* gener = new AliGenLightNuclei();
AliGenPythia* pythia = new AliGenPythia(-1);
pythia->SetCollisionSystem("p+", "p+");
pythia->SetEnergyCMS(7000);
gener->UsePerEventRates();
gener->AddGenerator(pythia, "PYTHIA", 1);
gener->SetCoalescenceMomentum(0.200); // default (GeV/c)
Function Members (Methods)
public:
protected:
| virtual void | TObject::DoError(int level, const char* location, const char* fmt, va_list va) const |
| virtual void | AliGenerator::KeepTrack(Int_t itrack) |
| void | TObject::MakeZombie() |
| virtual void | AliGenerator::PushTrack(Int_t done, Int_t parent, Int_t pdg, Float_t* pmom, Float_t* vpos, Float_t* polar, Float_t tof, TMCProcess mech, Int_t& ntr, Float_t weight = 1, Int_t is = 0) |
| virtual void | AliGenerator::PushTrack(Int_t done, Int_t parent, Int_t pdg, Double_t px, Double_t py, Double_t pz, Double_t e, Double_t vx, Double_t vy, Double_t vz, Double_t tof, Double_t polx, Double_t poly, Double_t polz, TMCProcess mech, Int_t& ntr, Float_t weight = 1, Int_t is = 0) |
| virtual void | AliGenerator::SetHighWaterMark(Int_t nt) |
private:
| AliGenLightNuclei(const AliGenLightNuclei& other) | |
| Bool_t | Coalescence(const TParticle* n1, const TParticle* n2) const |
| TParticle* | FindPartner(const TParticle* n0, const TList* nucleons, const TParticle* nx = 0) const |
| Int_t | GenerateDeuterons(const TList* protons, const TList* neutrons) |
| Int_t | GenerateHe3Nuclei(const TList* protons, const TList* neutrons) |
| Int_t | GenerateHe4Nuclei(const TList* protons, const TList* neutrons) |
| Int_t | GenerateTritons(const TList* protons, const TList* neutrons) |
| Double_t | GetPcm(Double_t p1x, Double_t p1y, Double_t p1z, Double_t m1, Double_t p2x, Double_t p2y, Double_t p2z, Double_t m2) const |
| Double_t | GetS(Double_t p1x, Double_t p1y, Double_t p1z, Double_t m1, Double_t p2x, Double_t p2y, Double_t p2z, Double_t m2) const |
| AliGenLightNuclei& | operator=(const AliGenLightNuclei& other) |
| void | PushDeuteron(TParticle* parent1, TParticle* parent2) |
| void | PushHe3Nucleus(TParticle* parent1, TParticle* parent2, TParticle* parent3) |
| void | PushHe4Nucleus(TParticle* parent1, TParticle* parent2, TParticle* parent3, TParticle* parent4) |
| void | PushNucleus(Int_t pdg, Double_t mass, TParticle* parent1, TParticle* parent2, TParticle* parent3 = 0, TParticle* parent4 = 0) |
| void | PushTriton(TParticle* parent1, TParticle* parent2, TParticle* parent3) |
Data Members
public:
| enum { | kDeuteron | |
| kAntiDeuteron | ||
| kTriton | ||
| kAntiTriton | ||
| kHe3Nucleus | ||
| kAntiHe3Nucleus | ||
| kAlpha | ||
| kAntiAlpha | ||
| kCluster | ||
| }; | ||
| enum AliGenerator::[unnamed] { | kThetaRange | |
| kVertexRange | ||
| kPhiRange | ||
| kPtRange | ||
| kYRange | ||
| kMomentumRange | ||
| kEtaRange | ||
| }; | ||
| enum TObject::EStatusBits { | kCanDelete | |
| kMustCleanup | ||
| kObjInCanvas | ||
| kIsReferenced | ||
| kHasUUID | ||
| kCannotPick | ||
| kNoContextMenu | ||
| kInvalidObject | ||
| }; | ||
| enum TObject::[unnamed] { | kIsOnHeap | |
| kNotDeleted | ||
| kZombie | ||
| kBitMask | ||
| kSingleKey | ||
| kOverwrite | ||
| kWriteDelete | ||
| }; |
protected:
| Int_t | AliGenerator::fAProjectile | Projectile A |
| Int_t | AliGenerator::fATarget | Target A |
| Int_t | AliGenerator::fAnalog | Flag for anolog or pt-weighted generation |
| Float_t | AliGenerator::fChildWeight | ChildWeight |
| AliCollisionGeometry* | AliGenerator::fCollisionGeometry | !Collision geometry |
| AliGenerator* | AliGenerator::fContainer | ! Local pointer to container |
| Float_t | AliGenerator::fCutVertexZ | Vertex cut in units of sigma_z |
| Float_t | AliGenerator::fEnergyCMS | Centre of mass energy |
| TList* | AliGenCocktail::fEntries | List of Generators |
| Float_t | AliGenerator::fEvPlane | the event plane |
| AliGenCocktailEventHeader* | AliGenCocktail::fHeader | !Header container |
| TGenerator* | AliGenerator::fMCEvGen | !Pointer to the generator |
| Int_t | AliGenCocktail::fNGenerators | Number of generators booked |
| TString | TNamed::fName | object identifier |
| Int_t | AliGenerator::fNpart | Maximum number of particles per event |
| TArrayF | AliGenerator::fOrigin | Origin of event |
| TArrayF | AliGenerator::fOsigma | Sigma of the Origin of event |
| Float_t | AliGenerator::fPMax | Minimum momentum of generation in GeV/c |
| Float_t | AliGenerator::fPMin | Minimum momentum of generation in GeV/c |
| Float_t | AliGenerator::fParentWeight | Parent Weight |
| Float_t | AliGenerator::fPhiMax | Maximum phi of generation in radians |
| Float_t | AliGenerator::fPhiMin | Minimum phi of generation in radians |
| Float_t | AliGenerator::fPileUpTimeWindow | Time window for pile-up events |
| TArrayF | AliGenCocktail::fProb | Probability of an event (if fSRandom == kTRUE) |
| TString | AliGenerator::fProjectile | Projectile |
| Float_t | AliGenerator::fPtMax | Maximum transverse momentum |
| Float_t | AliGenerator::fPtMin | Minimum transverse momentum |
| TRandom* | AliRndm::fRandom | ! Pointer to the random number generator |
| Bool_t | AliGenCocktail::fSRandom | Flag to select random generator from list |
| AliStack* | AliGenerator::fStack | ! Local pointer to stack |
| TString | AliGenerator::fTarget | Target |
| Float_t | AliGenerator::fThetaMax | Maximum theta of generation in radians |
| Float_t | AliGenerator::fThetaMin | Minimum theta of generation in radians |
| Float_t | AliGenerator::fTime | Event time smeared around time0 origin using sigma vertex |
| Float_t | AliGenerator::fTimeOrigin | Time0 origin in a run or event sample |
| TString | TNamed::fTitle | object title |
| Float_t | AliGenCocktail::fTotalRate | Total rate of the generators booked |
| Int_t | AliGenerator::fTrackIt | if 1, Track final state particles |
| Bool_t | AliGenCocktail::fUsePerEventRate | Flag to generate the events according to the rate per event |
| TArrayF | AliGenerator::fVMax | Minimum Decaylength |
| TArrayF | AliGenerator::fVMin | Minimum Decaylength |
| TArrayF | AliGenerator::fVertex | ! Vertex of current event |
| AliVertexGenerator* | AliGenerator::fVertexGenerator | ! Generator for the vertex |
| VertexSmear_t | AliGenerator::fVertexSmear | Vertex Smearing mode |
| VertexSource_t | AliGenerator::fVertexSource | Vertex source (internal/external) |
| Float_t | AliGenerator::fYMax | Maximum rapidity |
| Float_t | AliGenerator::fYMin | Minimum rapidity |
| Int_t | AliGenerator::fZProjectile | Projectile Z |
| Int_t | AliGenerator::fZTarget | Target Z |
| TObjLink* | AliGenCocktail::flnk1 | ! Iterator for first generator |
| TObjLink* | AliGenCocktail::flnk2 | ! Iterator for second generator |
Class Charts
| Inheritance Chart: | ||||||||||||||||||
|
Function documentation
void Generate()
delegate the particle generation to the cocktail and modify the stack adding the light nuclei
Bool_t Coalescence(const TParticle* n1, const TParticle* n2) const
returns true if the nucleons are inside of an sphere of radius p0 (assume the nucleons are in the same place e.g. PYTHIA, PHOJET,...)
TParticle* FindPartner(const TParticle* n0, const TList* nucleons, const TParticle* nx = 0) const
find the first nucleon partner within a sphere of radius p0 centered at n0 and exclude nucleon nx
Int_t GenerateDeuterons(const TList* protons, const TList* neutrons)
a deuteron is generated from a pair of p-n nucleons (the center of the sphere is one of the nucleons)
Int_t GenerateTritons(const TList* protons, const TList* neutrons)
a triton is generated from a cluster of p-n-n nucleons with same momentum (triangular configuration)
Int_t GenerateHe3Nuclei(const TList* protons, const TList* neutrons)
a He3 nucleus is generated from a cluster of p-n-p nucleons with same momentum (triangular configuration)
Int_t GenerateHe4Nuclei(const TList* protons, const TList* neutrons)
a He4 nucleus is generated from a cluster of p-n-p-n nucleons with same momentum (tetrahedron configuration)
void PushTriton(TParticle* parent1, TParticle* parent2, TParticle* parent3)
push a triton to the particle stack
void PushHe3Nucleus(TParticle* parent1, TParticle* parent2, TParticle* parent3)
push a He3 nucleus to the particle stack
void PushHe4Nucleus(TParticle* parent1, TParticle* parent2, TParticle* parent3, TParticle* parent4)
push a He4 nucleus to the particle stack
void PushNucleus(Int_t pdg, Double_t mass, TParticle* parent1, TParticle* parent2, TParticle* parent3 = 0, TParticle* parent4 = 0)
push a nucleus to the stack and tag the parents with the kCluster status code
Double_t GetS(Double_t p1x, Double_t p1y, Double_t p1z, Double_t m1, Double_t p2x, Double_t p2y, Double_t p2z, Double_t m2) const
square of the center of mass energy
Double_t GetPcm(Double_t p1x, Double_t p1y, Double_t p1z, Double_t m1, Double_t p2x, Double_t p2y, Double_t p2z, Double_t m2) const
momentum in the CM frame for 2 particles
AliGenLightNuclei& operator=(const AliGenLightNuclei& other)