PandAna.Demos.pi0_spectra Namespace Reference

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

Function Documentation

def PandAna.Demos.pi0_spectra.gaussian	(		x,
			x0,
			a,
			stdev,
			o
	)

Definition at line 145 of file pi0_spectra.py.

    def gaussian(x, x0, a, stdev, o):
        return a * np.exp( - ((x - x0) / stdev) ** 2 / 2) + o

def PandAna.Demos.pi0_spectra.kContain

(

tables

)

Definition at line 42 of file pi0_spectra.py.

References PandAna.Demos.pi0_spectra.kContain.

def kContain(tables):
    df = tables['rec.vtx.elastic']
    return (df['vtx.x'] < 180) & \
        (df['vtx.x'] > -180) & \
        (df['vtx.y'] < 180) & \
        (df['vtx.y'] > -180) & \
        (df['vtx.z'] < 1000) & \
        (df['vtx.z'] > 50)

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 kGammaCut(tables):
    df = tables['rec.vtx.elastic.fuzzyk.png.cvnpart']['photonid']
    return (df > 0.75).groupby(level=KL).agg(np.all)

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.

def kMass(tables):
    # Note: We could leave this check out, but you would get a warning about taking
    # the sqrt of negative numbers at the end (it won't crash like cafana).
    # dataframes can support NaNs just fine.
    check = tables['rec.vtx.elastic.fuzzyk']['npng'] == 2

    df = tables['rec.vtx.elastic.fuzzyk.png'][check]
    x = df['dir.x']
    y = df['dir.y']
    z = df['dir.z']
    
    # Compute the length of the dir vector and then normalize
    l = np.sqrt(x*x+y*y+z*z)
    x = x/l
    y = y/l
    z = z/l

    # compute the dot product
    dot = x.groupby(level=KL).prod()+y.groupby(level=KL).prod()+z.groupby(level=KL).prod()

    # multiply the energy of all prongs in each event together
    EProd = df['calE'].groupby(level=KL).prod()

    # return a dataframe with a single column of the invariant mass
    deadscale = 0.8747
    return 1000*deadscale*np.sqrt(2*EProd*(1-dot))
    # note NaNs can be removed by (df == df)

Variable Documentation

PandAna.Demos.pi0_spectra._

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.

PandAna.Demos.pi0_spectra.bkg = spectrum(tablesMC, cutBkg, kMass)

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.

PandAna.Demos.pi0_spectra.centers = (bins[:-1] + bins[1:])/2

Definition at line 148 of file pi0_spectra.py.

PandAna.Demos.pi0_spectra.color

Definition at line 167 of file pi0_spectra.py.

PandAna.Demos.pi0_spectra.cutBkg = cutTot&~kTruePi0

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.

PandAna.Demos.pi0_spectra.intbkg = bkg.integral(POT=POT)

Definition at line 134 of file pi0_spectra.py.

PandAna.Demos.pi0_spectra.inttot = tot.integral(POT=POT)

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

Initial value:

= Cut(lambda tables:
                        (tables['rec.vtx.elastic.fuzzyk.png']['maxplanecont'] > 4).\
                        groupby(level=KL).agg(np.all))

Definition at line 56 of file pi0_spectra.py.

PandAna.Demos.pi0_spectra.kPlaneGap

Initial value:

= Cut(lambda tables:
                    (tables['rec.vtx.elastic.fuzzyk.png']['maxplanegap'] > 1).\
                    groupby(level=KL).agg(np.all))

Definition at line 52 of file pi0_spectra.py.

PandAna.Demos.pi0_spectra.kTruePi0

Initial value:

= Cut(lambda tables:
                    tables['rec.sand.nue']['npi0'] > 0)

Definition at line 61 of file pi0_spectra.py.

PandAna.Demos.pi0_spectra.kTwoProng

Initial value:

= Cut(lambda tables:
                (tables['rec.vtx.elastic.fuzzyk']['npng'] == 2).groupby(level=KL).agg(np.any))

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.

tuple PandAna.Demos.pi0_spectra.pur = (inttot - intbkg)/inttot

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.

PandAna.Demos.pi0_spectra.tablesData = loader(filesData, limit=100)

Definition at line 107 of file pi0_spectra.py.

PandAna.Demos.pi0_spectra.tablesMC = loader(filesMC, limit=100)

Definition at line 103 of file pi0_spectra.py.

PandAna.Demos.pi0_spectra.tot = spectrum(tablesMC, cutTot, kMass)

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.