Initialization of TVector3 data member is corrected
[u/mrichter/AliRoot.git] / PHOS / AliPHOSClusterizerv1.h
CommitLineData
e9a88722 1#ifndef ALIPHOSCLUSTERIZERV1_H
d15a28e7 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
6ad0bfa0 6/* $Id$ */
7
702ab87e 8/* History of cvs commits:
9 *
a28333c5 10 * $Log: AliPHOSClusterizerv1.h,v $
11 * Revision 1.54 2007/10/01 20:24:08 kharlov
12 * Memory leaks fixed
13 *
e68222ce 14 * Revision 1.53 2007/08/28 12:55:07 policheh
15 * Loaders removed from the reconstruction code (C.Cheshkov)
16 *
9a2cdbdf 17 * Revision 1.52 2007/08/07 14:16:00 kharlov
18 * Quality assurance added (Yves Schutz)
19 *
ddd1a39c 20 * Revision 1.51 2007/04/11 11:55:45 policheh
21 * SetDistancesToBadChannels() added.
22 *
af0b0570 23 * Revision 1.50 2007/03/28 19:18:15 kharlov
24 * RecPoints recalculation in TSM removed
25 *
8c1fb709 26 * Revision 1.49 2007/03/06 06:51:27 kharlov
27 * Calculation of cluster properties dep. on vertex posponed to TrackSegmentMaker
28 *
a978bcea 29 * Revision 1.48 2006/08/30 16:12:52 kharlov
30 * Reconstruction of raw data from beam test 2006 (B.Polichtchouk)
31 *
02e48b45 32 * Revision 1.47 2006/08/25 16:56:30 kharlov
33 * Compliance with Effective C++
34 *
0378398c 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 *
f5eaa851 38 * Revision 1.45 2006/04/29 20:26:46 hristov
39 * Separate EMC and CPV calibration (Yu.Kharlov)
40 *
e95226ae 41 * Revision 1.44 2005/09/02 14:32:07 kharlov
42 * Calibration of raw data
43 *
44ae287e 44 * Revision 1.43 2005/05/28 14:19:04 schutz
45 * Compilation warnings fixed by T.P.
46 *
702ab87e 47 */
48
b2a60966 49//_________________________________________________________________________
50// Implementation version 1 of the clusterization algorithm
baef0810 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
b2a60966 55//
88cb7938 56//*-- Author: Yves Schutz (SUBATECH)
d15a28e7 57
58// --- ROOT system ---
ba898748 59class TClonesArray ;
a978bcea 60class TVector3 ;
d15a28e7 61// --- Standard library ---
62
63// --- AliRoot header files ---
64
65#include "AliPHOSClusterizer.h"
9a1398dd 66class AliPHOSEmcRecPoint ;
67class AliPHOSDigit ;
68class AliPHOSDigitizer ;
69class AliPHOSGeometry ;
d15a28e7 70
71class AliPHOSClusterizerv1 : public AliPHOSClusterizer {
72
73public:
74
0378398c 75 AliPHOSClusterizerv1() ;
9a2cdbdf 76 AliPHOSClusterizerv1(AliPHOSGeometry *geom);
9688c1dd 77 virtual ~AliPHOSClusterizerv1() ;
7e88424f 78
79 void InitParameters() ;
d15a28e7 80
1f96ab2a 81 virtual Int_t AreNeighbours(AliPHOSDigit * d1, AliPHOSDigit * d2)const ;
9a1398dd 82 // Checks if digits are in neighbour cells
3758d9fc 83
baef0810 84 virtual void GetNumberOfClustersFound(int * numb )const{ numb[0] = fNumberOfEmcClusters ;
85 numb[1] = fNumberOfCpvClusters ; }
9a1398dd 86
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 GetCpvClusteringThreshold()const{ return fCpvClusteringThreshold; }
91 virtual Float_t GetCpvLocalMaxCut()const { return fCpvLocMaxCut;}
92 virtual Float_t GetCpvLogWeight()const { return fW0CPV;}
d55c9739 93 virtual Float_t GetEcoreRadius()const { return fEcoreRadius;}
9a2cdbdf 94 // virtual const char * GetRecPointsBranch() const{ return GetName() ;}
9a1398dd 95
9a2cdbdf 96 virtual void Digits2Clusters(Option_t *option);
9a1398dd 97
702ab87e 98 void Print(const Option_t * = "")const ;
9a1398dd 99
d72dfbc3 100 virtual void SetEmcClusteringThreshold(Float_t cluth) { fEmcClusteringThreshold = cluth ; }
9a1398dd 101 virtual void SetEmcLocalMaxCut(Float_t cut) { fEmcLocMaxCut = cut ; }
102 virtual void SetEmcLogWeight(Float_t w) { fW0 = w ; }
ed4205d8 103 virtual void SetCpvClusteringThreshold(Float_t cluth) { fCpvClusteringThreshold = cluth ; }
9a1398dd 104 virtual void SetCpvLocalMaxCut(Float_t cut) { fCpvLocMaxCut = cut ; }
105 virtual void SetCpvLogWeight(Float_t w) { fW0CPV = w ; }
ba898748 106 virtual void SetUnfolding(Bool_t toUnfold = kTRUE ) { fToUnfold = toUnfold ;}
d55c9739 107 virtual void SetCoreRadius(Float_t coreRadius) { fEcoreRadius = coreRadius ;}
ba898748 108 //Switch to "on flyght" mode, without writing to TreeR and file
5d0435dd 109 void SetWriting(Bool_t toWrite = kFALSE){fWrite = toWrite;}
a978bcea 110 static Double_t ShowerShape(Double_t x, Double_t z) ; // Shape of EM shower used in unfolding;
9a1398dd 111 //class member function (not object member function)
7b7c1533 112 static void UnfoldingChiSquare(Int_t & nPar, Double_t * Grad, Double_t & fret, Double_t * x, Int_t iflag) ;
ba898748 113 // Chi^2 of the fit. Should be static to be passed to MINUIT
9a2cdbdf 114 // void Unload() ;
88cb7938 115 virtual const char * Version() const { return "clu-v1"; }
9a1398dd 116
1f96ab2a 117protected:
118
88cb7938 119 void WriteRecPoints() ;
1f96ab2a 120 virtual void MakeClusters( ) ;
121 virtual Bool_t IsInEmc (AliPHOSDigit * digit)const ; // Tells if id digit is in EMC
122 virtual Bool_t IsInCpv (AliPHOSDigit * digit)const ; // Tells if id digit is in CPV
ba898748 123 void CleanDigits(TClonesArray * digits) ;
af0b0570 124 void SetDistancesToBadChannels();
6f47f50d 125 virtual Float_t Calibrate(Float_t amp, Int_t absId) const ; // Tranforms ADC counts to energy
126 virtual Float_t CalibrateT(Float_t amp, Int_t absId) const ; //Tranforms Sample counts to sec.
67637327 127 Bool_t CheckTimeGate(Float_t t1, Float_t amp1, Float_t t2, Float_t amp2)const ; //Checks if time difference is reasonable
5024fc61 128
d15a28e7 129private:
9a2cdbdf 130 AliPHOSClusterizerv1(const AliPHOSClusterizerv1 & clu) ;
131 AliPHOSClusterizerv1 & operator = (const AliPHOSClusterizerv1 & obj);
7b7c1533 132
a0636361 133 Bool_t FindFit(AliPHOSEmcRecPoint * emcRP, AliPHOSDigit ** MaxAt, Float_t * maxAtEnergy,
baef0810 134 Int_t NPar, Float_t * FitParametres) const; //Used in UnfoldClusters, calls TMinuit
fbf811ec 135 void Init() ;
9a1398dd 136
9a1398dd 137 virtual void MakeUnfolding() ;
9a1398dd 138 void UnfoldCluster(AliPHOSEmcRecPoint * iniEmc,Int_t Nmax,
a0636361 139 AliPHOSDigit ** maxAt,Float_t * maxAtEnergy ) ; //Unfolds cluster using TMinuit package
9a1398dd 140 void PrintRecPoints(Option_t * option) ;
141
142private:
143
92f521a9 144 Bool_t fDefaultInit; //! Says if the task was created by defaut ctor (only parameters are initialized)
02e48b45 145 Int_t fEmcCrystals ; // number of EMC cristals in PHOS
3758d9fc 146
9a1398dd 147 Bool_t fToUnfold ; // To perform unfolding
ba898748 148 Bool_t fWrite ; // Write RecPoints to TreeR
d07a44a8 149 Bool_t fDigitsUsed[53760]; //Mark digits as already used in cluster (EMC:5*56*64 ; CPV: 5*56*128)
02e48b45 150
151 Int_t fNumberOfEmcClusters ; // number of EMC clusters found
9688c1dd 152 Int_t fNumberOfCpvClusters ; // number of CPV clusters found
7b7c1533 153
ba898748 154 Float_t fEmcClusteringThreshold ; // minimum energy to start EMC cluster
155 Float_t fCpvClusteringThreshold ; // minimum energy to start CPV cluster
9a1398dd 156 Float_t fEmcLocMaxCut ; // minimum energy difference to distinguish local maxima in a cluster
b2a60966 157 Float_t fW0 ; // logarithmic weight for the cluster center of gravity calculation
9a1398dd 158 Float_t fCpvLocMaxCut ; // minimum energy difference to distinguish local maxima in a CPV cluster
d72dfbc3 159 Float_t fW0CPV ; // logarithmic weight for the CPV cluster center of gravity calculation
67637327 160 Float_t fTimeGateLowAmp ; // Threshold for good/bad time measurement
161 Float_t fTimeGateLow ; // Time difference between cells with low amplitude
162 Float_t fTimeGateHigh ; // Time difference between cells with good time measurement
163
d55c9739 164 Float_t fEcoreRadius ; // Radius within which the core energy is calculated, in cm
f5eaa851 165
67637327 166 ClassDef(AliPHOSClusterizerv1,8) // Clusterizer implementation version 1
d15a28e7 167
168};
169
170#endif // AliPHOSCLUSTERIZERV1_H