--- quicklook/trigger/src/TriggerScan.cpp 2006/05/11 10:07:39 1.1 +++ quicklook/trigger/src/TriggerScan.cpp 2006/06/26 13:28:11 1.2 @@ -1,9 +1,11 @@ -/** +/* * TriggerScan * Author Nagni * Modified by G. De Rosa * Version 1.2 * Date 27 Apr 2006 + * Version 1.4 + * Date 26 Jun 2006: Inserted dependences on OrbitalTime and modified histogram ranges. * * Description: * Describe the performance of the Trigger. @@ -38,7 +40,6 @@ } void TriggerScan(TString base, TString outDirectory = "", TString format = ""){ - TLatex *t=new TLatex(); std::stringstream sst; if (outDirectory == "") outDirectory = base.Data(); @@ -70,59 +71,7 @@ } /* Histo DEF */ - - /* - * PTM Plane Counters histograms - */ - TCanvas *PMTCanvas = new TCanvas("PMTCanvas","PMTCanvas", 1280, 1024); - PMTCanvas->Divide(1,3); - //---------- - sst.str(""); - sst << "Histogram of PMT Plane Counters (S11 AND S12) based on " << base.Data(); - TH1S *pmt1Histo = new TH1S("PMT1Histo", sst.str().c_str(), 20, 0., 300.); - //---------- - sst.str(""); - sst << "Histogram of PMT Plane Counters (S21 AND S22) based on " << base.Data(); - TH1S *pmt2Histo = new TH1S("PMT2Histo", sst.str().c_str(), 20, 0., 300.); - //---------- - sst.str(""); - sst << "Histogram of PMT Plane Counters (S31 AND S32) based on " << base.Data(); - TH1S *pmt3Histo = new TH1S("PMT3Histo", sst.str().c_str(), 20, 0., 300.); - - /* - * Trigger Rate Counters histograms - */ - TCanvas *TrigRate1Canvas = new TCanvas("TrigRate1Canvas","TrigRate1Canvas", 1280, 1024); - TrigRate1Canvas->Divide(1,3); - - TCanvas *TrigRate2Canvas = new TCanvas("TrigRate2Canvas","TrigRate2Canvas", 1280, 1024); - TrigRate2Canvas->Divide(1,3); - //---------- - sst.str(""); - sst << "Histogram of Trigger Rate Counters (S11 OR S12) AND (S21 OR S22) AND (S31 OR S32) based on " << base.Data(); - TH1S *trRate1Histo = new TH1S("trRate1Histo", sst.str().c_str(), 20, 0., 5.); - //---------- - sst.str(""); - sst << "Histogram of Trigger Rate Counters (S11 AND S12) AND (S21 AND S22) AND (S31 AND S32) based on " << base.Data(); - TH1S *trRate2Histo = new TH1S("trRate2Histo", sst.str().c_str(), 20, 0., 5.); - //--------- - sst.str(""); - sst << "Histogram of Trigger Rate Counters (S21 OR S22) AND (S31 OR S32) based on " << base.Data(); - TH1S *trRate3Histo = new TH1S("trRate3Histo", sst.str().c_str(), 20, 0., 5.); - //--------- - sst.str(""); - sst << "Histogram of Trigger Rate Counters (S21 AND S22) AND (S31 AND S32) based on " << base.Data(); - TH1S *trRate4Histo = new TH1S("trRate4Histo", sst.str().c_str(), 20, 0., 5.); - //--------- - sst.str(""); - sst << "Histogram of Trigger Rate Counters S12 AND (S21 OR S22) based on " << base.Data(); - TH1S *trRate5Histo = new TH1S("trRate5Histo", sst.str().c_str(), 20, 0., 5.); - //--------- - sst.str(""); - sst << "Histogram of Trigger Rate Counters (S11 OR S12) AND (S31 OR S32) based on " << base.Data(); - TH1S *trRate6Histo = new TH1S("trRate6Histo", sst.str().c_str(), 20, 0., 5.); - - + /* * S4 and Calorimeter Counters histograms */ @@ -131,7 +80,7 @@ //---------- sst.str(""); sst << "S4 Counters based on " << base.Data(); - TH1S *S4Cal1Histo = new TH1S("S4Cal1Histo", sst.str().c_str(), 20, 0, 20); + TH1S *S4Cal1Histo = new TH1S("S4Cal1Histo", sst.str().c_str(), 20, 0, 100); //---------- sst.str(""); sst << "Calorimeter Counters based on " << base.Data(); @@ -158,21 +107,21 @@ //---------- sst.str(""); sst << "Histogram of S3 pattern trigger " << base.Data(); - TH1S *PatternTrig2Histo = new TH1S("PatternTrig2Histo", sst.str().c_str(), 16, 0, 16); + TH1S *PatternTrig2Histo = new TH1S("PatternTrig2Histo", sst.str().c_str(), 12, 0, 12); TCanvas *PatternTrig3 = new TCanvas("PatternTrig3","PatternTrig3", 1280, 1024); //---------- sst.str(""); sst << "Histogram of S2 pattern trigger " << base.Data(); - TH1S *PatternTrig3Histo = new TH1S("PatternTrig3Histo", sst.str().c_str(), 16, 0, 16); + TH1S *PatternTrig3Histo = new TH1S("PatternTrig3Histo", sst.str().c_str(), 8, 0, 8); TCanvas *PatternTrig4 = new TCanvas("PatternTrig4","PatternTrig4", 1280, 1024); //---------- sst.str(""); sst << "Histogram of S12 pattern trigger " << base.Data(); - TH1S *PatternTrig4Histo = new TH1S("PatternTrig4Histo", sst.str().c_str(), 16, 0, 16); + TH1S *PatternTrig4Histo = new TH1S("PatternTrig4Histo", sst.str().c_str(), 12, 0, 12); TCanvas *PatternTrig5 = new TCanvas("PatternTrig5","PatternTrig5", 1280, 1024); @@ -180,9 +129,7 @@ sst.str(""); sst << "Histogram of S11 pattern trigger " << base.Data(); TH1S *PatternTrig5Histo = new TH1S("PatternTrig5Histo", sst.str().c_str(), 16, 0, 16); - - - + int i = 0; Int_t size = nevents; @@ -190,7 +137,18 @@ Double_t *evCounter = new Double_t[size]; Double_t *deadTime = new Double_t[size]; Double_t *liveTime = new Double_t[size]; - + // + Double_t *trigrate0 = new Double_t[size]; + Double_t *trigrate1 = new Double_t[size]; + Double_t *trigrate2 = new Double_t[size]; + Double_t *trigrate3 = new Double_t[size]; + Double_t *trigrate4 = new Double_t[size]; + Double_t *trigrate5 = new Double_t[size]; + // + Double_t *pmt1 = new Double_t[size]; + Double_t *pmt2 = new Double_t[size]; + Double_t *pmt3 = new Double_t[size]; + // while(i < nevents){ HeadBr->GetEntry(i); TrigBr->GetEntry(i); @@ -198,31 +156,37 @@ /* time vector fill*/ ph = eh->GetPscuHeader(); time[i] = ph->GetOrbitalTime(); - + /* eventCounter vector fill*/ evCounter[i] = triggerEvent->evcount; - - /* live time vector fill*/ - liveTime[i] = 0.16 * triggerEvent->dltime[0]; - - /* dead time vector fill*/ - deadTime[i] = 0.01 * triggerEvent->dltime[1]; - + + if( evCounter[i] > 1 ){ + + /* live time vector fill*/ + liveTime[i] = 0.16 * triggerEvent->dltime[0]; + + /* dead time vector fill*/ + deadTime[i] = 0.01 * triggerEvent->dltime[1]; + + } + gStyle->SetStatY(0.85); - /* PMT Plane Counters fill*/ + /////////////////////// + + /* PMT Plane Counters */ - pmt1Histo->Fill((1/0.06)*triggerEvent->pmtpl[0]); - pmt2Histo->Fill((1/0.06)*triggerEvent->pmtpl[1]); - pmt3Histo->Fill((1/0.06)*triggerEvent->pmtpl[2]); - - /* Trigger Rate Counters fill*/ - trRate1Histo->Fill((1/4.)*triggerEvent->trigrate[0]); - trRate2Histo->Fill((1/4.)*triggerEvent->trigrate[1]); - trRate3Histo->Fill((1/4.)*triggerEvent->trigrate[2]); - trRate4Histo->Fill((1/4.)*triggerEvent->trigrate[3]); - trRate5Histo->Fill((1/4.)*triggerEvent->trigrate[4]); - trRate6Histo->Fill((1/4.)*triggerEvent->trigrate[5]); + pmt1[i] = 1.*triggerEvent->pmtpl[0]; + pmt2[i] = 1.*triggerEvent->pmtpl[1]; + pmt3[i] = 1.*triggerEvent->pmtpl[2]; + + /* Trigger Rate Counters*/ + trigrate0[i] = (1/4.)*triggerEvent->trigrate[0]; + trigrate1[i] = (1/4.)*triggerEvent->trigrate[1]; + trigrate2[i] = (1/4.)*triggerEvent->trigrate[2]; + trigrate3[i] = (1/4.)*triggerEvent->trigrate[3]; + trigrate4[i] = (1/4.)*triggerEvent->trigrate[4]; + trigrate5[i] = (1/4.)*triggerEvent->trigrate[5]; /* S4 and Calorimeter Counters fill*/ S4Cal1Histo->Fill(triggerEvent->s4calcount[0]); @@ -238,9 +202,6 @@ int trte1=(int)triggerEvent->patterntrig[1]; PatternTrig1Histo->AddBinContent(1,bit(trte1,0)); - PatternTrig1Histo->AddBinContent(2,bit(trte1,1)); - PatternTrig1Histo->AddBinContent(3,bit(trte1,2)); - int trte2=(int)triggerEvent->patterntrig[2]; @@ -339,6 +300,8 @@ sst << "Event Counter"; TGraph *evCounterGraph = new TGraph(nevents, time, evCounter); evCounterGraph->SetTitle(sst.str().c_str()); + evCounterGraph->GetXaxis()->SetTitle("OrbitalTime"); + evCounterGraph->GetXaxis()->CenterTitle(); evCounterGraph->Draw("APL"); /* @@ -355,12 +318,12 @@ sst << "Dead Time"; DeadTimeGraph->SetTitle(sst.str().c_str()); DeadTimeGraph->SetFillColor(2); - DeadTimeGraph->GetXaxis()->SetTitle("OnBoardTime"); + DeadTimeGraph->GetXaxis()->SetTitle("OrbitalTime"); DeadTimeGraph->GetXaxis()->CenterTitle(); DeadTimeGraph->GetYaxis()->SetTitle("DeadTime(ms)"); DeadTimeGraph->GetYaxis()->CenterTitle(); DeadTimeGraph->SetMaximum(100); - DeadTimeGraph->Draw("APL"); + DeadTimeGraph->Draw("AP"); //---------- sst.str(""); sst << "Graph of Live Time vs OBT based on " << base.Data(); @@ -370,44 +333,164 @@ sst << "Live Time"; LiveTimeGraph->SetTitle(sst.str().c_str()); LiveTimeGraph->SetFillColor(2); - LiveTimeGraph->GetXaxis()->SetTitle("OnBoardTime"); + LiveTimeGraph->GetXaxis()->SetTitle("OrbitalTime"); LiveTimeGraph->GetXaxis()->CenterTitle(); LiveTimeGraph->GetYaxis()->SetTitle("LiveTime(ms)"); LiveTimeGraph->GetYaxis()->CenterTitle(); - LiveTimeGraph->Draw("APL"); - - /* PMT Plane Counters draw*/ - + LiveTimeGraph->Draw("AP"); + + /* + * PMT Plane Counters Graph + */ + + TCanvas *PMTCanvas = new TCanvas("PMTCanvas","PMTCanvas", 1280, 1024); + PMTCanvas->Divide(1,3); + //---------- + sst.str(""); + sst << "Graph of PMT Plane Counters (S11 AND S12) vs OBT based on " << base.Data(); + TGraph *PMT1Graph = new TGraph(nevents, time, pmt1); PMTCanvas->cd(1); - pmt1Histo->GetXaxis()->SetTitle("Hz"); - pmt1Histo->Draw(); + sst.str(""); + sst << "PMT Plane Counters (S11 AND S12)"; + PMT1Graph->SetTitle(sst.str().c_str()); + PMT1Graph->SetFillColor(2); + PMT1Graph->GetXaxis()->SetTitle("OrbitalTime"); + PMT1Graph->GetXaxis()->CenterTitle(); + PMT1Graph->GetYaxis()->SetTitle("PMT Plane Counters"); + PMT1Graph->GetYaxis()->CenterTitle(); + PMT1Graph->Draw("AP"); + //---------- + sst.str(""); + sst << "Graph PMT Plane Counters (S21 AND S22) vs OBT based on " << base.Data(); + TGraph *PMT2Graph = new TGraph(nevents, time, pmt2); PMTCanvas->cd(2); - pmt2Histo->GetXaxis()->SetTitle("Hz"); - pmt2Histo->Draw(); + sst.str(""); + sst << "PMT Plane Counters (S21 AND S22)"; + PMT2Graph->SetTitle(sst.str().c_str()); + PMT2Graph->SetFillColor(2); + PMT2Graph->GetXaxis()->SetTitle("OrbitalTime"); + PMT2Graph->GetXaxis()->CenterTitle(); + PMT2Graph->GetYaxis()->SetTitle("PMT Plane Counters"); + PMT2Graph->GetYaxis()->CenterTitle(); + PMT2Graph->Draw("AP"); + //---------- + sst.str(""); + sst << "Graph PMT Plane Counters (S31 AND S32) vs OBT based on " << base.Data(); + TGraph *PMT3Graph = new TGraph(nevents, time, pmt3); PMTCanvas->cd(3); - pmt3Histo->GetXaxis()->SetTitle("Hz"); - pmt3Histo->Draw(); - - /* Trigger Rate Counters draw*/ + sst.str(""); + sst << "PMT Plane Counters (S31 AND S32)"; + PMT3Graph->SetTitle(sst.str().c_str()); + PMT3Graph->SetFillColor(2); + PMT3Graph->GetXaxis()->SetTitle("OrbitalTime"); + PMT3Graph->GetXaxis()->CenterTitle(); + PMT3Graph->GetYaxis()->SetTitle("PMT Plane Rate Counters"); + PMT3Graph->GetYaxis()->CenterTitle(); + PMT3Graph->Draw("AP"); + // + //---------- + + /* + * Trigger Rate Counters draw + */ + + TCanvas *TrigRate1Canvas = new TCanvas("TrigRate1Canvas","TrigRate1Canvas", 1280, 1024); + TrigRate1Canvas->Divide(1,3); + //---------- + sst.str(""); + sst << "Graph of Trigger Rate Counters (S11 OR S12) AND (S21 OR S22) AND (S31 OR S32) vs OBT based on " << base.Data(); + + TGraph *TrigRate0Graph = new TGraph(nevents, time, trigrate0); TrigRate1Canvas->cd(1); - trRate1Histo->GetXaxis()->SetTitle("Hz"); - trRate1Histo->Draw(); + sst.str(""); + sst << "Trigger Rate Counters (S11 OR S12) AND (S21 OR S22) AND (S31 OR S32)"; + TrigRate0Graph->SetTitle(sst.str().c_str()); + TrigRate0Graph->SetFillColor(2); + TrigRate0Graph->GetXaxis()->SetTitle("OrbitalTime"); + TrigRate0Graph->GetXaxis()->CenterTitle(); + TrigRate0Graph->GetYaxis()->SetTitle("Trigger Rate Counters"); + TrigRate0Graph->GetYaxis()->CenterTitle(); + TrigRate0Graph->Draw("AP"); + //---------- + sst.str(""); + sst << "Graph Trigger Rate Counters (S11 AND S12) AND (S21 AND S22) AND (S31 AND S32) vs OBT based on " << base.Data(); + TGraph *TrigRate1Graph = new TGraph(nevents, time, trigrate1); TrigRate1Canvas->cd(2); - trRate2Histo->GetXaxis()->SetTitle("Hz"); - trRate2Histo->Draw(); + sst.str(""); + sst << "Trigger Rate Counters (S11 AND S12) AND (S21 AND S22) AND (S31 AND S32)"; + TrigRate1Graph->SetTitle(sst.str().c_str()); + TrigRate1Graph->SetFillColor(2); + TrigRate1Graph->GetXaxis()->SetTitle("OrbitalTime"); + TrigRate1Graph->GetXaxis()->CenterTitle(); + TrigRate1Graph->GetYaxis()->SetTitle("Trigger Rate Counters"); + TrigRate1Graph->GetYaxis()->CenterTitle(); + TrigRate1Graph->Draw("AP"); + //---------- + sst.str(""); + sst << "Graph Trigger Rate Counters (S21 OR S22) AND (S31 OR S32) vs OBT based on " << base.Data(); + TGraph *TrigRate2Graph = new TGraph(nevents, time, trigrate2); TrigRate1Canvas->cd(3); - trRate3Histo->GetXaxis()->SetTitle("Hz"); - trRate3Histo->Draw(); + sst.str(""); + sst << "Trigger Rate Counters (S21 OR S22) AND (S31 OR S32) "; + TrigRate2Graph->SetTitle(sst.str().c_str()); + TrigRate2Graph->SetFillColor(2); + TrigRate2Graph->GetXaxis()->SetTitle("OrbitalTime"); + TrigRate2Graph->GetXaxis()->CenterTitle(); + TrigRate2Graph->GetYaxis()->SetTitle("Trigger Rate Counters"); + TrigRate2Graph->GetYaxis()->CenterTitle(); + TrigRate2Graph->Draw("AP"); + // + //---------- + // + TCanvas *TrigRate2Canvas = new TCanvas("TrigRate2Canvas","TrigRate2Canvas", 1280, 1024); + TrigRate2Canvas->Divide(1,3); + //---------- + sst.str(""); + sst << "Graph of Trigger Rate Counters (S21 AND S22) AND (S31 AND S32) vs OBT based on " << base.Data(); + + TGraph *TrigRate3Graph = new TGraph(nevents, time, trigrate3); TrigRate2Canvas->cd(1); - trRate4Histo->GetXaxis()->SetTitle("Hz"); - trRate4Histo->Draw(); + sst.str(""); + sst << "Trigger Rate Counters (S21 AND S22) AND (S31 AND S32)"; + TrigRate3Graph->SetTitle(sst.str().c_str()); + TrigRate3Graph->SetFillColor(2); + TrigRate3Graph->GetXaxis()->SetTitle("OrbitalTime"); + TrigRate3Graph->GetXaxis()->CenterTitle(); + TrigRate3Graph->GetYaxis()->SetTitle("Trigger Rate Counters"); + TrigRate3Graph->GetYaxis()->CenterTitle(); + TrigRate3Graph->Draw("AP"); + //---------- + sst.str(""); + sst << "Graph Trigger Rate Counters S12 AND (S21 AND S22) vs OBT based on " << base.Data(); + TGraph *TrigRate4Graph = new TGraph(nevents, time, trigrate4); TrigRate2Canvas->cd(2); - trRate5Histo->GetXaxis()->SetTitle("Hz"); - trRate5Histo->Draw(); + sst.str(""); + sst << "Trigger Rate Counters S12 AND (S21 AND S22)"; + TrigRate4Graph->SetTitle(sst.str().c_str()); + TrigRate4Graph->SetFillColor(2); + TrigRate4Graph->GetXaxis()->SetTitle("OrbitalTime"); + TrigRate4Graph->GetXaxis()->CenterTitle(); + TrigRate4Graph->GetYaxis()->SetTitle("Trigger Rate Counters"); + TrigRate4Graph->GetYaxis()->CenterTitle(); + TrigRate4Graph->Draw("AP"); + //---------- + sst.str(""); + sst << "Graph Trigger Rate Counters (S11 OR S12) AND (S31 OR S32) vs OBT based on " << base.Data(); + TGraph *TrigRate5Graph = new TGraph(nevents, time, trigrate5); TrigRate2Canvas->cd(3); - trRate6Histo->GetXaxis()->SetTitle("Hz"); - trRate6Histo->Draw(); - + sst.str(""); + sst << "Trigger Rate Counters (S11 OR S12) AND (S31 OR S32)"; + TrigRate5Graph->SetTitle(sst.str().c_str()); + TrigRate5Graph->SetFillColor(2); + TrigRate5Graph->GetXaxis()->SetTitle("OrbitalTime"); + TrigRate5Graph->GetXaxis()->CenterTitle(); + TrigRate5Graph->GetYaxis()->SetTitle("Trigger Rate Counters"); + TrigRate5Graph->GetYaxis()->CenterTitle(); + TrigRate5Graph->Draw("AP"); + //---------- + + gStyle->SetOptStat("mr"); + /* S4 and Calorimeter Counters draw*/ S4CalCanvas->cd(1); S4Cal1Histo->GetXaxis()->SetTitle("Hz"); @@ -418,477 +501,299 @@ /* Pattern Trigger draw*/ - gStyle->SetOptStat(1101); + // gStyle->SetOptStat(1101); + gStyle->SetOptStat("mr"); + PatternTrig5->cd(); - PatternTrig5Histo->Draw(""); + PatternTrig5Histo->SetMinimum(0.); + PatternTrig5Histo->Draw(""); PatternTrig4->cd(); + PatternTrig4Histo->SetMinimum(0.); PatternTrig4Histo->Draw(""); PatternTrig3->cd(); + PatternTrig3Histo->SetMinimum(0.); PatternTrig3Histo->Draw(""); PatternTrig2->cd(); + PatternTrig2Histo->SetMinimum(0.); PatternTrig2Histo->Draw(""); PatternTrig0->cd(); + PatternTrig0Histo->SetMinimum(0.); PatternTrig0Histo->Draw(""); PatternTrig1->cd(); + PatternTrig1Histo->SetMinimum(0.); PatternTrig1Histo->Draw(""); - gStyle->SetOptStat(1111); + + gStyle->SetOptStat("mr"); gStyle->SetStatH(0.3); gStyle->SetStatW(0.3); - + + TH1 *htemp; + TCanvas *RatePMT1 = new TCanvas("RatePMT1", "Rate PMT S11A", 1); RatePMT1->Divide (3,3); //---------- RatePMT1->cd(1); PhysicsTr->Draw("pmtcount1[0]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S111A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S111A"); //---------- RatePMT1->cd (2); PhysicsTr->Draw("pmtcount1[1]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S112A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S112A"); //---------- RatePMT1->cd (3); PhysicsTr->Draw("pmtcount1[2]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S113A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S113A"); //---------- RatePMT1->cd(4); PhysicsTr->Draw("pmtcount1[3]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S114A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S114A"); //---------- RatePMT1->cd(5); PhysicsTr->Draw("pmtcount1[4]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S115A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S115A"); //---------- RatePMT1->cd(6); PhysicsTr->Draw("pmtcount1[5]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S116A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S116A"); //---------- RatePMT1->cd(7); PhysicsTr->Draw("pmtcount1[6]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S117A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S117A"); //---------- RatePMT1->cd(8); PhysicsTr->Draw("pmtcount1[7]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S118A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S118A"); TCanvas *RatePMT2 = new TCanvas("RatePMT2", "Rate PMT S11B", 1); RatePMT2->Divide (3,3); //---------- RatePMT2->cd(1); PhysicsTr-> Draw("pmtcount1[8]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S111B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S111B"); //---------- RatePMT2->cd(2); PhysicsTr->Draw("pmtcount1[9]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S112B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S112B"); //---------- RatePMT2->cd(3); PhysicsTr->Draw("pmtcount1[10]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S113B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S113B"); //---------- RatePMT2->cd(4); PhysicsTr->Draw("pmtcount1[11]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S114B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S114B"); //---------- RatePMT2->cd(5); PhysicsTr-> Draw("pmtcount1[12]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S115B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S115B"); //---------- RatePMT2->cd(6); PhysicsTr-> Draw("pmtcount1[13]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S116B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S116B"); //---------- RatePMT2->cd(7); PhysicsTr-> Draw("pmtcount1[14]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S117B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S117B"); //---------- RatePMT2->cd(8); PhysicsTr->Draw("pmtcount1[15]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S118B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S118B"); TCanvas *RatePMT3 = new TCanvas("RatePMT3", "Rate PMT S12", 1); RatePMT3->Divide (3,4); //---------- RatePMT3->cd(1); PhysicsTr-> Draw("pmtcount1[16]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S121A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S121A"); //---------- RatePMT3->cd(2); PhysicsTr-> Draw("pmtcount1[17]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S122A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S122A"); //---------- RatePMT3->cd (3); PhysicsTr-> Draw("pmtcount1[18]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S123A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S123A"); //---------- RatePMT3->cd(4); PhysicsTr->Draw("pmtcount1[19]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S124A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S124A"); //---------- RatePMT3->cd(5); PhysicsTr-> Draw("pmtcount1[20]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S125A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S125A"); //---------- RatePMT3->cd(6); PhysicsTr-> Draw("pmtcount1[21]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S126A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S126A"); //---------- RatePMT3->cd(7); PhysicsTr->Draw("pmtcount1[22]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S121B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S121B"); //---------- RatePMT3->cd(8); PhysicsTr-> Draw("pmtcount1[23]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S122B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S122B"); //---------- RatePMT3->cd (9); PhysicsTr->Draw("pmtcount2[0]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S123B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S123B"); //---------- RatePMT3->cd (10); PhysicsTr->Draw("pmtcount2[1]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S124B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S124B"); //---------- RatePMT3->cd (11); PhysicsTr->Draw("pmtcount2[2]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S125B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S125B"); //---------- RatePMT3->cd (12); PhysicsTr->Draw("pmtcount2[3]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S126B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S126B"); TCanvas *RatePMT4 = new TCanvas("RatePMT4", "Rate PMT S2", 1); RatePMT4->Divide (3,3); //---------- RatePMT4->cd (1); PhysicsTr-> Draw("pmtcount2[4]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S211A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S211A"); //---------- RatePMT4->cd (2); PhysicsTr-> Draw("pmtcount2[5]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S212A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S212A"); //---------- RatePMT4->cd (3); PhysicsTr-> Draw("pmtcount2[6]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S211B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S211B"); //---------- RatePMT4->cd (4); PhysicsTr-> Draw("pmtcount2[7]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S212B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S212B"); //---------- RatePMT4->cd (5); PhysicsTr-> Draw("pmtcount2[8]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S221A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S221A"); //---------- RatePMT4->cd (6); PhysicsTr-> Draw("pmtcount2[9]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S222A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S222A"); //---------- RatePMT4->cd (7); PhysicsTr-> Draw("pmtcount2[10]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S221B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S221B"); //---------- RatePMT4->cd (8); PhysicsTr->Draw("pmtcount2[11]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S222B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S222B"); TCanvas *RatePMT5 = new TCanvas("RatePMT5", "Rate PMT S3", 1); RatePMT5->Divide (3,4); //---------- RatePMT5->cd (1); PhysicsTr-> Draw("pmtcount2[12]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S311A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S311A"); //---------- RatePMT5->cd (2); PhysicsTr-> Draw("pmtcount2[13]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S312A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S312A"); //---------- RatePMT5->cd (3); PhysicsTr-> Draw("pmtcount2[14]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S313A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S313A"); //---------- RatePMT5->cd (4); PhysicsTr-> Draw("pmtcount2[15]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S311B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S311B"); //---------- RatePMT5->cd (5); PhysicsTr-> Draw("pmtcount2[16]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S312B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S312B"); //---------- RatePMT5->cd (6); PhysicsTr-> Draw("pmtcount2[17]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S313B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S313B"); //---------- RatePMT5->cd (7); PhysicsTr-> Draw("pmtcount2[18]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S321A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S321A"); //---------- RatePMT5->cd (8); PhysicsTr-> Draw("pmtcount2[19]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S322A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S322A"); //---------- RatePMT5->cd (9); PhysicsTr-> Draw("pmtcount2[20]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S323A"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S323A"); //---------- RatePMT5->cd (10); PhysicsTr-> Draw("pmtcount2[21]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S321B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S321B"); //---------- RatePMT5->cd (11); PhysicsTr-> Draw("pmtcount2[22]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S322B"); + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S322B"); //---------- + gStyle->SetOptStat("mr"); + gStyle->SetStatH(0.3); + gStyle->SetStatW(0.3); RatePMT5->cd (12); PhysicsTr-> Draw("pmtcount2[23]"); - t=new TLatex(); - t->SetTextFont(32); - t->SetTextColor(1); - t->SetTextAlign(12); - t->SetTextSize(0.12); - t->DrawLatex(5.,20.,"S323B"); - + htemp = (TH1*)gPad->GetPrimitive("htemp"); + htemp->SetTitle("S323B"); + // print the output const char *saveas = format; if ( !strcmp(saveas,"ps") ) { @@ -903,9 +808,12 @@ PMTCanvas->Print(sst.str().c_str()); TrigRate1Canvas->Print(sst.str().c_str()); TrigRate2Canvas->Print(sst.str().c_str()); + + gStyle->SetOptStat("mr"); S4CalCanvas->Print(sst.str().c_str()); - gStyle->SetOptStat(1101); + // gStyle->SetOptStat(1101); + gStyle->SetOptStat("mr"); PatternTrig5->Print(sst.str().c_str()); PatternTrig4->Print(sst.str().c_str()); @@ -913,10 +821,12 @@ PatternTrig2->Print(sst.str().c_str()); PatternTrig0->Print(sst.str().c_str()); PatternTrig1->Print(sst.str().c_str()); + gStyle->SetOptStat("mr"); RatePMT1->Print(sst.str().c_str()); RatePMT2->Print(sst.str().c_str()); RatePMT3->Print(sst.str().c_str()); RatePMT4->Print(sst.str().c_str()); + // sst.str(""); sst << outDirectory.Data() << filename.Data() << "TriggerScan." << saveas << ")"; RatePMT5->Print(sst.str().c_str()); @@ -949,7 +859,9 @@ sst << outDirectory.Data() << filename.Data() << "TriggerScan6." << saveas; S4CalCanvas->SaveAs(sst.str().c_str()); - gStyle->SetOptStat(1101); + // gStyle->SetOptStat(1101); + gStyle->SetOptStat("mr"); + sst.str(""); sst << outDirectory.Data() << filename.Data() << "TriggerScan7." << saveas; @@ -975,6 +887,8 @@ sst << outDirectory.Data() << filename.Data() << "TriggerScan12." << saveas; PatternTrig1->SaveAs(sst.str().c_str()); + gStyle->SetOptStat("mr"); + sst.str(""); sst << outDirectory.Data() << filename.Data() << "TriggerScan13." << saveas; RatePMT1->SaveAs(sst.str().c_str()); @@ -990,6 +904,9 @@ sst.str(""); sst << outDirectory.Data() << filename.Data() << "TriggerScan16." << saveas; RatePMT4->SaveAs(sst.str().c_str()); + + // gStyle->SetOptStat(1111); + gStyle->SetOptStat("mr"); sst.str(""); sst << outDirectory.Data() << filename.Data() << "TriggerScan17." << saveas;