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