]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - EMCAL/AliEMCALRecPoint.h
To create back-to-back partons we need to call the function pyjoin from Pythia.
[u/mrichter/AliRoot.git] / EMCAL / AliEMCALRecPoint.h
index faac51f1f2ee3769b800a7512a0d0dc977b3b051..26e61b82cff6b0a00b14b8d7e21f1ac850e3465b 100644 (file)
@@ -47,17 +47,19 @@ class AliEMCALRecPoint : public AliCluster {
   virtual void    AddDigit(AliDigitNew &) const { 
     Fatal("AddDigit", "use AddDigit(AliEMCALDigit & digit, Float_t Energy )") ; 
   }
-  virtual void    AddDigit(AliEMCALDigit & digit, Float_t Energy); 
-  virtual Int_t   Compare(const TObject * obj) const;   
-  virtual Int_t   DistancetoPrimitive(Int_t px, Int_t py);
+  virtual void    AddDigit(AliEMCALDigit & digit, Float_t Energy, Bool_t shared); 
+  virtual Int_t   Compare(const TObject * obj) const;  
+  //virtual Int_t   DistancetoPrimitive(Int_t px, Int_t py);//Not used, remove?
   virtual void    Draw(Option_t * option="") ;
-  virtual void    ExecuteEvent(Int_t event, Int_t, Int_t) ;
+  //virtual void    ExecuteEvent(Int_t event, Int_t, Int_t) ;//Not used, remove?
 
   virtual void    SetClusterType(Int_t ver) { fClusterType = ver; }
   virtual Int_t   GetClusterType() const { return fClusterType; }
 
-  virtual void    EvalAll(Float_t logWeight, TClonesArray * digits);
-  virtual void    EvalLocalPosition(Float_t logWeight, TClonesArray * digits);
+  virtual void    EvalAll(Float_t logWeight, TClonesArray * digits, const Bool_t justClusters);
+  virtual void    EvalLocalPosition (Float_t logWeight, TClonesArray * digits);
+  virtual void    EvalGlobalPosition(Float_t logWeight, TClonesArray * digits);
+
   virtual void    EvalPrimaries(TClonesArray * digits) ;
   virtual void    EvalParents(TClonesArray * digits) ;
 
@@ -68,21 +70,21 @@ class AliEMCALRecPoint : public AliCluster {
   void   EvalLocalPositionFit(Double_t deff, Double_t w0, Double_t phiSlope,TClonesArray * digits);
   Bool_t EvalLocalPosition2(TClonesArray *digits, TArrayD &ed);
 
-  static  Bool_t  EvalLocalPositionFromDigits(const Double_t esum, const Double_t deff, const Double_t w0, 
+  Bool_t EvalLocalPositionFromDigits(const Double_t esum, const Double_t deff, const Double_t w0, 
                                               TClonesArray *digits, TArrayD &ed, TVector3 &locPos);
-  static  Bool_t  EvalLocalPositionFromDigits(TClonesArray *digits, TArrayD &ed, TVector3 &locPos);
+  Bool_t EvalLocalPositionFromDigits(TClonesArray *digits, TArrayD &ed, TVector3 &locPos);
   static  void    GetDeffW0(const Double_t esum, Double_t &deff,  Double_t &w0);
 
-  virtual void    GetGlobalPosition(TVector3 & gpos, TMatrixF & gmat) const; // return global position (x, y, z) in ALICE
+  //virtual void    GetGlobalPosition(TVector3 & gpos, TMatrixF & gmat) const; // return global position (x, y, z) in ALICE
   virtual void    GetGlobalPosition(TVector3 & gpos) const; // return global position (x, y, z) in ALICE
-  virtual void    GetLocalPosition(TVector3 & lpos) const;  // return local position  (x, y, z) in EMCAL SM
+  virtual void    GetLocalPosition(TVector3 & lpos)  const;  // return local position  (x, y, z) in EMCAL SM
   virtual Int_t * GetPrimaries(Int_t & number) const {number = fMulTrack ; 
                                                       return fTracksList ; }
-  virtual Int_t * GetParents(Int_t & number) const {number = fMulParent ; 
+  virtual Int_t * GetParents(Int_t & number)   const {number = fMulParent ; 
                                                       return fParentsList ; }
 
-  virtual Int_t   GetDigitsMultiplicity(void) const { return fMulDigit ; }
-  Int_t           GetIndexInList() const { return fIndexInList ; }
+  virtual Int_t   GetDigitsMultiplicity(void)  const { return fMulDigit ; }
+  Int_t           GetIndexInList()             const { return fIndexInList ; }
 
   Float_t         GetCoreEnergy()const {return fCoreEnergy ;}
   virtual Float_t GetDispersion()const {return fDispersion ;}
@@ -90,13 +92,14 @@ class AliEMCALRecPoint : public AliCluster {
   
   Float_t *   GetEnergiesList() const {return fEnergyList ;}       // gets the list of energies making this recpoint
   Double_t    GetPointEnergy() const;                              // gets point energy (sum of energy list)
-  Float_t *   GetTimeList() const {return fTimeList ;}       // gets the list of digit times in this recpoint
+  Float_t *   GetTimeList() const {return fTimeList ;}             // gets the list of digit times in this recpoint
   Float_t     GetMaximalEnergy(void) const ;                       // get the highest energy in the cluster
+  Int_t       GetMaximalEnergyIndex(void) const ;                  // get the index of highest energy digit
   Int_t       GetMaximumMultiplicity() const {return fMaxDigit ;}  // gets the maximum number of digits allowed
   Int_t       GetMultiplicity(void) const { return fMulDigit ; }   // gets the number of digits making this recpoint
-  Int_t       GetMultiplicityAtLevel(Float_t level) const ;  // computes multiplicity of digits with 
+  Int_t       GetMultiplicityAtLevel(Float_t level) const ;        // computes multiplicity of digits with 
   Int_t *     GetAbsId() const {return fAbsIdList;}
-  Int_t       GetAbsId(int i) const {if(i>=0 && i<fMulDigit)return fAbsIdList[i]; else return -1;}
+  Int_t       GetAbsId(Int_t i) const {if(i>=0 && i<fMulDigit)return fAbsIdList[i]; else return -1;}
   Int_t       GetAbsIdMaxDigit() const {return GetAbsId(fDigitIndMax);}
   Int_t       GetIndMaxDigit() const {return fDigitIndMax;}
   void        SetIndMaxDigit(const Int_t ind) {fDigitIndMax = ind;}
@@ -111,20 +114,25 @@ class AliEMCALRecPoint : public AliCluster {
   // Number of local maxima found in cluster in unfolding:
   // 0: no unfolding
   //-1: unfolding failed
-  Short_t     GetNExMax(void) const {return fNExMax ;}  // Number of maxima found in cluster in unfolding
-  void        SetNExMax(Int_t nmax=1){fNExMax = static_cast<Short_t>(nmax) ;}
-  Int_t       GetPrimaryIndex() const  ;
-  Float_t     GetTime(void) const{return  fTime ; }
-  virtual Bool_t  IsEmc(void)const { return kTRUE ;  }
-  virtual Bool_t  IsSortable() const { 
+  Short_t     GetNExMax(void)     const { return fNExMax ;}  // Number of maxima found in cluster in unfolding
+  void        SetNExMax(Int_t nmax=1)   {fNExMax = static_cast<Short_t>(nmax) ;}
+       
+  Int_t       GetPrimaryIndex()   const  ;
+       
+  Float_t     GetTime(void)       const { return  fTime ; }
+       
+  Bool_t      SharedCluster(void) const { return  fSharedCluster ; }
+  void        SetSharedCluster(Bool_t s){ fSharedCluster = s ; }
+       
+  virtual Bool_t  IsEmc(void)     const { return kTRUE ;  }
+  virtual Bool_t  IsSortable()    const { 
     // tells that this is a sortable object
     return kTRUE ; 
   }  
   virtual void    Paint(Option_t * option="");
   virtual void    Print(Option_t * option="") const ; 
   
-  static Double_t TmaxInCm(const Double_t e=0.0, const Int_t key=0);
+  Double_t TmaxInCm(const Double_t e=0.0, const Int_t key=0);
 
   Float_t     GetDistanceToBadTower() const {return fDistToBadTower;}
   void        EvalDistanceToBadChannels(AliCaloCalibPedestal* caloped);
@@ -146,8 +154,8 @@ private:
          Float_t fAmp ;            // summed amplitude of digits   
          Int_t   fIndexInList ;    // the index of this RecPoint in the
                                    // list stored in TreeR (to be set by analysis)
-         TVector3    fLocPos ;     // local position in the sub-detector coordinate
-         TMatrixF *  fLocPosM ;    // covariance matrix ;
+         TVector3    fGlobPos ;    // global position
+         TVector3    fLocPos ;     // local  position in the sub-detector coordinate
          Int_t       fMaxDigit ;   //! max initial size of digits array (not saved)
          Int_t       fMulDigit ;   // total multiplicity of digits       
          Int_t       fMaxTrack ;   //! max initial size of tracks array (not saved)
@@ -155,7 +163,7 @@ private:
          Int_t *     fDigitsList ; //[fMulDigit] list of digit's indexes from which the point was reconstructed
          Int_t *     fTracksList ; //[fMulTrack] list of tracks to which the point was assigned
 
-         Int_t   fClusterType;    // type of cluster stored: pseudocluster or v1
+         Int_t   fClusterType;       // type of cluster stored: v1
          Float_t fCoreEnergy ;       // energy in a shower core 
          Float_t fLambda[2] ;        // shower ellipse axes
          Float_t fDispersion ;       // shower dispersion
@@ -170,11 +178,11 @@ private:
          Int_t fMaxParent;           // Maximum number of parents allowed
          Int_t * fParentsList;       // [fMulParent] list of the parents of the digits
          Float_t * fDEParentsList;   // [fMulParent] list of the parents of the digits
-         Int_t   fSuperModuleNumber; // number identifying supermodule containing recpoint
-         // Aug 16, 2007
+         Int_t   fSuperModuleNumber; // number identifying supermodule containing recpoint, reference is cell with maximum energy.
          Int_t   fDigitIndMax;       // Index of digit with max energy in array fAbsIdList
-         Float_t fDistToBadTower;  // Distance to nearest bad tower
-
+         Float_t fDistToBadTower;    // Distance to nearest bad tower
+         Bool_t  fSharedCluster;     // States if cluster is shared by 2 SuperModules in same phi rack (0,1), (2,3) ... (10,11).
+       
   ClassDef(AliEMCALRecPoint,12) // RecPoint for EMCAL (Base Class)
  
 };