21 "$NOVA_PROD/mc/FA14-10-28/genie/fd/caf/010000/1000000/*_fhc_nonswap_*";
23 "$NOVA_PROD/mc/FA14-10-28/genie/fd/caf/010000/1000000/*_fhc_swap_*";
25 "$NOVA_PROD/mc/FA14-10-28/genie/fd/caf/010000/1000000/*_rhc_nonswap_*";
27 "$NOVA_PROD/mc/FA14-10-28/genie/fd/caf/010000/1000000/*_rhc_swap_*";
54 const double pot = 18e20;
55 const double potrhc = 18e20;
65 calcR.
SetTh23(TMath::Pi()*.2141);
106 std::vector<const IExperiment*> exptsR;
108 exptsR.push_back(&exptCCR);
109 exptsR.push_back(&exptQER);
110 exptsR.push_back(&exptCCRHCR);
111 exptsR.push_back(&exptQERHCR);
121 std::vector<const IExperiment*> exptsB;
123 exptsB.push_back(&exptCCB);
124 exptsB.push_back(&exptQEB);
125 exptsB.push_back(&exptCCRHCB);
126 exptsB.push_back(&exptQERHCB);
131 TGaxis::SetMaxDigits(3);
135 const double xmin = -.15;
136 const double xmax = .25;
137 const double ymin = -1.5e-3;
138 const double ymax = 1
e-3;
140 TCanvas*
c1 =
new TCanvas(
"c1",
"cptPlot", 1024, 768);
151 &kFitDeltaCPTSinSq2Theta23B, res, xmin, xmax,
152 &kFitDeltaCPTDmSq32NormB, res, ymin, ymax,
153 {&kFitSigmaCPTSinSq2Theta23B,&kFitSigmaCPTDmSq32NormB}
162 TLine *
xaxis =
new TLine(xmin,0,xmax,0);
163 xaxis->SetLineColor(kGray);
164 xaxis->SetLineStyle(7);
167 TLine *
yaxis =
new TLine(0,ymin,0,ymax);
168 yaxis->SetLineColor(kGray);
169 yaxis->SetLineStyle(7);
181 &kFitDeltaCPTSinSq2Theta23R, res, xmin, xmax,
182 &kFitDeltaCPTDmSq32NormR, res, ymin, ymax,
183 {&kFitSigmaCPTSinSq2Theta23R,&kFitSigmaCPTDmSq32NormR}
192 c1->SaveAs(
"CPT-sesitivity.pdf");
void SetDmsq21(const double &dmsq21) override
Cuts and Vars for the 2020 FD DiF Study.
std::map< std::string, double > xmax
void SetDmsq32(const double &dmsq32) override
TH2 * Gaussian68Percent2D(const FrequentistSurface &s)
Up-value surface for 68% confidence in 2D in gaussian approximation.
TH2 * Gaussian2Sigma2D(const FrequentistSurface &s)
Up-value surface for 2 sigma confidence in 2D in gaussian approximation.
const FitSinSq2Theta23CPT kFitSinSq2Theta23bar(osc::ENuSign::kNuBar)
Representation of a spectrum in any variable, with associated POT.
void SetTh12(const double &th12) override
Log-likelihood scan across two parameters.
void SetTh13(const double &th13) override
void DrawBestFit(Color_t color, Int_t marker=kFullCircle) const
Draw the best fit point.
void SetRho(double rho) override
const Binning kNumuEnergyBinning
Spectrum FakeData(double pot) const
Synonymous with AsimovData(). Retained for compatibility.
const FitDmSq32CPTHierarchy kFitDmSq32nuNorm(osc::ENuSign::kNu, false)
TH2 * Gaussian3Sigma2D(const FrequentistSurface &s)
Up-value surface for 3 sigma confidence in 2D in gaussian approximation.
const FitSinSq2Theta23CPT kFitSinSq2Theta23nu(osc::ENuSign::kNu)
virtual void Go() override
Load all the registered spectra.
void SetTh23(const double &th23) override
Combine multiple component experiments.
Collaborates with Spectrum and OscillatableSpectrum to fill spectra from CAF files.
void SetL(double L) override
const Var kQEE
Energy estimator for quasielastic CC events.
const FitDmSq32CPTHierarchy kFitDmSq32barNorm(osc::ENuSign::kNuBar, false)
Compare a single data spectrum to the MC + cosmics expectation.