PWGPP/ZDC/PlotZDCQAtrain.C

   1 #if !defined(__CINT__) || defined(__MAKECINT__)
   2
   3 #include <stdio.h>
   4 #include <stdlib.h>
   5 #include <TROOT.h>
   6 #include <Riostream.h>
   7 #include <TClassTable.h>
   8 #include <TStyle.h>
   9 #include <TMath.h>
  10 #include <TFile.h>
  11 #include <TCanvas.h>
  12 #include <TH1.h>
  13 #include <TH2.h>
  14 #include <TProfile.h>
  15 #include <TLine.h>
  16 #include <TGrid.h>
  17 #include <TBits.h>
  18 #include <TChain.h>
  19 #include <TNtuple.h>
  20 #include <TTree.h>
  21 #include <TBranch.h>
  22 #include <TFileMerger.h>
  23 #include <TGridResult.h>
  24 #include <TSystem.h>
  25
  26 #endif
  27
  28 void PlotZDCQAtrain(TString option="local",
  29                     Int_t nRun=0,
  30                     TString period="LHC11h",
  31                     TString recoPass="pass1_HLT",
  32                     TString qaTrain="QA",
  33                     TString fileName="QAresults.root"){
  34
  35   gStyle->SetOptStat(111);
  36   //  gStyle->SetOptTitle(0);
  37
  38   TFile *f;
  39   TString path;
  40   Int_t year=2011;
  41   if(period.Contains("LHC10")) year=2010;
  42   else if(period.Contains("LHC11")) year=2011;
  43
  44   if(option.Contains("local")){
  45     char fnamerun[20];
  46     sprintf(fnamerun,"QAresults%d.root",nRun);
  47     f = new TFile(fnamerun);
  48     printf("Opened file %s\n",f->GetName());
  49   }
  50   else{
  51     TGrid::Connect("alien:");
  52     char filenameAlien[250];
  53     sprintf(filenameAlien,"alien:///alice/data/%d/%s/%09d/ESDs/%s/%s",year,period.Data(),nRun,recoPass.Data(),fileName.Data());
  54     f = TFile::Open(filenameAlien);
  55     printf("Opened file %s\n",f->GetName());
  56   }
  57
  58   TDirectoryFile* df=(TDirectoryFile*)f->Get("ZDC_Performance");
  59   if(!df){
  60     printf("ZDC_Performance MISSING -> Exit\n");
  61     return;
  62   }
  63   TList* l=(TList*)df->Get("QAZDCHists");
  64   if(!df){
  65     printf("QAZDCHists TList MISSING -> Exit\n");
  66     return;
  67   }
  68
  69   TH2F *fDebunch       = (TH2F*)l->FindObject("fDebunch");
  70   TH1F *fhTDCZNSum     = (TH1F*)l->FindObject("fhTDCZNSum");
  71   TH1F *fhTDCZNDiff    = (TH1F*)l->FindObject("fhTDCZNDiff");
  72   TH1F *fhZNCSpectrum  = (TH1F*)l->FindObject("fhZNCSpectrum");
  73   TH1F *fhZNASpectrum  = (TH1F*)l->FindObject("fhZNASpectrum");
  74   TH1F *fhZPCSpectrum  = (TH1F*)l->FindObject("fhZPCSpectrum");
  75   TH1F *fhZPASpectrum  = (TH1F*)l->FindObject("fhZPASpectrum");
  76   TH1F *fhZEM1Spectrum = (TH1F*)l->FindObject("fhZEM1Spectrum");
  77   TH1F *fhZEM2Spectrum = (TH1F*)l->FindObject("fhZEM2Spectrum");
  78   TH1F *fhZNCpmc   = (TH1F*)l->FindObject("fhZNCpmc");
  79   TH1F *fhZNApmc   = (TH1F*)l->FindObject("fhZNApmc");
  80   TH1F *fhZPCpmc   = (TH1F*)l->FindObject("fhZPCpmc");
  81   TH1F *fhZPApmc   = (TH1F*)l->FindObject("fhZPApmc");
  82   TH2F *fhZNCCentroid  = (TH2F*)l->FindObject("fhZNCCentroid");
  83   TH2F *fhZNACentroid  = (TH2F*)l->FindObject("fhZNACentroid");
  84   TH1F *fhZNCemd   = (TH1F*)l->FindObject("fhZNCemd");
  85   TH1F *fhZNAemd   = (TH1F*)l->FindObject("fhZNAemd");
  86   TH1F *fhPMCZNCemd    = (TH1F*)l->FindObject("fhPMCZNCemd");
  87   TH1F *fhPMCZNAemd    = (TH1F*)l->FindObject("fhPMCZNAemd");
  88
  89   /*TH1D *hxZNC = fhZNCCentroid->ProjectionX("hxZNC");
  90   TH1D *hyZNC = fhZNCCentroid->ProjectionY("hyZNC");
  91   TH1D *hxZNA = fhZNACentroid->ProjectionX("hxZNA");
  92   TH1D *hyZNA = fhZNACentroid->ProjectionY("hyZNA");*/
  93
  94
  95   gStyle->SetPalette(1);
  96
  97   char nam[25];
  98   sprintf(nam,"QAhistos%d.root",nRun);
  99   TFile *fileout = TFile::Open(nam,"RECREATE");
 100         fDebunch->Write();
 101         fhTDCZNSum->Write();
 102         fhTDCZNDiff ->Write();
 103         fhZNCSpectrum->Write();
 104         fhZNASpectrum ->Write();
 105         fhZPCSpectrum ->Write();
 106         fhZPASpectrum ->Write();
 107         fhZEM1Spectrum->Write();
 108         fhZEM2Spectrum ->Write();
 109         fhZNCpmc->Write();
 110         fhZNApmc->Write();
 111         fhZPCpmc->Write();
 112         fhZPApmc->Write();
 113         fhZNCCentroid->Write();
 114         fhZNACentroid->Write();
 115         fhZNCemd->Write();
 116         fhZNAemd->Write();
 117         fhPMCZNCemd->Write();
 118         fhPMCZNAemd->Write();
 119         fhZNCCentroid->Write();
 120
 121   fileout->Close();
 122 }
 123