| d15a28e7 |
1 | #ifndef ALIPHOSEMCRECPOINT_H |
| 2 | #define ALIPHOSEMCRECPOINT_H |
| 3 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * |
| 4 | * See cxx source for full Copyright notice */ |
| 5 | |
| 6ad0bfa0 |
6 | /* $Id$ */ |
| 7 | |
| 702ab87e |
8 | /* History of cvs commits: |
| 9 | * |
| 10 | * $Log$ |
| c307c629 |
11 | * Revision 1.36 2007/04/05 10:18:58 policheh |
| 12 | * Introduced distance to nearest bad crystal. |
| 13 | * |
| dce76bbb |
14 | * Revision 1.35 2007/03/06 06:47:28 kharlov |
| 15 | * DP:Possibility to use actual vertex position added |
| 16 | * |
| 91daaf24 |
17 | * Revision 1.34 2005/05/28 14:19:04 schutz |
| 18 | * Compilation warnings fixed by T.P. |
| 19 | * |
| 702ab87e |
20 | */ |
| 21 | |
| b2a60966 |
22 | //_________________________________________________________________________ |
| 23 | // RecPoint implementation for PHOS-EMC |
| 24 | // An EmcRecPoint is a cluster of digits |
| 25 | // |
| 26 | //*-- Author: Dmitri Peressounko (RRC KI & SUBATECH) |
| d15a28e7 |
27 | |
| 28 | // --- ROOT system --- |
| 29 | |
| d15a28e7 |
30 | // --- Standard library --- |
| 31 | |
| 32 | // --- AliRoot header files --- |
| 33 | |
| e957fea8 |
34 | //#include "AliPHOSDigit.h" |
| d15a28e7 |
35 | #include "AliPHOSRecPoint.h" |
| 36 | |
| 37 | class AliPHOSEmcRecPoint : public AliPHOSRecPoint { |
| 38 | |
| 39 | public: |
| 40 | |
| 7932f811 |
41 | AliPHOSEmcRecPoint() ; |
| 73a68ccb |
42 | AliPHOSEmcRecPoint(const char * opt) ; |
| 55fe9d13 |
43 | AliPHOSEmcRecPoint(const AliPHOSEmcRecPoint & rp) ; |
| 6c370def |
44 | |
| 83974468 |
45 | virtual ~AliPHOSEmcRecPoint() ; |
| b2a60966 |
46 | |
| 702ab87e |
47 | //This virtual function has signature different from AliPHOSRecPoint::AddDigit |
| 48 | //it hides, not overrides. using - declaration should fix the problem, at least for |
| 49 | //g++ |
| 50 | using AliPHOSRecPoint::AddDigit; |
| 7932f811 |
51 | virtual void AddDigit(AliPHOSDigit & digit, Float_t Energy) ; // add a digit to the digits list |
| 2a941f4e |
52 | Int_t Compare(const TObject * obj) const; // method for sorting |
| 7932f811 |
53 | |
| 91daaf24 |
54 | virtual void EvalAll(Float_t logWeight, TClonesArray * digits) ; //Those tasks which can be done without vertex |
| 55 | virtual void EvalAll(Float_t logWeight, TVector3 &vtx, TClonesArray * digits) ; |
| 9688c1dd |
56 | |
| 702ab87e |
57 | //in base class this functions is non-const |
| 58 | virtual void ExecuteEvent(Int_t event, Int_t px, Int_t py) /*const*/; |
| 7932f811 |
59 | |
| 60 | Float_t GetCoreEnergy()const {return fCoreEnergy ;} |
| 61 | virtual Float_t GetDispersion()const {return fDispersion ;} |
| 62 | virtual void GetElipsAxis(Float_t * lambda)const { lambda[0] = fLambda[0] ; |
| 63 | lambda[1] = fLambda[1] ; } |
| 64 | Float_t * GetEnergiesList() const {return fEnergyList ;} // gets the list of energies making this recpoint |
| ad8cfaf4 |
65 | Float_t GetMaximalEnergy(void) const ; // get the highest energy in the cluster |
| 66 | Int_t GetMaximumMultiplicity() const {return fMaxDigit ;} // gets the maximum number of digits allowed |
| 67 | Int_t GetMultiplicity(void) const { return fMulDigit ; } // gets the number of digits making this recpoint |
| fc7e2f43 |
68 | Int_t GetMultiplicityAtLevel(Float_t level) const ; // computes multiplicity of digits with |
| ad8cfaf4 |
69 | // energy above relative level |
| 0bc3b8ed |
70 | Short_t GetNExMax(void) const {return fNExMax ;} // Number of maxima found in cluster in unfolding: |
| 71 | // 0: was no unfolging |
| 72 | //-1: unfolding failed |
| ce2a9a95 |
73 | void SetNExMax(Int_t nmax=1){fNExMax = static_cast<Short_t>(nmax) ;} |
| a0636361 |
74 | virtual Int_t GetNumberOfLocalMax(AliPHOSDigit ** maxAt, Float_t * maxAtEnergy, |
| 7932f811 |
75 | Float_t locMaxCut,TClonesArray * digits ) const ; |
| 76 | // searches for the local maxima |
| 092b50ba |
77 | //returns number of local maxima in parent cluster or -2 if unfolding failed |
| 9688c1dd |
78 | Float_t GetTime(void) const{return fTime ; } |
| ad8cfaf4 |
79 | Bool_t IsEmc(void) const { return kTRUE ; } // true if the recpoint is in EMC |
| 80 | Bool_t IsSortable() const {return kTRUE ; } // says that emcrecpoints are sortable objects |
| a8c47ab6 |
81 | void Print(Option_t *)const ; |
| 092b50ba |
82 | void Purify(Float_t threshold) ; //Removes digits below threshold |
| d15a28e7 |
83 | |
| ce2a9a95 |
84 | Float_t GetM2x() const {return fM2x; } // Get second X-moment |
| 85 | Float_t GetM2z() const {return fM2z; } // Get second Z-moment |
| 86 | Float_t GetM3x() const {return fM3x; } // Get third X-moment |
| 87 | Float_t GetM4z() const {return fM4z; } // Get forth Z-moment |
| 88 | Float_t GetPhixe() const {return fPhixe;} // Get angle between center gravity and eigen vector |
| 89 | |
| dce76bbb |
90 | Float_t GetDistanceToBadCrystal() const {return fDistToBadCrystal;} |
| 91 | void SetDistanceToBadCrystal(Float_t dist) {fDistToBadCrystal=dist;} |
| 92 | |
| 8c140292 |
93 | AliPHOSEmcRecPoint & operator = (const AliPHOSEmcRecPoint & /*rvalue*/) { return *this ; } |
| 6c370def |
94 | |
| 7932f811 |
95 | protected: |
| 9688c1dd |
96 | void EvalCoreEnergy(Float_t logWeight,TClonesArray * digits) ; |
| 91daaf24 |
97 | virtual void EvalLocalPosition(Float_t logWeight, TVector3 &vtx, TClonesArray * digits, TVector3 &vInc) ;// computes the position in the PHOS module |
| 98 | virtual void EvalDispersion(Float_t logWeight, TClonesArray * digits, TVector3 &vInc) ; // computes the dispersion of the shower |
| 99 | virtual void EvalElipsAxis(Float_t logWeight, TClonesArray * digits, TVector3 &vInc ); // computes the axis of shower ellipsoide |
| 100 | void EvalMoments(Float_t logWeight, TClonesArray * digits, TVector3 &vInc ); // computes shower moments |
| c307c629 |
101 | virtual void EvalPrimaries(TClonesArray * digits) ; |
| 9688c1dd |
102 | void EvalTime( TClonesArray * digits ); |
| 7932f811 |
103 | virtual Bool_t AreNeighbours(AliPHOSDigit * digit1, AliPHOSDigit * digit2 ) const ; |
| d15a28e7 |
104 | |
| baef0810 |
105 | Float_t fCoreEnergy ; // energy in a shower core |
| 106 | Float_t fLambda[2] ; // shower ellipse axes |
| 107 | Float_t fDispersion ; // shower dispersion |
| 108 | Float_t *fEnergyList ; //[fMulDigit] energy of digits |
| 9688c1dd |
109 | Float_t fTime ; // Time of the digit with maximal energy deposition |
| ce2a9a95 |
110 | Short_t fNExMax ; // number of (Ex-)maxima before unfolding |
| 111 | |
| 112 | Float_t fM2x; // Second moment along X axis |
| 113 | Float_t fM2z; // Second moment along Z axis |
| 114 | Float_t fM3x; // Third moment along X axis |
| 115 | Float_t fM4z; // Forth moment along Z axis |
| 116 | Float_t fPhixe; // Angle between center-gravity vector and eigen vector |
| dce76bbb |
117 | Float_t fDistToBadCrystal; // Distance to nearest bad crystal |
| ce2a9a95 |
118 | |
| 119 | Int_t fDebug; //! debug level (0 - no output) |
| 83974468 |
120 | |
| dce76bbb |
121 | ClassDef(AliPHOSEmcRecPoint,3) // EMC RecPoint (cluster) |
| d15a28e7 |
122 | |
| 123 | }; |
| 124 | |
| 125 | #endif // AliPHOSEMCRECPOINT_H |
git://git.uio.no / u / mrichter / AliRoot.git / blame