ZDC/MakeZDCRecoParam.C

   1 #if !defined(__CINT__) || defined(__MAKECINT__)
   2
   3 #include <Riostream.h>
   4 #include <TSystem.h>
   5 #include <TFile.h>
   6 #include <TH2F.h>
   7 #include <TH1D.h>
   8 #include "STEER/AliCDBEntry.h"
   9 #include "STEER/AliCDBGrid.h"
  10 #include "STEER/AliCDBId.h"
  11 #include "STEER/AliCDBLocal.h"
  12 #include "STEER/AliCDBManager.h"
  13 #include "STEER/AliCDBMetaData.h"
  14 #include "STEER/AliCDBPath.h"
  15 #include "STEER/AliCDBRunRange.h"
  16 #include "STEER/AliCDBStorage.h"
  17 #include "ZDC/AliZDC.h"
  18 #include "ZDC/AliZDCv3.h"
  19 #include "ZDC/AliZDCRecoParam.h"
  20 #include "ZDC/AliZDCRecoParampp.h"
  21 #include "ZDC/AliZDCRecoParamPbPb.h"
  22
  23 #endif
  24
  25 void MakeZDCRecoParam(Int_t type=0){
  26 //========================================================================
  27 //
  28 // Steering macro to create and store in OCDB
  29 //       ZDC reconstruction parameters
  30 //
  31 // Contact: chiara.oppedisano@to.infn.it
  32 //
  33 //========================================================================
  34
  35   AliCDBManager* cdb = AliCDBManager::Instance();
  36   //if(!cdb->IsDefaultStorageSet()) cdb->SetDefaultStorage("local://OCDB");
  37   cdb->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
  38
  39   if(type==0){
  40     AliZDCRecoParampp* zdcRecoParam = new AliZDCRecoParampp();
  41     // save in CDB storage
  42     AliCDBMetaData *md= new AliCDBMetaData();
  43     md->SetResponsible("Chiara Oppedisano");
  44     md->SetComment("p-p reconstruction parameters for ZDC");
  45     md->SetAliRootVersion(gSystem->Getenv("ARVERSION"));
  46     md->SetObjectClassName("AliZDCRecoParampp");
  47     AliCDBId id("ZDC/Calib/RecoParampp",0,AliCDBRunRange::Infinity());
  48     cdb->GetDefaultStorage()->Put(zdcRecoParam,id, md);
  49   }
  50   else if(type==1){
  51     //
  52     TH1::AddDirectory(0);
  53     TH2::AddDirectory(0);
  54     //
  55     TFile * fileHistos = TFile::Open("$ALICE_ROOT/ZDC/GlauberMCHistos.root");
  56     fileHistos->cd();
  57     //
  58     TH2F *hZDCvsZEM = (TH2F*) fileHistos->Get("hZDCvsZEM");
  59     hZDCvsZEM->SetDirectory(0);
  60     //
  61     TH2F *hZDCCvsZEM = (TH2F*) fileHistos->Get("hZDCCvsZEM");
  62     hZDCCvsZEM->SetDirectory(0);
  63     //
  64     TH2F *hZDCAvsZEM = (TH2F*) fileHistos->Get("hZDCAvsZEM");
  65     hZDCAvsZEM->SetDirectory(0);
  66     //
  67     TH1D* hDist = (TH1D*) fileHistos->Get("hDist");
  68     hDist->SetDirectory(0);
  69     //
  70     TH1D* hbDist = (TH1D*) fileHistos->Get("hbDist");
  71     hbDist->SetDirectory(0);
  72
  73     AliZDCRecoParamPbPb* zdcRecoParam =
  74         new AliZDCRecoParamPbPb(hZDCvsZEM, hZDCCvsZEM,
  75                 hZDCAvsZEM, hDist, hbDist, 0.1);
  76     //zdcRecoParam->Print("");
  77
  78     // save in CDB storage
  79     AliCDBMetaData *md= new AliCDBMetaData();
  80     md->SetResponsible("Chiara Oppedisano");
  81     md->SetComment("A-A reconstruction parameters for ZDC");
  82     md->SetAliRootVersion(gSystem->Getenv("ARVERSION"));
  83     md->SetObjectClassName("AliZDCRecoParamPbPb");
  84     AliCDBId id("ZDC/Calib/RecoParamPbPb", 0, AliCDBRunRange::Infinity());
  85     cdb->GetDefaultStorage()->Put(zdcRecoParam, id, md);
  86     //
  87     fileHistos->Close();
  88   }
  89   else{
  90     printf("Event type not implemented\n");
  91     return;
  92   }
  93
  94 }