Event weight stored in the generator event header,
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 20 Sep 2010 08:55:06 +0000 (08:55 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 20 Sep 2010 08:55:06 +0000 (08:55 +0000)
STEER/AliGenEventHeader.cxx
STEER/AliGenEventHeader.h
TEPEMGEN/AliGenEpEmv1.cxx
TEPEMGEN/AliGenEpEmv1.h

index db9c6c1..e5092b7 100644 (file)
@@ -31,7 +31,8 @@ ClassImp(AliGenEventHeader)
 AliGenEventHeader::AliGenEventHeader():
   fNProduced(-1),
   fVertex(3),
-  fInteractionTime(0.)
+  fInteractionTime(0.),
+  fEventWeight(1.)
 {
   //
   // Constructor
@@ -43,7 +44,8 @@ AliGenEventHeader::AliGenEventHeader(const char * name):
   TNamed(name, "Event Header"),
   fNProduced(-1),
   fVertex(3),
-  fInteractionTime(0.)
+  fInteractionTime(0.),
+  fEventWeight(1.)
 {
   //
   // Constructor
index 6b060da..9731879 100644 (file)
@@ -23,21 +23,23 @@ class AliGenEventHeader : public TNamed
   AliGenEventHeader();
   virtual ~AliGenEventHeader() {}
   // Getters
-  virtual Int_t           NProduced()  const       {return fNProduced;}
-  virtual void            PrimaryVertex(TArrayF &o) const;
+  virtual Int_t           NProduced()       const  {return fNProduced;}
   virtual Float_t         InteractionTime() const  {return fInteractionTime;}
+  virtual Float_t         EventWeight()     const  {return fEventWeight;}
+  virtual void            PrimaryVertex(TArrayF &o) const;
   // Setters
-  virtual void   SetNProduced(Int_t nprod)         {fNProduced=nprod;}
+  virtual void   SetNProduced(Int_t nprod)         {fNProduced = nprod;}
   virtual void   SetPrimaryVertex(const TArrayF &o);
   virtual void   SetInteractionTime(Float_t t)     {fInteractionTime = t;}
+  virtual void   SetEventWeight(Float_t w)         {fEventWeight = w;}
         
          
 protected:
   Int_t     fNProduced;                 // Number stable or undecayed particles
   TArrayF   fVertex;                    // Primary Vertex Position
   Float_t   fInteractionTime;           // Time of the interaction
-  
-  ClassDef(AliGenEventHeader, 3)        // Event header for primary event
+  Float_t   fEventWeight;               // Event weight
+  ClassDef(AliGenEventHeader, 4)        // Event header for primary event
 };
 
 #endif
index f85b703..3aa0f17 100644 (file)
@@ -83,7 +83,8 @@ AliGenEpEmv1::AliGenEpEmv1():
   fMass(0),
   fEpEmGen(0),
   fDebug(0),
-  fEvent(0)
+  fEvent(0),
+  fHeader(AliGenEventHeader())
 {
   // Default constructor
   // Avoid zero pt
@@ -184,5 +185,7 @@ void AliGenEpEmv1::Generate()
     printf("=====> AliGenEpEmv1::Generate(): \n   Event %d, sigma=%f +- %f kb\n",
           fEvent,fEpEmGen->GetXsection(),fEpEmGen->GetDsection());
   }
+  fHeader.SetEventWeight(weight);
+  AddHeader(&fHeader);
 }
 
index 26c6047..86cca3d 100644 (file)
@@ -12,6 +12,8 @@
 // 9 October 2002
 
 #include "AliGenMC.h"
+#include "AliGenEventHeader.h"
+
 class TEpEmGen;
 
 //-------------------------------------------------------------
@@ -33,6 +35,7 @@ public:
   TEpEmGen * fEpEmGen; // e+e- generator
   Int_t      fDebug;   // debug level
   Int_t      fEvent;   // internal event number
+  AliGenEventHeader fHeader; // Event Header
 ClassDef(AliGenEpEmv1,1) // Generator of single e+e- pair production in PbPb ultra-peripheral collisions
 };
 #endif