BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
EmcRecShowerPosLog.cxx
Go to the documentation of this file.
1//
2// Logarithmic weighted position attribute reconstruction
3//
4// Created by Zhe Wang 2004, 5, 16
5//
6#include "EmcRec/EmcRecShowerPosLog.h"
7#include "EmcRec/EmcRecParameter.h"
8
9#include <fstream>
10#include <iostream>
11
12#include "CLHEP/Matrix/Matrix.h"
13#include "EmcRecGeoSvc/IEmcRecGeoSvc.h"
14#include "GaudiKernel/Bootstrap.h"
15#include "GaudiKernel/ISvcLocator.h"
16
17
18
19void EmcRecShowerPosLog::Position( RecEmcShower& aShower ) {
20 RecEmcFractionMap::const_iterator cit;
21 HepPoint3D pos( 0, 0, 0 );
22 HepPoint3D possum( 0, 0, 0 );
23 double w, wsum = 0;
24
25 // cout<<"EmcRecShowerPosLog::Position Begin"<<endl;
26
27 IEmcRecGeoSvc* iGeoSvc;
28 ISvcLocator* svcLocator = Gaudi::svcLocator();
29 StatusCode sc = svcLocator->service( "EmcRecGeoSvc", iGeoSvc );
30 if ( sc != StatusCode::SUCCESS )
31 {
32 // cout<<"Error: Can't get EmcRecGeoSvc"<<endl;
33 }
34
35 EmcRecParameter& Para = EmcRecParameter::GetInstance();
36
37 double offset = Para.LogPosOffset();
38 // cout<<offset<<endl;
39
40 double e5x5 = -99999;
41 if ( Para.PositionMode2() == "all" )
42 {
43 double etot = aShower.getEAll();
44 for ( cit = aShower.Begin(); cit != aShower.End(); cit++ )
45 {
46 w = offset + log( ( cit->second.getEnergy() / etot ) * cit->second.getFraction() );
47 if ( w > 0 )
48 {
49 wsum += w;
50 pos = iGeoSvc->GetCFrontCenter( cit->second.getCellId() );
51 possum += pos * w;
52 }
53 }
54 }
55 else if ( Para.PositionMode2() == "3x3" )
56 {
57 double e3x3 = aShower.e3x3();
58 RecEmcFractionMap fracMap3x3 = aShower.getFractionMap3x3();
59 for ( cit = fracMap3x3.begin(); cit != fracMap3x3.end(); cit++ )
60 {
61 w = offset + log( ( cit->second.getEnergy() / e3x3 ) * cit->second.getFraction() );
62 if ( w > 0 )
63 {
64 wsum += w;
65 pos = iGeoSvc->GetCFrontCenter( cit->second.getCellId() );
66 possum += pos * w;
67 }
68 }
69 }
70 else
71 {
72 e5x5 = aShower.e5x5();
73 RecEmcFractionMap fracMap5x5 = aShower.getFractionMap5x5();
74 for ( cit = fracMap5x5.begin(); cit != fracMap5x5.end(); cit++ )
75 {
76 w = offset + log( ( cit->second.getEnergy() / e5x5 ) * cit->second.getFraction() );
77 if ( w > 0 )
78 {
79 wsum += w;
80 pos = iGeoSvc->GetCFrontCenter( cit->second.getCellId() );
81 possum += pos * w;
82 }
83 }
84 }
85 if ( wsum <= 0 )
86 {
87 // cout<<"[[[[[[[[[[[[[[["<<wsum<<endl;
88 }
89 pos = possum / wsum;
90
91 // aShower.setPositionNoCor(pos);
92 //----------------------------------------------------------------------
93 // position correction
94 RecEmcID id = aShower.getShowerId();
95 const unsigned int module = EmcID::barrel_ec( id );
96 const unsigned int thetaModule = EmcID::theta_module( id );
97 const unsigned int phiModule = EmcID::phi_module( id );
98 HepPoint3D posCorr( pos );
99
100 if ( module == 1 )
101 { // barrel
102
103 double IRShift = 0, parTheta[5], parPhi[5];
104
105 unsigned int theModule;
106 if ( thetaModule > 21 )
107 {
108 theModule = 43 - thetaModule;
109 id = EmcID::crystal_id( module, theModule, phiModule );
110 posCorr.setZ( -posCorr.z() );
111 }
112 else { theModule = thetaModule; }
113
114 if ( theModule == 21 ) { IRShift = 0; }
115 else if ( theModule == 20 ) { IRShift = 2.5; }
116 else { IRShift = 5.; }
117 HepPoint3D IR( 0, 0, IRShift );
118
119 // MethodMode=0 old parameters
120 // MethodMode=1 new parameters
121 if ( Para.MethodMode() == 0 )
122 { // using old parameter given by hemiao
123
124 for ( int i = 0; i < 5; i++ )
125 {
126 if ( theModule == 21 )
127 {
128 double par[5] = { 51.97, -0.1209, 6.175, -0.1431, -0.005497 };
129 parTheta[i] = par[i];
130 }
131 else if ( theModule == 20 )
132 {
133 double par[5] = { 56.31, -0.1135, 6.312, -1.216, -0.02978 };
134 parTheta[i] = par[i];
135 }
136 else
137 {
138 double par[5] = { 10.65, -0.2257, 2.216, -0.1562, -0.05552 };
139 parTheta[i] = par[i];
140 }
141 }
142 }
143 else if ( Para.MethodMode() == 1 )
144 { // using new parameter given by liuchunxiu
145
146 for ( int i = 0; i < 5; i++ )
147 {
148
149 if ( e5x5 > 1.0 ) parTheta[i] = Para.BarrLogThetaPara( theModule, i );
150 else if ( e5x5 <= 1.0 && e5x5 > 0.5 )
151 parTheta[i] = Para.BarrLogThetaPara( theModule + 22, i );
152 else if ( e5x5 <= 0.5 ) parTheta[i] = Para.BarrLogThetaPara( theModule + 44, i );
153
154 if ( e5x5 > 1.0 ) parPhi[i] = Para.BarrLogPhiPara( theModule, i );
155 else if ( e5x5 <= 1.0 && e5x5 > 0.5 )
156 parPhi[i] = Para.BarrLogPhiPara( theModule + 22, i );
157 else if ( e5x5 <= 0.5 ) parPhi[i] = Para.BarrLogPhiPara( theModule + 44, i );
158 }
159 }
160
161 HepPoint3D center01, center23; // center of point0,1 and point2,3 in crystal
162 center01 = ( iGeoSvc->GetCrystalPoint( id, 0 ) + iGeoSvc->GetCrystalPoint( id, 1 ) ) / 2;
163 center23 = ( iGeoSvc->GetCrystalPoint( id, 2 ) + iGeoSvc->GetCrystalPoint( id, 3 ) ) / 2;
164
165 double theta01, theta23, length, d;
166 theta01 = ( center01 - IR ).theta();
167 theta23 = ( center23 - IR ).theta();
168 length = ( center01 - IR ).mag();
169 d = length * tan( theta01 - theta23 ); // length in x direction
170
171 double xIn, xOut, deltaTheta;
172 xIn = length * tan( theta01 - ( posCorr - IR ).theta() ) - d / 2;
173 xOut = xIn - ( parTheta[0] * atan( xIn * parTheta[1] - parTheta[4] ) + parTheta[2] * xIn +
174 parTheta[3] );
175 deltaTheta = atan( ( xOut + d / 2 ) / length ) - atan( ( xIn + d / 2 ) / length );
176
177 // obtained by photon, not used, just for comparison
178 // double parPhi[5] = { 51.1, -0.1499, 7.62, 0.1301, 0.005581 };
179 double yIn, deltaY, deltaPhi;
180
181 // 1.843=1.5+0.343 degree
182 if ( phiModule == 0 ) { yIn = posCorr.phi() * 180. / CLHEP::pi - phiModule * 3. - 1.843; }
183 else if ( phiModule == 119 )
184 { yIn = posCorr.phi() * 180. / CLHEP::pi + 360. - phiModule * 3. - 1.843; }
185 else
186 {
187 yIn = posCorr.phi() < 0
188 ? posCorr.phi() * 180. / CLHEP::pi + 360. - phiModule * 3. - 1.843
189 : posCorr.phi() * 180. / CLHEP::pi - phiModule * 3. - 1.843;
190 }
191
192 // deltaY=yIn-Ymc => Ymc=yIn-deltaY
193 deltaY = parPhi[0] * atan( yIn * parPhi[1] - parPhi[4] ) + parPhi[2] * yIn + parPhi[3];
194 deltaPhi = deltaY * CLHEP::pi / 180.;
195
196 HepPoint3D rotateVector( 0, 1, 0 ); // y axis
197 rotateVector.rotateZ( center01.phi() );
198 posCorr.setZ( posCorr.z() - IRShift );
199 posCorr.rotate( -deltaTheta, rotateVector );
200 posCorr.setZ( posCorr.z() + IRShift );
201 // phi parameter is gotten by the average of e+ e- peak
202 if ( Para.MethodMode() == 0 ) { posCorr.rotateZ( -0.002994 ); }
203 else if ( Para.MethodMode() == 1 ) { posCorr.rotateZ( -deltaPhi ); }
204
205 if ( thetaModule > 21 ) { posCorr.setZ( -posCorr.z() ); }
206
207 double lengthemc;
208 lengthemc = abs( posCorr.z() / cos( posCorr.theta() ) );
209 // for Data
210 double posDataCorPar;
211 if ( Para.DataMode() == 1 )
212 {
213 posDataCorPar = Para.BarrPosDataCor( thetaModule, phiModule );
214 posCorr.setTheta( posCorr.theta() - posDataCorPar / lengthemc );
215 }
216
217 // for MC
218 double posMCCorPar;
219 if ( Para.DataMode() == 0 )
220 {
221 posMCCorPar = Para.BarrPosMCCor( thetaModule, phiModule );
222 posCorr.setTheta( posCorr.theta() - posMCCorPar / lengthemc );
223 }
224 }
225 else
226 { // endcap position corretion
227 if ( Para.MethodMode() == 1 )
228 {
229
230 double IRShift = 0, parTheta[5], parPhi[5];
231
232 if ( module == 0 )
233 { // east endcap
234
235 IRShift = 10;
236
237 for ( int i = 0; i < 5; i++ )
238 {
239
240 if ( e5x5 > 1.0 ) parTheta[i] = Para.EastLogThetaPara( thetaModule, i );
241 else if ( e5x5 <= 1.0 && e5x5 > 0.5 )
242 parTheta[i] = Para.EastLogThetaPara( thetaModule + 6, i );
243 else if ( e5x5 <= 0.5 ) parTheta[i] = Para.EastLogThetaPara( thetaModule + 6, i );
244
245 if ( e5x5 > 1.0 ) parPhi[i] = Para.EastLogPhiPara( 0, i );
246 else if ( e5x5 <= 1.0 && e5x5 > 0.5 ) parPhi[i] = Para.EastLogPhiPara( 1, i );
247 else if ( e5x5 <= 0.5 ) parPhi[i] = Para.EastLogPhiPara( 2, i );
248 }
249 }
250 else if ( module == 2 )
251 { // west endcap
252 IRShift = -10;
253 for ( int i = 0; i < 5; i++ )
254 {
255
256 if ( e5x5 > 1.0 ) parTheta[i] = Para.WestLogThetaPara( thetaModule, i );
257 else if ( e5x5 <= 1.0 && e5x5 > 0.5 )
258 parTheta[i] = Para.WestLogThetaPara( thetaModule + 6, i );
259 else if ( e5x5 <= 0.5 ) parTheta[i] = Para.WestLogThetaPara( thetaModule + 6, i );
260
261 if ( e5x5 > 1.0 ) parPhi[i] = Para.WestLogPhiPara( 0, i );
262 else if ( e5x5 <= 1.0 && e5x5 > 0.5 ) parPhi[i] = Para.WestLogPhiPara( 1, i );
263 else if ( e5x5 <= 0.5 ) parPhi[i] = Para.WestLogPhiPara( 2, i );
264 }
265 }
266
267 HepPoint3D IR( 0, 0, IRShift );
268
269 double theta12 = -9999, theta03 = -9999, theta23 = -9999, theta14 = -9999, delta = -9999;
270 double phi01 = -9999, phi23 = -9999, phi12 = -9999, phi34 = -9999, deltaphi = -9999;
271 double xIn, deltaX, deltaTheta, yIn, deltaY, deltaPhi;
272
273 double posphi = posCorr.phi();
274 if ( phiModule == 0 ) { posphi = posphi; }
275 else if ( ( ( thetaModule == 0 || thetaModule == 1 ) && phiModule == 63 ) ||
276 ( ( thetaModule == 2 || thetaModule == 3 ) && phiModule == 79 ) ||
277 ( ( thetaModule == 4 || thetaModule == 5 ) && phiModule == 95 ) )
278 { posphi = posphi + 2 * CLHEP::pi; }
279 else { posphi = posphi < 0 ? posphi + 2 * CLHEP::pi : posphi; }
280
281 HepPoint3D center, center01, center23, center12, center03, center34, center14;
282
283 if ( ( thetaModule == 1 &&
284 ( ( phiModule + 3 ) % 4 == 0 || ( phiModule + 2 ) % 4 == 0 ) ) ||
285 ( thetaModule == 3 && ( ( phiModule + 4 ) % 5 == 0 || ( phiModule + 3 ) % 5 == 0 ||
286 ( phiModule + 2 ) % 5 == 0 ) ) ) // pentagon
287 {
288
289 center23 =
290 ( iGeoSvc->GetCrystalPoint( id, 2 ) + iGeoSvc->GetCrystalPoint( id, 3 ) ) / 2.0;
291 center14 =
292 ( iGeoSvc->GetCrystalPoint( id, 1 ) + iGeoSvc->GetCrystalPoint( id, 4 ) ) / 2.0;
293 center = center23;
294 theta23 = ( center23 - IR ).theta();
295 theta14 = ( center14 - IR ).theta();
296 delta = theta14 - theta23;
297 xIn = ( ( ( posCorr - IR ).theta() - theta23 ) - delta * 0.5 ) / delta;
298
299 center12 =
300 ( iGeoSvc->GetCrystalPoint( id, 1 ) + iGeoSvc->GetCrystalPoint( id, 2 ) ) / 2.0;
301 center34 =
302 ( iGeoSvc->GetCrystalPoint( id, 3 ) + iGeoSvc->GetCrystalPoint( id, 4 ) ) / 2.0;
303 phi12 = center12.phi();
304 phi34 = center34.phi();
305
306 phi12 = phi12 < 0 ? phi12 + 2 * CLHEP::pi : phi12;
307 phi34 = phi34 < 0 ? phi34 + 2 * CLHEP::pi : phi34;
308 deltaphi = phi34 - phi12;
309 yIn = ( ( posphi - phi12 ) - deltaphi * 0.5 ) / deltaphi;
310 }
311 else
312 {
313
314 center12 =
315 ( iGeoSvc->GetCrystalPoint( id, 1 ) + iGeoSvc->GetCrystalPoint( id, 2 ) ) / 2.0;
316 center03 =
317 ( iGeoSvc->GetCrystalPoint( id, 0 ) + iGeoSvc->GetCrystalPoint( id, 3 ) ) / 2.0;
318 center = center12;
319
320 theta12 = ( center12 - IR ).theta();
321 theta03 = ( center03 - IR ).theta();
322 delta = theta03 - theta12;
323 xIn = ( ( ( posCorr - IR ).theta() - theta12 ) - delta * 0.5 ) / delta;
324
325 center01 =
326 ( iGeoSvc->GetCrystalPoint( id, 0 ) + iGeoSvc->GetCrystalPoint( id, 1 ) ) / 2.0;
327 center23 =
328 ( iGeoSvc->GetCrystalPoint( id, 2 ) + iGeoSvc->GetCrystalPoint( id, 3 ) ) / 2.0;
329 phi01 = center01.phi();
330 phi23 = center23.phi();
331 phi01 = phi01 < 0 ? phi01 + 2 * CLHEP::pi : phi01;
332 phi23 = phi23 < 0 ? phi23 + 2 * CLHEP::pi : phi23;
333 deltaphi = phi23 - phi01;
334 yIn = ( ( posphi - phi01 ) - deltaphi * 0.5 ) / deltaphi;
335 }
336 // deltaX=xIn-Xext => Xext=xIn-deltaX
337 deltaX = parTheta[0] * atan( xIn * parTheta[1] - parTheta[4] ) + parTheta[2] * xIn +
338 parTheta[3];
339 deltaTheta = deltaX * delta;
340 // deltaY=yIn-Ymc => Ymc=yIn-deltaY
341 deltaY = parPhi[0] * atan( yIn * parPhi[1] - parPhi[4] ) + parPhi[2] * yIn + parPhi[3];
342 deltaPhi = deltaY * deltaphi;
343
344 HepPoint3D rotateVector( 0, 1, 0 ); // y axis
345 rotateVector.rotateZ( center.phi() );
346 posCorr.setZ( posCorr.z() - IRShift );
347 posCorr.rotate( -deltaTheta, rotateVector );
348 posCorr.setZ( posCorr.z() + IRShift );
349 posCorr.rotateZ( -deltaPhi );
350
351 double lengthemc;
352 lengthemc = abs( posCorr.z() / cos( posCorr.theta() ) );
353 double posDataCorPar = 0.0;
354 if ( Para.DataMode() == 1 )
355 {
356 if ( module == 0 ) posDataCorPar = Para.EastPosDataCor( thetaModule, phiModule );
357 if ( module == 2 ) posDataCorPar = Para.WestPosDataCor( thetaModule, phiModule );
358 posCorr.setTheta( posCorr.theta() - posDataCorPar / lengthemc );
359 }
360 }
361 }
362
363 // PosCorr=0 without position correction
364 // PosCorr=1 with position correction
365 if ( Para.PosCorr() == 0 ) { aShower.setPosition( pos ); }
366
367 if ( Para.PosCorr() == 1 ) { aShower.setPosition( posCorr ); }
368 ////////////////////////////
369 if ( aShower.energy() < 1e-5 ) return;
370
371 double r, theta, phi;
372 double dtheta, dphi, dx, dy, dz;
373
374 r = posCorr.mag();
375 theta = posCorr.theta();
376 phi = posCorr.phi();
377 // cout<<"energy: "<<aShower.energy()<<" theta: "<<theta<<" phi: "<<phi<<endl;
378
379 if ( aShower.energy() > 0 )
380 {
381 dtheta = Para.SigTheta( 0 ) / sqrt( aShower.energy() ) + Para.SigTheta( 1 );
382 dphi = Para.SigPhi( 0 ) / sqrt( aShower.energy() ) + Para.SigPhi( 1 );
383 }
384 else
385 {
386 dtheta = 0;
387 dphi = 0;
388 }
389 // cout<<"dtheta: "<<dtheta<<" dphi: "<<dphi<<endl;
390
391 dx = fabs( iGeoSvc->GetBarrelR() * sin( phi ) * dphi );
392 dy = fabs( iGeoSvc->GetBarrelR() * cos( phi ) * dphi );
393 if ( theta > 0 )
394 dz = fabs( iGeoSvc->GetBarrelR() * dtheta / ( sin( theta ) * sin( theta ) ) );
395 else dz = 0;
396 // cout<<" dx: "<<dx<<" dy: "<<dy<<" dz: "<<dz<<endl;
397
398 aShower.setDtheta( dtheta );
399 aShower.setDphi( dphi );
400
401 //----------------------------------------------------------------------
402 // Error matrix
403 HepSymMatrix matrix( 3 ); // error matrix of r, theta, phi
404 matrix[0][0] = 0;
405 matrix[1][1] = dtheta * dtheta;
406 matrix[2][2] = dphi * dphi;
407 matrix[0][1] = 0;
408 matrix[0][2] = 0;
409 matrix[1][2] = 0;
410 // cout<<"matrix:\n"<<matrix<<endl;
411
412 HepMatrix matrix1( 3, 3 ), matrix2( 3, 3 );
413 matrix1[0][0] = sin( theta ) * cos( phi );
414 matrix1[0][1] = r * cos( theta ) * cos( phi );
415 matrix1[0][2] = -r * sin( theta ) * sin( phi );
416 matrix1[1][0] = sin( theta ) * sin( phi );
417 matrix1[1][1] = r * cos( theta ) * sin( phi );
418 matrix1[1][2] = r * sin( theta ) * cos( phi );
419 matrix1[2][0] = cos( theta );
420 matrix1[2][1] = -r * sin( theta );
421 matrix1[2][2] = 0;
422 // cout<<"matrix1:\n"<<matrix1<<endl;
423
424 for ( int i = 0; i < 3; i++ )
425 {
426 for ( int j = 0; j < 3; j++ ) { matrix2[i][j] = matrix1[j][i]; }
427 }
428 // cout<<"matrix2:\n"<<matrix2<<endl;
429
430 HepMatrix matrix4( 3, 3 );
431 matrix4 = matrix1 * matrix * matrix2;
432 // cout<<"matrix4:\n"<<matrix4<<endl;
433
434 HepSymMatrix matrix5( 3 ); // error matrix of x, y, z
435 for ( int i = 0; i < 3; i++ )
436 {
437 for ( int j = 0; j <= i; j++ ) { matrix5[i][j] = matrix4[i][j]; }
438 }
439 // cout<<"matrix5:\n"<<matrix5<<endl;
440
441 aShower.setErrorMatrix( matrix5 );
442
443 if ( matrix5[0][0] > 0 ) dx = sqrt( matrix5[0][0] );
444 if ( matrix5[1][1] > 0 ) dy = sqrt( matrix5[1][1] );
445 if ( matrix5[2][2] > 0 ) dz = sqrt( matrix5[2][2] );
446}
HepPoint3D
HepGeom::Point3D< double > HepPoint3D
Definition CosmicGenerator.h:36
etot
Double_t etot
Definition DataBase/tau_mode.c:8
RecEmcFractionMap
map< RecEmcID, RecEmcFraction, less< RecEmcID > > RecEmcFractionMap
Definition Emc/EmcRecEventModel/include/EmcRecEventModel/RecEmcFraction.h:55
w
double w
Definition EvtPhokharaDef.hh:28
tan
double tan(const BesAngle a)
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcGeom/BesAngle.h:189
sin
double sin(const BesAngle a)
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcGeom/BesAngle.h:185
cos
double cos(const BesAngle a)
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcGeom/BesAngle.h:187
DstEmcShower::setPosition
void setPosition(const HepPoint3D &pos)
Definition Event/DstEvent/include/DstEvent/DstEmcShower.h:62
DstEmcShower::setErrorMatrix
void setErrorMatrix(const HepSymMatrix &error)
Definition Event/DstEvent/include/DstEvent/DstEmcShower.h:75
EmcID::crystal_id
static Identifier crystal_id(const unsigned int barrel_ec, const unsigned int theta_module, const unsigned int phi_module)
For a single crystal.
Definition EmcID.cxx:63
EmcID::barrel_ec
static unsigned int barrel_ec(const Identifier &id)
Values of different levels (failure returns 0).
Definition EmcID.cxx:36
EmcID::theta_module
static unsigned int theta_module(const Identifier &id)
Definition EmcID.cxx:41
EmcID::phi_module
static unsigned int phi_module(const Identifier &id)
Definition EmcID.cxx:46
EmcRecParameter::LogPosOffset
double LogPosOffset() const
Definition EmcRecParameter.cxx:516
EmcRecParameter::BarrPosMCCor
double BarrPosMCCor(int ntheta, int nphi) const
Definition EmcRecParameter.cxx:641
EmcRecParameter::GetInstance
static EmcRecParameter & GetInstance()
Definition EmcRecParameter.cxx:494
EmcRecParameter::EastPosDataCor
double EastPosDataCor(int ntheta, int nphi) const
Definition EmcRecParameter.cxx:637
EmcRecParameter::WestPosDataCor
double WestPosDataCor(int ntheta, int nphi) const
Definition EmcRecParameter.cxx:633
EmcRecParameter::BarrLogPhiPara
double BarrLogPhiPara(int n, int m) const
Definition EmcRecParameter.cxx:553
EmcRecParameter::PosCorr
double PosCorr() const
Definition EmcRecParameter.cxx:524
EmcRecParameter::MethodMode
double MethodMode() const
Definition EmcRecParameter.cxx:522
EmcRecParameter::SigTheta
double SigTheta(int n) const
Definition EmcRecParameter.cxx:537
EmcRecParameter::WestLogThetaPara
double WestLogThetaPara(int n, int m) const
Definition EmcRecParameter.cxx:589
EmcRecParameter::BarrPosDataCor
double BarrPosDataCor(int ntheta, int nphi) const
Definition EmcRecParameter.cxx:629
EmcRecParameter::BarrLogThetaPara
double BarrLogThetaPara(int n, int m) const
Definition EmcRecParameter.cxx:549
EmcRecParameter::WestLogPhiPara
double WestLogPhiPara(int n, int m) const
Definition EmcRecParameter.cxx:595
EmcRecParameter::EastLogPhiPara
double EastLogPhiPara(int n, int m) const
Definition EmcRecParameter.cxx:593
EmcRecParameter::EastLogThetaPara
double EastLogThetaPara(int n, int m) const
Definition EmcRecParameter.cxx:585
EmcRecParameter::DataMode
double DataMode() const
Definition EmcRecParameter.cxx:526
EmcRecParameter::SigPhi
double SigPhi(int n) const
Definition EmcRecParameter.cxx:539
EmcRecShowerPosLog::Position
virtual void Position(RecEmcShower &aShower)
Definition EmcRecShowerPosLog.cxx:19
IEmcRecGeoSvc::GetBarrelR
virtual double GetBarrelR() const =0
IEmcRecGeoSvc::GetCFrontCenter
virtual HepPoint3D GetCFrontCenter(const Identifier &id) const =0
IEmcRecGeoSvc::GetCrystalPoint
virtual HepPoint3D GetCrystalPoint(const Identifier &id, const int i) const =0
RecEmcShower::End
RecEmcFractionMap::const_iterator End() const
Definition RecEmcShower.cxx:114
RecEmcShower::getFractionMap5x5
RecEmcFractionMap getFractionMap5x5() const
Definition RecEmcShower.cxx:206
RecEmcShower::Begin
RecEmcFractionMap::const_iterator Begin() const
Definition RecEmcShower.cxx:112
RecEmcShower::getFractionMap3x3
RecEmcFractionMap getFractionMap3x3() const
Definition RecEmcShower.cxx:183