Track matching and PID parameters added to AliEMCALRecParam

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

diff --git a/EMCAL/macros/RecParamDB/AliEMCALSetRecParamCDB.C b/EMCAL/macros/RecParamDB/AliEMCALSetRecParamCDB.C
index 21e9c2b94cc1a10189ad1ad664f1d159bf17790b..868b89182584998d1906ca937e185ae032f8e723 100644 (file)
--- a/EMCAL/macros/RecParamDB/AliEMCALSetRecParamCDB.C
+++ b/EMCAL/macros/RecParamDB/AliEMCALSetRecParamCDB.C
@@ -54,10 +54,116 @@ void SetRecParam()
   // Create reconstruction parameter object and set parameter values
   
   AliEMCALRecParam *recParamDB = new AliEMCALRecParam();
+  //Clusterization
   recParamDB->SetClusteringThreshold(0.5);
   recParamDB->SetW0(4.5);
   recParamDB->SetMinECut(0.45);
 
+  //Track matching
+  recParamDB->SetTrkCutX(6.0);
+  recParamDB->SetTrkCutY(6.0);
+  recParamDB->SetTrkCutZ(6.0);
+  recParamDB->SetTrkCutR(10.0);
+  recParamDB->SetTrkCutAlphaMin(-50.0);
+  recParamDB->SetTrkCutAlphaMin( 50.0);
+  recParamDB->SetTrkCutAngle(10000.0);      // i.e. exclude this for the moment
+
+  //PID
+  recParamDB->SetGamma(0,0,  0.038022);
+  recParamDB->SetGamma(0,1, -0.0001883);
+  recParamDB->SetGamma(0,2,  5.449e-06);
+  
+  recParamDB->SetGamma(1,0,  0.207313);
+  recParamDB->SetGamma(1,1, -0.000978);
+  recParamDB->SetGamma(1,2,  0.00001634);
+  
+  recParamDB->SetGamma(2,0,  0.043364);
+  recParamDB->SetGamma(2,1, -0.0002048);
+  recParamDB->SetGamma(2,2,  8.661e-06);
+  recParamDB->SetGamma(2,3, -1.353e-07);
+  
+  recParamDB->SetGamma(3,0,  0.265004);
+  recParamDB->SetGamma(3,1,  0.061298);
+  recParamDB->SetGamma(3,2, -0.003203);
+  recParamDB->SetGamma(3,3,  4.73e-05);
+  
+  recParamDB->SetGamma(4,0,  0.243579);
+  recParamDB->SetGamma(4,1, -1.614e-05);
+  
+  recParamDB->SetGamma(5,0,  0.002942);
+  recParamDB->SetGamma(5,1, -3.976e-05);
+  
+  recParamDB->SetHadron(0,0,  0.011945 / 3.);
+  recParamDB->SetHadron(0,1,  0.000386 / 3.);
+  recParamDB->SetHadron(0,2, -0.000014 / 3.);
+  recParamDB->SetHadron(0,3,  1.336e-07 / 3.);
+  
+  recParamDB->SetHadron(1,0,  0.496544);
+  recParamDB->SetHadron(1,1, -0.003226);
+  recParamDB->SetHadron(1,2,  0.00001678);
+  
+  recParamDB->SetHadron(2,0,  0.144838);
+  recParamDB->SetHadron(2,1, -0.002954);
+  recParamDB->SetHadron(2,2,  0.00008754);
+  recParamDB->SetHadron(2,3, -7.587e-07);
+  
+  recParamDB->SetHadron(3,0,  1.264461 / 7.);
+  recParamDB->SetHadron(3,1,  0.002097 / 7.);
+  
+  recParamDB->SetHadron(4,0,  0.261950);
+  recParamDB->SetHadron(4,1, -0.001078);
+  recParamDB->SetHadron(4,2,  0.00003237);
+  recParamDB->SetHadron(4,3, -3.241e-07);
+  recParamDB->SetHadron(4,4,  0.);
+  recParamDB->SetHadron(4,5,  0.);
+  recParamDB->SetHadron(5,0,  0.010317);
+  recParamDB->SetHadron(5,1,  0.);
+  recParamDB->SetHadron(5,2,  0.);
+  recParamDB->SetHadron(5,3,  0.);
+  recParamDB->SetHadron(5,4,  0.);
+  recParamDB->SetHadron(5,5,  0.);
+  
+  recParamDB->SetPiZero5to10(0,0, 0.009138);
+  recParamDB->SetPiZero5to10(0,1, 0.0006377);
+  
+  recParamDB->SetPiZero5to10(1,0, 0.08);
+  
+  recParamDB->SetPiZero5to10(2,0, -0.061119);
+  recParamDB->SetPiZero5to10(2,1,  0.019013);
+  
+  recParamDB->SetPiZero5to10(3,0,  0.2);
+  
+  recParamDB->SetPiZero5to10(4,0,  0.252044);
+  recParamDB->SetPiZero5to10(4,1, -0.002315);
+  
+  recParamDB->SetPiZero5to10(5,0,  0.002942);
+  recParamDB->SetPiZero5to10(5,1, -3.976e-05);
+  
+  recParamDB->SetPiZero10to60(0,0,  0.009138);
+  recParamDB->SetPiZero10to60(0,1,  0.0006377);
+  
+  recParamDB->SetPiZero10to60(1,0,  1.272837);
+  recParamDB->SetPiZero10to60(1,1, -0.069708);
+  recParamDB->SetPiZero10to60(1,2,  0.001568);
+  recParamDB->SetPiZero10to60(1,3, -1.162e-05);
+  
+  recParamDB->SetPiZero10to60(2,0,  0.139703);
+  recParamDB->SetPiZero10to60(2,1,  0.003687);
+  recParamDB->SetPiZero10to60(2,2, -0.000568);
+  recParamDB->SetPiZero10to60(2,3,  1.498e-05);
+  recParamDB->SetPiZero10to60(2,4, -1.174e-07);
+  
+  recParamDB->SetPiZero10to60(3,0, -0.826367);
+  recParamDB->SetPiZero10to60(3,1,  0.096951);
+  recParamDB->SetPiZero10to60(3,2, -0.002215);
+  recParamDB->SetPiZero10to60(3,3,  2.523e-05);
+  
+  recParamDB->SetPiZero10to60(4,0,  0.249890);
+  recParamDB->SetPiZero10to60(4,1, -0.000063);
+  
+  recParamDB->SetPiZero10to60(5,0,  0.002942);
+  recParamDB->SetPiZero10to60(5,1, -3.976e-05);
+  
   // Store calibration data into database
   
   AliCDBMetaData md;
@@ -65,7 +171,7 @@ void SetRecParam()
   md.SetBeamPeriod(beamPeriod);
   md.SetResponsible("Yuri Kharlov");
   
-  AliCDBId id("EMCAL/RecParam/Data",firstRun,lastRun);
+  AliCDBId id("EMCAL/Config/RecParam",firstRun,lastRun);
 
   AliCDBManager* man = AliCDBManager::Instance();  
   AliCDBStorage* loc = man->GetStorage(DBFolder.Data());
@@ -86,9 +192,9 @@ void GetRecParam()
   AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT");
   AliCDBManager::Instance()->SetSpecificStorage("EMCAL/*",DBFolder.Data());
 
-  AliCDBEntry* cdbEntry = AliCDBManager::Instance()->Get("EMCAL/RecParam/Data",runNumber);
+  AliCDBEntry* cdbEntry = AliCDBManager::Instance()->Get("EMCAL/Config/RecParam/",runNumber);
   if (cdbEntry == 0) {
-    cerr << "No CDBEntry found at path "<<DBFolder.Data()<<"/"<<"EMCAL/RecParam/Data"<<endl;
+    cerr << "No CDBEntry found at path "<<DBFolder.Data()<<"/"<<"EMCAL/Config/RecParam/"<<endl;
     return;
   }