**************************************************************************/
- //_________________________________________________________________________
+//_________________________________________________________________________
// Class for the management by the CPV reconstruction.
-//
+////
//*-- Author : Boris Polichtchouk (IHEP, Protvino) 6 Mar 2001
//
// --- ROOT system ---
// --- AliRoot header files ---
#include "AliPHOSRecCpvManager.h"
-#include "AliPHOSGetter.h"
+#include "AliPHOSGeometry.h"
+#include "AliPHOSLoader.h"
ClassImp(AliPHOSRecCpvManager)
// Put a comment here
}
-Float_t AliPHOSRecCpvManager::Dispersion(Float_t etot, Float_t ai, Float_t ei) const
+Float_t AliPHOSRecCpvManager::Dispersion(Float_t etot, Float_t ai) const
{
//"Dispresion" of energy deposition in the cell.
// etot is the total shower energy, ai is the
return kConst*ai*(1.-ai/etot);
}
-Float_t AliPHOSRecCpvManager::OneGamChi2(Float_t ai, Float_t ei, Float_t etot, Float_t& Gi)
+Float_t AliPHOSRecCpvManager::OneGamChi2(Float_t ai, Float_t ei, Float_t etot, Float_t& Gi) const
{
//"Chi2" for one cell.
// etot is the total "shower" energy, ai is the
//xi and yi are the distances along x and y from reference point
// to the pad center.
- AliPHOSGetter * gime = AliPHOSGetter::GetInstance() ;
- const AliPHOSGeometry* geom = gime->PHOSGeometry();
+ AliPHOSGeometry * geom = AliPHOSLoader::GetPHOSGeometry();
+
Float_t celZ = geom->GetPadSizeZ();
Float_t celY = geom->GetPadSizePhi();
- TMath::ATan(x*y/(7*kb*TMath::Sqrt((7*kb)*(7*kb) + x*x+y*y)))
+ TMath::ATan(x*y/(9*kb*TMath::Sqrt((9*kb)*(9*kb) + x*x+y*y)));
- Float_t fcml = ka*fff/6.2831853071796;
+ Float_t fcml = ka*fff/TMath::TwoPi();
// Info("Fcml", "fcml: %f", fcml) ;
return fcml;
- y*(1.-x*x/skv)*7*kb*TMath::Sqrt(skv)/((7*kb)*(7*kb)*skv+x*x*y*y)
+ y*(1.-x*x/skv)*9*kb*TMath::Sqrt(skv)/((9*kb)*(9*kb)*skv+x*x*y*y);
- Float_t grad = ka*gradient/6.2831853071796;
+ Float_t grad = ka*gradient/TMath::TwoPi();
return grad;
}
{
// Put a comment here
- const Float_t a = 1.0;
+ const Float_t ka = 1.0;
const Float_t kb = 0.70;
Float_t skv = kb*kb + x*x + y*y;
- x*(1.-y*y/skv)*7*kb*TMath::Sqrt(skv)/((7*kb)*(7*kb)*skv+x*x*y*y)
+ x*(1.-y*y/skv)*9*kb*TMath::Sqrt(skv)/((9*kb)*(9*kb)*skv+x*x*y*y);
- Float_t grad = a*gradient/6.2831853071796;
+ Float_t grad = ka*gradient/TMath::TwoPi();
return grad;
}