possiblity to read ESD friends and TPC outer check (Jacek)
[u/mrichter/AliRoot.git] / EMCAL / AliEMCALRecParam.h
CommitLineData
3a8be91c 1#ifndef ALIEMCALRECPARAM_H
2#define ALIEMCALRECPARAM_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//-----------------------------------------------------------------------------
9// Container of EMCAL reconstruction parameters
10// The purpose of this object is to store it to OCDB
b4133f05 11// and retrieve it in AliEMCALClusterizerv1, AliEMCALPID,
12// AliEMCALTracker and use it to configure AliEMCALRawUtils
13//
14//
3a8be91c 15// Author: Yuri Kharlov
16//-----------------------------------------------------------------------------
17
18// --- ROOT system ---
19
413e6b81 20#include "AliDetectorRecoParam.h"
225cd96d 21#include "AliLog.h"
3a8be91c 22
413e6b81 23class AliEMCALRecParam : public AliDetectorRecoParam
3a8be91c 24{
25public:
26
27 AliEMCALRecParam() ;
413e6b81 28 AliEMCALRecParam(const AliEMCALRecParam& recParam);
29 AliEMCALRecParam& operator = (const AliEMCALRecParam& recParam);
3a8be91c 30 virtual ~AliEMCALRecParam() {}
8ba062b1 31
225cd96d 32 //Clustering (Unfolding : Cynthia)
c47157cd 33 Float_t GetClusteringThreshold() const {return fClusteringThreshold;}
34 Float_t GetW0 () const {return fW0 ;}
35 Float_t GetMinECut () const {return fMinECut ;}
225cd96d 36 Float_t GetLocMaxCut () const {return fLocMaxCut ;}
37 Bool_t GetUnfold () const {return fUnfold ;}
3a8be91c 38 void SetClusteringThreshold(Float_t thrsh) {fClusteringThreshold = thrsh;}
39 void SetW0 (Float_t w0) {fW0 = w0 ;}
40 void SetMinECut (Float_t minEcut) {fMinECut = minEcut ;}
225cd96d 41 void SetLocMaxCut (Float_t locMaxCut) {fLocMaxCut = locMaxCut ;}
0561e246 42 void SetUnfold (Bool_t unfold) {fUnfold = unfold ; if(fUnfold) AliWarning("Cluster Unfolding ON. Implementing only for eta=0 case!!!");}
8ba062b1 43
44 //PID (Guenole)
45 Double_t GetGamma(Int_t i, Int_t j) const {return fGamma[i][j];}
46 Double_t GetHadron(Int_t i, Int_t j) const {return fHadron[i][j];}
47 Double_t GetPiZero5to10(Int_t i, Int_t j) const {return fPiZero5to10[i][j];}
48 Double_t GetPiZero10to60(Int_t i, Int_t j) const {return fPiZero10to60[i][j];}
49
50 void SetGamma(Int_t i, Int_t j,Double_t param ) {fGamma[i][j]=param;}
51 void SetHadron(Int_t i, Int_t j,Double_t param ) {fHadron[i][j]=param;}
52 void SetPiZero5to10(Int_t i, Int_t j,Double_t param) {fPiZero5to10[i][j]=param;}
53 void SetPiZero10to60(Int_t i, Int_t j,Double_t param) {fPiZero10to60[i][j]=param;}
54
55 //Track Matching (Alberto)
56 /* track matching cut setters */
57 void SetTrkCutX(Double_t value) {fTrkCutX = value;}
58 void SetTrkCutY(Double_t value) {fTrkCutY = value;}
59 void SetTrkCutZ(Double_t value) {fTrkCutZ = value;}
60 void SetTrkCutR(Double_t value) {fTrkCutR = value;}
61 void SetTrkCutAlphaMin(Double_t value) {fTrkCutAlphaMin = value;}
62 void SetTrkCutAlphaMax(Double_t value) {fTrkCutAlphaMax = value;}
63 void SetTrkCutAngle(Double_t value) {fTrkCutAngle = value;}
dcd86c5d 64 void SetTrkCutNITS(Double_t value) {fTrkCutNITS = value;}
65 void SetTrkCutNTPC(Double_t value) {fTrkCutNTPC = value;}
8ba062b1 66 /* track matching cut getters */
67 Double_t GetTrkCutX() const {return fTrkCutX;}
68 Double_t GetTrkCutY() const {return fTrkCutY;}
69 Double_t GetTrkCutZ() const {return fTrkCutZ;}
70 Double_t GetTrkCutR() const {return fTrkCutR;}
71 Double_t GetTrkCutAlphaMin() const {return fTrkCutAlphaMin;}
72 Double_t GetTrkCutAlphaMax() const {return fTrkCutAlphaMax;}
73 Double_t GetTrkCutAngle() const {return fTrkCutAngle;}
dcd86c5d 74 Double_t GetTrkCutNITS() const {return fTrkCutNITS;}
75 Double_t GetTrkCutNTPC() const {return fTrkCutNTPC;}
8ba062b1 76
b4133f05 77 //Raw signal fitting (Jenn)
78 /* raw signal setters */
79 void SetHighLowGainFactor(Double_t value) {fHighLowGainFactor = value;}
80 void SetOrderParameter(Int_t value) {fOrderParameter = value;}
81 void SetTau(Double_t value) {fTau = value;}
82 void SetNoiseThreshold(Int_t value) {fNoiseThreshold = value;}
83 void SetNPedSamples(Int_t value) {fNPedSamples = value;}
84 /* raw signal getters */
85 Double_t GetHighLowGainFactor() const {return fHighLowGainFactor;}
86 Int_t GetOrderParameter() const {return fOrderParameter;}
87 Double_t GetTau() const {return fTau;}
88 Int_t GetNoiseThreshold() const {return fNoiseThreshold;}
89 Int_t GetNPedSamples() const {return fNPedSamples;}
90
8ba062b1 91 virtual void Print(Option_t * option="") const ;
92
413e6b81 93 static AliEMCALRecParam* GetDefaultParameters();
f5fc991a 94 static AliEMCALRecParam* GetLowFluxParam();
95 static AliEMCALRecParam* GetHighFluxParam();
96
feedcab9 97 static const TObjArray* GetMappings();
3a8be91c 98
99private:
8ba062b1 100 //Clustering
3a8be91c 101 Float_t fClusteringThreshold ; // minimum energy to seed a EC digit in a cluster
102 Float_t fW0 ; // logarithmic weight for the cluster center of gravity calculation
103 Float_t fMinECut; // Minimum energy for a digit to be a member of a cluster
225cd96d 104 Bool_t fUnfold; // flag to perform cluster unfolding
105 Float_t fLocMaxCut; // minimum energy difference to consider local maxima in a cluster
3a8be91c 106
8ba062b1 107 //PID (Guenole)
108 Double_t fGamma[6][6]; // Parameter to Compute PID
109 Double_t fHadron[6][6]; // Parameter to Compute PID
110 Double_t fPiZero5to10[6][6]; // Parameter to Compute PID
111 Double_t fPiZero10to60[6][6]; // Parameter to Compute PID
112
113 //Track-Matching (Alberto)
114 Double_t fTrkCutX; // X-difference cut for track matching
115 Double_t fTrkCutY; // Y-difference cut for track matching
116 Double_t fTrkCutZ; // Z-difference cut for track matching
117 Double_t fTrkCutR; // cut on allowed track-cluster distance
118 Double_t fTrkCutAlphaMin; // cut on 'alpha' parameter for track matching (min)
119 Double_t fTrkCutAlphaMax; // cut on 'alpha' parameter for track matching (min)
120 Double_t fTrkCutAngle; // cut on relative angle between different track points for track matching
dcd86c5d 121 Double_t fTrkCutNITS; // Number of ITS hits for track matching
122 Double_t fTrkCutNTPC; // Number of TPC hits for track matching
123
b4133f05 124 //Raw signal fitting parameters (Jenn)
125 Double_t fHighLowGainFactor; //gain factor to convert between high and low gain
126 Int_t fOrderParameter; //order parameter for raw signal fit
127 Double_t fTau; //decay constant for raw signal fit
128 Int_t fNoiseThreshold; //threshold to consider signal or noise
129 Int_t fNPedSamples; //number of time samples to use in pedestal calculation
130
feedcab9 131 static TObjArray* fgkMaps; // ALTRO mappings for RCU0..RCUX
132
413e6b81 133 ClassDef(AliEMCALRecParam,6) // Reconstruction parameters
3a8be91c 134
135} ;
136
137#endif // ALIEMCALRECPARAM_H
8ba062b1 138