cana_time_peak.C File Reference

#include <iostream>
#include <cmath>
#include "TArrow.h"
#include "TCanvas.h"
#include "TChain.h"
#include "TCut.h"
#include "TH1F.h"
#include "TMath.h"
#include "TMinuit.h"
#include "TStyle.h"
#include "TText.h"
#include "TLegend.h"
#include "TLatex.h"

Go to the source code of this file.

Functions
TCut	run ("run>=15071&&run<=15131")
TCut	activity ("nhit>=25&&nhitx>=3&&nhity>=3")
TCut	y1cut ("y1<150")
TCut	x1cut ("abs(x1)<100")
TCut	z1cut ("280<z1 && z1<1200")
TCut	fiducial (x1cut &&y1cut &&z1cut)
TCut	cosnumi ("cosnumi>0.95")
TCut	canacuts (run &&fiducial &&activity &&cosnumi)
TCut	intime ("tave>=217.0 && tave<=227.8")
TCut	outtime ("(tave>0.0&&tave<200.0) || (tave>250.0&&tave<450.0)")
void	NDOSLabel ()
void	Preliminary ()
	Put NOvA Preliminary on plots. More...
void	CenterTitles (TH1 *histo)
void	nhit_proj (TChain &ch)
void	nhit_plot ()
void	y1_proj (TChain &ch)
void	y1_plot ()
void	x1_proj (TChain &ch)
void	x1_plot ()
void	z1_proj (TChain &ch)
void	z1_plot ()
void	cosnumi_proj (TChain &ch)
void	cosnumi_plot ()
void	time_wide_proj (TChain &ch)
void	time_wide_plot ()
void	time_zoom_proj (TChain &ch)
double	spill_profile (double t)
void	make_template (double sigma)
void	make_time_profile (double b, double f, double t0, double sigma)
double	chi2 ()
void	fcn (int &, double *, double &f, double *p, int)
void	minimize ()
void	time_zoom_plot ()
void	proj_all (TChain &ch)
void	plot_all ()

Variables
const char *	label = "NDOS Data, 09/12/13 - 09/27/13"
double	t1 = 217.0
double	t2 = 227.8
double	intimet = (t2-t1)
double	outtimet = 400.0
TH1F *	nhith
TH1F *	y1h
TH1F *	x1h
TH1F *	z1h
TH1F *	cosnumiinh
TH1F *	cosnumiouth
TH1F *	tavewideh
TH1F *	tavezoomh
int	ntbin_per_usec = 20
double	zoom_tbin_sz = 1.0/(float)ntbin_per_usec
double	zoom_tlo = floor(t1) - 5
double	zoom_thi = floor(t2) + 6
int	zoom_nbin = ntbin_per_usec*(zoom_thi-zoom_tlo)
static const int	nbatch = 6
static const double	dt_batch = 81*18.83E-3
static const double	dt_gap = 5*18.83E-3
static const double	dt_cycle = dt_batch + dt_gap
static const double	spill_length = nbatch*dt_batch + (nbatch-1)*dt_gap
TH1F *	template_spillh = 0
double	template_tlo
double	template_thi
int	template_nbin
TH1F *	time_profileh

Function Documentation

TCut activity

(

"nhit>=25&&nhitx>=3&&nhity>=3"

)

Referenced by cosnumi_proj(), x1_proj(), y1_proj(), and z1_proj().

TCut canacuts

(

run &&fiducial &&activity &&

cosnumi

)

Referenced by time_wide_proj(), and time_zoom_proj().

void CenterTitles

(

TH1 *

histo

)

Definition at line 77 of file cana_time_peak.C.

Referenced by cosnumi_plot(), nhit_plot(), time_wide_plot(), time_zoom_plot(), x1_plot(), y1_plot(), and z1_plot().

{
  histo->GetXaxis()->CenterTitle();
  histo->GetYaxis()->CenterTitle();
  histo->GetZaxis()->CenterTitle();  
}

double chi2

(

)

Definition at line 532 of file cana_time_peak.C.

References MECModelEnuComparisons::i, test_ParserArtEvents::log, nui, sum, tavezoomh, time_profileh, and zoom_nbin.

Referenced by upmuana::UpMuRecoAna::analyze(), upmuana::UpMuAnalysis::analyze(), cmf::ChiSqrCalculator::BinCNPChiSqr(), cmf::ChiSqrCalculator::BinPoissonChiSqr(), ana::CalcChi2(), CalcChi2(), Calculate1DChi2_simple(), Calculate2DChi2(), CalculateENuChi2(), CalculateQ2Chi2(), calculateWrongSignNumuQ1(), calculateWrongSignNumuQ2(), calculateWrongSignNumuQ3(), calculateWrongSignNumuQ4(), vdt::Segment::Chi2EndPoint(), ana::OscCovMxExperiment::ChiSq(), ana::CovMxExperiment::ChiSq(), cmf::ChiSqrCalculator::ChiSqrCNP(), cmf::ChiSqrCalculator::ChiSqrPoisson(), DedxComp(), DedxCompBirk(), DedxGamma(), align::SplitTracks::DemingRegFit(), doUnfolding(), FCContour(), fcn(), FindAngleFromChi2(), cosrej::CosRejFxs::fitchisq(), numusand::NumuSandFxs::fitchisq(), bpfit::DimuonFitter::FitView(), util::LinFit(), novaddt::TrackFit::LinFit(), geo::LinFit(), novaddt::TrackFit::LinFitMinDperp(), geo::LinFitMinDperp(), makeMatrixElementSurface(), nuebar_signif(), PionDedxComp(), ProtonDedxComp(), ProtonDedxCompBirk(), ana::test_best_spline(), and novaddt::TrackFit::UtilLinFit().

{
  int i;
  double sum = 0;
  double ni, nui;
  for (i=1; i<=zoom_nbin; ++i) {
    nui = time_profileh->GetBinContent(i);
    ni  = tavezoomh->    GetBinContent(i);
    sum += nui;
    if (ni>0.0) {
      sum += ni*(log(ni/nui)-1);
    }
  }
  return 2.0*sum;
}

TCut cosnumi	(	"	cosnumi,
		0.95"
	)

void cosnumi_plot

(

)

Definition at line 297 of file cana_time_peak.C.

References a1, a2, CenterTitles(), cosnumiinh, cosnumiouth, Draw(), NDOSLabel(), Preliminary(), std::sqrt(), x1, and y1.

Referenced by plot_all().

{
  gStyle->SetOptStat(0);
  
  TCanvas* cosnumic = new TCanvas("cosnumic","cosnumic",800,600);

  cosnumiinh->SetMarkerStyle(20);
  cosnumiouth->SetLineColor(kRed);
  cosnumiouth->SetLineWidth(2);

  cosnumiouth->SetMaximum(cosnumiinh->GetMaximum()+
                          1.5*sqrt(cosnumiinh->GetMaximum()));

  CenterTitles(cosnumiouth);
  cosnumiouth->Draw("hist");
  cosnumiinh-> Draw("same,e");

  Preliminary();

  double x1 = cosnumiouth->GetBinCenter(4);
  double y1 = cosnumiouth->GetBinContent(4);
  
  TLine* a1 = new TLine(x1+0.08, 1.12*y1, x1+0.03, 1.01*y1);
  a1->SetLineWidth(1);
  a1->SetLineColor(kRed);
  a1->Draw();

  TText* tx1 = new TText(x1+0.09, 1.14*y1,"out of time tracks");
  tx1->SetTextColor(kRed);
  tx1->Draw();

  x1 = cosnumiinh->GetBinCenter(18);
  y1 = cosnumiinh->GetBinContent(18);

  TLine* a2 = new TLine(x1-0.04, 1.01*y1, x1-0.08, 1.12*y1);
  a2->SetLineWidth(1);
  a2->SetLineColor(kBlack);
  a2->Draw();

  TText* tx2 = new TText(x1-0.06, 1.20*y1,"in time tracks");
  tx2->SetTextAlign(31);
  tx2->SetTextColor(kBlack);
  tx2->Draw();

  NDOSLabel();

  cosnumic->SaveAs("cana_cosnumic.pdf");
  cosnumic->SaveAs("cana_cosnumic.png");
}

void cosnumi_proj

(

TChain &

ch

)

Definition at line 284 of file cana_time_peak.C.

References activity(), cosnumiinh, cosnumiouth, fiducial(), intime(), intimet, outtime(), and outtimet.

Referenced by proj_all().

{
  cosnumiinh = new TH1F("cosnumiinh",
                        ";cos #theta_{NuMI}; Events / 0.05",
                        20, 0.0, 1.0);
  cosnumiouth = new TH1F("cosnumiouth",
                         ";cos #theta_{NuMI}; Events / 0.05",
                         20, 0.0, 1.0);
  ch.Project("cosnumiinh", "cosnumi",run&&activity&&fiducial&&intime);
  ch.Project("cosnumiouth","cosnumi",run&&activity&&fiducial&&outtime);

  cosnumiouth->Scale(intimet/outtimet);
}

void fcn	(	int &	,
		double *	,
		double &	f,
		double *	p,
		int
	)

Definition at line 550 of file cana_time_peak.C.

References chi2(), and make_time_profile().

Referenced by minimize().

{
  make_time_profile(p[0], p[1], p[2], p[3]);
  f = chi2();
}

TCut fiducial

(

x1cut &&y1cut &&

z1cut

)

Referenced by cosnumi_proj(), tree_maker_simple(), and useBDTG().

TCut intime

(

"tave>=217.0 && tave<=227.8"

)

Referenced by CalcRWithSystsNus17(), cosnumi_proj(), getHists_FNEX(), GetNueData2017(), GetNueData2018(), GetNueData2019(), ana::GetNueData2020(), GetNumuData2019(), ana::GetNumuData2020(), make_nus17_fc_surfs(), make_nus_fc_surfs(), makeMatrixElementSurface(), plot_datamcpred(), plot_predictions(), PlotNus17PredSystsData(), and PlotNusSensAna01().

void make_template

(

double

sigma

)

Definition at line 448 of file cana_time_peak.C.

References stan::math::floor(), Gaus(), MECModelEnuComparisons::i, calib::j, spill_length, spill_profile(), template_nbin, template_spillh, template_thi, template_tlo, w, and zoom_tbin_sz.

Referenced by make_time_profile().

{
  if (template_spillh!=0) {
    delete template_spillh;
    template_spillh = 0;
  }
  if (template_spillh==0) {
    template_tlo =
      zoom_tbin_sz*floor(-5.0*sigma/zoom_tbin_sz);
    template_thi =
      zoom_tbin_sz*
      floor((spill_length+5.0*sigma)/zoom_tbin_sz+1);
    template_nbin = 10*(template_thi-template_tlo)/zoom_tbin_sz;
    template_spillh = new TH1F("tempalte_spillh",";usec",
                               template_nbin,
                               template_tlo,
                               template_thi);
  }
  
  int i, j;
  double ti, tj;
  double s, w;
  for (i=1; i<=template_nbin; ++i) {
    ti = template_spillh->GetBinCenter(i);
    s = spill_profile(ti);
    for (j=1; j<template_nbin; ++j) {
      tj = template_spillh->GetBinCenter(j);
      w = TMath::Gaus(ti,tj,sigma);
      
      template_spillh->Fill(tj, w*s);
    }
  }
  template_spillh->Scale(template_spillh->GetMaximum());
}

void make_time_profile	(	double	b,
		double	f,
		double	t0,
		double	sigma
	)

Definition at line 488 of file cana_time_peak.C.

References MECModelEnuComparisons::i, make_template(), sigma(), confusionMatrixTree::t, template_nbin, template_spillh, time_profileh, zoom_nbin, zoom_thi, and zoom_tlo.

Referenced by fcn(), and minimize().

{
  int i;
  
  static int ifirst = 1;
  if (ifirst) {
    time_profileh = new TH1F("time_profileh",
                             ";usec",
                             zoom_nbin,
                             zoom_tlo,
                             zoom_thi);
    ifirst = 0;
  }
  time_profileh->Reset();

  static double last_sigma = 0.0;
  if (last_sigma != sigma) {
    last_sigma = sigma;
    make_template(sigma);
  }
  
  double s, t;
  //
  // Put the background in
  //
  for (i=1; i<=time_profileh->GetNbinsX(); ++i) {
    t = time_profileh->GetBinCenter(i);
    time_profileh->Fill(t, b);
  }
  //
  // Put the signal in
  //
  for (i=1; i<=template_nbin; ++i) {
    t = template_spillh->GetBinCenter(i);
    s = template_spillh->GetBinContent(i);
    time_profileh->Fill(t+t0, f*s);
  }
}

void minimize

(

)

Definition at line 558 of file cana_time_peak.C.

References PandAna.Demos.pi0_spectra::bkg, om::cerr, e, allTimeWatchdog::endl, fcn(), stan::math::fmin(), MECModelEnuComparisons::i, make_time_profile(), std::sqrt(), getGoodRuns4SAM::t0, tavezoomh, febshutoff_auto::val, submit_syst::x, submit_syst::y, and zoom_nbin.

Referenced by proj_all().

{
  int ierr;
  TMinuit minuit;
  minuit.SetPrintLevel(-99);
  minuit.SetFCN(fcn);
  minuit.mnparm(0, "bkg",  1e-3, 1e-4, 0, 10.0, ierr);
  minuit.mnparm(1, "sf",   1e-3, 1e-4, 0, 10.0, ierr);
  minuit.mnparm(2, "t0",   1e-3, 1e-4, 0, 10.0, ierr);
  minuit.mnparm(3, "sigt", 1e-3, 1e-4, 0, 10.0, ierr);
  
  double chi2best    = 99E99;
  double p[4];
  double sigp[4];
  
  int i;
  double t0;
  double x[zoom_nbin];
  double y[zoom_nbin];
  for (i=0; i<zoom_nbin; ++i) {
    t0 = tavezoomh->GetBinCenter(i+1);

    minuit.mnparm(2, "t0", t0, 1e-3, 0.0, 500.0, ierr);
    minuit.FixParameter(2);

    minuit.mnparm(3, "sigt", 0.512/sqrt(12)/sqrt(25), 1e-3, 0.0, 500.0, ierr);
    minuit.FixParameter(3);
    
    minuit.Migrad();
    
    double fmin, fedm, errdef;
    int npari, nparx, istat;
    minuit.mnstat(fmin,fedm,errdef,npari,nparx,istat);
    x[i] = t0;
    y[i] = fmin;
    std::cerr << i << "\t" << t0 << "\t" << fmin << std::endl;
    if (fmin<chi2best) {
      chi2best = fmin;
      
      double val, err, xlolim, xuplim;
      int iuint;

      static TString bkg("bkg");
      minuit.mnpout(0, bkg, val, err, xlolim, xuplim, iuint);
      p[0]    = val;
      sigp[1] = err;

      static TString sf("sf");
      minuit.mnpout(1, sf, val, err, xlolim, xuplim, iuint);
      p[1]    = val;
      sigp[1] = err;
      
      static TString t0s("t0");
      minuit.mnpout(2, t0s, val, err, xlolim, xuplim, iuint);
      p[2]    = val;
      sigp[2] = err;

      static TString sigt("sigt");
      minuit.mnpout(3, sigt, val, err, xlolim, xuplim, iuint);
      p[3]    = val;
      sigp[3] = err;
    }
  }
  make_time_profile(p[0], p[1], p[2], p[3]);

  std::cerr << "chi2 = " << chi2best << std::endl;
  std::cerr << "Best b   = " << p[0] <<" +/- " << sigp[0] << std::endl;
  std::cerr << "Best sf  = " << p[1] <<" +/- " << sigp[1] << std::endl;
  std::cerr << "Best t0  = " << p[2] <<" +/- " << sigp[2] << std::endl;
  std::cerr << "Best sigt= " << p[3] <<" +/- " << sigp[3] << std::endl;
}

void NDOSLabel

(

)

Definition at line 56 of file cana_time_peak.C.

References label.

Referenced by cosnumi_plot(), nhit_plot(), time_wide_plot(), time_zoom_plot(), x1_plot(), y1_plot(), and z1_plot().

{
  TLatex* tx = new TLatex(0.1, .93, label);
  tx->SetTextColor(kBlack);
  tx->SetNDC();
  tx->SetTextSize(1/30.);
  tx->SetTextAlign(12);
  tx->Draw();
}

void nhit_plot

(

)

Definition at line 95 of file cana_time_peak.C.

References CenterTitles(), std::log10(), NDOSLabel(), nhith, cet::pow(), and Preliminary().

Referenced by plot_all().

{
  gStyle->SetOptStat(0);

  TCanvas* nhitc = new TCanvas("nhitc","nhitc",800,600);
  nhitc->SetLogy();
  
  nhith->SetFillColor(kGreen-10);
  CenterTitles(nhith);
  nhith->Draw();
  Preliminary();

  TLine* nhitl = new TLine(25.0, 0, 25.0, nhith->GetMaximum());
  nhitl->SetLineStyle(3);
  nhitl->Draw();
  double ya1 = 0.75*log10(nhith->GetMaximum());
  double ya2 = 0.75*log10(nhith->GetMaximum())+0.9;
  double ya3 = pow(10,0.5*(ya1+ya2));
  ya1 = pow(10,ya1);
  ya2 = pow(10,ya2);
  TArrow* nhita1 = new TArrow(25.0, ya1, 30.0, ya1, 0.03);
  TArrow* nhita2 = new TArrow(25.0, ya2, 30.0, ya2, 0.03);
  nhita1->Draw();
  nhita2->Draw();

  TText* tx2 = new TText(30.0, ya3, "selected");
  tx2->SetTextFont(42);
  tx2->SetTextAlign(12);
  tx2->SetTextSize(0.035);
  tx2->Draw();

  NDOSLabel();

  nhitc->SaveAs("cana_nhit.pdf");
  nhitc->SaveAs("cana_nhit.png");
}

void nhit_proj

(

TChain &

ch

)

Definition at line 88 of file cana_time_peak.C.

References nhith.

Referenced by proj_all().

{
  nhith = new TH1F("nhith",";N_{hit};Events",101,-0.5,100.5);
  nhith->SetFillColor(kGreen-10);
  ch.Project("nhith","nhit", run);
}

TCut outtime

(

"(tave>0.0&&tave<200.0) || (tave>250.0&&tave<450.0)"

)

Referenced by cosnumi_proj(), getHists_FNEX(), GetNueCosmics2017(), GetNueCosmics2018(), and GetNueCosmicsFuture().

void plot_all

(

)

Definition at line 687 of file cana_time_peak.C.

References cosnumi_plot(), nhit_plot(), time_wide_plot(), time_zoom_plot(), x1_plot(), y1_plot(), and z1_plot().

{
  nhit_plot();
  y1_plot();
  x1_plot();
  z1_plot();
  cosnumi_plot();
  time_wide_plot();
  time_zoom_plot();
}

void Preliminary

(

)

Put NOvA Preliminary on plots.

Definition at line 67 of file cana_time_peak.C.

References prelim.

Referenced by cosnumi_plot(), nhit_plot(), time_wide_plot(), time_zoom_plot(), x1_plot(), y1_plot(), and z1_plot().

{
  TLatex* prelim = new TLatex(.9, .93, "NOvA Preliminary");
  prelim->SetTextColor(kBlue);
  prelim->SetNDC();
  prelim->SetTextSize(1/30.);
  prelim->SetTextAlign(32);
  prelim->Draw();
}

void proj_all

(

TChain &

ch

)

Definition at line 663 of file cana_time_peak.C.

References cosnumi_proj(), om::cout, allTimeWatchdog::endl, minimize(), nhit_proj(), time_wide_proj(), time_zoom_proj(), x1_proj(), y1_proj(), and z1_proj().

{

  std::cout<<"starting proj"<<std::endl;
  nhit_proj(ch);
  std::cout<<"nhit done"<<std::endl;
  y1_proj(ch);
  std::cout<<"y1 done"<<std::endl;
  x1_proj(ch);
  std::cout<<"x1 done"<<std::endl;
  z1_proj(ch);
  std::cout<<"z1 done"<<std::endl;
  cosnumi_proj(ch);
  std::cout<<"cosnumi done"<<std::endl;
  time_wide_proj(ch);
  std::cout<<"time wide done"<<std::endl;
  time_zoom_proj(ch);
  std::cout<<"time zoom done"<<std::endl;
  minimize();
  std::cout<<"minimize done"<<std::endl;
}

TCut run

(

"run>=15071&&run<=15131"

)

double spill_profile

(

double

t

)

Definition at line 424 of file cana_time_peak.C.

References dt_batch, dt_cycle, spill_length, and confusionMatrixTree::t.

Referenced by make_template().

{
  if (t<0.0)          return 0.0;
  if (t>spill_length) return 0.0;
  
  // Rewind cycle by cycle
  double ttmp = t;
  while (ttmp>dt_cycle) ttmp -= dt_cycle;
  
  // Check it we are inside the "on" part of the cycle
  if (ttmp<0.0)      return 0.0;
  if (ttmp>dt_batch) return 0.0;

  return 1.0;
}

void time_wide_plot

(

)

Definition at line 357 of file cana_time_peak.C.

References b, CenterTitles(), om::cout, allTimeWatchdog::endl, check_time_usage::float, stan::math::floor(), NDOSLabel(), Preliminary(), std::sqrt(), t1, t2, and tavewideh.

Referenced by plot_all().

{  
  gStyle->SetOptStat(0);
  
  TCanvas* cana_tave_wide = 
    new TCanvas("cana_tave_wide","cana_tave_wide",800,600);
  tavewideh->SetFillColor(kGreen-10);
  CenterTitles(tavewideh);
  tavewideh->SetMaximum(floor(tavewideh->GetMaximum())+1);
  tavewideh->Draw();
  Preliminary();
  
  unsigned int i1 = tavewideh->FindBin(t1);
  unsigned int i2 = tavewideh->FindBin(t2);
  double splusb = tavewideh->Integral(i1,i2);
  double b = 
    tavewideh->Integral(i1-25,i1-1) + 
    tavewideh->Integral(i2+1,i2+25);
  b *= (float)(i2-i1+1)/50.0;
  double fom = (splusb-b)/sqrt(b);
  std::cout << "FOM: " << fom << std::endl;

  NDOSLabel();

  /*
  TText* tx4 = new TText(240, 15, "Best fit for NuMI beam spill window");
  tx4->SetTextFont(32);
  tx4->SetTextSize(0.035);
  tx4->SetTextColor(kRed);
  tx4->Draw();
  */

  cana_tave_wide->SaveAs("cana_tave_wide.pdf");
  cana_tave_wide->SaveAs("cana_tave_wide.png");
}

void time_wide_proj

(

TChain &

ch

)

Definition at line 350 of file cana_time_peak.C.

References canacuts(), and tavewideh.

Referenced by proj_all().

{
  tavewideh = new TH1F("tavewideh",
                       ";NDOS event time [#musec];Events / 2 #musec",
                       250,0,500);
  ch.Project("tavewideh", "tave", canacuts);
}

void time_zoom_plot

(

)

Definition at line 632 of file cana_time_peak.C.

References CenterTitles(), NDOSLabel(), Preliminary(), tavezoomh, and time_profileh.

Referenced by plot_all().

{
  TCanvas* time_zoom_drawc = 
    new TCanvas("time_zoom_drawc",
                "time_zoom_drawc",800,600);
  tavezoomh->SetFillColor(kGreen-10);
  CenterTitles(tavezoomh);
  tavezoomh->SetLineWidth(1);
  tavezoomh->Draw();
  Preliminary();

  if (time_profileh) {
    time_profileh->SetLineWidth(2);
    time_profileh->SetLineColor(kRed);
    time_profileh->Draw("same");
  }

  TLegend *tleg = new TLegend(0.15,0.75,0.43,0.85);
  tleg->SetTextFont(42);
  tleg->SetTextSize(0.04);
  tleg->AddEntry(time_profileh,"Fit to NuMI Beam Structure");
  tleg->Draw();
  
  NDOSLabel();

  time_zoom_drawc->SaveAs("cana_time_zoom.png");
  time_zoom_drawc->SaveAs("cana_time_zoom.pdf");
}

void time_zoom_proj

(

TChain &

ch

)

Definition at line 401 of file cana_time_peak.C.

References canacuts(), tavezoomh, zoom_nbin, zoom_thi, and zoom_tlo.

Referenced by proj_all().

{  
  gStyle->SetOptStat(0);
  
  tavezoomh = new TH1F("tavezoomh",
                       ";NDOS event time [#musec];Events / 50 nsec",
                       zoom_nbin,
                       zoom_tlo,
                       zoom_thi);
  ch.Project("tavezoomh", "tave", canacuts);
}

void x1_plot

(

)

Definition at line 185 of file cana_time_peak.C.

References CenterTitles(), NDOSLabel(), Preliminary(), and x1h.

Referenced by plot_all().

{
  gStyle->SetOptStat(0);
  
  TCanvas* x1c = new TCanvas("x1c","x1c",800,600);
  
  x1h->SetFillColor(kGreen-10);
  CenterTitles(x1h);
  x1h->Draw();
  Preliminary();
  
  TLine* x1l1 = new TLine(-100.0, 0, -100.0, x1h->GetMaximum());
  TLine* x1l2 = new TLine( 100.0, 0,  100.0, x1h->GetMaximum());
  x1l1->SetLineStyle(3);
  x1l2->SetLineStyle(3);
  x1l1->Draw();
  x1l2->Draw();
  double ya1 = 0.66*x1h->GetMaximum();
  double ya2 = 0.60*x1h->GetMaximum();
  double ya3 = 0.5*(ya1+ya2);
  
  TArrow* x1a1 = new TArrow(-100.0, ya1, -80.0, ya1, 0.02);
  TArrow* x1a2 = new TArrow(-100.0, ya2, -80.0, ya2, 0.02);
  x1a1->Draw();
  x1a2->Draw();

  TText* tx2 = new TText(0.0,ya3,"selected");
  tx2->SetTextFont(42);
  tx2->SetTextAlign(22);
  tx2->SetTextSize(0.035);
  tx2->Draw();

  TArrow* x1a3 = new TArrow(100.0, ya1, 80.0, ya1, 0.025);
  TArrow* x1a4 = new TArrow(100.0, ya2, 80.0, ya2, 0.025);
  x1a3->Draw();
  x1a4->Draw();

  NDOSLabel();

  x1c->SaveAs("cana_x1c.pdf");
  x1c->SaveAs("cana_x1c.png");
}

void x1_proj

(

TChain &

ch

)

Definition at line 179 of file cana_time_peak.C.

References activity(), x1h, and y1cut().

Referenced by proj_all().

{
  x1h = new TH1F("x1h",";upper end x location (cm); Events / cm",
                 500,-250,250);
  ch.Project("x1h","x1", run&&activity&&y1cut);
}

TCut x1cut

(

)

Referenced by z1_proj().

void y1_plot

(

)

Definition at line 141 of file cana_time_peak.C.

References CenterTitles(), NDOSLabel(), Preliminary(), and y1h.

Referenced by plot_all().

{
  gStyle->SetOptStat(0);

  TCanvas* y1c = new TCanvas("y1c","y1c",800,600);
  
  y1h->SetFillColor(kGreen-10);
  CenterTitles(y1h);
  y1h->Draw();
  Preliminary();

  TLine* y1l = new TLine(150.0, 0, 150.0, y1h->GetMaximum());
  y1l->SetLineStyle(3);
  y1l->Draw();
  double ya1 = 0.66*y1h->GetMaximum();
  double ya2 = 0.60*y1h->GetMaximum();
  double ya3 = 0.5*(ya1+ya2);
  
  TArrow* y1a1 = new TArrow(150.0, ya1, 120.0, ya1, 0.03);
  TArrow* y1a2 = new TArrow(150.0, ya2, 120.0, ya2, 0.03);
  y1a1->Draw();
  y1a2->Draw();

  NDOSLabel();

  TText* tx2 = new TText(110.0,ya3,"selected");
  tx2->SetTextAlign(32);
  tx2->SetTextFont(42);
  tx2->SetTextSize(0.035);
  tx2->Draw();

  y1c->SaveAs("cana_y1c.pdf");
  y1c->SaveAs("cana_y1c.png");
}

void y1_proj

(

TChain &

ch

)

Definition at line 135 of file cana_time_peak.C.

References activity(), and y1h.

Referenced by proj_all().

{
  y1h = new TH1F("y1h",";upper end y location (cm); Events / cm",
                       500,-250,250);
  ch.Project("y1h","y1", run&&activity);
}

TCut y1cut

(

)

Referenced by x1_proj(), and z1_proj().

void z1_plot

(

)

Definition at line 237 of file cana_time_peak.C.

References CenterTitles(), NDOSLabel(), Preliminary(), and z1h.

Referenced by plot_all().

{
  gStyle->SetOptStat(0);

  TCanvas* z1c = new TCanvas("z1c","z1c",800,600);

  z1h->SetFillColor(kGreen-10);
  CenterTitles(z1h);
  z1h->Draw();
  Preliminary();

  TLine* z1l1 = new TLine( 280.0, 0,  280.0, z1h->GetMaximum());
  TLine* z1l2 = new TLine(1200.0, 0, 1200.0, z1h->GetMaximum());
  z1l1->SetLineStyle(3);
  z1l2->SetLineStyle(3);
  z1l1->Draw();
  z1l2->Draw();
  double ya1 = 0.66*z1h->GetMaximum();
  double ya2 = 0.60*z1h->GetMaximum();
  double ya3 = 0.5*(ya1+ya2);
  
  TArrow* z1a1 = new TArrow(280.0, ya1, 330.0, ya1, 0.025);
  TArrow* z1a2 = new TArrow(280.0, ya2, 330.0, ya2, 0.025);
  z1a1->Draw();
  z1a2->Draw();

  TText* tx2 = new TText(700.0,ya3,"selected");
  tx2->SetTextFont(42);
  tx2->SetTextAlign(22);
  tx2->SetTextSize(0.035);
  tx2->Draw();

  TArrow* z1a3 = new TArrow(1200.0, ya1, 1150.0, ya1, 0.025);
  TArrow* z1a4 = new TArrow(1200.0, ya2, 1150.0, ya2, 0.025);
  z1a3->Draw();
  z1a4->Draw();

  NDOSLabel();

  z1c->SaveAs("cana_z1c.pdf");
  z1c->SaveAs("cana_z1c.png");
}

void z1_proj

(

TChain &

ch

)

Definition at line 231 of file cana_time_peak.C.

References activity(), x1cut(), y1cut(), and z1h.

Referenced by proj_all().

{
  z1h = new TH1F("z1h",";upper end z location (cm); Events / 10 cm",
                 140,0,1400);
  ch.Project("z1h","z1",run&&activity&&y1cut&&x1cut);
}

TCut z1cut

(

)

Variable Documentation

TH1F* cosnumiinh

Definition at line 282 of file cana_time_peak.C.

Referenced by cosnumi_plot(), and cosnumi_proj().

TH1F* cosnumiouth

Definition at line 283 of file cana_time_peak.C.

Referenced by cosnumi_plot(), and cosnumi_proj().

const double dt_batch = 81*18.83E-3

static

Definition at line 419 of file cana_time_peak.C.

Referenced by spill_profile().

const double dt_cycle = dt_batch + dt_gap

static

Definition at line 421 of file cana_time_peak.C.

Referenced by spill_profile().

const double dt_gap = 5*18.83E-3

static

Definition at line 420 of file cana_time_peak.C.

double intimet = (t2-t1)

Definition at line 53 of file cana_time_peak.C.

Referenced by cosnumi_proj().

const char* label = "NDOS Data, 09/12/13 - 09/27/13"

Definition at line 34 of file cana_time_peak.C.

Referenced by cvn::RegPixelMap::Add(), cvn::PixelMap::Add(), Ana2017_sb_opening_macro(), AnaResultsSpectra(), cheat::TestTrackIds::beginJob(), dprf::ChannelPlots::beginRun(), BestFitToDataRelease(), BicountEllipse_dCP(), bin_composition_pie_chart(), BiprobEllipse_dCP(), dbutils::gui::BitmapEditor::BitmapEditor(), caf_numu_fd_validation_data(), caf_numu_fd_validation_MC(), caf_numu_fd_validation_MC_no_tau(), CalculateXSec(), dbutils::gui::ConfigSelectionDialog::ConfigSelectionDialog(), ana::CountingExperimentErrorBarChart(), cvn::RegPixelMapProducer::CreateMapGivenBoundary(), cvn::PixelMapProducer::CreateMapGivenBoundary(), cvn::RegPixelMapProducer::CreateMapGivenShowerVertex(), cvn::RegPixelMapProducer::CreateMapGivenVertex(), data_over_mc_profile(), DataVtxDistributionAna(), DataVtxDistributionLoad(), demoFitContours(), demoFitSlices(), DirNameToProbe(), DirNameToTarget(), DrawDeltaChiSqrWithContours(), DrawHCLabel(), drawLongTerm(), DrawResBin(), drawTimePlots(), drawVsPOT(), efficiency(), efficiencySA(), skim::CellHitSkimmer::EvaluateSlices(), ExtendedAxesAna(), fake_future_data(), FD_plots(), FitParamEffectsAna(), FitSystEffectsAna(), genie_contours(), getBatch(), getBNBPlots(), art::DataViewImpl::getByLabel(), GetGENIEShiftLabels(), getTimePeakPlots(), GetXSecShiftLabels(), goodness_of_fit(), GraphsToDataRelease(), dprf::ChannelPlots::HRFebToStringLabel(), rb::IsFiltered(), joint_fit_2017_slices(), joint_fit_2018_contours(), joint_fit_2018_slices(), joint_fit_2019_bestfit(), joint_fit_2019_contours(), joint_fit_2019_slices(), joint_fit_future_bestfit_univ(), joint_fit_future_contour_univ(), jointsensitivity(), CalibAnaPlot::LegendLabel(), ana::NumuSyst::LoadFrom(), ana::NCSyst::LoadFrom(), ana::LoadFrom(), ana::NueSyst::LoadFrom(), ana::FluxReweight::LoadFrom(), ana::NDOscillatableSpectrum::LoadFrom(), ana::ReweightableSpectrum::LoadFrom(), ana::OscillatableSpectrum::LoadFrom(), ana::SystMakerShift::LoadFrom(), ana::GenericSystematicDef< SRType >::LoadFrom(), ana::TruthReweight::LoadFrom(), ana::RecoReweight::LoadFrom(), ana::SystMaker::LoadFrom(), ana::Spectrum::LoadFrom(), ana::LoadFromFile(), make_pi0_xcheck(), make_pid(), calib::MakeAttenuationProfiles::MakeAttenuationProfiles(), MakePlots(), MakeSummaryTable(), NDOSLabel(), cvn::RegPixelMap::NInput(), cvn::PixelMap::NInput(), nuint09_1pi1(), nuint09_1pi2(), nuint09_1pi3(), nuint09_1pi4(), nuint09_coh1(), nuint09_coh2(), nuint09_coh3(), nuint09_qel1(), nuint09_qel2(), nuint09_qel3(), nuint09_qel4(), nuint09_qel5(), nuint09_qel6(), genie::flux::operator<<(), ana::CutOptimization::OptimizedSigmaOverSigma(), evdb::ParameterSetEdit::ParameterSetEdit(), pi0_xcheck(), plot_joint_fit_2020_slices(), plot_pid(), PlotEvents(), ana::PlotPurEff(), ana::PlotSpectra(), fxwgt::FluxWeightCalculator::produce(), reco_minus_true_panels(), skim::CellHitSkimmer::reconfigure(), run_joint_fit_2020_bestfit(), run_joint_fit_2020_contours(), run_joint_fit_2020_slices(), ana::SaveTo(), sensitivity2018(), sensitivity2020(), set_labels_1D(), set_labels_2D(), ana::SpectrumHandler::SetLoader(), evd::SliceButtonBar::SliceButtonBar(), SpectrumParamEffectsAna(), TheoryEllipse_dCP(), timingPeak(), Tutorial2019Fit(), Tutorial2019FitContours(), Tutorial2019FitSlices(), om::WatchListBox::UpdateCurrentList(), validationscript(), cvn::CVNEventDump::WriteMapTH2(), cvn::CVNCosmicEventDump::WriteMapTH2(), calib::MakeAttenuationProfiles::writeResults(), FDPredictionHandler::~FDPredictionHandler(), and NDPredictionHandler::~NDPredictionHandler().

const int nbatch = 6

static

Definition at line 418 of file cana_time_peak.C.

TH1F* nhith

Definition at line 87 of file cana_time_peak.C.

Referenced by make_nhitc(), nhit_plot(), and nhit_proj().

int ntbin_per_usec = 20

Definition at line 396 of file cana_time_peak.C.

double outtimet = 400.0

Definition at line 54 of file cana_time_peak.C.

Referenced by cosnumi_proj().

const double spill_length = nbatch*dt_batch + (nbatch-1)*dt_gap

static

Definition at line 422 of file cana_time_peak.C.

Referenced by make_template(), and spill_profile().

double t1 = 217.0

Definition at line 49 of file cana_time_peak.C.

Referenced by time_wide_plot().

double t2 = 227.8

Definition at line 50 of file cana_time_peak.C.

Referenced by art::detail::CanBeAggregated< std::array< T, N > >::aggregate(), singlekaon_xsec::Amatrix_NN(), genie::AlamSimoAtharVacasSKPXSec2014::Amatrix_NN(), singlekaon_xsec::Amatrix_NP(), genie::AlamSimoAtharVacasSKPXSec2014::Amatrix_NP(), singlekaon_xsec::Amatrix_PP(), genie::AlamSimoAtharVacasSKPXSec2014::Amatrix_PP(), zcl::SMMTriggerAna::analyze(), beamlinereco::ToFPositionRecoAnalysis::analyze(), beamlinereco::LEHitFinder< T >::BackwardFindingOfHitStart(), beamlinereco::CFDHitFinder< T >::BackwardFindingOfHitStart(), genie::BBA05ELFormFactorsModel::BBA05Fit(), novaddt::HoughPoint::calc(), novaddt::SlowMMTrigger::filter(), beamlinereco::CFDHitFinder< T >::FindHitWidth(), genie::HAIntranuke::Inelastic(), genie::HAIntranuke2018::Inelastic(), shadow::NOVASLocate::JulianDate(), runh::BadDataFilter::LoadBadTimeRanges(), main(), evtsum::EventSummary::MakeOutput(), genie::StrumiaVissaniIBDPXSec::MtxElm(), osc::OscCalc::P_internal_me(), osc::OscCalc::P_internal_mt(), osc::OscCalc::P_internal_te(), art::pointersEqual(), zcl::SPCluster::produce(), nl::TimeSeriesData::RangeTime(), locator::RndmLocatorTable::RndmLocatorTable(), evd::GeometryDrawer::ShowCellAndPlaneNumbers(), osc::SolveCubic(), sn::SortTrackEndsByTime(), and time_wide_plot().

TH1F* tavewideh

Definition at line 349 of file cana_time_peak.C.

Referenced by time_wide_plot(), and time_wide_proj().

TH1F* tavezoomh

Definition at line 395 of file cana_time_peak.C.

Referenced by chi2(), minimize(), time_zoom_plot(), and time_zoom_proj().

int template_nbin

Definition at line 447 of file cana_time_peak.C.

Referenced by make_template(), and make_time_profile().

TH1F* template_spillh = 0

Definition at line 444 of file cana_time_peak.C.

Referenced by make_template(), and make_time_profile().

double template_thi

Definition at line 446 of file cana_time_peak.C.

Referenced by make_template().

double template_tlo

Definition at line 445 of file cana_time_peak.C.

Referenced by make_template().

TH1F* time_profileh

Definition at line 487 of file cana_time_peak.C.

Referenced by chi2(), make_time_profile(), and time_zoom_plot().

TH1F* x1h

Definition at line 178 of file cana_time_peak.C.

Referenced by x1_plot(), and x1_proj().

TH1F* y1h

Definition at line 134 of file cana_time_peak.C.

Referenced by y1_plot(), and y1_proj().

TH1F* z1h

Definition at line 230 of file cana_time_peak.C.

Referenced by z1_plot(), and z1_proj().

int zoom_nbin = ntbin_per_usec*(zoom_thi-zoom_tlo)

Definition at line 400 of file cana_time_peak.C.

Referenced by chi2(), make_time_profile(), minimize(), and time_zoom_proj().

double zoom_tbin_sz = 1.0/(float)ntbin_per_usec

Definition at line 397 of file cana_time_peak.C.

Referenced by make_template().

double zoom_thi = floor(t2) + 6

Definition at line 399 of file cana_time_peak.C.

Referenced by make_time_profile(), and time_zoom_proj().

double zoom_tlo = floor(t1) - 5

Definition at line 398 of file cana_time_peak.C.

Referenced by make_time_profile(), and time_zoom_proj().