12 #include "CAFAna/Analysis/CovMxSurface.h" 14 #include "CAFAna/covmx/Utilities.h" 33 std::vector<covmx::Sample> samples;
34 std::vector<bool> unoscillated(samples.size(),
false);
47 for (
size_t i = 0;
i < samples.size(); ++
i) {
48 std::cout <<
"Sample: " << samples[
i].GetTag() <<
", disabling oscillations? " << unoscillated[
i] <<
std::endl;
62 std::vector<IPrediction*> preds;
66 std::vector<Spectrum*>
cosmic;
67 for (
auto sample : samples)
71 std::vector<Spectrum*>
data;
77 for (
size_t i = 0;
i < preds.size(); ++
i) {
81 data.back()->OverrideLivetime(livetimeFD);
82 if (cosmic[
i]) *data[
i] += *cosmic[
i];
98 for (
int i = 2;
i <= 4; ++
i) calc->
SetDm(
i, 0);
99 expt->DisableOscillations(unoscillated, noosccalc);
101 std::vector<const IChiSExperiment*> expts;
103 expts.push_back(expt);
107 std::vector<const IFitVar*> fitVars;
108 std::map<const IFitVar*, std::vector<double> > seedValues;
109 std::vector<double> fit_seeds;
118 TH2D* varHists[fitVars.size()];
119 unsigned int padWidth = (*std::max_element(fitVars.begin(), fitVars.end(),
121 return lhs->
ShortName().length() < rhs->ShortName().length();
122 }))->ShortName().length();
142 TFile*
f = TFile::Open(
"surface_joint.root",
"recreate");
143 surf.SaveTo(f->mkdir(
"surface"));
size_t JobNumber()
What's the process number for a grid job?
Represent the binning of a Spectrum's x-axis.
Cuts and Vars for the 2020 FD DiF Study.
void SetNus20Params(osc::OscCalcSterile *calc, std::string type="3flav")
bool RunningOnGrid()
Is this a grid (condor) job?
A simple Gaussian constraint on an arbitrary IFitVar.
Adapt the PMNS_Sterile calculator to standard interface.
const FitDmSq32Sterile kFitDmSq32Sterile
const double kAna2018SensitivityFHCNDPOT
static SystShifts Nominal()
const XML_Char const XML_Char * data
osc::OscCalcSterile * DefaultSterileCalc(int nflavors)
Create a sterile calculator with default assumptions for all parameters.
Sum up livetimes from individual cosmic triggers.
const FitSinSq2Theta13Sterile kFitSinSq2Theta13Sterile
Compare a single data spectrum to the MC + cosmics expectation.
const FitDelta13InPiUnitsSterile kFitDelta13InPiUnitsSterile
std::vector< float > Spectrum
Combine multiple component experiments.
void SetDm(int i, double dm)
const std::string & ShortName() const
IPrediction * LoadPrediction(std::string detector, bool rhc=false, std::string syst_type="all")
Function to load prediction object.
const FitSinSqTheta23Sterile kFitSinSqTheta23Sterile
std::unique_ptr< Spectrum > LoadCosmic(covmx::Sample sample, bool cvmfs=true)
Get cosmics for a given sample.
Interface definition for fittable variables.
const double kAna2018FHCPOT
Prevent histograms being added to the current directory.
static Binning Simple(int n, double lo, double hi, const std::vector< std::string > &labels={})
const double kAna2018FHCLivetime