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 **************************************************************************/
17 //_________________________________________________________________________
18 // Class for the management by the CPV reconstruction.
20 //*-- Author : Boris Polichtchouk (IHEP, Protvino) 6 Mar 2001
22 // --- ROOT system ---
24 // --- Standard library ---
26 // --- AliRoot header files ---
28 #include "AliPHOSRecCpvManager.h"
29 #include "AliPHOSGeometry.h"
30 #include "AliPHOSLoader.h"
32 ClassImp(AliPHOSRecCpvManager)
34 //____________________________________________________________________________
36 AliPHOSRecCpvManager::AliPHOSRecCpvManager()
40 fOneGamChisqCut = 3.; // If Chi2/dof > fOneGamChisqCut, split point.
42 fOneGamInitialStep = 0.00005;
44 fOneGamStepMin = 0.0005;
45 fOneGamNumOfIterations = 50;
47 fTwoGamInitialStep = 0.00005;
50 fTwoGamStepMin = 0.00005;
51 fTwoGamNumOfIterations = 50;
53 // fThr0 = 0.0285; // Min. energy of rec. point. If E<fThr0, point deleted.
54 // fSqdCut = 3.; // Min. distance (in cm) between two rec points.
59 SetTitle("Cpv Reconstruction Manager");
62 AliPHOSRecCpvManager::~AliPHOSRecCpvManager(void)
67 Float_t AliPHOSRecCpvManager::Dispersion(Float_t etot, Float_t ai) const
69 //"Dispresion" of energy deposition in the cell.
70 // etot is the total shower energy, ai is the
71 // calculated cell response,
72 // ei is the measured cell response.
74 const Float_t kConst = 1.5;
75 return kConst*ai*(1.-ai/etot);
78 Float_t AliPHOSRecCpvManager::OneGamChi2(Float_t ai, Float_t ei, Float_t etot, Float_t& Gi) const
80 //"Chi2" for one cell.
81 // etot is the total "shower" energy, ai is the
82 // calculated cell response,
83 // ei is the measured cell response.
85 const Float_t kConst = 1.5;
88 Float_t d = kConst*ai*(1.-ai/etot);
91 Gi = dd*(2.- dd*kConst*(1.-2.*ai/etot));
93 Info("OneGamChi2", " OneGamChi2 (ai,ei,etot,&Gi,chi2) %f %f %f %f %f", ai, ei, etot, Gi, da*da/d );
99 Float_t AliPHOSRecCpvManager::TwoGamChi2(Float_t ai, Float_t ei, Float_t etot, Float_t& gi) const
101 // Put a comment here
103 const Float_t kConst = 1.5;
105 Float_t da = ai - ei;
106 Float_t d = kConst*ai*(1.-ai/etot);
109 gi = dd*(2.- dd*kConst*(1.-2.*ai/etot));
115 void AliPHOSRecCpvManager::AG(Float_t ei, Float_t xi, Float_t yi, Float_t& ai, Float_t& gxi, Float_t& gyi )
117 //Calculates amplitude (ai) and gradients (gxi, gyi) of CPV pad response.
118 //Integrated response (total "shower energy") is e,
119 //xi and yi are the distances along x and y from reference point
120 // to the pad center.
122 AliPHOSGeometry * geom = AliPHOSLoader::GetPHOSGeometry();
124 Float_t celZ = geom->GetPadSizeZ();
125 Float_t celY = geom->GetPadSizePhi();
127 // // Info("AG", "celZ: %f celY: %f", celZ, celY) ;
129 Float_t dx = celZ/2.;
130 Float_t dy = celY/2.;
132 // // Float_t x = xi*celZ;
133 // // Float_t y = yi*celZ;
140 Float_t a = Fcml(x+dx,y+dy) - Fcml(x+dx,y-dy) - Fcml(x-dx,y+dy) + Fcml(x-dx,y-dy);
144 Float_t gx = GradX(x+dx,y+dy) - GradX(x+dx,y-dy) - GradX(x-dx,y+dy) + GradX(x-dx,y-dy);
147 Float_t gy = GradY(x+dx,y+dy) - GradY(x+dx,y-dy) - GradY(x-dx,y+dy) + GradY(x-dx,y-dy);
152 Float_t AliPHOSRecCpvManager::Fcml(Float_t x, Float_t y)
154 //Cumulative function
156 const Float_t ka = 1.0;
157 const Float_t kb = 0.70;
159 Float_t fff = TMath::ATan(x*y/( kb*TMath::Sqrt( (kb*kb) + x*x+y*y)))
160 - TMath::ATan(x*y/(3*kb*TMath::Sqrt((3*kb)*(3*kb) + x*x+y*y)))
161 + TMath::ATan(x*y/(5*kb*TMath::Sqrt((5*kb)*(5*kb) + x*x+y*y)))
162 - TMath::ATan(x*y/(7*kb*TMath::Sqrt((7*kb)*(7*kb) + x*x+y*y)))
163 + TMath::ATan(x*y/(9*kb*TMath::Sqrt((9*kb)*(9*kb) + x*x+y*y)));
165 Float_t fcml = ka*fff/TMath::TwoPi();
166 // Info("Fcml", "fcml: %f", fcml) ;
172 Float_t AliPHOSRecCpvManager::GradX(Float_t x, Float_t y)
174 // Put a comment here
176 const Float_t ka = 1.0;
177 const Float_t kb = 0.70;
179 Float_t skv = kb*kb + x*x + y*y;
181 Float_t gradient = y*(1.-x*x/skv)* kb*TMath::Sqrt(skv)/( kb*kb*skv+x*x*y*y)
182 - y*(1.-x*x/skv)*3*kb*TMath::Sqrt(skv)/((3*kb)*(3*kb)*skv+x*x*y*y)
183 + y*(1.-x*x/skv)*5*kb*TMath::Sqrt(skv)/((5*kb)*(5*kb)*skv+x*x*y*y)
184 - y*(1.-x*x/skv)*7*kb*TMath::Sqrt(skv)/((7*kb)*(7*kb)*skv+x*x*y*y)
185 + y*(1.-x*x/skv)*9*kb*TMath::Sqrt(skv)/((9*kb)*(9*kb)*skv+x*x*y*y);
187 Float_t grad = ka*gradient/TMath::TwoPi();
192 Float_t AliPHOSRecCpvManager::GradY(Float_t x, Float_t y)
194 // Put a comment here
196 const Float_t ka = 1.0;
197 const Float_t kb = 0.70;
199 Float_t skv = kb*kb + x*x + y*y;
200 Float_t gradient = x*(1.-y*y/skv)* kb*TMath::Sqrt(skv)/( kb*kb*skv+x*x*y*y)
201 - x*(1.-y*y/skv)*3*kb*TMath::Sqrt(skv)/((3*kb)*(3*kb)*skv+x*x*y*y)
202 + x*(1.-y*y/skv)*5*kb*TMath::Sqrt(skv)/((5*kb)*(5*kb)*skv+x*x*y*y)
203 - x*(1.-y*y/skv)*7*kb*TMath::Sqrt(skv)/((7*kb)*(7*kb)*skv+x*x*y*y)
204 + x*(1.-y*y/skv)*9*kb*TMath::Sqrt(skv)/((9*kb)*(9*kb)*skv+x*x*y*y);
206 Float_t grad = ka*gradient/TMath::TwoPi();