I added some methods that return squared distances in order to avoid sqrt calculation...

[u/mrichter/AliRoot.git] / STEER / AliGenPythiaEventHeader.cxx

diff --git a/STEER/AliGenPythiaEventHeader.cxx b/STEER/AliGenPythiaEventHeader.cxx
index a203bfbc5cb14e6ea4f347fc968e496d6c044330..e5be04ab321fcb95bd4362a76ec7eabf2124c020 100644 (file)
--- a/STEER/AliGenPythiaEventHeader.cxx
+++ b/STEER/AliGenPythiaEventHeader.cxx
@@ -26,10 +26,17 @@ AliGenPythiaEventHeader::AliGenPythiaEventHeader():
     fNUQJets(0), 
     fXJet(-1.),
     fYJet(-1.),
+    fInMediumLength(0.),
+    fImpactParameter(-1),
     fPtHard(0.) 
 {
 // Default Constructor
     for (Int_t i = 0; i < 4; i++) fZquench[i] = 0.;
+    for (Int_t i = 0; i < 4; i++)
+      for (Int_t j = 0; j < 10; j++) {
+       fJets[i][j] = 0.;
+       fUQJets[i][j] = 0.;
+      }
 }
 
 AliGenPythiaEventHeader::AliGenPythiaEventHeader(const char* name):
@@ -40,10 +47,17 @@ AliGenPythiaEventHeader::AliGenPythiaEventHeader(const char* name):
     fNUQJets(0), 
     fXJet(-1.),
     fYJet(-1.),
+    fInMediumLength(0.),
+    fImpactParameter(-1),
     fPtHard(0.) 
 {
 // Constructor
     for (Int_t i = 0; i < 4; i++) fZquench[i] = 0.;
+    for (Int_t i = 0; i < 4; i++)
+      for (Int_t j = 0; j < 10; j++) {
+       fJets[i][j] = 0.;
+       fUQJets[i][j] = 0.;
+      }
 }
 
 void AliGenPythiaEventHeader::AddJet(Float_t px, Float_t py, Float_t pz, Float_t e)