| | 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 --- |
| | 13 | class TVector3 ; |
| | 14 | |
| | 15 | // --- Standard library --- |
| | 16 | |
| | 17 | // --- AliRoot header files --- |
| | 18 | |
| | 19 | #include "AliRecPoint.h" |
| | 20 | #include "AliEMCALDigit.h" |
| | 21 | |
| | 22 | class AliEMCALRecPoint : public AliRecPoint { |
| | 23 | |
| | 24 | public: |
| | 25 | |
| | 26 | typedef TObjArray RecPointsList ; |
| | 27 | |
| | 28 | AliEMCALRecPoint() ; // ctor |
| | 29 | AliEMCALRecPoint(const char * opt) ; // ctor |
| | 30 | AliEMCALRecPoint(const AliEMCALRecPoint & rp):AliRecPoint(rp) { Fatal("cpy ctor", "not implemented") ; } |
| | 31 | |
| | 32 | virtual ~AliEMCALRecPoint(); |
| | 33 | virtual void AddDigit(AliDigitNew &){ Fatal("AddDigit", "use AddDigit(AliEMCALDigit & digit, Float_t Energy )") ; } |
| | 34 | virtual void AddDigit(AliEMCALDigit & digit, Float_t Energy); |
| | 35 | virtual Int_t Compare(const TObject * obj) const; |
| | 36 | virtual Int_t DistancetoPrimitive(Int_t px, Int_t py); |
| | 37 | virtual void Draw(Option_t * option="") ; |
| | 38 | virtual void ExecuteEvent(Int_t event, Int_t, Int_t) ; |
| | 39 | |
| | 40 | virtual void EvalAll(Float_t logWeight, TClonesArray * digits); |
| | 41 | virtual void EvalLocalPosition(Float_t logWeight, TClonesArray * digits) ; |
| | 42 | virtual void EvalPrimaries(TClonesArray * digits) ; |
| | 43 | virtual void EvalParents(TClonesArray * digits) ; |
| | 44 | |
| | 45 | // virtual void GetGlobalPosition(TVector3 & gpos, TMatrix & /*gmat*/) const; // return global position in ALICE |
| | 46 | virtual void GetGlobalPosition(TVector3 & gpos) const; // return global position (x, y, z) in ALICE |
| | 47 | virtual void GetLocalPosition(TVector3 & lpos) const; // return local position (eta, phi, r) in EMCAL |
| | 48 | virtual Int_t * GetPrimaries(Int_t & number) const {number = fMulTrack ; |
| | 49 | return fTracksList ; } |
| | 50 | virtual Int_t * GetParents(Int_t & number) const {number = fMulParent ; |
| | 51 | return fParentsList ; } |
| | 52 | Float_t GetCoreEnergy()const {return fCoreEnergy ;} |
| | 53 | virtual Float_t GetDispersion()const {return fDispersion ;} |
| | 54 | virtual void GetElipsAxis(Float_t * lambda)const {lambda[0] = fLambda[0]; lambda[1] = fLambda[1];}; |
| | 55 | |
| | 56 | Float_t * GetEnergiesList() const {return fEnergyList ;} // gets the list of energies making this recpoint |
| | 57 | Float_t GetMaximalEnergy(void) const ; // get the highest energy in the cluster |
| | 58 | Int_t GetMaximumMultiplicity() const {return fMaxDigit ;} // gets the maximum number of digits allowed |
| | 59 | Int_t GetMultiplicity(void) const { return fMulDigit ; } // gets the number of digits making this recpoint |
| | 60 | Int_t GetMultiplicityAtLevel(Float_t level) const ; // computes multiplicity of digits with |
| | 61 | // energy above relative level |
| | 62 | virtual Int_t GetNumberOfLocalMax(AliEMCALDigit ** maxAt, Float_t * maxAtEnergy, |
| | 63 | Float_t locMaxCut,TClonesArray * digits ) const ; |
| | 64 | // searches for the local maxima |
| | 65 | Float_t GetTime(void) const{return fTime ; } |
| | 66 | |
| | 67 | virtual Bool_t IsEmc(void)const { return kTRUE ; } |
| | 68 | virtual Bool_t IsSortable() const { |
| | 69 | // tells that this is a sortable object |
| | 70 | return kTRUE ; |
| | 71 | } |
| | 72 | virtual void Paint(Option_t * option=""); |
| | 73 | virtual void Print(Option_t * option="") const ; |
| | 74 | |
| | 75 | AliEMCALRecPoint & operator = (const AliEMCALRecPoint & ) { |
| | 76 | Fatal("operator =", "not implemented") ; |
| | 77 | return *this ; |
| | 78 | } |
| | 79 | |
| | 80 | protected: |
| | 81 | void EvalCoreEnergy(Float_t logWeight,TClonesArray * digits) ; |
| | 82 | virtual void EvalDispersion(Float_t logWeight,TClonesArray * digits) ; // computes the dispersion of the shower |
| | 83 | virtual void EvalElipsAxis(Float_t logWeight, TClonesArray * digits ); // computes the axis of shower ellipsoide |
| | 84 | void EvalTime( TClonesArray * digits ); |
| | 85 | virtual Bool_t AreNeighbours(AliEMCALDigit * digit1, AliEMCALDigit * digit2 ) const; |
| | 86 | Float_t ThetaToEta(Float_t arg) const; //Converts Theta (Radians) to Eta(Radians) |
| | 87 | Float_t EtaToTheta(Float_t arg) const; //Converts Eta (Radians) to Theta(Radians) |
| | 88 | |
| | 89 | Float_t fCoreEnergy ; // energy in a shower core |
| | 90 | Float_t fLambda[2] ; // shower ellipse axes |
| | 91 | Float_t fDispersion ; // shower dispersion |
| | 92 | Float_t *fEnergyList ; //[fMulDigit] energy of digits |
| | 93 | Float_t fTime ; // Time of the digit with maximal energy deposition |
| | 94 | Float_t fCoreRadius; // The radius in which the core energy is evaluated |
| | 95 | Int_t fMulParent; // Multiplicity of the parents |
| | 96 | Int_t fMaxParent; // Maximum number of parents allowed |
| | 97 | Int_t * fParentsList; // [fMulParent] list of the parents of the digits |
| | 98 | |
| | 99 | ClassDef(AliEMCALRecPoint,6) // RecPoint for EMCAL (Base Class) |
| | 100 | |
| | 101 | }; |
| | 102 | |
| | 103 | #endif // AliEMCALRECPOINT_H |
git://git.uio.no / u / mrichter / AliRoot.git / blame_incremental