Improved common vertex handling.

[u/mrichter/AliRoot.git] / STEER / AliGenerator.h

diff --git a/STEER/AliGenerator.h b/STEER/AliGenerator.h
index 4b866b31ed660a6b033db1f34755c8dd10cd6a10..b2362690cd496797937fe90eef58cee6cd88fb5f 100644 (file)
--- a/STEER/AliGenerator.h
+++ b/STEER/AliGenerator.h
@@ -12,13 +12,22 @@
 //                                                       //
 ///////////////////////////////////////////////////////////
 
-class TGenerator;
+#include <TArrayF.h>
+#include <TLorentzVector.h>
+#include <TMCProcess.h>
 
-#include "TLorentzVector.h"
-#include "TArrayF.h"
 #include "AliRndm.h"
 
-typedef enum { kNoSmear, kPerEvent, kPerTrack } VertexSmear_t;
+class AliCollisionGeometry;
+class AliStack;
+
+class TGenerator;
+
+
+
+typedef enum { kNoSmear, kPerEvent, kPerTrack }  VertexSmear_t;
+typedef enum { kExternal, kInternal, kContainer} VertexSource_t;
+
 
 class AliGenerator : public TNamed, public AliRndm
 {
@@ -34,7 +43,7 @@ class AliGenerator : public TNamed, public AliRndm
     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=100);
+    virtual void SetPtRange(Float_t ptmin=0, Float_t ptmax=1.e10);
     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,
@@ -48,8 +57,15 @@ class AliGenerator : public TNamed, public AliRndm
     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 source = kInternal) {fVertexSource = source;}
+    virtual void SetVertex(Float_t vx, Float_t vy, Float_t vz)
+       {fVertex[0] = vx; fVertex[1] = vy; fVertex[2] = vz;}
     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;}
 
@@ -61,9 +77,29 @@ class AliGenerator : public TNamed, public AliRndm
        {ox=fOrigin.At(0);oy=fOrigin.At(1);oz=fOrigin.At(2);}
     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
+  // Stack 
+    void SetStack (AliStack *stack) {fStack = stack;}
+    AliStack* GetStack(){return fStack;}
+  // Collision Geometry
+    virtual Bool_t ProvidesCollisionGeometry() const {return kFALSE;}
+    virtual Bool_t NeedsCollisionGeometry()    const {return kFALSE;}    
+    virtual AliCollisionGeometry* CollisionGeometry() {return fCollisionGeometry;}
+    virtual void SetCollisionGeometry(AliCollisionGeometry* geom) {fCollisionGeometry = geom;}
+ 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, TMCProcess mech, Int_t &ntr,
+                               Float_t weight = 1, Int_t is = 0);
+    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,
+                      TMCProcess mech, Int_t &ntr, Float_t weight = 1, Int_t is = 0);
+    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
@@ -79,18 +115,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
 
-    ClassDef(AliGenerator,1) // Base class for event generators
+    AliStack*   fStack;      //! Local pointer to stack
+    AliCollisionGeometry* fCollisionGeometry; //!Collision geometry
+    /*************************************************************************/
+    enum {kThetaRange    = BIT(14),
+         kVertexRange   = BIT(15),
+         kPhiRange      = BIT(16),
+         kPtRange       = BIT(17),
+         kYRange        = BIT(18),
+         kMomentumRange = BIT(19)     
+    };
+    ClassDef(AliGenerator,3) // Base class for event generators
 };
 
 #endif
+
+
+
+
+
+
+
+
+
+
+
+
+
+