]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - STEER/AliGenerator.h
Debug staff moved to AliRICH
[u/mrichter/AliRoot.git] / STEER / AliGenerator.h
index d773e9dbe867f364913727f5b92cfa87d0e82685..e9c41dcd64e3fc75a6c44c36843f8e9ce5a689a0 100644 (file)
 
 #include "TLorentzVector.h"
 #include "TArrayF.h"
-#include "TGenerator.h"
+#include "AliRndm.h"
+#include "AliMCProcess.h"
+class TGenerator;
+class AliStack;
+
 
 typedef enum { kNoSmear, kPerEvent, kPerTrack } VertexSmear_t;
+typedef enum { kExternal, kInternal}            VertexSource_t;
+
 
-class AliGenerator : public TNamed
+class AliGenerator : public TNamed, public AliRndm
 {
 
  public:
@@ -32,7 +38,7 @@ class AliGenerator : public TNamed
     virtual void SetOrigin(const TLorentzVector &o);
     virtual void SetSigma(Float_t sx, Float_t sy, Float_t sz);
     virtual void SetMomentumRange(Float_t pmin=0, Float_t pmax=1.e10);
-    virtual void SetPtRange(Float_t ptmin=0, Float_t ptmax=1.e10);
+    virtual void SetPtRange(Float_t ptmin=0, Float_t ptmax=100);
     virtual void SetPhiRange(Float_t phimin=-180., Float_t phimax=180);
     virtual void SetYRange(Float_t ymin=-100, Float_t ymax=100);
     virtual void SetVRange(Float_t vxmin, Float_t vxmax,
@@ -46,8 +52,14 @@ class AliGenerator : public TNamed
     virtual void SetChildWeight(Float_t wgt)  {fChildWeight=wgt;}    
     virtual void SetAnalog(Int_t flag=1) {fAnalog=flag;}       
     virtual void SetVertexSmear(VertexSmear_t smear) {fVertexSmear = smear;}
+    virtual void SetCutVertexZ(Float_t cut=999999.) {fCutVertexZ = cut;}
+    virtual void SetVertexSource(VertexSource_t smear) {fVertexSource = kInternal;}    
     virtual void SetTrackingFlag(Int_t flag=1) {fTrackIt=flag;}
-           
+    void Vertex();
+    void VertexExternal();
+    virtual void VertexInternal();
+    virtual void FinishRun(){;}
+    
     virtual void SetMC(TGenerator *theMC) 
        {if (!fgMCEvGen) fgMCEvGen =theMC;}
 
@@ -60,8 +72,24 @@ class AliGenerator : public TNamed
     virtual void GetOrigin(TLorentzVector &o) const
        {o[0]=fOrigin.At(0);o[1]=fOrigin.At(1);o[2]=fOrigin.At(2);o[3]=0;}
 
-protected:
-    static TGenerator* fgMCEvGen; // Pointer to the generator
+    void SetStack (AliStack *stack) {fStack = stack;}
+    AliStack* GetStack(){return fStack;}
+// Comminication with stack
+ protected:
+    virtual  void  SetTrack(Int_t done, Int_t parent, Int_t pdg,
+                               Float_t *pmom, Float_t *vpos, Float_t *polar,
+                               Float_t tof, AliMCProcess mech, Int_t &ntr,
+                               Float_t weight=1);
+    virtual  void  SetTrack(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,
+                      AliMCProcess mech, Int_t &ntr, Float_t weight=1);
+    virtual void   KeepTrack(Int_t itrack); 
+    virtual void   SetHighWaterMark(Int_t nt);
+    
+ protected:
+    static  TGenerator* fgMCEvGen; // Pointer to the generator
     Float_t     fThetaMin;     //Minimum theta of generation in radians
     Float_t     fThetaMax;     //Maximum theta of generation in radians
     Float_t     fPhiMin;       //Minimum phi of generation in radians
@@ -77,18 +105,41 @@ protected:
     Int_t       fNpart;        //Maximum number of particles per event
     Float_t     fParentWeight; //Parent Weight
     Float_t     fChildWeight;  //ChildWeight
-    Int_t       fTrackit;      // Track the generated final state particle if 1
-    Int_t       fAnalog;       //Flaf for anolog or pt-weighted generation
+    Int_t       fAnalog;       //Flag for anolog or pt-weighted generation
    //
-    VertexSmear_t     fVertexSmear; //Vertex Smearing mode
-    Int_t       fTrackIt;    // if 1 Track final state particles 
-    TArrayF     fOrigin;     //Origin of event
-    TArrayF     fOsigma;     //Sigma of the Origin of event
-
-    enum {kThetaRange=1, kVertexRange=2, kPhiRange=4, kPtRange=8,
-         kYRange=32, kMomentumRange=16};
+    VertexSmear_t     fVertexSmear;  //Vertex Smearing mode
+    VertexSource_t    fVertexSource; //Vertex source (internal/external)
+    Float_t     fCutVertexZ;    // Vertex cut in units of sigma_z
+    Int_t       fTrackIt;    // if 1, Track final state particles 
+    TArrayF     fOrigin;     // Origin of event
+    TArrayF     fOsigma;     // Sigma of the Origin of event
+    TArrayF     fVertex;     //! Vertex of current event
+    TArrayF     fEventVertex;   //!The current event vertex
 
-    ClassDef(AliGenerator,1) // Base class for event generators
+    AliStack*   fStack;      //! Local pointer to stack
+    /*************************************************************************/
+    enum {kThetaRange    = BIT(14),
+         kVertexRange   = BIT(15),
+         kPhiRange      = BIT(16),
+         kPtRange       = BIT(17),
+         kYRange        = BIT(18),
+         kMomentumRange = BIT(19)     
+    };
+    ClassDef(AliGenerator,2) // Base class for event generators
 };
 
 #endif
+
+
+
+
+
+
+
+
+
+
+
+
+
+