Possibility to set threshold in SZ runs

[u/mrichter/AliRoot.git] / PHOS / AliPHOSRecoParam.cxx

diff --git a/PHOS/AliPHOSRecoParam.cxx b/PHOS/AliPHOSRecoParam.cxx
index 1d6981cc6292ce0ea09278e8bbadff097dc38104..d12b9d677a5c2ef110de56ca9d6cfd763220ea2d 100644 (file)
--- a/PHOS/AliPHOSRecoParam.cxx
+++ b/PHOS/AliPHOSRecoParam.cxx
@@ -20,38 +20,63 @@
 // Author: Boris Polichtchouk.
 
 // --- AliRoot header files ---
 // Author: Boris Polichtchouk.
 
 // --- AliRoot header files ---

+#include "AliCDBManager.h"
+#include "AliCDBEntry.h"
+#include "AliLog.h"

 #include "AliPHOSRecoParam.h"
 
 ClassImp(AliPHOSRecoParam)
 
 #include "AliPHOSRecoParam.h"
 
 ClassImp(AliPHOSRecoParam)
 

+TObjArray* AliPHOSRecoParam::fgkMaps =0; //ALTRO mappings
+

 //-----------------------------------------------------------------------------
 AliPHOSRecoParam::AliPHOSRecoParam() :
 //-----------------------------------------------------------------------------
 AliPHOSRecoParam::AliPHOSRecoParam() :

-  TNamed(),
-  fClusteringThreshold(9999),
-  fLocMaxCut(9999),
-  fMinE(9999),
-  fW0(9999),
-  fSampleQualityCut(1.),
-  fEcoreRadius(3.),
-  fSubtractPedestals(kTRUE),
-  fUnfold(kTRUE),
-  fDecoderVersion("")
+  AliDetectorRecoParam(),
+  fEMCClusteringThreshold(0.2),
+  fEMCLocMaxCut(0.03),
+  fEMCRawDigitThreshold(2),
+  fEMCMinE(0.012),
+  fEMCW0(4.5),
+  fEMCSampleQualityCut(1.),
+  fEMCEcoreRadius(3.),
+  fEMCEcore2ESD(kFALSE),
+  fEMCSubtractPedestals(kTRUE),
+  fEMCUnfold(kTRUE),
+  fEMCEnergyCorrectionOn(kTRUE),
+  fEMCDecoderVersion(""),
+  fGlobalAltroOffset(0),
+  fGlobalAltroThreshold(0),
+  fCPVClusteringThreshold(0.0),
+  fCPVLocMaxCut(0.03),
+  fCPVMinE(0.0),
+  fCPVW0(4.0),
+  fCPVUnfold(kTRUE)

 {
   //Default constructor.
 }
 
 //-----------------------------------------------------------------------------
 AliPHOSRecoParam::AliPHOSRecoParam(const AliPHOSRecoParam& ):
 {
   //Default constructor.
 }
 
 //-----------------------------------------------------------------------------
 AliPHOSRecoParam::AliPHOSRecoParam(const AliPHOSRecoParam& ):

-  TNamed(),
-  fClusteringThreshold(9999),
-  fLocMaxCut(9999),
-  fMinE(9999),
-  fW0(9999),
-  fSampleQualityCut(1.),
-  fEcoreRadius(3.),
-  fSubtractPedestals(kTRUE),
-  fUnfold(kTRUE),
-  fDecoderVersion("")
+  AliDetectorRecoParam(),
+  fEMCClusteringThreshold(0.2),
+  fEMCLocMaxCut(0.03),
+  fEMCRawDigitThreshold(2),
+  fEMCMinE(0.012),
+  fEMCW0(4.5),
+  fEMCSampleQualityCut(1.),
+  fEMCEcoreRadius(3.),
+  fEMCEcore2ESD(kFALSE),
+  fEMCSubtractPedestals(kTRUE),
+  fEMCUnfold(kTRUE),
+  fEMCEnergyCorrectionOn(kTRUE),
+  fEMCDecoderVersion(""),
+  fGlobalAltroOffset(0),
+  fGlobalAltroThreshold(0),
+  fCPVClusteringThreshold(0.0),
+  fCPVLocMaxCut(0.03),
+  fCPVMinE(0.0),
+  fCPVW0(4.0),
+  fCPVUnfold(kTRUE)

 {
   //Copy constructor.
 }
 {
   //Copy constructor.
 }


@@ -62,16 +87,92 @@ AliPHOSRecoParam& AliPHOSRecoParam::operator = (const AliPHOSRecoParam& recoPara
   //Assignment operator.
 
   if(this != &recoParam) {
   //Assignment operator.
 
   if(this != &recoParam) {

-    fClusteringThreshold = recoParam.fClusteringThreshold;
-    fLocMaxCut           = recoParam.fLocMaxCut;
-    fMinE                = recoParam.fMinE;
-    fW0                  = recoParam.fW0;
-    fSampleQualityCut    = recoParam.fSampleQualityCut ;
-    fEcoreRadius         = recoParam.fEcoreRadius ;
-    fSubtractPedestals   = recoParam.fSubtractPedestals;
-    fUnfold              = recoParam.fUnfold;
-    fDecoderVersion      = recoParam.fDecoderVersion ;
+    fEMCClusteringThreshold = recoParam.fEMCClusteringThreshold;
+    fEMCLocMaxCut           = recoParam.fEMCLocMaxCut;
+    fEMCRawDigitThreshold   = recoParam.fEMCRawDigitThreshold;
+    fEMCMinE                = recoParam.fEMCMinE;
+    fEMCW0                  = recoParam.fEMCW0;
+    fEMCSampleQualityCut    = recoParam.fEMCSampleQualityCut;
+    fEMCEcoreRadius         = recoParam.fEMCEcoreRadius;
+    fEMCEcore2ESD           = recoParam.fEMCEcore2ESD;
+    fEMCSubtractPedestals   = recoParam.fEMCSubtractPedestals;
+    fEMCUnfold              = recoParam.fEMCUnfold;
+    fEMCEnergyCorrectionOn  = recoParam.fEMCEnergyCorrectionOn;
+    fEMCDecoderVersion      = recoParam.fEMCDecoderVersion;
+    fGlobalAltroOffset      = recoParam.fGlobalAltroOffset;
+    fGlobalAltroThreshold   = recoParam.fGlobalAltroThreshold;
+    fCPVClusteringThreshold = recoParam.fCPVClusteringThreshold;
+    fCPVLocMaxCut           = recoParam.fCPVLocMaxCut;
+    fCPVMinE                = recoParam.fCPVMinE;
+    fCPVW0                  = recoParam.fCPVW0;
+    fCPVUnfold              = recoParam.fCPVUnfold;

   }
 
   return *this;
 }
   }
 
   return *this;
 }

+
+//-----------------------------------------------------------------------------
+void AliPHOSRecoParam::Print(Option_t * /*option*/) const
+{
+  AliDebug(2,Form("PHOS reconstruction parameters:\n"
+                 "\tEMCClusteringThreshold = %f\n"
+                 "\tEMCLocMaxCut           = %f\n"
+                 "\tEMCRawDigitThreshold   = %f\n"
+                 "\tEMCMinE                = %f\n"
+                 "\tEMCW0                  = %f\n"
+                 "\tEMCSampleQualityCut    = %f\n"
+                 "\tEMCEcoreRadius         = %f\n"
+                 "\tEMCEcore2ESD           = %d\n"
+                 "\tEMCSubtractPedestals   = %d\n"
+                 "\tEMCUnfold              = %d\n"
+                 "\tEMCEnergyCorrectionOn  = %d\n"
+                 "\tEMCDecoderVersion      = \"%s\"\n"
+                 "\tGlobalAltroOffset      = %d",
+                 "\tGlobalAltroThreshold   = %d",
+                 fEMCClusteringThreshold,
+                 fEMCLocMaxCut,
+                 fEMCRawDigitThreshold,
+                 fEMCMinE,
+                 fEMCW0,
+                 fEMCSampleQualityCut,
+                 fEMCEcoreRadius,
+                 fEMCEcore2ESD,
+                 fEMCSubtractPedestals,
+                 fEMCUnfold,
+                 fEMCEnergyCorrectionOn,
+                 fEMCDecoderVersion.Data(),
+                 fGlobalAltroOffset,
+                 fGlobalAltroThreshold));
+
+}
+
+//-----------------------------------------------------------------------------
+AliPHOSRecoParam* AliPHOSRecoParam::GetDefaultParameters()
+{
+  //Default parameters for the reconstruction
+
+  AliPHOSRecoParam* params = new AliPHOSRecoParam();
+  return params;
+}
+
+//-----------------------------------------------------------------------------
+const TObjArray* AliPHOSRecoParam::GetMappings()
+{
+  //Returns array of AliAltroMappings for RCU0..RCU3.
+  //If not found, read it from OCDB.
+
+  //Quick check as follows:
+  //  root [0] AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT");
+  //  root [1] AliCDBManager::Instance()->SetRun(1);
+  //  root [2] TObjArray* maps = AliPHOSRecoParam::GetMappings();
+  //  root [3] maps->Print();
+  
+  if(fgkMaps) return fgkMaps;
+  
+  AliCDBEntry* entry = AliCDBManager::Instance()->Get("PHOS/Calib/Mapping");
+  if(entry)
+    fgkMaps = (TObjArray*)entry->GetObject();
+  
+  return fgkMaps;
+  
+}