UPdated for new README*txt files.

[u/mrichter/AliRoot.git] / ZDC / AliZDCReco.cxx

diff --git a/ZDC/AliZDCReco.cxx b/ZDC/AliZDCReco.cxx
index 431d980e2329cb215e4ca94df77dc619782d2225..62a86bf46b0706e74fa973d3574c6e9379450964 100644 (file)
--- a/ZDC/AliZDCReco.cxx
+++ b/ZDC/AliZDCReco.cxx
@@ -27,32 +27,145 @@
 
 ClassImp(AliZDCReco)
   
+
 //_____________________________________________________________________________
-AliZDCReco::AliZDCReco(Float_t ezn1, Float_t ezp1, Float_t ezdc1, Float_t ezem,
-       Float_t ezn2, Float_t ezp2, Float_t ezdc2,
-       Int_t detspnLeft, Int_t detsppLeft, Int_t detspnRight, Int_t detsppRight,
-       Int_t trspn, Int_t trspp, Int_t trsp, Int_t part, Float_t b)
+AliZDCReco::AliZDCReco() :
+       
+  TObject(),
+  fZN1Energy(0),
+  fZP1Energy(0),
+  fZN2Energy(0),
+  fZP2Energy(0),
+  //
+  fZEM1signal(0),
+  fZEM2signal(0),
+  //
+  fNDetSpecNLeft(0),
+  fNDetSpecPLeft(0),
+  fNDetSpecNRight(0),
+  fNDetSpecPRight(0),
+  fNTrueSpecNLeft(0),
+  fNTrueSpecPLeft(0),
+  fNTrueSpecLeft(0),
+  fNTrueSpecNRight(0),
+  fNTrueSpecPRight(0),
+  fNTrueSpecRight(0),
+  fNPartLeft(0),
+  fNPartRight(0),
+  fImpPar(0)
+
 { 
   //
-  // Standard constructor
+  // Default constructor
   //
-  fZN1energy  = ezn1;
-  fZP1energy  = ezp1;
-  fZDC1energy = ezdc1;
-  fZN2energy  = ezn2;
-  fZP2energy  = ezp2;
-  fZDC2energy = ezdc2;
-  fZEMenergy = ezem;
-  fNDetSpecNLeft = detspnLeft;
-  fNDetSpecPLeft = detsppLeft;
-  fNDetSpecNRight = detspnRight;
-  fNDetSpecPRight = detsppRight;
-  fNTrueSpecN = trspn;
-  fNTrueSpecP = trspp;
-  fNTrueSpec = trsp;
-  fNPart     = part;
-  fImpPar    = b;
+  for(Int_t i=0; i<5; i++){
+     fZN1EnTow[i] = 0;
+     fZP1EnTow[i] = 0;
+     fZN2EnTow[i] = 0;
+     fZP2EnTow[i] = 0;
+     fZN1SigLowRes[i] = 0;
+     fZP1SigLowRes[i] = 0;
+     fZN2SigLowRes[i] = 0;
+     fZP2SigLowRes[i] = 0;
+     
+  }
+}
   
+
+//_____________________________________________________________________________
+AliZDCReco::AliZDCReco(Float_t ezn1, Float_t ezp1, Float_t ezn2, Float_t ezp2,  
+            //
+            Float_t* ezn1tow, Float_t* ezp1tow,
+            Float_t* ezn2tow, Float_t* ezp2tow, 
+            Float_t* ezn1siglr, Float_t* ezp1siglr,
+            Float_t* ezn2siglr, Float_t* ezp2siglr,
+            Float_t ezem1, Float_t ezem2,
+            //    
+            Int_t detspnLeft,  Int_t detsppLeft, Int_t detspnRight,
+            Int_t detsppRight,  Int_t trspnLeft, Int_t trsppLeft, 
+            Int_t trspLeft, Int_t partLeft, Int_t trspnRight, 
+            Int_t trsppRight, Int_t trspRight, Int_t partRight,  
+            Float_t b) :
+       
+  TObject(),
+  fZN1Energy(ezn1),
+  fZP1Energy(ezp1),
+  fZN2Energy(ezn2),
+  fZP2Energy(ezp2),
+  //
+  fZEM1signal(ezem1),
+  fZEM2signal(ezem2),
+  //
+  fNDetSpecNLeft(detspnLeft),
+  fNDetSpecPLeft(detsppLeft),
+  fNDetSpecNRight(detspnRight),
+  fNDetSpecPRight(detsppRight),
+  fNTrueSpecNLeft(trspnLeft),
+  fNTrueSpecPLeft(trsppLeft),
+  fNTrueSpecLeft(trspLeft),
+  fNTrueSpecNRight(trspnRight),
+  fNTrueSpecPRight(trsppRight),
+  fNTrueSpecRight(trspRight),
+  fNPartLeft(partLeft),
+  fNPartRight(partRight),
+  fImpPar(b)
+
+{ 
+  //
+  // Constructor
+  //
+  for(Int_t j=0; j<5; j++){
+     fZN1EnTow[j] =  ezn1tow[j];
+     fZP1EnTow[j] =  ezp1tow[j];
+     fZN2EnTow[j] =  ezn2tow[j];
+     fZP2EnTow[j] =  ezp2tow[j];
+     fZN1SigLowRes[j] = ezn1siglr[j];
+     fZP1SigLowRes[j] = ezp1siglr[j];
+     fZN2SigLowRes[j] = ezn2siglr[j];
+     fZP2SigLowRes[j] = ezp2siglr[j];
+  }
+  
+}
+
+//______________________________________________________________________________
+AliZDCReco::AliZDCReco(const AliZDCReco &oldreco) :
+
+  TObject()
+{
+  // Copy constructor
+
+  fZN1Energy  = oldreco.GetZN1Energy();
+  fZP1Energy  = oldreco.GetZP1Energy();            
+  fZN2Energy  = oldreco.GetZN2Energy();           
+  fZP2Energy  = oldreco.GetZP2Energy();           
+  //
+  for(Int_t i=0; i<5; i++){      
+     fZN1EnTow[i]  = oldreco.GetZN1EnTow(i);
+     fZP1EnTow[i]  = oldreco.GetZP1EnTow(i);
+     fZN2EnTow[i]  = oldreco.GetZN2EnTow(i);
+     fZP2EnTow[i]  = oldreco.GetZP2EnTow(i);
+     fZN1SigLowRes[i] = oldreco.GetZN1SigLowRes(i);
+     fZP1SigLowRes[i] = oldreco.GetZP1SigLowRes(i);
+     fZN2SigLowRes[i] = oldreco.GetZN2SigLowRes(i);
+     fZP2SigLowRes[i] = oldreco.GetZP2SigLowRes(i);
+  }
+  //
+  fZEM1signal = oldreco.GetZEM1signal();       
+  fZEM2signal = oldreco.GetZEM2signal();       
+  //   
+  fNDetSpecNLeft = oldreco.GetNDetSpecNLeft(); 
+  fNDetSpecPLeft = oldreco.GetNDetSpecPLeft(); 
+  fNDetSpecNRight = oldreco.GetNDetSpecNRight();       
+  fNDetSpecPRight = oldreco.GetNDetSpecPRight();       
+  fNTrueSpecNLeft = oldreco.GetNTrueSpecNLeft();         
+  fNTrueSpecPLeft = oldreco.GetNTrueSpecPLeft();         
+  fNTrueSpecLeft = oldreco.GetNTrueSpecLeft();   
+  fNTrueSpecNRight = oldreco.GetNTrueSpecNRight();       
+  fNTrueSpecPRight = oldreco.GetNTrueSpecPRight();       
+  fNTrueSpecRight = oldreco.GetNTrueSpecRight();         
+  fNPartLeft = oldreco.GetNPartLeft();                  
+  fNPartRight = oldreco.GetNPartRight();                        
+  fImpPar = oldreco.GetImpPar();                        
 }
 
 //______________________________________________________________________________
@@ -60,9 +173,15 @@ void AliZDCReco::Print(Option_t *) const {
   //
   // Printing Reconstruction Parameters
   //
-  printf("     ---   Reconstruction -> EZN = %f TeV, EZP = %f TeV, EZDC = %f TeV,"
-        " EZEM = %f GeV \n             NDetSpecN = %d, NDetSpecP = %d, Nspecn = %d,"
-        " Nspecp = %d, Npart = %d, b = %f fm.\n ", 
-        fZN1energy,fZP1energy,fZDC1energy,fZEMenergy,fNDetSpecNLeft,
-        fNDetSpecPLeft,fNTrueSpecN,fNTrueSpecP,fNPart,fImpPar);
+  printf(" \t ---   Reconstruction -> EZN1 = %f TeV, EZP1 = %f TeV,  EZEM1 = %f GeV ,  EZEM2 = %f GeV \n "             
+        "EZN2 = %f TeV, EZP2 = %f TeV \n"
+        " \t NDetSpecNLeft = %d, NDetSpecPLeft = %d, NspecnLeft = %d,"
+        " NspecpLeft = %d, NpartLeft = %d"
+        " \t NDetSpecNRight = %d, NDetSpecPRight = %d, NspecnRight = %d,"
+        " NspecpRight = %d, NpartRight = %d"
+        " \t b = %f fm\n ", 
+        fZN1Energy,fZP1Energy,fZEM1signal,fZEM2signal, fZN2Energy, fZP2Energy,
+        fNDetSpecNLeft,fNDetSpecPLeft,fNTrueSpecNLeft,fNTrueSpecPLeft,fNPartLeft,
+        fNDetSpecNRight,fNDetSpecPRight,fNTrueSpecNRight,fNTrueSpecPRight,fNPartRight,
+        fImpPar);
 }