PWGPP/ZDC/macros/DrawPerformanceZDCQAMatch.C

   1 #if !defined(__CINT__) || defined(__MAKECINT__)
   2 #include <stdio.h>
   3 #include <stdlib.h>
   4 #include <TROOT.h>
   5 #include <Riostream.h>
   6 #include <TClassTable.h>
   7 #include <TStyle.h>
   8 #include <TMath.h>
   9 #include <TFile.h>
  10 #include <TCanvas.h>
  11 #include <TH1.h>
  12 #include <TH2.h>
  13 #include <TProfile.h>
  14 #include <TLine.h>
  15 #include <TGrid.h>
  16 #include <TBits.h>
  17 #include <TChain.h>
  18 #include <TNtuple.h>
  19 #include <TTree.h>
  20 #include <TBranch.h>
  21 #include <TFileMerger.h>
  22 #include <TGridResult.h>
  23 #include <TSystem.h>
  24 #include <TGaxis.h>
  25 #include <TRandom.h>
  26 #include <TLegend.h>
  27 #endif
  28
  29 Int_t DrawPerformanceZDCQAMatch(const char* inFile){
  30
  31 // Draw control histograms and generate output pictures
  32
  33 gSystem->Load("libSTAT");
  34 gSystem->Load("libANALYSIS");
  35 gSystem->Load("libANALYSISalice");
  36 gSystem->Load("libANALYSIScalib");
  37 gSystem->Load("libCORRFW");
  38 gSystem->Load("libANALYSISalice");
  39 gSystem->Load("libANALYSIScalib");
  40 gSystem->Load("libTender");
  41 gSystem->Load("libPWGPP");
  42
  43 gROOT->Reset();
  44 gROOT->SetStyle("Plain");
  45 gStyle->SetPalette(1);
  46 gStyle->SetOptStat(0);
  47 gStyle->SetTitleSize(0.025);
  48 TH1::AddDirectory(kFALSE);
  49
  50 TFile * fin = TFile::Open(inFile);
  51 TTree * ttree = (TTree*) fin->Get("trending");  //in
  52 TTree * tree = new TTree("tree","tree");        //out (to be summed)
  53
  54 if (!ttree){
  55   Printf("Invalid trending tree.");
  56   return 2;
  57 }
  58
  59 Int_t nRuns=ttree->GetEntries();
  60 TList list;
  61
  62 /*set graphic style*/
  63 gStyle->SetCanvasColor(kWhite);
  64 gStyle->SetFrameFillColor(kWhite);
  65 gStyle->SetFrameBorderMode(0);
  66 gStyle->SetCanvasBorderMode(0);
  67 gStyle->SetTitleFillColor(kWhite);
  68 gStyle->SetTitleBorderSize(0);
  69 gStyle->SetTitleFont(42);
  70 gStyle->SetTitleX(0.5);
  71 gStyle->SetTitleAlign(23);
  72 gStyle->SetTextFont(42);
  73 gStyle->SetStatColor(kWhite);
  74 gStyle->SetStatBorderSize(1);
  75 gStyle->SetOptStat(0);
  76 gStyle->SetTickLength(0.02,"y");
  77 gStyle->SetLabelSize(0.02,"xyz");
  78 gStyle->SetLabelOffset(0.03,"xyz");
  79
  80 char  outfilename[200]= "trending.root";
  81 TString plotDir(".");
  82
  83 legend = new TLegend(0.9,0.1,1.0,0.9);
  84 legend->SetFillColor(kWhite);
  85
  86 Int_t runNumber=0;
  87 Double_t ZNC_mean=0;
  88 Double_t ZNA_mean=0;
  89 Double_t ZPA_mean=0;
  90 Double_t ZPC_mean=0;
  91 Double_t ZEM1_mean=0;
  92 Double_t ZEM2_mean=0;
  93 Double_t ZNC_XCent=0;
  94 Double_t ZNC_YCent=0;
  95 Double_t ZNA_XCent=0;
  96 Double_t ZNA_YCent=0;
  97 Double_t ZNC_XCent_err=0;
  98 Double_t ZNC_YCent_err=0;
  99 Double_t ZNA_XCent_err=0;
 100 Double_t ZNA_YCent_err=0;
 101 Double_t ZN_TDC_Sum=0;
 102 Double_t ZN_TDC_Diff=0;
 103 Double_t ZN_TDC_Sum_err=0;
 104 Double_t ZN_TDC_Diff_err=0;
 105
 106 TH1D *hZNCpmcSpectrum = new TH1D("hZNCpmcSpectrum","hZNCpmcSpectrum",200.,0.,2000.);
 107 TH1D *hZNApmcSpectrum = new TH1D("hZNApmcSpectrum","hZNApmcSpectrum",200.,0.,2000.);
 108 TH1D *hZPCpmcSpectrum = new TH1D("hZPCpmcSpectrum","hZPCpmcSpectrum",200.,0.,2000.);
 109 TH1D *hZPApmcSpectrum = new TH1D("hZPApmcSpectrum","hZPApmcSpectrum",200.,0.,2000.);
 110 TH1D *hZEM1Spectrum = new TH1D("hZEM1Spectrum","hZEM1Spectrum",200.,0.,2000.);
 111 TH1D *hZEM2Spectrum = new TH1D("hZEM2Spectrum","hZEM2Spectrum",200.,0.,2000.);
 112
 113 ttree->SetBranchAddress("run",&runNumber);
 114 ttree->SetBranchAddress("ZNC_mean_value",&ZNC_mean);
 115 ttree->SetBranchAddress("ZNA_mean_value",&ZNA_mean);
 116 ttree->SetBranchAddress("ZPC_mean_value",&ZPC_mean);
 117 ttree->SetBranchAddress("ZPA_mean_value",&ZPA_mean);
 118 ttree->SetBranchAddress("ZEM1_mean_value",&ZEM1_mean);
 119 ttree->SetBranchAddress("ZEM2_mean_value",&ZEM2_mean);
 120 ttree->SetBranchAddress("ZNC_X_Centroid",&ZNC_XCent);
 121 ttree->SetBranchAddress("ZNC_Y_Centroid",&ZNC_YCent);
 122 ttree->SetBranchAddress("ZNA_X_Centroid",&ZNA_XCent);
 123 ttree->SetBranchAddress("ZNA_Y_Centroid",&ZNA_YCent);
 124 ttree->SetBranchAddress("ZNC_X_Centroid_Err",&ZNC_XCent_err);
 125 ttree->SetBranchAddress("ZNC_Y_Centroid_Err",&ZNC_YCent_err);
 126 ttree->SetBranchAddress("ZNA_X_Centroid_Err",&ZNA_XCent_err);
 127 ttree->SetBranchAddress("ZNA_Y_Centroid_Err",&ZNA_YCent_err);
 128 ttree->SetBranchAddress("ZN_TDC_Sum",&ZN_TDC_Sum);
 129 ttree->SetBranchAddress("ZN_TDC_Diff",&ZN_TDC_Diff);
 130 ttree->SetBranchAddress("ZN_TDC_Sum_Err",&ZN_TDC_Sum_err);
 131 ttree->SetBranchAddress("ZN_TDC_Diff_Err",&ZN_TDC_Diff_err);
 132
 133 TH1D *hznc = new TH1D("hznc","ZNC TC mean signal",3,-1,1);
 134 hznc->GetXaxis()->SetRangeUser(-1.,1.);
 135 hznc->SetDrawOption("EP");
 136 hznc->SetMarkerStyle(20);
 137 hznc->SetMarkerColor(kRed);
 138 hznc->SetLineColor(kRed);
 139
 140 TH1D *hzna = new TH1D("hzna","ZNA TC mean signal",3,-1,1);
 141 hzna->GetXaxis()->SetRangeUser(-1.,1.);
 142 hzna->SetDrawOption("EP");
 143 hzna->SetMarkerStyle(20);
 144 hzna->SetMarkerColor(kRed);
 145 hzna->SetLineColor(kRed);
 146
 147 TH1D *hzpc = new TH1D("hzpc","ZPC TC mean signal",3,-1,1);
 148 hzpc->GetXaxis()->SetRangeUser(-1.,1.);
 149 hzpc->SetDrawOption("EP");
 150 hzpc->SetMarkerStyle(20);
 151 hzpc->SetMarkerColor(kRed);
 152 hzpc->SetLineColor(kRed);
 153
 154 TH1D *hzpa = new TH1D("hzpa","ZPA TC mean signal",3,-1,1);
 155 hzpa->GetXaxis()->SetRangeUser(-1.,1.);
 156 hzpa->SetDrawOption("EP");
 157 hzpa->SetMarkerStyle(20);
 158 hzpa->SetMarkerColor(kRed);
 159 hzpa->SetLineColor(kRed);
 160
 161 TH1D *hzem1 = new TH1D("hzem1","ZEM1 mean signal",3,-1,1);
 162 hzem1->GetXaxis()->SetRangeUser(-1.,1.);
 163 hzem1->SetDrawOption("EP");
 164 hzem1->SetMarkerStyle(20);
 165 hzem1->SetMarkerColor(kRed);
 166 hzem1->SetLineColor(kRed);
 167
 168 TH1D *hzem2 = new TH1D("hzem2","ZEM2 mean signal",3,-1,1);
 169 hzem2->GetXaxis()->SetRangeUser(-1.,1.);
 170 hzem2->SetDrawOption("EP");
 171 hzem2->SetMarkerStyle(20);
 172 hzem2->SetMarkerColor(kRed);
 173 hzem2->SetLineColor(kRed);
 174
 175 TH1D *hzna_Xcentroid = new TH1D("hzna_Xcentroid","ZNA X centroid",3,-1,1);
 176 hzna_Xcentroid->GetXaxis()->SetRangeUser(-1.,1.);
 177 hzna_Xcentroid->SetDrawOption("EP");
 178 hzna_Xcentroid->SetMarkerStyle(20);
 179 hzna_Xcentroid->SetMarkerColor(kRed);
 180 hzna_Xcentroid->SetLineColor(kRed);
 181
 182 TH1D *hzna_Ycentroid = new TH1D("hzna_Ycentroid","ZNA Y centroid",3,-1,1);
 183 hzna_Ycentroid->GetXaxis()->SetRangeUser(-1.,1.);
 184 hzna_Ycentroid->SetDrawOption("EP");
 185 hzna_Ycentroid->SetMarkerStyle(20);
 186 hzna_Ycentroid->SetMarkerColor(kRed);
 187 hzna_Ycentroid->SetLineColor(kRed);
 188
 189 TH1D *hznc_Xcentroid = new TH1D("hznc_Xcentroid","ZNC X centroid",3,-1,1);
 190 hznc_Xcentroid->GetXaxis()->SetRangeUser(-1.,1.);
 191 hznc_Xcentroid->SetDrawOption("EP");
 192 hznc_Xcentroid->SetMarkerStyle(20);
 193 hznc_Xcentroid->SetMarkerColor(kRed);
 194 hznc_Xcentroid->SetLineColor(kRed);
 195
 196 TH1D *hznc_Ycentroid = new TH1D("hznc_Ycentroid","ZNC Y centroid",3,-1,1);
 197 hznc_Ycentroid->GetXaxis()->SetRangeUser(-1.,1.);
 198 hznc_Ycentroid->SetDrawOption("EP");
 199 hznc_Ycentroid->SetMarkerStyle(20);
 200 hznc_Ycentroid->SetMarkerColor(kRed);
 201 hznc_Ycentroid->SetLineColor(kRed);
 202
 203 TH1D *hzn_TDC_Sum = new TH1D("hzn_TDC_Sum","ZNC TDC + ZNA TDC",3,-1,1);
 204 hzn_TDC_Sum->GetXaxis()->SetRangeUser(-1.,1.);
 205 hzn_TDC_Sum->SetDrawOption("EP");
 206 hzn_TDC_Sum->SetMarkerStyle(20);
 207 hzn_TDC_Sum->SetMarkerColor(kRed);
 208 hzn_TDC_Sum->SetLineColor(kRed);
 209
 210 TH1D *hzn_TDC_Diff = new TH1D("hzn_TDC_Diff","ZNC TDC - ZNA TDC",3,-1,1);
 211 hzn_TDC_Diff->GetXaxis()->SetRangeUser(-1.,1.);
 212 hzn_TDC_Diff->SetDrawOption("EP");
 213 hzn_TDC_Diff->SetMarkerStyle(20);
 214 hzn_TDC_Diff->SetMarkerColor(kRed);
 215 hzn_TDC_Diff->SetLineColor(kRed);
 216
 217 char runlabel[10];
 218
 219 for (Int_t irun=0;irun<nRuns;irun++){
 220   ttree->GetEntry(irun);
 221  }
 222
 223 sprintf(runlabel,"%i",runNumber);
 224
 225 //----------------------------------------------------------------------
 226 //spectrum vs run
 227 //----------------------------------------------------------------------
 228
 229 hZNCpmcSpectrum = dynamic_cast<TH1D*> (fin->Get("fhZNCpmc"));
 230 hZNCpmcSpectrum->GetXaxis()->SetRangeUser(0.,2000);
 231 hZNCpmcSpectrum->Scale(1./hZNCpmcSpectrum->GetEntries());
 232 hZNCpmcSpectrum->SetLineColor(kRed);
 233 hZNCpmcSpectrum->SetLineWidth(2);
 234 hZNCpmcSpectrum->SetTitle("ZNC spectrum");
 235 hZNCpmcSpectrum->SetXTitle("ZNC signal (ADC ch.)");
 236
 237 hZNApmcSpectrum = dynamic_cast<TH1D*> (fin->Get("fhZNApmc"));
 238 hZNApmcSpectrum->GetXaxis()->SetRangeUser(0.,2000);
 239 hZNApmcSpectrum->Scale(1./hZNApmcSpectrum->GetEntries());
 240 hZNApmcSpectrum->SetLineColor(kRed);
 241 hZNApmcSpectrum->SetLineWidth(2);
 242 hZNApmcSpectrum->SetTitle("ZNA spectrum");
 243 hZNApmcSpectrum->SetXTitle("ZNA signal (ADC ch.)");
 244
 245 hZPCpmcSpectrum = dynamic_cast<TH1D*> (fin->Get("fhZPCpmc"));
 246 hZPCpmcSpectrum->GetXaxis()->SetRangeUser(0.,2000);
 247 hZPCpmcSpectrum->Scale(1./hZPCpmcSpectrum->GetEntries());
 248 hZPCpmcSpectrum->SetLineColor(kRed);
 249 hZPCpmcSpectrum->SetLineWidth(2);
 250 hZPCpmcSpectrum->SetTitle("ZPC spectrum");
 251 hZPCpmcSpectrum->SetXTitle("ZPC signal (ADC ch.)");
 252
 253 hZPApmcSpectrum = dynamic_cast<TH1D*> (fin->Get("fhZPApmc"));
 254 hZPApmcSpectrum->GetXaxis()->SetRangeUser(0.,2000);
 255 hZPApmcSpectrum->Scale(1./hZPApmcSpectrum->GetEntries());
 256 hZPApmcSpectrum->SetLineColor(kRed);
 257 hZPApmcSpectrum->SetLineWidth(2);
 258 hZPApmcSpectrum->SetTitle("ZPA spectrum");
 259 hZPApmcSpectrum->SetXTitle("ZPA signal (ADC ch.)");
 260
 261 hZEM1Spectrum = dynamic_cast<TH1D*> (fin->Get("fhZEM1Spectrum"));
 262 hZEM1Spectrum->GetXaxis()->SetRangeUser(0.,2000);
 263 hZEM1Spectrum->Scale(1./hZEM1Spectrum->GetEntries());
 264 hZEM1Spectrum->SetLineColor(kRed);
 265 hZEM1Spectrum->SetLineWidth(2);
 266 hZEM1Spectrum->SetTitle("ZEM1 spectrum");
 267 hZEM1Spectrum->SetXTitle("ZEM1 signal (ADC ch.)");
 268
 269 hZEM2Spectrum = dynamic_cast<TH1D*> (fin->Get("fhZEM2Spectrum"));
 270 hZEM2Spectrum->GetXaxis()->SetRangeUser(0.,2000);
 271 hZEM2Spectrum->Scale(1./hZEM2Spectrum->GetEntries());
 272 hZEM2Spectrum->SetLineColor(kRed);
 273 hZEM2Spectrum->SetLineWidth(2);
 274 hZEM2Spectrum->SetTitle("ZEM2 spectrum");
 275 hZEM2Spectrum->SetXTitle("ZEM2 signal (ADC ch.)");
 276
 277 //----------------------------------------------------------------------
 278 //variables vs run
 279 //----------------------------------------------------------------------
 280 hzna->SetBinContent(2,ZNA_mean);
 281 hzna->GetXaxis()->SetBinLabel(2,runlabel);
 282 hzna->GetXaxis()->SetLabelSize(0.05);
 283
 284 hzpc->SetBinContent(2,ZPC_mean);
 285 hzpc->GetXaxis()->SetBinLabel(2,runlabel);
 286 hzpc->GetXaxis()->SetLabelSize(0.05);
 287
 288 hzpa->SetBinContent(2,ZPA_mean);
 289 hzpa->GetXaxis()->SetBinLabel(2,runlabel);
 290 hzpa->GetXaxis()->SetLabelSize(0.05);
 291
 292 hzem1->SetBinContent(2,ZEM1_mean);
 293 hzem1->GetXaxis()->SetBinLabel(2,runlabel);
 294 hzem1->GetXaxis()->SetLabelSize(0.05);
 295
 296 hzem2->SetBinContent(2,ZEM2_mean);
 297 hzem2->GetXaxis()->SetBinLabel(2,runlabel);
 298 hzem2->GetXaxis()->SetLabelSize(0.05);
 299
 300 hzna_Xcentroid->SetBinContent(2,ZNA_XCent);
 301 hzna_Xcentroid->SetBinError(2,ZNA_XCent_err);
 302 hzna_Xcentroid->GetXaxis()->SetBinLabel(2,runlabel);
 303 hzna_Xcentroid->GetXaxis()->SetLabelSize(0.05);
 304 hzna_Xcentroid->GetYaxis()->SetTitle("(cm)");
 305
 306 hzna_Ycentroid->SetBinContent(2,ZNA_YCent);
 307 hzna_Ycentroid->SetBinError(2,ZNA_YCent_err);
 308 hzna_Ycentroid->GetXaxis()->SetBinLabel(2,runlabel);
 309 hzna_Ycentroid->GetXaxis()->SetLabelSize(0.05);
 310 hzna_Ycentroid->GetYaxis()->SetTitle("(cm)");
 311
 312 hznc_Xcentroid->SetBinContent(2,ZNC_XCent);
 313 hznc_Xcentroid->SetBinError(2,ZNC_XCent_err);
 314 hznc_Xcentroid->GetXaxis()->SetBinLabel(2,runlabel);
 315 hznc_Xcentroid->GetXaxis()->SetLabelSize(0.05);
 316 hznc_Xcentroid->GetYaxis()->SetTitle("(cm)");
 317
 318 hznc_Ycentroid->SetBinContent(2,ZNC_YCent);
 319 hznc_Ycentroid->SetBinError(2,ZNC_YCent_err);
 320 hznc_Ycentroid->GetXaxis()->SetBinLabel(2,runlabel);
 321 hznc_Ycentroid->GetXaxis()->SetLabelSize(0.05);
 322 hznc_Ycentroid->GetYaxis()->SetTitle("(cm)");
 323
 324 hzn_TDC_Sum->SetBinContent(2,ZN_TDC_Sum);
 325 hzn_TDC_Sum->SetBinError(2,ZN_TDC_Sum_err);
 326 hzn_TDC_Sum->GetXaxis()->SetBinLabel(2,runlabel);
 327 hzn_TDC_Sum->GetXaxis()->SetLabelSize(0.05);
 328 hzn_TDC_Sum->GetYaxis()->SetTitle("(ns)");
 329
 330 hzn_TDC_Diff->SetBinContent(2,ZN_TDC_Diff);
 331 hzn_TDC_Diff->SetBinError(2,ZN_TDC_Diff_err);
 332 hzn_TDC_Diff->GetXaxis()->SetBinLabel(2,runlabel);
 333 hzn_TDC_Diff->GetXaxis()->SetLabelSize(0.05);
 334 hzn_TDC_Diff->GetYaxis()->SetTitle("(ns)");
 335
 336 //----------------------------------------------------------------------
 337 //spectra
 338 //----------------------------------------------------------------------
 339
 340 TCanvas* cZNC_Spectra = new TCanvas("cZNC_Spectra","cZNC_Spectra",50,50,1200,900); gPad->SetLogy();
 341 hZNCpmcSpectrum->Draw();
 342 cZNC_Spectra->Print(Form("%s/cZNC_Spectra.png",plotDir.Data()));
 343
 344 TCanvas* cZNA_Spectra = new TCanvas("cZNA_Spectra","cZNA_Spectra",50,50,1200,900); gPad->SetLogy();
 345 hZNApmcSpectrum->Draw();
 346 cZNA_Spectra->Print(Form("%s/cZNA_Spectra.png",plotDir.Data()));
 347
 348 TCanvas* cZPC_Spectra = new TCanvas("cZPC_Spectra","cZPC_Spectra",50,50,1200,900); gPad->SetLogy();
 349 hZPCpmcSpectrum->Draw();
 350 cZPC_Spectra->Print(Form("%s/cZPC_Spectra.png",plotDir.Data()));
 351
 352 TCanvas* cZPA_Spectra = new TCanvas("cZPA_Spectra","cZPA_Spectra",50,50,1200,900); gPad->SetLogy();
 353 hZPApmcSpectrum->Draw();
 354 cZPA_Spectra->Print(Form("%s/cZPA_Spectra.png",plotDir.Data()));
 355
 356 TCanvas* cZEM1_Spectra = new TCanvas("cZEM1_Spectra","cZEM1_Spectra",50,50,1200,900); gPad->SetLogy();
 357 hZEM1Spectrum->Draw();
 358 cZEM1_Spectra->Print(Form("%s/cZEM1_Spectra.png",plotDir.Data()));
 359
 360 TCanvas* cZEM2_Spectra = new TCanvas("cZEM2_Spectra","cZEM2_Spectra",50,50,1200,900); gPad->SetLogy();
 361 hZEM2Spectrum->Draw();
 362 cZEM2_Spectra->Print(Form("%s/cZEM2_Spectra.png",plotDir.Data()));
 363
 364
 365 //---------------------------------------------------------------------------------------------------
 366 //variables
 367 //---------------------------------------------------------------------------------------------------
 368
 369 //---------------------------------------------------------------------------------------------------
 370 //means
 371 //---------------------------------------------------------------------------------------------------
 372 TCanvas* cZNC_Mean_Values = new TCanvas("cZNC_Mean_Values","cZNC_Mean_Values", 50,50,750,550);
 373 hznc->Draw("ep");
 374 hznc->SetBinContent(2,ZNC_mean);
 375 hznc->GetXaxis()->SetBinLabel(2,runlabel);
 376 hznc->GetXaxis()->SetLabelSize(0.05);
 377 cZNC_Mean_Values->Print(Form("%s/cZNC_Mean_Values.png",plotDir.Data()));
 378
 379 TCanvas* cZNA_Mean_Values = new TCanvas("cZNA_Mean_Values","cZNA_Mean_Values", 50,50,750,550);
 380 hzna->Draw("ep");
 381 cZNA_Mean_Values->Print(Form("%s/cZNA_Mean_Values.png",plotDir.Data()));
 382
 383 TCanvas* cZPC_Mean_Values = new TCanvas("cZPC_Mean_Values","cZPC_Mean_Values", 50,50,750,550);
 384 hzpc->Draw("ep");
 385 cZPC_Mean_Values->Print(Form("%s/cZPC_Mean_Values.png",plotDir.Data()));
 386
 387 TCanvas* cZPA_Mean_Values = new TCanvas("cZPA_Mean_Values","cZPA_Mean_Values", 50,50,750,550);
 388 hzpa->Draw("ep");
 389 cZPA_Mean_Values->Print(Form("%s/cZPA_Mean_Values.png",plotDir.Data()));
 390
 391 TCanvas* cZEM1_Mean_Values = new TCanvas("cZEM1_Mean_Values","cZEM1_Mean_Values", 50,50,750,550);
 392 hzem1->Draw("ep");
 393 cZEM1_Mean_Values->Print(Form("%s/cZEM1_Mean_Values.png",plotDir.Data()));
 394
 395 TCanvas* cZEM2_Mean_Values = new TCanvas("cZEM2_Mean_Values","cZEM2_Mean_Values", 50,50,750,550);
 396 hzem2->Draw("ep");
 397 cZEM2_Mean_Values->Print(Form("%s/cZEM2_Mean_Values.png",plotDir.Data()));
 398
 399 //---------------------------------------------------------------------------------------------------
 400 //centroids
 401 //---------------------------------------------------------------------------------------------------
 402 TCanvas* cZNA_X_centroid = new TCanvas("cZNA_X_centroid","cZNA_X_centroid", 50,50,750,550);
 403 hzna_Xcentroid->Draw();
 404 cZNA_X_centroid->Print(Form("%s/cZNA_X_centroid.png",plotDir.Data()));
 405
 406 TCanvas* cZNA_Y_centroid = new TCanvas("cZNA_Y_centroid","cZNA_Y_centroid", 50,50,750,550);
 407 hzna_Ycentroid->Draw();
 408 cZNA_Y_centroid->Print(Form("%s/cZNA_Y_centroid.png",plotDir.Data()));
 409
 410 TCanvas* cZNC_X_centroid = new TCanvas("cZNC_X_centroid","cZNC_X_centroid", 50,50,750,550);
 411 hznc_Xcentroid->Draw();
 412 cZNC_X_centroid->Print(Form("%s/cZNC_X_centroid.png",plotDir.Data()));
 413
 414 TCanvas* cZNC_Y_centroid = new TCanvas("cZNC_Y_centroid","cZNC_Y_centroid", 50,50,750,550);
 415 hznc_Ycentroid->Draw();
 416 cZNC_Y_centroid->Print(Form("%s/cZNC_Y_centroid.png",plotDir.Data()));
 417
 418 //---------------------------------------------------------------------------------
 419 //timing
 420 //---------------------------------------------------------------------------------
 421 TCanvas* cTimingSum = new TCanvas("cTimingSum","cTimingSum",50,50,750,550);
 422 hzn_TDC_Sum->Draw();
 423 cTimingSum->Print(Form("%s/cTimingSum.png",plotDir.Data()));
 424
 425 TCanvas* cTimingDiff = new TCanvas("cTimingDiff","cTimingDiff",50,50,750,550);
 426 hzn_TDC_Diff->Draw();
 427 cTimingDiff->Print(Form("%s/cTimingDiff.png",plotDir.Data()));
 428
 429 //----------------------------------------------------------------------
 430 //out
 431 //----------------------------------------------------------------------
 432 tree->Branch("run",&runNumber,"runNumber/I");
 433 tree->Branch("ZNC_mean_value",&ZNC_mean,"ZNC_mean/D");
 434 tree->Branch("ZNA_mean_value",&ZNA_mean,"ZNA_mean/D");
 435 tree->Branch("ZPC_mean_value",&ZPC_mean,"ZPC_mean/D");
 436 tree->Branch("ZPA_mean_value",&ZPA_mean,"ZPA_mean/D");
 437 tree->Branch("ZEM1_mean_value",&ZEM1_mean,"ZEM1_mean/D");
 438 tree->Branch("ZEM2_mean_value",&ZEM2_mean,"ZEM2_mean/D");
 439 tree->Branch("ZNC_X_Centroid",&ZNC_XCent,"ZNC_XCent/D");
 440 tree->Branch("ZNC_Y_Centroid",&ZNC_YCent,"ZNC_YCent/D");
 441 tree->Branch("ZNA_X_Centroid",&ZNA_XCent,"ZNA_XCent/D");
 442 tree->Branch("ZNA_Y_Centroid",&ZNA_YCent,"ZNA_YCent/D");
 443 tree->Branch("ZNC_X_Centroid_Err",&ZNC_XCent_err,"ZNC_XCent_err/D");
 444 tree->Branch("ZNC_Y_Centroid_Err",&ZNC_YCent_err,"ZNC_YCent_err/D");
 445 tree->Branch("ZNA_X_Centroid_Err",&ZNA_XCent_err,"ZNA_XCent_err/D");
 446 tree->Branch("ZNA_Y_Centroid_Err",&ZNA_YCent_err,"ZNA_YCent_err/D");
 447 tree->Branch("ZN_TDC_Sum",&ZN_TDC_Sum,"ZN_TDC_Sum/D");
 448 tree->Branch("ZN_TDC_Diff",&ZN_TDC_Diff,"ZN_TDC_Diff/D");
 449 tree->Branch("ZN_TDC_Sum_Err",&ZN_TDC_Sum_err,"ZN_TDC_Sum_err/D");
 450 tree->Branch("ZN_TDC_Diff_Err",&ZN_TDC_Diff_err,"ZN_TDC_Diff_err/D");
 451 tree->Fill();
 452
 453 list.Add(cZNC_Spectra);
 454 list.Add(cZNA_Spectra);
 455 list.Add(cZPC_Spectra);
 456 list.Add(cZPA_Spectra);
 457 list.Add(cZEM1_Spectra);
 458 list.Add(cZEM2_Spectra);
 459 list.Add(hznc);
 460 list.Add(hzna);
 461 list.Add(hzpc);
 462 list.Add(hzpa);
 463 list.Add(hzem1);
 464 list.Add(hzem2);
 465 list.Add(hzna_Xcentroid);
 466 list.Add(hzna_Ycentroid);
 467 list.Add(hznc_Xcentroid);
 468 list.Add(hznc_Ycentroid);
 469 list.Add(hzn_TDC_Sum);
 470 list.Add(hzn_TDC_Diff);
 471
 472 TFile * fout=new TFile(outfilename,"recreate");
 473 fout->cd();
 474 list.Write();
 475 tree.Write();
 476 fout->Close();
 477
 478 return 0;
 479
 480 }