[u/mrichter/AliRoot.git] / ZDC / ZDCQAtrending.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>
#include <TLegend.h>

#endif

void MakePlots(TString ntupleFileName);

void ZDCQAtrending(TString period,
		   TString recoPass,
		   TString qaTrain = "QA",
		   Bool_t useOnlyMerged = kTRUE,
		   Int_t firstRun = 0,
		   Int_t lastRun = 999999999,
		   TString runListFile = "",
		   TString fileName = "QAresults.root"){

  gStyle->SetOptStat(0);

  TGrid::Connect("alien:");
  Int_t year = 0;
  if(period.Contains("LHC09")) year = 2009;
  else if(period.Contains("LHC10")) year = 2010;
  else if(period.Contains("LHC11")) year = 2011;
  
  Bool_t useExternalList = kFALSE;
  Int_t runList[10000];
  Int_t totRuns = 0;
  if(runListFile.Length()>0){
    if(!gSystem->Exec(Form("ls -l %s > /dev/null 2>&1",runListFile.Data()))){
      printf("Use Run List from %s  --- runs to be analyzed:\n",runListFile.Data());
      useExternalList = kTRUE;
      FILE* rfil = fopen(runListFile.Data(),"r");
      Int_t nrun;
      while(!feof(rfil)){
	int stat = fscanf(rfil,"%d, ",&nrun);
	if(feof(rfil)) break;
	runList[totRuns++] = nrun;
      }
      for(Int_t ir = 0; ir<totRuns; ir++){
	printf("%d\n",runList[ir]);
      }
    }
    else{
      printf("File with run list does not exist\n");
    }
  }

  TString outFilNam = Form("ZDCtrend_%s_%s_%s.root",period.Data(),recoPass.Data(),qaTrain.Data());

 
  const Int_t nVariables = 27;
  TNtuple* ntzdc = new TNtuple("ntzdc","ZDC trending",
  "nrun:meanZNC:meanZPC:meanZNA:meanZPA:meanZEM1:meanZEM2:emeanZNC:emeanZPC:emeanZNA:emeanZPA:emeanZEM1:emeanZEM2:"
  "ZNClg:ZNAlg:eZNClg:eZNAlg:pmcZNClg:pmcZNAlg:epmcZNClg:epmcZNAlg:xZNC:yZNC:xZNA:yZNA:tdcSum:tdcDiff");
  Float_t xnt[nVariables];

  TBits* readRun = new TBits(999999);
  readRun->ResetAllBits();
  if(!useExternalList){
    if(!gSystem->Exec(Form("ls -l %s > /dev/null 2>&1",outFilNam.Data()))){
      TFile* oldfil = new TFile(outFilNam.Data());
      TNtuple* ntmp = (TNtuple*)oldfil->Get("ntzdc");
      Bool_t isOK = kFALSE;
      if(ntmp){
	if(ntmp->GetNvar() == ntzdc->GetNvar()){
	  isOK = kTRUE;
	  TObjArray* arr1 = (TObjArray*)ntzdc->GetListOfBranches();
	  TObjArray* arr2 = (TObjArray*)ntmp->GetListOfBranches();
	  for(Int_t iV = 0; iV<ntmp->GetNvar(); iV++){
	    TString vnam1 = arr1->At(iV)->GetName();
	    TString vnam2 = arr2->At(iV)->GetName();
	    if(vnam1 != vnam2) isOK = kFALSE;
	    ntmp->SetBranchAddress(vnam2.Data(),&xnt[iV]);
	  }
	  if(isOK){
	    for(Int_t nE = 0; nE<ntmp->GetEntries(); nE++){
	      ntmp->GetEvent(nE);
	      Int_t theRun = (Int_t)(xnt[0]+0.0001);
	      readRun->SetBitNumber(theRun);
	      ntzdc->Fill(xnt);
	    }
	  }
	}
      }
      if(!isOK){
	printf("Ntuple in local file not OK -> will be recreated\n");
      }
      oldfil->Close();
      delete oldfil;
    }
  }

  if(!gGrid||!gGrid->IsConnected()) {
    printf("gGrid not found! exit macro\n");
    return;
  }

  TString  path = Form("/alice/data/%d/%s/",year,period.Data());
  TGridResult *gr  =  gGrid->Query(path,fileName);
  Int_t nFiles  =  gr->GetEntries();
  printf(" --->%d merged QAresults.root files found in %s\n", nFiles,path.Data());
  if(nFiles < 1) return;

  Int_t nAnalyzedFiles = 0;
  if(nFiles > 1){
    for(Int_t iFil  =  0; iFil <nFiles ; iFil++) { 
      TString fileNameLong = Form("%s",gr->GetKey(iFil,"turl"));
      if(!fileNameLong.Contains(recoPass.Data())) continue;
// Commented Sep. 2011
//      if(!fileNameLong.Contains(qaTrain.Data())) continue;
//      if(fileNameLong.Contains("TRD") || fileNameLong.Contains("EMCAL")) continue;
      TString runNumber = fileNameLong;
      runNumber.ReplaceAll(Form("alien:///alice/data/%d/%s/",year,period.Data()),"");
      runNumber.Remove(9,runNumber.Sizeof());
   
      Int_t iRun = atoi(runNumber.Data());
      //printf("  runNumber: %s -> iRun %d\n",runNumber.Data(), iRun);
      if(useExternalList){
	Bool_t keepRun = kFALSE;
	for(Int_t ir = 0; ir<totRuns; ir++){
	  if(iRun == runList[ir]){
	    keepRun = kTRUE;
	    break;
	  }
	}
	if(!keepRun) continue;
      }
      if(readRun->TestBitNumber(iRun)){ 
	printf("Run %d already in local ntuple -> skipping it\n",iRun);
	continue;
      }
      if(iRun<firstRun) continue;
      if(iRun>lastRun) continue;    

printf("  fileNameLog: %s\n",fileNameLong.Data());
      if(useOnlyMerged){
	TString isMerged = fileNameLong;
	//
	if(!isMerged.Contains(recoPass.Data())) continue;
	if(isMerged.Contains("Stage")) continue;
printf("  ->  isMerged: %s\n",isMerged.Data());
	isMerged.Remove(0,isMerged.Sizeof()-20); 
printf("  ->  isMerged: %s\n",isMerged.Data());
	if(!isMerged.Contains(qaTrain) || !isMerged.Contains("QAresults.root")) continue;
printf("isMerged %s\n",isMerged.Data());
	
      }
      printf("Open File %s  \n",fileNameLong.Data());
      

      TFile* f = TFile::Open(fileNameLong.Data());  

      TDirectoryFile* df = (TDirectoryFile*)f->Get("ZDC_Performance");
      if(!df){
	printf("Run %d ZDC_Performance MISSING -> Exit\n",iRun);
	continue;
      }
      TList* l = (TList*)df->Get("QAZDCHists");
      if(!df){
	printf("Run %d QAZDCHists TList MISSING -> Exit\n",iRun);
	continue;
      }  

      nAnalyzedFiles++;
      if(!useOnlyMerged) readRun->SetBitNumber(iRun);

      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");
      
      Int_t index = 0;
      xnt[index++] = (Float_t)iRun;
      xnt[index++] = fhZNCSpectrum->GetMean();
      xnt[index++] = fhZPCSpectrum->GetMean();
      xnt[index++] = fhZNASpectrum->GetMean();
      xnt[index++] = fhZPASpectrum->GetMean();
      xnt[index++] = fhZEM1Spectrum->GetMean();
      xnt[index++] = fhZEM2Spectrum->GetMean();
      if(fhZNCSpectrum->GetEntries()>0)  xnt[index++] = fhZNCSpectrum->GetMean()/fhZNCSpectrum->GetEntries();
      if(fhZPCSpectrum->GetEntries()>0)  xnt[index++] = fhZPCSpectrum->GetMean()/fhZPCSpectrum->GetEntries();
      if(fhZNASpectrum->GetEntries()>0)  xnt[index++] = fhZNASpectrum->GetMean()/fhZNASpectrum->GetEntries();
      if(fhZPASpectrum->GetEntries()>0)  xnt[index++] = fhZPASpectrum->GetMean()/fhZPASpectrum->GetEntries();
      if(fhZEM1Spectrum->GetEntries()>0) xnt[index++] = fhZEM1Spectrum->GetMean()/fhZEM1Spectrum->GetEntries();
      if(fhZEM2Spectrum->GetEntries()>0) xnt[index++] = fhZEM2Spectrum->GetMean()/fhZEM2Spectrum->GetEntries();
      xnt[index++] = fhZNCemd->GetMean();
      xnt[index++] = fhZNAemd->GetMean();
      if(fhZNCemd->GetEntries()>0) xnt[index++] = fhZNCemd->GetMean()/fhZNCemd->GetEntries();
      if(fhZNAemd->GetEntries()>0) xnt[index++] = fhZNAemd->GetMean()/fhZNAemd->GetEntries();
      xnt[index++] = fhPMCZNCemd->GetMean();
      xnt[index++] = fhPMCZNAemd->GetMean();
      if(fhPMCZNCemd->GetEntries()>0) xnt[index++] = fhPMCZNCemd->GetMean()/fhPMCZNCemd->GetEntries();
      if(fhPMCZNAemd->GetEntries()>0) xnt[index++] = fhPMCZNAemd->GetMean()/fhPMCZNAemd->GetEntries();
      xnt[index++] = hxZNC->GetMean();
      xnt[index++] = hyZNC->GetMean();
      xnt[index++] = hxZNA->GetMean();
      xnt[index++] = hyZNA->GetMean();
      xnt[index++] = fhTDCZNSum->GetMean();
      xnt[index++] = fhTDCZNDiff->GetMean();
      ntzdc->Fill(xnt);
    }
  }
  printf("Number of analyzed files  =  %d\n",nAnalyzedFiles);

  if(nAnalyzedFiles>0){
    TFile* outfil = new TFile(outFilNam.Data(),"recreate");
    outfil->cd();
    ntzdc->Write();
    outfil->Close();
    
    MakePlots(outFilNam);
  }
}

void MakePlots(TString ntupleFileName){
  TFile* fil = new TFile(ntupleFileName.Data(),"read");
  if(!fil){
    printf("File with ntuple does not exist\n");
    return;
  }
  TNtuple* ntzdc = (TNtuple*)fil->Get("ntzdc");

  Float_t nrun;
  Float_t meanZNC,meanZPC,meanZNA,meanZPA,meanZEM1,meanZEM2;
  Float_t emeanZNC,emeanZPC,emeanZNA,emeanZPA,emeanZEM1,emeanZEM2;
  Float_t ZNClg,ZNAlg,pmcZNClg,pmcZNAlg;
  Float_t eZNClg,eZNAlg,epmcZNClg,epmcZNAlg;
  Float_t xZNC,yZNC,xZNA,yZNA,tdcSum,tdcDiff;
  
  ntzdc->SetBranchAddress("nrun",&nrun);
  ntzdc->SetBranchAddress("meanZNC",&meanZNC);
  ntzdc->SetBranchAddress("meanZPC",&meanZPC);
  ntzdc->SetBranchAddress("meanZNA",&meanZNA);
  ntzdc->SetBranchAddress("meanZPA",&meanZPA);
  ntzdc->SetBranchAddress("meanZEM1",&meanZEM1);
  ntzdc->SetBranchAddress("meanZEM2",&meanZEM2);
  ntzdc->SetBranchAddress("emeanZNC",&emeanZNC);
  ntzdc->SetBranchAddress("emeanZPC",&emeanZPC);
  ntzdc->SetBranchAddress("emeanZNA",&emeanZNA);
  ntzdc->SetBranchAddress("emeanZPA",&emeanZPA);
  ntzdc->SetBranchAddress("emeanZEM1",&emeanZEM1);
  ntzdc->SetBranchAddress("emeanZEM2",&emeanZEM2);
  ntzdc->SetBranchAddress("ZNClg",&ZNClg);
  ntzdc->SetBranchAddress("ZNAlg",&ZNAlg);
  ntzdc->SetBranchAddress("eZNClg",&eZNClg);
  ntzdc->SetBranchAddress("eZNAlg",&eZNAlg);
  ntzdc->SetBranchAddress("pmcZNClg",&pmcZNClg);
  ntzdc->SetBranchAddress("pmcZNAlg",&pmcZNAlg);
  ntzdc->SetBranchAddress("epmcZNClg",&epmcZNClg);
  ntzdc->SetBranchAddress("epmcZNAlg",&epmcZNAlg);
  ntzdc->SetBranchAddress("xZNC",&xZNC);
  ntzdc->SetBranchAddress("yZNC",&yZNC);
  ntzdc->SetBranchAddress("xZNA",&xZNA);
  ntzdc->SetBranchAddress("yZNA",&yZNA);
  ntzdc->SetBranchAddress("tdcSum",&tdcSum);
  ntzdc->SetBranchAddress("tdcDiff",&tdcDiff);

  TH1F *hznc      = new TH1F("hznc","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
  TH1F *hzna      = new TH1F("hzna","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
  TH1F *hzpc      = new TH1F("hzpc","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
  TH1F *hzpa      = new TH1F("hzpa","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
  TH1F *hzem1     = new TH1F("hzem1","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
  TH1F *hzem2     = new TH1F("hzem2","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
  TH1F *hznclg    = new TH1F("hznclg","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
  TH1F *hznalg    = new TH1F("hznalg","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
  TH1F *hpmcznclg = new TH1F("hpmcznclg","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
  TH1F *hpmcznalg = new TH1F("hpmcznalg","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
  TH1F *hxznc     = new TH1F("hxznc","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
  TH1F *hyznc     = new TH1F("hyznc","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
  TH1F *hxzna     = new TH1F("hxzna","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
  TH1F *hyzna     = new TH1F("hyzna","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
  TH1F *htdcsum   = new TH1F("htdcsum","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
  TH1F *htdcdiff  = new TH1F("htdcdiff","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());

  for(Int_t i = 0; i<ntzdc->GetEntries();i++){
    ntzdc->GetEvent(i);
    //
    hznc->SetBinContent(i+1, meanZNC);
    hznc->SetBinError(i+1, emeanZNC);
    hznc->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
    hznc->GetXaxis()->SetTitle("RUN #");
    hznc->GetYaxis()->SetTitle("ZNC mean signal");
    hzpc->SetBinContent(i+1, meanZPC);
    hzpc->SetBinError(i+1, emeanZPC);
    hzpc->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
    hzpc->GetXaxis()->SetTitle("RUN #");
    hzpc->GetYaxis()->SetTitle("ZPC mean signal");
    hzna->SetBinContent(i+1, meanZNA);
    hzna->SetBinError(i+1, emeanZNA);
    hzna->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
    hzna->GetXaxis()->SetTitle("RUN #");
    hzna->GetYaxis()->SetTitle("ZNA mean signal");
    hzpa->SetBinContent(i+1, meanZPA);
    hzpa->SetBinError(i+1, emeanZPA);
    hzpa->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
    hzpa->GetXaxis()->SetTitle("RUN #");
    hzpa->GetYaxis()->SetTitle("ZPA mean signal");
    hzem1->SetBinContent(i+1, meanZEM1);
    hzem1->SetBinError(i+1, emeanZEM1);
    hzem1->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
    hzem1->GetXaxis()->SetTitle("RUN #");
    hzem1->GetYaxis()->SetTitle("ZEM1 mean signal");
    hzem2->SetBinContent(i+1, meanZEM2);
    hzem2->SetBinError(i+1, emeanZEM2);
    hzem2->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
    hzem2->GetXaxis()->SetTitle("RUN #");
    hzem2->GetYaxis()->SetTitle("ZEM1 mean signal");
    hznclg->SetBinContent(i+1, ZNClg);
    hznclg->SetBinError(i+1, eZNClg);
    hznclg->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
    hznclg->GetXaxis()->SetTitle("RUN #");
    hznclg->GetYaxis()->SetTitle("ZNC LG mean signal");
    hznalg->SetBinContent(i+1, ZNAlg);
    hznalg->SetBinError(i+1, eZNAlg);
    hznalg->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
    hznalg->GetXaxis()->SetTitle("RUN #");
    hznalg->GetYaxis()->SetTitle("ZNA LG mean signal");
    hpmcznclg->SetBinContent(i+1, pmcZNClg);
    hpmcznclg->SetBinError(i+1, epmcZNClg);
    hpmcznclg->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
    hpmcznclg->GetXaxis()->SetTitle("RUN #");
    hpmcznclg->GetYaxis()->SetTitle("ZNC LG PMC mean signal");
    hpmcznalg->SetBinContent(i+1, pmcZNAlg);
    hpmcznalg->SetBinError(i+1, epmcZNAlg);
    hpmcznalg->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
    hpmcznalg->GetXaxis()->SetTitle("RUN #");
    hpmcznalg->GetYaxis()->SetTitle("ZNC LG PMC mean signal");
    hxznc->SetBinContent(i+1, xZNC);
    hxznc->SetBinError(i+1, 0.1*xZNC);
    hxznc->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
    hxznc->GetXaxis()->SetTitle("RUN #");
    hxznc->GetYaxis()->SetTitle("X_{ZN} (cm)");
    hyznc->SetBinContent(i+1, yZNC);
    hyznc->SetBinError(i+1, 0.1*yZNC);
    hyznc->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
    hyznc->GetXaxis()->SetTitle("RUN #");
    hyznc->GetYaxis()->SetTitle("Y_{ZN} (cm)");
    hxzna->SetBinContent(i+1, xZNA);
    hxzna->SetBinError(i+1, 0.1*xZNA);
    hxzna->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
    hxzna->GetXaxis()->SetTitle("RUN #");
    hxzna->GetYaxis()->SetTitle("X_{ZN} (cm)");
    hyzna->SetBinContent(i+1, yZNA);
    hyzna->SetBinError(i+1, 0.1*yZNA);
    hyzna->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
    hyzna->GetXaxis()->SetTitle("RUN #");
    hyzna->GetYaxis()->SetTitle("Y_{ZN} (cm)");
    htdcsum->SetBinContent(i+1, tdcSum);
    htdcsum->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
    htdcsum->GetXaxis()->SetTitle("RUN #");
    htdcsum->GetYaxis()->SetTitle("TDC Sum (ns)");
    htdcdiff->SetBinContent(i+1, tdcDiff);
    htdcdiff->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
    htdcdiff->GetXaxis()->SetTitle("RUN #");
    htdcdiff->GetYaxis()->SetTitle("TDC Diff (ns)");
  }

  
  TCanvas *c1 = new TCanvas("c1", "Mean value ZNs", 0, 0, 1200, 1000);
  c1->Divide(1,2);
  //
  c1->cd(1); 
  hznc->SetMarkerColor(kAzure+6); hznc->SetLineColor(kAzure+6); 
  hznc->SetMarkerStyle(21);
  hznc->Draw(""); 
  hzna->SetMarkerColor(kPink-2);  hzna->SetLineColor(kPink-2);
  hzna->SetMarkerStyle(20);
  hzna->Draw("SAME"); 
  //
  TLegend *l1 = new TLegend(0.44,0.18,0.54,0.32);
  l1->SetFillColor(kWhite);
  l1->AddEntry(hznc," ZNC " ,"P");
  l1->AddEntry(hzna," ZNA " ,"P");
  l1->Draw("");
  //
  c1->cd(2);
  hznalg->SetMarkerColor(kBlue+1); hznalg->SetLineColor(kBlue+1); 
  hznalg->SetMarkerStyle(20);  hznalg->SetMinimum(0);
  hznalg->Draw("");
  hznclg->SetMarkerColor(kPink+5); hznclg->SetLineColor(kPink+5); 
  hznclg->SetMarkerStyle(21);
  hznclg->Draw("SAME"); 
  //
  TLegend *l2 = new TLegend(0.44,0.18,0.54,0.32);
  l2->SetFillColor(kWhite);
  l2->AddEntry(hznc," ZNC LG " ,"P");
  l2->AddEntry(hzna," ZNA LG " ,"P");
  l2->Draw("");
  

  TCanvas *c1b = new TCanvas("c1b", "Mean value ZEMs ZPs", 200, 0, 1200, 1000);
  c1b->Divide(1,2);
  c1b->cd(1); 
  hzpc->SetMarkerColor(kAzure+6); hzpc->SetLineColor(kAzure+6); 
  hzpc->SetMarkerStyle(21);
  hzpc->Draw(""); 
  hzpa->SetMarkerColor(kPink-2);  hzpa->SetLineColor(kPink-2);
  hzpa->SetMarkerStyle(20);
  hzpa->Draw("SAME"); 
  //
  TLegend *l3 = new TLegend(0.44,0.18,0.54,0.32);
  l3->SetFillColor(kWhite);
  l3->AddEntry(hzpc," ZPC " ,"P");
  l3->AddEntry(hzpa," ZPA " ,"P");
  l3->Draw("");
  
  c1b->cd(2);
  hzem1->SetMarkerColor(kTeal-7); hzem1->SetLineColor(kTeal-7);
  hzem1->SetMarkerStyle(29);
  hzem2->SetMarkerColor(kTeal+5); hzem2->SetLineColor(kTeal+5);
  hzem2->SetMarkerStyle(30);
  hzem2->Draw(""); 
  hzem1->Draw("SAME"); 
  //
  TLegend *l4 = new TLegend(0.44,0.18,0.54,0.32);
  l4->SetFillColor(kWhite);
  l4->AddEntry(hzem1," ZEM1 " ,"P");
  l4->AddEntry(hzem2," ZEM2 " ,"P");
  l4->Draw("");
 
  
  /*TCanvas *c2 = new TCanvas("c2", "ZN centroids", 400, 400, 1200, 1000);
  c2->Divide(1,2);
  //
  c2->cd(1);
  hxznc->SetMarkerColor(kAzure+5); hxznc->SetLineColor(kAzure+5);
  hxznc->SetMarkerStyle(21); 
  hxzna->SetMarkerColor(kPink+5); hxzna->SetLineColor(kPink+5);
  hxzna->SetMarkerStyle(20); 
  hxznc->Draw(""); 
  hxzna->Draw("SAME");
  c2->cd(2);
  hyznc->SetMarkerColor(kAzure+5); hyznc->SetLineColor(kAzure+5);
  hyznc->SetMarkerStyle(21); 
  hyzna->SetMarkerColor(kPink+5); hyzna->SetLineColor(kPink+5);
  hyzna->SetMarkerStyle(20); 
  hyznc->Draw(""); 
  hyzna->Draw("SAME");*/

  
  TCanvas *c3 = new TCanvas("c3", "LG signals", 600, 0, 1200, 1000);
  c3->Divide(1,2);
  //
  c3->cd(1);
  hznalg->SetMarkerColor(kBlue+1); hznalg->SetLineColor(kBlue+1); 
  hznalg->SetMarkerStyle(20);
  hznalg->Draw(""); 
  hznclg->SetMarkerColor(kPink+5); hznclg->SetLineColor(kPink+5); 
  hznclg->SetMarkerStyle(21);
  hznclg->Draw("SAME"); 
  TLegend *l7 = new TLegend(0.44,0.18,0.54,0.32);
  l7->SetFillColor(kWhite);
  l7->AddEntry(hznclg," ZNC LG " ,"P");
  l7->AddEntry(hznalg," ZNA LG " ,"P");
  l7->Draw("");
  //
  c3->cd(2);
  hpmcznalg->SetMarkerColor(kAzure+7); hpmcznalg->SetLineColor(kAzure+7); 
  hpmcznalg->SetMarkerStyle(20); hpmcznalg->SetMinimum(0);
  hpmcznalg->Draw(""); 
  hpmcznclg->SetMarkerColor(kPink+6); hpmcznclg->SetLineColor(kPink+6); 
  hpmcznclg->SetMarkerStyle(21);
  hpmcznclg->Draw("SAME"); 
  TLegend *l8 = new TLegend(0.44,0.15,0.58,0.32);
  l8->SetFillColor(kWhite);
  l8->AddEntry(hpmcznclg," ZNC LG PMC" ,"P");
  l8->AddEntry(hpmcznalg," ZNA LG PMC" ,"P");
  l8->Draw("");
}
Commit	Line	Data
4eb3f0ea	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	#include <TLegend.h>
	26
	27	#endif
	28
	29	void MakePlots(TString ntupleFileName);
	30
	31	void ZDCQAtrending(TString period,
	32	TString recoPass,
	33	TString qaTrain = "QA",
	34	Bool_t useOnlyMerged = kTRUE,
	35	Int_t firstRun = 0,
	36	Int_t lastRun = 999999999,
	37	TString runListFile = "",
	38	TString fileName = "QAresults.root"){
	39
	40	gStyle->SetOptStat(0);
	41
	42	TGrid::Connect("alien:");
	43	Int_t year = 0;
	44	if(period.Contains("LHC09")) year = 2009;
	45	else if(period.Contains("LHC10")) year = 2010;
	46	else if(period.Contains("LHC11")) year = 2011;
	47
	48	Bool_t useExternalList = kFALSE;
	49	Int_t runList[10000];
f286e6ab	50	Int_t totRuns = 0;
4eb3f0ea	51	if(runListFile.Length()>0){
	52	if(!gSystem->Exec(Form("ls -l %s > /dev/null 2>&1",runListFile.Data()))){
	53	printf("Use Run List from %s --- runs to be analyzed:\n",runListFile.Data());
	54	useExternalList = kTRUE;
	55	FILE* rfil = fopen(runListFile.Data(),"r");
	56	Int_t nrun;
	57	while(!feof(rfil)){
f286e6ab	58	int stat = fscanf(rfil,"%d, ",&nrun);
4eb3f0ea	59	if(feof(rfil)) break;
	60	runList[totRuns++] = nrun;
	61	}
	62	for(Int_t ir = 0; ir<totRuns; ir++){
	63	printf("%d\n",runList[ir]);
	64	}
f286e6ab	65	}
f286e6ab	66	else{
4eb3f0ea	67	printf("File with run list does not exist\n");
	68	}
	69	}
	70
	71	TString outFilNam = Form("ZDCtrend_%s_%s_%s.root",period.Data(),recoPass.Data(),qaTrain.Data());
	72
	73
	74	const Int_t nVariables = 27;
	75	TNtuple* ntzdc = new TNtuple("ntzdc","ZDC trending",
	76	"nrun:meanZNC:meanZPC:meanZNA:meanZPA:meanZEM1:meanZEM2:emeanZNC:emeanZPC:emeanZNA:emeanZPA:emeanZEM1:emeanZEM2:"
	77	"ZNClg:ZNAlg:eZNClg:eZNAlg:pmcZNClg:pmcZNAlg:epmcZNClg:epmcZNAlg:xZNC:yZNC:xZNA:yZNA:tdcSum:tdcDiff");
	78	Float_t xnt[nVariables];
	79
	80	TBits* readRun = new TBits(999999);
	81	readRun->ResetAllBits();
	82	if(!useExternalList){
	83	if(!gSystem->Exec(Form("ls -l %s > /dev/null 2>&1",outFilNam.Data()))){
	84	TFile* oldfil = new TFile(outFilNam.Data());
	85	TNtuple* ntmp = (TNtuple*)oldfil->Get("ntzdc");
	86	Bool_t isOK = kFALSE;
	87	if(ntmp){
	88	if(ntmp->GetNvar() == ntzdc->GetNvar()){
	89	isOK = kTRUE;
	90	TObjArray* arr1 = (TObjArray*)ntzdc->GetListOfBranches();
	91	TObjArray* arr2 = (TObjArray*)ntmp->GetListOfBranches();
	92	for(Int_t iV = 0; iV<ntmp->GetNvar(); iV++){
	93	TString vnam1 = arr1->At(iV)->GetName();
	94	TString vnam2 = arr2->At(iV)->GetName();
	95	if(vnam1 != vnam2) isOK = kFALSE;
	96	ntmp->SetBranchAddress(vnam2.Data(),&xnt[iV]);
	97	}
	98	if(isOK){
	99	for(Int_t nE = 0; nE<ntmp->GetEntries(); nE++){
	100	ntmp->GetEvent(nE);
	101	Int_t theRun = (Int_t)(xnt[0]+0.0001);
	102	readRun->SetBitNumber(theRun);
	103	ntzdc->Fill(xnt);
	104	}
	105	}
	106	}
	107	}
	108	if(!isOK){
	109	printf("Ntuple in local file not OK -> will be recreated\n");
	110	}
	111	oldfil->Close();
	112	delete oldfil;
	113	}
	114	}
	115
	116	if(!gGrid\|\|!gGrid->IsConnected()) {
	117	printf("gGrid not found! exit macro\n");
	118	return;
	119	}
	120
	121	TString path = Form("/alice/data/%d/%s/",year,period.Data());
	122	TGridResult *gr = gGrid->Query(path,fileName);
	123	Int_t nFiles = gr->GetEntries();
f286e6ab	124	printf(" --->%d merged QAresults.root files found in %s\n", nFiles,path.Data());
f286e6ab	125	if(nFiles < 1) return;
4eb3f0ea	126
	127	Int_t nAnalyzedFiles = 0;
	128	if(nFiles > 1){
	129	for(Int_t iFil = 0; iFil <nFiles ; iFil++) {
	130	TString fileNameLong = Form("%s",gr->GetKey(iFil,"turl"));
	131	if(!fileNameLong.Contains(recoPass.Data())) continue;
f286e6ab	132	// Commented Sep. 2011
	133	// if(!fileNameLong.Contains(qaTrain.Data())) continue;
	134	// if(fileNameLong.Contains("TRD") \|\| fileNameLong.Contains("EMCAL")) continue;
4eb3f0ea	135	TString runNumber = fileNameLong;
	136	runNumber.ReplaceAll(Form("alien:///alice/data/%d/%s/",year,period.Data()),"");
	137	runNumber.Remove(9,runNumber.Sizeof());
	138
	139	Int_t iRun = atoi(runNumber.Data());
f286e6ab	140	//printf(" runNumber: %s -> iRun %d\n",runNumber.Data(), iRun);
4eb3f0ea	141	if(useExternalList){
	142	Bool_t keepRun = kFALSE;
	143	for(Int_t ir = 0; ir<totRuns; ir++){
	144	if(iRun == runList[ir]){
	145	keepRun = kTRUE;
	146	break;
	147	}
	148	}
	149	if(!keepRun) continue;
	150	}
	151	if(readRun->TestBitNumber(iRun)){
	152	printf("Run %d already in local ntuple -> skipping it\n",iRun);
	153	continue;
	154	}
	155	if(iRun<firstRun) continue;
	156	if(iRun>lastRun) continue;
	157
f286e6ab	158	printf(" fileNameLog: %s\n",fileNameLong.Data());
4eb3f0ea	159	if(useOnlyMerged){
4eb3f0ea	160	TString isMerged = fileNameLong;
f286e6ab	161	//
	162	if(!isMerged.Contains(recoPass.Data())) continue;
	163	if(isMerged.Contains("Stage")) continue;
	164	printf(" -> isMerged: %s\n",isMerged.Data());
	165	isMerged.Remove(0,isMerged.Sizeof()-20);
	166	printf(" -> isMerged: %s\n",isMerged.Data());
	167	if(!isMerged.Contains(qaTrain) \|\| !isMerged.Contains("QAresults.root")) continue;
	168	printf("isMerged %s\n",isMerged.Data());
	169
4eb3f0ea	170	}
f286e6ab	171	printf("Open File %s \n",fileNameLong.Data());
4eb3f0ea	172
	173
	174	TFile* f = TFile::Open(fileNameLong.Data());
	175
	176	TDirectoryFile* df = (TDirectoryFile*)f->Get("ZDC_Performance");
	177	if(!df){
	178	printf("Run %d ZDC_Performance MISSING -> Exit\n",iRun);
	179	continue;
	180	}
	181	TList* l = (TList*)df->Get("QAZDCHists");
	182	if(!df){
	183	printf("Run %d QAZDCHists TList MISSING -> Exit\n",iRun);
	184	continue;
	185	}
	186
	187	nAnalyzedFiles++;
	188	if(!useOnlyMerged) readRun->SetBitNumber(iRun);
	189
	190	TH1F fhTDCZNSum = (TH1F)l->FindObject("fhTDCZNSum");
	191	TH1F fhTDCZNDiff = (TH1F)l->FindObject("fhTDCZNDiff");
	192	TH1F fhZNCSpectrum = (TH1F)l->FindObject("fhZNCSpectrum");
	193	TH1F fhZNASpectrum = (TH1F)l->FindObject("fhZNASpectrum");
	194	TH1F fhZPCSpectrum = (TH1F)l->FindObject("fhZPCSpectrum");
	195	TH1F fhZPASpectrum = (TH1F)l->FindObject("fhZPASpectrum");
	196	TH1F fhZEM1Spectrum = (TH1F)l->FindObject("fhZEM1Spectrum");
	197	TH1F fhZEM2Spectrum = (TH1F)l->FindObject("fhZEM2Spectrum");
	198	/TH1F fhZNCpmc = (TH1F*)l->FindObject("fhZNCpmc");
	199	TH1F fhZNApmc = (TH1F)l->FindObject("fhZNApmc");
	200	TH1F fhZPCpmc = (TH1F)l->FindObject("fhZPCpmc");
	201	TH1F fhZPApmc = (TH1F)l->FindObject("fhZPApmc");*/
	202	TH2F fhZNCCentroid = (TH2F)l->FindObject("fhZNCCentroid");
	203	TH2F fhZNACentroid = (TH2F)l->FindObject("fhZNACentroid");
	204	TH1F fhZNCemd = (TH1F)l->FindObject("fhZNCemd");
	205	TH1F fhZNAemd = (TH1F)l->FindObject("fhZNAemd");
	206	TH1F fhPMCZNCemd = (TH1F)l->FindObject("fhPMCZNCemd");
	207	TH1F fhPMCZNAemd = (TH1F)l->FindObject("fhPMCZNAemd");
	208
	209	TH1D *hxZNC = fhZNCCentroid->ProjectionX("hxZNC");
	210	TH1D *hyZNC = fhZNCCentroid->ProjectionY("hyZNC");
	211	TH1D *hxZNA = fhZNACentroid->ProjectionX("hxZNA");
	212	TH1D *hyZNA = fhZNACentroid->ProjectionY("hyZNA");
	213
	214	Int_t index = 0;
	215	xnt[index++] = (Float_t)iRun;
	216	xnt[index++] = fhZNCSpectrum->GetMean();
	217	xnt[index++] = fhZPCSpectrum->GetMean();
	218	xnt[index++] = fhZNASpectrum->GetMean();
	219	xnt[index++] = fhZPASpectrum->GetMean();
	220	xnt[index++] = fhZEM1Spectrum->GetMean();
	221	xnt[index++] = fhZEM2Spectrum->GetMean();
	222	if(fhZNCSpectrum->GetEntries()>0) xnt[index++] = fhZNCSpectrum->GetMean()/fhZNCSpectrum->GetEntries();
	223	if(fhZPCSpectrum->GetEntries()>0) xnt[index++] = fhZPCSpectrum->GetMean()/fhZPCSpectrum->GetEntries();
	224	if(fhZNASpectrum->GetEntries()>0) xnt[index++] = fhZNASpectrum->GetMean()/fhZNASpectrum->GetEntries();
	225	if(fhZPASpectrum->GetEntries()>0) xnt[index++] = fhZPASpectrum->GetMean()/fhZPASpectrum->GetEntries();
	226	if(fhZEM1Spectrum->GetEntries()>0) xnt[index++] = fhZEM1Spectrum->GetMean()/fhZEM1Spectrum->GetEntries();
	227	if(fhZEM2Spectrum->GetEntries()>0) xnt[index++] = fhZEM2Spectrum->GetMean()/fhZEM2Spectrum->GetEntries();
	228	xnt[index++] = fhZNCemd->GetMean();
	229	xnt[index++] = fhZNAemd->GetMean();
	230	if(fhZNCemd->GetEntries()>0) xnt[index++] = fhZNCemd->GetMean()/fhZNCemd->GetEntries();
	231	if(fhZNAemd->GetEntries()>0) xnt[index++] = fhZNAemd->GetMean()/fhZNAemd->GetEntries();
	232	xnt[index++] = fhPMCZNCemd->GetMean();
	233	xnt[index++] = fhPMCZNAemd->GetMean();
	234	if(fhPMCZNCemd->GetEntries()>0) xnt[index++] = fhPMCZNCemd->GetMean()/fhPMCZNCemd->GetEntries();
	235	if(fhPMCZNAemd->GetEntries()>0) xnt[index++] = fhPMCZNAemd->GetMean()/fhPMCZNAemd->GetEntries();
236	xnt[index++] = hxZNC->GetMean();
237	xnt[index++] = hyZNC->GetMean();
238	xnt[index++] = hxZNA->GetMean();
239	xnt[index++] = hyZNA->GetMean();
240	xnt[index++] = fhTDCZNSum->GetMean();
241	xnt[index++] = fhTDCZNDiff->GetMean();
242	ntzdc->Fill(xnt);
243	}
244	}
245	printf("Number of analyzed files = %d\n",nAnalyzedFiles);
246
247	if(nAnalyzedFiles>0){
248	TFile* outfil = new TFile(outFilNam.Data(),"recreate");
249	outfil->cd();
250	ntzdc->Write();
251	outfil->Close();
252
253	MakePlots(outFilNam);
254	}
255	}
256
257	void MakePlots(TString ntupleFileName){
258	TFile* fil = new TFile(ntupleFileName.Data(),"read");
259	if(!fil){
260	printf("File with ntuple does not exist\n");
261	return;
262	}
263	TNtuple* ntzdc = (TNtuple*)fil->Get("ntzdc");
264
265	Float_t nrun;
266	Float_t meanZNC,meanZPC,meanZNA,meanZPA,meanZEM1,meanZEM2;
267	Float_t emeanZNC,emeanZPC,emeanZNA,emeanZPA,emeanZEM1,emeanZEM2;
268	Float_t ZNClg,ZNAlg,pmcZNClg,pmcZNAlg;
269	Float_t eZNClg,eZNAlg,epmcZNClg,epmcZNAlg;
270	Float_t xZNC,yZNC,xZNA,yZNA,tdcSum,tdcDiff;
271
272	ntzdc->SetBranchAddress("nrun",&nrun);
273	ntzdc->SetBranchAddress("meanZNC",&meanZNC);
274	ntzdc->SetBranchAddress("meanZPC",&meanZPC);
275	ntzdc->SetBranchAddress("meanZNA",&meanZNA);
276	ntzdc->SetBranchAddress("meanZPA",&meanZPA);
277	ntzdc->SetBranchAddress("meanZEM1",&meanZEM1);
278	ntzdc->SetBranchAddress("meanZEM2",&meanZEM2);
279	ntzdc->SetBranchAddress("emeanZNC",&emeanZNC);
280	ntzdc->SetBranchAddress("emeanZPC",&emeanZPC);
281	ntzdc->SetBranchAddress("emeanZNA",&emeanZNA);
282	ntzdc->SetBranchAddress("emeanZPA",&emeanZPA);
283	ntzdc->SetBranchAddress("emeanZEM1",&emeanZEM1);
284	ntzdc->SetBranchAddress("emeanZEM2",&emeanZEM2);
285	ntzdc->SetBranchAddress("ZNClg",&ZNClg);
286	ntzdc->SetBranchAddress("ZNAlg",&ZNAlg);
287	ntzdc->SetBranchAddress("eZNClg",&eZNClg);
288	ntzdc->SetBranchAddress("eZNAlg",&eZNAlg);
289	ntzdc->SetBranchAddress("pmcZNClg",&pmcZNClg);
290	ntzdc->SetBranchAddress("pmcZNAlg",&pmcZNAlg);
291	ntzdc->SetBranchAddress("epmcZNClg",&epmcZNClg);
292	ntzdc->SetBranchAddress("epmcZNAlg",&epmcZNAlg);
293	ntzdc->SetBranchAddress("xZNC",&xZNC);
294	ntzdc->SetBranchAddress("yZNC",&yZNC);
295	ntzdc->SetBranchAddress("xZNA",&xZNA);
296	ntzdc->SetBranchAddress("yZNA",&yZNA);
297	ntzdc->SetBranchAddress("tdcSum",&tdcSum);
298	ntzdc->SetBranchAddress("tdcDiff",&tdcDiff);
299
300	TH1F *hznc = new TH1F("hznc","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
301	TH1F *hzna = new TH1F("hzna","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
302	TH1F *hzpc = new TH1F("hzpc","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
303	TH1F *hzpa = new TH1F("hzpa","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
304	TH1F *hzem1 = new TH1F("hzem1","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
305	TH1F *hzem2 = new TH1F("hzem2","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
306	TH1F *hznclg = new TH1F("hznclg","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
307	TH1F *hznalg = new TH1F("hznalg","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
308	TH1F *hpmcznclg = new TH1F("hpmcznclg","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
309	TH1F *hpmcznalg = new TH1F("hpmcznalg","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
310	TH1F *hxznc = new TH1F("hxznc","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
311	TH1F *hyznc = new TH1F("hyznc","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
312	TH1F *hxzna = new TH1F("hxzna","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
313	TH1F *hyzna = new TH1F("hyzna","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
314	TH1F *htdcsum = new TH1F("htdcsum","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
315	TH1F *htdcdiff = new TH1F("htdcdiff","",(Int_t)ntzdc->GetEntries(),0.,ntzdc->GetEntries());
316
317	for(Int_t i = 0; i<ntzdc->GetEntries();i++){
318	ntzdc->GetEvent(i);
319	//
320	hznc->SetBinContent(i+1, meanZNC);
321	hznc->SetBinError(i+1, emeanZNC);
322	hznc->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
323	hznc->GetXaxis()->SetTitle("RUN #");
324	hznc->GetYaxis()->SetTitle("ZNC mean signal");
325	hzpc->SetBinContent(i+1, meanZPC);
326	hzpc->SetBinError(i+1, emeanZPC);
327	hzpc->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
328	hzpc->GetXaxis()->SetTitle("RUN #");
329	hzpc->GetYaxis()->SetTitle("ZPC mean signal");
330	hzna->SetBinContent(i+1, meanZNA);
331	hzna->SetBinError(i+1, emeanZNA);
332	hzna->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
333	hzna->GetXaxis()->SetTitle("RUN #");
334	hzna->GetYaxis()->SetTitle("ZNA mean signal");
335	hzpa->SetBinContent(i+1, meanZPA);
336	hzpa->SetBinError(i+1, emeanZPA);
337	hzpa->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
338	hzpa->GetXaxis()->SetTitle("RUN #");
339	hzpa->GetYaxis()->SetTitle("ZPA mean signal");
340	hzem1->SetBinContent(i+1, meanZEM1);
341	hzem1->SetBinError(i+1, emeanZEM1);
342	hzem1->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
343	hzem1->GetXaxis()->SetTitle("RUN #");
344	hzem1->GetYaxis()->SetTitle("ZEM1 mean signal");
345	hzem2->SetBinContent(i+1, meanZEM2);
346	hzem2->SetBinError(i+1, emeanZEM2);
347	hzem2->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
348	hzem2->GetXaxis()->SetTitle("RUN #");
349	hzem2->GetYaxis()->SetTitle("ZEM1 mean signal");
350	hznclg->SetBinContent(i+1, ZNClg);
351	hznclg->SetBinError(i+1, eZNClg);
352	hznclg->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
353	hznclg->GetXaxis()->SetTitle("RUN #");
354	hznclg->GetYaxis()->SetTitle("ZNC LG mean signal");
355	hznalg->SetBinContent(i+1, ZNAlg);
356	hznalg->SetBinError(i+1, eZNAlg);
357	hznalg->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
358	hznalg->GetXaxis()->SetTitle("RUN #");
359	hznalg->GetYaxis()->SetTitle("ZNA LG mean signal");
360	hpmcznclg->SetBinContent(i+1, pmcZNClg);
361	hpmcznclg->SetBinError(i+1, epmcZNClg);
362	hpmcznclg->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
363	hpmcznclg->GetXaxis()->SetTitle("RUN #");
364	hpmcznclg->GetYaxis()->SetTitle("ZNC LG PMC mean signal");
365	hpmcznalg->SetBinContent(i+1, pmcZNAlg);
366	hpmcznalg->SetBinError(i+1, epmcZNAlg);
367	hpmcznalg->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
368	hpmcznalg->GetXaxis()->SetTitle("RUN #");
369	hpmcznalg->GetYaxis()->SetTitle("ZNC LG PMC mean signal");
370	hxznc->SetBinContent(i+1, xZNC);
371	hxznc->SetBinError(i+1, 0.1*xZNC);
372	hxznc->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
373	hxznc->GetXaxis()->SetTitle("RUN #");
374	hxznc->GetYaxis()->SetTitle("X_{ZN} (cm)");
375	hyznc->SetBinContent(i+1, yZNC);
376	hyznc->SetBinError(i+1, 0.1*yZNC);
377	hyznc->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
378	hyznc->GetXaxis()->SetTitle("RUN #");
379	hyznc->GetYaxis()->SetTitle("Y_{ZN} (cm)");
380	hxzna->SetBinContent(i+1, xZNA);
381	hxzna->SetBinError(i+1, 0.1*xZNA);
382	hxzna->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
383	hxzna->GetXaxis()->SetTitle("RUN #");
384	hxzna->GetYaxis()->SetTitle("X_{ZN} (cm)");
385	hyzna->SetBinContent(i+1, yZNA);
386	hyzna->SetBinError(i+1, 0.1*yZNA);
387	hyzna->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
388	hyzna->GetXaxis()->SetTitle("RUN #");
389	hyzna->GetYaxis()->SetTitle("Y_{ZN} (cm)");
390	htdcsum->SetBinContent(i+1, tdcSum);
391	htdcsum->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
392	htdcsum->GetXaxis()->SetTitle("RUN #");
393	htdcsum->GetYaxis()->SetTitle("TDC Sum (ns)");
394	htdcdiff->SetBinContent(i+1, tdcDiff);
395	htdcdiff->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
396	htdcdiff->GetXaxis()->SetTitle("RUN #");
397	htdcdiff->GetYaxis()->SetTitle("TDC Diff (ns)");
398	}
399
f286e6ab	400
4eb3f0ea	401	TCanvas *c1 = new TCanvas("c1", "Mean value ZNs", 0, 0, 1200, 1000);
	402	c1->Divide(1,2);
	403	//
	404	c1->cd(1);
	405	hznc->SetMarkerColor(kAzure+6); hznc->SetLineColor(kAzure+6);
	406	hznc->SetMarkerStyle(21);
	407	hznc->Draw("");
	408	hzna->SetMarkerColor(kPink-2); hzna->SetLineColor(kPink-2);
	409	hzna->SetMarkerStyle(20);
	410	hzna->Draw("SAME");
	411	//
	412	TLegend *l1 = new TLegend(0.44,0.18,0.54,0.32);
	413	l1->SetFillColor(kWhite);
	414	l1->AddEntry(hznc," ZNC " ,"P");
	415	l1->AddEntry(hzna," ZNA " ,"P");
	416	l1->Draw("");
	417	//
	418	c1->cd(2);
	419	hznalg->SetMarkerColor(kBlue+1); hznalg->SetLineColor(kBlue+1);
	420	hznalg->SetMarkerStyle(20); hznalg->SetMinimum(0);
	421	hznalg->Draw("");
	422	hznclg->SetMarkerColor(kPink+5); hznclg->SetLineColor(kPink+5);
	423	hznclg->SetMarkerStyle(21);
	424	hznclg->Draw("SAME");
	425	//
	426	TLegend *l2 = new TLegend(0.44,0.18,0.54,0.32);
	427	l2->SetFillColor(kWhite);
	428	l2->AddEntry(hznc," ZNC LG " ,"P");
	429	l2->AddEntry(hzna," ZNA LG " ,"P");
	430	l2->Draw("");
	431
	432
	433	TCanvas *c1b = new TCanvas("c1b", "Mean value ZEMs ZPs", 200, 0, 1200, 1000);
	434	c1b->Divide(1,2);
	435	c1b->cd(1);
	436	hzpc->SetMarkerColor(kAzure+6); hzpc->SetLineColor(kAzure+6);
	437	hzpc->SetMarkerStyle(21);
	438	hzpc->Draw("");
	439	hzpa->SetMarkerColor(kPink-2); hzpa->SetLineColor(kPink-2);
	440	hzpa->SetMarkerStyle(20);
	441	hzpa->Draw("SAME");
	442	//
	443	TLegend *l3 = new TLegend(0.44,0.18,0.54,0.32);
	444	l3->SetFillColor(kWhite);
	445	l3->AddEntry(hzpc," ZPC " ,"P");
	446	l3->AddEntry(hzpa," ZPA " ,"P");
	447	l3->Draw("");
	448
	449	c1b->cd(2);
	450	hzem1->SetMarkerColor(kTeal-7); hzem1->SetLineColor(kTeal-7);
	451	hzem1->SetMarkerStyle(29);
	452	hzem2->SetMarkerColor(kTeal+5); hzem2->SetLineColor(kTeal+5);
	453	hzem2->SetMarkerStyle(30);
	454	hzem2->Draw("");
	455	hzem1->Draw("SAME");
	456	//
	457	TLegend *l4 = new TLegend(0.44,0.18,0.54,0.32);
	458	l4->SetFillColor(kWhite);
	459	l4->AddEntry(hzem1," ZEM1 " ,"P");
	460	l4->AddEntry(hzem2," ZEM2 " ,"P");
	461	l4->Draw("");
	462
	463
	464	/TCanvas c2 = new TCanvas("c2", "ZN centroids", 400, 400, 1200, 1000);
465	c2->Divide(1,2);
466	//
467	c2->cd(1);
468	hxznc->SetMarkerColor(kAzure+5); hxznc->SetLineColor(kAzure+5);
469	hxznc->SetMarkerStyle(21);
470	hxzna->SetMarkerColor(kPink+5); hxzna->SetLineColor(kPink+5);
471	hxzna->SetMarkerStyle(20);
472	hxznc->Draw("");
473	hxzna->Draw("SAME");
474	c2->cd(2);
475	hyznc->SetMarkerColor(kAzure+5); hyznc->SetLineColor(kAzure+5);
476	hyznc->SetMarkerStyle(21);
477	hyzna->SetMarkerColor(kPink+5); hyzna->SetLineColor(kPink+5);
478	hyzna->SetMarkerStyle(20);
479	hyznc->Draw("");
480	hyzna->Draw("SAME");*/
481
482
483	TCanvas *c3 = new TCanvas("c3", "LG signals", 600, 0, 1200, 1000);
484	c3->Divide(1,2);
485	//
486	c3->cd(1);
487	hznalg->SetMarkerColor(kBlue+1); hznalg->SetLineColor(kBlue+1);
488	hznalg->SetMarkerStyle(20);
489	hznalg->Draw("");
490	hznclg->SetMarkerColor(kPink+5); hznclg->SetLineColor(kPink+5);
491	hznclg->SetMarkerStyle(21);
492	hznclg->Draw("SAME");
493	TLegend *l7 = new TLegend(0.44,0.18,0.54,0.32);
494	l7->SetFillColor(kWhite);
495	l7->AddEntry(hznclg," ZNC LG " ,"P");
496	l7->AddEntry(hznalg," ZNA LG " ,"P");
497	l7->Draw("");
498	//
499	c3->cd(2);
500	hpmcznalg->SetMarkerColor(kAzure+7); hpmcznalg->SetLineColor(kAzure+7);
501	hpmcznalg->SetMarkerStyle(20); hpmcznalg->SetMinimum(0);
502	hpmcznalg->Draw("");
503	hpmcznclg->SetMarkerColor(kPink+6); hpmcznclg->SetLineColor(kPink+6);
504	hpmcznclg->SetMarkerStyle(21);
505	hpmcznclg->Draw("SAME");
506	TLegend *l8 = new TLegend(0.44,0.15,0.58,0.32);
507	l8->SetFillColor(kWhite);
508	l8->AddEntry(hpmcznclg," ZNC LG PMC" ,"P");
509	l8->AddEntry(hpmcznalg," ZNA LG PMC" ,"P");
510	l8->Draw("");
511	}