Classes | Typedefs | Enumerations | Functions | Variables
cmf Namespace Reference

Module to combine a set of results into a single file. More...

Classes

class  AnalyzeWeighters
 
class  BeamSyst
 
class  CAFToEventList
 
class  CalibSyst
 
class  CalibSystHistMaker
 
struct  ChiSqrCalcGridPoint
 
class  CombineResults
 
struct  ConstraintInfo
 
class  ContourMaker
 
class  CovarianceBinUtility
 
class  CovarianceFitHelper
 
class  CovarianceMatrixDecorrelator
 
class  CovarianceMatrixFitter
 
class  CovarianceMatrixMaker
 
class  CovarianceTools
 
struct  CovFitIteration
 
struct  DataVars
 
class  DataVarVals
 Base container for the MC related Vars that constitute an event. More...
 
class  Event
 
class  EventId
 
class  EventList
 
class  EventListManipulator
 
class  EventListToTextFile
 
class  FakePointResultMaker
 
class  FakeUniverse
 
class  GeniePCASyst
 
struct  GENIEVar
 
struct  GENIEVars
 
struct  GridPoint
 
struct  GridPointResults
 
class  InputPoint
 
class  MCVarVals
 Base container for the MC related Vars that constitute an event. More...
 
class  MedianSensitivityPointsMaker
 
struct  MetaData
 
struct  MetaDataLite
 
struct  NPInfo
 
class  NuisanceParameters
 
class  ParameterConstraints
 
struct  ParameterDetHasher
 
class  PlotPoint
 
class  PlotUtilities
 
class  PointResult
 
struct  POTLiveTime
 
class  RandomUniverses
 
class  ReadFitResults
 
struct  SelectOutput
 
class  SensitivityPointsMaker
 
struct  shiftconv
 
class  ShifterAndWeighter
 
struct  SpillSummary
 
class  StarPlot
 
struct  SubMatrixData
 
class  SystematicSigmas
 
struct  TruthVars
 
struct  WeightVars
 

Typedefs

typedef std::unique_ptr< EventEventPtr
 
typedef std::vector< EventPtrEventContainer
 
typedef std::map< cmf::MetaData, cmf::EventListEventListMap
 
typedef enum cmf::weight_type WeightType_t
 
typedef std::pair< std::string, novadaq::cnv::DetIdParameterDet
 
typedef std::unordered_map< std::string, double > ParameterMap
 
typedef std::unordered_map< ParameterDet, double, ParameterDetHasherParameterDetMap
 
typedef enum cmf::beam_type BeamType_t
 
typedef enum cmf::chisqr_type ChiSqrCalc_t
 
typedef enum cmf::sel_type SelectionType_t
 
typedef enum cmf::interaction_type InteractionType_t
 
typedef enum cmf::osc_params OscParm_t
 
typedef enum cmf::osc_type OscType_t
 
typedef enum cmf::calc_type CalcType_t
 
typedef enum cmf::file_type FileType_t
 
typedef std::map< cmf::ChiSqrCalcGridPoint, cmf::GridPointResultsGridPointResultsMap
 
typedef std::map< std::pair< double, double >, std::set< cmf::PointResult > > PointResMap
 
typedef std::map< std::pair< double, double >, double > PointMap
 
typedef enum cmf::dataMC_type DataMC_t
 
typedef std::map< cmf::MetaData, double > ExposureMap
 
typedef std::map< cmf::MetaData, std::vector< TH1F * > > EnergyHistMap
 

Enumerations

enum  weight_type { kSystWgtOnly =0, kOscWgtOnly, kSystAndOscWgt, kAllWgt }
 
enum  beam_type { kFHC = 0, kRHC, k0HC, kUnknownBeam }
 
enum  chisqr_type {
  kCovMat = 0, kPoisson, kCNP, kGauss,
  kUnknownChiSqr
}
 
enum  sel_type {
  kNuMuSelection = 0, kNuESelection, kNuESelectionLowPID, kNuESelectionMidPID,
  kNuESelectionHighPID, kNuESelectionPeripheral, kNuMuSelectionQ1, kNuMuSelectionQ2,
  kNuMuSelectionQ3, kNuMuSelectionQ4, kNCSelection, kUnknownSelection
}
 
enum  interaction_type {
  kUnknownInteraction = 0, kNuMuCC, kNuMuBarCC, kNuECC,
  kNuEBarCC, kNuTauCC, kNuTauBarCC, kNC,
  kCosmicMuon, kRockMuon
}
 
enum  osc_params {
  kL = 0, kRho, kDmsq21, kDmsq32,
  kTh12, kTh13, kTh23, kdCP,
  kTh24, kTh34, kDmsq41, kNumOscParams
}
 
enum  osc_type {
  kMCmm = 0, kMCmmbar, kMCme, kMCmebar,
  kMCmt, kMCmtbar, kMCee, kMCeebar,
  kMCem, kMCembar, kMCtm, kMCtmbar,
  kMCte, kMCtebar, kMCNC, kDataEvent,
  kCosmicEvent, kBadOscType, kNumOscTypes, kOscTypeNoneSet
}
 
enum  calc_type { kPMNS =0, kNSI, kSterile }
 
enum  file_type {
  kBeam =0, kSwap, kTauSwap, kDataFile,
  kCosmicBackgroundFile, kRockFluxSwap, kRockNonSwap, kNumFileTypes,
  kBadFileType, kUnknownFileType
}
 
enum  dataMC_type { kData = 0, kMC, kFakeData, kBoth }
 

Functions

std::ostream & operator<< (std::ostream &o, cmf::Event const &event)
 
void SerializeEvents (TTree *eventTree, cmf::MetaData const &md, cmf::EventList const &events)
 
void SerializeEventListMap (std::string const &dirName, cmf::EventListMap const &listMap)
 
std::ostream & operator<< (std::ostream &o, cmf::MCVarVals const &mcvv)
 
std::ostream & operator<< (std::ostream &o, cmf::EventId const &evid)
 
std::ostream & operator<< (std::ostream &o, cmf::DataVarVals const &dvv)
 
bool isSumSmallXSecJoint2017Key (unsigned char const &varkey)
 
unsigned char VarNameToKey (std::string const &name)
 
std::string KeyToVarName (unsigned char const &key)
 
std::string KeyToLatex (unsigned char const &varkey)
 
bool IsMCVar (unsigned char const &key)
 
bool IsMCVar (std::string const &name)
 
bool IsDataVar (unsigned char const &key)
 
bool IsDataVar (std::string const &name)
 
std::string ParameterDetToString (cmf::ParameterDet const &pd)
 
std::ostream & operator<< (std::ostream &o, cmf::InputPoint const &ip)
 
std::ostream & operator<< (std::ostream &o, cmf::PointResult const &pr)
 
static long DetectorTypeKey (novadaq::cnv::DetId const &det)
 
static long EpochKey (int const &ep)
 
static long FileTypeKey (cmf::FileType_t const &ft)
 
static long BeamTypeKey (cmf::BeamType_t const &bt)
 
static long SelectionTypeKey (cmf::SelectionType_t const &st)
 
static long InteractionTypeKey (cmf::InteractionType_t const &it)
 
static long DetectorBeamSelectionTypesToKey (novadaq::cnv::DetId const &det, cmf::BeamType_t const &bt, cmf::SelectionType_t const &sel)
 
static novadaq::cnv::DetId KeyToDetectorType (long const &key)
 
static int KeyToEpoch (long const &key)
 
static cmf::FileType_t KeyToFileType (long const &key)
 
static cmf::BeamType_t KeyToBeamType (long const &key)
 
static cmf::SelectionType_t KeyToSelectionType (long const &key)
 
static cmf::InteractionType_t KeyToInteractionType (long const &key)
 
static cmf::SelectionType_t StringToSelectionType (std::string const &str)
 
static cmf::BeamType_t StringToBeamType (std::string const &str)
 
static novadaq::cnv::DetId StringToDetectorType (std::string const &str)
 
static cmf::ChiSqrCalc_t StringToChiSqrType (std::string const &str)
 
long StringToKey (std::string const &str)
 
int EpochStringToInt (std::string const &epochStr)
 
std::ostream & operator<< (std::ostream &o, SpillSummary const &ss)
 
std::ostream & operator<< (std::ostream &o, cmf::SelectOutput const &so)
 
double FitFunction_ForROOT (const double *params)
 
template<class T >
static bool readSingleResult (art::Results const &r, std::string const &module, std::string const &instance, std::vector< T > &resvec)
 
template<class T >
static bool readVectorResults (art::Results const &r, std::string const &module, std::string const &instance, std::vector< T > &resvec)
 

Variables

static ShifterAndWeightergShifterAndWeighter = nullptr
 
static const unsigned char kUnknownVar = 0
 
static const unsigned char kNuE_CVN = 1
 
static const unsigned char kNuE_NumMichel = 2
 
static const unsigned char kHad_RecoE = 3
 
static const unsigned char kLep_RecoE = 4
 
static const unsigned char kNu_RecoE = 5
 
static const unsigned char kLep_RecoE_MCFrac = 6
 
static const unsigned char kRecoQ2 = 7
 
static const unsigned char kFake_Weight = 8
 
static const unsigned char kTrueE = 9
 
static const unsigned char kTruePDG = 10
 
static const unsigned char kTrueCCNC = 11
 
static const unsigned char kTrueIntType = 12
 
static const unsigned char kTruePDGOrig = 13
 
static const unsigned char kTrueParentPDG = 14
 
static const unsigned char kTrueParentPT = 15
 
static const unsigned char kTrueParentPZ = 16
 
static const unsigned char kTrueParentDecay = 17
 
static const unsigned char kTrueParentTargetPDG = 18
 
static const unsigned char kTrueHitNuc = 19
 
static const unsigned char kTrueIntMode = 20
 
static const unsigned char kXSecCVPPFX_Weight = 21
 
static const unsigned char kRPACCQE_Weight = 22
 
static const unsigned char kRPARES_Weight = 23
 
static const unsigned char kDISvnCC1pi_Weight = 24
 
static const unsigned char kEmpiricalMEC_Weight = 25
 
static const unsigned char kEmpiricalMECtoGENIEQE_Weight = 26
 
static const unsigned char kEmpiricalMECtoGENIERES_Weight = 27
 
static const unsigned char kOtherDIS_Weight = 28
 
static const unsigned char kMinus2Sigma = 29
 
static const unsigned char kMinus1Sigma = 30
 
static const unsigned char kPlus1Sigma = 31
 
static const unsigned char kPlus2Sigma = 32
 
static const unsigned char kMaCCQE = 33
 
static const unsigned char kMaCCRES = 34
 
static const unsigned char kMaNCRES = 35
 
static const unsigned char kMvCCRES = 36
 
static const unsigned char kRPACCQEshapeEnh = 37
 
static const unsigned char kRPACCQEshapeSupp = 38
 
static const unsigned char kFormZone = 39
 
static const unsigned char kFrInel_pi = 40
 
static const unsigned char kFrCEx_N = 41
 
static const unsigned char kFrElas_N = 42
 
static const unsigned char kFrAbs_N = 43
 
static const unsigned char kCCQEPauliSupViaKF = 44
 
static const unsigned char kSumSmallXSecJoint2017 = 45
 
static const unsigned char kMaNCEL = 46
 
static const unsigned char kEtaNCEL = 47
 
static const unsigned char kVecCCQE_shape = 48
 
static const unsigned char kMvNCRES = 49
 
static const unsigned char kMaCOHpi = 50
 
static const unsigned char kR0COHpi = 51
 
static const unsigned char kAhtBY = 52
 
static const unsigned char kBhtBY = 53
 
static const unsigned char kCV1uBY = 54
 
static const unsigned char kCV2uBY = 55
 
static const unsigned char kNormDISCC = 56
 
static const unsigned char kRnubarnuCC = 57
 
static const unsigned char kDISNuclMod = 58
 
static const unsigned char kAGKY_xF1pi = 59
 
static const unsigned char kAGKY_pT1pi = 60
 
static const unsigned char kMFP_pi = 61
 
static const unsigned char kMFP_N = 62
 
static const unsigned char kFrCEx_pi = 63
 
static const unsigned char kFrElas_pi = 64
 
static const unsigned char kFrAbs_pi = 65
 
static const unsigned char kFrPiProd_pi = 66
 
static const unsigned char kFrInel_N = 67
 
static const unsigned char kFrPiProd_N = 68
 
static const unsigned char kCCQEMomDistroFGtoSF = 69
 
static const unsigned char kBR1gamma = 70
 
static const unsigned char kBR1eta = 71
 
static const unsigned char kTheta_Delta2Npi = 72
 
static const unsigned char kLep_RecoE_UBL = 73
 
static const unsigned char kHad_RecoE_UBL = 74
 
static const unsigned char kNu_RecoE_UBL = 75
 
static const unsigned char kNu_RecoE_Quantiles_NuMu = 76
 
static const unsigned char kNu_RecoE_NuE_Peripheral = 77
 
const std::string cBeamType_Strings [4]
 
const std::string cBeamType_LatexStrings [4]
 
const std::string cChiSqrCalcType_Strings [5]
 
const std::string cSelectionType_Strings [12]
 
const std::string cSelectionType_LatexStrings [12]
 
const std::string cInteractionType_LatexStrings [10]
 
const std::string cInteractionType_Strings [10]
 
const double cOscParams_Scales [kNumOscParams]
 
const std::string cOscParams_LatexScales [kNumOscParams]
 
const std::string cOscParams_Strings [kNumOscParams]
 
const std::string cOscParams_Strings_Latex [kNumOscParams]
 
const std::string cOscType_Strings [20]
 
const std::string cOscType_LatexStrings [20]
 
const int cOscType_Colors [20]
 
const int cOscType_FillStyles [20]
 
const int cOscType_Markers [20]
 
const std::string cFileTypeStrings [10]
 
const std::string cFileTypeStrings_Latex [10]
 
const int kNumFCAlpha = 2
 
const double kFCAlpha [2] ={0.68, 0.90}
 
const int kFCAlphaColors [2] = {4, 2}
 
const std::string kFCAlphaStrings [2] = {"68% CL","90% CL"}
 
static NuisanceParametersgSystSigmaInstance = nullptr
 
static ParameterConstraintsgParConstInstance = nullptr
 
static SystematicSigmasgSystSigmaInstance = nullptr
 
static CovarianceFitHelpergCovFitHelp = nullptr
 
std::unique_ptr< cmf::EventListMapgMCEventLists
 
std::unique_ptr< cmf::InputPointgInputPoint
 
std::unique_ptr< std::vector< double > > gFFVals
 
std::unique_ptr< std::vector< std::string > > gFFNames
 
ROOT::Math::Minimizer * gFFMin
 
cmf::ChiSqrCalc_t gChiSqrCalc
 
double gMinChiSqr
 
int gMinStatus
 
static const double gPI = 3.141592653589793238462643383279502884
 
static PlotUtilitiesgPlotUtil = nullptr
 
static CovarianceBinUtilitygCovBinUtil = nullptr
 

Detailed Description

Module to combine a set of results into a single file.

Create a set of Covariance matrices to be used in fits.

Use covariance matrices in fits.

Module to do a fit of multiple random data draws at a given point in parameter space.

Module to do a fit of a single point in parameter space.

Author
brebe.nosp@m.l@fn.nosp@m.al.go.nosp@m.v

Typedef Documentation

typedef std::map<cmf::MetaData, std::vector<TH1F*> > cmf::EnergyHistMap

Definition at line 41 of file EventListManipulator.h.

typedef std::vector<EventPtr> cmf::EventContainer

Definition at line 96 of file Event.h.

Definition at line 148 of file Event.h.

typedef std::unique_ptr<Event> cmf::EventPtr

Definition at line 95 of file Event.h.

typedef std::map<cmf::MetaData, double > cmf::ExposureMap

Definition at line 40 of file EventListManipulator.h.

Definition at line 154 of file PointResult.h.

typedef std::pair<std::string, novadaq::cnv::DetId> cmf::ParameterDet

Definition at line 38 of file Constants.h.

typedef std::unordered_map<ParameterDet, double, ParameterDetHasher> cmf::ParameterDetMap

Definition at line 63 of file Constants.h.

typedef std::unordered_map<std::string, double> cmf::ParameterMap

Definition at line 62 of file Constants.h.

typedef std::map< cmf::GridPoint, double > cmf::PointMap

Definition at line 41 of file CMFReadFitResults_plugin.cc.

typedef std::map< std::pair< double, double >, std::set< cmf::PointResult > > cmf::PointResMap

Definition at line 40 of file CMFReadFitResults_plugin.cc.

Enumeration Type Documentation

Enumerator
kFHC 
kRHC 
k0HC 
kUnknownBeam 

Definition at line 68 of file Constants.h.

68  {
69  kFHC = 0,
70  kRHC,
71  k0HC,
73  } BeamType_t;
enum fnex::beam_type BeamType_t
Enumerator
kPMNS 
kNSI 
kSterile 

Definition at line 333 of file Constants.h.

333  {
334  kPMNS=0,
335  kNSI,
336  kSterile
337  } CalcType_t;
enum cmf::calc_type CalcType_t
Enumerator
kCovMat 
kPoisson 
kCNP 
kGauss 
kUnknownChiSqr 

Definition at line 85 of file Constants.h.

85  {
86  kCovMat = 0,
87  kPoisson,
88  kCNP,
89  kGauss,
91  } ChiSqrCalc_t;
enum cmf::chisqr_type ChiSqrCalc_t
Enumerator
kData 
kMC 
kFakeData 
kBoth 

Definition at line 33 of file EventListManipulator.h.

33  {
34  kData = 0,
35  kMC,
36  kFakeData,
37  kBoth
38  } DataMC_t;
const Color_t kMC
enum fnex::dataMC_type DataMC_t
Enumerator
kBeam 
kSwap 
kTauSwap 
kDataFile 
kCosmicBackgroundFile 
kRockFluxSwap 
kRockNonSwap 
kNumFileTypes 
kBadFileType 
kUnknownFileType 

Definition at line 367 of file Constants.h.

Enumerator
kUnknownInteraction 
kNuMuCC 
kNuMuBarCC 
kNuECC 
kNuEBarCC 
kNuTauCC 
kNuTauBarCC 
kNC 
kCosmicMuon 
kRockMuon 

Definition at line 144 of file Constants.h.

Enumerator
kL 
kRho 
kDmsq21 
kDmsq32 
kTh12 
kTh13 
kTh23 
kdCP 
kTh24 
kTh34 
kDmsq41 
kNumOscParams 

Definition at line 180 of file Constants.h.

180  {
181  kL = 0,
182  kRho,
183  kDmsq21,
184  kDmsq32,
185  kTh12,
186  kTh13,
187  kTh23,
188  kdCP,
189  kTh24,
190  kTh34,
191  kDmsq41,
193  } OscParm_t;
enum fnex::osc_params OscParm_t
Enumerator
kMCmm 
kMCmmbar 
kMCme 
kMCmebar 
kMCmt 
kMCmtbar 
kMCee 
kMCeebar 
kMCem 
kMCembar 
kMCtm 
kMCtmbar 
kMCte 
kMCtebar 
kMCNC 
kDataEvent 
kCosmicEvent 
kBadOscType 
kNumOscTypes 
kOscTypeNoneSet 

Definition at line 244 of file Constants.h.

Enumerator
kNuMuSelection 
kNuESelection 
kNuESelectionLowPID 
kNuESelectionMidPID 
kNuESelectionHighPID 
kNuESelectionPeripheral 
kNuMuSelectionQ1 
kNuMuSelectionQ2 
kNuMuSelectionQ3 
kNuMuSelectionQ4 
kNCSelection 
kUnknownSelection 

Definition at line 100 of file Constants.h.

Enumerator
kSystWgtOnly 
kOscWgtOnly 
kSystAndOscWgt 
kAllWgt 

Definition at line 32 of file ShifterAndWeighter.h.

Function Documentation

static long cmf::BeamTypeKey ( cmf::BeamType_t const &  bt)
static
static long cmf::DetectorBeamSelectionTypesToKey ( novadaq::cnv::DetId const &  det,
cmf::BeamType_t const &  bt,
cmf::SelectionType_t const &  sel 
)
inlinestatic

Definition at line 42 of file StaticFuncs.h.

References BeamTypeKey(), DetectorTypeKey(), and SelectionTypeKey().

Referenced by cmf::CovarianceBinUtility::CovarianceBinUtility(), cmf::CalibSystHistMaker::CreateHistograms(), cmf::MetaData::DetectorBeamSelectionKey(), cmf::CovarianceFitHelper::InitializeEventLists(), InteractionTypeKey(), and cmf::CovarianceFitHelper::MakeDataMCCanv().

45 {
47 }
static long SelectionTypeKey(cmf::SelectionType_t const &st)
Definition: StaticFuncs.h:20
static long DetectorTypeKey(novadaq::cnv::DetId const &det)
Definition: StaticFuncs.h:16
static long BeamTypeKey(cmf::BeamType_t const &bt)
Definition: StaticFuncs.h:19
static long cmf::DetectorTypeKey ( novadaq::cnv::DetId const &  det)
static
static long cmf::EpochKey ( int const &  ep)
static
int cmf::EpochStringToInt ( std::string const &  epochStr)

Definition at line 72 of file Structs.cxx.

References LOG_DEBUG, and LOG_WARNING.

Referenced by cmf::MetaData::MetaData(), and cmf::MetaDataLite::MetaDataLite().

72  {
73  size_t first = epochStr.find_first_of("abcdefghijklmnopqrstuvwxyz");
74  if(first != std::string::npos){
75  LOG_WARNING("MetaData")
76  << "use of subperiods is deprecated, only full periods are used";
77 
78  //auto alpha = epochStr.substr(first, first + 1);
79  //return (std::stol(epochStr.substr(0, first).c_str(), nullptr, 0) * 1000 + static_cast<int>(*(alpha.c_str())));
80  }
81 
82  LOG_DEBUG("MetaData")
83  << epochStr
84  << " "
85  << std::stol(epochStr.substr(0, 1).c_str(), nullptr, 0);
86 
87  return std::stol(epochStr.substr(0, 1).c_str(), nullptr, 0);
88  }
#define LOG_DEBUG(stream)
Definition: Messenger.h:149
#define LOG_WARNING(category)
static long cmf::FileTypeKey ( cmf::FileType_t const &  ft)
static
double cmf::FitFunction_ForROOT ( const double *  params)

Definition at line 48 of file CovarianceFitHelper.cxx.

References cmf::CovarianceFitHelper::AddIteration(), cmf::CovarianceFitHelper::ChiSqr(), cmf::CovarianceFitHelper::FillMCSpectrum(), MECModelEnuComparisons::i, cmf::CovarianceFitHelper::Instance(), cmf::ShifterAndWeighter::Instance(), novadaq::cnv::kFARDET, novadaq::cnv::kNEARDET, LOG_DEBUG, LOG_VERBATIM, make_pair(), cet::pow(), cmf::ShifterAndWeighter::SetCurrentVals(), and std::sin().

Referenced by cmf::CovarianceFitHelper::InitializeMinimizer().

49  {
50  // debug
51  LOG_VERBATIM("CovarianceFit") << "Processing guess: ";
52 
53  for(size_t i = 0; i < gFFVals->size(); ++i){
54  if((*gFFNames)[i] == "Th23")
55  LOG_VERBATIM("CovarianceFit")
56  << "\t"
57  << (*gFFNames)[i]
58  << " : "
59  << params[i]
60  << "; sin^2(Th23) : "
61  << std::pow(std::sin(params[i]), 2.);
62  else if((*gFFNames)[i] == "dCP")
63  LOG_VERBATIM("CovarianceFit")
64  << "\t"
65  << (*gFFNames)[i]
66  << " : "
67  << params[i]
68  << " = "
69  << params[i] / gPI
70  << " (pi)";
71  else
72  LOG_VERBATIM("CovarianceFit")
73  << "\t"
74  << (*gFFNames)[i]
75  << " : "
76  << params[i];
77 
78  }
79 
80  // Grab the parameter space location vector to edit in the InputPoint
81  ParameterDetMap pParameterSpaceLocations = gInputPoint->ParameterSpaceLocation();
82 
83  // Start by looping through condensed name (key) vector
84  // We are only adjusting the parameter space values, don't bother
85  // with the systematics
86 
87  size_t keys_size = gFFNames->size();
88  for(size_t i = 0; i < keys_size; ++i) {
89 
92 
93  LOG_DEBUG("CovarianceMatrixFitter")
94  << "fill parameter space location "
95  << i
96  << " "
97  << gFFNames->at(i)
98  << " "
99  << params[i];
100 
101  if(pParameterSpaceLocations.count(tempPairFar) > 0){
102  pParameterSpaceLocations[tempPairFar] = params[i];
103  }
104  if(pParameterSpaceLocations.count(tempPairNear) > 0){
105  pParameterSpaceLocations[tempPairNear] = params[i];
106  }
107  }
108 
109  // Make the new input point and then set the ShifterAndWeighter to use it
110  gInputPoint = std::make_unique<cmf::InputPoint>(pParameterSpaceLocations,
111  gInputPoint->SystematicPulls(),
112  gInputPoint->FixedParameters(),
113  gInputPoint->FixedSystematics());
114 
115  // Now update the ShifterAndWeighter to check the goodness of fit
117 
118  LOG_DEBUG("CovarianceMatrixFitter")
119  << "fill the MC spectrum with the new values";
120 
121  // fill the MC spectrum
123 
124  LOG_DEBUG("CovarianceMatrixFitter")
125  << "get the new chi^2";
126 
127  // Get Chisq
129  gChiSqrCalc);
130 
131  LOG_DEBUG("CovarianceMatrixFitter")
132  << "add to the iterations";
133 
134  // store this iteration of the fitter
136  *gInputPoint,
137  gChiSqrCalc);
138 
139  LOG_VERBATIM("CovarianceFit")
140  << "\t"
141  << "FitValue = "
142  << chisq;
143 
144  // Return the fit value; that's what we're looking to minimize
145  return chisq;
146 
147  }
void FillMCSpectrum(cmf::EventListMap const &eventLists)
#define LOG_DEBUG(stream)
Definition: Messenger.h:149
void SetCurrentVals(cmf::InputPoint const &curPoint)
constexpr T pow(T x)
Definition: pow.h:75
void AddIteration(double const &chi2, cmf::InputPoint const &curPoint, cmf::ChiSqrCalc_t const &chiSqrCalc=cmf::kCovMat)
std::pair< Spectrum *, CheatDecomp * > make_pair(SpectrumLoaderBase &loader_data, SpectrumLoaderBase &loader_mc, HistAxis *axis, Cut *cut, const SystShifts &shift, const Var &wei)
Definition: DataMCLoad.C:335
double ChiSqr(cmf::InputPoint const &curPoint, cmf::ChiSqrCalc_t const &chiSqrType=cmf::kCovMat)
Far Detector at Ash River, MN.
std::pair< std::string, novadaq::cnv::DetId > ParameterDet
Definition: Constants.h:38
std::unordered_map< ParameterDet, double, ParameterDetHasher > ParameterDetMap
Definition: Constants.h:63
Near Detector in the NuMI cavern.
static ShifterAndWeighter * Instance()
static CovarianceFitHelper * Instance()
T sin(T number)
Definition: d0nt_math.hpp:132
std::unique_ptr< cmf::EventListMap > gMCEventLists
static const double gPI
std::unique_ptr< std::vector< double > > gFFVals
#define LOG_VERBATIM(category)
cmf::ChiSqrCalc_t gChiSqrCalc
std::unique_ptr< cmf::InputPoint > gInputPoint
std::unique_ptr< std::vector< std::string > > gFFNames
static long cmf::InteractionTypeKey ( cmf::InteractionType_t const &  it)
static
bool cmf::IsDataVar ( unsigned char const &  key)
inline

Definition at line 526 of file VarVals.h.

References kFake_Weight.

Referenced by IsDataVar().

527  {
528  return (key > 0 && key <= cmf::kFake_Weight);
529  }
static const unsigned char kFake_Weight
Definition: VarVals.h:38
bool cmf::IsDataVar ( std::string const &  name)
inline

Definition at line 533 of file VarVals.h.

References IsDataVar(), and VarNameToKey().

534  {
536  }
const XML_Char * name
Definition: expat.h:151
bool IsDataVar(unsigned char const &key)
Definition: VarVals.h:526
unsigned char VarNameToKey(std::string const &name)
Definition: VarVals.h:192
bool cmf::IsMCVar ( unsigned char const &  key)
inline

Definition at line 513 of file VarVals.h.

References kFake_Weight.

Referenced by IsMCVar().

514  {
515  return (key > cmf::kFake_Weight);
516  }
static const unsigned char kFake_Weight
Definition: VarVals.h:38
bool cmf::IsMCVar ( std::string const &  name)
inline

Definition at line 520 of file VarVals.h.

References IsMCVar(), and VarNameToKey().

521  {
523  }
const XML_Char * name
Definition: expat.h:151
bool IsMCVar(unsigned char const &key)
Definition: VarVals.h:513
unsigned char VarNameToKey(std::string const &name)
Definition: VarVals.h:192
bool cmf::isSumSmallXSecJoint2017Key ( unsigned char const &  varkey)
inline

Definition at line 150 of file VarVals.h.

References kTheta_Delta2Npi.

151  {
152  return varkey == kSumSmallXSecJoint2017 ||
153  varkey == kMaNCEL ||
154  varkey == kEtaNCEL ||
155  varkey == kVecCCQE_shape ||
156  varkey == kMvNCRES ||
157  varkey == kMaCOHpi ||
158  varkey == kR0COHpi ||
159  varkey == kAhtBY ||
160  varkey == kBhtBY ||
161  varkey == kCV1uBY ||
162  varkey == kCV2uBY ||
163  varkey == kNormDISCC ||
164  varkey == kRnubarnuCC ||
165  varkey == kDISNuclMod ||
166  varkey == kAGKY_xF1pi ||
167  varkey == kAGKY_pT1pi ||
168  varkey == kFormZone ||
169  varkey == kMFP_pi ||
170  varkey == kMFP_N ||
171  varkey == kFrCEx_pi ||
172  varkey == kFrElas_pi ||
173  varkey == kFrInel_pi ||
174  varkey == kFrAbs_pi ||
175  varkey == kFrPiProd_pi ||
176  varkey == kFrCEx_N ||
177  varkey == kFrElas_N ||
178  varkey == kFrInel_N ||
179  varkey == kFrAbs_N ||
180  varkey == kFrPiProd_N ||
181  varkey == kCCQEPauliSupViaKF ||
182  varkey == kCCQEMomDistroFGtoSF ||
183  varkey == kBR1gamma ||
184  varkey == kBR1eta ||
185  varkey == kTheta_Delta2Npi;
186 
187  }
static const unsigned char kEtaNCEL
Definition: VarVals.h:86
static const unsigned char kTheta_Delta2Npi
Definition: VarVals.h:111
static const unsigned char kRnubarnuCC
Definition: VarVals.h:96
static const unsigned char kMaCOHpi
Definition: VarVals.h:89
static const unsigned char kFrElas_N
Definition: VarVals.h:79
static const unsigned char kCCQEMomDistroFGtoSF
Definition: VarVals.h:108
static const unsigned char kFrPiProd_N
Definition: VarVals.h:107
static const unsigned char kFrElas_pi
Definition: VarVals.h:103
static const unsigned char kMvNCRES
Definition: VarVals.h:88
static const unsigned char kVecCCQE_shape
Definition: VarVals.h:87
static const unsigned char kBR1gamma
Definition: VarVals.h:109
static const unsigned char kFormZone
Definition: VarVals.h:76
static const unsigned char kFrAbs_N
Definition: VarVals.h:80
static const unsigned char kMaNCEL
Definition: VarVals.h:85
static const unsigned char kBhtBY
Definition: VarVals.h:92
static const unsigned char kNormDISCC
Definition: VarVals.h:95
static const unsigned char kDISNuclMod
Definition: VarVals.h:97
static const unsigned char kAGKY_pT1pi
Definition: VarVals.h:99
static const unsigned char kCV1uBY
Definition: VarVals.h:93
static const unsigned char kCCQEPauliSupViaKF
Definition: VarVals.h:81
static const unsigned char kFrAbs_pi
Definition: VarVals.h:104
static const unsigned char kR0COHpi
Definition: VarVals.h:90
static const unsigned char kFrPiProd_pi
Definition: VarVals.h:105
static const unsigned char kSumSmallXSecJoint2017
Definition: VarVals.h:84
static const unsigned char kFrInel_pi
Definition: VarVals.h:77
static const unsigned char kMFP_N
Definition: VarVals.h:101
static const unsigned char kBR1eta
Definition: VarVals.h:110
static const unsigned char kCV2uBY
Definition: VarVals.h:94
static const unsigned char kMFP_pi
Definition: VarVals.h:100
static const unsigned char kAhtBY
Definition: VarVals.h:91
static const unsigned char kAGKY_xF1pi
Definition: VarVals.h:98
static const unsigned char kFrCEx_N
Definition: VarVals.h:78
static const unsigned char kFrCEx_pi
Definition: VarVals.h:102
static const unsigned char kFrInel_N
Definition: VarVals.h:106
static cmf::BeamType_t cmf::KeyToBeamType ( long const &  key)
inlinestatic

Definition at line 116 of file StaticFuncs.h.

References DetectorTypeKey(), EpochKey(), mkDefs::fileType, FileTypeKey(), KeyToDetectorType(), KeyToEpoch(), KeyToFileType(), novadaq::cnv::kUNKNOWN_DET, kUnknownBeam, kUnknownFileType, and LOG_DEBUG.

Referenced by InteractionTypeKey(), KeyToInteractionType(), KeyToSelectionType(), and cmf::CovarianceFitHelper::MakeEnergySpectraFromBins().

117 {
118  auto detType = cmf::KeyToDetectorType(key);
120  auto epoch = cmf::KeyToEpoch(key);
121 
122  long tmpKey = key - cmf::EpochKey(epoch);
123 
124  if(detType != novadaq::cnv::kUNKNOWN_DET) tmpKey -= cmf::DetectorTypeKey(detType);
126 
127  if(tmpKey < 10000) return cmf::kUnknownBeam;
128 
129  LOG_DEBUG("StaticFuncs")
130  << "KeyToBeam: "
131  << key
132  << " "
133  << detType
134  << "/"
135  << cmf::DetectorTypeKey(detType)
136  << " "
137  << fileType
138  << "/"
140  << " "
141  << tmpKey
142  << " "
143  << (tmpKey / 10000) - 1;
144 
145 
146  return (cmf::BeamType_t)((tmpKey / 10000) - 1);
147 }
#define LOG_DEBUG(stream)
Definition: Messenger.h:149
static novadaq::cnv::DetId KeyToDetectorType(long const &key)
Definition: StaticFuncs.h:50
fileType
Definition: mkDefs.py:112
static long DetectorTypeKey(novadaq::cnv::DetId const &det)
Definition: StaticFuncs.h:16
static cmf::FileType_t KeyToFileType(long const &key)
Definition: StaticFuncs.h:89
static int KeyToEpoch(long const &key)
Definition: StaticFuncs.h:64
static long FileTypeKey(cmf::FileType_t const &ft)
Definition: StaticFuncs.h:18
enum cmf::beam_type BeamType_t
static long EpochKey(int const &ep)
Definition: StaticFuncs.h:17
static novadaq::cnv::DetId cmf::KeyToDetectorType ( long const &  key)
inlinestatic

Definition at line 50 of file StaticFuncs.h.

References novadaq::cnv::kUNKNOWN_DET, and LOG_DEBUG.

Referenced by InteractionTypeKey(), KeyToBeamType(), KeyToEpoch(), KeyToFileType(), KeyToInteractionType(), KeyToSelectionType(), cmf::CovarianceFitHelper::MakeEnergySpectraFromBins(), and cmf::CovarianceBinUtility::SelectionHighEdges().

51 {
52  if(key < 1000000000) return novadaq::cnv::kUNKNOWN_DET;
53 
54  LOG_DEBUG("StaticFuncs")
55  << "KeyToDetector: "
56  << key
57  << " "
58  << (key / 1000000000);
59 
60  return (novadaq::cnv::DetId)(key / 1000000000);
61 }
#define LOG_DEBUG(stream)
Definition: Messenger.h:149
static int cmf::KeyToEpoch ( long const &  key)
inlinestatic

Definition at line 64 of file StaticFuncs.h.

References DetectorTypeKey(), findDuplicateFiles::key, KeyToDetectorType(), novadaq::cnv::kUNKNOWN_DET, and LOG_DEBUG.

Referenced by InteractionTypeKey(), KeyToBeamType(), KeyToFileType(), KeyToInteractionType(), and KeyToSelectionType().

65 {
66  auto detType = cmf::KeyToDetectorType(key);
67 
68  long tmpKey = key;
69  if(detType != novadaq::cnv::kUNKNOWN_DET) tmpKey -= cmf::DetectorTypeKey(detType);
70 
71  if(tmpKey < 10000000) return 0;
72 
73  LOG_DEBUG("StaticFuncs")
74  << "KeyToEpoch: "
75  << key
76  << " "
77  << detType
78  << "/"
79  << cmf::DetectorTypeKey(detType)
80  << " "
81  << tmpKey
82  << " "
83  << (tmpKey / 10000000);
84 
85  return (tmpKey / 10000000);
86 }
#define LOG_DEBUG(stream)
Definition: Messenger.h:149
static novadaq::cnv::DetId KeyToDetectorType(long const &key)
Definition: StaticFuncs.h:50
static long DetectorTypeKey(novadaq::cnv::DetId const &det)
Definition: StaticFuncs.h:16
static cmf::FileType_t cmf::KeyToFileType ( long const &  key)
inlinestatic

Definition at line 89 of file StaticFuncs.h.

References DetectorTypeKey(), EpochKey(), KeyToDetectorType(), KeyToEpoch(), novadaq::cnv::kUNKNOWN_DET, kUnknownFileType, and LOG_DEBUG.

Referenced by InteractionTypeKey(), KeyToBeamType(), KeyToInteractionType(), and KeyToSelectionType().

90 {
91  auto detType = cmf::KeyToDetectorType(key);
92  auto epoch = cmf::KeyToEpoch(key);
93 
94  // if the epoch is not part of the key, KeyToEpoch returns 0
95  long tmpKey = key - cmf::EpochKey(epoch);
96  if(detType != novadaq::cnv::kUNKNOWN_DET) tmpKey -= cmf::DetectorTypeKey(detType);
97 
98  if(tmpKey < 100000) return cmf::kUnknownFileType;
99 
100  LOG_DEBUG("StaticFuncs")
101  << "KeyToFile: "
102  << key
103  << " "
104  << detType
105  << "/"
106  << cmf::DetectorTypeKey(detType)
107  << " "
108  << tmpKey
109  << " "
110  << (tmpKey / 100000) - 1;
111 
112  return (cmf::FileType_t)((tmpKey / 100000) - 1);
113 }
#define LOG_DEBUG(stream)
Definition: Messenger.h:149
static novadaq::cnv::DetId KeyToDetectorType(long const &key)
Definition: StaticFuncs.h:50
static long DetectorTypeKey(novadaq::cnv::DetId const &det)
Definition: StaticFuncs.h:16
static int KeyToEpoch(long const &key)
Definition: StaticFuncs.h:64
static long EpochKey(int const &ep)
Definition: StaticFuncs.h:17
enum cmf::file_type FileType_t
static cmf::InteractionType_t cmf::KeyToInteractionType ( long const &  key)
inlinestatic

Definition at line 188 of file StaticFuncs.h.

References BeamTypeKey(), DetectorTypeKey(), EpochKey(), mkDefs::fileType, FileTypeKey(), KeyToBeamType(), KeyToDetectorType(), KeyToEpoch(), KeyToFileType(), KeyToSelectionType(), novadaq::cnv::kUNKNOWN_DET, kUnknownBeam, kUnknownFileType, kUnknownInteraction, kUnknownSelection, and SelectionTypeKey().

Referenced by InteractionTypeKey().

189 {
190  auto detType = cmf::KeyToDetectorType(key);
191  auto beamType = cmf::KeyToBeamType(key);
193  auto selType = cmf::KeyToSelectionType(key);
194  auto epoch = cmf::KeyToEpoch(key);
195 
196  long tmpKey = key - cmf::EpochKey(epoch);
197  if(detType != novadaq::cnv::kUNKNOWN_DET) tmpKey -= cmf::DetectorTypeKey(detType);
198  if(beamType != cmf::kUnknownBeam ) tmpKey -= cmf::BeamTypeKey(beamType);
200  if(selType != cmf::kUnknownSelection ) tmpKey -= cmf::SelectionTypeKey(selType);
201 
202  if(tmpKey < 1) return cmf::kUnknownInteraction;
203 
204  return (cmf::InteractionType_t)(tmpKey);
205 }
static novadaq::cnv::DetId KeyToDetectorType(long const &key)
Definition: StaticFuncs.h:50
enum cmf::interaction_type InteractionType_t
static long SelectionTypeKey(cmf::SelectionType_t const &st)
Definition: StaticFuncs.h:20
static cmf::BeamType_t KeyToBeamType(long const &key)
Definition: StaticFuncs.h:116
fileType
Definition: mkDefs.py:112
static long DetectorTypeKey(novadaq::cnv::DetId const &det)
Definition: StaticFuncs.h:16
static cmf::FileType_t KeyToFileType(long const &key)
Definition: StaticFuncs.h:89
static int KeyToEpoch(long const &key)
Definition: StaticFuncs.h:64
static long FileTypeKey(cmf::FileType_t const &ft)
Definition: StaticFuncs.h:18
static cmf::SelectionType_t KeyToSelectionType(long const &key)
Definition: StaticFuncs.h:150
static long BeamTypeKey(cmf::BeamType_t const &bt)
Definition: StaticFuncs.h:19
static long EpochKey(int const &ep)
Definition: StaticFuncs.h:17
std::string cmf::KeyToLatex ( unsigned char const &  varkey)
inline

Definition at line 502 of file VarVals.h.

References KeyToVarName(), kNu_RecoE, and latex().

503  {
505 
506  if(varkey == cmf::kNu_RecoE) latex = "E_{#nu} (GeV)";
507 
508  return latex;
509 
510  }
static const unsigned char kNu_RecoE
Definition: VarVals.h:35
::xsd::cxx::tree::string< char, simple_type > string
Definition: Database.h:154
std::string KeyToVarName(unsigned char const &key)
Definition: VarVals.h:354
void latex(double x, double y, std::string txt, double ang=0, int align=22, double size=2/30.)
static cmf::SelectionType_t cmf::KeyToSelectionType ( long const &  key)
inlinestatic

Definition at line 150 of file StaticFuncs.h.

References BeamTypeKey(), DetectorTypeKey(), EpochKey(), mkDefs::fileType, FileTypeKey(), KeyToBeamType(), KeyToDetectorType(), KeyToEpoch(), KeyToFileType(), novadaq::cnv::kUNKNOWN_DET, kUnknownBeam, kUnknownFileType, kUnknownSelection, and LOG_DEBUG.

Referenced by InteractionTypeKey(), KeyToInteractionType(), cmf::CovarianceFitHelper::MakeEnergySpectraFromBins(), cmf::CovarianceFitHelper::MakeResultPlots(), and cmf::CovarianceBinUtility::SelectionHighEdges().

151 {
152  auto detType = cmf::KeyToDetectorType(key);
153  auto beamType = cmf::KeyToBeamType(key);
155  auto epoch = cmf::KeyToEpoch(key);
156 
157  long tmpKey = key - cmf::EpochKey(epoch);
158  if(detType != novadaq::cnv::kUNKNOWN_DET) tmpKey -= cmf::DetectorTypeKey(detType);
159  if(beamType != cmf::kUnknownBeam ) tmpKey -= cmf::BeamTypeKey(beamType);
161 
162  if(tmpKey < 100) return cmf::kUnknownSelection;
163 
164  LOG_DEBUG("StaticFuncs")
165  << "KeyToSelection: "
166  << key
167  << " "
168  << detType
169  << "/"
170  << cmf::DetectorTypeKey(detType)
171  << " "
172  << beamType
173  << "/"
174  << cmf::BeamTypeKey(beamType)
175  << " "
176  << fileType
177  << "/"
179  << " "
180  << tmpKey
181  << " "
182  << (tmpKey / 100) - 1;
183 
184  return (cmf::SelectionType_t)((tmpKey / 100) - 1);
185 }
#define LOG_DEBUG(stream)
Definition: Messenger.h:149
static novadaq::cnv::DetId KeyToDetectorType(long const &key)
Definition: StaticFuncs.h:50
static cmf::BeamType_t KeyToBeamType(long const &key)
Definition: StaticFuncs.h:116
enum cmf::sel_type SelectionType_t
fileType
Definition: mkDefs.py:112
static long DetectorTypeKey(novadaq::cnv::DetId const &det)
Definition: StaticFuncs.h:16
static cmf::FileType_t KeyToFileType(long const &key)
Definition: StaticFuncs.h:89
static int KeyToEpoch(long const &key)
Definition: StaticFuncs.h:64
static long FileTypeKey(cmf::FileType_t const &ft)
Definition: StaticFuncs.h:18
static long BeamTypeKey(cmf::BeamType_t const &bt)
Definition: StaticFuncs.h:19
static long EpochKey(int const &ep)
Definition: StaticFuncs.h:17
std::string cmf::KeyToVarName ( unsigned char const &  key)
inline

Definition at line 354 of file VarVals.h.

Referenced by cmf::EventListManipulator::FillEventList(), cmf::EventListManipulator::FillHistogram(), KeyToLatex(), operator<<(), SerializeEvents(), cmf::DataVarVals::set_val_at(), and cmf::ShifterAndWeighter::TotalShift().

355  {
356  // else if(key == kNuE_LID ) return "NuE_LID";
357  // else if(key == kNuE_LEM ) return "NuE_LEM";
358  if (key == kNuE_CVN ) return "NuE_CVN";
359  else if(key == kNuE_NumMichel ) return "NuE_NumMichel";
360  else if(key == kHad_RecoE ) return "Had_RecoE";
361  else if(key == kHad_RecoE_UBL ) return "Had_RecoE_UBL";
362  else if(key == kLep_RecoE ) return "Lep_RecoE";
363  else if(key == kLep_RecoE_UBL ) return "Lep_RecoE_UBL";
364  else if(key == kNu_RecoE ) return "Nu_RecoE";
365  else if(key == kNu_RecoE_UBL ) return "Nu_RecoE_UBL";
366  else if(key == kNu_RecoE_Quantiles_NuMu ) return "Nu_RecoE_Quantiles_NuMu";
367  else if(key == kNu_RecoE_NuE_Peripheral ) return "Nu_RecoE_NuE_Peripheral";
368  else if(key == kLep_RecoE_MCFrac ) return "Lep_RecoE_MCFrac";
369  else if(key == kRecoQ2 ) return "Reco_Q2";
370  else if(key == kFake_Weight ) return "Fake_Weight";
371 
372  else if(key == kTrueE ) return "TrueE";
373  else if(key == kTruePDG ) return "TruePDG";
374  else if(key == kTrueCCNC ) return "TrueCCNC";
375  else if(key == kTrueIntType ) return "TrueIntType";
376  else if(key == kTruePDGOrig ) return "TruePDGOrig";
377  else if(key == kTrueParentPDG ) return "TrueParentPDG";
378  else if(key == kTrueParentPT ) return "TrueParentPT";
379  else if(key == kTrueParentPZ ) return "TrueParentPZ";
380  else if(key == kTrueParentDecay ) return "TrueParentDecay";
381  else if(key == kTrueParentTargetPDG) return "TrueParentTargetPDG";
382  else if(key == kTrueHitNuc ) return "TrueHitNuc";
383  else if(key == kTrueIntMode ) return "TrueIntMode";
384 
385  else if(key == kXSecCVPPFX_Weight ) return "XsecCVWgt";
386  else if(key == kRPACCQE_Weight ) return "RPACCQEWgt";
387  else if(key == kRPARES_Weight ) return "RPARESWgt";
388  else if(key == kDISvnCC1pi_Weight ) return "DISvnCC1piWgt";
389  else if(key == kEmpiricalMEC_Weight ) return "EmpiricalMECWgt";
390  else if(key == kEmpiricalMECtoGENIEQE_Weight ) return "EmpiricalMECtoGENIEQEWgt";
391  else if(key == kEmpiricalMECtoGENIERES_Weight) return "EmpiricalMECtoGENIERESWgt";
392  else if(key == kOtherDIS_Weight ) return "OtherDISWgt";
393 
394  else if(key == kMinus2Sigma ) return "-2sigma";
395  else if(key == kMinus1Sigma ) return "-1sigma";
396  else if(key == kPlus1Sigma ) return "+1sigma";
397  else if(key == kPlus2Sigma ) return "+2sigma";
398 
399  //GENIE knobs
400  //else if( isSumSmallXSecJoint2017Key(key) ) return "SumSmallXSecJoint2017";
401 
402  // NC elastic parameters
403  else if(key == kMaNCEL ) return "MaNCEL";
404  else if(key == kEtaNCEL ) return "EtaNCEL";
405 
406  // CCQE parameters
407  //else if(key == kNormCCQE ) return "NormCCQE";
408  //else if(key == kNormCCQE_EnuDep ) return "NormCCQE_EnuDep";
409  //else if(key == kMaCCQE_shape ) return "MaCCQE_shape";
410  else if(key == kMaCCQE ) return "MaCCQE";
411  else if(key == kVecCCQE_shape ) return "VecCCQE_shape";
412 
413  // Resonant production: CC
414  //else if(key == kNormCCRES ) return "NormCCRES";
415  //else if(key == kMaCCRES_shape ) return "MaCCRES_shape";
416  //else if(key == kMvCCRES_shape ) return "MvCCRES_shape";
417  else if(key == kMaCCRES ) return "MaCCRES";
418  else if(key == kMvCCRES ) return "MvCCRES";
419 
420  // Resonant production: NC
421  //else if(key == kNormNCRES ) return "NormNCRES";
422  //else if(key == kMaNCRES_shape ) return "MaNCRES_shape";
423  //else if(key == kMvNCRES_shape ) return "MvNCRES_shape";
424  else if(key == kMaNCRES ) return "MaNCRES";
425  else if(key == kMvNCRES ) return "MvNCRES";
426 
427  // Coherent pion production
428  else if(key == kMaCOHpi ) return "MaCOHpi";
429  else if(key == kR0COHpi ) return "R0COHpi";
430 
431  // Non-resonant background
432  //else if(key == kRvpCC1pi ) return "RvpCC1pi";
433  //else if(key == kRvpCC2pi ) return "RvpCC2pi";
434  //else if(key == kRvpNC1pi ) return "RvpNC1pi";
435  //else if(key == kRvpNC2pi ) return "RvpNC2pi";
436  //else if(key == kRvnCC1pi ) return "RvnCC1pi";
437  //else if(key == kRvnCC2pi ) return "RvnCC2pi";
438  //else if(key == kRvnNC1pi ) return "RvnNC1pi";
439  //else if(key == kRvnNC2pi ) return "RvnNC2pi";
440  //else if(key == kRvbarpCC1pi ) return "RvbarpCC1pi";
441  //else if(key == kRvbarpCC2pi ) return "RvbarpCC2pi";
442  //else if(key == kRvbarpNC1pi ) return "RvbarpNC1pi";
443  //else if(key == kRvbarpNC2pi ) return "RvbarpNC2pi";
444  //else if(key == kRvbarnCC1pi ) return "RvbarnCC1pi";
445  //else if(key == kRvbarnCC2pi ) return "RvbarnCC2pi";
446  //else if(key == kRvbarnNC1pi ) return "RvbarnNC1pi";
447  //else if(key == kRvbarnNC2pi ) return "RvbarnNC2pi";
448 
449  // DIS tweaking parameters
450  else if(key == kAhtBY ) return "AhtBY";
451  else if(key == kBhtBY ) return "BhtBY";
452  else if(key == kCV1uBY ) return "CV1uBY";
453  else if(key == kCV2uBY ) return "CV2uBY";
454  //else if(key == kAhtBYshape ) return "AhtBYshape";
455  //else if(key == kBhtBYshape ) return "BhtBYshape";
456  //else if(key == kCV1uBYshape ) return "CV1uBYshape";
457  //else if(key == kCV2uBYshape ) return "CV2uBYshape";
458  else if(key == kNormDISCC ) return "NormDISCC";
459  else if(key == kRnubarnuCC ) return "RnubarnuCC";
460  else if(key == kDISNuclMod ) return "DISNuclMod";
461 
462  //else if(key == kNCGENIE ) return "NC";
463 
464  // Hadronization (free nucleon target)
465  else if(key == kAGKY_xF1pi ) return "AGKY_xF1pi";
466  else if(key == kAGKY_pT1pi ) return "AGKY_pT1pi";
467 
468  // Medium-effects to hadronization
469  else if(key == kFormZone ) return "FormZone";
470 
471  // Intranuclear rescattering parameters.
472  else if(key == kMFP_pi ) return "MFP_pi";
473  else if(key == kMFP_N ) return "MFP_N";
474  else if(key == kFrCEx_pi ) return "FrCEx_pi";
475  else if(key == kFrElas_pi ) return "FrElas_pi";
476  else if(key == kFrInel_pi ) return "FrInel_pi";
477  else if(key == kFrAbs_pi ) return "FrAbs_pi";
478  else if(key == kFrPiProd_pi ) return "FrPiProd_pi";
479  else if(key == kFrCEx_N ) return "FrCEx_N";
480  else if(key == kFrElas_N ) return "FrElas_N";
481  else if(key == kFrInel_N ) return "FrInel_N";
482  else if(key == kFrAbs_N ) return "FrAbs_N";
483  else if(key == kFrPiProd_N ) return "FrPiProd_N";
484 
485  // Nuclear model
486  else if(key == kCCQEPauliSupViaKF ) return "CCQEPauliSupViaKF";
487  else if(key == kCCQEMomDistroFGtoSF) return "CCQEMomDistroFGtoSF";
488 
489  // Resonance decays
490  else if(key == kBR1gamma ) return "BR1gamma";
491  else if(key == kBR1eta ) return "BR1eta";
492  else if(key == kTheta_Delta2Npi ) return "Theta_Delta2Npi";
493 
494  // 2017 RPA shape systematics
495  else if(key == kRPACCQEshapeEnh ) return "RPACCQEshapeEnh";
496  else if(key == kRPACCQEshapeSupp ) return "RPACCQEshapeSupp";
497 
498  return "Unknown";
499  }
static const unsigned char kNu_RecoE
Definition: VarVals.h:36
static const unsigned char kEtaNCEL
Definition: VarVals.h:86
static const unsigned char kNuE_CVN
Definition: VarVals.h:32
static const unsigned char kTheta_Delta2Npi
Definition: VarVals.h:111
static const unsigned char kRnubarnuCC
Definition: VarVals.h:96
static const unsigned char kMaCOHpi
Definition: VarVals.h:89
static const unsigned char kFrElas_N
Definition: VarVals.h:79
static const unsigned char kMinus1Sigma
Definition: VarVals.h:66
static const unsigned char kEmpiricalMEC_Weight
Definition: VarVals.h:60
static const unsigned char kCCQEMomDistroFGtoSF
Definition: VarVals.h:108
static const unsigned char kTrueParentPT
Definition: VarVals.h:48
static const unsigned char kOtherDIS_Weight
Definition: VarVals.h:63
static const unsigned char kFrPiProd_N
Definition: VarVals.h:107
static const unsigned char kNu_RecoE_NuE_Peripheral
Definition: VarVals.h:119
static const unsigned char kMvCCRES
Definition: VarVals.h:73
static const unsigned char kDISvnCC1pi_Weight
Definition: VarVals.h:59
static const unsigned char kMaNCRES
Definition: VarVals.h:72
static const unsigned char kEmpiricalMECtoGENIERES_Weight
Definition: VarVals.h:62
static const unsigned char kFrElas_pi
Definition: VarVals.h:103
static const unsigned char kMvNCRES
Definition: VarVals.h:88
static const unsigned char kVecCCQE_shape
Definition: VarVals.h:87
static const unsigned char kRecoQ2
Definition: VarVals.h:38
static const unsigned char kRPACCQEshapeSupp
Definition: VarVals.h:75
static const unsigned char kXSecCVPPFX_Weight
Definition: VarVals.h:56
static const unsigned char kNu_RecoE_Quantiles_NuMu
Definition: VarVals.h:118
static const unsigned char kTrueHitNuc
Definition: VarVals.h:52
static const unsigned char kBR1gamma
Definition: VarVals.h:109
static const unsigned char kFormZone
Definition: VarVals.h:76
static const unsigned char kFake_Weight
Definition: VarVals.h:39
static const unsigned char kFrAbs_N
Definition: VarVals.h:80
static const unsigned char kTrueParentPDG
Definition: VarVals.h:47
static const unsigned char kMaNCEL
Definition: VarVals.h:85
static const unsigned char kRPARES_Weight
Definition: VarVals.h:58
static const unsigned char kBhtBY
Definition: VarVals.h:92
static const unsigned char kLep_RecoE_UBL
Definition: VarVals.h:115
static const unsigned char kTruePDGOrig
Definition: VarVals.h:46
static const unsigned char kNormDISCC
Definition: VarVals.h:95
static const unsigned char kTrueCCNC
Definition: VarVals.h:44
static const unsigned char kDISNuclMod
Definition: VarVals.h:97
static const unsigned char kAGKY_pT1pi
Definition: VarVals.h:99
static const unsigned char kPlus1Sigma
Definition: VarVals.h:67
static const unsigned char kCV1uBY
Definition: VarVals.h:93
static const unsigned char kLep_RecoE_MCFrac
Definition: VarVals.h:37
static const unsigned char kMaCCQE
Definition: VarVals.h:70
static const unsigned char kNuE_NumMichel
Definition: VarVals.h:33
static const unsigned char kCCQEPauliSupViaKF
Definition: VarVals.h:81
static const unsigned char kTrueIntType
Definition: VarVals.h:45
static const unsigned char kFrAbs_pi
Definition: VarVals.h:104
static const unsigned char kTrueParentDecay
Definition: VarVals.h:50
static const unsigned char kR0COHpi
Definition: VarVals.h:90
static const unsigned char kFrPiProd_pi
Definition: VarVals.h:105
static const unsigned char kTrueParentTargetPDG
Definition: VarVals.h:51
static const unsigned char kFrInel_pi
Definition: VarVals.h:77
static const unsigned char kMFP_N
Definition: VarVals.h:101
static const unsigned char kBR1eta
Definition: VarVals.h:110
static const unsigned char kEmpiricalMECtoGENIEQE_Weight
Definition: VarVals.h:61
static const unsigned char kCV2uBY
Definition: VarVals.h:94
static const unsigned char kRPACCQE_Weight
Definition: VarVals.h:57
static const unsigned char kMFP_pi
Definition: VarVals.h:100
static const unsigned char kMinus2Sigma
Definition: VarVals.h:65
static const unsigned char kTrueE
Definition: VarVals.h:42
static const unsigned char kPlus2Sigma
Definition: VarVals.h:68
static const unsigned char kAhtBY
Definition: VarVals.h:91
static const unsigned char kAGKY_xF1pi
Definition: VarVals.h:98
static const unsigned char kHad_RecoE
Definition: VarVals.h:34
static const unsigned char kTrueIntMode
Definition: VarVals.h:53
static const unsigned char kFrCEx_N
Definition: VarVals.h:78
static const unsigned char kFrCEx_pi
Definition: VarVals.h:102
static const unsigned char kRPACCQEshapeEnh
Definition: VarVals.h:74
static const unsigned char kMaCCRES
Definition: VarVals.h:71
static const unsigned char kLep_RecoE
Definition: VarVals.h:35
static const unsigned char kNu_RecoE_UBL
Definition: VarVals.h:117
static const unsigned char kTrueParentPZ
Definition: VarVals.h:49
static const unsigned char kHad_RecoE_UBL
Definition: VarVals.h:116
static const unsigned char kFrInel_N
Definition: VarVals.h:106
static const unsigned char kTruePDG
Definition: VarVals.h:43
std::ostream& cmf::operator<< ( std::ostream &  o,
cmf::Event const &  event 
)

Definition at line 27 of file Event.cxx.

References cmf::Event::DataVals(), cmf::Event::DataValsFilled(), allTimeWatchdog::endl, cmf::Event::EventID(), cmf::Event::EventIDFilled(), cmf::Event::MCVals(), cmf::Event::MCValsFilled(), and msf_helper::o.

Referenced by cmf::DataVarVals::DataVarVals(), and cmf::MCVarVals::MCVarVals().

28  {
29  o << "Event Variables are: \n";
30 
31  if(event.EventIDFilled() ) o << event.EventID() ;
32  if(event.DataValsFilled()) o << event.DataVals();
33  if(event.MCValsFilled() ) o << event.MCVals() ;
34 
35  o << std::endl;
36 
37  return o;
38  }
std::ostream& cmf::operator<< ( std::ostream &  o,
cmf::SelectOutput const &  so 
)

Definition at line 80 of file CMFEventListToTextFile_module.cc.

References cmf::SelectOutput::cvn, cmf::SelectOutput::cycle, cmf::SelectOutput::energy, cmf::SelectOutput::event, msf_helper::o, cmf::SelectOutput::run, cmf::SelectOutput::selName, cmf::SelectOutput::slice, cmf::SelectOutput::subrun, and cmf::SelectOutput::wgt.

81  {
82  o << so.selName
83  << " "
84  << so.run
85  << " "
86  << so.subrun
87  << " "
88  << so.event
89  << " "
90  << so.slice
91  << " "
92  << so.cycle
93  << " "
94  << so.energy
95  << " "
96  << so.cvn
97  << " "
98  << so.wgt;
99 
100  return o;
101  }
std::ostream& cmf::operator<< ( std::ostream &  o,
SpillSummary const &  ss 
)
inline

Definition at line 153 of file Structs.h.

References cmf::SpillSummary::goodPOT, cmf::SpillSummary::liveTime, cmf::SpillSummary::numGoodSpills, msf_helper::o, cmf::SpillSummary::totalNumSpills, and cmf::SpillSummary::totalPOT.

154  {
155  o << "good POT: "
156  << ss.goodPOT
157  << " total POT: "
158  << ss.totalPOT
159  << " live time: "
160  << ss.liveTime
161  << " good spills: "
162  << ss.numGoodSpills
163  << " total spills: "
164  << ss.totalNumSpills;
165 
166  return o;
167  }
Float_t ss
Definition: plot.C:24
std::ostream& cmf::operator<< ( std::ostream &  o,
cmf::PointResult const &  pr 
)

Definition at line 159 of file PointResult.cxx.

References cmf::PointResult::ChiSqr(), allTimeWatchdog::endl, novadaq::cnv::kFARDET, novadaq::cnv::kNEARDET, msf_helper::o, cmf::PointResult::ParameterSpaceLocation(), cet::pow(), art::right(), std::sin(), and cmf::PointResult::SystematicPulls().

160  {
161  o << "PointResult:" << std::endl
162  << "\tChiSqr = " << std::setw(5) << std::right << pr.ChiSqr() << std::endl;
163 
164  o << "PointResult FD:" << std::endl;
165  for(auto itr : pr.ParameterSpaceLocation(novadaq::cnv::kFARDET) ){
166  if(itr.first=="dCP")
167  o << "\t" << std::setw(20) << itr.first << " " << itr.second << " = " << itr.second/TMath::Pi() << "(pi)" << std::endl;
168  else if(itr.first=="Th23")
169  o << "\t" << std::setw(20) << itr.first << " " << itr.second << ": sin^2(Th23) : " << std::pow(std::sin(itr.second), 2.) << std::endl;
170  else o << "\t" << std::setw(20) << itr.first << " " << itr.second << std::endl;
171  }
172  for(auto itr : pr.SystematicPulls(novadaq::cnv::kFARDET) ){
173  o << "\t" << std::setw(20) << itr.first << " " << itr.second << std::endl;
174  }
175 
176  o << "PointResult ND:" << std::endl;
177  for(auto itr : pr.ParameterSpaceLocation(novadaq::cnv::kNEARDET) ){
178  if(itr.first=="dCP")
179  o << "\t" << std::setw(20) << itr.first << " " << itr.second << " = " << itr.second/TMath::Pi() << "(pi)" << std::endl;
180  else if(itr.first=="Th23")
181  o << "\t" << std::setw(20) << itr.first << " " << itr.second << ": sin^2(Th23) : " << std::pow(std::sin(itr.second), 2.) << std::endl;
182  else o << "\t" << std::setw(20) << itr.first << " " << itr.second << std::endl;
183  }
184  for(auto itr : pr.SystematicPulls(novadaq::cnv::kNEARDET) ){
185  o << "\t" << std::setw(20) << itr.first << " " << itr.second << std::endl;
186  }
187 
188  return o;
189  }
constexpr auto const & right(const_AssnsIter< L, R, D, Dir > const &a, const_AssnsIter< L, R, D, Dir > const &b)
Definition: AssnsIter.h:112
constexpr T pow(T x)
Definition: pow.h:75
Far Detector at Ash River, MN.
Near Detector in the NuMI cavern.
T sin(T number)
Definition: d0nt_math.hpp:132
std::ostream& cmf::operator<< ( std::ostream &  o,
cmf::MCVarVals const &  mcvv 
)

Definition at line 235 of file VarVals.cxx.

References allTimeWatchdog::endl, kDISvnCC1pi_Weight, kEmpiricalMEC_Weight, kEmpiricalMECtoGENIEQE_Weight, kEmpiricalMECtoGENIERES_Weight, KeyToVarName(), kOtherDIS_Weight, kRPACCQE_Weight, kRPARES_Weight, kTrueCCNC, kTrueE, kTrueHitNuc, kTrueIntMode, kTrueIntType, kTrueParentDecay, kTrueParentPDG, kTrueParentPT, kTrueParentPZ, kTrueParentTargetPDG, kTruePDG, kTruePDGOrig, kXSecCVPPFX_Weight, msf_helper::o, and cmf::MCVarVals::val_at().

236  {
237  o << "Event MC related variables are: "
238  << "\n\t" << cmf::KeyToVarName(cmf::kTrueE) << " " << mcvv.val_at(cmf::kTrueE )
239  << "\n\t" << cmf::KeyToVarName(cmf::kTruePDG) << " " << mcvv.val_at(cmf::kTruePDG )
240  << "\n\t" << cmf::KeyToVarName(cmf::kTrueCCNC) << " " << mcvv.val_at(cmf::kTrueCCNC )
241  << "\n\t" << cmf::KeyToVarName(cmf::kTrueIntType) << " " << mcvv.val_at(cmf::kTrueIntType )
242  << "\n\t" << cmf::KeyToVarName(cmf::kTruePDGOrig) << " " << mcvv.val_at(cmf::kTruePDGOrig )
243  << "\n\t" << cmf::KeyToVarName(cmf::kTrueParentPDG) << " " << mcvv.val_at(cmf::kTrueParentPDG )
244  << "\n\t" << cmf::KeyToVarName(cmf::kTrueParentPT) << " " << mcvv.val_at(cmf::kTrueParentPT )
245  << "\n\t" << cmf::KeyToVarName(cmf::kTrueParentPZ) << " " << mcvv.val_at(cmf::kTrueParentPZ )
246  << "\n\t" << cmf::KeyToVarName(cmf::kTrueParentDecay) << " " << mcvv.val_at(cmf::kTrueParentDecay )
247  << "\n\t" << cmf::KeyToVarName(cmf::kTrueParentTargetPDG) << " " << mcvv.val_at(cmf::kTrueParentTargetPDG)
248  << "\n\t" << cmf::KeyToVarName(cmf::kTrueHitNuc) << " " << mcvv.val_at(cmf::kTrueHitNuc )
249  << "\n\t" << cmf::KeyToVarName(cmf::kTrueIntMode) << " " << mcvv.val_at(cmf::kTrueIntMode )
250  << "\n\t" << cmf::KeyToVarName(cmf::kXSecCVPPFX_Weight) << " " << mcvv.val_at(cmf::kXSecCVPPFX_Weight )
251  << "\n\t" << cmf::KeyToVarName(cmf::kRPACCQE_Weight) << " " << mcvv.val_at(cmf::kRPACCQE_Weight )
252  << "\n\t" << cmf::KeyToVarName(cmf::kRPARES_Weight) << " " << mcvv.val_at(cmf::kRPARES_Weight )
253  << "\n\t" << cmf::KeyToVarName(cmf::kDISvnCC1pi_Weight) << " " << mcvv.val_at(cmf::kDISvnCC1pi_Weight )
254  << "\n\t" << cmf::KeyToVarName(cmf::kEmpiricalMEC_Weight) << " " << mcvv.val_at(cmf::kEmpiricalMEC_Weight)
257  << "\n\t" << cmf::KeyToVarName(cmf::kOtherDIS_Weight) << " " << mcvv.val_at(cmf::kOtherDIS_Weight )
258  << "\n\t and too many GENIE knobs to list"
259  << std::endl;
260 
261  return o;
262  }
static const unsigned char kTruePDGOrig
Definition: VarVals.h:45
static const unsigned char kTruePDG
Definition: VarVals.h:42
static const unsigned char kEmpiricalMEC_Weight
Definition: VarVals.h:59
static const unsigned char kXSecCVPPFX_Weight
Definition: VarVals.h:55
static const unsigned char kRPACCQE_Weight
Definition: VarVals.h:56
static const unsigned char kTrueE
Definition: VarVals.h:41
static const unsigned char kDISvnCC1pi_Weight
Definition: VarVals.h:58
static const unsigned char kTrueParentPZ
Definition: VarVals.h:48
static const unsigned char kTrueParentPDG
Definition: VarVals.h:46
static const unsigned char kTrueIntType
Definition: VarVals.h:44
static const unsigned char kTrueIntMode
Definition: VarVals.h:52
static const unsigned char kTrueParentPT
Definition: VarVals.h:47
static const unsigned char kTrueParentTargetPDG
Definition: VarVals.h:50
static const unsigned char kOtherDIS_Weight
Definition: VarVals.h:62
static const unsigned char kEmpiricalMECtoGENIEQE_Weight
Definition: VarVals.h:60
std::string KeyToVarName(unsigned char const &key)
Definition: VarVals.h:354
static const unsigned char kRPARES_Weight
Definition: VarVals.h:57
static const unsigned char kTrueCCNC
Definition: VarVals.h:43
static const unsigned char kEmpiricalMECtoGENIERES_Weight
Definition: VarVals.h:61
static const unsigned char kTrueParentDecay
Definition: VarVals.h:49
static const unsigned char kTrueHitNuc
Definition: VarVals.h:51
std::ostream& cmf::operator<< ( std::ostream &  o,
cmf::EventId const &  evid 
)

Definition at line 295 of file VarVals.cxx.

References cmf::EventId::cycle, cmf::EventId::event, msf_helper::o, cmf::EventId::run, cmf::EventId::slice, and cmf::EventId::subrun.

296  {
297  o << "run "
298  << evid.run
299  << "\tsubrun "
300  << evid.subrun
301  << "\tevent "
302  << evid.event
303  << "\tcycle "
304  << evid.cycle
305  << "\tslice "
306  << evid.slice;
307 
308  return o;
309  }
std::ostream& cmf::operator<< ( std::ostream &  o,
cmf::DataVarVals const &  dvv 
)

Definition at line 426 of file VarVals.cxx.

References allTimeWatchdog::endl, cmf::MetaData::isMC, KeyToVarName(), kFake_Weight, kHad_RecoE, kLep_RecoE, kLep_RecoE_MCFrac, kNu_RecoE, kNuE_CVN, kNuE_NumMichel, kNuESelection, kNuMuSelection, msf_helper::o, cmf::MetaData::selectionType, and cmf::DataVarVals::val_at().

427  {
428  cmf::MetaData numumd;
429  cmf::MetaData nuemd;
430 
431  numumd.isMC = false;
432  nuemd .isMC = false;
435 
436  o << "Event Data related variables are: "
437  << "\n\t" << cmf::KeyToVarName(cmf::kNuE_CVN) << " " << dvv.val_at(cmf::kNuE_CVN, nuemd )
438  << "\n\t" << cmf::KeyToVarName(cmf::kNuE_NumMichel) << " " << dvv.val_at(cmf::kNuE_NumMichel, nuemd )
439  << "\n\t" << cmf::KeyToVarName(cmf::kHad_RecoE) << " " << dvv.val_at(cmf::kHad_RecoE, nuemd )
440  << "\n\t" << cmf::KeyToVarName(cmf::kLep_RecoE) << " " << dvv.val_at(cmf::kLep_RecoE, nuemd )
441  << "\n\t" << cmf::KeyToVarName(cmf::kLep_RecoE_MCFrac) << " " << dvv.val_at(cmf::kLep_RecoE_MCFrac, numumd)
442  << "\n\t" << cmf::KeyToVarName(cmf::kNu_RecoE ) << " nue " << dvv.val_at(cmf::kNu_RecoE, nuemd )
443  << "\n\t" << cmf::KeyToVarName(cmf::kNu_RecoE ) << " numu " << dvv.val_at(cmf::kNu_RecoE, numumd)
444  << "\n\t" << cmf::KeyToVarName(cmf::kFake_Weight) << " " << dvv.val_at(cmf::kFake_Weight, nuemd )
445  << std::endl;
446 
447  return o;
448  }
static const unsigned char kFake_Weight
Definition: VarVals.h:38
static const unsigned char kNu_RecoE
Definition: VarVals.h:35
bool isMC
Definition: Structs.h:64
static const unsigned char kNuE_CVN
Definition: VarVals.h:31
cmf::SelectionType_t selectionType
Definition: Structs.h:67
static const unsigned char kLep_RecoE
Definition: VarVals.h:34
static const unsigned char kNuE_NumMichel
Definition: VarVals.h:32
static const unsigned char kLep_RecoE_MCFrac
Definition: VarVals.h:36
std::string KeyToVarName(unsigned char const &key)
Definition: VarVals.h:354
static const unsigned char kHad_RecoE
Definition: VarVals.h:33
std::ostream& cmf::operator<< ( std::ostream &  o,
cmf::InputPoint const &  ip 
)

Definition at line 463 of file InputPoint.cxx.

References allTimeWatchdog::endl, cmf::InputPoint::FixedParameters(), cmf::InputPoint::FixedSystematics(), novadaq::cnv::kFARDET, novadaq::cnv::kNEARDET, msf_helper::o, cmf::InputPoint::ParameterSpaceLocation(), cet::pow(), std::sin(), and cmf::InputPoint::SystematicPulls().

464  {
465 
466  o << "InputPoint ND:" << std::endl;
467  for(auto const& itr : ip.ParameterSpaceLocation(novadaq::cnv::kNEARDET) ){
468  if (itr.first=="dCP")
469  o << "\t" << std::setw(20) << itr.first << " " << std::setw(10) << itr.second << " = " << itr.second/TMath::Pi() << "(pi)" << std::endl;
470  else if(itr.first=="Th23")
471  o << "\t" << std::setw(20) << itr.first << " " << std::setw(10) << itr.second << ": sin^2(Th23) : " << std::pow(std::sin(itr.second), 2.) << std::endl;
472  else
473  o << "\t" << std::setw(20) << itr.first << " " << std::setw(10) << itr.second << std::endl;
474  }
475  for(auto const& itr : ip.SystematicPulls(novadaq::cnv::kNEARDET) ){
476  o << "\t"<< std::setw(10) << itr.first << " " << std::setw(10) << itr.second << std::endl;
477  }
478  o << "Fixed:" << std::endl;
479  for(auto const& itr: ip.FixedParameters(novadaq::cnv::kNEARDET) )
480  o << "\t"<< std::setw(10) << itr << " " << std::endl;
481  for(auto const& itr: ip.FixedSystematics(novadaq::cnv::kNEARDET) )
482  o << "\t"<< std::setw(10) << itr << " " << std::endl;
483 
484  o << "\nInputPoint FD:" << std::endl;
485  for(auto const& itr : ip.ParameterSpaceLocation(novadaq::cnv::kFARDET) ){
486  if (itr.first=="dCP")
487  o << "\t" << std::setw(20) << itr.first << " " << std::setw(10) << itr.second << " = " << itr.second/TMath::Pi() << "(pi)" << std::endl;
488  else if(itr.first=="Th23")
489  o << "\t" << std::setw(20) << itr.first << " " << std::setw(10) << itr.second << ": sin^2(Th23) : " << std::pow(std::sin(itr.second), 2.) << std::endl;
490  else
491  o << "\t"<< std::setw(10) << itr.first << " " << std::setw(10) << itr.second << std::endl;
492  }
493  for(auto const& itr : ip.SystematicPulls(novadaq::cnv::kFARDET) ){
494  o << "\t"<< std::setw(10) << itr.first
495  << " " << std::setw(10) << itr.second << std::endl;
496  }
497  o << "Fixed:" << std::endl;
498  for(auto const& itr: ip.FixedParameters(novadaq::cnv::kFARDET) )
499  o << std::setw(10)<< itr << " " << std::endl;
500  for(auto const& itr: ip.FixedSystematics(novadaq::cnv::kFARDET) )
501  o << std::setw(10)<< itr << " " << std::endl;
502 
503  return o;
504  }
constexpr T pow(T x)
Definition: pow.h:75
TString ip
Definition: loadincs.C:5
Far Detector at Ash River, MN.
Near Detector in the NuMI cavern.
T sin(T number)
Definition: d0nt_math.hpp:132
std::string cmf::ParameterDetToString ( cmf::ParameterDet const &  pd)
inline

Definition at line 412 of file Constants.h.

References novadaq::cnv::kNEARDET.

413 {
414  std::string retStr(pd.first);
415 
416  if (pd.second == novadaq::cnv::kNEARDET) retStr += "ND";
417  else if(pd.second == novadaq::cnv::kNEARDET) retStr += "FD";
418 
419  return retStr;
420 }
Near Detector in the NuMI cavern.
::xsd::cxx::tree::string< char, simple_type > string
Definition: Database.h:154
template<class T >
static bool cmf::readSingleResult ( art::Results const &  r,
std::string const &  module,
std::string const &  instance,
std::vector< T > &  resvec 
)
static
template<class T >
static bool cmf::readVectorResults ( art::Results const &  r,
std::string const &  module,
std::string const &  instance,
std::vector< T > &  resvec 
)
static
static long cmf::SelectionTypeKey ( cmf::SelectionType_t const &  st)
static
void cmf::SerializeEventListMap ( std::string const &  dirName,
cmf::EventListMap const &  listMap 
)

Definition at line 107 of file Event.cxx.

References dir, fetch_tb_beamline_files::md, art::TFileDirectory::mkdir(), SerializeEvents(), ss, and cmf::MetaData::ToString().

Referenced by cmf::EventList::AddEvent(), and cmf::CAFToEventList::endJob().

109  {
110 
111  // make a directory in the output file, and the act of making it puts us
112  // into that directory of the file
114  auto dir = tfs->mkdir(dirName.c_str());
115 
118 
119  // first the metadata tree
120  // leave this pointer dangling because the fileService created & manages it
121  TTree *metadataTree = dir.make<TTree>( "metadata", "Event list metadata" );
122  metadataTree->Branch("metadata", &md);
123  metadataTree->Branch("spillsummary", &ss);
124 
125  for(auto const& itr : listMap){
126 
127  md = itr.first;
128  ss = itr.second.ListSpillSummary();
129 
130  metadataTree->Fill();
131 
132  // same principle as above for the event tree
133  TTree *eventTree = dir.make<TTree>(md.ToString().c_str(), "CovarianceMatrixFit event records");
134 
135  cmf::SerializeEvents(eventTree,
136  md,
137  itr.second);
138 
139  } // end loop over the lists
140  }
TFileDirectory mkdir(std::string const &dir, std::string const &descr="")
Float_t ss
Definition: plot.C:24
void SerializeEvents(TTree *eventTree, cmf::MetaData const &md, cmf::EventList const &events)
Definition: Event.cxx:41
std::string ToString() const
Definition: Structs.cxx:138
TDirectory * dir
Definition: macro.C:5
void cmf::SerializeEvents ( TTree *  eventTree,
cmf::MetaData const &  md,
cmf::EventList const &  events 
)

Definition at line 41 of file Event.cxx.

References cmf::DataVars::DataVars(), cmf::MetaData::isMC, kCCQEPauliSupViaKF, findDuplicateFiles::key, KeyToVarName(), kMaCCQE, LOG_DEBUG, cmf::TruthVars::TruthVars(), and cmf::WeightVars::WeightVars().

Referenced by cmf::EventList::AddEvent(), and SerializeEventListMap().

44  {
45 
46  cmf::DataVars dataVars;
47  cmf::TruthVars truthVars;
48  cmf::WeightVars weightVars;
49  std::vector<cmf::GENIEVar> genieVars(cmf::kCCQEPauliSupViaKF);
50  cmf::EventId eventId;
51 
52  std::string eventidst("run/I:subrun/I:event/I:cycle/I:slice/I");
53  std::string datavarst("fNue_CVN/F:fNuE_NumMichel/F:fHad_RecoE/F:fLep_RecoE/F:fLep_RecoE_MCFrac/F:fRecoQ2/F:fFake_Weight/F");
54  std::string weightvst("fXSecCVPPFX_Weight/F:fRPACCQE_Weight/F:fRPARES_Weight/F:fDISvnCC1pi_Weight/F:fEmpiricalMEC_Weight/F:fEmpiricalMECtoGENIEQE_Weight/F:fEmpiricalMECtoGENIERES_Weight/F:fOtherDIS_Weight/F");
55  std::string truthvrst("fTrueE/F:fTruePDG/F:fTrueCCNC/F:fTrueIntType/F:fTruePDGOrig/F:fTrueParentPDG/F:fTrueNucHit/F:fTrueParentPT/F:fTrueParentPZ/F:fTrueParentDecay/F:fTrueParentTargetPDG/F:fTrueIntMode/F");
56  // truthvrst. append("fTrueParentPT/F:fTrueParentPZ/F:fTrueParentDecay/F:fTrueParentTargetPDG/F");
57  std::string genievrst("fMinus2Sigma/F:fMinus1Sigma/F:fPlus1Sigma/F:fPlus2Sigma/F");
58 
59  eventTree->Branch("eventId" , &eventId , eventidst.c_str());
60  eventTree->Branch("dataVars", &dataVars, datavarst.c_str());
61 
62  if(md.isMC){
63  eventTree->Branch("truthVars" , &truthVars , truthvrst.c_str());
64  eventTree->Branch("weightVars", &weightVars, weightvst.c_str());
65 
66  //for(auto key = cmf::kMaNCEL; key < cmf::kOtherGENIE + 1; ++key){
67  //for(auto key = cmf::kMaCCQE; key < cmf::kSumSmallXSecJoint2017 + 1; ++key){
68  for(auto key = cmf::kMaCCQE; key < cmf::kCCQEPauliSupViaKF + 1; ++key){
69  LOG_DEBUG("Event")
70  << "For branch "
72  << ", get genieVars["
73  << key - cmf::kMaCCQE
74  << "]";
75  eventTree->Branch(cmf::KeyToVarName(key).c_str(), &genieVars[key - cmf::kMaCCQE], genievrst.c_str());
76  }
77  } // end if MC
78 
79  for(auto & itr : events){
80  eventId = itr->EventID();
81  dataVars = itr->DataVals().DataVars();
82 
83  if(md.isMC){
84  truthVars = itr->MCVals().TruthVars();
85  weightVars = itr->MCVals().WeightVars();
86 
87  // for(auto key = cmf::kMaNCEL; key < cmf::kOtherGENIE + 1; ++key)
88  //for(auto key = cmf::kMaCCQE; key < cmf::kSumSmallXSecJoint2017 + 1; ++key){
89  for(auto key = cmf::kMaCCQE; key < cmf::kCCQEPauliSupViaKF + 1; ++key){
90  LOG_DEBUG("Event")
91  << "For branch "
93  << ", get genieVars["
94  << key - cmf::kMaCCQE
95  << "]";
96  genieVars[key - cmf::kMaCCQE] = itr->MCVals().GENIEVars().GENIEVarByKey(key);
97  }
98  }
99 
100  eventTree->Fill();
101  }
102  }
#define LOG_DEBUG(stream)
Definition: Messenger.h:149
static const unsigned char kCCQEPauliSupViaKF
Definition: VarVals.h:80
::xsd::cxx::tree::string< char, simple_type > string
Definition: Database.h:154
std::string KeyToVarName(unsigned char const &key)
Definition: VarVals.h:354
static const unsigned char kMaCCQE
Definition: VarVals.h:69
void events(int which)
Definition: Cana.C:52
static cmf::BeamType_t cmf::StringToBeamType ( std::string const &  str)
inlinestatic

Definition at line 223 of file StaticFuncs.h.

References kFHC, kRHC, and kUnknownBeam.

Referenced by cmf::CovarianceFitHelper::InitializeEventLists(), InteractionTypeKey(), and cmf::RandomUniverses::reconfigure().

224 {
225  if (str == "RHC") return cmf::kRHC;
226  else if(str == "FHC") return cmf::kFHC;
227 
228  return cmf::kUnknownBeam;
229 }
static cmf::ChiSqrCalc_t cmf::StringToChiSqrType ( std::string const &  str)
static
static novadaq::cnv::DetId cmf::StringToDetectorType ( std::string const &  str)
inlinestatic

Definition at line 232 of file StaticFuncs.h.

References kCNP, kCovMat, novadaq::cnv::kFARDET, kGauss, novadaq::cnv::kNEARDET, kPoisson, novadaq::cnv::kUNKNOWN_DET, kUnknownChiSqr, runNovaSAM::str, and StringToChiSqrType().

Referenced by cmf::CovarianceFitHelper::InitializeEventLists(), InteractionTypeKey(), and cmf::RandomUniverses::reconfigure().

233 {
234  if (str.find("Near") != std::string::npos) return novadaq::cnv::kNEARDET;
235  else if(str.find("Far") != std::string::npos) return novadaq::cnv::kFARDET;
236 
238 }
Far Detector at Ash River, MN.
Near Detector in the NuMI cavern.
long cmf::StringToKey ( std::string const &  str)

Definition at line 15 of file Structs.cxx.

References DetectorTypeKey(), EpochKey(), FileTypeKey(), InteractionTypeKey(), kBeam, kCosmicMuon, kDataFile, findDuplicateFiles::key, novadaq::cnv::kFARDET, kNC, novadaq::cnv::kNEARDET, kNuEBarCC, kNuECC, kNuESelection, kNuESelectionHighPID, kNuESelectionLowPID, kNuESelectionMidPID, kNuESelectionPeripheral, kNuMuBarCC, kNuMuCC, kNuMuSelection, kNuMuSelectionQ1, kNuMuSelectionQ2, kNuMuSelectionQ3, kNuMuSelectionQ4, kNuTauBarCC, kNuTauCC, kRockFluxSwap, kRockNonSwap, kSwap, kTauSwap, and SelectionTypeKey().

Referenced by cmf::EventListManipulator::EpochEventCap(), and cmf::EventListManipulator::reconfigure().

16  {
17  long key = 0;
18 
19  // check for the detector
20  if (str.find("Near") != std::string::npos) key += cmf::DetectorTypeKey(novadaq::cnv::kNEARDET);
21  else if(str.find("Far" ) != std::string::npos) key += cmf::DetectorTypeKey(novadaq::cnv::kFARDET );
22 
23  // check the epoch
24  if (str.find("P1") != std::string::npos) key += cmf::EpochKey(1);
25  else if(str.find("P2") != std::string::npos) key += cmf::EpochKey(2);
26  else if(str.find("P3") != std::string::npos) key += cmf::EpochKey(3);
27  else if(str.find("P4") != std::string::npos) key += cmf::EpochKey(4);
28  else if(str.find("P5") != std::string::npos) key += cmf::EpochKey(5);
29  else if(str.find("P6") != std::string::npos) key += cmf::EpochKey(6);
30  else if(str.find("P7") != std::string::npos) key += cmf::EpochKey(7);
31  else if(str.find("P8") != std::string::npos) key += cmf::EpochKey(8);
32  else if(str.find("P9") != std::string::npos) key += cmf::EpochKey(9);
33 
34  // check the file type
35  if (str.find("Beam" ) != std::string::npos) key += cmf::FileTypeKey(cmf::kBeam );
36  else if(str.find("TauSwap" ) != std::string::npos) key += cmf::FileTypeKey(cmf::kTauSwap );
37  else if(str.find("FluxSwap" ) != std::string::npos) key += cmf::FileTypeKey(cmf::kSwap );
38  else if(str.find("RockFluxSwap") != std::string::npos) key += cmf::FileTypeKey(cmf::kRockFluxSwap);
39  else if(str.find("RockNonSwap" ) != std::string::npos) key += cmf::FileTypeKey(cmf::kRockNonSwap );
40  else if(str.find("Data" ) != std::string::npos) key += cmf::FileTypeKey(cmf::kDataFile );
41 
42  // check for a selection
43  if (str.find("NuMuSelQ1" ) != std::string::npos) key += cmf::SelectionTypeKey(cmf::kNuMuSelectionQ1 );
44  else if(str.find("NuMuSelQ2" ) != std::string::npos) key += cmf::SelectionTypeKey(cmf::kNuMuSelectionQ2 );
45  else if(str.find("NuMuSelQ3" ) != std::string::npos) key += cmf::SelectionTypeKey(cmf::kNuMuSelectionQ3 );
46  else if(str.find("NuMuSelQ4" ) != std::string::npos) key += cmf::SelectionTypeKey(cmf::kNuMuSelectionQ4 );
47  else if(str.find("NuMuSel" ) != std::string::npos) key += cmf::SelectionTypeKey(cmf::kNuMuSelection );
48  else if(str.find("NuESel_AllPID" ) != std::string::npos) key += cmf::SelectionTypeKey(cmf::kNuESelection );
49  else if(str.find("NuESel_LowPID" ) != std::string::npos) key += cmf::SelectionTypeKey(cmf::kNuESelectionLowPID );
50  else if(str.find("NuESel_MidPID" ) != std::string::npos) key += cmf::SelectionTypeKey(cmf::kNuESelectionMidPID );
51  else if(str.find("NuESel_HighPID" ) != std::string::npos) key += cmf::SelectionTypeKey(cmf::kNuESelectionHighPID );
52  else if(str.find("NuESel_Peripheral") != std::string::npos) key += cmf::SelectionTypeKey(cmf::kNuESelectionPeripheral);
53 
54  // check for an interaction type
55  if (str.find("NuMuCC" ) != std::string::npos) key += cmf::InteractionTypeKey(cmf::kNuMuCC );
56  else if(str.find("NuMuBarCC" ) != std::string::npos) key += cmf::InteractionTypeKey(cmf::kNuMuBarCC );
57  else if(str.find("NuECC" ) != std::string::npos) key += cmf::InteractionTypeKey(cmf::kNuECC );
58  else if(str.find("NuEBarCC" ) != std::string::npos) key += cmf::InteractionTypeKey(cmf::kNuEBarCC );
59  else if(str.find("NuTauCC" ) != std::string::npos) key += cmf::InteractionTypeKey(cmf::kNuTauCC );
60  else if(str.find("NuTauBarCC") != std::string::npos) key += cmf::InteractionTypeKey(cmf::kNuTauBarCC);
61  else if(str.find("NC" ) != std::string::npos) key += cmf::InteractionTypeKey(cmf::kNC );
62  else if(str.find("CosmicMuon") != std::string::npos) key += cmf::InteractionTypeKey(cmf::kCosmicMuon);
63 
64  // check if this is data or MC, multiply by +/- 1 depending on the answer
65  if (str.find("MC") != std::string::npos) key *= 1;
66  else if(str.find("Data") != std::string::npos) key *= -1;
67 
68  return key;
69  }
static long SelectionTypeKey(cmf::SelectionType_t const &st)
Definition: StaticFuncs.h:20
Far Detector at Ash River, MN.
static long DetectorTypeKey(novadaq::cnv::DetId const &det)
Definition: StaticFuncs.h:16
Near Detector in the NuMI cavern.
static long InteractionTypeKey(cmf::InteractionType_t const &it)
Definition: StaticFuncs.h:21
static long FileTypeKey(cmf::FileType_t const &ft)
Definition: StaticFuncs.h:18
static long EpochKey(int const &ep)
Definition: StaticFuncs.h:17
static cmf::SelectionType_t cmf::StringToSelectionType ( std::string const &  str)
inlinestatic

Definition at line 208 of file StaticFuncs.h.

References kNCSelection, kNuESelectionHighPID, kNuESelectionLowPID, kNuESelectionPeripheral, kNuMuSelectionQ1, kNuMuSelectionQ2, kNuMuSelectionQ3, kNuMuSelectionQ4, and kUnknownSelection.

Referenced by cmf::CovarianceFitHelper::InitializeEventLists(), InteractionTypeKey(), and cmf::RandomUniverses::reconfigure().

209 {
210  if (str == "NuESel_LowPID" ) return cmf::kNuESelectionLowPID;
211  else if(str == "NuESel_HighPID" ) return cmf::kNuESelectionHighPID;
212  else if(str == "NuESel_Peripheral") return cmf::kNuESelectionPeripheral;
213  else if(str == "NuMuSelQ1" ) return cmf::kNuMuSelectionQ1;
214  else if(str == "NuMuSelQ2" ) return cmf::kNuMuSelectionQ2;
215  else if(str == "NuMuSelQ3" ) return cmf::kNuMuSelectionQ3;
216  else if(str == "NuMuSelQ4" ) return cmf::kNuMuSelectionQ4;
217  else if(str == "NCSel" ) return cmf::kNCSelection;
218 
219  return cmf::kUnknownSelection;
220 }
unsigned char cmf::VarNameToKey ( std::string const &  name)
inline

Definition at line 192 of file VarVals.h.

References kAGKY_pT1pi, kAGKY_xF1pi, kAhtBY, kBhtBY, kBR1eta, kBR1gamma, kCCQEMomDistroFGtoSF, kCCQEPauliSupViaKF, kCV1uBY, kCV2uBY, kDISNuclMod, kDISvnCC1pi_Weight, kEmpiricalMEC_Weight, kEmpiricalMECtoGENIEQE_Weight, kEmpiricalMECtoGENIERES_Weight, kEtaNCEL, kFake_Weight, kFormZone, kFrAbs_N, kFrAbs_pi, kFrCEx_N, kFrCEx_pi, kFrElas_N, kFrElas_pi, kFrInel_N, kFrInel_pi, kFrPiProd_N, kFrPiProd_pi, kHad_RecoE, kHad_RecoE_UBL, kLep_RecoE, kLep_RecoE_MCFrac, kLep_RecoE_UBL, kMaCCQE, kMaCCRES, kMaCOHpi, kMaNCEL, kMaNCRES, kMFP_N, kMFP_pi, kMinus1Sigma, kMinus2Sigma, kMvCCRES, kMvNCRES, kNormDISCC, kNu_RecoE, kNu_RecoE_NuE_Peripheral, kNu_RecoE_Quantiles_NuMu, kNu_RecoE_UBL, kNuE_CVN, kNuE_NumMichel, kOtherDIS_Weight, kPlus1Sigma, kPlus2Sigma, kR0COHpi, kRecoQ2, kRnubarnuCC, kRPACCQE_Weight, kRPACCQEshapeEnh, kRPACCQEshapeSupp, kRPARES_Weight, kTheta_Delta2Npi, kTrueCCNC, kTrueE, kTrueHitNuc, kTrueIntMode, kTrueIntType, kTrueParentDecay, kTrueParentPDG, kTrueParentPT, kTrueParentPZ, kTrueParentTargetPDG, kTruePDG, kTruePDGOrig, kVecCCQE_shape, kXSecCVPPFX_Weight, and LOG_DEBUG.

Referenced by cmf::MCVarVals::GENIEWgts(), IsDataVar(), IsMCVar(), cmf::MCVarVals::set_GENIEWgt(), cmf::MCVarVals::set_val_at(), cmf::DataVarVals::set_val_at(), cmf::MCVarVals::SetGENIEVarVal(), cmf::MCVarVals::val_at(), and cmf::DataVarVals::val_at().

193  {
194  if (name == "NuE_CVN" ) return kNuE_CVN;
195  else if(name == "NuE_NumMichel" ) return kNuE_NumMichel;
196  else if(name == "Had_RecoE" ) return kHad_RecoE;
197  else if(name == "Had_RecoE_UBL" ) return kHad_RecoE_UBL;
198  else if(name == "Lep_RecoE" ) return kLep_RecoE;
199  else if(name == "Lep_RecoE_UBL" ) return kLep_RecoE_UBL;
200  else if(name == "Nu_RecoE" ) return kNu_RecoE;
201  else if(name == "Nu_RecoE_UBL" ) return kNu_RecoE_UBL;
202  else if(name == "Nu_RecoE_Quantiles_NuMu" ) return kNu_RecoE_Quantiles_NuMu;
203  else if(name == "Nu_RecoE_NuE_Peripheral" ) return kNu_RecoE_NuE_Peripheral;
204  else if(name == "Lep_RecoE_MCFrac" ) return kLep_RecoE_MCFrac;
205  else if(name == "Reco_Q2" ) return kRecoQ2;
206  else if(name == "Fake_Weight" ) return kFake_Weight;
207 
208  else if(name == "TrueE" ) return kTrueE;
209  else if(name == "TruePDG" ) return kTruePDG;
210  else if(name == "TrueCCNC" ) return kTrueCCNC;
211  else if(name == "TrueIntType" ) return kTrueIntType;
212  else if(name == "TruePDGOrig" ) return kTruePDGOrig;
213  else if(name == "TrueParentPDG" ) return kTrueParentPDG;
214  else if(name == "TrueParentPT" ) return kTrueParentPT;
215  else if(name == "TrueParentPZ" ) return kTrueParentPZ;
216  else if(name == "TrueParentDecay" ) return kTrueParentDecay;
217  else if(name == "TrueParentTargetPDG" ) return kTrueParentTargetPDG;
218  else if(name == "TrueHitNuc" ) return kTrueHitNuc;
219  else if(name == "TrueIntMode" ) return kTrueIntMode;
220 
221  else if(name == "XsecCVWgt" ) return kXSecCVPPFX_Weight;
222  else if(name == "RPACCQEWgt" ) return kRPACCQE_Weight;
223  else if(name == "RPARESWgt" ) return kRPARES_Weight;
224  else if(name == "DISvnCC1piWgt" ) return kDISvnCC1pi_Weight;
225  else if(name == "EmpiricalMECWgt" ) return kEmpiricalMEC_Weight;
226  else if(name == "EmpiricalMECtoGENIEQEWgt" ) return kEmpiricalMECtoGENIEQE_Weight;
227  else if(name == "EmpiricalMECtoGENIERESWgt") return kEmpiricalMECtoGENIERES_Weight;
228  else if(name == "OtherDISWgt" ) return kOtherDIS_Weight;
229 
230  // the next two are for backwards compatibility with files using an
231  // old naming scheme for those variables
232  else if(name == "TuftsCC_Weight" ) return kXSecCVPPFX_Weight;
233  else if(name == "RPACCQE_Weight" ) return kRPACCQE_Weight;
234 
235  else if(name == "-2sigma" ) return kMinus2Sigma;
236  else if(name == "-1sigma" ) return kMinus1Sigma;
237  else if(name == "+1sigma" ) return kPlus1Sigma;
238  else if(name == "+2sigma" ) return kPlus2Sigma;
239 
240  //GENIE knobs
241  //else if( isSumSmallXSecJoint2017Name(name) ) return kSumSmallXSecJoint2017;
242  //else if( isUnusedGENIEName (name) ) return 0;
243 
244  // NC elastic parameters
245  else if(name == "MaNCEL" ) return kMaNCEL;
246  else if(name == "EtaNCEL" ) return kEtaNCEL;
247 
248  // CCQE parameters
249  //else if(name == "NormCCQE" ) return kNormCCQE;
250  //else if(name == "NormCCQE_EnuDep" ) return kNormCCQE_EnuDep;
251  //else if(name == "MaCCQE_shape" ) return kMaCCQE_shape;
252  else if(name == "MaCCQE" ) return kMaCCQE;
253  else if(name == "VecCCQE_shape" ) return kVecCCQE_shape;
254 
255  // Resonant production: CC
256  //else if(name == "NormCCRES" ) return kNormCCRES;
257  //else if(name == "MaCCRES_shape" ) return kMaCCRES_shape;
258  //else if(name == "MvCCRES_shape" ) return kMvCCRES_shape;
259  else if(name == "MaCCRES" ) return kMaCCRES;
260  else if(name == "MvCCRES" ) return kMvCCRES;
261 
262  // Resonant production: NC
263  //else if(name == "NormNCRES" ) return kNormNCRES;
264  //else if(name == "MaNCRES_shape" ) return kMaNCRES_shape;
265  //else if(name == "MvNCRES_shape" ) return kMvNCRES_shape;
266  else if(name == "MaNCRES" ) return kMaNCRES;
267  else if(name == "MvNCRES" ) return kMvNCRES;
268 
269  // Coherent pion production
270  else if(name == "MaCOHpi" ) return kMaCOHpi;
271  else if(name == "R0COHpi" ) return kR0COHpi;
272 
273  // Non-resonant background
274  //else if(name == "RvpCC1pi" ) return kRvpCC1pi;
275  //else if(name == "RvpCC2pi" ) return kRvpCC2pi;
276  //else if(name == "RvpNC1pi" ) return kRvpNC1pi;
277  //else if(name == "RvpNC2pi" ) return kRvpNC2pi;
278  //else if(name == "RvnCC1pi" ) return kRvnCC1pi;
279  //else if(name == "RvnCC2pi" ) return kRvnCC2pi;
280  //else if(name == "RvnNC1pi" ) return kRvnNC1pi;
281  //else if(name == "RvnNC2pi" ) return kRvnNC2pi;
282  //else if(name == "RvbarpCC1pi" ) return kRvbarpCC1pi;
283  //else if(name == "RvbarpCC2pi" ) return kRvbarpCC2pi;
284  //else if(name == "RvbarpNC1pi" ) return kRvbarpNC1pi;
285  //else if(name == "RvbarpNC2pi" ) return kRvbarpNC2pi;
286  //else if(name == "RvbarnCC1pi" ) return kRvbarnCC1pi;
287  //else if(name == "RvbarnCC2pi" ) return kRvbarnCC2pi;
288  //else if(name == "RvbarnNC1pi" ) return kRvbarnNC1pi;
289  //else if(name == "RvbarnNC2pi" ) return kRvbarnNC2pi;
290 
291  // DIS tweaking parameters
292  else if(name == "AhtBY" ) return kAhtBY;
293  else if(name == "BhtBY" ) return kBhtBY;
294  else if(name == "CV1uBY" ) return kCV1uBY;
295  else if(name == "CV2uBY" ) return kCV2uBY;
296  //else if(name == "AhtBYshape" ) return kAhtBYshape;
297  //else if(name == "BhtBYshape" ) return kBhtBYshape;
298  //else if(name == "CV1uBYshape" ) return kCV1uBYshape;
299  //else if(name == "CV2uBYshape" ) return kCV2uBYshape;
300  else if(name == "NormDISCC" ) return kNormDISCC;
301  else if(name == "RnubarnuCC" ) return kRnubarnuCC;
302  else if(name == "DISNuclMod" ) return kDISNuclMod;
303 
304  //else if(name == "NC" ) return kNCGENIE;
305 
306  // Hadronization (free nucleon target)
307  else if(name == "AGKY_xF1pi" ) return kAGKY_xF1pi;
308  else if(name == "AGKY_pT1pi" ) return kAGKY_pT1pi;
309 
310  // Medium-effects to hadronization
311  else if(name == "FormZone" ) return kFormZone;
312 
313  // Intranuclear rescattering parameters.
314  else if(name == "MFP_pi" ) return kMFP_pi;
315  else if(name == "MFP_N" ) return kMFP_N;
316  else if(name == "FrCEx_pi" ) return kFrCEx_pi;
317  else if(name == "FrElas_pi" ) return kFrElas_pi;
318  else if(name == "FrInel_pi" ) return kFrInel_pi;
319  else if(name == "FrAbs_pi" ) return kFrAbs_pi;
320  else if(name == "FrPiProd_pi" ) return kFrPiProd_pi;
321  else if(name == "FrCEx_N" ) return kFrCEx_N;
322  else if(name == "FrElas_N" ) return kFrElas_N;
323  else if(name == "FrInel_N" ) return kFrInel_N;
324  else if(name == "FrAbs_N" ) return kFrAbs_N;
325  else if(name == "FrPiProd_N" ) return kFrPiProd_N;
326 
327  // Nuclear model
328  else if(name == "CCQEPauliSupViaKF" ) return kCCQEPauliSupViaKF;
329  else if(name == "CCQEMomDistroFGtoSF" ) return kCCQEMomDistroFGtoSF;
330 
331  // Resonance decays
332  else if(name == "BR1gamma" ) return kBR1gamma;
333  else if(name == "BR1eta" ) return kBR1eta;
334  else if(name == "Theta_Delta2Npi" ) return kTheta_Delta2Npi;
335 
336  // 2017 RPA shape systematics
337  else if(name == "RPACCQEshapeEnh" ) return kRPACCQEshapeEnh;
338  else if(name == "RPACCQEshapeSupp" ) return kRPACCQEshapeSupp;
339 
340  // if we passed it the name Unknown intentionally, just return 0
341  else if(name == "Unknown" ) return 0;
342 
343  LOG_DEBUG("VarNameToKey")
344  << "Variable name "
345  << name
346  << " does not correspond to any key, returning 0";
347 
348  return 0;
349  }
#define LOG_DEBUG(stream)
Definition: Messenger.h:149
const XML_Char * name
Definition: expat.h:151
static const unsigned char kNu_RecoE
Definition: VarVals.h:36
static const unsigned char kEtaNCEL
Definition: VarVals.h:86
static const unsigned char kNuE_CVN
Definition: VarVals.h:32
static const unsigned char kTheta_Delta2Npi
Definition: VarVals.h:111
static const unsigned char kRnubarnuCC
Definition: VarVals.h:96
static const unsigned char kMaCOHpi
Definition: VarVals.h:89
static const unsigned char kFrElas_N
Definition: VarVals.h:79
static const unsigned char kMinus1Sigma
Definition: VarVals.h:66
static const unsigned char kEmpiricalMEC_Weight
Definition: VarVals.h:60
static const unsigned char kCCQEMomDistroFGtoSF
Definition: VarVals.h:108
static const unsigned char kTrueParentPT
Definition: VarVals.h:48
static const unsigned char kOtherDIS_Weight
Definition: VarVals.h:63
static const unsigned char kFrPiProd_N
Definition: VarVals.h:107
static const unsigned char kNu_RecoE_NuE_Peripheral
Definition: VarVals.h:119
static const unsigned char kMvCCRES
Definition: VarVals.h:73
static const unsigned char kDISvnCC1pi_Weight
Definition: VarVals.h:59
static const unsigned char kMaNCRES
Definition: VarVals.h:72
static const unsigned char kEmpiricalMECtoGENIERES_Weight
Definition: VarVals.h:62
static const unsigned char kFrElas_pi
Definition: VarVals.h:103
static const unsigned char kMvNCRES
Definition: VarVals.h:88
static const unsigned char kVecCCQE_shape
Definition: VarVals.h:87
static const unsigned char kRecoQ2
Definition: VarVals.h:38
static const unsigned char kRPACCQEshapeSupp
Definition: VarVals.h:75
static const unsigned char kXSecCVPPFX_Weight
Definition: VarVals.h:56
static const unsigned char kNu_RecoE_Quantiles_NuMu
Definition: VarVals.h:118
static const unsigned char kTrueHitNuc
Definition: VarVals.h:52
static const unsigned char kBR1gamma
Definition: VarVals.h:109
static const unsigned char kFormZone
Definition: VarVals.h:76
static const unsigned char kFake_Weight
Definition: VarVals.h:39
static const unsigned char kFrAbs_N
Definition: VarVals.h:80
static const unsigned char kTrueParentPDG
Definition: VarVals.h:47
static const unsigned char kMaNCEL
Definition: VarVals.h:85
static const unsigned char kRPARES_Weight
Definition: VarVals.h:58
static const unsigned char kBhtBY
Definition: VarVals.h:92
static const unsigned char kLep_RecoE_UBL
Definition: VarVals.h:115
static const unsigned char kTruePDGOrig
Definition: VarVals.h:46
static const unsigned char kNormDISCC
Definition: VarVals.h:95
static const unsigned char kTrueCCNC
Definition: VarVals.h:44
static const unsigned char kDISNuclMod
Definition: VarVals.h:97
static const unsigned char kAGKY_pT1pi
Definition: VarVals.h:99
static const unsigned char kPlus1Sigma
Definition: VarVals.h:67
static const unsigned char kCV1uBY
Definition: VarVals.h:93
static const unsigned char kLep_RecoE_MCFrac
Definition: VarVals.h:37
static const unsigned char kMaCCQE
Definition: VarVals.h:70
static const unsigned char kNuE_NumMichel
Definition: VarVals.h:33
static const unsigned char kCCQEPauliSupViaKF
Definition: VarVals.h:81
static const unsigned char kTrueIntType
Definition: VarVals.h:45
static const unsigned char kFrAbs_pi
Definition: VarVals.h:104
static const unsigned char kTrueParentDecay
Definition: VarVals.h:50
static const unsigned char kR0COHpi
Definition: VarVals.h:90
static const unsigned char kFrPiProd_pi
Definition: VarVals.h:105
static const unsigned char kTrueParentTargetPDG
Definition: VarVals.h:51
static const unsigned char kFrInel_pi
Definition: VarVals.h:77
static const unsigned char kMFP_N
Definition: VarVals.h:101
static const unsigned char kBR1eta
Definition: VarVals.h:110
static const unsigned char kEmpiricalMECtoGENIEQE_Weight
Definition: VarVals.h:61
static const unsigned char kCV2uBY
Definition: VarVals.h:94
static const unsigned char kRPACCQE_Weight
Definition: VarVals.h:57
static const unsigned char kMFP_pi
Definition: VarVals.h:100
static const unsigned char kMinus2Sigma
Definition: VarVals.h:65
static const unsigned char kTrueE
Definition: VarVals.h:42
static const unsigned char kPlus2Sigma
Definition: VarVals.h:68
static const unsigned char kAhtBY
Definition: VarVals.h:91
static const unsigned char kAGKY_xF1pi
Definition: VarVals.h:98
static const unsigned char kHad_RecoE
Definition: VarVals.h:34
static const unsigned char kTrueIntMode
Definition: VarVals.h:53
static const unsigned char kFrCEx_N
Definition: VarVals.h:78
static const unsigned char kFrCEx_pi
Definition: VarVals.h:102
static const unsigned char kRPACCQEshapeEnh
Definition: VarVals.h:74
static const unsigned char kMaCCRES
Definition: VarVals.h:71
static const unsigned char kLep_RecoE
Definition: VarVals.h:35
static const unsigned char kNu_RecoE_UBL
Definition: VarVals.h:117
static const unsigned char kTrueParentPZ
Definition: VarVals.h:49
static const unsigned char kHad_RecoE_UBL
Definition: VarVals.h:116
static const unsigned char kFrInel_N
Definition: VarVals.h:106
static const unsigned char kTruePDG
Definition: VarVals.h:43

Variable Documentation

const std::string cmf::cBeamType_LatexStrings[4]
Initial value:
= { "FHC",
"RHC",
"0HC"
"Unknown Beam" }

Definition at line 80 of file Constants.h.

Referenced by cmf::CovarianceFitHelper::Make1DSpectra(), cmf::CovarianceFitHelper::MakeEnergySpectraFromBins(), and cmf::MetaData::ToStringLatex().

const std::string cmf::cBeamType_Strings[4]
const std::string cmf::cChiSqrCalcType_Strings[5]
const std::string cmf::cFileTypeStrings[10]
Initial value:
= {"Beam",
"FluxSwap",
"TauSwap",
"Data",
"CosmicBackground",
"RockFluxSwap",
"RockNonSwap",
"num_file_types",
"bad_file_type",
"UnknownFileType"}

Definition at line 380 of file Constants.h.

Referenced by cmf::CovarianceMatrixMaker::beginJob(), cmf::MetaData::DetectorFilePeriodString(), cmf::CovarianceFitHelper::Make1DSpectra(), cmf::MetaData::OscType(), and cmf::MetaData::ToString().

const std::string cmf::cFileTypeStrings_Latex[10]
Initial value:
= {"Beam",
"Flux Swap",
"#tau Swap",
"Data",
"Cosmic #mu Background",
"Rock Flux Swap #mu Background",
"Rock Non Swap #mu Background",
"num_file_types",
"bad_file_type",
"none_set"}

Definition at line 391 of file Constants.h.

Referenced by cmf::MetaData::ToStringLatex().

const std::string cmf::cInteractionType_LatexStrings[10]
Initial value:
={"Unknown Interaction",
"#nu_{#mu} CC",
"#bar{#nu}_{#mu} CC",
"#nu_{e} CC",
"#bar{#nu}_{e} CC",
"#nu_{#tau} CC",
"#bar{#nu}_{#tau} CC",
"NC",
"Cosmic Muon",
"Rock Muon"}

Definition at line 157 of file Constants.h.

Referenced by cmf::MetaData::ToStringLatex().

const std::string cmf::cInteractionType_Strings[10]
Initial value:
={"UnknownInteraction",
"NuMuCC",
"NuMuBarCC",
"NuECC",
"NuEBarCC",
"NuTauCC",
"NuTauBarCC",
"NC",
"CosmicMuon",
"RockMuon"}

Definition at line 168 of file Constants.h.

Referenced by cmf::MetaData::OscType(), and cmf::MetaData::ToString().

const std::string cmf::cOscParams_LatexScales[kNumOscParams]
Initial value:
= {" (m)",
" (g/cc)",
" (#times10^{-5} eV^{2})",
" (#times10^{-3} eV^{2})",
"",
"",
"",
" (#pi rad)",
"",
"",
"(eV^{2})"}

Definition at line 207 of file Constants.h.

Referenced by cmf::PlotUtilities::Initialize().

const double cmf::cOscParams_Scales[kNumOscParams]
Initial value:
= {1.,
1.,
1.e5,
1.e3,
1.,
1.,
1.,
1. / 3.14159,
1.,
1.,
1.}

Definition at line 195 of file Constants.h.

Referenced by cmf::ContourMaker::AddPointToMap(), and cmf::PlotUtilities::Initialize().

const std::string cmf::cOscParams_Strings[kNumOscParams]
const std::string cmf::cOscParams_Strings_Latex[kNumOscParams]
Initial value:
= {"L",
"#rho",
"#Delta m^{2}_{21}",
"#Delta m^{2}_{32}",
"#theta_{12}",
"#theta_{13}",
"#theta_{23}",
"#delta_{CP}",
"#theta_{24}",
"#theta_{34}",
"#Delta m^{2}_{41}"}

Definition at line 231 of file Constants.h.

Referenced by cmf::PlotUtilities::Make1DPlot(), cmf::PlotUtilities::Make2DContours(), cmf::PlotUtilities::Make2DHist(), cmf::StarPlot::MakeStarPlot(), and cmf::PlotUtilities::StoreContourGraphs().

const int cmf::cOscType_Colors[20]
Initial value:
= {kGreen + 0, kGreen + 0,
kGreen + 2, kGreen + 2,
kGreen - 8, kGreen - 8,
kMagenta + 0, kMagenta + 0,
kMagenta + 2, kMagenta + 2,
kCyan + 0, kCyan + 0,
kCyan + 2, kCyan + 2,
kGray,
kBlack,
kYellow,
kPink,
kPink + 5,
kPink - 5}

Definition at line 299 of file Constants.h.

const int cmf::cOscType_FillStyles[20]
Initial value:
= {1001, 3001, 1001, 3001, 1001,3001,
1001, 3001, 1001, 3001,
1001, 3001, 1001, 3001,
1001,
1001,
1001,
1001,
0,
0}

Definition at line 313 of file Constants.h.

const std::string cmf::cOscType_LatexStrings[20]
Initial value:
={"MC #nu_{#mu} #rightarrow #nu_{#mu}",
"MC #bar{#nu}_{#mu} #rightarrow #bar{#nu}_{#mu}",
"MC #nu_{#mu} #rightarrow #nu_{e}",
"MC #bar{#nu}_{#mu} #rightarrow #bar{#nu}_{e}",
"MC #nu_{#mu} #rightarrow #nu_{#tau}",
"MC #bar{#nu}_{#mu} #rightarrow #bar{#nu}_{#tau}",
"MC #nu_{e} #rightarrow #nu_{e}",
"MC #bar{#nu}_{e} #rightarrow #bar{#nu}_{e}",
"MC #nu_{e} #rightarrow #nu_{#mu}",
"MC #bar{#nu}_{e} #rightarrow #bar{#nu}_{#mu}",
"MC #nu_{#tau} #rightarrow #nu_{#mu}",
"MC #bar{#nu}_{#tau} #rightarrow #bar{#nu}_{#mu}",
"MC #nu_{#tau} #rightarrow #nu_{e}",
"MC #bar{#nu}_{#tau} #rightarrow #bar{#nu}_{e}",
"MC NC",
"DATA",
"COSMIC",
"bad_osc_type",
"num_osc_types",
"none_set"}

Definition at line 278 of file Constants.h.

const int cmf::cOscType_Markers[20]
Initial value:
= {22, 22, 23, 23, 24, 24,
32, 32, 26, 26,
25, 25, 27, 27,
33,
33,
20,
21,
0,
0}

Definition at line 323 of file Constants.h.

const std::string cmf::cOscType_Strings[20]
Initial value:
= {"mc_mm", "mc_mmbar", "mc_me", "mc_mebar", "mc_mt", "mc_mtbar",
"mc_ee", "mc_eebar", "mc_em", "mc_embar",
"mc_tm", "mc_tmbar", "mc_te", "mc_tebar",
"mc_nc",
"data",
"cosmic",
"bad_osc_type",
"num_osc_types",
"none_set"}

Definition at line 268 of file Constants.h.

const std::string cmf::cSelectionType_LatexStrings[12]
Initial value:
= {"#nu_{#mu} Selection",
"#nu_{e} Selection",
"#nu_{e} Selection - Low PID",
"#nu_{e} Selection - Mid PID",
"#nu_{e} Selection - High PID",
"#nu_{e} Selection - Peripheral",
"#nu_{#mu} Selection - Quantile 1",
"#nu_{#mu} Selection - Quantile 2",
"#nu_{#mu} Selection - Quantile 3",
"#nu_{#mu} Selection - Quantile 4",
"NC Selection",
"Unknown Selection"}

Definition at line 128 of file Constants.h.

Referenced by cmf::CovarianceFitHelper::Make1DSpectra(), cmf::CovarianceFitHelper::MakeEnergySpectraFromBins(), and cmf::MetaData::ToStringLatex().

const std::string cmf::cSelectionType_Strings[12]
cmf::ChiSqrCalc_t cmf::gChiSqrCalc

Definition at line 41 of file CovarianceFitHelper.cxx.

CovarianceBinUtility* cmf::gCovBinUtil = nullptr
static

Definition at line 15 of file CovarianceBinUtility.cxx.

Referenced by cmf::CovarianceBinUtility::Instance().

CovarianceFitHelper* cmf::gCovFitHelp = nullptr
static

Definition at line 33 of file CovarianceFitHelper.cxx.

Referenced by cmf::CovarianceFitHelper::Instance().

ROOT::Math::Minimizer* cmf::gFFMin

Definition at line 40 of file CovarianceFitHelper.cxx.

std::unique_ptr<std::vector<std::string> > cmf::gFFNames

Definition at line 39 of file CovarianceFitHelper.cxx.

std::unique_ptr<std::vector<double> > cmf::gFFVals

Definition at line 38 of file CovarianceFitHelper.cxx.

std::unique_ptr<cmf::InputPoint> cmf::gInputPoint

Definition at line 37 of file CovarianceFitHelper.cxx.

Referenced by cmf::CovarianceFitHelper::Minimize().

std::unique_ptr<cmf::EventListMap> cmf::gMCEventLists

Definition at line 36 of file CovarianceFitHelper.cxx.

double cmf::gMinChiSqr
int cmf::gMinStatus

Definition at line 43 of file CovarianceFitHelper.cxx.

Referenced by cmf::CovarianceFitHelper::Result().

ParameterConstraints* cmf::gParConstInstance = nullptr
static

Definition at line 12 of file ParameterConstraints.cxx.

Referenced by cmf::ParameterConstraints::Instance().

const double cmf::gPI = 3.141592653589793238462643383279502884
static

Definition at line 44 of file CovarianceFitHelper.cxx.

PlotUtilities* cmf::gPlotUtil = nullptr
static

Definition at line 22 of file CMFPlotUtilities.cxx.

Referenced by cmf::PlotUtilities::Instance().

ShifterAndWeighter* cmf::gShifterAndWeighter = nullptr
static

Definition at line 451 of file ShifterAndWeighter.cxx.

Referenced by cmf::ShifterAndWeighter::Instance().

SystematicSigmas* cmf::gSystSigmaInstance = nullptr
static

Definition at line 11 of file SystematicSigmas.cxx.

NuisanceParameters* cmf::gSystSigmaInstance = nullptr
static
const unsigned char cmf::kAGKY_pT1pi = 60
static

Definition at line 98 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kAGKY_xF1pi = 59
static

Definition at line 97 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kAhtBY = 52
static

Definition at line 90 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kBhtBY = 53
static

Definition at line 91 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kBR1eta = 71
static

Definition at line 109 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kBR1gamma = 70
static

Definition at line 108 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kCCQEMomDistroFGtoSF = 69
static

Definition at line 107 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kCCQEPauliSupViaKF = 44
static
const unsigned char cmf::kCV1uBY = 54
static

Definition at line 92 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kCV2uBY = 55
static

Definition at line 93 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kDISNuclMod = 58
static

Definition at line 96 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kDISvnCC1pi_Weight = 24
static
const unsigned char cmf::kEmpiricalMEC_Weight = 25
static
const unsigned char cmf::kEmpiricalMECtoGENIEQE_Weight = 26
static
const unsigned char cmf::kEmpiricalMECtoGENIERES_Weight = 27
static
const unsigned char cmf::kEtaNCEL = 47
static

Definition at line 85 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kFake_Weight = 8
static
const double cmf::kFCAlpha[2] ={0.68, 0.90}

Definition at line 404 of file Constants.h.

const int cmf::kFCAlphaColors[2] = {4, 2}

Definition at line 405 of file Constants.h.

const std::string cmf::kFCAlphaStrings[2] = {"68% CL","90% CL"}

Definition at line 406 of file Constants.h.

const unsigned char cmf::kFormZone = 39
static
const unsigned char cmf::kFrAbs_N = 43
static
const unsigned char cmf::kFrAbs_pi = 65
static

Definition at line 103 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kFrCEx_N = 41
static
const unsigned char cmf::kFrCEx_pi = 63
static

Definition at line 101 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kFrElas_N = 42
static
const unsigned char cmf::kFrElas_pi = 64
static

Definition at line 102 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kFrInel_N = 67
static

Definition at line 105 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kFrInel_pi = 40
static
const unsigned char cmf::kFrPiProd_N = 68
static

Definition at line 106 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kFrPiProd_pi = 66
static

Definition at line 104 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kHad_RecoE = 3
static
const unsigned char cmf::kHad_RecoE_UBL = 74
static

Definition at line 115 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kLep_RecoE = 4
static
const unsigned char cmf::kLep_RecoE_MCFrac = 6
static
const unsigned char cmf::kLep_RecoE_UBL = 73
static

Definition at line 114 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kMaCCQE = 33
static
const unsigned char cmf::kMaCCRES = 34
static
const unsigned char cmf::kMaCOHpi = 50
static

Definition at line 88 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kMaNCEL = 46
static

Definition at line 84 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kMaNCRES = 35
static
const unsigned char cmf::kMFP_N = 62
static

Definition at line 100 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kMFP_pi = 61
static

Definition at line 99 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kMinus1Sigma = 30
static

Definition at line 65 of file VarVals.h.

Referenced by cmf::MCVarVals::SetGENIEVarVal(), and VarNameToKey().

const unsigned char cmf::kMinus2Sigma = 29
static

Definition at line 64 of file VarVals.h.

Referenced by cmf::MCVarVals::SetGENIEVarVal(), and VarNameToKey().

const unsigned char cmf::kMvCCRES = 36
static
const unsigned char cmf::kMvNCRES = 49
static

Definition at line 87 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kNormDISCC = 56
static

Definition at line 94 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kNu_RecoE = 5
static
const unsigned char cmf::kNu_RecoE_NuE_Peripheral = 77
static

Definition at line 118 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kNu_RecoE_Quantiles_NuMu = 76
static

Definition at line 117 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kNu_RecoE_UBL = 75
static

Definition at line 116 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kNuE_CVN = 1
static
const unsigned char cmf::kNuE_NumMichel = 2
static
const int cmf::kNumFCAlpha = 2

Definition at line 403 of file Constants.h.

const unsigned char cmf::kOtherDIS_Weight = 28
static
const unsigned char cmf::kPlus1Sigma = 31
static

Definition at line 66 of file VarVals.h.

Referenced by cmf::MCVarVals::SetGENIEVarVal(), and VarNameToKey().

const unsigned char cmf::kPlus2Sigma = 32
static

Definition at line 67 of file VarVals.h.

Referenced by cmf::MCVarVals::SetGENIEVarVal(), and VarNameToKey().

const unsigned char cmf::kR0COHpi = 51
static

Definition at line 89 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kRecoQ2 = 7
static
const unsigned char cmf::kRnubarnuCC = 57
static

Definition at line 95 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kRPACCQE_Weight = 22
static
const unsigned char cmf::kRPACCQEshapeEnh = 37
static
const unsigned char cmf::kRPACCQEshapeSupp = 38
static
const unsigned char cmf::kRPARES_Weight = 23
static
const unsigned char cmf::kSumSmallXSecJoint2017 = 45
static

Definition at line 83 of file VarVals.h.

const unsigned char cmf::kTheta_Delta2Npi = 72
static

Definition at line 110 of file VarVals.h.

Referenced by isSumSmallXSecJoint2017Key(), and VarNameToKey().

const unsigned char cmf::kTrueCCNC = 11
static
const unsigned char cmf::kTrueE = 9
static
const unsigned char cmf::kTrueHitNuc = 19
static
const unsigned char cmf::kTrueIntMode = 20
static
const unsigned char cmf::kTrueIntType = 12
static
const unsigned char cmf::kTrueParentDecay = 17
static
const unsigned char cmf::kTrueParentPDG = 14
static
const unsigned char cmf::kTrueParentPT = 15
static
const unsigned char cmf::kTrueParentPZ = 16
static
const unsigned char cmf::kTrueParentTargetPDG = 18
static
const unsigned char cmf::kTruePDG = 10
static
const unsigned char cmf::kTruePDGOrig = 13
static
const unsigned char cmf::kUnknownVar = 0
static

Definition at line 26 of file VarVals.h.

const unsigned char cmf::kVecCCQE_shape = 48
static

Definition at line 86 of file VarVals.h.

Referenced by VarNameToKey().

const unsigned char cmf::kXSecCVPPFX_Weight = 21
static