]> git.uio.no Git - u/mrichter/AliRoot.git/blame_incremental - EMCAL/AliEMCALRecPoint.h
git://git.uio.no / u / mrichter / AliRoot.git / blame_incremental
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (non-incremental) | history | HEAD
Full set of updated misalignment macros (Raffaele)
[u/mrichter/AliRoot.git] / EMCAL / AliEMCALRecPoint.h
... / ...
CommitLineData
1#ifndef ALIEMCALRECPOINT_H
2#define ALIEMCALRECPOINT_H
3/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
5//_________________________________________________________________________
6// Base Class for EMCAL Reconstructed Points
7// A recpoint being equivalent to a cluster in encal terminology
8//*-- Author: Yves Schutz (SUBATECH)
9//*-- Author: Dmitri Peressounko (RRC KI & SUBATECH)
10//*-- Author: Heather Gray (LBL): merged AliEMCALRecPoint and AliEMCALTowerRecPoint 02/04
11
12// --- ROOT system ---
13class TVector3 ;
14
15// --- Standard library ---
16
17// --- AliRoot header files ---
18
19#include "AliRecPoint.h"
20class AliEMCALDigit;
21class AliEMCALGeometry;
22
23class AliEMCALRecPoint : public AliRecPoint {
24
25 public:
26
27 typedef TObjArray RecPointsList ;
28
29 AliEMCALRecPoint() ; // ctor
30 AliEMCALRecPoint(const char * opt) ; // ctor
31 AliEMCALRecPoint(const AliEMCALRecPoint & rp);
32
33 virtual ~AliEMCALRecPoint();
34 virtual void AddDigit(AliDigitNew &){ Fatal("AddDigit", "use AddDigit(AliEMCALDigit & digit, Float_t Energy )") ; }
35 virtual void AddDigit(AliEMCALDigit & digit, const Float_t Energy);
36 virtual Int_t Compare(const TObject * obj) const;
37 virtual Int_t DistancetoPrimitive(Int_t px, Int_t py);
38 virtual void Draw(Option_t * option="") ;
39 virtual void ExecuteEvent(Int_t event, Int_t, Int_t) ;
40
41 virtual void SetClusterType(Int_t ver) { fClusterType = ver; }
42 virtual Int_t GetClusterType() const { return fClusterType; }
43
44 virtual void EvalAll(Float_t logWeight, TClonesArray * digits);
45 virtual void EvalLocalPosition(Float_t logWeight, TClonesArray * digits) ;
46 // void EvalLocalPositionSimple(TClonesArray *digits); // ??
47 virtual void EvalPrimaries(TClonesArray * digits) ;
48 virtual void EvalParents(TClonesArray * digits) ;
49
50 using AliRecPoint::GetGlobalPosition;
51 virtual void GetGlobalPosition(TVector3 & gpos) const; // return global position (x, y, z) in ALICE
52 virtual void GetLocalPosition(TVector3 & lpos) const; // return local position (x, y, z) in EMCAL SM
53 virtual Int_t * GetPrimaries(Int_t & number) const {number = fMulTrack ;
54 return fTracksList ; }
55 virtual Int_t * GetParents(Int_t & number) const {number = fMulParent ;
56 return fParentsList ; }
57 Float_t GetCoreEnergy()const {return fCoreEnergy ;}
58 virtual Float_t GetDispersion()const {return fDispersion ;}
59 virtual void GetElipsAxis(Float_t * lambda)const {lambda[0] = fLambda[0]; lambda[1] = fLambda[1];};
60
61 Float_t * GetEnergiesList() const {return fEnergyList ;} // gets the list of energies making this recpoint
62 Double_t GetPointEnergy() const; // gets point energy (sum of energy list)
63 Float_t * GetTimeList() const {return fTimeList ;} // gets the list of digit times in this recpoint
64 Float_t GetMaximalEnergy(void) const ; // get the highest energy in the cluster
65 Int_t GetMaximumMultiplicity() const {return fMaxDigit ;} // gets the maximum number of digits allowed
66 Int_t GetMultiplicity(void) const { return fMulDigit ; } // gets the number of digits making this recpoint
67 Int_t GetMultiplicityAtLevel(Float_t level) const ; // computes multiplicity of digits with
68 Int_t * GetAbsId() const {return fAbsIdList;}
69 Int_t GetAbsId(int i) const {if(i>=0 && i<fMulDigit)return fAbsIdList[i]; else return -1;}
70 // energy above relative level
71 virtual Int_t GetNumberOfLocalMax(AliEMCALDigit ** maxAt, Float_t * maxAtEnergy,
72 Float_t locMaxCut,TClonesArray * digits ) const ;
73 // searches for the local maxima
74
75 Int_t GetPrimaryIndex() const ;
76 Float_t GetTime(void) const{return fTime ; }
77
78 virtual Bool_t IsEmc(void)const { return kTRUE ; }
79 virtual Bool_t IsSortable() const {
80 // tells that this is a sortable object
81 return kTRUE ;
82 }
83 virtual void Paint(Option_t * option="");
84 virtual void Print(Option_t * option="") const ;
85
86 AliEMCALRecPoint & operator = (const AliEMCALRecPoint & ) {
87 Fatal("operator =", "not implemented") ;
88 return *this ;
89 }
90
91protected:
92 void EvalCoreEnergy(Float_t logWeight,TClonesArray * digits) ;
93 virtual void EvalDispersion(Float_t logWeight,TClonesArray * digits) ; // computes the dispersion of the shower
94 virtual void EvalElipsAxis(Float_t logWeight, TClonesArray * digits ); // computes the axis of shower ellipsoide
95 void EvalTime( TClonesArray * digits );
96 virtual Bool_t AreNeighbours(AliEMCALDigit * digit1, AliEMCALDigit * digit2 ) const;
97 Float_t ThetaToEta(Float_t arg) const; //Converts Theta (Radians) to Eta(Radians)
98 Float_t EtaToTheta(Float_t arg) const; //Converts Eta (Radians) to Theta(Radians)
99
100private:
101 AliEMCALGeometry* fGeomPtr; //! Pointer to geometry for utilities
102
103 Int_t fClusterType; // type of cluster stored:
104 // pseudocluster or v1
105 Float_t fCoreEnergy ; // energy in a shower core
106 Float_t fLambda[2] ; // shower ellipse axes
107 Float_t fDispersion ; // shower dispersio
108 Float_t *fEnergyList ; //[fMulDigit] energy of digits
109 Float_t *fTimeList ; //[fMulDigit] time of digits
110 Int_t *fAbsIdList; //[fMulDigit] absId of digits
111 Float_t fTime ; // Time of the digit with maximal energy deposition
112 Float_t fCoreRadius; // The radius in which the core energy is evaluated
113 Int_t fMulParent; // Multiplicity of the parents
114 Int_t fMaxParent; // Maximum number of parents allowed
115 Int_t * fParentsList; // [fMulParent] list of the parents of the digits
116 Int_t fSuperModuleNumber; // number identifying supermodule containing recpoint
117
118 ClassDef(AliEMCALRecPoint,7) // RecPoint for EMCAL (Base Class)
119
120};
121
122#endif // AliEMCALRECPOINT_H