tbpots.C File Reference

Go to the source code of this file.

Functions
double	clamp (const double low, const double x, const double high)
double	rounduptenths (const double x)
double	rounddowntenths (const double x)
void	pots (const char *const in, const char *out)
void	tbpots (const char *const more="")

Function Documentation

double clamp	(	const double	low,
		const double	x,
		const double	high
	)

Definition at line 1 of file tbpots.C.

References submit_syst::x.

{
  if(x < low) return low;
  if(x > high) return high;
  return x;
}

void pots	(	const char *const	in,
		const char *	out
	)

Definition at line 18 of file tbpots.C.

References a, b, bins, plot_validation_datamc::c, demo5::c1, MECModelEnuComparisons::i, infile, calib::j, kBlue, kRed, MECModelEnuComparisons::leg, extractScale::mean, submit_nova_art::mode, getGoodRuns4SAM::n, rounduptenths(), sum, confusionMatrixTree::t, tmp, shutoffs::week, and Munits::year.

Referenced by tbpots().

{
  const double textsize = 0.05;
  const unsigned int week = 7, year = 365, month = 30;

  TCanvas *c1 = new TCanvas("c1", "c1",0,0,800,500);
  c1->SetCanvasSize(1000,500);
  gStyle->SetOptStat(0);
  gStyle->SetOptTitle(0);
  c1->SetGrid(0,1);
  gStyle->SetGridStyle(kSolid);
  gStyle->SetGridColor(kGray);
  c1->SetLeftMargin(0.08);
  c1->SetBottomMargin(0.08);
  c1->SetTopMargin(0.025);
  c1->SetRightMargin(0.23);
  c1->SetTicky(1);
  c1->SetTickx(1);
  c1->SetFrameLineWidth(2);

  ifstream infile(in);
  
  vector<string> date;
  vector<double> uptimefrac, potdelivered, potrecorded, bins;

  struct datatmp{
    string s;
    double a, b, c;
  };
  datatmp tmp;

  int n = 0;
  while(infile >> tmp.s >> tmp.a >> tmp.b >> tmp.c){
    date        .push_back(tmp.s); 
    uptimefrac  .push_back(tmp.a);
    potrecorded .push_back(tmp.b);
    potdelivered.push_back(tmp.c);
    bins        .push_back(1 + n++);
  }

  double maxpots = 1;
  for(int i = 0; i < potdelivered.size(); i++)  
    if(potdelivered[i] > maxpots)
      maxpots = potdelivered[i];


  TGraph * delivered = new TGraph(date.size(), &bins[0], &potdelivered[0]);
  delivered->SetName("delivered");
  TGraph * recorded = new TGraph(date.size(), &bins[0], &potrecorded[0]);
  recorded->SetName("recorded");

  enum mode{weekmode, monthmode, yearmode } themode;
  if(bins.size() < month+1) themode = weekmode;
  else if(bins.size() < month*4) themode = monthmode;
  else themode = yearmode;
  const double highy = rounduptenths(maxpots);
  
  const int nbin = themode == weekmode?  (week*((bins.size() +  week - 1)/ week)+1):
                   themode ==monthmode? (month*((bins.size() + month - 1)/month)+1):
                                         (year*((bins.size() +  year - 1)/ year)+1);

  TH2D * dum = new TH2D("dum", "", nbin, 0, nbin, 1, 0, highy);

  int div =  themode == weekmode?  700 + (bins.size() + week - 1)/week:
    themode ==monthmode?  300 + (bins.size() + month- 1)/month:
                         1200 + (bins.size() + year - 1)/year;
  TGaxis * xa = new TGaxis(1, 0, nbin, 0, 0, nbin, div, "N");
  TGaxis * xat= new TGaxis(1, highy, nbin, highy, 0, nbin, div, "N-");
  xa->SetLabelSize(0);
  xat->SetLabelSize(0);

  for(int i = 0; i < bins.size(); i += (themode == weekmode? week :
                                        themode == monthmode?month: year))
    //      shift about half a bin right ---v
    dum->GetXaxis()->SetBinLabel(i+1, Form("   %s", date[i].c_str()));

  dum->GetXaxis()->SetTitle("");
  dum->GetYaxis()->SetTitle("10^{6} scintillator counts per day");
  dum->GetYaxis()->SetTitleSize(textsize);
  dum->GetXaxis()->SetTitleSize(textsize);
  dum->GetYaxis()->SetLabelSize(textsize);
  dum->GetXaxis()->SetLabelSize(textsize*1.5);
  dum->GetXaxis()->SetLabelOffset(0.01);
  dum->GetYaxis()->SetTitleOffset(0.70);
  dum->GetXaxis()->SetBit(TAxis::kLabelsHori);
  dum->GetYaxis()->CenterTitle(1);
  dum->GetYaxis()->SetDecimals();
  dum->GetYaxis()->SetTickSize(0.01);
  dum->GetXaxis()->CenterTitle(1);
  dum->GetXaxis()->SetTickSize(0);
  dum->Draw();
  xa->Draw();
  xat->Draw();

  delivered->SetMarkerStyle(kOpenDiamond);
  delivered->SetMarkerSize(themode == yearmode?0.3: 1.3);
  delivered->SetMarkerColor(kBlue+1);
  delivered->Draw("p");

  recorded->SetMarkerStyle(kMultiply);
  recorded->SetMarkerSize(themode == yearmode?0.3: 1.3);
  recorded->SetMarkerColor(kRed);
  recorded->Draw("p");

  TGraph * deliveredbyweek = dynamic_cast<TGraph*>(delivered->Clone("deliveredbyweek"));

  for(int i = 0; i < delivered->GetN(); i+=week){
    double sum = 0;
    int n = 0;
    for(int j = i; j < i+week && j < delivered->GetN(); j++){
      sum += delivered->GetY()[j];
      n++;
    }
    const double mean = sum/n;
    for(int j = i; j < i+week && j < delivered->GetN(); j++){

      // Make transitions sharper
      if(j == i) deliveredbyweek->GetX()[j] -= 0.25;
      if(j == i+week-1 || j == delivered->GetN()-1) deliveredbyweek->GetX()[j] += 0.25;

      deliveredbyweek->GetY()[j] = mean;
    }
  }

  deliveredbyweek->SetLineColor(kBlue+1);
  deliveredbyweek->SetLineWidth(themode == weekmode?3:2);
  deliveredbyweek->Draw("l"); 

  TGraph * recordedbyweek = dynamic_cast<TGraph*>(recorded->Clone("recordedbyweek"));

  double lastweeksdelivered = 0, lastweeksrecorded = 0;
  for(int i = 0; i < recorded->GetN(); i+=week){
    double sumdem = 0, sumnum = 0, sum = 0;
    double sumx = 0;
    int n = 0;
    for(int j = i; j < i+week && j < recorded->GetN(); j++){
      sumnum += recorded->GetY()[j];
      sumdem += delivered->GetY()[j];
      sumx += recorded->GetX()[j];
      n++;
    }
    const double meaneff = sumnum/sumdem;
    const double mean = sumnum/n;
    const double meanx = sumx/n;
    for(int j = i; j < i+week && j < recorded->GetN(); j++){
      // Make transitions sharper
      if(j == i) recordedbyweek->GetX()[j] -= 0.25;
      if(j == i+week-1 || j == recorded->GetN()-1) recordedbyweek->GetX()[j] += 0.25;

      recordedbyweek->GetY()[j] = mean;
    }

    if(n == 7){ // save last full week
      lastweeksdelivered = sumdem;
      lastweeksrecorded  = sumnum;
    }

    if(recorded->GetN() <= week*9){
      TText * t = meaneff!=meaneff?new TText(meanx, i%2?highy/8.:highy/6., ""):
                  new TText(meanx, i%2?highy/8.:highy/6., Form("%.1f%%", meaneff*100));
      t->SetTextSize(textsize);
      t->SetTextFont(42);
      t->SetTextAlign(22);
      t->SetTextColor(kRed);
      t->Draw();
    }
  }

  recordedbyweek->SetLineColor(kRed);
  recordedbyweek->SetLineWidth(themode == weekmode?3:2);
  recordedbyweek->Draw("l"); 

  c1->SaveAs(Form("%s-noleg.pdf", out));
  c1->SaveAs(Form("%s-noleg.eps", out));
  c1->SaveAs(Form("%s-noleg.png", out));

  const double legbound = 0.45;

  TLegend leg(0.77, legbound, 0.995, 0.81);
  leg.SetTextSize(textsize);
  leg.SetBorderSize(0);
  leg.SetMargin(0.19);
  leg.SetFillStyle(0);
  leg.AddEntry((TH1D*)NULL, "Test Beam", "");
  leg.AddEntry((TH1D*)NULL, "UTC days", "");
  leg.AddEntry(delivered, "Delivered daily", "p");
  leg.AddEntry(recorded, "Recorded daily", "p");
  leg.AddEntry(deliveredbyweek, "Delivered weekly", "l");
  leg.AddEntry(recordedbyweek, "Recorded weekly", "l");
  leg.Draw();

  TLegend leg2(0.77, 0.24, 0.995, legbound-0.03);
  leg2.SetTextSize(textsize);
  leg2.SetBorderSize(0);
  leg2.SetMargin(0.04);
  leg2.SetFillStyle(0);
  leg2.AddEntry((TH1D*)NULL, "Last full week", "");
  leg2.AddEntry((TH1D*)NULL, Form("%.2f#times10^{6} delivered",lastweeksdelivered), "");
  leg2.AddEntry((TH1D*)NULL, Form("%.2f#times10^{6} recorded", lastweeksrecorded), "");
  leg2.Draw();

  c1->SaveAs(Form("%s.pdf", out));
  c1->SaveAs(Form("%s.eps", out));
  c1->SaveAs(Form("%s.png", out));
}

double rounddowntenths

(

const double

x

)

Definition at line 13 of file tbpots.C.

References makeTrainCVSamples::int.

{
  return int(x*10)/10.;
}

double rounduptenths

(

const double

x

)

Definition at line 8 of file tbpots.C.

References makeTrainCVSamples::int.

Referenced by pots().

{
  return int(x*10 + 1)/10.;
}

void tbpots

(

const char *const

more = ""

)

Definition at line 224 of file tbpots.C.

References pots().

{
  pots(Form("tb_for_plots%s.dat", more), Form("tbpots%s", more));
}