Change reco param default and OCDB settings, PbPb clusterizer is v2 and minimum cell...

[u/mrichter/AliRoot.git] / EMCAL / macros / RecParamDB / AliEMCALSetRecParamCDB.C

diff --git a/EMCAL/macros/RecParamDB/AliEMCALSetRecParamCDB.C b/EMCAL/macros/RecParamDB/AliEMCALSetRecParamCDB.C
index 95a2c69394f700490805d730e9047a9d97ba8bcb..91801ae97735c2c4e63af2221498ad093952f046 100644 (file)
--- a/EMCAL/macros/RecParamDB/AliEMCALSetRecParamCDB.C
+++ b/EMCAL/macros/RecParamDB/AliEMCALSetRecParamCDB.C
@@ -16,7 +16,7 @@
 #endif
 
 
-void AliEMCALSetRecParamCDB(AliRecoParam::EventSpecie_t default = AliRecoParam::kDefault)
+void AliEMCALSetRecParamCDB(AliRecoParam::EventSpecie_t default = AliRecoParam::kLowMult)
 {
   
   // Create an object AliEMCALRecParam and store it to OCDB
@@ -37,6 +37,7 @@ void AliEMCALSetRecParamCDB(AliRecoParam::EventSpecie_t default = AliRecoParam::
     //Add to the recParamArray
     recParamDB->SetEventSpecie(AliRecoParam::kDefault);
     recParamArray->AddLast(recParamDB);
+         
   }
   
   //Add other options here, if desired, for
@@ -126,24 +127,29 @@ AliEMCALRecParam* GetHighMultiplicityParameters()
   //AliEMCALRecParam *recParamDB = AliEMCALRecParam::GetDefaultParameters();
   
   AliEMCALRecParam* params =  AliEMCALRecParam::GetDefaultParameters();
-  
   //Clusterization
-  params->SetClusteringThreshold(0.5);
-  params->SetW0(4.5);
-  params->SetMinECut(0.45);
+
+  // params->SetClusterizerFlag(AliEMCALRecParam::kClusterizerv1);
+  // params->SetClusteringThreshold(0.5);
+  // params->SetMinECut(0.45);
+
+  params->SetClusterizerFlag(AliEMCALRecParam::kClusterizerv2);
+  params->SetClusteringThreshold(0.1); // 100 MeV                                             
+  params->SetMinECut(0.05);  //50 MeV    
+   
   params->SetUnfold(kFALSE);
-  params->SetLocMaxCut(0.03);
-  
+
+  params->SetW0(4.5);
+
+  params->SetTimeCut(250e-9);//250 ns
+  params->SetTimeMin(425e-9);//425 ns
+  params->SetTimeMax(825e-9);//825 ns
+
   //Track matching
-  params->SetTrkCutX(6.0);
-  params->SetTrkCutY(6.0);
-  params->SetTrkCutZ(6.0);
-  params->SetTrkCutR(10.0);
-  params->SetTrkCutAlphaMin(-50.0);
-  params->SetTrkCutAlphaMax( 50.0);
-  params->SetTrkCutNITS(3.0);
+
+  params->SetTrkCutNITS(1.0);
   params->SetTrkCutNTPC(20.0);
-  params->SetTrkCutAngle(10000.0);      // i.e. exclude this for the moment
+  params->SetExtrapolateStep(20.);
   
   //PID
        
@@ -286,7 +292,11 @@ AliEMCALRecParam* GetHighMultiplicityParameters()
   params->SetTau(2.35);
   params->SetNoiseThreshold(3);
   params->SetNPedSamples(5);
-  
+  params->SetRemoveBadChannels(kFALSE);
+  params->SetFittingAlgorithm(0);//(AliCaloConstants::kStandard);  
+  params->SetFALTROUsage(kTRUE); 
+  params->SetLEDFit(kFALSE);   
+
   return params ;
 }      
 
@@ -298,8 +308,24 @@ AliEMCALRecParam* GetLowMultiplicityParameters()
   //AliEMCALRecParam *recParamDB = AliEMCALRecParam::GetLowFluxParam();
   
   AliEMCALRecParam* params =  AliEMCALRecParam::GetDefaultParameters();
+  //params->SetClusterizerFlag(AliEMCALRecParam::kClusterizerNxN);
+  params->SetClusterizerFlag(AliEMCALRecParam::kClusterizerv1);
   params->SetClusteringThreshold(0.1); // 100 MeV                                             
-  params->SetMinECut(0.01);  //10 MeV          
+  params->SetMinECut(0.05);  //50 MeV          
+
+  params->SetUnfold(kFALSE);
+
+  params->SetW0(4.5);
+
+  params->SetTimeCut(250e-9);//250 ns
+  params->SetTimeMin(425e-9);//425 ns
+  params->SetTimeMax(825e-9);//825 ns
+
+  // Track Matching
+  
+  params->SetTrkCutNITS(1.0);
+  params->SetTrkCutNTPC(20.0);
+  params->SetExtrapolateStep(20.);  
   
   //PID parameters for pp  implemented 
   // as a first step, all array elements are initialized to 0.0
@@ -437,7 +463,18 @@ AliEMCALRecParam* GetLowMultiplicityParameters()
   params->SetHadronEnergyProb(2,  2.956727e-01);
   params->SetHadronEnergyProb(3, -3.051022e+01);
   params->SetHadronEnergyProb(4, -6.036931e-02);
-  
+
+  // raw signal fitting
+  params->SetHighLowGainFactor(16.);
+  params->SetOrderParameter(2);
+  params->SetTau(2.35);
+  params->SetNoiseThreshold(3);
+  params->SetNPedSamples(5);
+  params->SetRemoveBadChannels(kFALSE);
+  params->SetFittingAlgorithm(0);//(AliCaloConstants::kStandard);  
+  params->SetFALTROUsage(kTRUE); 
+  params->SetLEDFit(kFALSE);   
+
   return params;
   
 }