1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
18 //_________________________________________________________________________
19 // Class for the management by the Emc reconstruction.
21 //*-- Author : Boris Polichtchouk (IHEP, Protvino) 6 Mar 2001
23 // --- ROOT system ---
27 // --- Standard library ---
29 // --- AliRoot header files ---
31 #include "AliPHOSRecEmcManager.h"
33 ClassImp(AliPHOSRecEmcManager)
35 //____________________________________________________________________________
37 AliPHOSRecEmcManager::AliPHOSRecEmcManager():
38 fOneGamChisqCut(1.3f),
39 fOneGamInitialStep(0.00005f),
41 fOneGamStepMin(0.0005f),
42 fOneGamNumOfIterations(50),
43 fTwoGamInitialStep(0.00005f),
46 fTwoGamStepMin(0.00005),
47 fTwoGamNumOfIterations(50),
52 SetTitle("Emc Reconstruction Manager");
55 AliPHOSRecEmcManager::~AliPHOSRecEmcManager(void) {}
57 Float_t AliPHOSRecEmcManager::Dispersion(Float_t ei) const
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.
67 Float_t AliPHOSRecEmcManager::OneGamChi2(Float_t ai, Float_t ei, Float_t fi, Float_t& gi) const
69 // Chi2 used in OneGam (one-gamma fitting).
70 // gi is d(Chi2)/d(ai).
74 Float_t d = ei; // we assume that sigma(E) = sqrt(E)
81 Float_t AliPHOSRecEmcManager::TwoGamChi2(Float_t ai, Float_t ei, Float_t fi, Float_t& gi) const
86 Float_t d = ei; // we assume that sigma(E) = sqrt(E)
93 void AliPHOSRecEmcManager::AG(Float_t ei, Float_t xi, Float_t yi, Float_t& ai, Float_t& gxi, Float_t& gyi )
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
99 //Shape of the shower is from PHOS TDR.
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) ) ) ;
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))));
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))));