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 */
8 /* History of cvs commits:
10 * $Log: AliPHOSClusterizerv1.h,v $
11 * Revision 1.54 2007/10/01 20:24:08 kharlov
14 * Revision 1.53 2007/08/28 12:55:07 policheh
15 * Loaders removed from the reconstruction code (C.Cheshkov)
17 * Revision 1.52 2007/08/07 14:16:00 kharlov
18 * Quality assurance added (Yves Schutz)
20 * Revision 1.51 2007/04/11 11:55:45 policheh
21 * SetDistancesToBadChannels() added.
23 * Revision 1.50 2007/03/28 19:18:15 kharlov
24 * RecPoints recalculation in TSM removed
26 * Revision 1.49 2007/03/06 06:51:27 kharlov
27 * Calculation of cluster properties dep. on vertex posponed to TrackSegmentMaker
29 * Revision 1.48 2006/08/30 16:12:52 kharlov
30 * Reconstruction of raw data from beam test 2006 (B.Polichtchouk)
32 * Revision 1.47 2006/08/25 16:56:30 kharlov
33 * Compliance with Effective C++
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
38 * Revision 1.45 2006/04/29 20:26:46 hristov
39 * Separate EMC and CPV calibration (Yu.Kharlov)
41 * Revision 1.44 2005/09/02 14:32:07 kharlov
42 * Calibration of raw data
44 * Revision 1.43 2005/05/28 14:19:04 schutz
45 * Compilation warnings fixed by T.P.
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
56 //*-- Author: Yves Schutz (SUBATECH)
58 // --- ROOT system ---
61 // --- Standard library ---
63 // --- AliRoot header files ---
65 #include "AliPHOSClusterizer.h"
66 class AliPHOSEmcRecPoint ;
68 class AliPHOSDigitizer ;
69 class AliPHOSGeometry ;
71 class AliPHOSClusterizerv1 : public AliPHOSClusterizer {
75 AliPHOSClusterizerv1() ;
76 AliPHOSClusterizerv1(AliPHOSGeometry *geom);
77 virtual ~AliPHOSClusterizerv1() ;
79 void InitParameters() ;
81 virtual Int_t AreNeighbours(AliPHOSDigit * d1, AliPHOSDigit * d2)const ;
82 // Checks if digits are in neighbour cells
84 virtual void GetNumberOfClustersFound(int * numb )const{ numb[0] = fNumberOfEmcClusters ;
85 numb[1] = fNumberOfCpvClusters ; }
87 virtual Float_t GetEmcClusteringThreshold()const{ return fEmcClusteringThreshold;}
88 virtual Float_t GetEmcLocalMaxCut()const { return fEmcLocMaxCut;}
89 virtual Float_t GetEmcLogWeight()const { return fW0;}
90 virtual Float_t GetEmcTimeGate() const { return fEmcTimeGate ; }
91 virtual Float_t GetCpvClusteringThreshold()const{ return fCpvClusteringThreshold; }
92 virtual Float_t GetCpvLocalMaxCut()const { return fCpvLocMaxCut;}
93 virtual Float_t GetCpvLogWeight()const { return fW0CPV;}
94 virtual Float_t GetEcoreRadius()const { return fEcoreRadius;}
95 // virtual const char * GetRecPointsBranch() const{ return GetName() ;}
97 virtual void Digits2Clusters(Option_t *option);
99 void Print(const Option_t * = "")const ;
101 virtual void SetEmcClusteringThreshold(Float_t cluth) { fEmcClusteringThreshold = cluth ; }
102 virtual void SetEmcLocalMaxCut(Float_t cut) { fEmcLocMaxCut = cut ; }
103 virtual void SetEmcLogWeight(Float_t w) { fW0 = w ; }
104 virtual void SetEmcTimeGate(Float_t gate) { fEmcTimeGate = gate ;}
105 virtual void SetCpvClusteringThreshold(Float_t cluth) { fCpvClusteringThreshold = cluth ; }
106 virtual void SetCpvLocalMaxCut(Float_t cut) { fCpvLocMaxCut = cut ; }
107 virtual void SetCpvLogWeight(Float_t w) { fW0CPV = w ; }
108 virtual void SetUnfolding(Bool_t toUnfold = kTRUE ) { fToUnfold = toUnfold ;}
109 virtual void SetCoreRadius(Float_t coreRadius) { fEcoreRadius = coreRadius ;}
110 //Switch to "on flyght" mode, without writing to TreeR and file
111 void SetWriting(Bool_t toWrite = kFALSE){fWrite = toWrite;}
112 static Double_t ShowerShape(Double_t x, Double_t z) ; // Shape of EM shower used in unfolding;
113 //class member function (not object member function)
114 static void UnfoldingChiSquare(Int_t & nPar, Double_t * Grad, Double_t & fret, Double_t * x, Int_t iflag) ;
115 // Chi^2 of the fit. Should be static to be passed to MINUIT
117 virtual const char * Version() const { return "clu-v1"; }
121 void WriteRecPoints() ;
122 virtual void MakeClusters( ) ;
123 virtual Bool_t IsInEmc (AliPHOSDigit * digit)const ; // Tells if id digit is in EMC
124 virtual Bool_t IsInCpv (AliPHOSDigit * digit)const ; // Tells if id digit is in CPV
125 void CleanDigits(TClonesArray * digits) ;
126 void SetDistancesToBadChannels();
127 virtual Float_t Calibrate(Float_t amp, Int_t absId) ; // Tranforms ADC counts to energy
130 AliPHOSClusterizerv1(const AliPHOSClusterizerv1 & clu) ;
131 AliPHOSClusterizerv1 & operator = (const AliPHOSClusterizerv1 & obj);
133 Bool_t FindFit(AliPHOSEmcRecPoint * emcRP, AliPHOSDigit ** MaxAt, Float_t * maxAtEnergy,
134 Int_t NPar, Float_t * FitParametres) const; //Used in UnfoldClusters, calls TMinuit
137 virtual void MakeUnfolding() ;
138 void UnfoldCluster(AliPHOSEmcRecPoint * iniEmc,Int_t Nmax,
139 AliPHOSDigit ** maxAt,Float_t * maxAtEnergy ) ; //Unfolds cluster using TMinuit package
140 void PrintRecPoints(Option_t * option) ;
144 Bool_t fDefaultInit; //! Says if the task was created by defaut ctor (only parameters are initialized)
145 Int_t fEmcCrystals ; // number of EMC cristals in PHOS
147 Bool_t fToUnfold ; // To perform unfolding
148 Bool_t fWrite ; // Write RecPoints to TreeR
149 Bool_t fDigitsUsed[53760]; //Mark digits as already used in cluster (EMC:5*56*64 ; CPV: 5*56*128)
151 Int_t fNumberOfEmcClusters ; // number of EMC clusters found
152 Int_t fNumberOfCpvClusters ; // number of CPV clusters found
154 Float_t fEmcClusteringThreshold ; // minimum energy to start EMC cluster
155 Float_t fCpvClusteringThreshold ; // minimum energy to start CPV cluster
156 Float_t fEmcLocMaxCut ; // minimum energy difference to distinguish local maxima in a cluster
157 Float_t fW0 ; // logarithmic weight for the cluster center of gravity calculation
158 Float_t fCpvLocMaxCut ; // minimum energy difference to distinguish local maxima in a CPV cluster
159 Float_t fW0CPV ; // logarithmic weight for the CPV cluster center of gravity calculation
160 Float_t fEmcTimeGate ; // Maximum time difference between the digits in ont EMC cluster
161 Float_t fEcoreRadius ; // Radius within which the core energy is calculated, in cm
163 ClassDef(AliPHOSClusterizerv1,7) // Clusterizer implementation version 1
167 #endif // AliPHOSCLUSTERIZERV1_H