syst. check that reads the MC information for the case of TPC-only
[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 fi, Float_t& gi) const
68 {
69   // Chi2 used in OneGam (one-gamma fitting).
70   // gi is d(Chi2)/d(ai).
71
72   fi = 0 ; 
73   Float_t da = ai - ei;
74   Float_t d = ei; // we assume that sigma(E) = sqrt(E)
75   gi = 2.*(ai-ei)/d;
76
77   return da*da/d;
78
79 }
80
81 Float_t AliPHOSRecEmcManager::TwoGamChi2(Float_t ai, Float_t ei, Float_t fi, Float_t& gi) const
82 {
83   // calculates chi^2
84   fi = 0 ; 
85   Float_t da = ai - ei;
86   Float_t d = ei; // we assume that sigma(E) = sqrt(E)
87   gi = 2.*(ai-ei)/d;
88
89   return da*da/d;
90
91 }
92
93 void AliPHOSRecEmcManager::AG(Float_t ei, Float_t xi, Float_t yi, Float_t& ai, Float_t& gxi, Float_t& gyi )
94 {
95   //Calculates amplitude (ai) and gradients (gxi, gyi) of CPV pad response.
96   //Integrated response (total "shower energy") is E, 
97   //xi and yi are the distances along x and y from reference point 
98   // to the pad center.
99   //Shape of the shower is from PHOS TDR.
100
101
102   Float_t r = TMath::Sqrt(xi*xi + yi*yi);
103   Float_t r4    = r*r*r*r ;
104   Float_t r295  = TMath::Power(r, 2.95) ;
105   Float_t shape = ei*TMath::Exp( -r4 * (1. / (2.32 + 0.26 * r4) + 0.0316 / (1 + 0.0652 * r295) ) ) ;
106   ai = shape;
107
108
109   //d(shape)/d(xi)
110   gxi = (-(TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),2)*
111        ((-0.006077944*xi*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),
112              0.4750000000000001))/
113           TMath::Power(1 + 0.0652*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),1.475),2) - 
114          (1.04*xi*(TMath::Power(xi,2) + TMath::Power(yi,2)))/
115           TMath::Power(2.32 + 0.26*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),2),2))) - 
116     4*xi*(TMath::Power(xi,2) + TMath::Power(yi,2))*
117      (0.0316/(1 + 0.0652*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),1.475)) + 
118        1./(2.32 + 0.26*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),2))))/
119   TMath::Power(TMath::E(),TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),2)*
120     (0.0316/(1 + 0.0652*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),1.475)) + 
121      1./(2.32 + 0.26*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),2))));
122   
123   gxi = gxi*ei;
124
125   //d(shape)/d(yi)
126   gyi = (-(TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),2)*
127        ((-0.006077944*yi*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),
128              0.4750000000000001))/
129           TMath::Power(1 + 0.0652*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),1.475),2) - 
130          (1.04*yi*(TMath::Power(xi,2) + TMath::Power(yi,2)))/
131           TMath::Power(2.32 + 0.26*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),2),2))) - 
132     4*yi*(TMath::Power(xi,2) + TMath::Power(yi,2))*
133      (0.0316/(1 + 0.0652*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),1.475)) + 
134        1./(2.32 + 0.26*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),2))))/
135   TMath::Power(TMath::E(),TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),2)*
136     (0.0316/(1 + 0.0652*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),1.475)) + 
137      1./(2.32 + 0.26*TMath::Power(TMath::Power(xi,2) + TMath::Power(yi,2),2))));
138
139
140   gyi = gyi*ei;
141
142 }
143
144
145
146
147
148
149