Adjusting reco-parameters for warm PHOS to reduce ESD/AOD size
authorkharlov <kharlov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 10 Sep 2008 17:57:55 +0000 (17:57 +0000)
committerkharlov <kharlov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 10 Sep 2008 17:57:55 +0000 (17:57 +0000)
PHOS/AliPHOSClusterizerv1.cxx
PHOS/AliPHOSRecoParam.cxx
PHOS/AliPHOSRecoParam.h
PHOS/AliPHOSReconstructor.cxx
PHOS/Calib/RecoParam/Run0_999999999_v0_s0.root
PHOS/MakePHOSRecoParam.C

index 1ed11b2..5f5d221 100644 (file)
@@ -384,8 +384,10 @@ void AliPHOSClusterizerv1::InitParameters()
   const AliPHOSRecoParam* recoParam = AliPHOSReconstructor::GetRecoParam();
   if(!recoParam) AliFatal("Reconstruction parameters are not set!");
 
-  fCpvClusteringThreshold  = recoParam->GetEMCClusteringThreshold();
-  fEmcClusteringThreshold  = recoParam->GetCPVClusteringThreshold();
+  recoParam->Print();
+
+  fEmcClusteringThreshold  = recoParam->GetEMCClusteringThreshold();
+  fCpvClusteringThreshold  = recoParam->GetCPVClusteringThreshold();
   
   fEmcLocMaxCut            = recoParam->GetEMCLocalMaxCut();
   fCpvLocMaxCut            = recoParam->GetCPVLocalMaxCut();
index 9b5aed0..592546f 100644 (file)
@@ -22,6 +22,7 @@
 // --- AliRoot header files ---
 #include "AliCDBManager.h"
 #include "AliCDBEntry.h"
+#include "AliLog.h"
 #include "AliPHOSRecoParam.h"
 
 ClassImp(AliPHOSRecoParam)
@@ -41,6 +42,7 @@ AliPHOSRecoParam::AliPHOSRecoParam() :
   fEMCSubtractPedestals(kTRUE),
   fEMCUnfold(kTRUE),
   fEMCDecoderVersion(""),
+  fGlobalAltroOffset(0),
   fCPVClusteringThreshold(0.0),
   fCPVLocMaxCut(0.03),
   fCPVMinE(0.0),
@@ -63,6 +65,7 @@ AliPHOSRecoParam::AliPHOSRecoParam(const AliPHOSRecoParam& ):
   fEMCSubtractPedestals(kTRUE),
   fEMCUnfold(kTRUE),
   fEMCDecoderVersion(""),
+  fGlobalAltroOffset(0),
   fCPVClusteringThreshold(0.0),
   fCPVLocMaxCut(0.03),
   fCPVMinE(0.0),
@@ -88,6 +91,7 @@ AliPHOSRecoParam& AliPHOSRecoParam::operator = (const AliPHOSRecoParam& recoPara
     fEMCSubtractPedestals   = recoParam.fEMCSubtractPedestals;
     fEMCUnfold              = recoParam.fEMCUnfold;
     fEMCDecoderVersion      = recoParam.fEMCDecoderVersion;
+    fGlobalAltroOffset      = recoParam.fGlobalAltroOffset;
     fCPVClusteringThreshold = recoParam.fCPVClusteringThreshold;
     fCPVLocMaxCut           = recoParam.fCPVLocMaxCut;
     fCPVMinE                = recoParam.fCPVMinE;
@@ -99,6 +103,35 @@ AliPHOSRecoParam& AliPHOSRecoParam::operator = (const AliPHOSRecoParam& recoPara
 }
 
 //-----------------------------------------------------------------------------
+void AliPHOSRecoParam::Print(Option_t *option) const
+{
+  AliDebug(2,Form("PHOS reconstruction parameters:\n"
+                 "\tEMCClusteringThreshold = %f\n"
+                 "\tEMCLocMaxCut           = %f\n"
+                 "\tEMCMinE                = %f\n"
+                 "\tEMCW0                  = %f\n"
+                 "\tEMCSampleQualityCut    = %f\n"
+                 "\tEMCEcoreRadius         = %f\n"
+                 "\tEMCEcore2ESD           = %d\n"
+                 "\tEMCSubtractPedestals   = %d\n"
+                 "\tEMCUnfold              = %d\n"
+                 "\tEMCDecoderVersion      = \"%s\"\n"
+                 "\tGlobalAltroOffset      = %d",
+                 fEMCClusteringThreshold,
+                 fEMCLocMaxCut,
+                 fEMCMinE,
+                 fEMCW0,
+                 fEMCSampleQualityCut,
+                 fEMCEcoreRadius,
+                 fEMCEcore2ESD,
+                 fEMCSubtractPedestals,
+                 fEMCUnfold,
+                 fEMCDecoderVersion.Data(),
+                 fGlobalAltroOffset));
+
+}
+
+//-----------------------------------------------------------------------------
 AliPHOSRecoParam* AliPHOSRecoParam::GetDefaultParameters()
 {
   //Default parameters for the reconstruction
index 7dd8781..55e373b 100644 (file)
@@ -55,6 +55,8 @@ public:
   void SetCPVLogWeight(Float_t w)                    { fCPVW0                 =w;       }
   void SetCPVUnfolding(Bool_t toUnfold=kFALSE)       { fCPVUnfold            =toUnfold;}
 
+  virtual void Print(const Option_t *option="") const;
+
   static AliPHOSRecoParam* GetDefaultParameters();
   static const  TObjArray* GetMappings();
 
index c1fcddf..e601566 100644 (file)
@@ -223,7 +223,7 @@ void AliPHOSReconstructor::FillESD(TTree* digitsTree, TTree* clustersTree,
   // Add to CaloCells only EMC digits with non-zero energy 
   for (Int_t idig = 0 ; idig < nDigits ; idig++) {
     const AliPHOSDigit * dig = (const AliPHOSDigit*)fgDigitsArray->At(idig);
-    if(dig->GetId() <= knEMC && dig->GetEnergy() > 0 ){
+    if(dig->GetId() <= knEMC && dig->GetEnergy() > GetRecoParam()->GetEMCMinE() ){
       //printf("i %d; id %d; amp %f; time %e\n",
       //idignew,dig->GetId(),dig->GetEnergy(), dig->GetTime());
       phsCells.SetCell(idignew,dig->GetId(), dig->GetEnergy(), dig->GetTime());   
@@ -334,7 +334,6 @@ void  AliPHOSReconstructor::ConvertDigits(AliRawReader* rawReader, TTree* digits
   // Converts raw data to
   // PHOS digits
   // Works on a single-event basis
-
   rawReader->Reset() ; 
 
   AliPHOSRawDecoder * dc ;
index aac6195..a4154af 100644 (file)
Binary files a/PHOS/Calib/RecoParam/Run0_999999999_v0_s0.root and b/PHOS/Calib/RecoParam/Run0_999999999_v0_s0.root differ
index 1a554a2..e160bfc 100644 (file)
@@ -19,17 +19,18 @@ void MakePHOSRecoParam(AliRecoParam::EventSpecie_t default=AliRecoParam::kDefaul
 
   {
     AliPHOSRecoParam * phosRecoParam = AliPHOSRecoParam::GetDefaultParameters();
-    phosRecoParam->SetEventSpecie(AliRecoParam::kCosmic);
-    phosRecoParam->SetEMCSubtractPedestals(kTRUE);
-    phosRecoParam->SetEMCMinE(0.01);
-    phosRecoParam->SetEMCClusteringThreshold(0.02);
+    phosRecoParam->SetEventSpecie(AliRecoParam::kDefault);
+    phosRecoParam->SetEMCSubtractPedestals(kFALSE);
+    phosRecoParam->SetEMCMinE(0.045);
+    phosRecoParam->SetEMCClusteringThreshold(0.20);
     phosRecoParam->SetEMCDecoderVersion("v1");
     recoParamArray->AddLast(phosRecoParam);
   }
   {
     AliPHOSRecoParam * phosRecoParam = AliPHOSRecoParam::GetDefaultParameters();
     recoParamArray->AddLast(phosRecoParam);
-    phosRecoParam->SetEventSpecie(AliRecoParam::kDefault);
+    phosRecoParam->SetEMCClusteringThreshold(0.02008);
+    phosRecoParam->SetEventSpecie(AliRecoParam::kCosmic);
   }
 
   // Set the default