Use CalcMass()

[u/mrichter/AliRoot.git] / EVGEN / AliGenHijing.h

diff --git a/EVGEN/AliGenHijing.h b/EVGEN/AliGenHijing.h
index edf2834fada006ebf6f2b90c3bda5127b3b712ab..8eeb3415108d7cbf81c0ffb4c5bb5f57711dbd51 100644 (file)
--- a/EVGEN/AliGenHijing.h
+++ b/EVGEN/AliGenHijing.h
@@ -34,20 +34,12 @@ class AliGenHijing : public AliGenMC
     virtual void    SetEnergyCMS(Float_t energy=5500) {fEnergyCMS=energy;}
     virtual void    SetReferenceFrame(TString frame="CMS")
        {fFrame=frame;}
-    virtual void    SetProjectile(TString proj="A",
-                                 Int_t a=208, Int_t z=82)
-       {fProjectile = proj;
-       fAProjectile = a;
-       fZProjectile = z;}    
-    virtual void    SetTarget(TString tar="A",
-                             Int_t a=208, Int_t z=82)
-       {fTarget = tar;
-       fATarget = a;
-       fZTarget =z;}    
-    virtual void    SetImpactParameterRange(Float_t bmin = 0, Float_t bmax =15.)
-       { fMinImpactParam=bmin;
-       fMaxImpactParam=bmax;
-       }
+    virtual void    SetProjectile(TString proj="A", Int_t a=208, Int_t z=82)
+       {fProjectile = proj; fAProjectile = a; fZProjectile = z;}    
+    virtual void    SetTarget(TString tar="A", Int_t a=208, Int_t z=82)
+       {fTarget = tar; fATarget = a; fZTarget = z;}    
+    virtual void    SetImpactParameterRange(Float_t bmin = 0, Float_t bmax = 15.)
+       {fMinImpactParam=bmin; fMaxImpactParam=bmax;}
     virtual void    KeepFullEvent();
     virtual void    SetJetQuenching(Int_t flag=1)     {fQuench     = flag;}
     virtual void    SetShadowing(Int_t flag=1)        {fShadowing  = flag;}
@@ -58,20 +50,44 @@ class AliGenHijing : public AliGenMC
     virtual void    SetSelectAll(Int_t flag=0)        {fSelectAll  = flag;}
     virtual void    SetRadiation(Int_t flag=3)        {fRadiation  = flag;}    
     virtual void    SetSpectators(Int_t spects=1)     {fSpectators = spects;}
-    virtual void    SetPtMinJet(Float_t ptmin)        {fPtMinJet   = ptmin;}
+    virtual void    SetPtJet(Float_t ptmin)           {fPtMinJet   = ptmin;}
+    virtual void    SetSimpleJets(Int_t flag=0)       {fSimpleJet  = flag;}
+    virtual void    SetNoGammas(Int_t flag=0)         {fNoGammas   = flag;}
+           
     virtual void    SetJetEtaRange(Float_t etamin = -20., Float_t etamax = 20.)
        {fEtaMinJet = etamin; fEtaMaxJet = etamax;}
-    virtual void    SetBoostLHC(Int_t flag=0)         {fLHC        = flag;}
-    
-           
-    AliGenHijing &  operator=(const AliGenHijing & rhs);
+    virtual void    SetJetPhiRange(Float_t phimin = -180., Float_t phimax = 180.)
+       {fPhiMinJet = TMath::Pi()*phimin/180.; fPhiMaxJet = TMath::Pi()*phimax/180.;}
+    virtual void    SetBoostLHC(Int_t flag = 0)         {fLHC        = flag;}
+// Getters
+    virtual Float_t GetEnergyCMS() {return fEnergyCMS;}
+    virtual TString GetReferenceFrame() {return fFrame;}
+    virtual void    GetProjectile(TString& tar, Int_t& a, Int_t& z)
+       {tar = fProjectile; a = fAProjectile; z = fZProjectile;}    
+    virtual void    GetTarget(TString& tar, Int_t& a, Int_t& z)
+       {tar = fTarget; a = fATarget; z = fZTarget;}    
+    virtual void    GetImpactParameterRange(Float_t& bmin, Float_t& bmax)
+       {bmin = fMinImpactParam; bmax = fMaxImpactParam;}
+    virtual Int_t   GetJetQuenching()                    {return fQuench;}
+    virtual Int_t   GetShadowing()                       {return fShadowing;}
+    virtual Int_t   GetTrigger(Int_t flag=kNoTrigger)    {return fTrigger;}
+    virtual Int_t   GetFlavor(Int_t flag=0)              {return fFlavor;}
+    virtual Int_t   GetRadiation(Int_t flag=3)           {return fRadiation;}    
+    virtual Int_t   GetSpectators(Int_t spects=1)        {return fSpectators;}
+    virtual Float_t GetPtJet(Float_t ptmin)              {return fPtMinJet;}
+    virtual void    GetJetEtaRange(Float_t& etamin, Float_t& etamax)
+       {etamin = fEtaMinJet; etamax = fEtaMaxJet;}
+    virtual void    GetJetPhiRange(Float_t& phimin, Float_t& phimax)
+       {phimin = fPhiMinJet*180./TMath::Pi(); phimax = fPhiMaxJet*180./TMath::Pi();}
+
 // Physics Routines        
     virtual void EvaluateCrossSections();
     virtual void Boost(TClonesArray* particles);
     virtual TGraph* CrossSection()     {return fDsigmaDb;}
     virtual TGraph* BinaryCollisions() {return fDnDb;}
     virtual Bool_t  CheckTrigger();
-    
+//
+    AliGenHijing &  operator=(const AliGenHijing & rhs);
  protected:
     Bool_t SelectFlavor(Int_t pid);
     void   MakeHeader();
@@ -109,7 +125,11 @@ class AliGenHijing : public AliGenMC
     Float_t     fPtMinJet;       // Minimum Pt of triggered Jet
     Float_t     fEtaMinJet;      // Minimum eta of triggered Jet
     Float_t     fEtaMaxJet;      // Maximum eta of triggered Jet
+    Float_t     fPhiMinJet;      // At least one of triggered Jets must be in this
+    Float_t     fPhiMaxJet;      // phi range
     Int_t       fRadiation;      // Flag to switch on/off initial and final state radiation
+    Int_t       fSimpleJet;      // Flag to produce simple tiggered jet topology
+    Int_t       fNoGammas;       // Don't write gammas if flag "on"
     
 // ZDC proposal (by Chiara) to store num. of SPECTATORS protons and neutrons
     Int_t      fSpecn;          // Num. of spectator neutrons
@@ -124,7 +144,7 @@ class AliGenHijing : public AliGenMC
     // check if stable
     Bool_t Stable(TParticle*  particle);
     
-    ClassDef(AliGenHijing,1) // AliGenerator interface to Hijing
+    ClassDef(AliGenHijing,3) // AliGenerator interface to Hijing
 };
 #endif