minor fixes for cluster unfolding
authorjklay <jklay@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 27 May 2008 19:35:59 +0000 (19:35 +0000)
committerjklay <jklay@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 27 May 2008 19:35:59 +0000 (19:35 +0000)
EMCAL/AliEMCALClusterizerv1.cxx
EMCAL/AliEMCALRecParam.h
EMCAL/AliEMCALRecPoint.cxx

index 19ff314..a7e4fff 100644 (file)
@@ -270,7 +270,7 @@ Bool_t AliEMCALClusterizerv1::FindFit(AliEMCALRecPoint * RecPoint, AliEMCALDigit
   Double_t p2 = 0.0 ;
 
   gMinuit->mnexcm("SET STR", &p2, 0, ierflg) ;   // force TMinuit to reduce function calls
-  gMinuit->mnexcm("SET GRA", &p1, 1, ierflg) ;   // force TMinuit to use my gradient
+  //  gMinuit->mnexcm("SET GRA", &p1, 1, ierflg) ;   // force TMinuit to use my gradient
   gMinuit->SetMaxIterations(5);
   gMinuit->mnexcm("SET NOW", &p2 , 0, ierflg) ;  // No Warnings
   gMinuit->mnexcm("MIGRAD", &p0, 0, ierflg) ;    // minimize
@@ -355,6 +355,7 @@ void AliEMCALClusterizerv1::InitParameters()
     fECAW0                  = recParam->GetW0();
     fMinECut                = recParam->GetMinECut();    
     fToUnfold               = recParam->GetUnfold();
+    if(fToUnfold) AliWarning("Cluster Unfolding ON. Implementing only for eta=0 case!!!"); 
     fECALocMaxCut           = recParam->GetLocMaxCut();
 
     AliDebug(1,Form("Reconstruction parameters: fECAClusteringThreshold=%.3f, fECAW=%.3f, fMinECut=%.3f, fToUnfold=%d, fECALocMaxCut=%.3f",
index 2fb0b45..1e2b4dd 100644 (file)
@@ -37,7 +37,7 @@ public:
   void SetW0                 (Float_t w0)      {fW0 = w0                    ;}
   void SetMinECut            (Float_t minEcut) {fMinECut = minEcut          ;}
   void SetLocMaxCut          (Float_t locMaxCut) {fLocMaxCut = locMaxCut    ;}
-  void SetUnfold             (Bool_t unfold)     {fUnfold = unfold          ; if(fUnfold) AliInfo("Cluster Unfolding ON. Implementing only for eta=0 case!!!");}
+  void SetUnfold             (Bool_t unfold)     {fUnfold = unfold          ; if(fUnfold) AliWarning("Cluster Unfolding ON. Implementing only for eta=0 case!!!");}
 
   //PID (Guenole)
   Double_t GetGamma(Int_t i, Int_t j) const    {return fGamma[i][j];} 
index 398194c..b4d94fc 100644 (file)
@@ -55,7 +55,7 @@ AliEMCALRecPoint::AliEMCALRecPoint()
     fMulTrack(0), fDigitsList(0), fTracksList(0),
     fClusterType(-1), fCoreEnergy(0), fDispersion(0),
     fEnergyList(0), fTimeList(0), fAbsIdList(0),
-    fTime(0.), fCoreRadius(10),  //HG check this
+    fTime(0.), fNExMax(0), fCoreRadius(10),  //HG check this 
     fDETracksList(0), fMulParent(0), fMaxParent(0),
     fParentsList(0), fDEParentsList(0), fSuperModuleNumber(0),
     fDigitIndMax(-1)
@@ -77,7 +77,7 @@ AliEMCALRecPoint::AliEMCALRecPoint(const char *)
     fDigitsList(new Int_t[fMaxDigit]), fTracksList(new Int_t[fMaxTrack]),
     fClusterType(-1), fCoreEnergy(0), fDispersion(0),
     fEnergyList(new Float_t[fMaxDigit]), fTimeList(new Float_t[fMaxDigit]), 
-    fAbsIdList(new Int_t[fMaxDigit]), fTime(-1.), fCoreRadius(10),
+    fAbsIdList(new Int_t[fMaxDigit]), fTime(-1.), fNExMax(0), fCoreRadius(10),
     fDETracksList(new Float_t[fMaxTrack]), fMulParent(0), fMaxParent(1000),
     fParentsList(new Int_t[fMaxParent]), fDEParentsList(new Float_t[fMaxParent]),
     fSuperModuleNumber(0), fDigitIndMax(-1)
@@ -106,7 +106,7 @@ AliEMCALRecPoint::AliEMCALRecPoint(const AliEMCALRecPoint & rp)
     fClusterType(rp.fClusterType), fCoreEnergy(rp.fCoreEnergy), 
     fDispersion(rp.fDispersion),
     fEnergyList(new Float_t[rp.fMaxDigit]), fTimeList(new Float_t[rp.fMaxDigit]), 
-    fAbsIdList(new Int_t[rp.fMaxDigit]), fTime(rp.fTime), fCoreRadius(rp.fCoreRadius),
+    fAbsIdList(new Int_t[rp.fMaxDigit]), fTime(rp.fTime), fNExMax(rp.fNExMax),fCoreRadius(rp.fCoreRadius),
     fDETracksList(new Float_t[rp.fMaxTrack]), fMulParent(rp.fMulParent), 
     fMaxParent(rp.fMaxParent), fParentsList(new Int_t[rp.fMaxParent]), 
     fDEParentsList(new Float_t[rp.fMaxParent]),
@@ -179,6 +179,7 @@ AliEMCALRecPoint& AliEMCALRecPoint::operator= (const AliEMCALRecPoint &rp)
     fAbsIdList[i] = rp.fAbsIdList[i];
   }
   fTime = rp.fTime;
+  fNExMax = rp.fNExMax;
   fCoreRadius = rp.fCoreRadius;
   for(Int_t i = 0; i < fMaxTrack; i++) fDETracksList[i] = rp.fDETracksList[i];
   fMulParent = rp.fMulParent;