Macro for a manual check of the automatic fit on pedestal data

diff --git a/ZDC/CheckPedestal.C b/ZDC/CheckPedestal.C
new file mode 100644 (file)
index 0000000..d1764ad
--- /dev/null
+++ b/ZDC/CheckPedestal.C
@@ -0,0 +1,196 @@
+#if !defined(__CINT__) || defined(__MAKECINT__)
+
+#include <TROOT.h>
+#include <TStyle.h>
+#include <Riostream.h>
+#include "TClassTable.h"
+#include <TStopwatch.h>
+#include <TDatime.h>
+#include <TClassTable.h>
+#include <TH1.h>
+#include <TH2.h>
+#include <TF1.h>
+#include <TProfile.h>
+#include <TFunction.h>
+#include <TCanvas.h>
+#include <TFile.h>
+
+#endif
+
+void CheckPedestal( Int_t optPlot = 1, const char * histoFName = "ZDCPedHisto.root")
+{
+
+  int const kNChannels = 24;
+  //
+  TFile * file = new TFile(histoFName, "READ");
+  file->cd();
+  TH1F::AddDirectory(0);
+  //
+  TH1F *hPedhg[kNChannels], *hPedOutOfTimehg[kNChannels];
+  TH2F *hPedCorrhg[kNChannels];
+  TH1F *hPedlg[kNChannels], *hPedOutOfTimelg[kNChannels];
+  TH2F *hPedCorrlg[kNChannels];
+  //
+  char namhist1hg[50], namhist2hg[50], namhist3hg[50];
+  char namhist1lg[50], namhist2lg[50], namhist3lg[50];
+  for(Int_t j=0; j<kNChannels; j++){
+     if(j<=4){ // ZNC
+       sprintf(namhist1hg,"PedZNChg_%d",j);
+       sprintf(namhist2hg,"PedZNChgOutOfTime_%d",j);
+       sprintf(namhist3hg,"PedCorrZNChg_%d",j);
+       //
+       sprintf(namhist1lg,"PedZNClg_%d",j);
+       sprintf(namhist2lg,"PedZNClgOutOfTime_%d",j);
+       sprintf(namhist3lg,"PedCorrZNClg_%d",j);
+     }
+     else if(j>=5 && j<=9){ // ZPC
+       sprintf(namhist1hg,"PedZPChg_%d",j-5);
+       sprintf(namhist2hg,"PedZPChgOutOfTime_%d",j-5);
+       sprintf(namhist3hg,"PedCorrZPChg_%d",j-5);
+       //
+       sprintf(namhist1lg,"PedZPClg_%d",j-5);
+       sprintf(namhist2lg,"PedZPClgOutOfTime_%d",j-5);
+       sprintf(namhist3lg,"PedCorrZPClg_%d",j-5);       
+     }
+     else if(j==10 || j==11){ // ZEM
+       sprintf(namhist1hg,"PedZEMhg_%d",j-10);
+       sprintf(namhist2hg,"PedZEMhgOutOfTime_%d",j-10);
+       sprintf(namhist3hg,"PedCorrZEMhg_%d",j-10);
+       //
+       sprintf(namhist1lg,"PedZEMlg_%d",j-10);
+       sprintf(namhist2lg,"PedZEMlgOutOfTime_%d",j-10);
+       sprintf(namhist3lg,"PedCorrZEMlg_%d",j-10);
+     }
+     else if(j>=12 && j<=16){ // ZNA
+       sprintf(namhist1hg,"PedZNAhg_%d",j-12);
+       sprintf(namhist2hg,"PedZNAhgOutOfTime_%d",j-12);
+       sprintf(namhist3hg,"PedCorrZNAhg_%d",j-12);
+       //
+       sprintf(namhist1lg,"PedZNAlg_%d",j-12);
+       sprintf(namhist2lg,"PedZNAlgOutOfTime_%d",j-12);
+       sprintf(namhist3lg,"PedCorrZNAlg_%d",j-12);
+     }
+     else if(j>=17 && j<=21){ // ZPA
+       sprintf(namhist1hg,"PedZPAhg_%d",j-17);
+       sprintf(namhist2hg,"PedZPAhgOutOfTime_%d",j-17);
+       sprintf(namhist3hg,"PedCorrZPAhg_%d",j-17);
+       //
+       sprintf(namhist1lg,"PedZPAlg_%d",j-17);
+       sprintf(namhist2lg,"PedZPAlgOutOfTime_%d",j-17);
+       sprintf(namhist3lg,"PedCorrZPAlg_%d",j-17);
+     }
+     else if(j==22 || j==24){ //Reference PMs
+       sprintf(namhist1hg,"PedRefhg_%d",j-22);
+       sprintf(namhist2hg,"PedRefhgOutOfTime_%d",j-22);
+       sprintf(namhist3hg,"PedCorrRefhg_%d",j-22);
+       //
+       sprintf(namhist1lg,"PedReflg_%d",j-22);
+       sprintf(namhist2lg,"PedReflgOutOfTime_%d",j-22);
+       sprintf(namhist3lg,"PedCorrReflg_%d",j-22);
+     }
+     // --- High gain chain histos
+     hPedhg[j] = (TH1F*) file->Get(namhist1hg);
+     hPedOutOfTimehg[j] = (TH1F*) file->Get(namhist2hg);
+     hPedCorrhg[j] =  (TH2F*) file->Get(namhist3hg);
+     // --- Low gain chain histos
+     hPedlg[j] = (TH1F*) file->Get(namhist1lg);
+     hPedOutOfTimelg[j] =(TH1F*) file->Get(namhist2lg);
+     hPedCorrlg[j] = (TH2F*) file->Get(namhist3lg);
+  }
+
+ if(optPlot==1){
+  // Plot the retrieved histos
+  //***********************************************************
+  // #### ROOT initialization
+  gROOT->Reset();
+  gStyle->SetCanvasColor(10);
+  gStyle->SetFrameFillColor(10);
+  gStyle->SetOptTitle(0);
+  gStyle->SetOptStat(1111);
+  gStyle->SetOptFit(111);
+  gStyle->SetTitleTextColor(9);
+  gStyle->SetStatTextColor(4);
+  gStyle->SetLineColor(1);
+  gStyle->SetPalette(1);
+  //***********************************************************
+  TCanvas *c6 = new TCanvas("c6","Side C correlations",0,200,1000,800);
+  c6->Divide(5,4);
+  for(Int_t t=0; t<10; t++){
+    c6->cd(t+1);
+    hPedCorrhg[t]->Draw();
+    c6->cd(t+11);
+    hPedCorrlg[t]->Draw();
+  }
+  c6->Print("CorrSideC.ps");
+  //
+  TCanvas *c7 = new TCanvas("c7","Side A correlations",300,200,1000,800);
+  c7->Divide(5,4);
+  for(Int_t t=0; t<10; t++){
+    c7->cd(t+1);
+    hPedCorrhg[t+12]->Draw();
+    c7->cd(t+11);
+    hPedCorrlg[t+12]->Draw();
+  }
+  c7->Print("CorrSideA.ps");
+  //
+  TCanvas *c8 = new TCanvas("c8","ZEM correlations",400,200,400,400);
+  c8->Divide(2,2);
+  for(Int_t t=0; t<2; t++){
+    c8->cd(t+1);
+    hPedCorrhg[t+10]->Draw();
+    c8->cd(t+3);
+    hPedCorrlg[t+10]->Draw();
+  }
+  c8->Print("CorrZEM.ps");
+  //***********************************************************
+  TCanvas *c1 = new TCanvas("c1","ZNC pedestals",0,0,1000,400);
+  c1->Divide(5,2);
+  for(Int_t y=0; y<5; y++){
+    c1->cd(y+1);
+    hPedhg[y]->Draw();
+    c1->cd(y+6);
+    hPedlg[y]->Draw();
+  }
+  c1->Print("ZNCPed.ps");
+  //
+  TCanvas *c2 = new TCanvas("c2","ZPC pedestals",300,0,1000,400);
+  c2->Divide(5,2);
+  for(Int_t y=0; y<5; y++){
+    c2->cd(y+1);
+    hPedhg[y+5]->Draw();
+    c2->cd(y+6);
+    hPedlg[y+5]->Draw();
+  }
+  c2->Print("ZPCPed.ps");
+  //
+  TCanvas *c3 = new TCanvas("c3","ZEM pedestals",400,0,400,400);
+  c3->Divide(2,2);
+  for(Int_t y=0; y<2; y++){
+    c3->cd(y+1);
+    hPedhg[y+10]->Draw();
+    c3->cd(y+3);
+    hPedlg[y+10]->Draw();
+  }
+  c3->Print("ZEMPed.ps");
+  //
+  TCanvas *c4 = new TCanvas("c4","ZNA pedestals",0,400,1000,400);
+  c4->Divide(5,2);
+  for(Int_t y=0; y<5; y++){
+    c4->cd(y+1);
+    hPedhg[y+12]->Draw();
+    c4->cd(y+6);
+    hPedlg[y+12]->Draw();
+  }
+  c4->Print("ZNAPed.ps");
+  //
+  TCanvas *c5 = new TCanvas("c5","ZPA pedestals",300,400,1000,400);
+  c5->Divide(5,2);
+  for(Int_t y=0; y<5; y++){
+    c5->cd(y+1);
+    hPedhg[y+17]->Draw();
+    c5->cd(y+6);
+    hPedlg[y+17]->Draw();
+  }
+  c5->Print("ZPAPed.ps");
+ }
+}