+ char outfilename[200];
+ sprintf(outfilename, "trend_%s",trendFileName);
+ TFile * fout=new TFile(outfilename,"recreate");
+ fout->cd();
+ lista.Write();
+ fout->Close();
+
+ TString plotDir(Form("PlotsTrending"));
+ gSystem->Exec(Form("mkdir %s",plotDir.Data()));
+
+ TCanvas* cAvDiffTimeVsRun = new TCanvas("cAvDiffTimeVsRun","cAvDiffTimeVsRun",50,50,750,550);
+ gPad->SetGridx();
+ gPad->SetGridy();
+ hAvDiffTimeVsRun->Draw();
+ cAvDiffTimeVsRun->Print(Form("%s/cAvDiffTimeVsRun.png",plotDir.Data()));
+
+ TCanvas* cPeakDiffTimeVsRun = new TCanvas("cPeakDiffTimeVsRun","cPeakDiffTimeVsRun", 50,50,750,550);
+ hPeakDiffTimeVsRun->Draw();
+ cPeakDiffTimeVsRun->Print(Form("%s/cPeakDiffTimeVsRun",plotDir.Data()));
+
+ TCanvas* cSpreadDiffTimeVsRun = new TCanvas("cSpreadDiffTimeVsRun","cSpreadDiffTimeVsRun", 50,50,750,550);
+ hSpreadDiffTimeVsRun->Draw();
+ cSpreadDiffTimeVsRun->Print(Form("%s/cSpreadDiffTimeVsRun",plotDir.Data()));
+
+ TCanvas* cAvTimeVsRun = new TCanvas("cAvTimeVsRun","cAvTimeVsRun", 50,50,750,550);
+ hAvTimeVsRun->Draw();
+ cAvTimeVsRun->Print(Form("%s/cAvTimeVsRun",plotDir.Data()));
+
+ TCanvas* cPeakTimeVsRun = new TCanvas("cPeakTimeVsRun","cPeakTimeVsRun", 50,50,750,550);
+ hPeakTimeVsRun->Draw();
+ cPeakTimeVsRun->Print(Form("%s/cPeakTimeVsRun",plotDir.Data()));
+
+ TCanvas* cSpreadTimeVsRun = new TCanvas("cSpreadTimeVsRun","cSpreadTimeVsRun", 50,50,750,550);
+ hSpreadTimeVsRun->Draw();
+ cSpreadTimeVsRun->Print(Form("%s/cSpreadTimeVsRun",plotDir.Data()));
+
+ TCanvas* cAvRawTimeVsRun = new TCanvas("cAvRawTimeVsRun","cAvRawTimeVsRun", 50,50,750,550);
+ hAvRawTimeVsRun->Draw();
+ cAvRawTimeVsRun->Print(Form("%s/cAvRawTimeVsRun",plotDir.Data()));
+
+ TCanvas* cPeakRawTimeVsRun = new TCanvas("cPeakRawTimeVsRun","cPeakRawTimeVsRun", 50,50,750,550);
+ hPeakRawTimeVsRun->Draw();
+ cPeakRawTimeVsRun->Print(Form("%s/cPeakRawTimeVsRun",plotDir.Data()));
+
+ TCanvas* cSpreadRawTimeVsRun = new TCanvas("cSpreadRawTimeVsRun","cSpreadRawTimeVsRun", 50,50,750,550);
+ hSpreadRawTimeVsRun->Draw();
+ cSpreadRawTimeVsRun->Print(Form("%s/cSpreadRawTimeVsRun",plotDir.Data()));
+
+ TCanvas* cAvTotVsRun = new TCanvas("cAvTotVsRun","cAvTotVsRun", 50,50,750,550);
+ hAvTotVsRun->Draw();
+ cAvTotVsRun->Print(Form("%s/cAvTotVsRun",plotDir.Data()));
+
+ TCanvas* cPeakTotVsRun = new TCanvas("cPeakTotVsRun","cPeakTotVsRun", 50,50,750,550);
+ hPeakTotVsRun->Draw();
+ cPeakTotVsRun->Print(Form("%s/cPeakTotVsRun",plotDir.Data()));
+
+ TCanvas* cSpreadTotVsRun = new TCanvas("cSpreadTotVsRun","cSpreadTotVsRun", 50,50,750,550);
+ hSpreadTotVsRun->Draw();
+ cSpreadTotVsRun->Print(Form("%s/cSpreadTotVsRun",plotDir.Data()));
+
+ TCanvas* cNegTimeRatioVsRun = new TCanvas("cNegTimeRatioVsRun","cNegTimeRatioVsRun", 50,50,750,550);
+ hNegTimeRatioVsRun->Draw();
+ cNegTimeRatioVsRun->Print(Form("%s/cNegTimeRatioVsRun",plotDir.Data()));
+
+ TCanvas* cOrphansRatioVsRun = new TCanvas("cOrphansRatioVsRun","cOrphansRatioVsRun", 50,50,750,550);
+ hOrphansRatioVsRun->Draw();
+ cOrphansRatioVsRun->Print(Form("%s/cOrphansRatioVsRun",plotDir.Data()));
+
+ TCanvas* cMeanLVsRun = new TCanvas("cMeanLVsRun","cMeanLVsRun", 50,50,750,550);
+ hMeanLVsRun->Draw();
+ cMeanLVsRun->Print(Form("%s/cMeanLVsRun",plotDir.Data()));
+
+ TCanvas* cNegLRatioVsRun = new TCanvas("cNegLRatioVsRun","cNegLRatioVsRun", 50,50,750,550);
+ hNegLRatioVsRun->Draw();
+ cNegLRatioVsRun->Print(Form("%s/cNegLRatioVsRun",plotDir.Data()));
+
+ TCanvas* cMatchEffVsRun = new TCanvas("cMatchEffVsRun","cMatchEffVsRun", 50,50,750,550);
+ hMatchEffVsRun->Draw();
+ cMatchEffVsRun->Print(Form("%s/cMatchEffVsRun",plotDir.Data()));
+
+ TCanvas* cPeakT0AVsRun = new TCanvas("cPeakT0AVsRun","cPeakT0AVsRun", 50,50,750,550);
+ hPeakT0AVsRun->Draw();
+ cPeakT0AVsRun->Print(Form("%s/cPeakT0AVsRun.png",plotDir.Data()));
+
+ TCanvas* cPeakT0CVsRun = new TCanvas("cPeakT0CVsRun","cPeakT0CVsRun", 50,50,750,550);
+ hPeakT0CVsRun->Draw();
+ cPeakT0CVsRun->Print(Form("%s/cPeakT0CVsRun.png",plotDir.Data()));
+
+ TCanvas* cPeakT0ACVsRun = new TCanvas("cPeakT0ACVsRun","cPeakT0ACVsRun", 50,50,750,550);
+ hPeakT0ACVsRun->Draw();
+ cPeakT0ACVsRun->Print(Form("%s/cPeakT0ACVsRun.png",plotDir.Data()));
+
+ TCanvas* cT0fillResVsRun = new TCanvas("cT0fillResVsRun","cT0fillResVsRun", 50,50,750,550);
+ hT0fillResVsRun->Draw();
+ cT0fillResVsRun->Print(Form("%s/cT0fillResVsRun.png",plotDir.Data()));
+
+
+ return 0;