Calibration macros directory

diff --git a/EMCAL/macros/CalibrationDB/AliEMCALDecalibrate.C b/EMCAL/macros/CalibrationDB/AliEMCALDecalibrate.C
new file mode 100644 (file)
index 0000000..eeb6017
--- /dev/null
+++ b/EMCAL/macros/CalibrationDB/AliEMCALDecalibrate.C
@@ -0,0 +1,41 @@
+
+// Digitize and decalibrate events assuming that SDigits 
+// have been already produced.
+// Decalibration coefficients are located in the local file
+// DeCalibDB/EMCAL/Calib/GainFactors_and_Pedestals/Run_xxx.root
+// Author: Boris Polichtchouk (Boris.Polichtchouk@cern.ch)
+// Modified from PHOS script for EMCAL by Gustavo Conesa
+
+void AliEMCALDecalibrate(Int_t nevents=3)
+{
+
+  //Load (de)calibration database into aliroot session
+  //and set it to AliEMCALGetter.
+  
+  AliEMCALCalibData* deCal  = (AliEMCALCalibData*)(AliCDBManager::Instance()
+    ->GetStorage("local://DeCalibDB")->Get("EMCAL/Calib/GainFactors_and_Pedestals",1)
+    ->GetObject());
+  
+  //Loader  
+  AliRunLoader* rl=0x0;
+  
+  rl = AliRunLoader::Open("galice.root",
+                         AliConfig::GetDefaultEventFolderName(),
+                         "update");
+  if (rl == 0x0)
+    {
+      gAlice->Fatal("Config.C","Can not instatiate the Run Loader");
+      return;
+    }
+
+  AliEMCALLoader *emcalLoader = dynamic_cast<AliEMCALLoader*>
+    (rl->GetDetectorLoader("EMCAL"));
+  //  rl->LoadSDigits("EMCAL");
+  emcalLoader->SetCalibData(deCal);
+
+  AliSimulation sim ; 
+  sim.SetRunGeneration(kFALSE) ;
+  sim.SetMakeSDigits("") ;
+  sim.SetMakeDigits("EMCAL") ;
+  sim.Run(nevents) ;  
+}

diff --git a/EMCAL/macros/CalibrationDB/AliEMCALLocalReconstruction.C b/EMCAL/macros/CalibrationDB/AliEMCALLocalReconstruction.C
new file mode 100644 (file)
index 0000000..5ca2477
--- /dev/null
+++ b/EMCAL/macros/CalibrationDB/AliEMCALLocalReconstruction.C
@@ -0,0 +1,17 @@
+
+
+//Run EMCAL clusterization using information from calibration database.
+// Author: Gustavo Conesa
+
+AliEMCALLocalReconstruction()
+{
+
+  // Open local calibration data base
+  AliCDBStorage* storage = AliCDBManager::Instance()->GetStorage("local://CalibDB");
+  
+  // Make clusterization, calibration parameters will be taken from CDB
+
+  AliReconstruction rec ; 
+  rec.SetRunLocalReconstruction("EMCAL") ;
+  rec.Run() ;
+}

diff --git a/EMCAL/macros/CalibrationDB/AliEMCALSetCDB.C b/EMCAL/macros/CalibrationDB/AliEMCALSetCDB.C
new file mode 100644 (file)
index 0000000..8bc8f96
--- /dev/null
+++ b/EMCAL/macros/CalibrationDB/AliEMCALSetCDB.C
@@ -0,0 +1,194 @@
+
+// Script to create calibration parameters and store them into CDB
+// Two sets of calibration parameters can be created:
+// 1) equal parameters
+// 2) randomly distributed parameters for decalibrated detector silumations
+// Modified from PHOS script for EMCAL by Gustavo Conesa
+
+#if !defined(__CINT__)
+#include "TControlBar.h"
+#include "TString.h"
+#include "TRandom.h"
+#include "TH2F.h"
+#include "TCanvas.h"
+
+#include "AliRun.h"
+#include "AliEMCALCalibData.h"
+#include "AliCDBMetaData.h"
+#include "AliCDBId.h"
+#include "AliCDBEntry.h"
+#include "AliCDBManager.h"
+#include "AliCDBStorage.h"
+#endif
+
+
+void AliEMCALSetCDB()
+{
+  TControlBar *menu = new TControlBar("vertical","EMCAL CDB");
+  menu->AddButton("Help to run EMCAL CDB","Help()",
+                 "Explains how to use EMCAL CDS menus");
+  menu->AddButton("Equal CC","SetCC(0)",
+                 "Set equal calibration coefficients");
+  menu->AddButton("Decalibrate","SetCC(1)",
+                 "Set random decalibration calibration coefficients");
+  menu->AddButton("Read equal CC","GetCC(0)",
+                 "Read initial equal calibration coefficients");
+  menu->AddButton("Read random CC","GetCC(1)",
+                 "Read random decalibration calibration coefficients");
+  menu->Show();
+}
+
+//------------------------------------------------------------------------
+void Help()
+{
+  char *string =
+    "\nSet calibration parameters and write them into ALICE CDB. Press button \"Equal CC\" to create equal pedestals and gain factors. Press button \"Decalibrate\" to create random pedestals and gain factors to imitate decalibrated detector\n";
+  printf(string);
+}
+
+//------------------------------------------------------------------------
+void SetCC(Int_t flag=0)
+{
+  // Writing calibration coefficients into the Calibration DB
+  // Arguments:
+  //   flag=0: all calibration coefficients are equal
+  //   flag=1: all calibration coefficients random (decalibration)
+  // Author: Boris Polishchuk (Boris.Polichtchouk at cern.ch)
+
+  TString DBFolder;
+  Int_t firstRun   =  0;
+  Int_t lastRun    = 10;
+  Int_t beamPeriod =  1;
+  char* objFormat  = "";
+
+  if      (flag == 0) {
+    DBFolder  ="local://InitCalibDB";
+    firstRun  =  0;
+    lastRun   =  0;
+    objFormat = "EMCAL initial gain factors and pedestals";
+  }
+  else if (flag == 1) {
+    DBFolder  ="local://DeCalibDB";
+    firstRun  =  0;
+    lastRun   = 12;
+    objFormat = "EMCAL random pedestals and ADC gain factors (12x48x24)";
+  }
+  
+  AliEMCALCalibData *calibda=new AliEMCALCalibData("EMCAL");
+  
+  Float_t fADCpedestal = 0.005;
+  Float_t fADCchanel   = 0.0015;
+  
+  TRandom rn;
+  Int_t nSMod  = 12;
+  Int_t nCol   = 48;
+  Int_t nRow   = 24;
+  Int_t nRow2  = 12; //Modules 11 and 12 are half modules
+
+  for(Int_t supermodule=1; supermodule < nSMod+1; supermodule++) {
+    for(Int_t column=1; column< nCol+1; column++) {
+      if(supermodule > 10)
+       nRow = nRow2;
+      for(Int_t row=1; row< nRow+1; row++) {
+       if (flag == 1) {
+         // Decalibration:
+         // Spread calibration coefficients uniformly with
+         // Cmax/Cmin = 5, (Cmax-Cmin)/2 = 0.0015
+         // and pedestals 0.005 +-10%
+         fADCchanel  =rn.Uniform(0.00075,0.00375);
+         fADCpedestal=rn.Uniform(0.0045,0.0055);
+       }
+       calibda->SetADCchannel (supermodule,column,row,fADCchanel);
+//     cout<<"Set SM: "<<supermodule<<" col "<<column<<" row "<<row
+//         <<" chan "<<fADCchanel<<endl;
+       calibda->SetADCpedestal(supermodule,column,row,fADCpedestal);
+      }
+    }
+  }
+
+  //Store calibration data into database
+  
+  AliCDBMetaData md;
+  md.SetComment(objFormat);
+  md.SetBeamPeriod(beamPeriod);
+  md.SetResponsible("Boris Polichtchouk");
+  
+  AliCDBId id("EMCAL/Calib/GainFactors_and_Pedestals",firstRun,lastRun);
+
+  AliCDBManager* man = AliCDBManager::Instance();  
+  AliCDBStorage* loc = man->GetStorage(DBFolder.Data());
+  loc->Put(calibda, id, &md);
+
+}
+
+//------------------------------------------------------------------------
+void GetCC(Int_t flag=0)
+{
+  // Read calibration coefficients into the Calibration DB
+  // Arguments:
+  //   flag=0: all calibration coefficients are equal
+  //   flag=1: all calibration coefficients random (decalibration)
+  // Author: Yuri.Kharlov at cern.ch
+
+  TString DBFolder;
+
+  if      (flag == 0) {
+    DBFolder  ="local://InitCalibDB";
+  }
+  else if (flag == 1) {
+    DBFolder  ="local://DeCalibDB";
+  }
+
+  AliEMCALCalibData* clb  = (AliEMCALCalibData*)
+    (AliCDBManager::Instance()
+     ->GetStorage(DBFolder.Data())
+     ->Get("EMCAL/Calib/GainFactors_and_Pedestals",
+          gAlice->GetRunNumber())->GetObject());
+
+  Int_t nSMod = 12;
+  Int_t nCol  = 48;
+  Int_t nRow  = 24;
+  Int_t nRow2 = 12; //Modules 11 and 12 are half modules
+
+  TH2F *hPed[nSMod], *hGain[nSMod];
+  TCanvas *cPed  = new TCanvas("cPed" ,"Pedestals"   , 10,10,400,800);
+  TCanvas *cGain = new TCanvas("cGain","Gain factors",410,10,400,800);
+  cPed ->Divide(1,12);
+  cGain->Divide(1,12);
+
+  for (Int_t supermodule=1; supermodule<=nSMod; supermodule++) {
+
+    if(supermodule > 10)
+      nRow = nRow2;
+
+    TString namePed="hPed";
+    namePed+=supermodule;
+    TString titlePed="Pedestals in supermodule ";
+    titlePed+=supermodule;
+    hPed[supermodule-1] = new TH2F(namePed.Data(),titlePed.Data(),
+                           nCol,1.,1.*nCol,nRow,1.,1.*nRow);
+
+    TString nameGain="hGain";
+    nameGain+=supermodule;
+    TString titleGain="Gain factors in supermodule ";
+    titleGain+=supermodule;
+    hGain[supermodule-1] = new TH2F(nameGain.Data(),titleGain.Data(),
+                                   nCol,1.,1.*nCol,nRow,1.,1.*nRow);
+    for (Int_t column=1; column<=nCol; column++) {
+      for (Int_t row=1; row<=nRow; row++) {
+       Float_t ped  = clb->GetADCpedestal(supermodule,column,row);
+       Float_t gain = clb->GetADCchannel (supermodule,column,row);
+       cout<<"Get SM: "<<supermodule<<" col "<<column<<" row "<<row
+           <<" chan "<<gain<<endl;
+       hPed[supermodule-1] ->SetBinContent(column,row,ped);
+       hGain[supermodule-1]->SetBinContent(column,row,gain);
+      }
+    }
+    cPed ->cd(supermodule);
+    hPed[supermodule-1] ->Draw("lego2");
+    cGain->cd(supermodule);
+    hGain[supermodule-1]->Draw("lego2");
+  }
+  cPed ->Print("pedestals.eps");
+  cGain->Print("gains.eps");
+}