coding convention

[u/mrichter/AliRoot.git] / PHOS / AliPHOSEmcRecPoint.h

diff --git a/PHOS/AliPHOSEmcRecPoint.h b/PHOS/AliPHOSEmcRecPoint.h
index e3761a9dcaf21ad9db3c4f17f2cd34dab1ff14dd..bdfe3c59c25128f02e637f3c72e7cea8b7971635 100644 (file)
--- a/PHOS/AliPHOSEmcRecPoint.h
+++ b/PHOS/AliPHOSEmcRecPoint.h
@@ -13,14 +13,11 @@
 
 // --- ROOT system ---
 
-//#include "TObject.h"
-#include "TArrayI.h"
- 
 // --- Standard library ---
 
 // --- AliRoot header files ---
 
-#include "AliPHOSDigit.h"
+//#include "AliPHOSDigit.h"
 #include "AliPHOSRecPoint.h"
 
 class AliPHOSEmcRecPoint : public AliPHOSRecPoint  {
@@ -29,11 +26,7 @@ public:
 
   AliPHOSEmcRecPoint() ;
   AliPHOSEmcRecPoint(const char * opt) ;
-  AliPHOSEmcRecPoint(const AliPHOSEmcRecPoint & rp) {
-    // cpy ctor requested by Coding Convention 
-    // but not yet needed
-    assert(0==1) ; 
-  } 
+  AliPHOSEmcRecPoint(const AliPHOSEmcRecPoint & rp) ; 
  
   virtual ~AliPHOSEmcRecPoint() ;  
 
@@ -54,6 +47,10 @@ public:
   Int_t       GetMultiplicity(void) const { return fMulDigit ; }   // gets the number of digits making this recpoint
   Int_t       GetMultiplicityAtLevel(const Float_t level) const ;  // computes multiplicity of digits with 
                                                                    // energy above relative level
+  Short_t     GetNExMax(void) const {return fNExMax ;}             // Number of maxima found in cluster in unfolding:
+                                                                   // 0: was no unfolging
+                                                                   //-1: unfolding failed
+  void        SetNExMax(Int_t nmax=1){fNExMax = static_cast<Short_t>(nmax) ;}
   virtual Int_t GetNumberOfLocalMax(AliPHOSDigit **  maxAt, Float_t * maxAtEnergy,
                                     Float_t locMaxCut,TClonesArray * digits ) const ; 
                                                                    // searches for the local maxima 
@@ -61,20 +58,23 @@ public:
   Float_t     GetTime(void) const{return  fTime ; } 
   Bool_t      IsEmc(void) const { return kTRUE ; }                 // true if the recpoint is in EMC
   Bool_t      IsSortable() const {return kTRUE ; }                 // says that emcrecpoints are sortable objects 
-  void        Print(Option_t * opt = "void")const ; 
+  void        Print(Option_t *)const ; 
   void        Purify(Float_t threshold) ;                          //Removes digits below threshold
 
-  AliPHOSEmcRecPoint & operator = (const AliPHOSEmcRecPoint & rvalue)  {
-    // assignement operator requested by coding convention but not needed
-    assert(0==1) ;
-    return *this ; 
-  }
+  Float_t     GetM2x()   const {return fM2x;  } // Get second X-moment
+  Float_t     GetM2z()   const {return fM2z;  } // Get second Z-moment
+  Float_t     GetM3x()   const {return fM3x;  } // Get third  X-moment
+  Float_t     GetM4z()   const {return fM4z;  } // Get forth  Z-moment
+  Float_t     GetPhixe() const {return fPhixe;} // Get angle between center gravity and eigen vector
+
+  AliPHOSEmcRecPoint & operator = (const AliPHOSEmcRecPoint & /*rvalue*/)  { return *this ; }
 
  protected:
           void  EvalCoreEnergy(Float_t logWeight,TClonesArray * digits) ;             
   virtual void  EvalLocalPosition(Float_t logWeight,TClonesArray * digits) ;// computes the position in the PHOS module 
   virtual void  EvalDispersion(Float_t logWeight,TClonesArray * digits) ;   // computes the dispersion of the shower
   virtual void  EvalElipsAxis(Float_t logWeight, TClonesArray * digits );   // computes the axis of shower ellipsoide
+          void  EvalMoments(Float_t logWeight, TClonesArray * digits );     // computes shower moments
           void  EvalTime( TClonesArray * digits );
   virtual Bool_t AreNeighbours(AliPHOSDigit * digit1, AliPHOSDigit * digit2 ) const ;
 
@@ -83,8 +83,17 @@ public:
   Float_t fDispersion ;       // shower dispersion
   Float_t *fEnergyList ;      //[fMulDigit] energy of digits
   Float_t fTime ;             // Time of the digit with maximal energy deposition
+  Short_t fNExMax ;           // number of (Ex-)maxima before unfolding
+
+  Float_t fM2x;               // Second moment along X axis
+  Float_t fM2z;               // Second moment along Z axis
+  Float_t fM3x;               // Third  moment along X axis
+  Float_t fM4z;               // Forth  moment along Z axis
+  Float_t fPhixe;             // Angle between center-gravity vector and eigen vector
+
+  Int_t fDebug;               //! debug level (0 - no output)
   
-  ClassDef(AliPHOSEmcRecPoint,1)  // EMC RecPoint (cluster)
+  ClassDef(AliPHOSEmcRecPoint,2)  // EMC RecPoint (cluster)
 
 };