message += fTowerLocMaxCut ;
message += "\n EMC Logarothmic weight = " ;
message += fW0 ;
- message += "\n CPV Clustering threshold = " ;
+ message += "\n Pre Shower Clustering threshold = " ;
message += fPreShoClusteringThreshold ;
- message += "\n CPV Local Maximum cut = " ;
+ message += "\n Pre Shower Local Maximum cut = " ;
message += fPreShoLocMaxCut ;
- message += "\n CPV Logarothmic weight = " ;
+ message += "\n Pre Shower Logarothmic weight = " ;
message += fW0CPV ;
if(fToUnfold)
message +="\nUnfolding on\n" ;
fAirGap = 5.0; // cm, air gap between EMCAL mother volume and
// active material.
fAlFrontThick = 3.18; // cm, Thickness of front Al layer
- fPbRadThickness = 0.5; // cm, Thickness of theh Pb radiators.
+ fPbRadThickness = 0.5; // cm, Thickness of the Pb radiators.
fPreShowerSintThick = 0.6; // cm, Thickness of the sintilator for the
// preshower part of the calorimeter
fFullShowerSintThick = 0.5; // cm, Thickness of the sintilator for the
kDeg2Rad = TMath::Pi() / static_cast<Double_t>(180) ;
cyl_radius = GetIPDistance()+ GetAirGap() ;
x = cyl_radius * TMath::Cos(phi * kDeg2Rad ) ;
- y = cyl_radius * TMath::Cos(phi * kDeg2Rad ) ;
+ y = cyl_radius * TMath::Sin(phi * kDeg2Rad ) ;
z = cyl_radius / TMath::Tan(theta * kDeg2Rad ) ;
return;
TFile * file = static_cast<TFile*>(gROOT->GetFile(fRecPointsFileName));
if(!file)
file = TFile::Open(fRecPointsFileName) ;
-
+ Info("ReadTreeR", "file=%s\n",fRecPointsFileName.Data());
// Get Digits Tree header from file
TString treeName("TreeR") ;
void AliEMCALTowerRecPoint::EvalDispersion(Float_t logWeight,TClonesArray * digits)
{
// Calculates the dispersion of the shower at the origine of the RecPoint
- printf("**************** EVAL Dispersion *****************") ;
Float_t d = 0. ;
Float_t wtot = 0. ;
Float_t w = TMath::Max(0.,logWeight+TMath::Log(fEnergyList[iDigit]/fAmp ) ) ;
d += w*((xi-x)*(xi-x) + (yi-y)*(yi-y)+ (zi-z)*(zi-z) ) ;
wtot+=w ;
- printf("xi=%f, x=%f\n yi=%f, y=%f\n zi=%f, z=%f\n phi=%f phii=%f theta=%f thetaii=%f\n\n",xi,x,yi,y,zi,z, fPhi, phii, fTheta, thetai) ;
}
if ( wtot > 0 )