20 #include "CAFAna/Core/Binning.h" 29 #include "TMultiGraph.h" 34 #include "TPaveText.h" 43 #include "Utilities/func/MathUtil.h" 47 #define location 0 // 0 for FNAL 1 for caltech 48 #define mcset 1 // 0 = ideal MC goes with pot like 18e20, or 1 = realistic conditions mc for 1st analysis POT 55 TLatex*
prelim =
new TLatex(.9, .95,
"NO#nuA Simulation");
56 prelim->SetTextColor(kGray+1);
58 prelim->SetTextSize(2/30.);
59 prelim->SetTextAlign(32);
66 TGaxis::SetMaxDigits(2);
68 gStyle->SetTitleOffset(.85,
"y");
69 gStyle->SetTitleOffset(.84,
"x");
70 gStyle->SetTitleSize(0.05,
"x");
77 fname =
"$NOVA_DATA/mc/forFA/fdgenie.nonswap.decaf.root";
78 fnameSwap =
"$NOVA_DATA/mc/forFA/fdgenie.swap.decaf.root";
81 fname =
"prod_decaf_S15-05-22a_fd_genie_fhc_nonswap_numu_contain";
82 fnameSwap =
"prod_decaf_S15-05-22a_fd_genie_fhc_fluxswap_numu_contain";
89 fname =
"prod_decaf_S15-05-22a_fd_genie_fhc_nonswap_ideal_numu_contain";
90 fnameSwap =
"prod_decaf_S15-05-22a_fd_genie_fhc_fluxswap_ideal_numu_contain";
93 fname =
"prod_decaf_S15-05-22a_fd_genie_fhc_nonswap_ideal_numu_contain";
94 fnameSwap =
"prod_decaf_S15-05-22a_fd_genie_fhc_fluxswap_ideal_numu_contain";
130 text.SetTextSize(0.04);
139 new TCanvas(
"comp1b",
"CC spectra .95");
142 ccu2->SetLineStyle(7);
143 ccu2->Draw(
"hist same");
144 gPad->Print(
"plots/comp1b.ps");
148 TH1F *hcc =
new TH1F(
"cc",
"cc", 10, 0, 1);
149 TH1F *hqe =
new TH1F(
"qe",
"qe", 10, 0, 1);
150 TH1F *huc =
new TH1F(
"uc",
"uc", 10, 0, 1);
151 TH1F *hall =
new TH1F(
"all",
"all", 10, 0, 1);
153 hcc->SetLineWidth(2);
154 hqe->SetLineWidth(2);
155 huc->SetLineWidth(2);
156 hall->SetLineWidth(2);
157 hcc->SetLineColor(2);
158 hqe->SetLineColor(4);
159 huc->SetLineColor(kMagenta);
160 hall->SetLineColor(1);
163 new TCanvas(
"contoursfhc1",
"contoursfhc1") ;
173 TLegend *
leg =
new TLegend(0.11,0.11,0.42,0.35);
174 leg->SetFillColor(10);
175 leg->SetFillStyle(0);
176 leg->SetBorderSize(0);
177 if(
pot==18e20) leg->SetHeader(
"18e20 POT, FHC, 14 kton: 90% CL");
178 else if(
pot==6e20) leg->SetHeader(
"6e20 POT, FHC, 14 kton: 90% CL");
179 else if(
pot==2e20) leg->SetHeader(
"2e20 POT, FHC, 14 kton: 90% CL");
180 else if(
pot==1e20) leg->SetHeader(
"1e20 POT, FHC, 14 kton: 90% CL");
182 leg->AddEntry(hall,
"Contained CC",
"L");
184 text.DrawLatex(0.3,0.85,
"NO#nuA #nu_{#mu} Sensitivity (14kton, 700kW)");
186 gPad->Print(
"plots/basic.1.ps");
187 gPad->Print(
"plots/basic.1.C");
191 new TCanvas(
"contoursfhc2",
"contoursfhc2") ;
200 text.DrawLatex(0.3,0.85,
"NO#nuA #nu_{#mu} Sensitivity (14kton, 700kW)");
204 gPad->Print(
"plots/basic.fhc.2.ps");
205 gPad->Print(
"plots/basic.fhc.2.C");
207 new TCanvas(
"contourscombine",
"contours90percent") ;
216 TH1F *solid =
new TH1F(
"90% CL",
"90% CL", 10, 0, 1);
217 TH1F *dashed =
new TH1F(
"68% CL",
"68% CL", 10, 0, 1);
219 solid->SetLineWidth(2);
220 dashed->SetLineWidth(2);
221 dashed->SetLineStyle(7);
222 solid->SetLineColor(1);
223 dashed->SetLineColor(1);
225 TLegend *leg4 =
new TLegend(0.12,0.12,0.5,0.4);
226 leg4->SetFillColor(10);
227 leg4->SetFillStyle(0);
228 leg4->SetHeader(
"18e20 POT FHC, 14 kton");
230 leg4->AddEntry(solid,
"2 #sigma",
"L");
231 leg4->AddEntry(dashed,
"1 #sigma",
"L");
232 leg4->SetBorderSize(0);
234 text.DrawLatex(0.3,0.92,
"NO#nuA #nu_{#mu} Sensitivity (14kton, 700kW)");
236 gPad->Print(
"plots/basic.overlay.ps");
237 gPad->Print(
"plots/basic.overlay.C");
void SetDmsq21(const T &dmsq21) override
Cuts and Vars for the 2020 FD DiF Study.
void SetTh12(const T &th12) override
void SetTh23(const T &th23) override
TH1D * ToTH1(double exposure, Color_t col=kBlack, Style_t style=kSolid, EExposureType expotype=kPOT, EBinType bintype=kBinContent) const
Histogram made from this Spectrum, scaled to some exposure.
Adapt the PMNS calculator to standard interface.
const FitSinSq2Theta23 kFitSinSq2Theta23
void SetL(double L) override
TH2 * Gaussian68Percent2D(const FrequentistSurface &s)
Up-value surface for 68% confidence in 2D in gaussian approximation.
virtual Spectrum PredictUnoscillated() const
TH2 * Gaussian2Sigma2D(const FrequentistSurface &s)
Up-value surface for 2 sigma confidence in 2D in gaussian approximation.
Representation of a spectrum in any variable, with associated POT.
Log-likelihood scan across two parameters.
void DrawBestFit(Color_t color, Int_t marker=kFullCircle) const
Draw the best fit point.
const Binning kNumuEnergyBinning
Spectrum FakeData(double pot) const
Synonymous with AsimovData(). Retained for compatibility.
void SetdCP(const T &dCP) override
void DrawContour(TH2 *fc, Style_t style, Color_t color, double minchi=-1)
void SetDmsq32(const T &dmsq32) override
virtual void Go() override
Load all the registered spectra.
void SetTh13(const T &th13) override
const ConstrainedFitVarWithPrior fitDmSq32Scaled_UniformPrior & kFitDmSq32Scaled
void SetRho(double rho) override
Collaborates with Spectrum and OscillatableSpectrum to fill spectra from CAF files.
TH2 * Gaussian90Percent2D(const FrequentistSurface &s)
Up-value surface for 90% confidence in 2D in gaussian approximation.
TH1 * DataMCComparison(const Spectrum &data, const Spectrum &mc, EBinType bintype)
Compare a single data spectrum to the MC + cosmics expectation.