TClonesArray->Clone() removed
[u/mrichter/AliRoot.git] / PHOS / AliPHOSClusterizerv1.h
1 #ifndef ALIPHOSCLUSTERIZERV1_H
2 #define ALIPHOSCLUSTERIZERV1_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4  * See cxx source for full Copyright notice                               */
5
6 /* $Id$ */
7
8 /* History of cvs commits:
9  *
10  * $Log: AliPHOSClusterizerv1.h,v $
11  * Revision 1.54  2007/10/01 20:24:08  kharlov
12  * Memory leaks fixed
13  *
14  * Revision 1.53  2007/08/28 12:55:07  policheh
15  * Loaders removed from the reconstruction code (C.Cheshkov)
16  *
17  * Revision 1.52  2007/08/07 14:16:00  kharlov
18  * Quality assurance added (Yves Schutz)
19  *
20  * Revision 1.51  2007/04/11 11:55:45  policheh
21  * SetDistancesToBadChannels() added.
22  *
23  * Revision 1.50  2007/03/28 19:18:15  kharlov
24  * RecPoints recalculation in TSM removed
25  *
26  * Revision 1.49  2007/03/06 06:51:27  kharlov
27  * Calculation of cluster properties dep. on vertex posponed to TrackSegmentMaker
28  *
29  * Revision 1.48  2006/08/30 16:12:52  kharlov
30  * Reconstruction of raw data from beam test 2006 (B.Polichtchouk)
31  *
32  * Revision 1.47  2006/08/25 16:56:30  kharlov
33  * Compliance with Effective C++
34  *
35  * Revision 1.46  2006/08/01 12:20:17  cvetan
36  * 1. Adding a possibility to read and reconstruct an old rcu formatted raw data. This is controlled by an option of AliReconstruction and AliPHOSReconstructor. 2. In case of raw data processing (without galice.root) create the default AliPHOSGeometry object. Most likely this should be moved to the CDB
37  *
38  * Revision 1.45  2006/04/29 20:26:46  hristov
39  * Separate EMC and CPV calibration (Yu.Kharlov)
40  *
41  * Revision 1.44  2005/09/02 14:32:07  kharlov
42  * Calibration of raw data
43  *
44  * Revision 1.43  2005/05/28 14:19:04  schutz
45  * Compilation warnings fixed by T.P.
46  *
47  */
48
49 //_________________________________________________________________________
50 //  Implementation version 1 of the clusterization algorithm                     
51 //  Performs clusterization (collects neighbouring active cells) and 
52 //  unfolding of the clusters with several local maxima.  
53 //  results are stored in TreeR#, branches PHOSEmcRP (EMC recPoints),
54 //  PHOSCpvRP (CPV RecPoints) and AliPHOSClusterizer
55 //
56 //*-- Author: Yves Schutz (SUBATECH)
57
58 // --- ROOT system ---
59 class TClonesArray ;
60 class TVector3 ;
61 // --- Standard library ---
62
63 // --- AliRoot header files ---
64
65 #include "AliPHOSClusterizer.h"
66 class AliPHOSEmcRecPoint ; 
67 class AliPHOSDigit ;
68 class AliPHOSDigitizer ;
69 class AliPHOSGeometry ;
70
71 class AliPHOSClusterizerv1 : public AliPHOSClusterizer {
72   
73 public:
74   
75   AliPHOSClusterizerv1() ;
76   AliPHOSClusterizerv1(AliPHOSGeometry *geom);
77   virtual ~AliPHOSClusterizerv1()  ;
78   
79   virtual Int_t   AreNeighbours(AliPHOSDigit * d1, AliPHOSDigit * d2)const ; 
80                                // Checks if digits are in neighbour cells 
81
82   virtual void    GetNumberOfClustersFound(int * numb )const{  numb[0] = fNumberOfEmcClusters ; 
83                                                                numb[1] = fNumberOfCpvClusters ; }
84
85   virtual Float_t GetEmcClusteringThreshold()const{ return fEmcClusteringThreshold;}
86   virtual Float_t GetEmcLocalMaxCut()const        { return fEmcLocMaxCut;} 
87   virtual Float_t GetEmcLogWeight()const          { return fW0;}  
88   virtual Float_t GetEmcTimeGate() const          { return fEmcTimeGate ; }
89   virtual Float_t GetCpvClusteringThreshold()const{ return fCpvClusteringThreshold;  } 
90   virtual Float_t GetCpvLocalMaxCut()const        { return fCpvLocMaxCut;} 
91   virtual Float_t GetCpvLogWeight()const          { return fW0CPV;}  
92   //  virtual const char *  GetRecPointsBranch() const{ return GetName() ;}
93
94   virtual void    Digits2Clusters(Option_t *option);
95
96   void Print(const Option_t * = "")const ;
97
98   virtual void SetEmcClusteringThreshold(Float_t cluth)  { fEmcClusteringThreshold = cluth ; }
99   virtual void SetEmcLocalMaxCut(Float_t cut)            { fEmcLocMaxCut = cut ; }
100   virtual void SetEmcLogWeight(Float_t w)                { fW0 = w ; }
101   virtual void SetEmcTimeGate(Float_t gate)              { fEmcTimeGate = gate ;}
102   virtual void SetCpvClusteringThreshold(Float_t cluth)  { fCpvClusteringThreshold = cluth ; }
103   virtual void SetCpvLocalMaxCut(Float_t cut)            { fCpvLocMaxCut = cut ; }
104   virtual void SetCpvLogWeight(Float_t w)                { fW0CPV = w ; }
105   virtual void SetUnfolding(Bool_t toUnfold = kTRUE )    { fToUnfold = toUnfold ;}
106   //Switch to "on flyght" mode, without writing to TreeR and file  
107   void SetWriting(Bool_t toWrite = kFALSE){fWrite = toWrite;} 
108   static Double_t ShowerShape(Double_t x, Double_t z) ; // Shape of EM shower used in unfolding; 
109                                             //class member function (not object member function)
110   static void UnfoldingChiSquare(Int_t & nPar, Double_t * Grad, Double_t & fret, Double_t * x, Int_t iflag)  ;
111                                             // Chi^2 of the fit. Should be static to be passed to MINUIT
112   //  void Unload() ; 
113   virtual const char * Version() const { return "clu-v1"; }  
114
115 protected:
116
117   void           WriteRecPoints() ;
118   virtual void   MakeClusters( ) ;            
119   virtual Bool_t IsInEmc (AliPHOSDigit * digit)const ;     // Tells if id digit is in EMC
120   virtual Bool_t IsInCpv (AliPHOSDigit * digit)const ;     // Tells if id digit is in CPV
121   void           CleanDigits(TClonesArray * digits) ;
122   void           SetDistancesToBadChannels();
123
124 private:
125   AliPHOSClusterizerv1(const AliPHOSClusterizerv1 & clu) ;
126   AliPHOSClusterizerv1 & operator = (const AliPHOSClusterizerv1 & obj);
127
128   Bool_t  FindFit(AliPHOSEmcRecPoint * emcRP, AliPHOSDigit ** MaxAt, Float_t * maxAtEnergy, 
129                   Int_t NPar, Float_t * FitParametres) const; //Used in UnfoldClusters, calls TMinuit
130   void    Init() ;
131   void    InitParameters() ;
132
133   virtual void   MakeUnfolding() ;
134   void           UnfoldCluster(AliPHOSEmcRecPoint * iniEmc,Int_t Nmax, 
135                        AliPHOSDigit ** maxAt,Float_t * maxAtEnergy ) ; //Unfolds cluster using TMinuit package
136   void           PrintRecPoints(Option_t * option) ;
137
138 private:
139
140   Bool_t  fDefaultInit;              //! Says if the task was created by defaut ctor (only parameters are initialized)
141   Int_t   fEmcCrystals ;             // number of EMC cristals in PHOS
142
143   Bool_t  fToUnfold ;                // To perform unfolding 
144   Bool_t  fWrite ;                   // Write RecPoints to TreeR  
145   Bool_t  fDigitsUsed[53760];        //Mark digits as already used in cluster (EMC:5*56*64 ; CPV: 5*56*128)
146  
147   Int_t   fNumberOfEmcClusters ;     // number of EMC clusters found
148   Int_t   fNumberOfCpvClusters ;     // number of CPV clusters found
149  
150   Float_t fEmcClusteringThreshold ;  // minimum energy to start EMC cluster
151   Float_t fCpvClusteringThreshold ;  // minimum energy to start CPV cluster
152   Float_t fEmcLocMaxCut ;            // minimum energy difference to distinguish local maxima in a cluster
153   Float_t fW0 ;                      // logarithmic weight for the cluster center of gravity calculation
154   Float_t fCpvLocMaxCut ;            // minimum energy difference to distinguish local maxima in a CPV cluster
155   Float_t fW0CPV ;                   // logarithmic weight for the CPV cluster center of gravity calculation
156   //  Int_t fRecPointsInRun ;            //! Total number of recpoints in one run
157   Float_t fEmcTimeGate ;             // Maximum time difference between the digits in ont EMC cluster
158
159
160   ClassDef(AliPHOSClusterizerv1,6)   // Clusterizer implementation version 1
161
162 };
163
164 #endif // AliPHOSCLUSTERIZERV1_H