Oscillation analysis framework, runs over CAF files outside of ART. More...
Namespaces | |
CosmicIdUtils | |
Current | |
Enumeration for interaction currents (CC/NC) | |
ExampleAnalysis | |
Flavors | |
Enumeration of neutrino transition modes. | |
SBNOsc | |
Sign | |
Enumeration for neutrino sign (neutrino/antineutrino) | |
Classes | |
class | TPCPurityInfoAna |
class | NuFitPenalizer |
class | Penalizer_GlbLike |
class | NuFitPenalizerCDR |
class | Penalizer_GlbLikeCDR |
class | Fitter |
Perform MINUIT fits in one or two dimensions. More... | |
class | FitAxis |
Collect information describing the axis of a fit variable. More... | |
class | GradientDescent |
class | MedianSurface |
class | OscCalcNoHash |
Helper for Surface::FillSurfacePoint. More... | |
class | Surface |
Log-likelihood scan across two parameters. More... | |
class | Binning |
Represent the binning of a Spectrum's x-axis. More... | |
class | _Cut |
Template for Cut and SpillCut. More... | |
class | InstallHandlers |
class | EnsembleRatio |
class | EnsembleSpectrum |
class | FileListSource |
Simple file source based on an explicit list provided by the user. More... | |
class | FileReducer |
Create smaller CAFs. More... | |
class | _HistAxis |
Collect information describing the x-axis of an analysis histogram. More... | |
class | HistCache |
Helper for Spectrum. More... | |
class | IFileSource |
Interface class for accessing ROOT files in sequence. More... | |
class | IFitVar |
Interface definition for fittable variables. More... | |
class | IConstrainedFitVar |
Base class for variables with constraints. Apply penalty outside range. More... | |
class | ISyst |
Encapsulate code to systematically shift a caf::StandardRecord. More... | |
class | Loaders |
Collection of SpectrumLoaders for many configurations. More... | |
class | MultiVar2DFunc |
Helper for MultiVar2D. More... | |
class | _MultiVar |
class | OscCalcSterileApprox |
class | OscCalcSterileApproxAdjustable |
class | OscCurve |
Transition probability for any one channel as a function of energy. More... | |
class | OscillatableSpectrum |
Spectrum with true energy information, allowing it to be oscillated More... | |
class | ProfilerSupport |
Support for the –prof commandline option. More... | |
class | Progress |
A simple ascii-art progress bar. More... | |
class | Ratio |
Represent the ratio between two spectra. More... | |
class | ReweightableSpectrum |
Spectrum with the value of a second variable, allowing for reweighting More... | |
class | SAMProjectSource |
Fetch files from a pre-existing SAM project. More... | |
class | SAMQuerySource |
File source based on a SAM query or dataset (definition) More... | |
class | Spectrum |
Representation of a spectrum in any variable, with associated POT. More... | |
struct | CompareByID |
class | CutVarCache |
Helper for HandleRecord. More... | |
class | SpectrumLoader |
Collaborates with Spectrum and OscillatableSpectrum to fill spectra from CAF files. More... | |
class | SpectrumLoaderBase |
Base class for the various types of spectrum loader. More... | |
class | NullLoader |
Dummy loader that doesn't load any files. More... | |
class | SystRegistry |
class | SystShifts |
Simple record of shifts applied to systematic parameters. More... | |
class | ThreadPool |
A very simple thread pool for use by Surface. More... | |
class | DontAddDirectory |
Prevent histograms being added to the current directory. More... | |
class | IFDHSilent |
ifdh calls between construction and destruction produce no output More... | |
class | FloatingExceptionOnNaN |
Alter floating-point exception flag. More... | |
class | LLPerBinFracSystErr |
class | Var2DFunc |
Helper for Var2D. More... | |
class | Var3DFunc |
Helper for Var3D. More... | |
class | _Var |
Most useful for combining weights. More... | |
class | WildcardSource |
File source based on a wildcard (glob) More... | |
class | CCFlavSel |
Helper for defining true CC event cuts. More... | |
class | NCFlavOrig |
Helper for defining true CC event cuts. More... | |
class | CountingExperiment |
Compare a data spectrum to MC expectation using only the event count. More... | |
class | GaussianConstraint |
A simple Gaussian constraint on an arbitrary IFitVar. More... | |
class | IExperiment |
Base class defining interface for experiments. More... | |
class | MultiExperiment |
Combine multiple component experiments. More... | |
class | MultiExperimentSBN |
Combine multiple component experiments. More... | |
class | RatioExperiment |
class | ReactorExperiment |
Very simple model allowing inclusion of reactor constraints. More... | |
class | SingleSampleExperiment |
Compare a single data spectrum to the MC + cosmics expectation. More... | |
class | SolarConstraints |
Constraints on the parameters and from solar experiments. More... | |
class | IExtrap |
Interface to extrapolation procedures. More... | |
class | TrivialExtrap |
"Extrapolation" that simply returns the MC prediction More... | |
class | IncrementalCholeskyDecomp |
class | IPrediction |
Standard interface to all prediction techniques. More... | |
class | PredictionExtrap |
Take the output of an extrapolation and oscillate it as required. More... | |
class | IPredictionGenerator |
Given loaders and an MC shift, Generate() generates an IPrediction. More... | |
class | PredictionIncDirt |
Prediction summing detector and dirt components. More... | |
class | PredictionInterp |
Implements systematic errors by interpolation between shifted templates. More... | |
class | PredictionLinFit |
Parameterize a collection of universes as a function of the syst knobs. More... | |
class | PredictionNoExtrap |
Prediction that just uses one detector's MC, with no extrapolation. More... | |
class | NoExtrapPredictionGenerator |
class | PredictionNoOsc |
Prediction that wraps a simple Spectrum. More... | |
class | NoOscPredictionGenerator |
class | PredictionSBNExtrap |
TODO comment. More... | |
class | SBNExtrapGenerator |
class | PredictionScaleComp |
Prediction broken down into arbitrary components whose scales can be varied independently. More... | |
class | BoosterFluxHadronSyst |
struct | BoosterFluxHadronSystVector |
class | EnergyScaleSyst |
class | NuMIFluxSyst |
class | UniverseWeight |
class | SBNWeightSyst |
class | SystComponentScale |
Uncertainty in the scale of a single component of the spectrum. More... | |
class | UniverseOracle |
class | FitTheta13 |
More... | |
class | FitSinSq2Theta13 |
More... | |
class | FitDeltaInPiUnits |
More... | |
class | FitTheta23 |
More... | |
class | FitSinSqTheta23 |
More... | |
class | FitSinSq2Theta23 |
More... | |
class | FitDmSq32 |
More... | |
class | FitDmSq32Scaled |
More... | |
class | FitTanSqTheta12 |
More... | |
class | FitSinSq2Theta12 |
More... | |
class | FitDmSq21 |
More... | |
class | FitRho |
More... | |
class | FitDmSq32Sterile |
More... | |
class | FitDmSq41Sterile |
More... | |
class | FitDmSq43Sterile |
More... | |
class | FitDelta13InPiUnitsSterile |
More... | |
class | FitDelta14InPiUnitsSterile |
More... | |
class | FitDelta24InPiUnitsSterile |
More... | |
class | FitTheta13Sterile |
More... | |
class | FitSinSqTheta13Sterile |
More... | |
class | FitTheta23Sterile |
More... | |
class | FitSinSqTheta23Sterile |
More... | |
class | FitTheta14Sterile |
More... | |
class | FitSinSqTheta14Sterile |
More... | |
class | FitSinSq2Theta14Sterile |
More... | |
class | FitTheta24Sterile |
More... | |
class | FitSinSqTheta24Sterile |
More... | |
class | FitSinSq2Theta24Sterile |
More... | |
class | FitTheta34Sterile |
More... | |
class | FitSinSqTheta34Sterile |
More... | |
class | FitSinSq2Theta34Sterile |
More... | |
class | FitTheta13InDegreesSterile |
More... | |
class | FitTheta23InDegreesSterile |
More... | |
class | FitTheta14InDegreesSterile |
More... | |
class | FitTheta24InDegreesSterile |
More... | |
class | FitTheta34InDegreesSterile |
More... | |
class | FitDmSqSterile |
More... | |
class | FluxTimesNuclei |
class | ApaCrossCosmicIdAlg |
class | CosmicIdAlg |
class | CpaCrossCosmicIdAlg |
class | CrtHitCosmicIdAlg |
class | CrtTrackCosmicIdAlg |
class | FiducialVolumeCosmicIdAlg |
class | GeometryCosmicIdAlg |
class | PandoraT0CosmicIdAlg |
class | StoppingParticleCosmicIdAlg |
class | PFPSliceValidation |
Typedefs | |
typedef double( | ExposureFunc_t )(const caf::SRSpill *spill) |
typedef _Cut< caf::SRSliceProxy > | SliceCut |
Representation of a cut (selection) to be applied to a caf::StandardRecord object. More... | |
typedef _Cut< caf::SRSliceProxy > | Cut |
typedef _Cut< caf::SRSpillProxy > | SpillCut |
Equivalent of Cut acting on caf::SRSpill. For use in spill-by-spill data quality cuts. More... | |
typedef _Cut < caf::SRSpillTruthBranch > | SpillTruthCut |
Cut designed to be used over the nuTree, ie all neutrinos, not just those that got slices. More... | |
typedef _HistAxis< Var > | HistAxis |
typedef _HistAxis< SpillVar > | SpillHistAxis |
typedef _MultiVar < caf::SRSliceProxy > | MultiVar |
typedef _MultiVar < caf::SRSpillProxy > | SpillMultiVar |
typedef _Var< caf::SRSliceProxy > | Var |
Representation of a variable to be retrieved from a caf::StandardRecord object. More... | |
typedef _Var< caf::SRSpillProxy > | SpillVar |
Equivalent of Var acting on caf::SRSpill. More... | |
Functions | |
void | ResetOscCalcToDefault (osc::IOscCalcAdjustable *calc) |
Reset calc to default assumptions for all parameters. More... | |
osc::IOscCalcAdjustable * | DefaultOscCalc () |
Create a new calc with default assumptions for all parameters. More... | |
void | ResetOscCalcToDefaultIH (osc::IOscCalcAdjustable *calc) |
osc::IOscCalcAdjustable * | DefaultOscCalcIH () |
void | ResetSterileCalcToDefault (osc::OscCalcSterile *calc) |
Reset calc to default assumptions for all parameters. More... | |
osc::OscCalcSterile * | DefaultSterileCalc (int nflavors) |
Create a sterile calc with default assumptions for all parameters. More... | |
osc::IOscCalcAdjustable * | NuFitOscCalc (int hie) |
bool | HasVar (std::vector< const IFitVar * > oscVars, std::string name) |
osc::IOscCalcAdjustable * | ThrownNuFitOscCalc (int hie, std::vector< const IFitVar * > oscVars) |
osc::IOscCalcAdjustable * | NuFitOscCalcPlusOneSigma (int hie) |
osc::IOscCalcAdjustable * | NuFitOscCalcCDR (int hie) |
osc::IOscCalcAdjustable * | ThrownNuFitOscCalcCDR (int hie) |
osc::IOscCalcAdjustable * | NuFitOscCalcCDRPlusOneSigma (int hie) |
void | MakeNMinusOneSpectra (SpectrumLoader &loader, const Cut &sigcut, const Cut &presel, const std::vector< HistAxis > &axes, const std::vector< double > &cut_pos, std::vector< Spectrum > &sigs, std::vector< Spectrum > &bkgs) |
Make a series of spectra leaving out one cut in turn. More... | |
double | FindOptimumCut (TH1 *hsig, TH1 *hbkg, double &best_fom) |
Search for optimum cut position given signal and background histograms. More... | |
double | OptimizeOneCut (const std::string &wildcard, double pot, const Cut &sigcut, const Cut &presel, const std::vector< HistAxis > &axes, std::vector< double > &cut_pos) |
Scan all cuts and update the one giving the largest FOM gain. More... | |
void | OptimizeCuts (const std::string &wildcard, double pot, const Cut &sigcut, const Cut &presel, const std::vector< HistAxis > &axes, std::vector< double > &cut_pos) |
Repeatedly invoke OptimizeOneCut until the FOM increase becomes small. More... | |
double | SimpleFOM (const Spectrum &obs, const Spectrum &unosc, double pot=0) |
Figure-of-merit with no systematics, for binned data. More... | |
TH1 * | Profile (const IExperiment *expt, osc::IOscCalcAdjustable *calc, const IFitVar *v, int nbinsx, double minx, double maxx, double minchi=-1, const std::vector< const IFitVar * > &profVars={}, const std::vector< const ISyst * > &profSysts={}, const std::map< const IFitVar *, std::vector< double >> &seedPts={}, const std::vector< SystShifts > &systsSeedPts={}, std::map< const IFitVar *, TGraph * > &profVarsMap=empty_vars_map, std::map< const ISyst *, TGraph * > &systsMap=empty_syst_map) |
scan in one variable, profiling over all others More... | |
TH1 * | Profile (const IExperiment *expt, osc::IOscCalcAdjustable *calc, const ISyst *s, int nbinsx, double minx, double maxx, double input_minchi, const std::vector< const IFitVar * > &profVars, const std::vector< const ISyst * > &profSysts, const std::map< const IFitVar *, std::vector< double >> &seedPts, const std::vector< SystShifts > &systSeedPts, std::map< const IFitVar *, TGraph * > &profVarsMap, std::map< const ISyst *, TGraph * > &profSystsMap) |
TH1 * | SqrtProfile (const IExperiment *expt, osc::IOscCalcAdjustable *calc, const IFitVar *v, int nbinsx, double minx, double maxx, double minchi=-1, std::vector< const IFitVar * > profVars={}, std::vector< const ISyst * > profSysts={}, const std::map< const IFitVar *, std::vector< double >> &seedPts={}, const std::vector< SystShifts > &systsSeedPts={}, std::map< const IFitVar *, TGraph * > &profVarsMap=empty_vars_map, std::map< const ISyst *, TGraph * > &systsMap=empty_syst_map) |
Forward to Profile but sqrt the result for a crude significance. More... | |
TH1 * | SqrtProfile (const IExperiment *expt, osc::IOscCalcAdjustable *calc, const ISyst *s, int nbinsx, double minx, double maxx, double minchi, std::vector< const IFitVar * > profVars, std::vector< const ISyst * > profSysts, const std::map< const IFitVar *, std::vector< double >> &seedPts, const std::vector< SystShifts > &systSeedPts, std::map< const IFitVar *, TGraph * > &profVarsMap, std::map< const ISyst *, TGraph * > &systsMap) |
TH1 * | Slice (const IExperiment *expt, osc::IOscCalcAdjustable *calc, const IFitVar *v, int nbinsx, double minx, double maxx, double minchi=-1) |
scan in one variable, holding all others constant More... | |
TH1 * | Slice (const IExperiment *expt, osc::IOscCalcAdjustable *calc, const ISyst *s, int nbinsx, double minx, double maxx, double minchi) |
TH1 * | SqrtSlice (const IExperiment *expt, osc::IOscCalcAdjustable *calc, const IFitVar *v, int nbinsx, double minx, double maxx, double minchi=-1) |
Forward to Slice but sqrt the result for a crude significance. More... | |
TH1 * | SqrtSlice (const IExperiment *expt, osc::IOscCalcAdjustable *calc, const ISyst *s, int nbinsx, double minx, double maxx, double minchi) |
TGraph * | FindValley (const IExperiment *expt, osc::IOscCalcAdjustable *calc, const IFitVar &scanVar, const IFitVar &fitVar, int nbinsx, double xmin, double xmax, const std::vector< const IFitVar * > &profVars={}, const std::vector< const ISyst * > &profSysts={}, const std::map< const IFitVar *, std::vector< double >> &seedPts={}, const std::vector< SystShifts > &systsSeedPts={}, bool transpose=false) |
Find the minimum in one variable as a function of another. More... | |
std::vector< double > | FindCurveCrossings (TH1 *h, double critVal) |
Intended for use on the output of Profile. More... | |
Fitter::Precision | operator| (Fitter::Precision a, Fitter::Precision b) |
TH1 * | DataMCComparison (const Spectrum &data, const Spectrum &mc) |
TH1 * | DataMCComparisonAreaNormalized (const Spectrum &data, const Spectrum &mc) |
TH1 * | DataMCComparisonComponents (const Spectrum &data, const IPrediction *mc, osc::IOscCalc *calc) |
Plot MC broken down into flavour components, overlayed with data. More... | |
TH1 * | GetMCSystTotal (const IPrediction *mc, osc::IOscCalc *calc, const SystShifts &shift, std::string hist_name, double pot, bool force1D) |
TH1 * | GetMCTotal (const IPrediction *mc, osc::IOscCalc *calc, std::string hist_name, double pot, bool force1D) |
std::vector< TH1 * > | GetMCComponents (const IPrediction *mc, osc::IOscCalc *calc, std::string hist_name, double pot, bool force1D) |
std::vector< TH1 * > | GetMCSystComponents (const IPrediction *mc, osc::IOscCalc *calc, const SystShifts &shift, std::string hist_name, double pot, bool force1D) |
std::vector< TH1 * > | GetMCTotalForSystShifts (const IPrediction *mc, osc::IOscCalc *calc, const ISyst *syst, std::string hist_base_name, double pot, bool force1D) |
void | DataMCRatio (const Spectrum &data, const IPrediction *mc, osc::IOscCalc *calc, double miny=0, double maxy=3) |
Plot data/MC ratio for the given spectrum. Normalize MC to Data by POT. More... | |
void | DataMCRatio (const Spectrum &data, const Spectrum &mc, double miny=0, double maxy=3) |
Plot data/MC ratio for the given spectrum. Normalize MC to Data by POT. More... | |
void | DataMCAreaNormalizedRatio (const Spectrum &data, const IPrediction *mc, osc::IOscCalc *calc, double miny=0, double maxy=3) |
Plot data/MC ratio for the given spectrum. Normalize MC to Data by area. More... | |
void | DataMCAreaNormalizedRatio (const Spectrum &data, const Spectrum &mc, double miny=0, double maxy=3) |
Plot data/MC ratio for the given spectrum. Normalize MC to Data by area. More... | |
void | RatioPlot (const Spectrum &data, const Spectrum &expected, const Spectrum &fit, double miny=0, double maxy=1.2) |
Plot data/expected, compared with fit/expected. More... | |
void | PlotWithSystErrorBand (IPrediction *pred, const std::vector< const ISyst * > &systs, osc::IOscCalc *calc, double pot, int col=-1, int errCol=-1, float headroom=1.3, bool newaxis=true) |
Plot prediction with +/-1sigma error band. More... | |
void | PlotWithSystErrorBand (const Spectrum &nominal, const std::vector< Spectrum > &upShifts, const std::vector< Spectrum > &downShifts, double pot, int col=-1, int errCol=-1, float headroom=1.3, bool newaxis=true) |
Plot prediction with error band. More... | |
THStack * | ToTHStack (const std::vector< std::pair< Spectrum, Color_t >> &s, double pot) |
Can call like ToTHStack({{h1, kRed}, {h2, kBlue}}, pot) More... | |
double | PointDistanceToBox (double x, double y, double x0, double y0, double x1, double y1) |
Helper for AutoPlaceLegend. More... | |
TLegend * | AutoPlaceLegend (double dx, double dy, double yPin=-1) |
Create a legend, maximizing distance from all histograms. More... | |
TGraphAsymmErrors * | GraphWithPoissonErrors (const TH1 *h, bool noErrorsXaxis=false, bool drawEmptyBins=true) |
Calculate statistical errors appropriate for small Poisson numbers. More... | |
TGraph * | ShadeBetweenHistograms (TH1 *hmin, TH1 *hmax) |
TGraphAsymmErrors * | ProfileQuantile (const TH2 *hist, const std::string &axis_name, const std::string &graph_name="", const std::pair< double, double > &quantile_divisions={0.159, 0.841}) |
Calculate profile with error bars corresponding to specified quantiles of a 2D distribution (by default, 68% coverage) More... | |
void | CountingExperimentErrorBarChart (const std::map< std::string, double > &systs, double statErr=0, bool bkgdOrSig=false, bool shortchart=false) |
Make a simple plot of relative size of different errors. More... | |
double | BinCenter (const TAxis *ax, int bin, bool islog) |
TH2 * | Flat (double level, const Surface &s) |
Helper function for the gaussian approximation surfaces. More... | |
TH2 * | Gaussian68Percent2D (const Surface &s) |
Up-value surface for 68% confidence in 2D in gaussian approximation. More... | |
TH2 * | Gaussian90Percent2D (const Surface &s) |
Up-value surface for 90% confidence in 2D in gaussian approximation. More... | |
TH2 * | Gaussian95Percent2D (const Surface &s) |
Up-value surface for 95% confidence in 2D in gaussian approximation. More... | |
TH2 * | Gaussian2Sigma2D (const Surface &s) |
Up-value surface for 2 sigma confidence in 2D in gaussian approximation. More... | |
TH2 * | Gaussian99Percent2D (const Surface &s) |
Up-value surface for 99% confidence in 2D in gaussian approximation. More... | |
TH2 * | Gaussian3Sigma2D (const Surface &s) |
Up-value surface for 3 sigma confidence in 2D in gaussian approximation. More... | |
TH2 * | Gaussian5Sigma2D (const Surface &s) |
Up-value surface for 5 sigma confidence in 2D in gaussian approximation. More... | |
TH2 * | Gaussian68Percent1D (const Surface &s) |
Up-value surface for 68% confidence in 1D in gaussian approximation. More... | |
TH2 * | Gaussian90Percent1D (const Surface &s) |
Up-value surface for 90% confidence in 1D in gaussian approximation. More... | |
TH2 * | Gaussian95Percent1D (const Surface &s) |
Up-value surface for 95% confidence in 1D in gaussian approximation. More... | |
TH2 * | Gaussian2Sigma1D (const Surface &s) |
Up-value surface for 2 sigma confidence in 1D in gaussian approximation. More... | |
TH2 * | Gaussian99Percent1D (const Surface &s) |
Up-value surface for 99% confidence in 1D in gaussian approximation. More... | |
TH2 * | Gaussian3Sigma1D (const Surface &s) |
Up-value surface for 3 sigma confidence in 1D in gaussian approximation. More... | |
TH2 * | Gaussian90Percent1D1Sided (const Surface &s) |
Up-value surface for 90% confidence in 1D in 1-sided gaussian approximation. More... | |
TH2 * | Gaussian95Percent1D1Sided (const Surface &s) |
Up-value surface for 95% confidence in 1D in 1-sided gaussian approximation. More... | |
TH2 * | Gaussian99Percent1D1Sided (const Surface &s) |
Up-value surface for 99% confidence in 1D in 1-sided gaussian approxiamtion. More... | |
TH2 * | Gaussian3Sigma1D1Sided (const Surface &s) |
Up-value surface for 3 sigma confidence in 1D in 1-sided gaussian approximation. More... | |
TH2 * | Gaussian5Sigma1D1Sided (const Surface &s) |
Up-value surface for 5 sigma confidence in 1D in 1-sided gaussian approximation. More... | |
Binning | TrueEnergyBins () |
Default true-energy bin edges. More... | |
Binning | TrueLOverTrueEBins () |
std::function< ExposureFunc_t > | CombineExposures (const std::function< ExposureFunc_t > &a, const std::function< ExposureFunc_t > &b) |
template<class T > | |
_Cut< T > | operator&& (const _Cut< T > &a, const _Cut< T > &b) |
template Cut | operator&&< caf::SRSliceProxy > (const Cut &a, const Cut &b) |
template SpillCut | operator&&< caf::SRSpillProxy > (const SpillCut &a, const SpillCut &b) |
template<class T > | |
_Cut< T > | operator|| (const _Cut< T > &a, const _Cut< T > &b) |
template Cut | operator||< caf::SRSliceProxy > (const Cut &a, const Cut &b) |
template SpillCut | operator||< caf::SRSpillProxy > (const SpillCut &a, const SpillCut &b) |
template<class T > | |
_Cut< T > | operator! (const _Cut< T > &a) |
template Cut | operator!< caf::SRSliceProxy > (const Cut &a) |
template SpillCut | operator!< caf::SRSpillProxy > (const SpillCut &a) |
template<class T > | |
_Cut< T > | operator> (const _Var< T > &v, double c) |
template<class T > | |
_Cut< T > | operator>= (const _Var< T > &v, double c) |
template<class T > | |
_Cut< T > | operator< (const _Var< T > &v, double c) |
template<class T > | |
_Cut< T > | operator<= (const _Var< T > &v, double c) |
template<class T > | |
_Cut< T > | operator== (const _Var< T > &v, double c) |
template<class T > | |
_Cut< T > | operator!= (const _Var< T > &v, double c) |
template<class T > | |
_Cut< T > | operator> (const _Var< T > &a, const _Var< T > &b) |
template<class T > | |
_Cut< T > | operator>= (const _Var< T > &a, const _Var< T > &b) |
template<class T > | |
_Cut< T > | operator== (const _Var< T > &a, const _Var< T > &b) |
template<class T > | |
_Cut< T > | operator> (double c, const _Var< T > &v) |
template<class T > | |
_Cut< T > | operator< (double c, const _Var< T > &v) |
template<class T > | |
_Cut< T > | operator>= (double c, const _Var< T > &v) |
template<class T > | |
_Cut< T > | operator<= (double c, const _Var< T > &v) |
template<class T > | |
_Cut< T > | operator!= (double c, const _Var< T > &v) |
template<class T > | |
_Cut< T > | operator< (const _Var< T > &a, const _Var< T > &b) |
template<class T > | |
_Cut< T > | operator<= (const _Var< T > &a, const _Var< T > &b) |
template<class T > | |
_Cut< T > | operator!= (const _Var< T > &a, const _Var< T > &b) |
template Cut | operator> (const Var &, double) |
template Cut | operator< (const Var &, double) |
template Cut | operator>= (const Var &, double) |
template Cut | operator<= (const Var &, double) |
template Cut | operator== (const Var &, double) |
template Cut | operator!= (const Var &, double) |
template Cut | operator> (const Var &, const Var &) |
template Cut | operator< (const Var &, const Var &) |
template Cut | operator>= (const Var &, const Var &) |
template Cut | operator<= (const Var &, const Var &) |
template Cut | operator== (const Var &, const Var &) |
template Cut | operator!= (const Var &, const Var &) |
template Cut | operator> (double, const Var &) |
template Cut | operator< (double, const Var &) |
template Cut | operator>= (double, const Var &) |
template Cut | operator<= (double, const Var &) |
const Cut | kNoCut ([](const caf::SRSliceProxy *){return true;}) |
The simplest possible cut: pass everything, used as a default. More... | |
const SpillCut | kNoSpillCut ([](const caf::SRSpillProxy *){return true;}) |
The simplest possible cut: pass everything, used as a default. More... | |
void | gdb_backtrace () |
void | handle_terminate () |
void | handle_signal (int sig, siginfo_t *, void *) |
EnsembleRatio | operator/ (const EnsembleSpectrum &lhs, const EnsembleSpectrum &rhs) |
void | DrawErrorBand (TH1 *nom, TGraphAsymmErrors *band, int bandCol, double alpha) |
EnsembleSpectrum | operator* (const EnsembleRatio &lhs, const EnsembleSpectrum &rhs) |
EnsembleSpectrum | operator/ (const EnsembleRatio &lhs, const EnsembleSpectrum &rhs) |
void | ClearTrueParticles (caf::StandardRecord *sr) |
std::vector< std::string > | GetGenieWeightNames () |
std::pair< double, double > | GetGenieDialLimits (const std::string &name) |
int | GetGenieIndex (const std::string &name, bool quiet) |
std::string | GetGenieWeightName (int index) |
double | GetGenieMin (int index) |
double | GetGenieMax (int index) |
fBins ({binsX, binsY}) | |
fVars ({varX, varY}) | |
template<> | |
std::unique_ptr< osc::IOscCalc > | LoadFrom< osc::IOscCalc > (TDirectory *dir) |
template<> | |
void | SaveTo (const osc::IOscCalc &x, TDirectory *dir) |
template<class T > | |
std::unique_ptr< T > | LoadFrom (TDirectory *dir) |
template<> | |
std::unique_ptr< IDecomp > | LoadFrom< IDecomp > (TDirectory *dir) |
template<> | |
std::unique_ptr< IExtrap > | LoadFrom< IExtrap > (TDirectory *dir) |
template<> | |
std::unique_ptr< IPrediction > | LoadFrom< IPrediction > (TDirectory *dir) |
template<> | |
std::unique_ptr< IExperiment > | LoadFrom< IExperiment > (TDirectory *dir) |
template<> | |
std::unique_ptr < ModularExtrapComponent > | LoadFrom< ModularExtrapComponent > (TDirectory *dir) |
template<> | |
std::unique_ptr< IBkgdEstimator > | LoadFrom< IBkgdEstimator > (TDirectory *dir) |
template<class T > | |
void | SaveTo (const T &x, TDirectory *dir) |
template<class T > | |
std::unique_ptr< T > | LoadFromFile (const std::string &fname, const std::string &label) |
template<class T > | |
void | SaveToFile (const T &x, const std::string &fname, const std::string &label) |
template<class T > | |
_MultiVar< T > | MultiVar2D (const _MultiVar< T > &a, const Binning &binsa, const _MultiVar< T > &b, const Binning &binsb) |
template MultiVar | MultiVar2D (const MultiVar &, const Binning &, const MultiVar &, const Binning &) |
double | Si (double x) |
double | AvgSinSq (double k, double a, double b) |
OscCalcSterileApproxAdjustable * | DefaultSterileApproxCalc (SterileOscAngles angles) |
const OscCalcSterileApprox * | DowncastToSterileApprox (const osc::IOscCalc *calc, bool allowFail) |
OscCalcSterileApprox * | DowncastToSterileApprox (osc::IOscCalc *calc, bool allowFail) |
SterileOscAngles | operator| (const SterileOscAngles a, const SterileOscAngles b) |
SterileOscAngles | operator& (const SterileOscAngles a, const SterileOscAngles b) |
const Var | kBaseline ([](const caf::SRSliceProxy *slc) -> double{return slc->truth.baseline *1e-3;}) |
fCachedHash (0) | |
Ratio | operator/ (const Spectrum &lhs, const Spectrum &rhs) |
void | ProjectionX (TH2D *from, TH1D *to) |
Helper for Unweighted. More... | |
void | ProjectionY (TH2D *from, TH1D *to) |
Helper for WeightingVariable. More... | |
Spectrum | operator* (const Ratio &lhs, const Spectrum &rhs) |
Spectrum | operator/ (const Ratio &lhs, const Spectrum &rhs) |
std::string | UniqueName () |
Return a different string each time, for creating histograms. More... | |
std::string | Experiment () |
$EXPERIMENT or a nice error message and abort More... | |
std::string | SAMExperiment () |
$SAM_EXPERIMENT or a nice error message and abort More... | |
std::unique_ptr< TMatrixD > | CalcCovMx (const std::vector< TArrayD * > &binSets, int firstBin=0, int lastBin=-1) |
Compute bin-to-bin covariance matrix from a collection of sets of bin contents. More... | |
double | LogLikelihood (double exp, double obs) |
The log-likelihood formula for a single bin. More... | |
double | LogLikelihood (const TH1 *exp, const TH1 *obs, bool useOverflow=false) |
The log-likelihood formula from the PDG. More... | |
double | Chi2CovMx (const TVectorD *e, const TVectorD *o, const TMatrixD *covmxinv) |
double | Chi2CovMx (const TH1 *e, const TH1 *o, const TMatrixD *covmxinv) |
TH2F * | ExpandedHistogram (const std::string &title, int nbinsx, double xmin, double xmax, bool xlog, int nbinsy, double ymin, double ymax, bool ylog) |
std::unique_ptr< TMatrixD > | SymmMxInverse (const TMatrixD &mx) |
TH1D * | MakeTH1D (const char *name, const char *title, const Binning &bins) |
TH2D * | MakeTH2D (const char *name, const char *title, const Binning &binsx, const Binning &binsy) |
TH2 * | ToTH2 (const Spectrum &s, double exposure, ana::EExposureType expotype, const Binning &binsx, const Binning &binsy, ana::EBinType bintype) |
TH2 * | ToTH2 (const Ratio &r, const Binning &binsx, const Binning &binsy) |
TH2 * | ToTH2Helper (std::unique_ptr< TH1 > h1, const Binning &binsx, const Binning &binsy, ana::EBinType bintype) |
TH3 * | ToTH3 (const Spectrum &s, double exposure, ana::EExposureType expotype, const Binning &binsx, const Binning &binsy, const Binning &binsz, ana::EBinType bintype) |
TH3 * | ToTH3 (const Ratio &r, const Binning &binsx, const Binning &binsy, const Binning &binsz) |
TH3 * | ToTH3Helper (std::unique_ptr< TH1 > h1, const Binning &binsx, const Binning &binsy, const Binning &binsz, ana::EBinType bintype) |
std::vector< std::string > | Wildcard (const std::string &wildcardString) |
std::string | FindCAFAnaDir () |
std::vector< std::string > | LoadFileList (const std::string &listfile) |
std::map< std::string, std::string > | GetCAFMetadata (TDirectory *dir) |
void | CombineMetadata (std::map< std::string, std::string > &base, const std::map< std::string, std::string > &add, std::set< std::string > &mask) |
void | WriteCAFMetadata (TDirectory *dir, const std::map< std::string, std::string > &meta) |
bool | RunningOnGrid () |
size_t | Stride (bool allow_default) |
size_t | Offset (bool allow_default) |
int | Limit () |
size_t | JobNumber () |
size_t | NumJobs () |
bool | AlmostEqual (double a, double b) |
std::string | pnfs2xrootd (std::string loc, bool unauth) |
void | EnsurePositiveDefinite (TH2 *mat) |
TH1 * | GetMaskHist (const Spectrum &s, double xmin, double xmax, double ymin, double ymax) |
double | FindQuantile (double frac, std::vector< double > &xs) |
template<class T > | |
_Var< T > | Var2D (const _Var< T > &a, const Binning &binsa, const _Var< T > &b, const Binning &binsb) |
Variable formed from two input variables. More... | |
template<class T > | |
_Var< T > | Var2D (const _Var< T > &a, int na, double a0, double a1, const _Var< T > &b, int nb, double b0, double b1) |
Variable formed from two input variables. More... | |
template Var | Var2D (const Var &, const Binning &, const Var &, const Binning &) |
template Var | Var2D (const Var &, int, double, double, const Var &, int, double, double) |
template<class T > | |
_Var< T > | Var3D (const _Var< T > &a, const Binning &binsa, const _Var< T > &b, const Binning &binsb, const _Var< T > &c, const Binning &binsc) |
This is just like a Var2D, but useful for 3D Spectra. More... | |
template<class T > | |
_Var< T > | Var3D (const _Var< T > &a, int na, double a0, double a1, const _Var< T > &b, int nb, double b0, double b1, const _Var< T > &c, int nc, double c0, double c1) |
This is just like a Var2D, but useful for 3D Spectra. More... | |
template Var | Var3D (const Var &, const Binning &, const Var &, const Binning &, const Var &, const Binning &) |
template Var | Var3D (const Var &, int, double, double, const Var &, int, double, double, const Var &, int, double, double) |
Var | Scaled (const Var &v, double s) |
Use to rescale another variable. More... | |
Var | Constant (double c) |
Use to weight events up and down by some factor. More... | |
Var | Sqrt (const Var &v) |
Use to take sqrt of a var. More... | |
template<class T > | |
_Var< T > | operator* (const _Var< T > &a, const _Var< T > &b) |
template<class T > | |
_Var< T > | operator/ (const _Var< T > &a, const _Var< T > &b) |
template<class T > | |
_Var< T > | operator+ (const _Var< T > &a, const _Var< T > &b) |
template<class T > | |
_Var< T > | operator- (const _Var< T > &a, const _Var< T > &b) |
template Var | operator* (const Var &, const Var &) |
template Var | operator/ (const Var &, const Var &) |
template Var | operator+ (const Var &, const Var &) |
template Var | operator- (const Var &, const Var &) |
const Var | kUnweighted ([](const caf::SRSliceProxy *){return 1;}) |
The simplest possible Var, always 1. Used as a default weight. More... | |
const SpillVar | kSpillUnweighted ([](const caf::SRSpillProxy *){return 1;}) |
const Cut | kIsNueApp (CCFlavSel(12, 14)) |
Select CC . More... | |
const Cut | kIsNumuCC (CCFlavSel(14, 14)) |
Select CC . More... | |
const Cut | kIsBeamNue (CCFlavSel(12, 12)) |
Select CC . More... | |
const Cut | kIsNumuApp (CCFlavSel(14, 12)) |
Select CC . More... | |
const Cut | kIsTauFromMu (CCFlavSel(16, 14)) |
Select CC . More... | |
const Cut | kIsTauFromE (CCFlavSel(16, 12)) |
Select CC . More... | |
const Cut | kIsNCFromNumu (NCFlavOrig(14)) |
const Cut | kIsNCFromNue (NCFlavOrig(12)) |
const ReactorExperiment * | DayaBayConstraint2014 () |
A ReactorExperiment initialized with the Nu2014 Daya Bay constraints. More... | |
const ReactorExperiment * | WorldReactorConstraint2015 () |
Weighted average of all experiments as of first nue paper writing. More... | |
const ReactorExperiment * | WorldReactorConstraint2016 () |
Updated value for SecondAna based on the latest PDG. More... | |
const ReactorExperiment * | WorldReactorConstraint2017 () |
Reactor constraint from PDG 2017 update. More... | |
const SolarConstraints | kSolarConstraintsPDG2017 (7.53e-5, 0.18e-5, 0.851, 0.020) |
fSplitBySign (mode==kSplitBySign) | |
const BoosterFluxHadronSyst * | GetBoosterFluxHadronSyst (unsigned int i) |
BoosterFluxHadronSystVector | GetBoosterFluxHadronSysts (unsigned int N) |
std::vector< const ISyst * > | GetEnergySysts () |
std::vector< const ISyst * > | GetBigEnergySysts () |
const NuMIFluxSyst * | GetNuMIFluxSyst (const std::string &dir, const std::string &prefix, const std::string &name) |
std::vector< const ISyst * > | GetNuMIHadronProductionFluxSysts () |
These are envelopes not real systs. TODO make clearer in naming. More... | |
std::vector< const ISyst * > | GetNuMIBeamlineFluxSysts () |
std::vector< const ISyst * > | GetNuMIPCAFluxSysts (unsigned int Npcs) |
std::vector< const ISyst * > | GetAllNuMIFluxSysts (unsigned int Npcs) |
Combination of all beamline systs plus Npcs hadron production components. More... | |
std::vector< std::string > | GetSBNGenieWeightNames () |
const std::vector< const ISyst * > & | GetSBNGenieWeightSysts () |
std::vector< std::string > | GetSBNBoosterFluxWeightNames () |
const std::vector< const ISyst * > & | GetSBNBoosterFluxWeightSysts () |
const std::vector< const ISyst * > & | GetSBNWeightSysts () |
Var | GetUniverseWeight (const std::string &psetName, int univIdx) |
caf::SRGlobal | GetSRGlobal () |
void | PrintSRGlobal (const caf::SRGlobal &global) |
Spectrum | UnfoldIterative (const Spectrum &reco, const ReweightableSpectrum &recoVsTrue, unsigned int nIterations) |
Spectrum | UnfoldSVD (const Spectrum &reco, const ReweightableSpectrum &recoVsTrue, unsigned int reg) |
Singular value decomposition unfolding using ROOT's TSVDDecomp. More... | |
Spectrum | UnfoldTikhonov (const Spectrum &reco, const ReweightableSpectrum &recoVsTrue, double regStrength) |
Unfolding using Tikhonov regularization (penalizing true spectra with large second derivatives) More... | |
bool | IsCCQEOnArgon (const caf::SRTrueInteractionProxy *nu, int pdg) |
bool | IsNCQEOnArgon (const caf::SRTrueInteractionProxy *nu, int pdg) |
Cut | IsCCQEOnArgonCut (int pdg) |
TString | thisCellColor (double weird) |
TString | fixLatexName (TString mystring) |
void | printTableHeader (int quantId=0) |
void | printTableFooter () |
void | printEventsLine (std::string cutname, float nue, float numu, float nc, float cos, float other, float eff, float pur) |
void | SplitCanvas2 (TCanvas *&c1, TPad *&pad1, TPad *&pad2) |
float | GetHistMax (std::vector< TH1 * > histos) |
void | PimpHist (TH1 *histo, Color_t color, Style_t linestyle, int linewidth, Style_t markerstyle=8, double markersize=8) |
void | FillWithDimColor (TH1 *h, bool usealpha=false, float dim=0.8) |
void | DrawComponentsLegend (TH1 *hnue, TH1 *hnumu, TH1 *hnc, TH1 *hcos, TH1 *hother) |
void | DrawSigBkgLegend (TH1 *h1, char *name1, TH1 *h2, char *name2) |
void | DrawSigBkgIntLegend (TH1 *h1, char *name1, double iSig, TH1 *h2, char *name2, double iBkg) |
void | DrawIntEffPurLegend (TH1D *g1, char *name1, TH1D *g2, char *name2) |
void | DrawPurLegend (TH1D *g1, char *name1) |
void | DrawEffPurLegend (TH1D *g1, char *name1, TH1D *g2, char *name2) |
void | DrawSigBkgIntText (TH1 *hsig, TH1 *hbkg, float textsize) |
TH1D * | EffOrPurHistogram (TH1 *hSelSignal, TH1 *hSelBack, TH1 *hSignal, bool geteff) |
TGraph * | EffOrPurGraph (TH1 *hSelSignal, TH1 *hSelBack, TH1 *hSignal, bool geteff) |
static bool | Icarus202208contained (const caf::SRTrackProxy &trk) |
bool | IsPrimary (const caf::SRTrueParticleProxy &p) |
Whether this is a primary particle or generated by a secondary interaction. More... | |
bool | HasBraggPeak (const caf::SRTrueParticleProxy &p) |
Whether this particle should have a bragg peak in the detector. More... | |
bool | IsGenie (const caf::SRTrueParticleProxy &p) |
Whether this particle was generated by genie (as opposed to geant or corsika) More... | |
bool | IsStable (const caf::SRTrueParticleProxy &p) |
Whether this is a stable particle as generated by genie. More... | |
const Cut | kIsAntiNu ([](const caf::SRSliceProxy *slc){if(slc->truth.index< 0) return false;return slc->truth.pdg< 0;}) |
const Cut | kIsCC ([](const caf::SRSliceProxy *slc){if(slc->truth.index< 0) return false;return(slc->truth.iscc==1);}) |
const Cut | kIsNC ([](const caf::SRSliceProxy *slc){if(slc->truth.index< 0) return false;return(slc->truth.isnc==1);}) |
const caf::SRShowerProxy * | LargestRecoShower (const caf::SRSliceProxy *slc) |
Pointer to largest reconstructed shower, or null pointer if none exists. More... | |
const caf::SRTrackProxy * | LongestRecoTrack (const caf::SRSliceProxy *slc) |
Pointer to longest reconstructed shower, or null pointer if none exists. More... | |
const caf::SRTrkChi2PIDProxy * | BestPlaneChi2PID (const caf::SRTrackProxy *trk) |
const caf::SRTrkChi2PIDProxy * | LongestTrackBestPlaneChi2PID (const caf::SRSliceProxy *slc) |
static bool | Icarus202208_proton_cut (const caf::SRTrackProxy &trk) |
const Var | kNuScore ([](const caf::SRSliceProxy *slc) -> double{return slc->nu_score;}) |
const Var | kSlcIsRecoNu ([](const caf::SRSliceProxy *slc) -> double{return!slc->is_clear_cosmic;}) |
Variables | |
const double | kNuFitDmsq21CV = 7.39e-5 |
const double | kNuFitTh12CV = 33.82 * TMath::Pi()/180 |
const double | kNuFitDmsq32CVNH = +2.525e-3 - kNuFitDmsq21CV |
const double | kNuFitTh23CVNH = 49.6 * TMath::Pi()/180 |
const double | kNuFitTh13CVNH = 8.61 * TMath::Pi()/180 |
const double | kNuFitdCPCVNH = 215 * TMath::Pi()/180 |
const double | kNuFitDmsq32CVIH = -2.512e-3 |
const double | kNuFitTh23CVIH = 49.8 * TMath::Pi()/180 |
const double | kNuFitTh13CVIH = 8.65 * TMath::Pi()/180 |
const double | kNuFitdCPCVIH = 284 * TMath::Pi()/180 |
const double | kNuFitDmsq21Err = ((8.01-6.79)/6)*1e-5 |
const double | kNuFitTh12Err = ((36.27-31.61)/6) * TMath::Pi()/180 |
const double | kNuFitDmsq32ErrNH = ((2.625-2.427)/6)*1e-3 |
const double | kNuFitTh23ErrNH = ((52.4-40.3)/6) * TMath::Pi()/180 |
const double | kNuFitTh13ErrNH = ((8.99-8.22)/6) * TMath::Pi()/180 |
const double | kNuFitDmsq32ErrIH = ((2.611-2.412)/6)*1e-3 |
const double | kNuFitTh23ErrIH = ((52.5-40.6)/6) * TMath::Pi()/180 |
const double | kNuFitTh13ErrIH = ((9.03-8.27)/6) * TMath::Pi()/180 |
const double | kBaseline = 1284.9 |
const double | kEarthDensity = 2.848 |
const double | kPOTnominal = 6.6e20 |
const double | kPOTuBoone = 1.3e21 |
const double | kBaselineSBND = 0.11 |
const double | kBaselineMicroBoone = 0.47 |
const double | kBaselineIcarus = 0.6 |
const int | kSBND = 0 |
const int | kMicroBoone = 1 |
const int | kICARUS = 2 |
const std::vector< double > | kBLs = {kBaselineSBND, kBaselineMicroBoone, kBaselineIcarus} |
static std::map< const IFitVar *, TGraph * > | empty_vars_map |
static std::map< const ISyst *, TGraph * > | empty_syst_map |
const Color_t | kTotalMCColor = kRed |
const Color_t | kTotalMCErrorBandColor = kRed-10 |
const Color_t | kNueSignalColor = kViolet-5 |
const Color_t | kTotalBackgroundColor = kAzure+2 |
const Color_t | kNCBackgroundColor = kAzure |
const Color_t | kBeamNueBackgroundColor = kPink+9 |
const Color_t | kCosmicBackgroundColor = kAzure+1 |
const Color_t | kNumuBackgroundColor = kGreen+2 |
const Color_t | kNormalHierarchyColor = kAzure+3 |
const Color_t | kInvertedHierarchyColor = kOrange+9 |
const Style_t | k90PercentConfidenceStyle = 9 |
const Style_t | k68PercentConfidenceStyle = 7 |
Dashed. More... | |
const Color_t | kPrimColorNH = kAzure+2 |
const Color_t | kSecoColorNH = kAzure+1 |
const Color_t | kPrimColorIH = kOrange+9 |
const Color_t | kSecoColorIH = kOrange+10 |
const Style_t | kFillStyleNH = 3013 |
const Style_t | kFillStyleIH = 3003 |
const Style_t | k1SigmaConfidenceStyle = 7 |
const Style_t | k2SigmaConfidenceStyle = 9 |
const Style_t | k3SigmaConfidenceStyle = 10 |
const Color_t | kCentralValueColorNH = kAzure+3 |
const Color_t | k1SigmaConfidenceColorNH = kAzure+3 |
const Color_t | k90PercConfidenceColorNH = kAzure+2 |
const Color_t | k2SigmaConfidenceColorNH = kAzure+2 |
const Color_t | k3SigmaConfidenceColorNH = kAzure+1 |
const Color_t | kCentralValueColorIH = kOrange+9 |
const Color_t | k1SigmaConfidenceColorIH = kOrange+9 |
const Color_t | k2SigmaConfidenceColorIH = kOrange+10 |
const Color_t | k3SigmaConfidenceColorIH = kOrange+8 |
const Color_t | k90PercConfidenceColorIH = kOrange+10 |
const Color_t | kCentralValueColor = kGray+2 |
const Style_t | kCentralValueStyle = kSolid |
const Float_t | kBlessedTitleFontSize = 0.08 |
const Float_t | kBlessedLabelFontSize = 0.06 |
const Binning | kTrueEnergyBins = TrueEnergyBins() |
Default true-energy bin edges. More... | |
const Binning | kTrueLOverEBins = Binning::Simple(100, 0, 2) |
bool | gIsException = false |
static InstallHandlers | gHandlerInstaller |
const Var | kTrueE ([](const caf::SRSliceProxy *slc) -> double{return slc->truth.E;}) |
const Var | kTrueLOverE = kBaseline / kTrueE |
static NullLoader | kNullLoader |
Dummy loader that doesn't load any files. More... | |
const SystShifts | kNoShift = SystShifts::Nominal() |
const Cut | kHasMatchedNu ([](const caf::SRSliceProxy *slc){return slc->truth.index >=0;}) |
const Cut | kIsNC ([](const caf::SRSliceProxy *slc){return kHasMatchedNu(slc)&&slc->truth.isnc;}) |
Is this a Neutral Current event? More... | |
const Cut | kIsCC ([](const caf::SRSliceProxy *slc){return kHasMatchedNu(slc)&&slc->truth.iscc;}) |
const Cut | kIsAntiNu ([](const caf::SRSliceProxy *slc){return kHasMatchedNu(slc)&&slc->truth.pdg< 0;}) |
Is this truly an antineutrino? More... | |
const EnergyScaleSyst | kEnergyScaleMuon (EnergyScaleSystTerm::kConstant, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kAll, 0.02,"EnergyScaleMuon","Correlated linear E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonSqrt (EnergyScaleSystTerm::kSqrt, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kAll, 0.02,"EnergyScaleMuonSqrt","Correlated sqrt E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonInvSqrt (EnergyScaleSystTerm::kInverseSqrt, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kAll, 0.02,"EnergyScaleMuonInvSqrt","Correlated inv sqrt E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonND (EnergyScaleSystTerm::kConstant, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kSBND, 0.02,"EnergyScaleMuonND","Uncorrelated SBND linear E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonSqrtND (EnergyScaleSystTerm::kSqrt, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kSBND, 0.02,"EnergyScaleMuonSqrtND","Uncorrelated SBND sqrt E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonInvSqrtND (EnergyScaleSystTerm::kInverseSqrt, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kSBND, 0.02,"EnergyScaleMuonInvSqrtND","Uncorrelated SBND inv sqrt E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonUB (EnergyScaleSystTerm::kConstant, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kMicroBooNE, 0.02,"EnergyScaleMuonUB","Uncorrelated MicroBooNE linear E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonSqrtUB (EnergyScaleSystTerm::kSqrt, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kMicroBooNE, 0.02,"EnergyScaleMuonSqrtUB","Uncorrelated MicroBooNE sqrt E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonInvSqrtUB (EnergyScaleSystTerm::kInverseSqrt, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kMicroBooNE, 0.02,"EnergyScaleMuonInvSqrtUB","Uncorrelated MicroBooNE inv sqrt E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonFD (EnergyScaleSystTerm::kConstant, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kICARUS, 0.02,"EnergyScaleMuonFD","Uncorrelated ICARUS linear E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonSqrtFD (EnergyScaleSystTerm::kSqrt, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kICARUS, 0.02,"EnergyScaleMuonSqrtFD","Uncorrelated ICARUS sqrt E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonInvSqrtFD (EnergyScaleSystTerm::kInverseSqrt, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kICARUS, 0.02,"EnergyScaleMuonInvSqrtFD","Uncorrelated ICARUS inv sqrt E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleHadron (EnergyScaleSystTerm::kConstant, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kAll, 0.05,"EnergyScaleHadron","Correlated linear E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronSqrt (EnergyScaleSystTerm::kSqrt, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kAll, 0.05,"EnergyScaleHadronSqrt","Correlated sqrt E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronInvSqrt (EnergyScaleSystTerm::kInverseSqrt, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kAll, 0.05,"EnergyScaleHadronInvSqrt","Correlated inv sqrt E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronND (EnergyScaleSystTerm::kConstant, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kSBND, 0.05,"EnergyScaleHadronND","Uncorrelated SBND linear E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronSqrtND (EnergyScaleSystTerm::kSqrt, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kSBND, 0.05,"EnergyScaleHadronSqrtND","Uncorrelated SBND sqrt E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronInvSqrtND (EnergyScaleSystTerm::kInverseSqrt, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kSBND, 0.05,"EnergyScaleHadronInvSqrtND","Uncorrelated SBND inv sqrt E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronUB (EnergyScaleSystTerm::kConstant, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kMicroBooNE, 0.05,"EnergyScaleHadronUB","Uncorrelated MicroBooNE linear E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronSqrtUB (EnergyScaleSystTerm::kSqrt, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kMicroBooNE, 0.05,"EnergyScaleHadronSqrtUB","Uncorrelated MicroBooNE sqrt E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronInvSqrtUB (EnergyScaleSystTerm::kInverseSqrt, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kMicroBooNE, 0.05,"EnergyScaleHadronInvSqrtUB","Uncorrelated MicroBooNE inv sqrt E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronFD (EnergyScaleSystTerm::kConstant, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kICARUS, 0.05,"EnergyScaleHadronFD","Uncorrelated ICARUS linear E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronSqrtFD (EnergyScaleSystTerm::kSqrt, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kICARUS, 0.05,"EnergyScaleHadronSqrtFD","Uncorrelated ICARUS sqrt E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronInvSqrtFD (EnergyScaleSystTerm::kInverseSqrt, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kICARUS, 0.05,"EnergyScaleHadronInvSqrtFD","Uncorrelated ICARUS inv sqrt E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleMuonBig (EnergyScaleSystTerm::kConstant, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kAll, 0.05,"EnergyScaleMuonBig","Correlated linear E_{#mu} scale ") |
const EnergyScaleSyst | kEnergyScaleMuonSqrtBig (EnergyScaleSystTerm::kSqrt, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kAll, 0.05,"EnergyScaleMuonSqrtBig","Correlated sqrt E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonInvSqrtBig (EnergyScaleSystTerm::kInverseSqrt, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kAll, 0.05,"EnergyScaleMuonInvSqrtBig","Correlated inv sqrt E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonNDBig (EnergyScaleSystTerm::kConstant, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kSBND, 0.05,"EnergyScaleMuonNDBig","Uncorrelated SBND linear E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonSqrtNDBig (EnergyScaleSystTerm::kSqrt, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kSBND, 0.05,"EnergyScaleMuonSqrtNDBig","Uncorrelated SBND sqrt E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonInvSqrtNDBig (EnergyScaleSystTerm::kInverseSqrt, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kSBND, 0.05,"EnergyScaleMuonInvSqrtNDBig","Uncorrelated SBND inv sqrt E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonUBBig (EnergyScaleSystTerm::kConstant, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kMicroBooNE, 0.05,"EnergyScaleMuonUBBig","Uncorrelated MicroBooNE linear E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonSqrtUBBig (EnergyScaleSystTerm::kSqrt, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kMicroBooNE, 0.05,"EnergyScaleMuonSqrtUBBig","Uncorrelated MicroBooNE sqrt E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonInvSqrtUBBig (EnergyScaleSystTerm::kInverseSqrt, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kMicroBooNE, 0.05,"EnergyScaleMuonInvSqrtUBBig","Uncorrelated MicroBooNE inv sqrt E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonFDBig (EnergyScaleSystTerm::kConstant, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kICARUS, 0.05,"EnergyScaleMuonFDBig","Uncorrelated ICARUS linear E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonSqrtFDBig (EnergyScaleSystTerm::kSqrt, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kICARUS, 0.05,"EnergyScaleMuonSqrtFDBig","Uncorrelated ICARUS sqrt E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleMuonInvSqrtFDBig (EnergyScaleSystTerm::kInverseSqrt, EnergyScaleSystParticle::kMuon, EnergyScaleSystDetector::kICARUS, 0.05,"EnergyScaleMuonInvSqrtFDBig","Uncorrelated ICARUS inv sqrt E_{#mu} scale") |
const EnergyScaleSyst | kEnergyScaleHadronBig (EnergyScaleSystTerm::kConstant, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kAll, 0.10,"EnergyScaleHadronBig","Correlated linear E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronSqrtBig (EnergyScaleSystTerm::kSqrt, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kAll, 0.10,"EnergyScaleHadronSqrtBig","Correlated sqrt E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronInvSqrtBig (EnergyScaleSystTerm::kInverseSqrt, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kAll, 0.10,"EnergyScaleHadronInvSqrtBig","Correlated inv sqrt E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronNDBig (EnergyScaleSystTerm::kConstant, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kSBND, 0.10,"EnergyScaleHadronNDBig","Uncorrelated SBND linear E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronSqrtNDBig (EnergyScaleSystTerm::kSqrt, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kSBND, 0.10,"EnergyScaleHadronSqrtNDBig","Uncorrelated SBND sqrt E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronInvSqrtNDBig (EnergyScaleSystTerm::kInverseSqrt, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kSBND, 0.10,"EnergyScaleHadronInvSqrtNDBig","Uncorrelated SBND inv sqrt E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronUBBig (EnergyScaleSystTerm::kConstant, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kMicroBooNE, 0.10,"EnergyScaleHadronUBBig","Uncorrelated MicroBooNE linear E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronSqrtUBBig (EnergyScaleSystTerm::kSqrt, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kMicroBooNE, 0.10,"EnergyScaleHadronSqrtUBBig","Uncorrelated MicroBooNE sqrt E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronInvSqrtUBBig (EnergyScaleSystTerm::kInverseSqrt, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kMicroBooNE, 0.10,"EnergyScaleHadronInvSqrtUBBig","Uncorrelated MicroBooNE inv sqrt E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronFDBig (EnergyScaleSystTerm::kConstant, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kICARUS, 0.10,"EnergyScaleHadronFDBig","Uncorrelated ICARUS linear E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronSqrtFDBig (EnergyScaleSystTerm::kSqrt, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kICARUS, 0.10,"EnergyScaleHadronSqrtFDBig","Uncorrelated ICARUS sqrt E_{had} scale") |
const EnergyScaleSyst | kEnergyScaleHadronInvSqrtFDBig (EnergyScaleSystTerm::kInverseSqrt, EnergyScaleSystParticle::kHadron, EnergyScaleSystDetector::kICARUS, 0.10,"EnergyScaleHadronInvSqrtFDBig","Uncorrelated ICARUS inv sqrt E_{had} scale") |
const FitTheta13 | kFitTheta13 = FitTheta13() |
More... | |
const FitSinSq2Theta13 | kFitSinSq2Theta13 = FitSinSq2Theta13() |
More... | |
const FitDeltaInPiUnits | kFitDeltaInPiUnits = FitDeltaInPiUnits() |
More... | |
const FitTheta23 | kFitTheta23 = FitTheta23() |
More... | |
const FitSinSqTheta23 | kFitSinSqTheta23 = FitSinSqTheta23() |
More... | |
const FitSinSq2Theta23 | kFitSinSq2Theta23 = FitSinSq2Theta23() |
More... | |
const FitDmSq32 | kFitDmSq32 = FitDmSq32() |
More... | |
const FitDmSq32Scaled | kFitDmSq32Scaled = FitDmSq32Scaled() |
More... | |
const FitTanSqTheta12 | kFitTanSqTheta12 = FitTanSqTheta12() |
More... | |
const FitSinSq2Theta12 | kFitSinSq2Theta12 = FitSinSq2Theta12() |
More... | |
const FitDmSq21 | kFitDmSq21 = FitDmSq21() |
More... | |
const FitRho | kFitRho = FitRho() |
More... | |
const FitDmSq32Sterile | kFitDmSq32Sterile = FitDmSq32Sterile() |
More... | |
const FitDmSq41Sterile | kFitDmSq41Sterile = FitDmSq41Sterile() |
More... | |
const FitDmSq43Sterile | kFitDmSq43Sterile = FitDmSq43Sterile() |
More... | |
const FitDelta13InPiUnitsSterile | kFitDelta13InPiUnitsSterile = FitDelta13InPiUnitsSterile() |
More... | |
const FitDelta14InPiUnitsSterile | kFitDelta14InPiUnitsSterile = FitDelta14InPiUnitsSterile() |
More... | |
const FitDelta24InPiUnitsSterile | kFitDelta24InPiUnitsSterile = FitDelta24InPiUnitsSterile() |
More... | |
const FitTheta13Sterile | kFitTheta13Sterile = FitTheta13Sterile() |
More... | |
const FitSinSqTheta13Sterile | kFitSinSqTheta13Sterile = FitSinSqTheta13Sterile() |
More... | |
const FitTheta23Sterile | kFitTheta23Sterile = FitTheta23Sterile() |
More... | |
const FitSinSqTheta23Sterile | kFitSinSqTheta23Sterile = FitSinSqTheta23Sterile() |
More... | |
const FitTheta14Sterile | kFitTheta14Sterile = FitTheta14Sterile() |
More... | |
const FitSinSqTheta14Sterile | kFitSinSqTheta14Sterile = FitSinSqTheta14Sterile() |
More... | |
const FitSinSq2Theta14Sterile | kFitSinSq2Theta14Sterile = FitSinSq2Theta14Sterile() |
More... | |
const FitTheta24Sterile | kFitTheta24Sterile = FitTheta24Sterile() |
More... | |
const FitSinSqTheta24Sterile | kFitSinSqTheta24Sterile = FitSinSqTheta24Sterile() |
More... | |
const FitSinSq2Theta24Sterile | kFitSinSq2Theta24Sterile = FitSinSq2Theta24Sterile() |
More... | |
const FitTheta34Sterile | kFitTheta34Sterile = FitTheta34Sterile() |
More... | |
const FitSinSqTheta34Sterile | kFitSinSqTheta34Sterile = FitSinSqTheta34Sterile() |
More... | |
const FitSinSq2Theta34Sterile | kFitSinSq2Theta34Sterile = FitSinSq2Theta34Sterile() |
More... | |
const FitTheta13InDegreesSterile | kFitTheta13InDegreesSterile = FitTheta13InDegreesSterile() |
More... | |
const FitTheta23InDegreesSterile | kFitTheta23InDegreesSterile = FitTheta23InDegreesSterile() |
More... | |
const FitTheta14InDegreesSterile | kFitTheta14InDegreesSterile = FitTheta14InDegreesSterile() |
More... | |
const FitTheta24InDegreesSterile | kFitTheta24InDegreesSterile = FitTheta24InDegreesSterile() |
More... | |
const FitTheta34InDegreesSterile | kFitTheta34InDegreesSterile = FitTheta34InDegreesSterile() |
More... | |
const FitDmSqSterile | kFitDmSqSterile = FitDmSqSterile() |
m^2 More... | |
f sin | theta_ |
const Var | kInvXSec ([](const caf::SRSliceProxy *sr){const double GeV2perm2=2.56819e31;return GeV2perm2/sr->truth.xsec;}) |
const SpillCut | kFirstEvents = kEvt < 10 |
const SpillCut | kFlashTrigger ([](const caf::SRSpillProxy *sr){return(sr->pass_flashtrig);}) |
const SpillCut | kCRTHitVetoND ([](const caf::SRSpillProxy *sr){for(auto const &crtHit:sr->crt_hits){return false;}return true;}) |
const SpillCut | kCRTHitVetoFD ([](const caf::SRSpillProxy *sr){for(auto const &crtHit:sr->crt_hits){auto thistime=crtHit.time-1600.;return false;}return true;}) |
const Cut | kFiducialVolumeND ([](const caf::SRSliceProxy *slc){return PtInVolAbsX(slc->vertex, fvndAbs);}) |
const Cut | kActiveVolumeND ([](const caf::SRSliceProxy *slc){return PtInVolAbsX(slc->vertex, avnd);}) |
const Cut | kFiducialVolumeFDCryo1 ([](const caf::SRSliceProxy *slc){return PtInVol(slc->vertex, fvfd_cryo1);}) |
const Cut | kFiducialVolumeFDCryo2 ([](const caf::SRSliceProxy *slc){return PtInVol(slc->vertex, fvfd_cryo2);}) |
const Cut | kActiveVolumeFDCryo1 ([](const caf::SRSliceProxy *slc){return PtInVol(slc->vertex, avfd_cryo1);}) |
const Cut | kActiveVolumeFDCryo2 ([](const caf::SRSliceProxy *slc){return PtInVol(slc->vertex, avfd_cryo2);}) |
const Cut | kSlcIsRecoNu ([](const caf::SRSliceProxy *slc){return!slc->is_clear_cosmic;}) |
const Cut | kSlcNuScoreCut ([](const caf::SRSliceProxy *slc){return(kSlcIsRecoNu(slc)&&slc->nu_score >0.4);}) |
const Cut | kSlcHasFlashMatch ([](const caf::SRSliceProxy *slc){return slc->fmatch.present;}) |
const Cut | kSlcFlashMatchCut ([](const caf::SRSliceProxy *slc){return(kSlcHasFlashMatch(slc)&&slc->fmatch.score >0 &&slc->fmatch.score< 6);}) |
const Cut | kContainedFD = kFiducialVolumeFDCryo1 || kFiducialVolumeFDCryo2 |
const Cut | kRecoShower ([](const caf::SRSliceProxy *slc){const int largestShwIdx(kLargestRecoShowerIdx(slc));return false;return(slc->reco.shw[largestShwIdx].bestplane_energy > 0.f &&slc->reco.shw[largestShwIdx].bestplane_dEdx > 0.f &&slc->reco.shw[largestShwIdx].conversion_gap > 0.f);}) |
const Cut | kNueBasicCut ([](const caf::SRSliceProxy *slc){const int largestShwIdx(kLargestRecoShowerIdx(slc));return false;return(slc->reco.shw[largestShwIdx].bestplane_energy > 0.2f &&slc->reco.shw[largestShwIdx].bestplane_dEdx< 3 &&slc->reco.shw[largestShwIdx].conversion_gap< 3);}) |
const Cut | kShowerEnergyCut ([](const caf::SRSliceProxy *slc){const int largestShwIdx(kLargestRecoShowerIdx(slc));return false;return(slc->reco.shw[largestShwIdx].bestplane_energy > 0.2f);}) |
const Cut | kShowerdEdxCut ([](const caf::SRSliceProxy *slc){const int largestShwIdx(kLargestRecoShowerIdx(slc));return false;return(slc->reco.shw[largestShwIdx].bestplane_dEdx< 3.625f);}) |
const Cut | kShowerConvGapCut ([](const caf::SRSliceProxy *slc){const int largestShwIdx(kLargestRecoShowerIdx(slc));return false;return(slc->reco.shw[largestShwIdx].conversion_gap< 3.25f);}) |
const Cut | kNueNumShowersCut ([](const caf::SRSliceProxy *slc){return((unsigned int) kRecoShowers_EnergyCut(slc)==1);}) |
const Cut | kNueHasTrackCut ([](const caf::SRSliceProxy *slc){return slc->reco.ntrk > 0;}) |
const Cut | kNueTrackContainmentCut ([](const caf::SRSliceProxy *slc){const int longestTrackIdx(kLongestTrackIdx(slc));return true;return PtInVol(slc->reco.trk[longestTrackIdx].end, fvndExit);}) |
const Cut | kNueTrackLenCut ([](const caf::SRSliceProxy *slc){return kLongestTrackLength(slc)< 110.f;}) |
const Cut | kNueMuonCut ([](const caf::SRSliceProxy *slc){return(kLongestTrackLength(slc)< 80.f||kLongestTrackChi2Muon(slc) > 30.f||kLongestTrackChi2Proton(slc)< 60.f);}) |
const Cut | kShowerDensityCut ([](const caf::SRSliceProxy *slc){const int largestShwIdx(kLargestRecoShowerIdx(slc));return false;return(slc->reco.shw[largestShwIdx].density > 4.5);}) |
const Cut | kShowerOpenAngleCut ([](const caf::SRSliceProxy *slc){return kRecoShower_OpenAngle(slc)< 12.f;}) |
const Cut | kShowerRazzleElectronCut = kRecoShowerRazzlePID == 11 |
const Cut | kShowerRazzleElectronScoreCut = kRecoShowerRazzleElectronScore > 0.8 |
const Cut | kNueLongestTrackDazzleMuonCut = kLongestTrackDazzlePID != 13 |
const Cut | kNueLongestTrackDazzleMuonScoreCut = kLongestTrackDazzleMuonScore < 0.8 |
const Cut | kNueAllTrackDazzleMuonCut ([](const caf::SRSliceProxy *slc){for(auto const &trk:slc->reco.trk){return false;}return true;}) |
const Cut | kNueAllTrackDazzleMuonScoreCut ([](const caf::SRSliceProxy *slc){for(auto const &trk:slc->reco.trk){return false;}return true;}) |
const Cut | kNueContainedND ([](const caf::SRSliceProxy *slc){const int largestShwIdx(kLargestRecoShowerIdx(slc));return false;double this_endx=slc->reco.shw[largestShwIdx].start.x+(slc->reco.shw[largestShwIdx].dir.x *slc->reco.shw[largestShwIdx].len);double this_endy=slc->reco.shw[largestShwIdx].start.y+(slc->reco.shw[largestShwIdx].dir.y *slc->reco.shw[largestShwIdx].len);double this_endz=slc->reco.shw[largestShwIdx].start.z+(slc->reco.shw[largestShwIdx].dir.z *slc->reco.shw[largestShwIdx].len);bool startx=(fvndAbs.xmin< std::abs(slc->reco.shw[largestShwIdx].start.x))&&(std::abs(slc->reco.shw[largestShwIdx].start.x)< fvndAbs.xmax);bool endx=(fvndAbs.xmin< std::abs(this_endx))&&(std::abs(this_endx)< fvndAbs.xmax);bool starty=(fvndAbs.ymin< slc->reco.shw[largestShwIdx].start.y)&&(slc->reco.shw[largestShwIdx].start.y< fvndAbs.ymax);bool endy=(fvndAbs.ymin< this_endy)&&(this_endy< fvndAbs.ymax);bool startz=(fvndAbs.zmin< slc->reco.shw[largestShwIdx].start.z)&&(slc->reco.shw[largestShwIdx].start.z< fvndAbs.zmax);bool endz=(fvndAbs.zmin< this_endz)&&(this_endz< fvndAbs.zmax);return(startx &&endx &&starty &&endy &&startz &&endz);}) |
const Cut | kNueContainedFD ([](const caf::SRSliceProxy *slc){const int largestShwIdx(kLargestRecoShowerIdx(slc));return false;double this_endx=slc->reco.shw[largestShwIdx].start.x+(slc->reco.shw[largestShwIdx].dir.x *slc->reco.shw[largestShwIdx].len);double this_endy=slc->reco.shw[largestShwIdx].start.y+(slc->reco.shw[largestShwIdx].dir.y *slc->reco.shw[largestShwIdx].len);double this_endz=slc->reco.shw[largestShwIdx].start.z+(slc->reco.shw[largestShwIdx].dir.z *slc->reco.shw[largestShwIdx].len);bool startx=(fvfd_cryo1.xmin< slc->reco.shw[largestShwIdx].start.x)&&(slc->reco.shw[largestShwIdx].start.x< fvfd_cryo1.xmax);bool endx=(fvfd_cryo1.xmin< this_endx)&&(this_endx< fvfd_cryo1.xmax);bool starty=(fvfd_cryo1.ymin< slc->reco.shw[largestShwIdx].start.y)&&(slc->reco.shw[largestShwIdx].start.y< fvfd_cryo1.ymax);bool endy=(fvfd_cryo1.ymin< this_endy)&&(this_endy< fvfd_cryo1.ymax);bool startz=(fvfd_cryo1.zmin< slc->reco.shw[largestShwIdx].start.z)&&(slc->reco.shw[largestShwIdx].start.z< fvfd_cryo1.zmax);bool endz=(fvfd_cryo1.zmin< this_endz)&&(this_endz< fvfd_cryo1.zmax);return(startx &&endx &&starty &&endy &&startz &&endz);}) |
const Cut | kNueMuonCutOLD |
const Cut | kNueMuonCutNEW |
const Cut | kPreNueSelND = kFiducialVolumeND && kSlcNuScoreCut && kSlcFlashMatchCut |
const Cut | kRecoNueSel = kRecoShower && kShowerEnergyCut |
const Cut | kFullNueSel = kNueTrackLenCut && kShowerConvGapCut && kShowerdEdxCut && kShowerDensityCut |
const Cut | kRazzleDazzleNueSel = kNueAllTrackDazzleMuonCut && kShowerRazzleElectronScoreCut |
const Cut | kNueFlashScoreFDCut = kSlcFlashMatchCut |
const Cut | kNuePandoraScoreFDCut = kSlcNuScoreCut |
const Cut | kRecoShowerFD = kRecoShower && kNueNumShowersCut && kShowerdEdxCut && kShowerConvGapCut && kNueTrackLenCut && kShowerDensityCut && kShowerEnergyCut |
const Cut | kNueFD = kContainedFD && kNueFlashScoreFDCut && kNuePandoraScoreFDCut && kRecoShowerFD |
const Cut | kNumuBasicQual ([](const caf::SRSliceProxy *slc){bool hastrk=(slc->reco.ntrk > 0);if(!hastrk) return hastrk;unsigned int muIdx=(unsigned int) kPrimMuonIdx(slc);double len=slc->reco.trk[muIdx].len;return(len > 0);}) |
const Cut | kHasFlashMatch ([](const caf::SRSliceProxy *slc){return(slc->fmatch.present);}) |
const Cut | kFlashMatchScore ([](const caf::SRSliceProxy *slc){return(slc->fmatch.time > 0 &&slc->fmatch.score > 6);}) |
const Cut | kFlashMatchNumuCut = kHasFlashMatch && kFlashMatchScore |
const Cut | kNumuTrkLen ([](const caf::SRSliceProxy *slc){bool hastrk=(slc->reco.ntrk > 0);if(!hastrk) return hastrk;unsigned int muIdx=(unsigned int) kPrimMuonIdx(slc);double len=slc->reco.trk[muIdx].len;return(len > 50);}) |
const Cut | kIsNuSlice = ( kTruthIndex >= 0.f ) |
const Cut | kIsCosmic = ( !kIsNuSlice ) |
const Cut | kIsNuMuCC ([](const caf::SRSliceProxy *slc){return(kIsNuSlice(slc)&&slc->truth.iscc &&(slc->truth.pdg==14||slc->truth.pdg==-14));}) |
const Cut | kIsNuOther = ( kIsNuSlice && !kIsNuMuCC ) |
const Cut | kCryo0 ([](const caf::SRSliceProxy *slc){return(!isnan(slc->vertex.x)&&slc->vertex.x< 0);}) |
const Cut | kTFiducial ([](const caf::SRSliceProxy *slc){return(!isnan(slc->truth.position.x)&&((slc->truth.position.x< -71.1-25 &&slc->truth.position.x >-369.33+25)||(slc->truth.position.x > 71.1+25 &&slc->truth.position.x< 369.33-25))&&!isnan(slc->truth.position.y)&&(slc->truth.position.y >-181.7+25 &&slc->truth.position.y< 134.8-25)&&!isnan(slc->truth.position.z)&&(slc->truth.position.z >-895.95+30 &&slc->truth.position.z< 895.95-50));}) |
const Cut | kRFiducial ([](const caf::SRSliceProxy *slc){return(!isnan(slc->vertex.x)&&((slc->vertex.x< -71.1-25 &&slc->vertex.x >-369.33+25)||(slc->vertex.x > 71.1+25 &&slc->vertex.x< 369.33-25))&&!isnan(slc->vertex.y)&&(slc->vertex.y >-181.7+25 &&slc->vertex.y< 134.8-25)&&!isnan(slc->vertex.z)&&(slc->vertex.z >-895.95+30 &&slc->vertex.z< 895.95-50));}) |
const Cut | kNotClearCosmic ([](const caf::SRSliceProxy *slc){return!slc->is_clear_cosmic;}) |
const Cut | kNuScore ([](const caf::SRSliceProxy *slc){return(!isnan(slc->nu_score)&&slc->nu_score > 0.4);}) |
const Cut | kFMScore ([](const caf::SRSliceProxy *slc){return(!isnan(slc->fmatch.score)&&slc->fmatch.score< 7.0);}) |
const Cut | kPTrack ([](const caf::SRSliceProxy *slc){return(kPTrackInd(slc) >=0);}) |
const Cut | kPTrackContained ([](const caf::SRSliceProxy *slc){int Ind=kPTrackInd(slc);bool Contained(false);if(Ind >=0){auto const &trk=slc->reco.trk.at(Ind);Contained=(!isnan(trk.end.x)&&(trk.end.x< -71.1-25 &&trk.end.x >-369.33+25)&&!isnan(trk.end.y)&&(trk.end.y >-181.7+25 &&trk.end.y< 134.8-25)&&!isnan(trk.end.z)&&(trk.end.z >-895.95+30 &&trk.end.z< 895.95-50));}return Contained;}) |
const Cut | kPTrackExiting ([](const caf::SRSliceProxy *slc){int Ind=kPTrackInd(slc);bool Exiting(false);if(Ind >=0){auto const &trk=slc->reco.trk.at(Ind);Exiting=!(!isnan(trk.end.x)&&(trk.end.x< -71.1-25 &&trk.end.x >-369.33+25)&&!isnan(trk.end.y)&&(trk.end.y >-181.7+25 &&trk.end.y< 134.8-25)&&!isnan(trk.end.z)&&(trk.end.z >-895.95+30 &&trk.end.z< 895.95-50));}return Exiting;}) |
const Cut | kNuMuCC_Cryo0 = ( kIsNuMuCC && kCryo0 ) |
const Cut | kCosmic_Cryo0 = ( kIsCosmic && kCryo0 ) |
const Cut | kNuOther_Cryo0 = ( kIsNuOther && kCryo0 ) |
const Cut | kNuMuCC_TFiducial = ( kNuMuCC_Cryo0 && kTFiducial ) |
const Cut | kNuMuCC_RFiducial = ( kNuMuCC_Cryo0 && kRFiducial ) |
const Cut | kCosmic_TFiducial = ( kCosmic_Cryo0 && kTFiducial ) |
const Cut | kCosmic_RFiducial = ( kCosmic_Cryo0 && kRFiducial ) |
const Cut | kNuOther_TFiducial = ( kNuOther_Cryo0 && kTFiducial ) |
const Cut | kNuOther_RFiducial = ( kNuOther_Cryo0 && kRFiducial ) |
const Cut | kNuMuCC_ClearCos = ( kNuMuCC_TFiducial && kNotClearCosmic ) |
const Cut | kCosmic_ClearCos = ( kCosmic_RFiducial && kNotClearCosmic ) |
const Cut | kNuOther_ClearCos = ( kNuOther_TFiducial && kNotClearCosmic ) |
const Cut | kNuMuCC_NuScore = ( kNuMuCC_ClearCos && kNuScore ) |
const Cut | kCosmic_NuScore = ( kCosmic_ClearCos && kNuScore ) |
const Cut | kNuOther_NuScore = ( kNuOther_ClearCos && kNuScore ) |
const Cut | kNuMuCC_FMScore = ( kNuMuCC_NuScore && kFMScore ) |
const Cut | kCosmic_FMScore = ( kCosmic_NuScore && kFMScore ) |
const Cut | kNuOther_FMScore = ( kNuOther_NuScore && kFMScore ) |
const Cut | kNuMuCC_PTrack = ( kNuMuCC_FMScore && kPTrack ) |
const Cut | kCosmic_PTrack = ( kCosmic_FMScore && kPTrack ) |
const Cut | kNuOther_PTrack = ( kNuOther_FMScore && kPTrack ) |
const Cut | kNuMuCC_FullSelection = ( kCryo0 && kRFiducial && kNotClearCosmic && kNuScore && kFMScore && kPTrack ) |
const SpillCut | kLongTrack ([](const caf::SRSpillProxy *sr){bool ProgCut=false;bool NuIsNuMuCC, IsMuon, Contained;for(auto const &nu:sr->mc.nu){NuIsNuMuCC=nu.iscc &&(nu.pdg==14||nu.pdg==-14);for(auto const &prim:nu.prim){IsMuon=(prim.pdg==13||prim.pdg==-13);Contained=(prim.contained==1);ProgCut=ProgCut||(NuIsNuMuCC &&IsMuon &&((Contained &&prim.length > 50.0)||(!Contained &&prim.length > 100.0)));}}return ProgCut;}) |
const Cut | kIcarus202208FMTimeCut = kFMTimeVar > 0 && kFMTimeVar < 1.8 |
const Cut | kIcarus202208FMScoreCut = kFMScoreVar < 9 |
const Cut | kIcarus202208LongTrackDirCut = kCRLongestTrackDirY > -0.91 |
const Cut | kIcarus202208FoundMuon = kIcarus202208MuonIdx >= 0 |
const Cut | kIcarus202208RecoFiducial ([](const caf::SRSliceProxy *slc){return(!isnan(slc->vertex.x)&&((slc->vertex.x< -71.1-25 &&slc->vertex.x >-369.33+25)||(slc->vertex.x > 71.1+25 &&slc->vertex.x< 369.33-25))&&!isnan(slc->vertex.y)&&(slc->vertex.y >-181.7+25 &&slc->vertex.y< 134.8-25)&&!isnan(slc->vertex.z)&&(slc->vertex.z >-895.95+30 &&slc->vertex.z< 895.95-50));}) |
const Cut | kIcarus202208NumuSelection = kIcarus202208RecoFiducial && kIcarus202208FMScoreCut && kIcarus202208FMTimeCut && kIcarus202208LongTrackDirCut && kIcarus202208FoundMuon |
const Cut | kIcarus202208NoPion = kIcarus202208NumPions == 0 |
const Cut | kIcarus202208ContainedHadrons ([](const caf::SRSliceProxy *slc){auto idx=kIcarus202208MuonIdx(slc);int muID=-1;if(idx >=0) muID=slc->reco.trk.at(idx).pfp.id;for(auto &trk:slc->reco.trk){if(!Icarus202208contained(trk)) return false;}return true;}) |
const Cut | kIcarus202208ContainedMuon ([](const caf::SRSliceProxy *slc){return kIcarus202208FoundMuon(slc)&&Icarus202208contained(slc->reco.trk.at(kIcarus202208MuonIdx(slc)));}) |
const Cut | kIcarus202208ContainedMuonAndHadrons = kIcarus202208ContainedMuon && kIcarus202208ContainedHadrons |
const Cut | kIcarus202208QELike = kIcarus202208NumuSelection && kIcarus202208NoPion && kIcarus202208ContainedHadrons |
const Cut | kIcarus202208QELikeContainedMuon = kIcarus202208QELike && kIcarus202208ContainedMuon |
const Cut | kInFV ([](const caf::SRSliceProxy *slc){return(slc->vertex.x >(-199.15+10)&&slc->vertex.x< (199.15-10)&&slc->vertex.y >(-200.+10)&&slc->vertex.y< (200.-10)&&slc->vertex.z >(0.0+10)&&slc->vertex.z< (500.-50));}) |
const Cut | kSlcFlashMatchTimeCut ([](const caf::SRSliceProxy *slc){if(std::isnan(slc->fmatch.time)) return false;bool InBeam=(slc->fmatch.time > 0.&&slc->fmatch.time< 1.800);return(InBeam);}) |
const Cut | kSlcFlashMatchScoreCut ([](const caf::SRSliceProxy *slc){return(slc->fmatch.score< 7);}) |
const Cut | kHasPrimaryMuonTrk ([](const caf::SRSliceProxy *slc){int ptrkid=kPrimaryMuonTrkIdx(slc);double ptrkrecop=kPrimaryMuonTrkP(slc);return(ptrkid!=-1 &&!isnan(ptrkrecop)&&ptrkrecop > 0.&&ptrkrecop< 7.5);}) |
const Cut | kCRTTrackAngleCut ([](const caf::SRSliceProxy *slc){int ptrkid=kPrimaryMuonTrkIdx(slc);const caf::SRTrackProxy &ptrk=slc->reco.trk[ptrkid];return(isnan(ptrk.crttrack.angle)||ptrk.crttrack.angle > 0.05);}) |
const Cut | kCRTHitDistanceCut ([](const caf::SRSliceProxy *slc){int ptrkid=kPrimaryMuonTrkIdx(slc);const caf::SRCRTHitMatchProxy &crthit=slc->reco.trk[ptrkid].crthit;if(std::isnan(crthit.hit.time)||std::isnan(crthit.distance)) return true;return crthit.distance > 5||(crthit.hit.time > 0.&&crthit.hit.time< 1.800);}) |
const Cut | kInBeamSpill |
Does the event fall inside the beam spill window. More... | |
const Cut | kInTimingSideband |
Is the event far from the start and end of the spill window. More... | |
const Cut | kTrueActiveVolumeND ([](const caf::SRSliceProxy *slc){return kHasNu(slc)&&PtInVolAbsX(slc->truth.position, avnd);}) |
const Cut | kTrueFiducialVolumeND ([](const caf::SRSliceProxy *slc){return kHasNu(slc)&&PtInVolAbsX(slc->truth.position, fvndAbs);}) |
const Cut | kTrueActiveVolumeFDCryo1 ([](const caf::SRSliceProxy *slc){return kHasNu(slc)&&PtInVol(slc->truth.position, avfd_cryo1);}) |
const Cut | kTrueActiveVolumeFDCryo2 ([](const caf::SRSliceProxy *slc){return kHasNu(slc)&&PtInVol(slc->truth.position, avfd_cryo2);}) |
const Cut | kTrueFiducialVolumeFDCryo1 ([](const caf::SRSliceProxy *slc){return kHasNu(slc)&&PtInVol(slc->truth.position, fvfd_cryo1);}) |
const Cut | kTrueFiducialVolumeFDCryo2 ([](const caf::SRSliceProxy *slc){return kHasNu(slc)&&PtInVol(slc->truth.position, fvfd_cryo2);}) |
const Cut | kIsNu ([](const caf::SRSliceProxy *slc){if(slc->truth.index< 0) return false;return slc->truth.pdg > 0;}) |
const Cut | kHasNu ([](const caf::SRSliceProxy *slc){return slc->truth.index >=0;}) |
const Cut | kIsNue ([](const caf::SRSliceProxy *slc){return slc->truth.index >=0 &&abs(slc->truth.pdg)==12;}) |
const Cut | kIsNumu ([](const caf::SRSliceProxy *slc){return slc->truth.index >=0 &&abs(slc->truth.pdg)==14;}) |
const Cut | kIsNutau ([](const caf::SRSliceProxy *slc){return slc->truth.index >=0 &&abs(slc->truth.pdg)==16;}) |
const Cut | kVtxDistMagCut ([](const caf::SRSliceProxy *slc){if(slc->truth.index< 0) return true;return(kTruthVtxDistMag(slc)< 1);}) |
const Cut | kSlcCompletenessCut ([](const caf::SRSliceProxy *slc){if(slc->truth.index< 0) return false;return(kCompletness(slc) > 0.5);}) |
const SpillCut | kIsCosmicSpill ([](const caf::SRSpillProxy *sr){return(sr->mc.nnu==0);}) |
const SpillCut | kIsSingleNuSpill ([](const caf::SRSpillProxy *sr){return(sr->mc.nnu< 2);}) |
const SpillCut | kIsNueSpill ([](const caf::SRSpillProxy *sr){if(kIsCosmicSpill(sr)||!kIsSingleNuSpill(sr)) return false;return(std::abs(sr->mc.nu[0].pdg)==12);}) |
const SpillCut | kIsNumuSpill ([](const caf::SRSpillProxy *sr){if(kIsCosmicSpill(sr)||!kIsSingleNuSpill(sr)) return false;return(std::abs(sr->mc.nu[0].pdg)==14);}) |
const SpillCut | kIsNutauSpill ([](const caf::SRSpillProxy *sr){if(kIsCosmicSpill(sr)||!kIsSingleNuSpill(sr)) return false;return(std::abs(sr->mc.nu[0].pdg)==16);}) |
const SpillCut | kIsCCSpill ([](const caf::SRSpillProxy *sr){if(kIsCosmicSpill(sr)||!kIsSingleNuSpill(sr)) return false;return((bool) sr->mc.nu[0].iscc);}) |
const SpillCut | kIsNCSpill ([](const caf::SRSpillProxy *sr){if(kIsCosmicSpill(sr)||!kIsSingleNuSpill(sr)) return false;return((bool) sr->mc.nu[0].isnc);}) |
const std::vector< double > | kLowEnergyEdges |
const std::vector< double > | kLowEnergyGeVEdges |
const Binning | kNueEnergyBinning = Binning::Simple( 20, 0., 5000.) |
const Binning | kNumuEnergyBinning = Binning::Simple( 20, 0., 5000.) |
const Binning | kLowEnergyBinning = Binning::Custom( kLowEnergyEdges ) |
const Binning | kLowEnergyGeVBinning = Binning::Custom( kLowEnergyGeVEdges ) |
const Binning | kPositionXNDBinning = Binning::Simple(40, avnd.xmin, avnd.xmax) |
const Binning | kPositionYNDBinning = Binning::Simple(40, avnd.ymin, avnd.ymax) |
const Binning | kPositionZNDBinning = Binning::Simple(50, avnd.zmin, avnd.zmax) |
const Binning | kPositionXFDBinning = Binning::Simple(47, avfd_cryo1.xmin, avfd_cryo1.xmax) |
const Binning | kPositionYFDBinning = Binning::Simple(35, avfd_cryo1.ymin, avfd_cryo1.ymax) |
const Binning | kPositionZFDBinning = Binning::Simple(90, avfd_cryo1.zmin, avfd_cryo1.zmax) |
const Binning | kCRTXFDBinning = Binning::Simple(60, crtfd.xmin, crtfd.xmax) |
const Binning | kCRTYFDBinning = Binning::Simple(100, crtfd.ymin, crtfd.ymax) |
const Binning | kCRTZFDBinning = Binning::Simple(130, crtfd.zmin, crtfd.zmax) |
const Binning | kBDTBinning = Binning::Simple(50,0, 1.0) |
const Var | kLargestRecoShowerIdx ([](const caf::SRSliceProxy *slc) -> int{int bestIdx(-1);double maxEnergy(-1);for(unsigned int i=0;i< slc->reco.nshw;i++){auto const &shw=slc->reco.shw[i];continue;if(shw.bestplane_energy > maxEnergy){bestIdx=i;maxEnergy=shw.bestplane_energy;}}return bestIdx;}) |
const Var | kRecoShower_BestEnergy ([](const caf::SRSliceProxy *slc) -> double{const caf::SRShowerProxy *shw=LargestRecoShower(slc);return shw?double(shw->bestplane_energy):-5;}) |
const Var | kRecoShower_TruePdg ([](const caf::SRSliceProxy *slc) -> int{const caf::SRShowerProxy *shw=LargestRecoShower(slc);return shw?int(shw->truth.p.pdg):-5;}) |
const Var | kRecoShower_BestdEdx ([](const caf::SRSliceProxy *slc) -> double{const caf::SRShowerProxy *shw=LargestRecoShower(slc);if(!shw||shw->bestplane_dEdx< 0) return-5.;return shw->bestplane_dEdx;}) |
const Var | kRecoShower_ConversionGap ([](const caf::SRSliceProxy *slc) -> double{const caf::SRShowerProxy *shw=LargestRecoShower(slc);return shw?double(shw->conversion_gap):-5.;}) |
const Var | kRecoShower_Density ([](const caf::SRSliceProxy *slc) -> double{const caf::SRShowerProxy *shw=LargestRecoShower(slc);return shw?double(shw->density):-5.;}) |
const Var | kRecoShower_Energy ([](const caf::SRSliceProxy *slc) -> double{const caf::SRShowerProxy *shw=LargestRecoShower(slc);return shw?double(shw->plane[1].energy):-5.;}) |
const Var | kRecoShower_Length ([](const caf::SRSliceProxy *slc) -> double{const caf::SRShowerProxy *shw=LargestRecoShower(slc);return shw?double(shw->len):-5.;}) |
const Var | kRecoShower_OpenAngle ([](const caf::SRSliceProxy *slc) -> double{const caf::SRShowerProxy *shw=LargestRecoShower(slc);return shw?180.*shw->open_angle/M_PI:-5.;}) |
const Var | kRecoShower_StartX ([](const caf::SRSliceProxy *slc) -> double{const caf::SRShowerProxy *shw=LargestRecoShower(slc);return shw?double(shw->start.x):-9999.;}) |
const Var | kRecoShower_StartY ([](const caf::SRSliceProxy *slc) -> double{const caf::SRShowerProxy *shw=LargestRecoShower(slc);return shw?double(shw->start.y):-9999.;}) |
const Var | kRecoShower_StartZ ([](const caf::SRSliceProxy *slc) -> double{const caf::SRShowerProxy *shw=LargestRecoShower(slc);return shw?double(shw->start.z):-9999.;}) |
const Var | kRecoShower_EndX ([](const caf::SRSliceProxy *slc) -> double{const caf::SRShowerProxy *shw=LargestRecoShower(slc);return shw?double(shw->start.x+shw->dir.x *shw->len):-9999.;}) |
const Var | kRecoShower_EndY ([](const caf::SRSliceProxy *slc) -> double{const caf::SRShowerProxy *shw=LargestRecoShower(slc);return shw?double(shw->start.y+shw->dir.y *shw->len):-9999.;}) |
const Var | kRecoShower_EndZ ([](const caf::SRSliceProxy *slc) -> double{const caf::SRShowerProxy *shw=LargestRecoShower(slc);return shw?double(shw->start.z+shw->dir.z *shw->len):-9999.;}) |
const Var | kRecoShower_densityGradient ([](const caf::SRSliceProxy *slc) -> double{const caf::SRShowerProxy *shw=LargestRecoShower(slc);return shw?double(shw->selVars.densityGradient):-5.;}) |
const Var | kRecoShower_densityGradientPower ([](const caf::SRSliceProxy *slc) -> double{const caf::SRShowerProxy *shw=LargestRecoShower(slc);return shw?double(shw->selVars.densityGradientPower):-5.;}) |
const Var | kRecoShower_trackLength ([](const caf::SRSliceProxy *slc) -> double{const caf::SRShowerProxy *shw=LargestRecoShower(slc);return shw?double(shw->selVars.trackLength):-5.;}) |
const Var | kRecoShower_trackWidth ([](const caf::SRSliceProxy *slc) -> double{const caf::SRShowerProxy *shw=LargestRecoShower(slc);return shw?double(shw->selVars.trackWidth):-5.;}) |
const Var | kRecoShowers_EnergyCut ([](const caf::SRSliceProxy *slc) -> unsigned{unsigned int counter(0);for(auto const &shw:slc->reco.shw){++counter;}return counter;}) |
const Var | kLongestTrackIdx ([](const caf::SRSliceProxy *slc) -> int{int bestIdx(-1);double maxLength(-1);for(unsigned int i=0;i< slc->reco.ntrk;i++){auto const &trk=slc->reco.trk[i];continue;if(trk.len > maxLength){bestIdx=i;maxLength=trk.len;}}return bestIdx;}) |
const Var | kLongestTrackTruePdg ([](const caf::SRSliceProxy *slc) -> int{const caf::SRTrackProxy *trk=LongestRecoTrack(slc);return trk?int(trk->truth.p.pdg):-5;}) |
const Var | kLongestTrackLength ([](const caf::SRSliceProxy *slc) -> double{const caf::SRTrackProxy *trk=LongestRecoTrack(slc);return trk?double(trk->len):-5.;}) |
const Var | kLongestTrackChi2Muon ([](const caf::SRSliceProxy *slc) -> double{const caf::SRTrkChi2PIDProxy *chi2=LongestTrackBestPlaneChi2PID(slc);return chi2?double(chi2->chi2_muon):-5.;}) |
const Var | kLongestTrackChi2Pion ([](const caf::SRSliceProxy *slc) -> double{const caf::SRTrkChi2PIDProxy *chi2=LongestTrackBestPlaneChi2PID(slc);return chi2?double(chi2->chi2_pion):-5.;}) |
const Var | kLongestTrackChi2Kaon ([](const caf::SRSliceProxy *slc) -> double{const caf::SRTrkChi2PIDProxy *chi2=LongestTrackBestPlaneChi2PID(slc);return chi2?double(chi2->chi2_kaon):-5.;}) |
const Var | kLongestTrackChi2Proton ([](const caf::SRSliceProxy *slc) -> double{const caf::SRTrkChi2PIDProxy *chi2=LongestTrackBestPlaneChi2PID(slc);return chi2?double(chi2->chi2_proton):-5.;}) |
const Var | kMuonTrackLength ([](const caf::SRSliceProxy *slc) -> double{const caf::SRTrackProxy *trk=LongestRecoTrack(slc);if(!trk) return-5;if(trk &&(kLongestTrackChi2Muon(slc)< 30.f &&kLongestTrackChi2Proton(slc) > 60.f)){return trk->len;}return-5;}) |
const Var | kLongestTrackDazzlePID ([](const caf::SRSliceProxy *slc) -> int{const caf::SRTrackProxy *trk=LongestRecoTrack(slc);return trk?(int) trk->dazzle.pdg:-5;}) |
const Var | kLongestTrackDazzleMuonScore ([](const caf::SRSliceProxy *slc) -> float{const caf::SRTrackProxy *trk=LongestRecoTrack(slc);return trk?(float) trk->dazzle.muonScore:-5.f;}) |
const Var | kRecoShowerRazzlePID ([](const caf::SRSliceProxy *slc) -> int{const caf::SRShowerProxy *shw=LargestRecoShower(slc);return shw?(int) shw->razzle.pdg:-5;}) |
const Var | kRecoShowerRazzleElectronScore ([](const caf::SRSliceProxy *slc) -> float{const caf::SRShowerProxy *shw=LargestRecoShower(slc);return shw?(float) shw->razzle.electronScore:-5.f;}) |
const Var | kPrimMuonIdx ([](const caf::SRSliceProxy *slc) -> double{if((int) slc->reco.ntrk==0) return-5.0;double best_idx=0;double best_len=-5.0;for(unsigned int trkIdx=0;trkIdx< slc->reco.ntrk;trkIdx++){auto &trk=slc->reco.trk[trkIdx];if(trk.chi2pid[2].pid_ndof< 0) return-5.0;bool isMuonLike=trk.chi2pid[2].chi2_pion > trk.chi2pid[2].chi2_muon;if(isMuonLike &&trk.len > best_len){best_len=trk.len;best_idx=trkIdx;}}return best_idx;}) |
const Var | kPrimTrkLen ([](const caf::SRSliceProxy *slc) -> double{double len=-5.0;if(slc->reco.ntrk > 0){int muIdx=(int) kPrimMuonIdx(slc);if(muIdx >=0){len=slc->reco.trk[muIdx].len;}}return len;}) |
const Var | kTruthIndex = SIMPLEVAR(truth.index) |
const Var | kPrimaryEnergy = SIMPLEVAR(truth.E) |
const Var | kMuMaxTrack ([](const caf::SRSliceProxy *slc) -> float{float len(-5.f);for(auto const &trk:slc->reco.trk){if((trk.truth.p.pdg==13||trk.truth.p.pdg==-13)&&trk.truth.p.length > len) len=trk.truth.p.length;}return len;}) |
const Var | kPTrackInd ([](const caf::SRSliceProxy *slc) -> int{float Longest(0);int PTrackInd(-1);for(std::size_t i(0);i< slc->reco.trk.size();++i){auto const &trk=slc->reco.trk.at(i);if(trk.bestplane==-1) continue;const float Atslc=std::hypot(slc->vertex.x-trk.start.x, slc->vertex.y-trk.start.y, slc->vertex.z-trk.start.z);const bool AtSlice=(Atslc< 10.0 &&trk.pfp.parent_is_primary);const float Chi2Proton=trk.chi2pid[trk.bestplane].chi2_proton;const float Chi2Muon=trk.chi2pid[trk.bestplane].chi2_muon;const bool Contained=(!isnan(trk.end.x)&&(trk.end.x< -71.1-25 &&trk.end.x >-369.33+25)&&!isnan(trk.end.y)&&(trk.end.y >-181.7+25 &&trk.end.y< 134.8-25)&&!isnan(trk.end.z)&&(trk.end.z >-895.95+30 &&trk.end.z< 895.95-50));const bool MaybeMuonExiting=(!Contained &&trk.len > 100);const bool MaybeMuonContained=(Contained &&Chi2Proton > 60 &&Chi2Muon< 30 &&trk.len > 50);{Longest=trk.len;PTrackInd=i;}}return PTrackInd;}) |
const Var | kRecoMuonP ([](const caf::SRSliceProxy *slc) -> float{float p(-5.f);{auto const &trk=slc->reco.trk.at(kPTrackInd(slc));const bool Contained=(!isnan(trk.end.x)&&(trk.end.x< -71.1-25 &&trk.end.x >-369.33+25)&&!isnan(trk.end.y)&&(trk.end.y >-181.7+25 &&trk.end.y< 134.8-25)&&!isnan(trk.end.z)&&(trk.end.z >-895.95+30 &&trk.end.z< 895.95-50));if(Contained) p=trk.rangeP.p_muon;else p=trk.mcsP.fwdP_muon;}return p;}) |
const Var | kTrueMuonP ([](const caf::SRSliceProxy *slc) -> float{float p(-5.f);{auto const &trk=slc->reco.trk.at(kPTrackInd(slc));p=std::hypot(trk.truth.p.genp.x, trk.truth.p.genp.y, trk.truth.p.genp.z);}return p;}) |
const Var | kCRLongestTrackDirY = SIMPLEVAR(nuid.crlongtrkdiry) |
const Var | kFMTimeVar = SIMPLEVAR(fmatch.time) |
const Var | kFMScoreVar = SIMPLEVAR(fmatch.score) |
const Var | kIcarus202208MuonIdx ([](const caf::SRSliceProxy *slc) -> int{float Longest(0);int PTrackInd(-1);for(std::size_t i(0);i< slc->reco.trk.size();++i){auto const &trk=slc->reco.trk.at(i);if(trk.bestplane==-1) continue;const float Atslc=std::hypot(slc->vertex.x-trk.start.x, slc->vertex.y-trk.start.y, slc->vertex.z-trk.start.z);const bool AtSlice=(Atslc< 10.0 &&trk.pfp.parent_is_primary);const float Chi2Proton=trk.chi2pid[trk.bestplane].chi2_proton;const float Chi2Muon=trk.chi2pid[trk.bestplane].chi2_muon;const bool Contained=(!isnan(trk.end.x)&&((trk.end.x< -71.1-5 &&trk.end.x >-369.33+5)||(trk.end.x< 71.1+5 &&trk.end.x >+369.33-5))&&!isnan(trk.end.y)&&(trk.end.y >-181.7+5 &&trk.end.y< 134.8-5)&&!isnan(trk.end.z)&&(trk.end.z >-895.95+5 &&trk.end.z< 895.95-5));const bool MaybeMuonExiting=(!Contained &&trk.len > 100);const bool MaybeMuonContained=(Contained &&Chi2Proton > 60 &&Chi2Muon< 30 &&trk.len > 50);if(AtSlice &&(MaybeMuonExiting||MaybeMuonContained)&&trk.len > Longest){Longest=trk.len;PTrackInd=i;}}return PTrackInd;}) |
const Var | kIcarus202208NumPions ([](const caf::SRSliceProxy *slc){int count=0;auto idx=kIcarus202208MuonIdx(slc);int muID=-1;if(idx >=0) muID=slc->reco.trk.at(idx).pfp.id;for(auto &trk:slc->reco.trk){++count;}return count;}) |
const Var | kPrimaryMuonTrkIdx ([](const caf::SRSliceProxy *slc){double longest=-1;int best_idx=-1;double dist=-1;bool atslc, contained, maybe_muon_exiting, maybe_muon_contained;float chi2_proton, chi2_muon;for(unsigned int trkidx=0;trkidx< slc->reco.trk.size();++trkidx){const caf::SRTrackProxy &trk=slc->reco.trk[trkidx];if(trk.bestplane==-1) continue;dist=sqrt((pow(trk.start.x-slc->vertex.x, 2))+(pow(trk.start.y-slc->vertex.y, 2))+(pow(trk.start.z-slc->vertex.z, 2)));atslc=dist< 10;if(!atslc||!trk.pfp.parent_is_primary) continue;contained=((-199.15+10)< trk.end.x &&trk.end.x< (199.15-10)&&(-200.+10)< trk.end.y &&trk.end.y< (200.-10)&&(0.0+10)< trk.end.z &&trk.end.z< (500.-50));chi2_proton=trk.chi2pid[trk.bestplane].chi2_proton;chi2_muon=trk.chi2pid[trk.bestplane].chi2_muon;maybe_muon_exiting=!contained &&trk.len > 100;maybe_muon_contained=contained &&chi2_proton > 60 &&chi2_muon< 30 &&trk.len > 50;if(!maybe_muon_contained &&!maybe_muon_exiting) continue;if(trk.len > longest){longest=trk.len;best_idx=trkidx;}}return best_idx;}) |
const Var | kPrimaryMuonTrkP ([](const caf::SRSliceProxy *slc){float recop(-5.f);bool contained(false);if(kPrimaryMuonTrkIdx(slc) >=0){auto const &ptrk=slc->reco.trk.at(kPrimaryMuonTrkIdx(slc));contained=((-199.15+10)< ptrk.end.x &&ptrk.end.x< (199.15-10)&&(-200.+10)< ptrk.end.y &&ptrk.end.y< (200.-10)&&(0.0+10)< ptrk.end.z &&ptrk.end.z< (500.-50));if(contained) recop=ptrk.rangeP.p_muon;else recop=ptrk.mcsP.fwdP_muon;}return recop;}) |
const Var | kCRTTrkTime ([](const caf::SRSliceProxy *slc) -> double{float crttrktime(-5.f);if(kPrimaryMuonTrkIdx(slc) >=0){int ptrkid=kPrimaryMuonTrkIdx(slc);const caf::SRTrackProxy &ptrk=slc->reco.trk[ptrkid];crttrktime=ptrk.crthit.hit.time;}return crttrktime;}) |
const Var | kCRTTrkAngle ([](const caf::SRSliceProxy *slc) -> double{float crttrkangle(-5.f);if(kPrimaryMuonTrkIdx(slc) >=0){int ptrkid=kPrimaryMuonTrkIdx(slc);const caf::SRTrackProxy &ptrk=slc->reco.trk[ptrkid];crttrkangle=ptrk.crttrack.angle;}return crttrkangle;}) |
const Var | kCRTHitDist ([](const caf::SRSliceProxy *slc) -> double{float crttrkdist(-5.f);if(kPrimaryMuonTrkIdx(slc) >=0){int ptrkid=kPrimaryMuonTrkIdx(slc);const caf::SRTrackProxy &ptrk=slc->reco.trk[ptrkid];crttrkdist=ptrk.crthit.distance;}return crttrkdist;}) |
const Var | kPrimaryMuonTrkLen ([](const caf::SRSliceProxy *slc) -> double{float ptrklen(-5.f);if(kPrimaryMuonTrkIdx(slc) >=0){int ptrkid=kPrimaryMuonTrkIdx(slc);const caf::SRTrackProxy &ptrk=slc->reco.trk[ptrkid];ptrklen=ptrk.len;}return ptrklen;}) |
const Var | kHasTruthMatch ([](const caf::SRSliceProxy *slc) -> double{return(slc->truth.index!=-1);}) |
const Var | kCompletness ([](const caf::SRSliceProxy *slc) -> double{return(kHasTruthMatch(slc)?(float) slc->tmatch.eff:-5.f);}) |
const Var | kTruthEnergy ([](const caf::SRSliceProxy *slc) -> double{return(kHasTruthMatch(slc)?(float) slc->truth.E:-5.f);}) |
const Var | kTruthVtxX ([](const caf::SRSliceProxy *slc) -> double{return(kHasTruthMatch(slc)?(float) slc->truth.position.x:-999.f);}) |
const Var | kTruthVtxY ([](const caf::SRSliceProxy *slc) -> double{return(kHasTruthMatch(slc)?(float) slc->truth.position.y:-999.f);}) |
const Var | kTruthVtxZ ([](const caf::SRSliceProxy *slc) -> double{return(kHasTruthMatch(slc)?(float) slc->truth.position.z:-999.f);}) |
const Var | kTruthVtxDistX ([](const caf::SRSliceProxy *slc) -> double{return(kHasTruthMatch(slc)?(float) std::abs(slc->truth.position.x-kSlcVtxX(slc)):-999.f);}) |
const Var | kTruthVtxDistY ([](const caf::SRSliceProxy *slc) -> double{return(kHasTruthMatch(slc)?(float) std::abs(slc->truth.position.y-kSlcVtxY(slc)):-999.f);}) |
const Var | kTruthVtxDistZ ([](const caf::SRSliceProxy *slc) -> double{return(kHasTruthMatch(slc)?(float) std::abs(slc->truth.position.z-kSlcVtxZ(slc)):-999.f);}) |
const Var | kTruthVtxDistMag ([](const caf::SRSliceProxy *slc) -> double{return(kHasTruthMatch(slc)?(float) std::hypot(kTruthVtxDistX(slc), kTruthVtxDistY(slc), kTruthVtxDistZ(slc)):-5.f);}) |
const SpillVar | kTruthNuEnergy ([](const caf::SRSpillProxy *sr) -> double{return(sr->mc.nnu!=1?-5.f:(float) sr->mc.nu[0].E);}) |
const SpillVar | kTruthLeptonEnergy ([](const caf::SRSpillProxy *sr) -> double{if(sr->mc.nnu!=1) return-5.f;for(auto const &prim:sr->mc.nu[0].prim){if(std::abs(prim.pdg)!=11 &&std::abs(prim.pdg)!=13) continue;return prim.startE;}return-5.f;}) |
const SpillVar | kRun = SIMPLESPILLVAR(hdr.run) |
const SpillVar | kEvt = SIMPLESPILLVAR(hdr.evt) |
const Var | kCounting = kUnweighted |
const SpillVar | kSpillCounting = kSpillUnweighted |
const SpillMultiVar | kCRTHitX ([](const caf::SRSpillProxy *sr){std::vector< double > positions;for(const auto &hit:sr->crt_hits){positions.push_back(hit.position.x);}return positions;}) |
const SpillMultiVar | kCRTHitY ([](const caf::SRSpillProxy *sr){std::vector< double > positions;for(const auto &hit:sr->crt_hits){positions.push_back(hit.position.y);}return positions;}) |
const SpillMultiVar | kCRTHitZ ([](const caf::SRSpillProxy *sr){std::vector< double > positions;for(const auto &hit:sr->crt_hits){positions.push_back(hit.position.z);}return positions;}) |
const SpillMultiVar | kCRTHitPE ([](const caf::SRSpillProxy *sr){std::vector< double > pes;for(const auto &hit:sr->crt_hits){pes.push_back(hit.pe);}return pes;}) |
const SpillMultiVar | kCRTHitTime ([](const caf::SRSpillProxy *sr){std::vector< double > times;for(const auto &hit:sr->crt_hits){times.push_back(hit.time);}return times;}) |
const SpillMultiVar | kCRTHitTimeFD ([](const caf::SRSpillProxy *sr){std::vector< double > times;for(const auto &hit:sr->crt_hits){auto thistime=hit.time-1600.;times.push_back(thistime);}return times;}) |
const Var | kSlcVtxX ([](const caf::SRSliceProxy *slc) -> double{return slc->vertex.x;}) |
const Var | kSlcVtxY ([](const caf::SRSliceProxy *slc) -> double{return slc->vertex.y;}) |
const Var | kSlcVtxZ ([](const caf::SRSliceProxy *slc) -> double{return slc->vertex.z;}) |
const Var | kSlcNuScore ([](const caf::SRSliceProxy *slc) -> double{return slc->nu_score;}) |
const Var | kSlcHasFlash ([](const caf::SRSliceProxy *slc) -> double{return slc->fmatch.present;}) |
const Var | kSlcFlashScore ([](const caf::SRSliceProxy *slc) -> double{return((bool) kSlcHasFlash(slc)?(float) slc->fmatch.score:-5.f);}) |
Oscillation analysis framework, runs over CAF files outside of ART.
Beam spill duration is 1.6 microseconds TO DO: Change numbers to reflect reality: Where the beam spill actually starts after the beam trigger?
All analysis code is defined in namespace "ana"
typedef _Cut<caf::SRSliceProxy> ana::Cut |
typedef double( ana::ExposureFunc_t)(const caf::SRSpill *spill) |
typedef _HistAxis<Var> ana::HistAxis |
Definition at line 62 of file HistAxis.h.
typedef _MultiVar<caf::SRSliceProxy> ana::MultiVar |
Definition at line 48 of file MultiVar.h.
typedef _Cut<caf::SRSliceProxy> ana::SliceCut |
Representation of a cut (selection) to be applied to a caf::StandardRecord object.
A Cut consists of a function, taking a StandardRecord and returning a boolean indicating if that event passes the cut.
Cut objects may be combined with the standard boolean operations && || and !
typedef _Cut<caf::SRSpillProxy> ana::SpillCut |
typedef _HistAxis<SpillVar> ana::SpillHistAxis |
Definition at line 63 of file HistAxis.h.
typedef _MultiVar<caf::SRSpillProxy> ana::SpillMultiVar |
Definition at line 49 of file MultiVar.h.
typedef _Cut<caf::SRSpillTruthBranch> ana::SpillTruthCut |
typedef _Var<caf::SRSpillProxy> ana::SpillVar |
typedef _Var<caf::SRSliceProxy> ana::Var |
Representation of a variable to be retrieved from a caf::StandardRecord object.
A Var consists of a function, taking a StandardRecord and returning the value of the variable (which may be some complicated function).
enum ana::DataSource |
Is this data-file representing beam spills or cosmic spills?
Enumerator | |
---|---|
kBeam | |
kCosmic |
Definition at line 28 of file SpectrumLoaderBase.h.
enum ana::EBinType |
Enumerator | |
---|---|
kBinContent |
Regular histogram. |
kBinDensity |
Divide bin contents by bin widths. |
Definition at line 36 of file sbnana/sbnana/CAFAna/Core/Utilities.h.
enum ana::EExposureType |
For use as an argument to Spectrum::ToTH1.
Enumerator | |
---|---|
kPOT | |
kLivetime |
Definition at line 44 of file sbnana/sbnana/CAFAna/Core/Utilities.h.
|
strong |
|
strong |
Enumerator | |
---|---|
kMuon | |
kHadron | |
kNeutron | |
kEM | |
kChargedHadron |
Definition at line 18 of file EnergySysts.h.
|
strong |
Enumerator | |
---|---|
kConstant | |
kSqrt | |
kInverseSqrt |
Definition at line 12 of file EnergySysts.h.
|
strong |
Desired match type in UniverseOracle::ClosestShiftIndex.
Enumerator | |
---|---|
kAbove | |
kBelow | |
kEither |
Definition at line 12 of file UniverseOracle.h.
|
strong |
Enumerator | |
---|---|
kNone | |
kSinSq2ThetaMuMu | |
kSinSq2ThetaMuE | |
kSinSq2ThetaEE |
Definition at line 7 of file OscCalcSterileApprox.h.
bool ana::AlmostEqual | ( | double | a, |
double | b | ||
) |
Definition at line 799 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
TLegend * ana::AutoPlaceLegend | ( | double | dx, |
double | dy, | ||
double | yPin = -1 |
||
) |
Create a legend, maximizing distance from all histograms.
dx | Width of the legend, fraction of pad width |
dy | Height of the legend, fraction of pad height |
yPin,height | in NDC (frac) to pin center of legend, i.e only move x |
Definition at line 510 of file Plots.cxx.
double ana::AvgSinSq | ( | double | k, |
double | a, | ||
double | b | ||
) |
Definition at line 145 of file OscCalcSterileApprox.cxx.
const caf::SRTrkChi2PIDProxy* ana::BestPlaneChi2PID | ( | const caf::SRTrackProxy * | trk | ) |
Definition at line 193 of file NueVars.cxx.
|
inline |
Definition at line 32 of file sbnana/sbnana/CAFAna/Analysis/Surface.cxx.
std::unique_ptr< TMatrixD > ana::CalcCovMx | ( | const std::vector< TArrayD * > & | binSets, |
int | firstBin = 0 , |
||
int | lastBin = -1 |
||
) |
Compute bin-to-bin covariance matrix from a collection of sets of bin contents.
binSets | Collection of sets of bins from which covariances should be calculated |
firstBin | The first bin that should be considered (inclusive) |
lastBin | The last bin that should be considered (inclusive). -1 means "last in set" |
Note TH1D is a child class of TArrayD – so you can pass a vector of TH1D* to this method.
Definition at line 118 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
double ana::Chi2CovMx | ( | const TVectorD * | e, |
const TVectorD * | o, | ||
const TMatrixD * | covmxinv | ||
) |
Definition at line 221 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
double ana::Chi2CovMx | ( | const TH1 * | e, |
const TH1 * | o, | ||
const TMatrixD * | covmxinv | ||
) |
Definition at line 230 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
void ana::ClearTrueParticles | ( | caf::StandardRecord * | sr | ) |
Definition at line 25 of file FileReducer.cxx.
std::function<ExposureFunc_t> ana::CombineExposures | ( | const std::function< ExposureFunc_t > & | a, |
const std::function< ExposureFunc_t > & | b | ||
) |
void ana::CombineMetadata | ( | std::map< std::string, std::string > & | base, |
const std::map< std::string, std::string > & | add, | ||
std::set< std::string > & | mask | ||
) |
Definition at line 612 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
Var ana::Constant | ( | double | c | ) |
void ana::CountingExperimentErrorBarChart | ( | const std::map< std::string, double > & | systs, |
double | statErr = 0 , |
||
bool | bkgdOrSig = false , |
||
bool | shortchart = false |
||
) |
Make a simple plot of relative size of different errors.
systs | Map from error label to size in percent |
statErr | Percentage statistical error. Optional. |
bkgdOrSig | Label axis for background or signal. Default (false) = bkgd |
void ana::DataMCAreaNormalizedRatio | ( | const Spectrum & | data, |
const IPrediction * | mc, | ||
osc::IOscCalc * | calc, | ||
double | miny, | ||
double | maxy | ||
) |
Plot data/MC ratio for the given spectrum. Normalize MC to Data by area.
Definition at line 328 of file Plots.cxx.
void ana::DataMCAreaNormalizedRatio | ( | const Spectrum & | data, |
const Spectrum & | mc, | ||
double | miny, | ||
double | maxy | ||
) |
Plot data/MC ratio for the given spectrum. Normalize MC to Data by area.
Definition at line 337 of file Plots.cxx.
TH1 * ana::DataMCComparison | ( | const Spectrum & | data, |
const Spectrum & | mc | ||
) |
TH1 * ana::DataMCComparisonAreaNormalized | ( | const Spectrum & | data, |
const Spectrum & | mc | ||
) |
TH1 * ana::DataMCComparisonComponents | ( | const Spectrum & | data, |
const IPrediction * | mc, | ||
osc::IOscCalc * | calc | ||
) |
Plot MC broken down into flavour components, overlayed with data.
Definition at line 94 of file Plots.cxx.
void ana::DataMCRatio | ( | const Spectrum & | data, |
const IPrediction * | mc, | ||
osc::IOscCalc * | calc, | ||
double | miny, | ||
double | maxy | ||
) |
Plot data/MC ratio for the given spectrum. Normalize MC to Data by POT.
Definition at line 299 of file Plots.cxx.
void ana::DataMCRatio | ( | const Spectrum & | data, |
const Spectrum & | mc, | ||
double | miny, | ||
double | maxy | ||
) |
Plot data/MC ratio for the given spectrum. Normalize MC to Data by POT.
Definition at line 308 of file Plots.cxx.
const ReactorExperiment * ana::DayaBayConstraint2014 | ( | ) |
A ReactorExperiment initialized with the Nu2014 Daya Bay constraints.
Definition at line 56 of file ReactorExperiment.cxx.
osc::IOscCalcAdjustable * ana::DefaultOscCalc | ( | ) |
Create a new calc with default assumptions for all parameters.
Definition at line 39 of file Calcs.cxx.
osc::IOscCalcAdjustable * ana::DefaultOscCalcIH | ( | ) |
Definition at line 58 of file Calcs.cxx.
OscCalcSterileApproxAdjustable * ana::DefaultSterileApproxCalc | ( | SterileOscAngles | angles | ) |
Definition at line 293 of file OscCalcSterileApprox.cxx.
osc::OscCalcSterile * ana::DefaultSterileCalc | ( | int | nflavors | ) |
Create a sterile calc with default assumptions for all parameters.
Definition at line 87 of file Calcs.cxx.
const OscCalcSterileApprox * ana::DowncastToSterileApprox | ( | const osc::IOscCalc * | calc, |
bool | allowFail | ||
) |
Definition at line 309 of file OscCalcSterileApprox.cxx.
OscCalcSterileApprox * ana::DowncastToSterileApprox | ( | osc::IOscCalc * | calc, |
bool | allowFail | ||
) |
Definition at line 326 of file OscCalcSterileApprox.cxx.
void ana::DrawComponentsLegend | ( | TH1 * | hnue, |
TH1 * | hnumu, | ||
TH1 * | hnc, | ||
TH1 * | hcos, | ||
TH1 * | hother | ||
) |
Definition at line 146 of file tools_nuesel_icarus.h.
void ana::DrawEffPurLegend | ( | TH1D * | g1, |
char * | name1, | ||
TH1D * | g2, | ||
char * | name2 | ||
) |
Definition at line 230 of file tools_nuesel_icarus.h.
void ana::DrawErrorBand | ( | TH1 * | nom, |
TGraphAsymmErrors * | band, | ||
int | bandCol, | ||
double | alpha | ||
) |
Definition at line 196 of file EnsembleSpectrum.cxx.
void ana::DrawIntEffPurLegend | ( | TH1D * | g1, |
char * | name1, | ||
TH1D * | g2, | ||
char * | name2 | ||
) |
Definition at line 199 of file tools_nuesel_icarus.h.
void ana::DrawPurLegend | ( | TH1D * | g1, |
char * | name1 | ||
) |
Definition at line 217 of file tools_nuesel_icarus.h.
void ana::DrawSigBkgIntLegend | ( | TH1 * | h1, |
char * | name1, | ||
double | iSig, | ||
TH1 * | h2, | ||
char * | name2, | ||
double | iBkg | ||
) |
Definition at line 180 of file tools_nuesel_icarus.h.
void ana::DrawSigBkgIntText | ( | TH1 * | hsig, |
TH1 * | hbkg, | ||
float | textsize | ||
) |
Definition at line 243 of file tools_nuesel_icarus.h.
void ana::DrawSigBkgLegend | ( | TH1 * | h1, |
char * | name1, | ||
TH1 * | h2, | ||
char * | name2 | ||
) |
Definition at line 166 of file tools_nuesel_icarus.h.
TGraph* ana::EffOrPurGraph | ( | TH1 * | hSelSignal, |
TH1 * | hSelBack, | ||
TH1 * | hSignal, | ||
bool | geteff | ||
) |
Definition at line 292 of file tools_nuesel_icarus.h.
TH1D* ana::EffOrPurHistogram | ( | TH1 * | hSelSignal, |
TH1 * | hSelBack, | ||
TH1 * | hSignal, | ||
bool | geteff | ||
) |
Definition at line 271 of file tools_nuesel_icarus.h.
void ana::EnsurePositiveDefinite | ( | TH2 * | mat | ) |
Definition at line 903 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
TH2F* ana::ExpandedHistogram | ( | const std::string & | title, |
int | nbinsx, | ||
double | xmin, | ||
double | xmax, | ||
bool | xlog, | ||
int | nbinsy, | ||
double | ymin, | ||
double | ymax, | ||
bool | ylog | ||
) |
Definition at line 243 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
std::string ana::Experiment | ( | ) |
$EXPERIMENT or a nice error message and abort
Definition at line 92 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
ana::fBins | ( | {binsX, binsY} | ) |
ana::fCachedHash | ( | 0 | ) |
Definition at line 42 of file OscillatableSpectrum.cxx.
void ana::FillWithDimColor | ( | TH1 * | h, |
bool | usealpha = false , |
||
float | dim = 0.8 |
||
) |
Definition at line 130 of file tools_nuesel_icarus.h.
std::string ana::FindCAFAnaDir | ( | ) |
Definition at line 562 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
std::vector< double > ana::FindCurveCrossings | ( | TH1 * | h, |
double | critVal | ||
) |
double ana::FindOptimumCut | ( | TH1 * | hsig, |
TH1 * | hbkg, | ||
double & | best_fom | ||
) |
Search for optimum cut position given signal and background histograms.
hsig | Signal distribution | |
hbkg | Background distribution | |
[out] | best_fom | Figure of Merit obtained by best cut |
Definition at line 37 of file CutOptimizer.cxx.
double ana::FindQuantile | ( | double | frac, |
std::vector< double > & | xs | ||
) |
Definition at line 974 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
TGraph * ana::FindValley | ( | const IExperiment * | expt, |
osc::IOscCalcAdjustable * | calc, | ||
const IFitVar & | scanVar, | ||
const IFitVar & | fitVar, | ||
int | nbinsx, | ||
double | xmin, | ||
double | xmax, | ||
const std::vector< const IFitVar * > & | profVars = {} , |
||
const std::vector< const ISyst * > & | profSysts = {} , |
||
const std::map< const IFitVar *, std::vector< double >> & | seedPts = {} , |
||
const std::vector< SystShifts > & | systsSeedPts = {} , |
||
bool | transpose = false |
||
) |
Find the minimum in one variable as a function of another.
transpose | plot scanVar on the y axis |
Definition at line 657 of file Fit.cxx.
TString ana::fixLatexName | ( | TString | mystring | ) |
Definition at line 35 of file tools_nuesel_icarus.h.
TH2* ana::Flat | ( | double | level, |
const Surface & | s | ||
) |
Helper function for the gaussian approximation surfaces.
Definition at line 525 of file sbnana/sbnana/CAFAna/Analysis/Surface.cxx.
ana::fSplitBySign | ( | mode | = = kSplitBySign | ) |
Definition at line 37 of file PredictionInterp.cxx.
ana::fVars | ( | {varX, varY} | ) |
Definition at line 39 of file HistAxis.cxx.
TH2 * ana::Gaussian2Sigma1D | ( | const Surface & | s | ) |
Up-value surface for 2 sigma confidence in 1D in gaussian approximation.
Definition at line 548 of file sbnana/sbnana/CAFAna/Analysis/Surface.cxx.
TH2 * ana::Gaussian2Sigma2D | ( | const Surface & | s | ) |
Up-value surface for 2 sigma confidence in 2D in gaussian approximation.
Definition at line 540 of file sbnana/sbnana/CAFAna/Analysis/Surface.cxx.
TH2 * ana::Gaussian3Sigma1D | ( | const Surface & | s | ) |
Up-value surface for 3 sigma confidence in 1D in gaussian approximation.
Definition at line 550 of file sbnana/sbnana/CAFAna/Analysis/Surface.cxx.
TH2 * ana::Gaussian3Sigma1D1Sided | ( | const Surface & | s | ) |
Up-value surface for 3 sigma confidence in 1D in 1-sided gaussian approximation.
Definition at line 555 of file sbnana/sbnana/CAFAna/Analysis/Surface.cxx.
TH2 * ana::Gaussian3Sigma2D | ( | const Surface & | s | ) |
Up-value surface for 3 sigma confidence in 2D in gaussian approximation.
Definition at line 542 of file sbnana/sbnana/CAFAna/Analysis/Surface.cxx.
TH2 * ana::Gaussian5Sigma1D1Sided | ( | const Surface & | s | ) |
Up-value surface for 5 sigma confidence in 1D in 1-sided gaussian approximation.
Definition at line 556 of file sbnana/sbnana/CAFAna/Analysis/Surface.cxx.
TH2 * ana::Gaussian5Sigma2D | ( | const Surface & | s | ) |
Up-value surface for 5 sigma confidence in 2D in gaussian approximation.
Definition at line 543 of file sbnana/sbnana/CAFAna/Analysis/Surface.cxx.
TH2 * ana::Gaussian68Percent1D | ( | const Surface & | s | ) |
Up-value surface for 68% confidence in 1D in gaussian approximation.
Definition at line 545 of file sbnana/sbnana/CAFAna/Analysis/Surface.cxx.
TH2 * ana::Gaussian68Percent2D | ( | const Surface & | s | ) |
Up-value surface for 68% confidence in 2D in gaussian approximation.
Definition at line 537 of file sbnana/sbnana/CAFAna/Analysis/Surface.cxx.
TH2 * ana::Gaussian90Percent1D | ( | const Surface & | s | ) |
Up-value surface for 90% confidence in 1D in gaussian approximation.
Definition at line 546 of file sbnana/sbnana/CAFAna/Analysis/Surface.cxx.
TH2 * ana::Gaussian90Percent1D1Sided | ( | const Surface & | s | ) |
Up-value surface for 90% confidence in 1D in 1-sided gaussian approximation.
Definition at line 552 of file sbnana/sbnana/CAFAna/Analysis/Surface.cxx.
TH2 * ana::Gaussian90Percent2D | ( | const Surface & | s | ) |
Up-value surface for 90% confidence in 2D in gaussian approximation.
Definition at line 538 of file sbnana/sbnana/CAFAna/Analysis/Surface.cxx.
TH2 * ana::Gaussian95Percent1D | ( | const Surface & | s | ) |
Up-value surface for 95% confidence in 1D in gaussian approximation.
Definition at line 547 of file sbnana/sbnana/CAFAna/Analysis/Surface.cxx.
TH2 * ana::Gaussian95Percent1D1Sided | ( | const Surface & | s | ) |
Up-value surface for 95% confidence in 1D in 1-sided gaussian approximation.
Definition at line 553 of file sbnana/sbnana/CAFAna/Analysis/Surface.cxx.
TH2 * ana::Gaussian95Percent2D | ( | const Surface & | s | ) |
Up-value surface for 95% confidence in 2D in gaussian approximation.
Definition at line 539 of file sbnana/sbnana/CAFAna/Analysis/Surface.cxx.
TH2 * ana::Gaussian99Percent1D | ( | const Surface & | s | ) |
Up-value surface for 99% confidence in 1D in gaussian approximation.
Definition at line 549 of file sbnana/sbnana/CAFAna/Analysis/Surface.cxx.
TH2 * ana::Gaussian99Percent1D1Sided | ( | const Surface & | s | ) |
Up-value surface for 99% confidence in 1D in 1-sided gaussian approxiamtion.
Definition at line 554 of file sbnana/sbnana/CAFAna/Analysis/Surface.cxx.
TH2 * ana::Gaussian99Percent2D | ( | const Surface & | s | ) |
Up-value surface for 99% confidence in 2D in gaussian approximation.
Definition at line 541 of file sbnana/sbnana/CAFAna/Analysis/Surface.cxx.
void ana::gdb_backtrace | ( | ) |
Definition at line 16 of file DebugHelpers.cxx.
std::vector< const ISyst * > ana::GetAllNuMIFluxSysts | ( | unsigned int | Npcs | ) |
Combination of all beamline systs plus Npcs hadron production components.
Definition at line 160 of file NuMIFluxSysts.cxx.
std::vector< const ISyst * > ana::GetBigEnergySysts | ( | ) |
Definition at line 185 of file EnergySysts.cxx.
const BoosterFluxHadronSyst * ana::GetBoosterFluxHadronSyst | ( | unsigned int | i | ) |
Definition at line 63 of file BoosterFluxSysts.cxx.
BoosterFluxHadronSystVector ana::GetBoosterFluxHadronSysts | ( | unsigned int | N | ) |
N | total number of systematics |
Definition at line 76 of file BoosterFluxSysts.cxx.
std::map<std::string, std::string> ana::GetCAFMetadata | ( | TDirectory * | dir | ) |
Definition at line 587 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
std::vector< const ISyst * > ana::GetEnergySysts | ( | ) |
Definition at line 128 of file EnergySysts.cxx.
std::pair< double, double > ana::GetGenieDialLimits | ( | const std::string & | name | ) |
Definition at line 33 of file GenieWeightList.cxx.
int ana::GetGenieIndex | ( | const std::string & | name, |
bool | quiet | ||
) |
Definition at line 121 of file GenieWeightList.cxx.
double ana::GetGenieMax | ( | int | index | ) |
Definition at line 148 of file GenieWeightList.cxx.
double ana::GetGenieMin | ( | int | index | ) |
Definition at line 143 of file GenieWeightList.cxx.
std::string ana::GetGenieWeightName | ( | int | index | ) |
Definition at line 138 of file GenieWeightList.cxx.
std::vector< std::string > ana::GetGenieWeightNames | ( | ) |
Definition at line 9 of file GenieWeightList.cxx.
float ana::GetHistMax | ( | std::vector< TH1 * > | histos | ) |
Definition at line 108 of file tools_nuesel_icarus.h.
TH1* ana::GetMaskHist | ( | const Spectrum & | s, |
double | xmin, | ||
double | xmax, | ||
double | ymin, | ||
double | ymax | ||
) |
Definition at line 936 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
std::vector< TH1 * > ana::GetMCComponents | ( | const IPrediction * | mc, |
osc::IOscCalc * | calc, | ||
std::string | hist_name, | ||
double | pot, | ||
bool | force1D = false |
||
) |
A vector of histograms for the MC components. Easier to manipulate elsewhere Not ideal as returned pointers probably won't get deleted... but very useful...
Definition at line 167 of file Plots.cxx.
std::vector< TH1 * > ana::GetMCSystComponents | ( | const IPrediction * | mc, |
osc::IOscCalc * | calc, | ||
const SystShifts & | shift, | ||
std::string | hist_name, | ||
double | pot, | ||
bool | force1D | ||
) |
Definition at line 176 of file Plots.cxx.
TH1 * ana::GetMCSystTotal | ( | const IPrediction * | mc, |
osc::IOscCalc * | calc, | ||
const SystShifts & | shift, | ||
std::string | hist_name, | ||
double | pot, | ||
bool | force1D | ||
) |
TH1 * ana::GetMCTotal | ( | const IPrediction * | mc, |
osc::IOscCalc * | calc, | ||
std::string | hist_name, | ||
double | pot, | ||
bool | force1D | ||
) |
Definition at line 158 of file Plots.cxx.
std::vector< TH1 * > ana::GetMCTotalForSystShifts | ( | const IPrediction * | mc, |
osc::IOscCalc * | calc, | ||
const ISyst * | syst, | ||
std::string | hist_base_name, | ||
double | pot, | ||
bool | force1D | ||
) |
Definition at line 278 of file Plots.cxx.
std::vector< const ISyst * > ana::GetNuMIBeamlineFluxSysts | ( | ) |
Definition at line 127 of file NuMIFluxSysts.cxx.
const NuMIFluxSyst * ana::GetNuMIFluxSyst | ( | const std::string & | dir, |
const std::string & | prefix, | ||
const std::string & | name | ||
) |
Definition at line 91 of file NuMIFluxSysts.cxx.
std::vector< const ISyst * > ana::GetNuMIHadronProductionFluxSysts | ( | ) |
These are envelopes not real systs. TODO make clearer in naming.
Definition at line 105 of file NuMIFluxSysts.cxx.
std::vector< const ISyst * > ana::GetNuMIPCAFluxSysts | ( | unsigned int | Npcs = 100 | ) |
Npcs | Number of principal components. High-indexed components should have negligible effect on the analysis. |
Definition at line 149 of file NuMIFluxSysts.cxx.
std::vector<std::string> ana::GetSBNBoosterFluxWeightNames | ( | ) |
Definition at line 172 of file SBNWeightSysts.cxx.
const std::vector< const ISyst * > & ana::GetSBNBoosterFluxWeightSysts | ( | ) |
Definition at line 188 of file SBNWeightSysts.cxx.
std::vector< std::string > ana::GetSBNGenieWeightNames | ( | ) |
Definition at line 113 of file SBNWeightSysts.cxx.
const std::vector< const ISyst * > & ana::GetSBNGenieWeightSysts | ( | ) |
Definition at line 159 of file SBNWeightSysts.cxx.
const std::vector< const ISyst * > & ana::GetSBNWeightSysts | ( | ) |
Definition at line 200 of file SBNWeightSysts.cxx.
caf::SRGlobal ana::GetSRGlobal | ( | ) |
Definition at line 22 of file UniverseOracle.cxx.
Var ana::GetUniverseWeight | ( | const std::string & | psetName, |
int | univIdx | ||
) |
Definition at line 27 of file SBNWeightSysts.h.
TGraphAsymmErrors * ana::GraphWithPoissonErrors | ( | const TH1 * | h, |
bool | noErrorsXaxis, | ||
bool | drawEmptyBins | ||
) |
void ana::handle_signal | ( | int | sig, |
siginfo_t * | , | ||
void * | |||
) |
Definition at line 34 of file DebugHelpers.cxx.
void ana::handle_terminate | ( | ) |
Definition at line 27 of file DebugHelpers.cxx.
bool ana::HasBraggPeak | ( | const caf::SRTrueParticleProxy & | p | ) |
Whether this particle should have a bragg peak in the detector.
Definition at line 11 of file TrueParticleHelpers.cxx.
bool ana::HasVar | ( | std::vector< const IFitVar * > | oscVars, |
std::string | name | ||
) |
Definition at line 39 of file CalcsNuFit.cxx.
|
static |
Definition at line 48 of file NumuVarsIcarus202208.cxx.
|
static |
Definition at line 9 of file NumuCutsIcarus202208.cxx.
bool ana::IsCCQEOnArgon | ( | const caf::SRTrueInteractionProxy * | nu, |
int | pdg | ||
) |
Definition at line 10 of file Flux.cxx.
Cut ana::IsCCQEOnArgonCut | ( | int | pdg | ) |
Definition at line 46 of file Flux.cxx.
bool ana::IsGenie | ( | const caf::SRTrueParticleProxy & | p | ) |
Whether this particle was generated by genie (as opposed to geant or corsika)
Definition at line 22 of file TrueParticleHelpers.cxx.
bool ana::IsNCQEOnArgon | ( | const caf::SRTrueInteractionProxy * | nu, |
int | pdg | ||
) |
Definition at line 22 of file Flux.cxx.
bool ana::IsPrimary | ( | const caf::SRTrueParticleProxy & | p | ) |
Whether this is a primary particle or generated by a secondary interaction.
Definition at line 7 of file TrueParticleHelpers.cxx.
bool ana::IsStable | ( | const caf::SRTrueParticleProxy & | p | ) |
Whether this is a stable particle as generated by genie.
Definition at line 26 of file TrueParticleHelpers.cxx.
size_t ana::JobNumber | ( | ) |
Definition at line 776 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
const Var ana::kBaseline | ( | [] (const caf::SRSliceProxy *slc) -> double{return slc->truth.baseline *1e-3;} | ) |
const Cut ana::kIsAntiNu | ( | [] (const caf::SRSliceProxy *slc){if(slc->truth.index< 0) return false;return slc->truth.pdg< 0;} | ) |
const Cut ana::kIsBeamNue | ( | CCFlavSel(12, 12) | ) |
Select CC .
const Cut ana::kIsCC | ( | [] (const caf::SRSliceProxy *slc){if(slc->truth.index< 0) return false;return(slc->truth.iscc==1);} | ) |
const Cut ana::kIsNC | ( | [] (const caf::SRSliceProxy *slc){if(slc->truth.index< 0) return false;return(slc->truth.isnc==1);} | ) |
const Cut ana::kIsNCFromNue | ( | NCFlavOrig(12) | ) |
const Cut ana::kIsNCFromNumu | ( | NCFlavOrig(14) | ) |
const Cut ana::kIsNueApp | ( | CCFlavSel(12, 14) | ) |
Select CC .
const Cut ana::kIsNumuApp | ( | CCFlavSel(14, 12) | ) |
Select CC .
const Cut ana::kIsNumuCC | ( | CCFlavSel(14, 14) | ) |
Select CC .
const Cut ana::kIsTauFromE | ( | CCFlavSel(16, 12) | ) |
Select CC .
const Cut ana::kIsTauFromMu | ( | CCFlavSel(16, 14) | ) |
Select CC .
const Cut ana::kNoCut | ( | [] (const caf::SRSliceProxy *){return true;} | ) |
The simplest possible cut: pass everything, used as a default.
const SpillCut ana::kNoSpillCut | ( | [] (const caf::SRSpillProxy *){return true;} | ) |
The simplest possible cut: pass everything, used as a default.
const Var ana::kNuScore | ( | [] (const caf::SRSliceProxy *slc) -> double{return slc->nu_score;} | ) |
const Var ana::kSlcIsRecoNu | ( | [] (const caf::SRSliceProxy *slc) -> double{return!slc->is_clear_cosmic;} | ) |
const SolarConstraints ana::kSolarConstraintsPDG2017 | ( | 7.53e- | 5, |
0.18e- | 5, | ||
0. | 851, | ||
0. | 020 | ||
) |
const SpillVar ana::kSpillUnweighted | ( | [] (const caf::SRSpillProxy *){return 1;} | ) |
const Var ana::kUnweighted | ( | [] (const caf::SRSliceProxy *){return 1;} | ) |
The simplest possible Var, always 1. Used as a default weight.
const caf::SRShowerProxy* ana::LargestRecoShower | ( | const caf::SRSliceProxy * | slc | ) |
Pointer to largest reconstructed shower, or null pointer if none exists.
Definition at line 26 of file NueVars.cxx.
int ana::Limit | ( | ) |
Definition at line 758 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
std::vector<std::string> ana::LoadFileList | ( | const std::string & | listfile | ) |
Definition at line 568 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
std::unique_ptr<T> ana::LoadFrom | ( | TDirectory * | dir | ) |
Definition at line 22 of file LoadFromFile.h.
std::unique_ptr<IBkgdEstimator> ana::LoadFrom< IBkgdEstimator > | ( | TDirectory * | dir | ) |
std::unique_ptr<IDecomp> ana::LoadFrom< IDecomp > | ( | TDirectory * | dir | ) |
std::unique_ptr< IExperiment > ana::LoadFrom< IExperiment > | ( | TDirectory * | dir | ) |
Definition at line 22 of file IExperiment.cxx.
std::unique_ptr< IExtrap > ana::LoadFrom< IExtrap > | ( | TDirectory * | dir | ) |
Definition at line 17 of file IExtrap.cxx.
std::unique_ptr< IPrediction > ana::LoadFrom< IPrediction > | ( | TDirectory * | dir | ) |
Definition at line 26 of file IPrediction.cxx.
std::unique_ptr<ModularExtrapComponent> ana::LoadFrom< ModularExtrapComponent > | ( | TDirectory * | dir | ) |
std::unique_ptr< osc::IOscCalc > ana::LoadFrom< osc::IOscCalc > | ( | TDirectory * | dir | ) |
Definition at line 22 of file LoadFromFile.cxx.
std::unique_ptr<T> ana::LoadFromFile | ( | const std::string & | fname, |
const std::string & | label | ||
) |
Definition at line 61 of file LoadFromFile.h.
double ana::LogLikelihood | ( | double | exp, |
double | obs | ||
) |
The log-likelihood formula for a single bin.
exp | Expected count |
obs | Observed count |
Handles zero observed or expected events correctly.
Definition at line 163 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
double ana::LogLikelihood | ( | const TH1 * | exp, |
const TH1 * | obs, | ||
bool | useOverflow = false |
||
) |
The log-likelihood formula from the PDG.
exp | The expected spectrum |
obs | The corresponding observed spectrum |
Includes underflow bin and an option for overflow bin (off by default) and handles zero observed or expected events correctly.
Definition at line 202 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
const caf::SRTrackProxy* ana::LongestRecoTrack | ( | const caf::SRSliceProxy * | slc | ) |
Pointer to longest reconstructed shower, or null pointer if none exists.
Definition at line 173 of file NueVars.cxx.
const caf::SRTrkChi2PIDProxy* ana::LongestTrackBestPlaneChi2PID | ( | const caf::SRSliceProxy * | slc | ) |
Definition at line 199 of file NueVars.cxx.
void ana::MakeNMinusOneSpectra | ( | SpectrumLoader & | loader, |
const Cut & | sigcut, | ||
const Cut & | presel, | ||
const std::vector< HistAxis > & | axes, | ||
const std::vector< double > & | cut_pos, | ||
std::vector< Spectrum > & | sigs, | ||
std::vector< Spectrum > & | bkgs | ||
) |
Make a series of spectra leaving out one cut in turn.
This function does not call Go() on the loader, allowing you to combine it with other spectra-filling if desired.
loader | SpectrumLoader to associate spectra with | |
sigcut | Definition of signal (inverse defines background) | |
presel | Cuts that will always be applied | |
axes | Variables that will be cut on | |
cut_pos | Cut position for each variable | |
[out] | sigs | Resulting signal spectra |
[out] | bkgs | Resulting background spectra |
Definition at line 10 of file CutOptimizer.cxx.
TH1D* ana::MakeTH1D | ( | const char * | name, |
const char * | title, | ||
const Binning & | bins | ||
) |
Definition at line 393 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
TH2D* ana::MakeTH2D | ( | const char * | name, |
const char * | title, | ||
const Binning & | binsx, | ||
const Binning & | binsy | ||
) |
Definition at line 399 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
_MultiVar< T > ana::MultiVar2D | ( | const _MultiVar< T > & | a, |
const Binning & | binsa, | ||
const _MultiVar< T > & | b, | ||
const Binning & | binsb | ||
) |
Definition at line 69 of file MultiVar.cxx.
template SpillMultiVar ana::MultiVar2D | ( | const MultiVar & | , |
const Binning & | , | ||
const MultiVar & | , | ||
const Binning & | |||
) |
osc::IOscCalcAdjustable * ana::NuFitOscCalc | ( | int | hie | ) |
Definition at line 12 of file CalcsNuFit.cxx.
osc::IOscCalcAdjustable * ana::NuFitOscCalcCDR | ( | int | hie | ) |
Definition at line 10 of file CalcsNuFit_cdr.cxx.
osc::IOscCalcAdjustable * ana::NuFitOscCalcCDRPlusOneSigma | ( | int | hie | ) |
Definition at line 73 of file CalcsNuFit_cdr.cxx.
osc::IOscCalcAdjustable * ana::NuFitOscCalcPlusOneSigma | ( | int | hie | ) |
Definition at line 98 of file CalcsNuFit.cxx.
size_t ana::NumJobs | ( | ) |
Definition at line 787 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
size_t ana::Offset | ( | bool | allow_default | ) |
Definition at line 734 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
_Cut< T > ana::operator! | ( | const _Cut< T > & | a | ) |
template Cut ana::operator!< caf::SRSliceProxy > | ( | const Cut & | a | ) |
template SpillCut ana::operator!< caf::SRSpillProxy > | ( | const SpillCut & | a | ) |
_Cut< T > ana::operator!= | ( | const _Var< T > & | v, |
double | c | ||
) |
_Cut< T > ana::operator!= | ( | double | c, |
const _Var< T > & | v | ||
) |
_Cut< T > ana::operator!= | ( | const _Var< T > & | a, |
const _Var< T > & | b | ||
) |
template Cut ana::operator!= | ( | const Var & | , |
double | |||
) |
template Cut ana::operator!= | ( | const Var & | , |
const Var & | |||
) |
|
inline |
Definition at line 97 of file OscCalcSterileApprox.h.
_Cut< T > ana::operator&& | ( | const _Cut< T > & | a, |
const _Cut< T > & | b | ||
) |
Definition at line 49 of file Cut.cxx.
template Cut ana::operator&&< caf::SRSliceProxy > | ( | const Cut & | a, |
const Cut & | b | ||
) |
template SpillCut ana::operator&&< caf::SRSpillProxy > | ( | const SpillCut & | a, |
const SpillCut & | b | ||
) |
|
inline |
Definition at line 79 of file EnsembleSpectrum.h.
_Var< T > ana::operator* | ( | const _Var< T > & | a, |
const _Var< T > & | b | ||
) |
template Var ana::operator* | ( | const Var & | , |
const Var & | |||
) |
|
inline |
Definition at line 353 of file Spectrum.h.
_Var< T > ana::operator+ | ( | const _Var< T > & | a, |
const _Var< T > & | b | ||
) |
template Var ana::operator+ | ( | const Var & | , |
const Var & | |||
) |
_Var< T > ana::operator- | ( | const _Var< T > & | a, |
const _Var< T > & | b | ||
) |
template Var ana::operator- | ( | const Var & | , |
const Var & | |||
) |
|
inline |
|
inline |
Definition at line 39 of file EnsembleRatio.h.
|
inline |
Definition at line 80 of file EnsembleSpectrum.h.
_Var< T > ana::operator/ | ( | const _Var< T > & | a, |
const _Var< T > & | b | ||
) |
template Var ana::operator/ | ( | const Var & | , |
const Var & | |||
) |
|
inline |
Definition at line 354 of file Spectrum.h.
_Cut< T > ana::operator< | ( | const _Var< T > & | v, |
double | c | ||
) |
_Cut< T > ana::operator< | ( | double | c, |
const _Var< T > & | v | ||
) |
_Cut< T > ana::operator< | ( | const _Var< T > & | a, |
const _Var< T > & | b | ||
) |
template Cut ana::operator< | ( | const Var & | , |
double | |||
) |
template Cut ana::operator< | ( | const Var & | , |
const Var & | |||
) |
template Cut ana::operator< | ( | double | , |
const Var & | |||
) |
_Cut< T > ana::operator<= | ( | const _Var< T > & | v, |
double | c | ||
) |
_Cut< T > ana::operator<= | ( | double | c, |
const _Var< T > & | v | ||
) |
_Cut< T > ana::operator<= | ( | const _Var< T > & | a, |
const _Var< T > & | b | ||
) |
template Cut ana::operator<= | ( | const Var & | , |
double | |||
) |
template Cut ana::operator<= | ( | const Var & | , |
const Var & | |||
) |
template Cut ana::operator<= | ( | double | , |
const Var & | |||
) |
_Cut< T > ana::operator== | ( | const _Var< T > & | v, |
double | c | ||
) |
_Cut< T > ana::operator== | ( | const _Var< T > & | a, |
const _Var< T > & | b | ||
) |
template Cut ana::operator== | ( | const Var & | , |
double | |||
) |
template Cut ana::operator== | ( | const Var & | , |
const Var & | |||
) |
_Cut< T > ana::operator> | ( | const _Var< T > & | v, |
double | c | ||
) |
_Cut< T > ana::operator> | ( | const _Var< T > & | a, |
const _Var< T > & | b | ||
) |
_Cut< T > ana::operator> | ( | double | c, |
const _Var< T > & | v | ||
) |
template Cut ana::operator> | ( | const Var & | , |
double | |||
) |
template Cut ana::operator> | ( | const Var & | , |
const Var & | |||
) |
template Cut ana::operator> | ( | double | , |
const Var & | |||
) |
_Cut< T > ana::operator>= | ( | const _Var< T > & | v, |
double | c | ||
) |
_Cut< T > ana::operator>= | ( | const _Var< T > & | a, |
const _Var< T > & | b | ||
) |
_Cut< T > ana::operator>= | ( | double | c, |
const _Var< T > & | v | ||
) |
template Cut ana::operator>= | ( | const Var & | , |
double | |||
) |
template Cut ana::operator>= | ( | const Var & | , |
const Var & | |||
) |
template Cut ana::operator>= | ( | double | , |
const Var & | |||
) |
|
inline |
Definition at line 90 of file OscCalcSterileApprox.h.
|
inline |
_Cut< T > ana::operator|| | ( | const _Cut< T > & | a, |
const _Cut< T > & | b | ||
) |
Definition at line 76 of file Cut.cxx.
template Cut ana::operator||< caf::SRSliceProxy > | ( | const Cut & | a, |
const Cut & | b | ||
) |
template SpillCut ana::operator||< caf::SRSpillProxy > | ( | const SpillCut & | a, |
const SpillCut & | b | ||
) |
void ana::OptimizeCuts | ( | const std::string & | wildcard, |
double | pot, | ||
const Cut & | sigcut, | ||
const Cut & | presel, | ||
const std::vector< HistAxis > & | axes, | ||
std::vector< double > & | cut_pos | ||
) |
Repeatedly invoke OptimizeOneCut until the FOM increase becomes small.
Definition at line 99 of file CutOptimizer.cxx.
double ana::OptimizeOneCut | ( | const std::string & | wildcard, |
double | pot, | ||
const Cut & | sigcut, | ||
const Cut & | presel, | ||
const std::vector< HistAxis > & | axes, | ||
std::vector< double > & | cut_pos | ||
) |
Scan all cuts and update the one giving the largest FOM gain.
wildcard | File name / wildcard / dataset name | |
pot | POT to scale to | |
sigcut | Definition of signal (inverse defines background) | |
presel | Cuts that will always be applied | |
axes | Variables that will be cut on | |
[out] | cut_pos | Starting cut position for each variable. One of these values will be updated |
Definition at line 58 of file CutOptimizer.cxx.
void ana::PimpHist | ( | TH1 * | histo, |
Color_t | color, | ||
Style_t | linestyle, | ||
int | linewidth, | ||
Style_t | markerstyle = 8 , |
||
double | markersize = 8 |
||
) |
Definition at line 119 of file tools_nuesel_icarus.h.
void ana::PlotWithSystErrorBand | ( | IPrediction * | pred, |
const std::vector< const ISyst * > & | systs, | ||
osc::IOscCalc * | calc, | ||
double | pot, | ||
int | col = -1 , |
||
int | errCol = -1 , |
||
float | headroom = 1.3 , |
||
bool | newaxis = true |
||
) |
Plot prediction with +/-1sigma error band.
When multiple systematics are used, the errors are the quadrature sum
pred | The prediction. Must be able to generate shifted spectra |
systs | List of systematics to include in error band |
calc | Oscillation parameters |
pot | POT to evaluate prediction at |
col | Color of the prediction, default kRed |
errCol | Color of the shading, default col-7 (kRed-7 is light red) |
headroom | Fraction of maximum bin for headroom, default 30% |
Definition at line 377 of file Plots.cxx.
void ana::PlotWithSystErrorBand | ( | const Spectrum & | nominal, |
const std::vector< Spectrum > & | upShifts, | ||
const std::vector< Spectrum > & | downShifts, | ||
double | pot, | ||
int | col = -1 , |
||
int | errCol = -1 , |
||
float | headroom = 1.3 , |
||
bool | newaxis = true |
||
) |
Plot prediction with error band.
When multiple systematics are used, the errors are the quadrature sum
nominal | Nominal spectrum |
upShifts | Vector of spectra which have + shifts |
downShifts | Vector of spectra which have - shifts, same order as + |
pot | POT to scale spectra to |
col | Color of the prediction, default kRed |
errCol | Color of the shading, default col-7(kRed-7 is light red) |
headroom | Fraction of maximum bin for headroom, default 30% |
Definition at line 404 of file Plots.cxx.
std::string ana::pnfs2xrootd | ( | std::string | loc, |
bool | unauth | ||
) |
Definition at line 808 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
double ana::PointDistanceToBox | ( | double | x, |
double | y, | ||
double | x0, | ||
double | y0, | ||
double | x1, | ||
double | y1 | ||
) |
void ana::printEventsLine | ( | std::string | cutname, |
float | nue, | ||
float | numu, | ||
float | nc, | ||
float | cos, | ||
float | other, | ||
float | eff, | ||
float | pur | ||
) |
Definition at line 67 of file tools_nuesel_icarus.h.
void ana::PrintSRGlobal | ( | const caf::SRGlobal & | global | ) |
Definition at line 54 of file UniverseOracle.cxx.
void ana::printTableFooter | ( | ) |
Definition at line 60 of file tools_nuesel_icarus.h.
void ana::printTableHeader | ( | int | quantId = 0 | ) |
Definition at line 45 of file tools_nuesel_icarus.h.
TH1 * ana::Profile | ( | const IExperiment * | expt, |
osc::IOscCalcAdjustable * | calc, | ||
const IFitVar * | v, | ||
int | nbinsx, | ||
double | minx, | ||
double | maxx, | ||
double | minchi = -1 , |
||
const std::vector< const IFitVar * > & | profVars = {} , |
||
const std::vector< const ISyst * > & | profSysts = {} , |
||
const std::map< const IFitVar *, std::vector< double >> & | seedPts = {} , |
||
const std::vector< SystShifts > & | systsSeedPts = {} , |
||
std::map< const IFitVar *, TGraph * > & | profVarsMap = empty_vars_map , |
||
std::map< const ISyst *, TGraph * > & | systsMap = empty_syst_map |
||
) |
scan in one variable, profiling over all others
expt | The experiment to retrieve chisq values from | |
calc | Initial values of all oscillation parameters | |
v | Scan over this variable | |
nbinsx | Binning | |
minx | Binning | |
maxx | Binning | |
minchi | Set non-default to force a chisq value to evaluate delta chisqs against. Useful for comparing two profiles. If set to zero it will not zero-adjust and the axis will be labelled without "delta" | |
profVars | Profile over these variables | |
profSysts | Profile over these systematics | |
seedPts | Set each var to each of the values. Try all combinations. Beware of combinatorical explosion... | |
systSeedPts | If non-empty, try fit starting at each of these | |
[out] | profVarsMap | Pass empty map. Returns best values of each var. |
[out] | systsMap | Pass empty map. Returns best values of each syst. |
Definition at line 390 of file Fit.cxx.
TH1 * ana::Profile | ( | const IExperiment * | expt, |
osc::IOscCalcAdjustable * | calc, | ||
const ISyst * | s, | ||
int | nbinsx, | ||
double | minx, | ||
double | maxx, | ||
double | input_minchi, | ||
const std::vector< const IFitVar * > & | profVars, | ||
const std::vector< const ISyst * > & | profSysts, | ||
const std::map< const IFitVar *, std::vector< double >> & | seedPts, | ||
const std::vector< SystShifts > & | systSeedPts, | ||
std::map< const IFitVar *, TGraph * > & | profVarsMap, | ||
std::map< const ISyst *, TGraph * > & | profSystsMap | ||
) |
Definition at line 477 of file Fit.cxx.
TGraphAsymmErrors * ana::ProfileQuantile | ( | const TH2 * | hist, |
const std::string & | axisName, | ||
const std::string & | graphName, | ||
const std::pair< double, double > & | quantileDivisions | ||
) |
Calculate profile with error bars corresponding to specified quantiles of a 2D distribution (by default, 68% coverage)
Definition at line 728 of file Plots.cxx.
void ana::ProjectionX | ( | TH2D * | from, |
TH1D * | to | ||
) |
void ana::ProjectionY | ( | TH2D * | from, |
TH1D * | to | ||
) |
void ana::RatioPlot | ( | const Spectrum & | data, |
const Spectrum & | expected, | ||
const Spectrum & | fit, | ||
double | miny, | ||
double | maxy | ||
) |
Plot data/expected, compared with fit/expected.
Definition at line 348 of file Plots.cxx.
void ana::ResetOscCalcToDefault | ( | osc::IOscCalcAdjustable * | calc | ) |
void ana::ResetOscCalcToDefaultIH | ( | osc::IOscCalcAdjustable * | calc | ) |
Definition at line 47 of file Calcs.cxx.
void ana::ResetSterileCalcToDefault | ( | osc::OscCalcSterile * | calc | ) |
Reset calc to default assumptions for all parameters.
Definition at line 66 of file Calcs.cxx.
bool ana::RunningOnGrid | ( | ) |
Definition at line 697 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
std::string ana::SAMExperiment | ( | ) |
$SAM_EXPERIMENT or a nice error message and abort
Definition at line 105 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
void ana::SaveTo | ( | const T & | x, |
TDirectory * | dir | ||
) |
Definition at line 52 of file LoadFromFile.h.
void ana::SaveTo | ( | const osc::IOscCalc & | x, |
TDirectory * | dir | ||
) |
Definition at line 91 of file LoadFromFile.cxx.
void ana::SaveToFile | ( | const T & | x, |
const std::string & | fname, | ||
const std::string & | label | ||
) |
Definition at line 75 of file LoadFromFile.h.
Var ana::Scaled | ( | const Var & | v, |
double | s | ||
) |
Use to rescale another variable.
Definition at line 159 of file Var.cxx.
TGraph * ana::ShadeBetweenHistograms | ( | TH1 * | hmin, |
TH1 * | hmax | ||
) |
double ana::Si | ( | double | x | ) |
Definition at line 138 of file OscCalcSterileApprox.cxx.
double ana::SimpleFOM | ( | const Spectrum & | obs, |
const Spectrum & | unosc, | ||
double | pot = 0 |
||
) |
Figure-of-merit with no systematics, for binned data.
obs | The observed data |
unosc | A spectrum representing the null hypothesis |
pot | POT to scale to. Leave at default to adopt POT from obs |
Definition at line 28 of file Fit.cxx.
TH1 * ana::Slice | ( | const IExperiment * | expt, |
osc::IOscCalcAdjustable * | calc, | ||
const IFitVar * | v, | ||
int | nbinsx, | ||
double | minx, | ||
double | maxx, | ||
double | minchi | ||
) |
scan in one variable, holding all others constant
Definition at line 611 of file Fit.cxx.
TH1 * ana::Slice | ( | const IExperiment * | expt, |
osc::IOscCalcAdjustable * | calc, | ||
const ISyst * | s, | ||
int | nbinsx, | ||
double | minx, | ||
double | maxx, | ||
double | minchi | ||
) |
Definition at line 620 of file Fit.cxx.
void ana::SplitCanvas2 | ( | TCanvas *& | c1, |
TPad *& | pad1, | ||
TPad *& | pad2 | ||
) |
Definition at line 83 of file tools_nuesel_icarus.h.
TH1 * ana::SqrtProfile | ( | const IExperiment * | expt, |
osc::IOscCalcAdjustable * | calc, | ||
const IFitVar * | v, | ||
int | nbinsx, | ||
double | minx, | ||
double | maxx, | ||
double | minchi, | ||
std::vector< const IFitVar * > | profVars, | ||
std::vector< const ISyst * > | profSysts, | ||
const std::map< const IFitVar *, std::vector< double >> & | seedPts, | ||
const std::vector< SystShifts > & | systSeedPts, | ||
std::map< const IFitVar *, TGraph * > & | profVarsMap, | ||
std::map< const ISyst *, TGraph * > & | systsMap | ||
) |
Forward to Profile but sqrt the result for a crude significance.
Definition at line 565 of file Fit.cxx.
TH1 * ana::SqrtProfile | ( | const IExperiment * | expt, |
osc::IOscCalcAdjustable * | calc, | ||
const ISyst * | s, | ||
int | nbinsx, | ||
double | minx, | ||
double | maxx, | ||
double | minchi, | ||
std::vector< const IFitVar * > | profVars, | ||
std::vector< const ISyst * > | profSysts, | ||
const std::map< const IFitVar *, std::vector< double >> & | seedPts, | ||
const std::vector< SystShifts > & | systSeedPts, | ||
std::map< const IFitVar *, TGraph * > & | profVarsMap, | ||
std::map< const ISyst *, TGraph * > & | systsMap | ||
) |
Definition at line 588 of file Fit.cxx.
TH1 * ana::SqrtSlice | ( | const IExperiment * | expt, |
osc::IOscCalcAdjustable * | calc, | ||
const IFitVar * | v, | ||
int | nbinsx, | ||
double | minx, | ||
double | maxx, | ||
double | minchi | ||
) |
TH1 * ana::SqrtSlice | ( | const IExperiment * | expt, |
osc::IOscCalcAdjustable * | calc, | ||
const ISyst * | s, | ||
int | nbinsx, | ||
double | minx, | ||
double | maxx, | ||
double | minchi | ||
) |
Definition at line 643 of file Fit.cxx.
size_t ana::Stride | ( | bool | allow_default | ) |
Definition at line 710 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
std::unique_ptr<TMatrixD> ana::SymmMxInverse | ( | const TMatrixD & | mx | ) |
TString ana::thisCellColor | ( | double | weird | ) |
Definition at line 26 of file tools_nuesel_icarus.h.
osc::IOscCalcAdjustable * ana::ThrownNuFitOscCalc | ( | int | hie, |
std::vector< const IFitVar * > | oscVars | ||
) |
Definition at line 45 of file CalcsNuFit.cxx.
osc::IOscCalcAdjustable * ana::ThrownNuFitOscCalcCDR | ( | int | hie | ) |
Definition at line 38 of file CalcsNuFit_cdr.cxx.
TH2* ana::ToTH2 | ( | const Spectrum & | s, |
double | exposure, | ||
ana::EExposureType | expotype, | ||
const Binning & | binsx, | ||
const Binning & | binsy, | ||
ana::EBinType | bintype | ||
) |
Definition at line 407 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
TH2* ana::ToTH2 | ( | const Ratio & | r, |
const Binning & | binsx, | ||
const Binning & | binsy | ||
) |
Definition at line 417 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
TH2* ana::ToTH2Helper | ( | std::unique_ptr< TH1 > | h1, |
const Binning & | binsx, | ||
const Binning & | binsy, | ||
ana::EBinType | bintype | ||
) |
Definition at line 427 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
TH3* ana::ToTH3 | ( | const Spectrum & | s, |
double | exposure, | ||
ana::EExposureType | expotype, | ||
const Binning & | binsx, | ||
const Binning & | binsy, | ||
const Binning & | binsz, | ||
ana::EBinType | bintype | ||
) |
Definition at line 454 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
TH3* ana::ToTH3 | ( | const Ratio & | r, |
const Binning & | binsx, | ||
const Binning & | binsy, | ||
const Binning & | binsz | ||
) |
Definition at line 467 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
TH3* ana::ToTH3Helper | ( | std::unique_ptr< TH1 > | h1, |
const Binning & | binsx, | ||
const Binning & | binsy, | ||
const Binning & | binsz, | ||
ana::EBinType | bintype | ||
) |
Definition at line 478 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
THStack * ana::ToTHStack | ( | const std::vector< std::pair< Spectrum, Color_t >> & | s, |
double | pot | ||
) |
Binning ana::TrueEnergyBins | ( | ) |
Default true-energy bin edges.
Definition at line 152 of file Binning.cxx.
Binning ana::TrueLOverTrueEBins | ( | ) |
Definition at line 185 of file Binning.cxx.
Spectrum ana::UnfoldIterative | ( | const Spectrum & | reco, |
const ReweightableSpectrum & | recoVsTrue, | ||
unsigned int | nIterations | ||
) |
Definition at line 8 of file UnfoldIterative.cxx.
Spectrum ana::UnfoldSVD | ( | const Spectrum & | reco, |
const ReweightableSpectrum & | recoVsTrue, | ||
unsigned int | reg | ||
) |
Singular value decomposition unfolding using ROOT's TSVDDecomp.
SVD decomposition requires the reco and true spectra to have the same number of bins (recoVsTrue must be square).
Lower values of reg are more regularized
Definition at line 12 of file UnfoldSVD.cxx.
Spectrum ana::UnfoldTikhonov | ( | const Spectrum & | reco, |
const ReweightableSpectrum & | recoVsTrue, | ||
double | regStrength | ||
) |
Unfolding using Tikhonov regularization (penalizing true spectra with large second derivatives)
Currently requires 1D reco and true axes, but that is not a fundamental restriction of the method and could be lifted in future versions.
Large regStrength corresponds to more regularization (biased towards smooth results), small strength suffers more from statistical fluctuations.
Definition at line 16 of file UnfoldTikhonov.cxx.
std::string ana::UniqueName | ( | ) |
Return a different string each time, for creating histograms.
Definition at line 32 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
_Var< T > ana::Var2D | ( | const _Var< T > & | a, |
const Binning & | binsa, | ||
const _Var< T > & | b, | ||
const Binning & | binsb | ||
) |
_Var< T > ana::Var2D | ( | const _Var< T > & | a, |
int | na, | ||
double | a0, | ||
double | a1, | ||
const _Var< T > & | b, | ||
int | nb, | ||
double | b0, | ||
double | b1 | ||
) |
Variable formed from two input variables.
The binning of each variable has to be given to allow conversion into a 1D binned representation.
Definition at line 117 of file Var.cxx.
template Var ana::Var2D | ( | const Var & | , |
const Binning & | , | ||
const Var & | , | ||
const Binning & | |||
) |
template Var ana::Var2D | ( | const Var & | , |
int | , | ||
double | , | ||
double | , | ||
const Var & | , | ||
int | , | ||
double | , | ||
double | |||
) |
_Var< T > ana::Var3D | ( | const _Var< T > & | a, |
const Binning & | binsa, | ||
const _Var< T > & | b, | ||
const Binning & | binsb, | ||
const _Var< T > & | c, | ||
const Binning & | binsc | ||
) |
_Var< T > ana::Var3D | ( | const _Var< T > & | a, |
int | na, | ||
double | a0, | ||
double | a1, | ||
const _Var< T > & | b, | ||
int | nb, | ||
double | b0, | ||
double | b1, | ||
const _Var< T > & | c, | ||
int | nc, | ||
double | c0, | ||
double | c1 | ||
) |
This is just like a Var2D, but useful for 3D Spectra.
Definition at line 142 of file Var.cxx.
template Var ana::Var3D | ( | const Var & | , |
const Binning & | , | ||
const Var & | , | ||
const Binning & | , | ||
const Var & | , | ||
const Binning & | |||
) |
template Var ana::Var3D | ( | const Var & | , |
int | , | ||
double | , | ||
double | , | ||
const Var & | , | ||
int | , | ||
double | , | ||
double | , | ||
const Var & | , | ||
int | , | ||
double | , | ||
double | |||
) |
std::vector<std::string> ana::Wildcard | ( | const std::string & | wildcardString | ) |
Definition at line 534 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
const ReactorExperiment * ana::WorldReactorConstraint2015 | ( | ) |
Weighted average of all experiments as of first nue paper writing.
Definition at line 64 of file ReactorExperiment.cxx.
const ReactorExperiment * ana::WorldReactorConstraint2016 | ( | ) |
Updated value for SecondAna based on the latest PDG.
Definition at line 74 of file ReactorExperiment.cxx.
const ReactorExperiment * ana::WorldReactorConstraint2017 | ( | ) |
Reactor constraint from PDG 2017 update.
Definition at line 81 of file ReactorExperiment.cxx.
void ana::WriteCAFMetadata | ( | TDirectory * | dir, |
const std::map< std::string, std::string > & | meta | ||
) |
Definition at line 674 of file sbnana/sbnana/CAFAna/Core/Utilities.cxx.
|
static |
Definition at line 89 of file DebugHelpers.cxx.
bool ana::gIsException = false |
Definition at line 25 of file DebugHelpers.cxx.
const Color_t ana::k1SigmaConfidenceColorIH = kOrange+9 |
Definition at line 54 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Color_t ana::k1SigmaConfidenceColorNH = kAzure+3 |
Definition at line 48 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Style_t ana::k1SigmaConfidenceStyle = 7 |
Definition at line 43 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Color_t ana::k2SigmaConfidenceColorIH = kOrange+10 |
Definition at line 55 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Color_t ana::k2SigmaConfidenceColorNH = kAzure+2 |
Definition at line 50 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Style_t ana::k2SigmaConfidenceStyle = 9 |
Definition at line 44 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Color_t ana::k3SigmaConfidenceColorIH = kOrange+8 |
Definition at line 56 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Color_t ana::k3SigmaConfidenceColorNH = kAzure+1 |
Definition at line 51 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Style_t ana::k3SigmaConfidenceStyle = 10 |
Definition at line 45 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Style_t ana::k68PercentConfidenceStyle = 7 |
Dashed.
Definition at line 31 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Color_t ana::k90PercConfidenceColorIH = kOrange+10 |
Definition at line 57 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Color_t ana::k90PercConfidenceColorNH = kAzure+2 |
Definition at line 49 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Style_t ana::k90PercentConfidenceStyle = 9 |
Definition at line 30 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Cut ana::kActiveVolumeFDCryo1 |
Definition at line 19 of file SBNAna/Cuts/Cuts.h.
const Cut ana::kActiveVolumeFDCryo2 |
Definition at line 20 of file SBNAna/Cuts/Cuts.h.
const Cut ana::kActiveVolumeND |
Definition at line 16 of file SBNAna/Cuts/Cuts.h.
const Var ana::kBaseline([](const caf::SRSliceProxy *slc) -> double{return slc->truth.baseline *1e-3;}) = 1284.9 |
Definition at line 41 of file CalcsNuFit.h.
const Binning ana::kBDTBinning = Binning::Simple(50,0, 1.0) |
Definition at line 40 of file Binnings.cxx.
const Color_t ana::kBeamNueBackgroundColor = kPink+9 |
Definition at line 18 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Float_t ana::kBlessedLabelFontSize = 0.06 |
Definition at line 63 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Float_t ana::kBlessedTitleFontSize = 0.08 |
Definition at line 62 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const std::vector<double> ana::kBLs = {kBaselineSBND, kBaselineMicroBoone, kBaselineIcarus} |
const Color_t ana::kCentralValueColor = kGray+2 |
Definition at line 59 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Color_t ana::kCentralValueColorIH = kOrange+9 |
Definition at line 53 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Color_t ana::kCentralValueColorNH = kAzure+3 |
Definition at line 47 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Style_t ana::kCentralValueStyle = kSolid |
Definition at line 60 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Var ana::kCompletness |
Definition at line 8 of file TruthVars.h.
const Cut ana::kContainedFD = kFiducialVolumeFDCryo1 || kFiducialVolumeFDCryo2 |
Definition at line 30 of file SBNAna/Cuts/Cuts.h.
const Cut ana::kCosmic_ClearCos = ( kCosmic_RFiducial && kNotClearCosmic ) |
Definition at line 112 of file NumuCutsIcarus202106.cxx.
Definition at line 94 of file NumuCutsIcarus202106.cxx.
const Cut ana::kCosmic_FMScore = ( kCosmic_NuScore && kFMScore ) |
Definition at line 124 of file NumuCutsIcarus202106.cxx.
const Cut ana::kCosmic_NuScore = ( kCosmic_ClearCos && kNuScore ) |
Definition at line 118 of file NumuCutsIcarus202106.cxx.
const Cut ana::kCosmic_PTrack = ( kCosmic_FMScore && kPTrack ) |
Definition at line 130 of file NumuCutsIcarus202106.cxx.
const Cut ana::kCosmic_RFiducial = ( kCosmic_Cryo0 && kRFiducial ) |
Definition at line 104 of file NumuCutsIcarus202106.cxx.
const Cut ana::kCosmic_TFiducial = ( kCosmic_Cryo0 && kTFiducial ) |
Definition at line 102 of file NumuCutsIcarus202106.cxx.
const Color_t ana::kCosmicBackgroundColor = kAzure+1 |
Definition at line 19 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Var ana::kCounting = kUnweighted |
Definition at line 11 of file NumuVarsIcarus202208.cxx.
const Var ana::kCRTHitDist |
Definition at line 16 of file NumuVarsSBND202106.h.
const Cut ana::kCRTHitDistanceCut |
Definition at line 13 of file NumuCutsSBND202106.h.
const SpillMultiVar ana::kCRTHitPE |
const SpillMultiVar ana::kCRTHitTime |
const SpillMultiVar ana::kCRTHitTimeFD |
const SpillCut ana::kCRTHitVetoFD |
Definition at line 14 of file SBNAna/Cuts/Cuts.h.
const SpillCut ana::kCRTHitVetoND |
Definition at line 13 of file SBNAna/Cuts/Cuts.h.
const SpillMultiVar ana::kCRTHitX |
const SpillMultiVar ana::kCRTHitY |
const SpillMultiVar ana::kCRTHitZ |
const Cut ana::kCRTTrackAngleCut |
Definition at line 12 of file NumuCutsSBND202106.h.
const Var ana::kCRTTrkAngle |
Definition at line 14 of file NumuVarsSBND202106.h.
const Var ana::kCRTTrkTime |
Definition at line 15 of file NumuVarsSBND202106.h.
const Binning ana::kCRTXFDBinning = Binning::Simple(60, crtfd.xmin, crtfd.xmax) |
Definition at line 36 of file Binnings.cxx.
const Binning ana::kCRTYFDBinning = Binning::Simple(100, crtfd.ymin, crtfd.ymax) |
Definition at line 37 of file Binnings.cxx.
const Binning ana::kCRTZFDBinning = Binning::Simple(130, crtfd.zmin, crtfd.zmax) |
Definition at line 38 of file Binnings.cxx.
const Cut ana::kCryo0 |
Definition at line 18 of file NumuCutsIcarus202106.h.
const double ana::kEarthDensity = 2.848 |
Definition at line 42 of file CalcsNuFit.h.
const EnergyScaleSyst ana::kEnergyScaleHadron |
Definition at line 64 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronBig |
Definition at line 96 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronFD |
Definition at line 76 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronFDBig |
Definition at line 108 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronInvSqrt |
Definition at line 66 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronInvSqrtBig |
Definition at line 98 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronInvSqrtFD |
Definition at line 78 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronInvSqrtFDBig |
Definition at line 110 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronInvSqrtND |
Definition at line 70 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronInvSqrtNDBig |
Definition at line 102 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronInvSqrtUB |
Definition at line 74 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronInvSqrtUBBig |
Definition at line 106 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronND |
Definition at line 68 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronNDBig |
Definition at line 100 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronSqrt |
Definition at line 65 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronSqrtBig |
Definition at line 97 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronSqrtFD |
Definition at line 77 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronSqrtFDBig |
Definition at line 109 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronSqrtND |
Definition at line 69 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronSqrtNDBig |
Definition at line 101 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronSqrtUB |
Definition at line 73 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronSqrtUBBig |
Definition at line 105 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronUB |
Definition at line 72 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleHadronUBBig |
Definition at line 104 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuon |
Definition at line 48 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonBig |
Definition at line 80 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonFD |
Definition at line 60 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonFDBig |
Definition at line 92 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonInvSqrt |
Definition at line 50 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonInvSqrtBig |
Definition at line 82 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonInvSqrtFD |
Definition at line 62 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonInvSqrtFDBig |
Definition at line 94 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonInvSqrtND |
Definition at line 54 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonInvSqrtNDBig |
Definition at line 86 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonInvSqrtUB |
Definition at line 58 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonInvSqrtUBBig |
Definition at line 90 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonND |
Definition at line 52 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonNDBig |
Definition at line 84 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonSqrt |
Definition at line 49 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonSqrtBig |
Definition at line 81 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonSqrtFD |
Definition at line 61 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonSqrtFDBig |
Definition at line 93 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonSqrtND |
Definition at line 53 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonSqrtNDBig |
Definition at line 85 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonSqrtUB |
Definition at line 57 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonSqrtUBBig |
Definition at line 89 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonUB |
Definition at line 56 of file EnergySysts.h.
const EnergyScaleSyst ana::kEnergyScaleMuonUBBig |
Definition at line 88 of file EnergySysts.h.
const SpillVar ana::kEvt = SIMPLESPILLVAR(hdr.evt) |
const Cut ana::kFiducialVolumeFDCryo1 |
Definition at line 21 of file SBNAna/Cuts/Cuts.h.
const Cut ana::kFiducialVolumeFDCryo2 |
Definition at line 22 of file SBNAna/Cuts/Cuts.h.
const Cut ana::kFiducialVolumeND |
Definition at line 17 of file SBNAna/Cuts/Cuts.h.
const Style_t ana::kFillStyleIH = 3003 |
Definition at line 41 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Style_t ana::kFillStyleNH = 3013 |
Definition at line 40 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const FitDelta13InPiUnitsSterile ana::kFitDelta13InPiUnitsSterile = FitDelta13InPiUnitsSterile() |
Definition at line 66 of file FitVarsSterile.h.
const FitDelta14InPiUnitsSterile ana::kFitDelta14InPiUnitsSterile = FitDelta14InPiUnitsSterile() |
Definition at line 81 of file FitVarsSterile.h.
const FitDelta24InPiUnitsSterile ana::kFitDelta24InPiUnitsSterile = FitDelta24InPiUnitsSterile() |
Definition at line 96 of file FitVarsSterile.h.
const FitDeltaInPiUnits ana::kFitDeltaInPiUnits = FitDeltaInPiUnits() |
const FitDmSq32Scaled ana::kFitDmSq32Scaled = FitDmSq32Scaled() |
const FitDmSq32Sterile ana::kFitDmSq32Sterile = FitDmSq32Sterile() |
Definition at line 21 of file FitVarsSterile.h.
const FitDmSq41Sterile ana::kFitDmSq41Sterile = FitDmSq41Sterile() |
Definition at line 36 of file FitVarsSterile.h.
const FitDmSq43Sterile ana::kFitDmSq43Sterile = FitDmSq43Sterile() |
Definition at line 51 of file FitVarsSterile.h.
const FitDmSqSterile ana::kFitDmSqSterile = FitDmSqSterile() |
m^2
Definition at line 23 of file FitVarsSterileApprox.h.
const FitSinSq2Theta12 ana::kFitSinSq2Theta12 = FitSinSq2Theta12() |
const FitSinSq2Theta13 ana::kFitSinSq2Theta13 = FitSinSq2Theta13() |
const FitSinSq2Theta14Sterile ana::kFitSinSq2Theta14Sterile = FitSinSq2Theta14Sterile() |
Definition at line 223 of file FitVarsSterile.h.
const FitSinSq2Theta23 ana::kFitSinSq2Theta23 = FitSinSq2Theta23() |
const FitSinSq2Theta24Sterile ana::kFitSinSq2Theta24Sterile = FitSinSq2Theta24Sterile() |
Definition at line 277 of file FitVarsSterile.h.
const FitSinSq2Theta34Sterile ana::kFitSinSq2Theta34Sterile = FitSinSq2Theta34Sterile() |
Definition at line 331 of file FitVarsSterile.h.
const FitSinSqTheta13Sterile ana::kFitSinSqTheta13Sterile = FitSinSqTheta13Sterile() |
Definition at line 132 of file FitVarsSterile.h.
const FitSinSqTheta14Sterile ana::kFitSinSqTheta14Sterile = FitSinSqTheta14Sterile() |
Definition at line 204 of file FitVarsSterile.h.
const FitSinSqTheta23 ana::kFitSinSqTheta23 = FitSinSqTheta23() |
const FitSinSqTheta23Sterile ana::kFitSinSqTheta23Sterile = FitSinSqTheta23Sterile() |
Definition at line 168 of file FitVarsSterile.h.
const FitSinSqTheta24Sterile ana::kFitSinSqTheta24Sterile = FitSinSqTheta24Sterile() |
Definition at line 259 of file FitVarsSterile.h.
const FitSinSqTheta34Sterile ana::kFitSinSqTheta34Sterile = FitSinSqTheta34Sterile() |
Definition at line 313 of file FitVarsSterile.h.
const FitTanSqTheta12 ana::kFitTanSqTheta12 = FitTanSqTheta12() |
const FitTheta13 ana::kFitTheta13 = FitTheta13() |
const FitTheta13InDegreesSterile ana::kFitTheta13InDegreesSterile = FitTheta13InDegreesSterile() |
Definition at line 349 of file FitVarsSterile.h.
const FitTheta13Sterile ana::kFitTheta13Sterile = FitTheta13Sterile() |
Definition at line 114 of file FitVarsSterile.h.
const FitTheta14InDegreesSterile ana::kFitTheta14InDegreesSterile = FitTheta14InDegreesSterile() |
Definition at line 385 of file FitVarsSterile.h.
const FitTheta14Sterile ana::kFitTheta14Sterile = FitTheta14Sterile() |
Definition at line 186 of file FitVarsSterile.h.
const FitTheta23 ana::kFitTheta23 = FitTheta23() |
const FitTheta23InDegreesSterile ana::kFitTheta23InDegreesSterile = FitTheta23InDegreesSterile() |
Definition at line 367 of file FitVarsSterile.h.
const FitTheta23Sterile ana::kFitTheta23Sterile = FitTheta23Sterile() |
Definition at line 150 of file FitVarsSterile.h.
const FitTheta24InDegreesSterile ana::kFitTheta24InDegreesSterile = FitTheta24InDegreesSterile() |
Definition at line 403 of file FitVarsSterile.h.
const FitTheta24Sterile ana::kFitTheta24Sterile = FitTheta24Sterile() |
Definition at line 241 of file FitVarsSterile.h.
const FitTheta34InDegreesSterile ana::kFitTheta34InDegreesSterile = FitTheta34InDegreesSterile() |
Definition at line 421 of file FitVarsSterile.h.
const FitTheta34Sterile ana::kFitTheta34Sterile = FitTheta34Sterile() |
Definition at line 295 of file FitVarsSterile.h.
const Cut ana::kFlashMatchNumuCut = kHasFlashMatch && kFlashMatchScore |
Definition at line 29 of file NumuCuts.cxx.
const Cut ana::kFlashMatchScore |
Definition at line 11 of file NumuCuts.h.
const SpillCut ana::kFlashTrigger |
Definition at line 11 of file SBNAna/Cuts/Cuts.h.
const Cut ana::kFMScore |
Definition at line 28 of file NumuCutsIcarus202106.h.
const Var ana::kFMScoreVar = SIMPLEVAR(fmatch.score) |
Definition at line 13 of file NumuVarsIcarus202208.cxx.
Definition at line 12 of file NumuVarsIcarus202208.cxx.
const Cut ana::kFullNueSel = kNueTrackLenCut && kShowerConvGapCut && kShowerdEdxCut && kShowerDensityCut |
const Cut ana::kHasFlashMatch |
Definition at line 9 of file NumuCuts.h.
const Cut ana::kHasMatchedNu |
Definition at line 9 of file CAFAna/Cuts/TruthCuts.h.
const Cut ana::kHasNu |
Definition at line 10 of file SBNAna/Cuts/TruthCuts.h.
const Cut ana::kHasPrimaryMuonTrk |
Definition at line 9 of file NumuCutsSBND202106.h.
const Var ana::kHasTruthMatch |
Definition at line 7 of file TruthVars.h.
const Cut ana::kIcarus202208ContainedHadrons |
Definition at line 16 of file NumuCutsIcarus202208.h.
const Cut ana::kIcarus202208ContainedMuon |
Definition at line 17 of file NumuCutsIcarus202208.h.
const Cut ana::kIcarus202208ContainedMuonAndHadrons = kIcarus202208ContainedMuon && kIcarus202208ContainedHadrons |
Definition at line 47 of file NumuCutsIcarus202208.cxx.
const Cut ana::kIcarus202208FMScoreCut = kFMScoreVar < 9 |
Definition at line 18 of file NumuCutsIcarus202208.cxx.
const Cut ana::kIcarus202208FMTimeCut = kFMTimeVar > 0 && kFMTimeVar < 1.8 |
Definition at line 17 of file NumuCutsIcarus202208.cxx.
const Cut ana::kIcarus202208FoundMuon = kIcarus202208MuonIdx >= 0 |
Definition at line 20 of file NumuCutsIcarus202208.cxx.
const Cut ana::kIcarus202208LongTrackDirCut = kCRLongestTrackDirY > -0.91 |
Definition at line 19 of file NumuCutsIcarus202208.cxx.
const Var ana::kIcarus202208MuonIdx |
Definition at line 11 of file NumuVarsIcarus202208.h.
const Cut ana::kIcarus202208NoPion = kIcarus202208NumPions == 0 |
Definition at line 33 of file NumuCutsIcarus202208.cxx.
const Var ana::kIcarus202208NumPions |
Definition at line 12 of file NumuVarsIcarus202208.h.
const Cut ana::kIcarus202208NumuSelection = kIcarus202208RecoFiducial && kIcarus202208FMScoreCut && kIcarus202208FMTimeCut && kIcarus202208LongTrackDirCut && kIcarus202208FoundMuon |
Definition at line 31 of file NumuCutsIcarus202208.cxx.
const Cut ana::kIcarus202208QELike = kIcarus202208NumuSelection && kIcarus202208NoPion && kIcarus202208ContainedHadrons |
Definition at line 48 of file NumuCutsIcarus202208.cxx.
const Cut ana::kIcarus202208QELikeContainedMuon = kIcarus202208QELike && kIcarus202208ContainedMuon |
Definition at line 49 of file NumuCutsIcarus202208.cxx.
const Cut ana::kIcarus202208RecoFiducial |
Definition at line 11 of file NumuCutsIcarus202208.h.
const Cut ana::kInBeamSpill |
Does the event fall inside the beam spill window.
const Cut ana::kInFV |
Definition at line 7 of file NumuCutsSBND202106.h.
const Cut ana::kInTimingSideband |
Is the event far from the start and end of the spill window.
const Color_t ana::kInvertedHierarchyColor = kOrange+9 |
Definition at line 28 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Var ana::kInvXSec([](const caf::SRSliceProxy *sr){const double GeV2perm2=2.56819e31;return GeV2perm2/sr->truth.xsec;}) |
const Cut ana::kIsAntiNu |
Is this truly an antineutrino?
Definition at line 70 of file CAFAna/Cuts/TruthCuts.h.
const Cut ana::kIsCC |
Definition at line 17 of file CAFAna/Cuts/TruthCuts.h.
const SpillCut ana::kIsCCSpill |
Definition at line 37 of file SBNAna/Cuts/TruthCuts.h.
const Cut ana::kIsCosmic = ( !kIsNuSlice ) |
Definition at line 10 of file NumuCutsIcarus202106.cxx.
const SpillCut ana::kIsCosmicSpill |
Definition at line 31 of file SBNAna/Cuts/TruthCuts.h.
const Cut ana::kIsNC |
Is this a Neutral Current event?
In this case the selection function is simple enough that we can include it inline as a lambda function.
Definition at line 15 of file CAFAna/Cuts/TruthCuts.h.
const SpillCut ana::kIsNCSpill |
Definition at line 38 of file SBNAna/Cuts/TruthCuts.h.
const Cut ana::kIsNu |
Definition at line 9 of file SBNAna/Cuts/TruthCuts.h.
const Cut ana::kIsNue |
Definition at line 12 of file SBNAna/Cuts/TruthCuts.h.
const SpillCut ana::kIsNueSpill |
Definition at line 33 of file SBNAna/Cuts/TruthCuts.h.
const Cut ana::kIsNumu |
Definition at line 13 of file SBNAna/Cuts/TruthCuts.h.
const Cut ana::kIsNuMuCC |
Definition at line 13 of file NumuCutsIcarus202106.h.
const SpillCut ana::kIsNumuSpill |
Definition at line 34 of file SBNAna/Cuts/TruthCuts.h.
const Cut ana::kIsNuOther = ( kIsNuSlice && !kIsNuMuCC ) |
Definition at line 16 of file NumuCutsIcarus202106.cxx.
const Cut ana::kIsNuSlice = ( kTruthIndex >= 0.f ) |
Definition at line 8 of file NumuCutsIcarus202106.cxx.
const Cut ana::kIsNutau |
Definition at line 14 of file SBNAna/Cuts/TruthCuts.h.
const SpillCut ana::kIsNutauSpill |
Definition at line 35 of file SBNAna/Cuts/TruthCuts.h.
const SpillCut ana::kIsSingleNuSpill |
Definition at line 32 of file SBNAna/Cuts/TruthCuts.h.
const SpillCut ana::kLongTrack |
Definition at line 83 of file NumuCutsIcarus202106.h.
const Binning ana::kLowEnergyBinning = Binning::Custom( kLowEnergyEdges ) |
Definition at line 25 of file Binnings.cxx.
const std::vector<double> ana::kLowEnergyEdges |
Definition at line 9 of file Binnings.cxx.
const Binning ana::kLowEnergyGeVBinning = Binning::Custom( kLowEnergyGeVEdges ) |
Definition at line 26 of file Binnings.cxx.
const std::vector<double> ana::kLowEnergyGeVEdges |
Definition at line 15 of file Binnings.cxx.
const Var ana::kMuMaxTrack |
Definition at line 11 of file NumuVarsIcarus202106.h.
const Color_t ana::kNCBackgroundColor = kAzure |
Definition at line 17 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Color_t ana::kNormalHierarchyColor = kAzure+3 |
Definition at line 27 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const SystShifts ana::kNoShift = SystShifts::Nominal() |
Definition at line 61 of file SystShifts.h.
const Cut ana::kNotClearCosmic |
Definition at line 24 of file NumuCutsIcarus202106.h.
const Binning ana::kNueEnergyBinning = Binning::Simple( 20, 0., 5000.) |
Definition at line 22 of file Binnings.cxx.
const Cut ana::kNueFD = kContainedFD && kNueFlashScoreFDCut && kNuePandoraScoreFDCut && kRecoShowerFD |
const Cut ana::kNueFlashScoreFDCut = kSlcFlashMatchCut |
const Cut ana::kNueLongestTrackDazzleMuonCut = kLongestTrackDazzlePID != 13 |
Definition at line 131 of file NueCuts.cxx.
const Cut ana::kNueLongestTrackDazzleMuonScoreCut = kLongestTrackDazzleMuonScore < 0.8 |
Definition at line 132 of file NueCuts.cxx.
const Cut ana::kNueMuonCutNEW |
const Cut ana::kNueMuonCutOLD |
const Cut ana::kNuePandoraScoreFDCut = kSlcNuScoreCut |
const Color_t ana::kNueSignalColor = kViolet-5 |
Definition at line 15 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const double ana::kNuFitdCPCVIH = 284 * TMath::Pi()/180 |
Definition at line 26 of file CalcsNuFit.h.
const double ana::kNuFitdCPCVNH = 215 * TMath::Pi()/180 |
Definition at line 21 of file CalcsNuFit.h.
const double ana::kNuFitDmsq21CV = 7.39e-5 |
Definition at line 14 of file CalcsNuFit.h.
const double ana::kNuFitDmsq21Err = ((8.01-6.79)/6)*1e-5 |
Definition at line 29 of file CalcsNuFit.h.
const double ana::kNuFitDmsq32CVIH = -2.512e-3 |
Definition at line 23 of file CalcsNuFit.h.
const double ana::kNuFitDmsq32CVNH = +2.525e-3 - kNuFitDmsq21CV |
Definition at line 18 of file CalcsNuFit.h.
const double ana::kNuFitDmsq32ErrIH = ((2.611-2.412)/6)*1e-3 |
Definition at line 36 of file CalcsNuFit.h.
const double ana::kNuFitDmsq32ErrNH = ((2.625-2.427)/6)*1e-3 |
Definition at line 32 of file CalcsNuFit.h.
const double ana::kNuFitTh12CV = 33.82 * TMath::Pi()/180 |
Definition at line 15 of file CalcsNuFit.h.
const double ana::kNuFitTh12Err = ((36.27-31.61)/6) * TMath::Pi()/180 |
Definition at line 30 of file CalcsNuFit.h.
const double ana::kNuFitTh13CVIH = 8.65 * TMath::Pi()/180 |
Definition at line 25 of file CalcsNuFit.h.
const double ana::kNuFitTh13CVNH = 8.61 * TMath::Pi()/180 |
Definition at line 20 of file CalcsNuFit.h.
const double ana::kNuFitTh13ErrIH = ((9.03-8.27)/6) * TMath::Pi()/180 |
Definition at line 38 of file CalcsNuFit.h.
const double ana::kNuFitTh13ErrNH = ((8.99-8.22)/6) * TMath::Pi()/180 |
Definition at line 34 of file CalcsNuFit.h.
const double ana::kNuFitTh23CVIH = 49.8 * TMath::Pi()/180 |
Definition at line 24 of file CalcsNuFit.h.
const double ana::kNuFitTh23CVNH = 49.6 * TMath::Pi()/180 |
Definition at line 19 of file CalcsNuFit.h.
const double ana::kNuFitTh23ErrIH = ((52.5-40.6)/6) * TMath::Pi()/180 |
Definition at line 37 of file CalcsNuFit.h.
const double ana::kNuFitTh23ErrNH = ((52.4-40.3)/6) * TMath::Pi()/180 |
Definition at line 33 of file CalcsNuFit.h.
|
static |
Dummy loader that doesn't load any files.
Useful when a loader is required for a component you want to ignore
Definition at line 269 of file SpectrumLoaderBase.h.
const Color_t ana::kNumuBackgroundColor = kGreen+2 |
This is the color for the numu CC background to the nue analysis. I doubt numu plots will ever use this, preferring to just show the total MC.
Definition at line 22 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Cut ana::kNumuBasicQual |
Definition at line 7 of file NumuCuts.h.
const Cut ana::kNuMuCC_ClearCos = ( kNuMuCC_TFiducial && kNotClearCosmic ) |
Definition at line 110 of file NumuCutsIcarus202106.cxx.
Definition at line 92 of file NumuCutsIcarus202106.cxx.
const Cut ana::kNuMuCC_FMScore = ( kNuMuCC_NuScore && kFMScore ) |
Definition at line 122 of file NumuCutsIcarus202106.cxx.
const Cut ana::kNuMuCC_FullSelection = ( kCryo0 && kRFiducial && kNotClearCosmic && kNuScore && kFMScore && kPTrack ) |
Definition at line 135 of file NumuCutsIcarus202106.cxx.
const Cut ana::kNuMuCC_NuScore = ( kNuMuCC_ClearCos && kNuScore ) |
Definition at line 116 of file NumuCutsIcarus202106.cxx.
const Cut ana::kNuMuCC_PTrack = ( kNuMuCC_FMScore && kPTrack ) |
Definition at line 128 of file NumuCutsIcarus202106.cxx.
const Cut ana::kNuMuCC_RFiducial = ( kNuMuCC_Cryo0 && kRFiducial ) |
Definition at line 100 of file NumuCutsIcarus202106.cxx.
const Cut ana::kNuMuCC_TFiducial = ( kNuMuCC_Cryo0 && kTFiducial ) |
Definition at line 98 of file NumuCutsIcarus202106.cxx.
const Binning ana::kNumuEnergyBinning = Binning::Simple( 20, 0., 5000.) |
Definition at line 23 of file Binnings.cxx.
const Cut ana::kNumuTrkLen |
Definition at line 15 of file NumuCuts.h.
const Cut ana::kNuOther_ClearCos = ( kNuOther_TFiducial && kNotClearCosmic ) |
Definition at line 114 of file NumuCutsIcarus202106.cxx.
const Cut ana::kNuOther_Cryo0 = ( kIsNuOther && kCryo0 ) |
Definition at line 96 of file NumuCutsIcarus202106.cxx.
const Cut ana::kNuOther_FMScore = ( kNuOther_NuScore && kFMScore ) |
Definition at line 126 of file NumuCutsIcarus202106.cxx.
const Cut ana::kNuOther_NuScore = ( kNuOther_ClearCos && kNuScore ) |
Definition at line 120 of file NumuCutsIcarus202106.cxx.
const Cut ana::kNuOther_PTrack = ( kNuOther_FMScore && kPTrack ) |
Definition at line 132 of file NumuCutsIcarus202106.cxx.
const Cut ana::kNuOther_RFiducial = ( kNuOther_Cryo0 && kRFiducial ) |
Definition at line 108 of file NumuCutsIcarus202106.cxx.
const Cut ana::kNuOther_TFiducial = ( kNuOther_Cryo0 && kTFiducial ) |
Definition at line 106 of file NumuCutsIcarus202106.cxx.
const Var ana::kNuScore |
Definition at line 26 of file NumuCutsIcarus202106.h.
const Binning ana::kPositionXFDBinning = Binning::Simple(47, avfd_cryo1.xmin, avfd_cryo1.xmax) |
Definition at line 32 of file Binnings.cxx.
const Binning ana::kPositionXNDBinning = Binning::Simple(40, avnd.xmin, avnd.xmax) |
Definition at line 29 of file Binnings.cxx.
const Binning ana::kPositionYFDBinning = Binning::Simple(35, avfd_cryo1.ymin, avfd_cryo1.ymax) |
Definition at line 33 of file Binnings.cxx.
const Binning ana::kPositionYNDBinning = Binning::Simple(40, avnd.ymin, avnd.ymax) |
Definition at line 30 of file Binnings.cxx.
const Binning ana::kPositionZFDBinning = Binning::Simple(90, avfd_cryo1.zmin, avfd_cryo1.zmax) |
Definition at line 34 of file Binnings.cxx.
const Binning ana::kPositionZNDBinning = Binning::Simple(50, avnd.zmin, avnd.zmax) |
Definition at line 31 of file Binnings.cxx.
const Cut ana::kPreNueSelND = kFiducialVolumeND && kSlcNuScoreCut && kSlcFlashMatchCut |
Definition at line 9 of file NumuVarsIcarus202106.cxx.
const Var ana::kPrimaryMuonTrkIdx |
Definition at line 11 of file NumuVarsSBND202106.h.
const Var ana::kPrimaryMuonTrkLen |
Definition at line 13 of file NumuVarsSBND202106.h.
const Var ana::kPrimaryMuonTrkP |
Definition at line 12 of file NumuVarsSBND202106.h.
const Color_t ana::kPrimColorIH = kOrange+9 |
Definition at line 37 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Color_t ana::kPrimColorNH = kAzure+2 |
Definition at line 34 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Var ana::kPrimMuonIdx |
Definition at line 11 of file NumuVars.h.
const Var ana::kPrimTrkLen |
Definition at line 9 of file NumuVars.h.
const Cut ana::kPTrack |
Definition at line 30 of file NumuCutsIcarus202106.h.
const Cut ana::kPTrackContained |
Definition at line 32 of file NumuCutsIcarus202106.h.
const Cut ana::kPTrackExiting |
Definition at line 34 of file NumuCutsIcarus202106.h.
const Var ana::kPTrackInd |
Definition at line 13 of file NumuVarsIcarus202106.h.
const Cut ana::kRazzleDazzleNueSel = kNueAllTrackDazzleMuonCut && kShowerRazzleElectronScoreCut |
const Var ana::kRecoMuonP |
Definition at line 15 of file NumuVarsIcarus202106.h.
const Cut ana::kRecoNueSel = kRecoShower && kShowerEnergyCut |
const Cut ana::kRecoShowerFD = kRecoShower && kNueNumShowersCut && kShowerdEdxCut && kShowerConvGapCut && kNueTrackLenCut && kShowerDensityCut && kShowerEnergyCut |
const Cut ana::kRFiducial |
Definition at line 22 of file NumuCutsIcarus202106.h.
const SpillVar ana::kRun = SIMPLESPILLVAR(hdr.run) |
const Color_t ana::kSecoColorIH = kOrange+10 |
Definition at line 38 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Color_t ana::kSecoColorNH = kAzure+1 |
Definition at line 35 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Cut ana::kShowerRazzleElectronCut = kRecoShowerRazzlePID == 11 |
Definition at line 128 of file NueCuts.cxx.
const Cut ana::kShowerRazzleElectronScoreCut = kRecoShowerRazzleElectronScore > 0.8 |
Definition at line 129 of file NueCuts.cxx.
const Cut ana::kSlcCompletenessCut |
Definition at line 29 of file SBNAna/Cuts/TruthCuts.h.
const Cut ana::kSlcFlashMatchCut |
Definition at line 28 of file SBNAna/Cuts/Cuts.h.
const Cut ana::kSlcFlashMatchScoreCut |
Definition at line 15 of file NumuCutsSBND202106.h.
const Cut ana::kSlcFlashMatchTimeCut |
Definition at line 14 of file NumuCutsSBND202106.h.
const Cut ana::kSlcHasFlashMatch |
Definition at line 27 of file SBNAna/Cuts/Cuts.h.
const Cut ana::kSlcIsRecoNu |
Definition at line 24 of file SBNAna/Cuts/Cuts.h.
const Cut ana::kSlcNuScoreCut |
Definition at line 25 of file SBNAna/Cuts/Cuts.h.
const SpillVar ana::kSpillCounting = kSpillUnweighted |
const Cut ana::kTFiducial |
Definition at line 20 of file NumuCutsIcarus202106.h.
const Color_t ana::kTotalBackgroundColor = kAzure+2 |
Definition at line 16 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Color_t ana::kTotalMCColor = kRed |
Definition at line 13 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Color_t ana::kTotalMCErrorBandColor = kRed-10 |
Definition at line 14 of file sbnana/sbnana/CAFAna/Analysis/Style.h.
const Cut ana::kTrueActiveVolumeFDCryo1 |
Definition at line 21 of file SBNAna/Cuts/TruthCuts.h.
const Cut ana::kTrueActiveVolumeFDCryo2 |
Definition at line 22 of file SBNAna/Cuts/TruthCuts.h.
const Cut ana::kTrueActiveVolumeND |
Definition at line 19 of file SBNAna/Cuts/TruthCuts.h.
const Var ana::kTrueE([](const caf::SRSliceProxy *slc) -> double{return slc->truth.E;}) |
const Binning ana::kTrueEnergyBins = TrueEnergyBins() |
const Cut ana::kTrueFiducialVolumeFDCryo1 |
Definition at line 24 of file SBNAna/Cuts/TruthCuts.h.
const Cut ana::kTrueFiducialVolumeFDCryo2 |
Definition at line 25 of file SBNAna/Cuts/TruthCuts.h.
const Cut ana::kTrueFiducialVolumeND |
Definition at line 20 of file SBNAna/Cuts/TruthCuts.h.
Definition at line 30 of file OscillatableSpectrum.cxx.
const Binning ana::kTrueLOverEBins = Binning::Simple(100, 0, 2) |
const Var ana::kTrueMuonP |
Definition at line 17 of file NumuVarsIcarus202106.h.
const Var ana::kTruthEnergy |
Definition at line 10 of file TruthVars.h.
Definition at line 7 of file NumuVarsIcarus202106.cxx.
const SpillVar ana::kTruthLeptonEnergy |
Definition at line 23 of file TruthVars.h.
const SpillVar ana::kTruthNuEnergy |
Definition at line 22 of file TruthVars.h.
const Var ana::kTruthVtxDistMag |
Definition at line 20 of file TruthVars.h.
const Var ana::kTruthVtxDistX |
Definition at line 16 of file TruthVars.h.
const Var ana::kTruthVtxDistY |
Definition at line 17 of file TruthVars.h.
const Var ana::kTruthVtxDistZ |
Definition at line 18 of file TruthVars.h.
const Var ana::kTruthVtxX |
Definition at line 12 of file TruthVars.h.
const Var ana::kTruthVtxY |
Definition at line 13 of file TruthVars.h.
const Var ana::kTruthVtxZ |
Definition at line 14 of file TruthVars.h.
const Cut ana::kVtxDistMagCut |
Definition at line 27 of file SBNAna/Cuts/TruthCuts.h.
f sin ana::theta_ |
Definition at line 27 of file FitVarsSterileApprox.h.