True npart interfaced to HIJING (header)
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 13 Sep 2012 16:09:30 +0000 (16:09 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 13 Sep 2012 16:09:30 +0000 (16:09 +0000)
HIJING/hijing/himain1.inc
STEER/STEERBase/AliGenHijingEventHeader.cxx
STEER/STEERBase/AliGenHijingEventHeader.h
THijing/AliGenHijing.cxx
THijing/Hcommon.h
THijing/THijing.cxx
THijing/THijing.h

index 4ff6f13..12a6dbb 100644 (file)
@@ -1,6 +1,6 @@
 * $Id$
         REAL EATT, BB
-       COMMON/HIMAIN1/NATT,EATT,JATT,NT,NP,N0,N01,N10,N11,BB
+       COMMON/HIMAIN1/NATT,EATT,JATT,NT,NP,N0,N01,N10,N11,BB,NPART
 #ifndef BLANKET_SAVE
        SAVE  /HIMAIN1/
 #endif
index 4873513..887b944 100644 (file)
@@ -21,6 +21,7 @@ ClassImp(AliGenHijingEventHeader)
 AliGenHijingEventHeader::AliGenHijingEventHeader():
     fTotalEnergy(0.),
     fTrials(0),
+    fNPart(0),
     fJet1(0., 0., 0., 0.),
     fJet2(0., 0., 0., 0.),
     fJetFsr1(0., 0., 0., 0.),
@@ -33,6 +34,7 @@ AliGenHijingEventHeader::AliGenHijingEventHeader(const char* name):
     AliGenEventHeader(name),
     fTotalEnergy(0.),
     fTrials(0),
+    fNPart(0),
     fJet1(0., 0., 0., 0.),
     fJet2(0., 0., 0., 0.),
     fJetFsr1(0., 0., 0., 0.),
index f92ea57..374699e 100644 (file)
@@ -20,8 +20,8 @@ class AliGenHijingEventHeader : public AliGenEventHeader, public AliCollisionGeo
   // Getters
   Float_t TotalEnergy() const {return fTotalEnergy;} 
   Int_t   Trials() const {return fTrials;}
-  
-         
+  Int_t   GetTrueNPart() const {return fNPart;}
+         
   // Setters
   void SetTotalEnergy(Float_t energy)  {fTotalEnergy=energy;}
   void SetJets(const TLorentzVector* jet1, const TLorentzVector* jet2,
@@ -31,11 +31,12 @@ class AliGenHijingEventHeader : public AliGenEventHeader, public AliCollisionGeo
               TLorentzVector& jet3, TLorentzVector& jet4) const  
       {jet1 = fJet1; jet2 = fJet2; jet3 = fJetFsr1; jet4 = fJetFsr2;}
   void SetTrials(Int_t trials) {fTrials = trials;}
-         
+  void SetTrueNPart(Int_t npart) {fNPart = npart;} 
 protected:
   Float_t fTotalEnergy;              // Total energy of produced particles
   Int_t   fTrials;                   // Number of trials to fulfill trigger condition
-  
+  Int_t   fNPart;                    // True number of participants 
   TLorentzVector  fJet1;             // 4-Momentum-Vector of first   triggered jet  
   TLorentzVector  fJet2;             // 4-Momentum-Vector of second  triggered jet     
   TLorentzVector  fJetFsr1;          // 4-Momentum-Vector of first   triggered jet  
index 915e8cf..aa1a28f 100644 (file)
@@ -563,6 +563,7 @@ void AliGenHijing::MakeHeader()
     fHeader.SetSpectators(fProjectileSpecn, fProjectileSpecp,
                          fTargetSpecn,fTargetSpecp);
     fHeader.SetReactionPlaneAngle(fHijing->GetHINT1(20));
+    fHeader.SetTrueNPart(fHijing->GetNPART());
 
 // 4-momentum vectors of the triggered jets.
 //
index 9797685..d763ec1 100644 (file)
@@ -321,7 +321,8 @@ typedef struct {
    Int_t      n01;
    Int_t      n10;
    Int_t      n11;
-   Float_t    bb;   
+   Float_t    bb;
+   Int_t      npart;
    } Himain1Common;
 
 #define HIMAIN1 COMMON_BLOCK(HIMAIN1,himain1)
index 5868245..a0f2fe1 100644 (file)
@@ -506,6 +506,13 @@ Int_t  THijing::GetNATT() const
 
 }
 
+Int_t  THijing::GetNPART() const
+{
+// Get true number of participants
+   return HIMAIN1.npart;
+
+}
+
 //______________________________________________________________________________
 Float_t  THijing::GetEATT() const
 {
index 58420aa..532d8b0 100644 (file)
@@ -82,6 +82,8 @@ public:
    
    virtual Int_t       GetNATT() const;
 
+   virtual Int_t       GetNPART() const;
+
    virtual Float_t     GetEATT() const;
 
    virtual Int_t       GetJATT() const;