Move option to reclusterizer with calibrated amplitudes from NxN to mother Clusterize...
authorgconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 26 Oct 2010 08:06:34 +0000 (08:06 +0000)
committergconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 26 Oct 2010 08:06:34 +0000 (08:06 +0000)
EMCAL/AliEMCALClusterizer.cxx
EMCAL/AliEMCALClusterizer.h
EMCAL/AliEMCALClusterizerNxN.cxx
EMCAL/AliEMCALClusterizerNxN.h

index ad7a237..f3cae20 100644 (file)
@@ -60,6 +60,9 @@ class AliCDBStorage;
 
 ClassImp(AliEMCALClusterizer)
 
+Bool_t AliEMCALClusterizer::fgkIsInputCalibrated = kFALSE;
+
+
 //____________________________________________________________________________
 AliEMCALClusterizer::AliEMCALClusterizer():
   fDigitsArr(NULL),
@@ -210,15 +213,30 @@ Float_t  AliEMCALClusterizer::Calibrate(const Float_t amp, const Float_t time, c
     }
     //Check if time is too large or too small, indication of a noisy channel, remove in this case
     if(time > fTimeMax || time < fTimeMin) return 0;
+    
+    if (AliEMCALClusterizer::fgkIsInputCalibrated == kTRUE)
+    {
+      AliDebug(10, Form("Input already calibrated!"));
+      return amp;
+    }
          
     fADCchannelECA  = fCalibData->GetADCchannel (iSupMod,ieta,iphi);
     fADCpedestalECA = fCalibData->GetADCpedestal(iSupMod,ieta,iphi);
     
+    
     return -fADCpedestalECA + amp * fADCchannelECA ;        
     
   }
-  else //Return energy with default parameters if calibration is not available
+  else{ //Return energy with default parameters if calibration is not available
+    
+    if (AliEMCALClusterizer::fgkIsInputCalibrated == kTRUE)
+    {
+      AliDebug(10, Form("Input already calibrated!"));
+      return amp;
+    }    
+    
     return -fADCpedestalECA + amp * fADCchannelECA ; 
+  }
   
 }
 
@@ -480,5 +498,14 @@ void AliEMCALClusterizer::SetOutput(TTree *clustersTree)
   fTreeR->Branch("EMCALECARP", "TObjArray", &fRecPoints, bufsize, split);
 }
 
+//___________________________________________________________________
+void   AliEMCALClusterizer::SetInputCalibrated(Bool_t val)
+{
+  //
+  // input is calibrated - the case when we run already on ESD
+  //
+  AliEMCALClusterizer::fgkIsInputCalibrated = val;
+}
+
 
 
index 438cd9d..4aa0956 100644 (file)
@@ -70,14 +70,16 @@ public:
   virtual void PrintRecPoints(Option_t * option);
   virtual void PrintRecoInfo();                        //*MENU*
 
-
+  static void     SetInputCalibrated(Bool_t val);
   
   virtual const char * Version() const {Warning("Version", "Not Defined") ; return 0 ; } 
 
 protected:
 
   virtual void MakeClusters() = 0;
-
+  
+  static Bool_t fgkIsInputCalibrated; // to enable reclusterization from ESD cells
+  
   TClonesArray *fDigitsArr; // Array with EMCAL digits
   TTree *fTreeR;            // Tree with output clusters
   TObjArray    *fRecPoints; // Array with EMCAL clusters
@@ -113,7 +115,7 @@ private:
   
   
   
-  ClassDef(AliEMCALClusterizer,3)  // Clusterization algorithm class 
+  ClassDef(AliEMCALClusterizer,4)  // Clusterization algorithm class 
 } ;
 
 #endif // AliEMCALCLUSTERIZER_H
index 8c6de2d..a4eadaa 100644 (file)
@@ -65,8 +65,6 @@
 
 ClassImp(AliEMCALClusterizerNxN)
 
-Bool_t AliEMCALClusterizerNxN::fgkIsInputCalibrated = kFALSE;
-
 //____________________________________________________________________________
 AliEMCALClusterizerNxN::AliEMCALClusterizerNxN()
   : AliEMCALClusterizer()
@@ -360,11 +358,3 @@ void AliEMCALClusterizerNxN::MakeClusters()
   AliDebug(1,Form("total no of clusters %d from %d digits",fNumberOfECAClusters,fDigitsArr->GetEntriesFast())); 
 }
 
-//___________________________________________________________________
-void   AliEMCALClusterizerNxN::SetInputCalibrated(Bool_t val)
-{
-  //
-  // input is calibrated - the case when we run already on ESD
-  //
-  AliEMCALClusterizerNxN::fgkIsInputCalibrated = val;
-}
index 4e5457e..007844d 100644 (file)
@@ -35,21 +35,17 @@ public:
 
   virtual void   Digits2Clusters(Option_t *option);                // Does the job
 
-  static void     SetInputCalibrated(Bool_t val);
-
   virtual const char * Version() const { return "clu-NxN" ; }  
   
 protected:
 
   virtual void   MakeClusters();            
 
-  static Bool_t fgkIsInputCalibrated; // to enable reclusterization from ESD cells
-
 private:
   AliEMCALClusterizerNxN(const AliEMCALClusterizerNxN &); //copy ctor
   AliEMCALClusterizerNxN & operator = (const AliEMCALClusterizerNxN &);
 
-   ClassDef(AliEMCALClusterizerNxN,1)   // Clusterizer implementation version 1
+   ClassDef(AliEMCALClusterizerNxN,2)   // Clusterizer implementation version 1
 
 };