implemented copy constructor (protected)
[u/mrichter/AliRoot.git] / EMCAL / AliEMCALRecParam.cxx
CommitLineData
3a8be91c 1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
16/* $Id$ */
17// --- AliRoot header files ---
18#include "AliEMCALRecParam.h"
19#include "AliLog.h"
20
21ClassImp(AliEMCALRecParam)
22
23//-----------------------------------------------------------------------------
24// Container of EMCAL reconstruction parameters
25// The purpose of this object is to store it to OCDB
26// and retrieve it in AliEMCALClusterizerv1
27// Author: Yuri Kharlov
28//-----------------------------------------------------------------------------
29
30AliEMCALRecParam::AliEMCALRecParam():
8ba062b1 31 fClusteringThreshold(0.5),fW0(4.5),fMinECut(0.45), //clustering
32 fTrkCutX(6.0), fTrkCutY(6.0), fTrkCutZ(6.0), fTrkCutR(10.0),//track matching
33 fTrkCutAlphaMin(-50.0), fTrkCutAlphaMax(50.0), fTrkCutAngle(10000.0) //track matching
c47157cd 34{
35 // default reco values
8ba062b1 36
37 //PID parameters (Guenole)
38
39 // as a first step, all array elements are initialized to 0.0
40 Int_t i, j;
41 for (i = 0; i < 6; i++) {
42 for (j = 0; j < 6; j++) {
43 fGamma[i][j] = fHadron[i][j] = fPiZero5to10[i][j] = fPiZero10to60[i][j] = 0.;
44 }
45 }
46
47 // then, only the ones which must be not zero are initialized
48 // while the others will remain to the value 0.0
49
50 fGamma[0][0] = 0.038022;
51 fGamma[0][1] = -0.0001883;
52 fGamma[0][2] = 5.449e-06;
53
54 fGamma[1][0] = 0.207313;
55 fGamma[1][1] = -0.000978;
56 fGamma[1][2] = 0.00001634;
57
58 fGamma[2][0] = 0.043364;
59 fGamma[2][1] = -0.0002048;
60 fGamma[2][2] = 8.661e-06;
61 fGamma[2][3] = -1.353e-07;
62
63 fGamma[3][0] = 0.265004;
64 fGamma[3][1] = 0.061298;
65 fGamma[3][2] = -0.003203;
66 fGamma[3][3] = 4.73e-05;
67
68 fGamma[4][0] = 0.243579;
69 fGamma[4][1] = -1.614e-05;
70
71 fGamma[5][0] = 0.002942;
72 fGamma[5][1] = -3.976e-05;
73
74 fHadron[0][0] = 0.011945 / 3.;
75 fHadron[0][1] = 0.000386 / 3.;
76 fHadron[0][2] = -0.000014 / 3.;
77 fHadron[0][3] = 1.336e-07 / 3.;
78
79 fHadron[1][0] = 0.496544;
80 fHadron[1][1] = -0.003226;
81 fHadron[1][2] = 0.00001678;
82
83 fHadron[2][0] = 0.144838;
84 fHadron[2][1] = -0.002954;
85 fHadron[2][2] = 0.00008754;
86 fHadron[2][3] = -7.587e-07;
87
88 fHadron[3][0] = 1.264461 / 7.;
89 fHadron[3][1] = 0.002097 / 7.;
90
91 fHadron[4][0] = 0.261950;
92 fHadron[4][1] = -0.001078;
93 fHadron[4][2] = 0.00003237;
94 fHadron[4][3] = -3.241e-07;
95 fHadron[4][4] = 0.;
96 fHadron[4][5] = 0.;
97 fHadron[5][0] = 0.010317;
98 fHadron[5][1] = 0.;
99 fHadron[5][2] = 0.;
100 fHadron[5][3] = 0.;
101 fHadron[5][4] = 0.;
102 fHadron[5][5] = 0.;
103
104 fPiZero5to10[0][0] = 0.009138;
105 fPiZero5to10[0][1] = 0.0006377;
106
107 fPiZero5to10[1][0] = 0.08;
108
109 fPiZero5to10[2][0] = -0.061119;
110 fPiZero5to10[2][1] = 0.019013;
111
112 fPiZero5to10[3][0] = 0.2;
113
114 fPiZero5to10[4][0] = 0.252044;
115 fPiZero5to10[4][1] = -0.002315;
116
117 fPiZero5to10[5][0] = 0.002942;
118 fPiZero5to10[5][1] = -3.976e-05;
119
120 fPiZero10to60[0][0] = 0.009138;
121 fPiZero10to60[0][1] = 0.0006377;
122
123 fPiZero10to60[1][0] = 1.272837;
124 fPiZero10to60[1][1] = -0.069708;
125 fPiZero10to60[1][2] = 0.001568;
126 fPiZero10to60[1][3] = -1.162e-05;
127
128 fPiZero10to60[2][0] = 0.139703;
129 fPiZero10to60[2][1] = 0.003687;
130 fPiZero10to60[2][2] = -0.000568;
131 fPiZero10to60[2][3] = 1.498e-05;
132 fPiZero10to60[2][4] = -1.174e-07;
133
134 fPiZero10to60[3][0] = -0.826367;
135 fPiZero10to60[3][1] = 0.096951;
136 fPiZero10to60[3][2] = -0.002215;
137 fPiZero10to60[3][3] = 2.523e-05;
138
139 fPiZero10to60[4][0] = 0.249890;
140 fPiZero10to60[4][1] = -0.000063;
141
142 fPiZero10to60[5][0] = 0.002942;
143 fPiZero10to60[5][1] = -3.976e-05;
144
145 //track matching
146
147
148
c47157cd 149}
3a8be91c 150
151//-----------------------------------------------------------------------------
152void AliEMCALRecParam::Print(Option_t *) const
153{
154 printf("AliEMCALRecParam::Print()\n");
155 // Print reconstruction parameters to stdout
8ba062b1 156 AliInfo(Form("Clusterization parameters :\n fClusteringThreshold=%.3f,\n fW0=%.3f,\n fMinECut=%.3f\n",
3a8be91c 157 fClusteringThreshold,fW0,fMinECut));
158
8ba062b1 159 AliInfo(Form("Track-matching cuts :\n x %f, y %f, z %f, R %f \n alphaMin %f, alphaMax %f, Angle %f\n", fTrkCutX, fTrkCutY, fTrkCutZ, fTrkCutR,fTrkCutAlphaMin,fTrkCutAlphaMax, fTrkCutAngle));
160
161 AliInfo(Form("PID parameters, Gamma :\n"));
162 for(Int_t i = 0; i < 6; i++){
163 for(Int_t j = 0; j < 6; j++){
164 printf(" %f, ", fGamma[i][j]);
165 }
166 printf("\n");
167 }
168
169 printf("\n");
170
171 AliInfo(Form("PID parameters, Hadron :\n"));
172 for(Int_t i = 0; i < 6; i++){
173 for(Int_t j = 0; j < 6; j++){
174 printf(" %f, ", fHadron[i][j]);
175 }
176 printf("\n");
177 }
178
179 printf("\n");
180
181 AliInfo(Form("PID parameters, Pi0zero5to10 :\n"));
182 for(Int_t i = 0; i < 6; i++){
183 for(Int_t j = 0; j < 6; j++){
184 printf(" %f, ", fPiZero5to10[i][j]);
185 }
186 printf("\n");
187 }
188
189 printf("\n");
190
191 AliInfo(Form("PID parameters, Pi0zero10to60 :\n"));
192 for(Int_t i = 0; i < 6; i++){
193 for(Int_t j = 0; j < 6; j++){
194 printf(" %f, ", fPiZero10to60[i][j]);
195 }
196 printf("\n");
197 }
198
199 printf("\n");
200
3a8be91c 201}
8ba062b1 202