Functions | |
def | kGammaCut (tables) |
Look for events where all prongs are photon like multiliners need to be done with def x(tables): blah return blah then x = Cut(x) More... | |
def | kContain (tables) |
def | kMass (tables) |
Computes the invariant mass of two prong events. More... | |
def | gaussian (x, x0, a, stdev, o) |
Variables | |
kTwoProng | |
Look for vertices with two prongs Note: We don't need to explicitly check if there is a vertex like in CAFAna. More... | |
kGammaCut = Cut(kGammaCut) | |
kContain = Cut(kContain) | |
kPlaneGap | |
kPlaneContig | |
kTruePi0 | |
kMass = Var(kMass) | |
start = time.time() | |
string | dMC = '/pnfs/nova/persistent/users/karlwarb/HDF5-Training-19-02-26/ND-GIBUU-FHC' |
list | filesMC = [os.path.join(dMC,f) for f in os.listdir(dMC) if 'h5caf.h5' in f] |
tablesMC = loader(filesMC, limit=100) | |
string | dData = '/pnfs/nova/persistent/users/karlwarb/HDF5-Training-19-02-26/ND-Data-FHC' |
list | filesData = [os.path.join(dData,f) for f in os.listdir(dData) if 'h5caf.h5' in f] |
tablesData = loader(filesData, limit=100) | |
cutTot = kTwoProng&kGammaCut&kContain&kPlaneContig&kPlaneGap | |
cutBkg = cutTot&~kTruePi0 | |
data = spectrum(tablesData, cutTot, kMass) | |
bkg = spectrum(tablesMC, cutBkg, kMass) | |
tot = spectrum(tablesMC, cutTot, kMass) | |
POT = data.POT() | |
inttot = tot.integral(POT=POT) | |
intbkg = bkg.integral(POT=POT) | |
tuple | pur = (inttot - intbkg)/inttot |
int | nbins = 8 |
tuple | range = (0,400) |
d | |
bins | |
m | |
_ | |
b | |
int | centers = (bins[:-1] + bins[1:])/2 |
dataparam | |
datacov | |
p0 | |
mcparam | |
mccov | |
dataerr = np.sqrt(np.diag(datacov)) | |
mcerr = np.sqrt(np.diag(mccov)) | |
string | datamu = 'Data $\mu$: ' |
string | datasi = 'Data $\sigma$: ' |
string | mcmu = 'MC $\mu$: ' |
string | mcsi = 'MC $\sigma$: ' |
figsize | |
weights | |
histtype | |
color | |
label | |
derr = poisson_conf_interval(d,'frequentist-confidence') | |
yerr | |
fmt | |
handles | |
labels | |
c | |
loc | |
fontsize | |
horizontalalignment | |
verticalalignment | |
transform | |
def PandAna.Demos.pi0_spectra.gaussian | ( | x, | |
x0, | |||
a, | |||
stdev, | |||
o | |||
) |
Definition at line 145 of file pi0_spectra.py.
def PandAna.Demos.pi0_spectra.kContain | ( | tables | ) |
def PandAna.Demos.pi0_spectra.kGammaCut | ( | tables | ) |
Look for events where all prongs are photon like multiliners need to be done with def x(tables): blah return blah then x = Cut(x)
Definition at line 36 of file pi0_spectra.py.
References PandAna.Demos.pi0_spectra.kGammaCut.
def PandAna.Demos.pi0_spectra.kMass | ( | tables | ) |
Computes the invariant mass of two prong events.
Definition at line 67 of file pi0_spectra.py.
References PandAna.Demos.pi0_spectra.kMass.
|
private |
Definition at line 142 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.b |
Definition at line 143 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.bins |
Definition at line 141 of file pi0_spectra.py.
Definition at line 116 of file pi0_spectra.py.
Referenced by CVNphoton(), CVNphotonHist(), FOMCalc(), getTimePeakPlots(), ana.graphAsymmErrorWithBkgScaled(), minimize(), TMVA::TMVAGlob.NormalizeHists(), nuebar_signif(), plot_quantbound_twosamples(), plot_quantile_boundaries_twosamples_2020(), plot_rationoosc(), proj_1D(), ana.RatioAsymmErrorWithBkg(), ana.RatioAsymmErrorWithBkgScaled(), ReMId(), ReMIdHist(), resolutionplotter(), and testfom().
PandAna.Demos.pi0_spectra.c |
Definition at line 179 of file pi0_spectra.py.
Definition at line 148 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.color |
Definition at line 167 of file pi0_spectra.py.
Definition at line 112 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.cutTot = kTwoProng&kGammaCut&kContain&kPlaneContig&kPlaneGap |
Definition at line 111 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.d |
Definition at line 141 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.data = spectrum(tablesData, cutTot, kMass) |
Definition at line 115 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.datacov |
Definition at line 151 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.dataerr = np.sqrt(np.diag(datacov)) |
Definition at line 154 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.datamu = 'Data $\mu$: ' |
Definition at line 157 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.dataparam |
Definition at line 151 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.datasi = 'Data $\sigma$: ' |
Definition at line 158 of file pi0_spectra.py.
string PandAna.Demos.pi0_spectra.dData = '/pnfs/nova/persistent/users/karlwarb/HDF5-Training-19-02-26/ND-Data-FHC' |
Definition at line 105 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.derr = poisson_conf_interval(d,'frequentist-confidence') |
Definition at line 171 of file pi0_spectra.py.
string PandAna.Demos.pi0_spectra.dMC = '/pnfs/nova/persistent/users/karlwarb/HDF5-Training-19-02-26/ND-GIBUU-FHC' |
Definition at line 101 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.figsize |
Definition at line 164 of file pi0_spectra.py.
list PandAna.Demos.pi0_spectra.filesData = [os.path.join(dData,f) for f in os.listdir(dData) if 'h5caf.h5' in f] |
Definition at line 106 of file pi0_spectra.py.
list PandAna.Demos.pi0_spectra.filesMC = [os.path.join(dMC,f) for f in os.listdir(dMC) if 'h5caf.h5' in f] |
Definition at line 102 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.fmt |
Definition at line 172 of file pi0_spectra.py.
Referenced by murem::MuonRemoveAna.analyze(), convert_mcnp_txt_to_root(), skim::SkimmerAna.FillNumuHistograms(), skim::CosmicBeamComparison.FillNumuHistograms(), bsf::BremShowerFilter.filter(), GetCommandLineArgs(), evd::SliceNavigator.GetProductsImplPart(), skim::EvaluatorNumuCCpi.Initialize(), skim::EvaluatorNumu.Initialize(), MakeFileName(), genie::utils::print.P3AsString(), genie::utils::print.P4AsShortString(), genie::utils::print.P4AsString(), lem::MergeMatches.produce(), trk::KalmanTrackMerge.produce(), lem::Preselection.produce(), lem::MakeMatches.produce(), dt::ViewMerger.produce(), lem::FindLEMMatches.produce(), lem::MakePID.produce(), slid::LIDBuilder.produce(), lem::LEM.produce(), dif::DiFShowerFinder.produce(), slid::SPIDBuilder.produce(), remid::RecoMuon.produce(), lem::LEMClient.produce(), upmuana::UpMuProb.produce(), murem::MuonRemove.produce(), gibuu::GiBUURegen.produce(), skim::NumuCCpiSkimmer.Skim(), skim::NumuSkimmer.Skim(), murem::MuonRemove.SliceSelection(), murem::MRCCAna.SliceSelection(), genie::Spline.Spline(), genie::utils::print.Vec3AsString(), genie::pathsegutils.Vec3AsString(), and genie::utils::print.X4AsString().
PandAna.Demos.pi0_spectra.fontsize |
Definition at line 188 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.handles |
Definition at line 178 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.histtype |
Definition at line 167 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.horizontalalignment |
Definition at line 188 of file pi0_spectra.py.
Definition at line 134 of file pi0_spectra.py.
Definition at line 133 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.kContain = Cut(kContain) |
Definition at line 50 of file pi0_spectra.py.
Referenced by PandAna.Demos.pi0_spectra.kContain().
PandAna.Demos.pi0_spectra.kGammaCut = Cut(kGammaCut) |
Definition at line 39 of file pi0_spectra.py.
Referenced by PandAna.Demos.pi0_spectra.kGammaCut().
PandAna.Demos.pi0_spectra.kMass = Var(kMass) |
Definition at line 94 of file pi0_spectra.py.
Referenced by PandAna.Demos.pi0_spectra.kMass().
PandAna.Demos.pi0_spectra.kPlaneContig |
Definition at line 56 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.kPlaneGap |
Definition at line 52 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.kTruePi0 |
Definition at line 61 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.kTwoProng |
Look for vertices with two prongs Note: We don't need to explicitly check if there is a vertex like in CAFAna.
Note 2: One liners can be done with lambda functions directly in the cut.
Definition at line 28 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.label |
Definition at line 167 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.labels |
Definition at line 178 of file pi0_spectra.py.
Referenced by Beam(), CalculateTotalCovariance(), ana.CompareNDDataMCFromVector(), ana.CompareNDDataOneMC(), ana.CompareNDDataTwoMC(), ana.CompareOneShiftPred(), ana.CompareTwoShiftPred(), ana.ComparisonTableNbins(), ana::PredictionCombinePeriods.ComponentCC(), demo_numuOnly(), draw_fn_coverage(), genie_syst_universe(), ana.GetNDComponents(), ana.GetNDDecompsFromDecomp(), ana::Binning.Labels(), ana::NDOscillatableSpectrum.LoadFrom(), ana::ReweightableSpectrum.LoadFrom(), ana::OscillatableSpectrum.LoadFrom(), ana::Spectrum.LoadFrom(), ana.MakeNueSystematicsFile(), plot_michel_final(), plot_ND_DataMC(), plot_ND_DataMC_energybinning(), plot_nd_spectra_2018(), ana::PredictionCombinePeriods.PredictComponentSyst(), ana.PrintRawEventCounts(), reco_minus_true_panels(), ana::PredictionExtendToPeripheral.ReduceHelper(), ana::PredictionExtendToPeripheral.ReduceHelperNC(), evdb::ScanFrame.ScanFrame(), ana::SpectrumHandler.SetLoader(), ana::Binning.SimpleHelper(), ana.sort_chisq_tot(), systematics_extrap_comp_from_pred_interp(), systematics_summary_from_pred_interp(), TableNueNDComponents(), TableNueNDDataMC(), TableNuePredictionFourBins(), ana::Multiverse.ToSpectrum(), ana::Multiverse.Transform(), ana::NOvARwgtSyst.TruthShift(), and ana.XAxisDeltaCPLabels().
PandAna.Demos.pi0_spectra.loc |
Definition at line 185 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.m |
Definition at line 142 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.mccov |
Definition at line 152 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.mcerr = np.sqrt(np.diag(mccov)) |
Definition at line 155 of file pi0_spectra.py.
Referenced by om::OnMonProd.OnMonProd().
PandAna.Demos.pi0_spectra.mcmu = 'MC $\mu$: ' |
Definition at line 159 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.mcparam |
Definition at line 152 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.mcsi = 'MC $\sigma$: ' |
Definition at line 160 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.nbins = 8 |
Definition at line 139 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.p0 |
Definition at line 151 of file pi0_spectra.py.
Referenced by remid::ReMIdTrain.analyze(), beamlinereco::CFDHitFinder< T >.BackwardFindingOfHitStart(), geo.ClampRayToDetectorVolume(), geo::GeometryBase.ClosestApproach(), convert(), genie::GHepRecord.EventGenerationMode(), fill_simple(), trk::KalmanTrack.FilterOnly(), beamlinereco::CFDHitFinder< T >.FindPeakValue(), bpfit::BreakPoint.FitTracks(), genie::utils::xml.GetXMLPathList(), genie::HAIntranuke.Inelastic(), genie::HAIntranuke2018.Inelastic(), genie::GRV98LO.Initialize(), ana.LowHadMuonEAct(), ana.LowHadMuonEActandCat(), ana.LowHadVisibleHadE(), make_cosmic_eff_table(), zcl::SMMCluster.MissC(), ana.MuonEAct(), ana::numubarccpi0.MuonEAct(), ana.MuonEActandCat(), ana::numubarccpi0.MuonEActandCat(), ana.MuonEActive(), ana.MuonEActiveandCatcher(), ana.MuonECat(), nuEEnergy(), cmf.NuEEnergy(), ana.NueRecoE_2017FDFit(), ana.NueRecoE_2018FHCFit(), ana.NueRecoE_2018RHCFit(), ana.NueRecoE_2020FHCFit(), ana.NueRecoE_2020FHCFit_2D3D(), ana.NueRecoE_2020RHCFit(), ana.NueRecoE_2020RHCFit_2D3D(), plot_cosmic_eff(), ana.RecoEFitRHC_poly2(), ana.RecoEFitRHC_poly2X(), ana.RecoEFitRHC_poly3(), ana.RecoEFitRHC_poly3X(), skim::ParametersNue.SANueRecoE(), bpfit::Path.ScrubOne(), genie::flux::GDk2NuFlux.SetFluxWindow(), genie::flux::GNuMIFlux.SetFluxWindow(), ana.TAKEPol4_bias(), ana.TAKEPol6_biasHM(), geo::GeometryTest.testFindPlanes(), numue::NumuEAna.TotalLengthInDetector(), cheat::RecoCheckAna.TotalLengthInDetector(), trk::KalmanTrackAna.TotalLengthInDetector(), dt::Chain.ToTrack(), ana.VisibleHadE(), ana::numubarccpi0.VisibleHadE(), and genie::flux::GDk2NuFluxXMLHelper.~GDk2NuFluxXMLHelper().
PandAna.Demos.pi0_spectra.POT = data.POT() |
Definition at line 123 of file pi0_spectra.py.
Definition at line 135 of file pi0_spectra.py.
Referenced by fuzz::FuzzyKValidate.analyze(), caf::Proxy< caf::SRCosmic >.CheckEquals(), caf::Proxy< caf::SRNeutrino >.CheckEquals(), caf::Proxy< caf::SRParticleTruth >.CheckEquals(), caf::Proxy< caf::SRMRCCParent >.CheckEquals(), flat::Flat< caf::SRCosmic >.Clear(), flat::Flat< caf::SRNeutrino >.Clear(), flat::Flat< caf::SRParticleTruth >.Clear(), flat::Flat< caf::SRMRCCParent >.Clear(), flat::Flat< caf::SRCosmic >.Fill(), flat::Flat< caf::SRNeutrino >.Fill(), flat::Flat< caf::SRParticleTruth >.Fill(), flat::Flat< caf::SRMRCCParent >.Fill(), FOMPlot(), murem::FindParent.MatchToOrigSlice(), caf::Proxy< caf::SRCosmic >.operator=(), caf::Proxy< caf::SRNeutrino >.operator=(), caf::Proxy< caf::SRParticleTruth >.operator=(), caf::Proxy< caf::SRMRCCParent >.operator=(), plot_uncertainty(), and ana.ProngTruthCut().
PandAna.Demos.pi0_spectra.range = (0,400) |
Definition at line 140 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.start = time.time() |
Definition at line 99 of file pi0_spectra.py.
Definition at line 107 of file pi0_spectra.py.
Definition at line 103 of file pi0_spectra.py.
Definition at line 117 of file pi0_spectra.py.
Referenced by lem::MakePID.AvgAndFit(), bin_composition_pie_chart(), ana::UnfoldTikhonov.ChiSq(), DefaultNueLegend(), ana.DefaultNueLegend(), DefaultNumuLegend(), draw_michel_plots(), chaninfo::CosmicEff.endJob(), genie::BaryonResonanceDecayer.EvolveDeltaBR(), CherenkovCalc.FastMCFill(), calib::AttenFit.FitQuality(), calib::AttenuationFit.FitQuality(), initialise(), ana::SpectrumComponents.LoadFrom(), make_cosmic_eff_table(), make_plots(), calib::AttenFit.MedianProfile(), calib::AttenuationFit.MedianProfile(), lem.MMapFileAtAddress(), genie::alvarezruso::AlvarezRusoCOHPiPDXSec.NuclearCurrent(), ana::UnfoldMaxEnt.PenaltyTerm(), PeripheralCuts(), HepGeom::BasicVector3D< float >.perp2(), plot_nd_spectra_2018(), plot_timing_resolution(), readTree(), ana::UnfoldTikhonov.RecoTruthPrediction(), genie::flux::GAstroFlux.SetRelNuPopulations(), test_nue2017Prediction(), testfom(), calib::AttenFit.TruncatedMeanProfile(), and calib::AttenuationFit.TruncatedMeanProfile().
PandAna.Demos.pi0_spectra.transform |
Definition at line 189 of file pi0_spectra.py.
Referenced by novaddt::HotMapMaker.AppendCold(), novaddt::HotMapMaker.AppendHot(), novaddt::HotMapMaker.AppendRate(), novaddt::Multiplet.ApplyCorrection(), ana::BayesianMarginal.BayesianMarginal(), ana.CheckSystsMatchWithNOvARwgt(), CreateFidSelection(), CreateRockBoxSelection(), HepGeom::Point3D< float >.distance(), HepGeom::Point3D< double >.distance(), drawSystsShiftingNDdata(), mcmc.ExptPtrs(), art.for_each_group(), generate_fd_fake_events(), ana::covmx::CovarianceMatrix.GetBinnings(), ana::covmx::CovarianceMatrix.GetNBins(), ana::covmx::CovarianceMatrix.GetNComponents(), getPredictions(), IOManager.IOManager(), novaddt::HotMapMaker.MakeMask(), HepGeom::Vector3D< float >.operator=(), HepGeom::Normal3D< float >.operator=(), HepGeom::Vector3D< double >.operator=(), HepGeom::Normal3D< double >.operator=(), FindNearlineConfig.ParseArguments(), plot_3flavor_withsysts(), ana::MCMCSamples.SaveTo(), HepGeom::Transform3D.setIdentity(), ana::Multiverse.Transform(), and cet.transform_all().
PandAna.Demos.pi0_spectra.verticalalignment |
Definition at line 188 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.weights |
Definition at line 167 of file pi0_spectra.py.
PandAna.Demos.pi0_spectra.yerr |
Definition at line 172 of file pi0_spectra.py.