25 for(
int i = 0;
i < 1000; ++
i){
26 const double vx = gRandom->Uniform(0, 1);
29 const double vy = gRandom->Uniform(0, 1);
35 for(
int i = 0;
i < 1000; ++
i){
36 const double va = gRandom->Uniform(0, 1);
37 alpha.SetValue(&osc, va);
39 const double vb = gRandom->Uniform(0, 1);
40 beta.SetValue(&osc, vb);
46 TH2*
axes =
new TH2F(
"",
"", 120, -.1, 1.1, 120, -.1, 1.1);
47 axes->GetXaxis()->SetTitle(alpha.LatexName().c_str());
48 axes->GetYaxis()->SetTitle(
beta.LatexName().c_str());
50 for(
int ix = 0; ix <= 20; ++ix){
51 TGraph*
g =
new TGraph;
52 for(
int iy = 0; iy <= 100; ++iy){
55 g->SetPoint(iy, alpha.GetValue(&osc),
beta.GetValue(&osc));
61 for(
int iy = 0; iy <= 20; ++iy){
62 TGraph*
g =
new TGraph;
63 for(
int ix = 0; ix <= 100; ++ix){
66 g->SetPoint(ix, alpha.GetValue(&osc),
beta.GetValue(&osc));
75 axes =
new TH2F(
"",
"", 120, -.1, 1.1, 120, -.1, 1.1);
76 axes->GetXaxis()->SetTitle(x.
LatexName().c_str());
77 axes->GetYaxis()->SetTitle(y.
LatexName().c_str());
80 for(
int ia = 0; ia <= 20; ++ia){
81 TGraph*
g =
new TGraph;
82 for(
int ib = 0; ib <= 100; ++ib){
83 alpha.SetValue(&osc, ia*.05);
84 beta.SetValue(&osc, ib*.01);
91 for(
int ib = 0; ib <= 20; ++ib){
92 TGraph*
g =
new TGraph;
93 for(
int ia = 0; ia <= 100; ++ia){
94 alpha.SetValue(&osc, ia*.01);
95 beta.SetValue(&osc, ib*.05);
const std::string & LatexName() const
Cuts and Vars for the 2020 FD DiF Study.
fvar< T > fabs(const fvar< T > &x)
FitVarsProductMarg FitVarSinSqTheta23SinSq2Theta13Marg()
Return a FitVarsProductMarg to marginalize over this if using kFitSinSqTheta23SinSq2Theta13.
void test_fitvarsproduct()
Optimized version of OscCalcPMNS.
FitVarsProduct FitVarSinSqTheta23SinSq2Theta13()
Return a FitVarsProduct for .
Oscillation probability calculators.
virtual void SetValue(osc::IOscCalcAdjustable *osc, double val) const =0
const ConstrainedFitVarWithPrior fitSsqTh23_UniformPriorSsqTh23 & kFitSinSqTheta23
assert(nhit_max >=nhit_nbins)
Interface definition for fittable variables.
const FitSinSq2Theta13 kFitSinSq2Theta13
virtual double GetValue(const osc::IOscCalcAdjustable *osc) const =0