Fix
[u/mrichter/AliRoot.git] / PHOS / AliPHOSRecEmcManager.cxx
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
18 //_________________________________________________________________________
19 // Class for the management by the Emc reconstruction.
20 ////                  
21 //*-- Author   : Boris Polichtchouk (IHEP, Protvino) 6 Mar 2001
22 //
23 // --- ROOT system ---
24
25 #include <TMath.h>
26
27 // --- Standard library ---
28
29 // --- AliRoot header files ---
30
31 #include "AliPHOSRecEmcManager.h"
32
33 ClassImp(AliPHOSRecEmcManager) 
34
35 //____________________________________________________________________________
36
37 AliPHOSRecEmcManager::AliPHOSRecEmcManager():
38   fOneGamChisqCut(1.3f),
39   fOneGamInitialStep(0.00005f),
40   fOneGamChisqMin(1.f),
41   fOneGamStepMin(0.0005f),
42   fOneGamNumOfIterations(50),
43   fTwoGamInitialStep(0.00005f),
44   fTwoGamChisqMin(1.f),
45   fTwoGamEmin(0.1f),
46   fTwoGamStepMin(0.00005),
47   fTwoGamNumOfIterations(50),
48   fThr0(0.f),
49   fSqdCut(0.f)
50 {
51   // default ctor
52   SetTitle("Emc Reconstruction Manager");
53 }
54
55 AliPHOSRecEmcManager::~AliPHOSRecEmcManager(void) {}
56
57 Float_t AliPHOSRecEmcManager::Dispersion(Float_t ei) const
58 {
59   //"Dispresion" of energy deposition in the cell.
60   // eTot is the total shower energy, ai is the
61   // calculated cell response,
62   // ei is the measured cell response.
63
64   return ei;
65 }
66
67 Float_t AliPHOSRecEmcManager::OneGamChi2(Float_t ai, Float_t ei, Float_t, Float_t& gi) const
68 {
69   // Chi2 used in OneGam (one-gamma fitting).
70   // gi is d(Chi2)/d(ai).
71
72   Float_t da = ai - ei;
73   Float_t d = ei; // we assume that sigma(E) = sqrt(E)
74   gi = 2.*(ai-ei)/d;
75
76   return da*da/d;
77
78 }
79
80 Float_t AliPHOSRecEmcManager::TwoGamChi2(Float_t ai, Float_t ei, Float_t, Float_t& gi) const
81 {
82   // calculates chi^2
83   Float_t da = ai - ei;
84   Float_t d = ei; // we assume that sigma(E) = sqrt(E)
85   gi = 2.*(ai-ei)/d;
86
87   return da*da/d;
88
89 }
90
91 void AliPHOSRecEmcManager::AG(Float_t ei, Float_t xi, Float_t yi, Float_t& ai, Float_t& gxi, Float_t& gyi )
92 {
93   //Calculates amplitude (ai) and gradients (gxi, gyi) of CPV pad response.
94   //Integrated response (total "shower energy") is E, 
95   //xi and yi are the distances along x and y from reference point 
96   // to the pad center.
97   //Shape of the shower is from PHOS TDR.
98
99
100   Float_t r = TMath::Sqrt(xi*xi + yi*yi);
101   Float_t r4    = r*r*r*r ;
102   Float_t r295  = TMath::Power(r, 2.95) ;
103   Float_t shape = ei*TMath::Exp( -r4 * (1. / (2.32 + 0.26 * r4) + 0.0316 / (1 + 0.0652 * r295) ) ) ;
104   ai = shape;
105
106
107   //d(shape)/d(xi)
108   gxi = (-(TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),2)*
109        ((-0.006077944*xi*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),
110              0.4750000000000001))/
111           TMath::Power(1 + 0.0652*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),1.475),2) - 
112          (1.04*xi*(TMath::Power(xi,2) + TMath::Power(yi,2)))/
113           TMath::Power(2.32 + 0.26*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),2),2))) - 
114     4*xi*(TMath::Power(xi,2) + TMath::Power(yi,2))*
115      (0.0316/(1 + 0.0652*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),1.475)) + 
116        1./(2.32 + 0.26*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),2))))/
117   TMath::Power(TMath::E(),TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),2)*
118     (0.0316/(1 + 0.0652*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),1.475)) + 
119      1./(2.32 + 0.26*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),2))));
120   
121   gxi = gxi*ei;
122
123   //d(shape)/d(yi)
124   gyi = (-(TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),2)*
125        ((-0.006077944*yi*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),
126              0.4750000000000001))/
127           TMath::Power(1 + 0.0652*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),1.475),2) - 
128          (1.04*yi*(TMath::Power(xi,2) + TMath::Power(yi,2)))/
129           TMath::Power(2.32 + 0.26*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),2),2))) - 
130     4*yi*(TMath::Power(xi,2) + TMath::Power(yi,2))*
131      (0.0316/(1 + 0.0652*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),1.475)) + 
132        1./(2.32 + 0.26*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),2))))/
133   TMath::Power(TMath::E(),TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),2)*
134     (0.0316/(1 + 0.0652*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),1.475)) + 
135      1./(2.32 + 0.26*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),2))));
136
137
138   gyi = gyi*ei;
139
140 }
141
142
143
144
145
146
147