6 line.DrawLine(x,ymin,x,ymax);
11 arr.DrawArrow(x,ymax,x+dx,ymax,0.01,
"|>");
18 leg->SetBorderSize(0);
19 leg->SetTextSize(0.04);
21 leg->SetY1(leg->GetY2()-leg->GetNRows()*0.05);
30 TLatex *
l =
new TLatex(x,y,text);
33 l->SetTextColor(kBlack);
43 TString lasttime =
TString::Format(
"Last updated on: %s", tnow.AsString());
44 TLatex *ltx =
MiscText(0.6,0.96,0.03,lasttime);
54 strcpy(str,instr.c_str());
60 for(
int i=0;
i<strlen(str);
i++){
64 strStripped[
c] = str[
i];
71 strStripped[
c] =
'\0';
78 int GMToffset = (t.Convert(
true) - t.Convert());
79 gSystem->Setenv(
"TZ",
"UTC");
82 int maxt = t.Convert();
83 int mint = TDatime(20140801,0).Convert();
89 if(ptype==0) dt *= 30;
91 if(ptype==3) dt = maxt - mint;
93 TString tsel =
TString::Format(
"&& lastsec > %d && lastsec < %d",maxt - dt + 2*GMToffset, maxt + 2*GMToffset);
97 TString tlab =
"Date (dd/mm/yy - Central Time)";
98 TString tform =
"%d/%m/%y";
99 TString suf =
"Month";
105 tlab =
"Date (dd/mm - Central Time)";
115 tlab =
"Central Time (hour:min)";
123 tbin = (maxt - mint)/30;
125 tlab =
"Date (dd/mm/yy - Central Time)";
131 TString metricstr =
"TMath::Max(TMath::Min(timingPeakStart * 1e-3,230.24),210.26)";
133 TString
title =
"Timing Peak Start (#mus)";
135 double xmin = 210.25;
136 double xmax = 230.25;
138 gStyle->SetTimeOffset(0);
139 gStyle->SetOptStat(0);
141 if(sel==
"") sel =
"1";
143 if(title==
"") title = metricstr;
145 TString selpre =
"((" + sel +
") && preliminary)";
147 TChain *
mytree =
new TChain(
"mytree");
148 mytree->Add(filenames);
152 TCanvas *
c1 =
new TCanvas(
"c1",
"c1",1000,600);
153 c1->SetBottomMargin(0.12);
154 c1->SetLeftMargin(0.12);
155 c1->SetRightMargin(0.07);
156 c1->SetTopMargin(0.07);
158 TLegend *
leg =
new TLegend(0.15,0.65,0.65,0.9);
159 TLegend *legpre =
new TLegend(0.67,0.65,0.95,0.9);
160 mytree->Draw(metricstr,sel,
"hist");
162 hAll = (TH1F*)(mytree->GetHistogram())->
Clone(
"");
163 if(xmin == -9999 || xmax == -9999){
165 xmin = hAll->GetXaxis()->GetXmin();
166 xmax = hAll->GetXaxis()->GetXmax();
176 if(
int(xmax-xmin)%nbins!=0){
178 xmax += (nbins -
int(xmax-xmin)%
nbins)/2;
179 xmin -= (nbins -
int(xmax-xmin)%
nbins);
183 nbins = (xmax -
xmin) / nbins;
188 hAll =
new TH1F(
"hAll",
"hAll",nbins,xmin,xmax);
189 mytree->Draw(metricstr+
">>hAll",sel,
"hist");
192 hAll->SetLineColor(1);
193 hAll->SetLineWidth(2);
194 hAll->SetTitle(
";"+title+
";Subruns");
195 hAll->GetXaxis()->SetTitleOffset(1);
196 hAll->GetYaxis()->SetTitleOffset(1);
197 hAll->GetXaxis()->CenterTitle();
198 hAll->GetYaxis()->CenterTitle();
199 hAll->GetXaxis()->SetNdivisions();
200 hAll->GetYaxis()->SetNdivisions();
201 hAll->SetMinimum(0.1);
202 hAll->SetStats(kFALSE);
203 double max = hAll->GetMaximum();
204 hAll->SetMaximum(
pow(10,1.5*(
log10(max)+1)-1));
205 leg->AddEntry(hAll,
"Good Subruns",
"f");
207 mytree->Draw(metricstr,
"!(pass_all && !preliminary) && "+sel,
"hist same");
208 mytree->GetHistogram()->SetLineColor(1);
209 mytree->GetHistogram()->SetLineWidth(2);
210 mytree->GetHistogram()->SetFillColor(
kOrange+7);
211 leg->AddEntry(mytree->GetHistogram(),
"Failed Slice",
"f");
213 mytree->Draw(metricstr,
"!(pass_db && pass_hits && passemptyspill && pass_timingpeak && pass_time && pass_duration && !preliminary) && "+sel,
"hist same");
214 mytree->GetHistogram()->SetLineColor(1);
215 mytree->GetHistogram()->SetLineWidth(2);
216 mytree->GetHistogram()->SetFillColor(
kGreen+1);
217 leg->AddEntry(mytree->GetHistogram(),
"Failed Diblock",
"f");
219 mytree->Draw(metricstr,
"!(pass_hits && passemptyspill && pass_timingpeak && pass_time && pass_duration && !preliminary) && "+sel,
"hist same");
220 mytree->GetHistogram()->SetLineColor(1);
221 mytree->GetHistogram()->SetLineWidth(2);
222 mytree->GetHistogram()->SetFillColor(
kBlue+2);
223 leg->AddEntry(mytree->GetHistogram(),
"Failed Hit Rate",
"f");
225 mytree->Draw(metricstr,
"!(passemptyspill && pass_timingpeak && pass_time && pass_duration && !preliminary) && "+sel,
"hist same");
226 mytree->GetHistogram()->SetLineColor(1);
227 mytree->GetHistogram()->SetLineWidth(2);
228 mytree->GetHistogram()->SetFillColor(
kBlue-9);
229 leg->AddEntry(mytree->GetHistogram(),
"Failed Empty Spill",
"f");
231 mytree->Draw(metricstr,
"!(pass_timingpeak && pass_time && pass_duration && !preliminary) && "+sel,
"hist same");
232 mytree->GetHistogram()->SetLineColor(1);
233 mytree->GetHistogram()->SetLineWidth(2);
235 mytree->GetHistogram()->SetFillColor(
kRed+1);
236 leg->AddEntry(mytree->GetHistogram(),
"Failed Timing Peak",
"f");
238 mytree->Draw(metricstr,
"!(pass_time && pass_duration && !preliminary) && "+sel,
"hist same");
239 mytree->GetHistogram()->SetLineColor(1);
240 mytree->GetHistogram()->SetLineWidth(2);
241 mytree->GetHistogram()->SetFillColor(kMagenta+2);
242 leg->AddEntry(mytree->GetHistogram(),
"Failed NuMI Live Time",
"f");
244 mytree->Draw(metricstr,
"!(pass_time && !preliminary) && "+sel,
"hist same");
245 mytree->GetHistogram()->SetLineColor(1);
246 mytree->GetHistogram()->SetLineWidth(2);
247 mytree->GetHistogram()->SetFillColor(kBlack);
248 leg->AddEntry(mytree->GetHistogram(),
"Failed Timestamp",
"f");
250 mytree->Draw(metricstr,selpre,
"hist same");
251 mytree->GetHistogram()->SetLineColor(1);
252 mytree->GetHistogram()->SetLineWidth(2);
253 mytree->GetHistogram()->SetFillColor(18);
254 legpre->AddEntry(mytree->GetHistogram(),
" ",
"f");
256 mytree->Draw(metricstr,
"!pass_all && "+selpre,
"hist same");
257 mytree->GetHistogram()->SetLineColor(1);
258 mytree->GetHistogram()->SetLineWidth(2);
259 mytree->GetHistogram()->SetFillColor(
kGreen-9);
260 legpre->AddEntry(mytree->GetHistogram(),
" ",
"f");
262 mytree->Draw(metricstr,
"!(pass_time && pass_duration) && "+selpre,
"hist same");
263 mytree->GetHistogram()->SetLineColor(1);
264 mytree->GetHistogram()->SetLineWidth(2);
265 mytree->GetHistogram()->SetFillColor(kMagenta-7);
266 legpre->AddEntry(mytree->GetHistogram(),
" ",
"f");
268 mytree->Draw(metricstr,
"!(pass_time) && "+selpre,
"hist same");
269 mytree->GetHistogram()->SetLineColor(1);
270 mytree->GetHistogram()->SetLineWidth(2);
271 mytree->GetHistogram()->SetFillColor(kGray+2);
272 legpre->AddEntry(mytree->GetHistogram(),
" Preliminary",
"f");
275 leg->SetY1(0.9-leg->GetNRows()*0.05);
277 legpre->SetNColumns(2);
279 legpre->SetColumnSeparation(-0.05);
283 MiscText(0.15,0.94,0.055,
"Near Detector");
285 DrawArrow(216.75,-1,0.1,0.7*hAll->GetBinContent(hAll->GetMaximumBin()));
286 DrawArrow(219.25,1,0.1,0.7*hAll->GetBinContent(hAll->GetMaximumBin()));
287 c1->SaveAs(
"pngfiles/NearDet-t00-P1GoodDataSelTimingPeakStart"+suf+
".png");
std::map< std::string, double > xmax
void StripSpecialChar(string &instr)
void UpdateStamp(int ptype, int maxt)
void GoodDataSelTimingPeakStart(int ptype=0)
void SetLeg(TLegend *leg)
TLatex * MiscText(float x, float y, float size, TString text)
fvar< T > floor(const fvar< T > &x)
void DrawArrow(double x, double dx, double ymin, double ymax)
void Format(TGraph *gr, int lcol, int lsty, int lwid, int mcol, int msty, double msiz)
T max(sqlite3 *const db, std::string const &table_name, std::string const &column_name)
fvar< T > ceil(const fvar< T > &x)