[u/mrichter/AliRoot.git] / PWGPP / ZDC / PlotZDCQAtrain.C

#if !defined(__CINT__) || defined(__MAKECINT__)

#include <stdio.h>
#include <stdlib.h>
#include <TROOT.h>
#include <Riostream.h>
#include <TClassTable.h>
#include <TStyle.h>
#include <TMath.h>
#include <TFile.h>
#include <TCanvas.h>
#include <TH1.h>
#include <TH2.h>
#include <TProfile.h>
#include <TLine.h>
#include <TGrid.h>
#include <TBits.h>
#include <TChain.h>
#include <TNtuple.h>
#include <TTree.h>
#include <TBranch.h>
#include <TFileMerger.h>
#include <TGridResult.h>
#include <TSystem.h>

#endif

void PlotZDCQAtrain(TString option="local",
		    Int_t nRun=0,
		    TString period="LHC11h",
		    TString recoPass="pass1_HLT",
		    TString qaTrain="QA",
		    TString fileName="QAresults.root"){

  gStyle->SetOptStat(111);
  //  gStyle->SetOptTitle(0);
  
  TFile *f;
  TString path;
  Int_t year=2011;
  if(period.Contains("LHC10")) year=2010;
  else if(period.Contains("LHC11")) year=2011;

  if(option.Contains("local")){
    char fnamerun[20];
    sprintf(fnamerun,"QAresults%d.root",nRun);
    f = new TFile(fnamerun);
    printf("Opened file %s\n",f->GetName());
  }
  else{
    TGrid::Connect("alien:");
    char filenameAlien[250];
    sprintf(filenameAlien,"alien:///alice/data/%d/%s/%09d/ESDs/%s/%s",year,period.Data(),nRun,recoPass.Data(),fileName.Data());    
    f = TFile::Open(filenameAlien);
    printf("Opened file %s\n",f->GetName());
  }

  TDirectoryFile* df=(TDirectoryFile*)f->Get("ZDC_Performance");
  if(!df){
    printf("ZDC_Performance MISSING -> Exit\n");
    return;
  }
  TList* l=(TList*)df->Get("QAZDCHists");
  if(!df){
    printf("QAZDCHists TList MISSING -> Exit\n");
    return;
  }

  TH2F *fDebunch       = (TH2F*)l->FindObject("fDebunch");
  TH1F *fhTDCZNSum     = (TH1F*)l->FindObject("fhTDCZNSum");	 
  TH1F *fhTDCZNDiff    = (TH1F*)l->FindObject("fhTDCZNDiff");	 
  TH1F *fhZNCSpectrum  = (TH1F*)l->FindObject("fhZNCSpectrum");  
  TH1F *fhZNASpectrum  = (TH1F*)l->FindObject("fhZNASpectrum");  
  TH1F *fhZPCSpectrum  = (TH1F*)l->FindObject("fhZPCSpectrum");  
  TH1F *fhZPASpectrum  = (TH1F*)l->FindObject("fhZPASpectrum");  
  TH1F *fhZEM1Spectrum = (TH1F*)l->FindObject("fhZEM1Spectrum"); 
  TH1F *fhZEM2Spectrum = (TH1F*)l->FindObject("fhZEM2Spectrum"); 
  TH1F *fhZNCpmc   = (TH1F*)l->FindObject("fhZNCpmc");   
  TH1F *fhZNApmc   = (TH1F*)l->FindObject("fhZNApmc");   
  TH1F *fhZPCpmc   = (TH1F*)l->FindObject("fhZPCpmc");   
  TH1F *fhZPApmc   = (TH1F*)l->FindObject("fhZPApmc");  
  TH2F *fhZNCCentroid  = (TH2F*)l->FindObject("fhZNCCentroid");  
  TH2F *fhZNACentroid  = (TH2F*)l->FindObject("fhZNACentroid");   
  TH1F *fhZNCemd   = (TH1F*)l->FindObject("fhZNCemd");   
  TH1F *fhZNAemd   = (TH1F*)l->FindObject("fhZNAemd");   
  TH1F *fhPMCZNCemd    = (TH1F*)l->FindObject("fhPMCZNCemd");	 
  TH1F *fhPMCZNAemd    = (TH1F*)l->FindObject("fhPMCZNAemd");	 
 
  /*TH1D *hxZNC = fhZNCCentroid->ProjectionX("hxZNC");
  TH1D *hyZNC = fhZNCCentroid->ProjectionY("hyZNC");
  TH1D *hxZNA = fhZNACentroid->ProjectionX("hxZNA");
  TH1D *hyZNA = fhZNACentroid->ProjectionY("hyZNA");*/


  gStyle->SetPalette(1);
  
  char nam[25];
  sprintf(nam,"QAhistos%d.root",nRun);
  TFile *fileout = TFile::Open(nam,"RECREATE");
        fDebunch->Write();
  	fhTDCZNSum->Write();    
	fhTDCZNDiff ->Write();   
	fhZNCSpectrum->Write();  
	fhZNASpectrum ->Write(); 
	fhZPCSpectrum ->Write(); 
	fhZPASpectrum ->Write(); 
	fhZEM1Spectrum->Write(); 
	fhZEM2Spectrum ->Write();
	fhZNCpmc->Write();
	fhZNApmc->Write();
	fhZPCpmc->Write();
	fhZPApmc->Write();
	fhZNCCentroid->Write();  
	fhZNACentroid->Write();  
	fhZNCemd->Write();
	fhZNAemd->Write();
	fhPMCZNCemd->Write();    
	fhPMCZNAemd->Write();
	fhZNCCentroid->Write();
  
  fileout->Close();
}
Commit	Line	Data
e3284025	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