Xdiff_gwt.C File Reference

Go to the source code of this file.

Functions
	for ()
	for (int ibin=1;ibin< hgenie->GetNbinsX()+1;++ibin) Xtrue[ibin-1]
gStyle	SetOptStat (0)
c1	cd ()
h2f	SetXTitle ("E_{#nu_{e}} (GeV)")
h2f	SetYTitle ("#nu_{e} CC #sigma (10^{-38} cm^{2}/nucleon)")
h2f	Draw ()
bin1_2sigma	SetLineColor (3)
bin1_2sigma	SetFillColor (3)
bin1_2sigma	Draw ("same")
bin1_1sigma	SetLineColor (5)
bin1_1sigma	SetFillColor (5)
cout<<"genie syst: "<< Xtrue_1p[0]<<" "<< Xtrue_1p[1]<<" "<< Xtrue_1p[2]<<" "<< Xtrue_1p[3]<< endl;hgenie->	SetLineWidth (3)
hgenie	SetLineColor (4)
hgenie	SetLineStyle (2)
hgenie	Draw ("HIST same")
GRdata_sys	SetLineWidth (2)
GRdata_sys	SetLineColor (2)
GRdata_sys	Draw ("P z same")
gStyle	SetEndErrorSize (4)
GRdata	SetMarkerStyle (8)
GRdata	SetMarkerSize (1.0)
GRdata	Draw ("P same")
leg	SetLineColor (0)
leg	SetFillColor (0)
leg	AddEntry (GRdata,"data","p")
leg	AddEntry (GRdata,"stat. uncert.","E")
leg	AddEntry (GRdata_sys,"stat. + syst. uncert.","E")
leg	AddEntry (hgenie,"GENIE","l")
leg	AddEntry (bin1_1sigma,"GENIE #pm 1#sigma","f")
leg	AddEntry (bin1_2sigma,"GENIE #pm 2#sigma","f")

Variables
double	T =2.68
double	pot_flux =1.
double	Ebins [11] ={0.}
TFile *	fgwt =new TFile("Etrue_gwt_total_ccnue.root")
TH1F *	gwt_1p =(TH1F*)fgwt->Get("ccnue_1p_total")
TH1F *	gwt_2p =(TH1F*)fgwt->Get("ccnue_2p_total")
TH1F *	gwt_1m =(TH1F*)fgwt->Get("ccnue_1m_total")
TH1F *	gwt_2m =(TH1F*)fgwt->Get("ccnue_2m_total")
vector< double >	Xtrue_1p
vector< double >	Xtrue_2p
vector< double >	Xtrue_1m
vector< double >	Xtrue_2m
TFile *	fEreco =new TFile("Ereco.root")
double	Nsub [10] ={0.}
double	Nsub_err [10] ={0.}
double	Nsub_dataerr [10] ={0.}
double	Nsub_bkgerr [10] ={0.}
double	Nbkg [10] ={0.}
double	Nbkg_syst [10] ={0.}
double	Nsig [10] ={0.}
TH1F *	hd =(TH1F*)fEreco->Get("Ereco_data")
TH1F *	hbkg =(TH1F*)fEreco->Get("Ereco_bkg")
TH1F *	hsig =(TH1F*)fEreco->Get("Ereco_sig")
double	flux_syst [10] ={0.,0.,0.101,0.114,0.107,0.121,0,0,0,0}
double	ehadron_syst [10] ={0.,0.,0.099,0.069,0.054,0.015,0,0,0,0}
double	seleff_syst [10] ={0.,0.,0.071,0.071,0.071,0.071,0,0,0,0}
double	bkg_sf [10] ={0,0,0.95,0.95,0.95,0.95,0,0,0,0}
double	bkg_sf_err [10] ={0,0,0.21,0.21,0.21,0.21,0,0,0,0}
double	bkg_sf_err_rel [10] ={0.}
double	unfold_scale [10] ={0.}
TFile *	funfold =new TFile("BinByBinUnfold.root")
TH1F *	hunfold =(TH1F*)funfold->Get("ratio")
double	ideff [10] ={0.}
TFile *	feff =new TFile("SelEff_nueCC.root")
TH1F *	hid =(TH1F*)feff->Get("toteff")
TFile *	fflux =new TFile("newflux.root")
double	flux [10] ={0.}
double	flux_weight [10] ={0,0,0.903,0.906,0.915,0.932,0,0,0,0}
TH1F *	hflux =(TH1F*)fflux->Get("nue_05_fid")
double	total_syst [10] ={0.}
double	Xdata [10] ={0.}
double	Xbins [10] ={0.}
double	Xbins_err [10] ={0.}
double	Xdata_err [10] ={0.}
double	Xdata_allerr [10] ={0.}
TCanvas *	c1 =new TCanvas("Xdiff_bdt_gwt","",700,600)
TH2F *	h2f =new TH2F("h2f","NOvA, 2.6 #times 10^{20} POT",10,1.001,2.999,10,0,5.49)
TBox *	bin1_2sigma =new TBox(1.,Xtrue_2m[0],1.5,Xtrue_2p[0])
TBox *	bin2_2sigma =new TBox(1.5,Xtrue_2m[1],2.0,Xtrue_2p[1])
TBox *	bin3_2sigma =new TBox(2.0,Xtrue_2m[2],2.5,Xtrue_2p[2])
TBox *	bin4_2sigma =new TBox(2.5,Xtrue_2m[3],3.0,Xtrue_2p[3])
TBox *	bin1_1sigma =new TBox(1.,Xtrue_1m[0],1.5,Xtrue_1p[0])
TBox *	bin2_1sigma =new TBox(1.5,Xtrue_1m[1],2.0,Xtrue_1p[1])
TBox *	bin3_1sigma =new TBox(2.0,Xtrue_1m[2],2.5,Xtrue_1p[2])
TBox *	bin4_1sigma =new TBox(2.5,Xtrue_1m[3],3.0,Xtrue_1p[3])
	GRdata_sys = new TGraphErrors(10,Xbins,Xdata,Xbins_err,Xdata_allerr)
	GRdata = new TGraphErrors(10,Xbins,Xdata,Xbins_err,Xdata_err)
TLegend *	leg = new TLegend(0.6,0.6,0.89,0.89)

Function Documentation

leg AddEntry	(	GRdata	,
		"data"	,
		"p"
	)

Referenced by demoStarPlots(), canMan::drawOverlay(), drawSensitivity(), energyResolution(), make_starplots(), MakeCanvas(), MakeComparisonPlots(), MakeSplitCans(), PlotHistProp(), and PlotMCOnly().

leg AddEntry	(	GRdata	,
		"stat. uncert."	,
		"E"
	)

leg AddEntry	(	GRdata_sys	,
		"stat. + syst. uncert."	,
		"E"
	)

leg AddEntry	(	hgenie	,
		"GENIE"	,
		"l"
	)

leg AddEntry	(	bin1_1sigma	,
		"GENIE #pm 1#sigma"	,
		"f"
	)

leg AddEntry	(	bin1_2sigma	,
		"GENIE #pm 2#sigma"	,
		"f"
	)

c1 cd

(

)

h2f Draw

(

)

fBaryonPT2pdfTwk2 Draw

(

"same"

)

hgenie Draw

(

"HIST same"

)

GRdata_sys Draw

(

"P z same"

)

nova_stat Draw

(

"P same"

)

for

(

)

Definition at line 9 of file Xdiff_gwt.C.

                   {0.};

for

(

int

ibin = 1; ibin<hgenie->GetNbinsX()+1; ++ibin

)

Definition at line 31 of file Xdiff_gwt.C.

References make_syst_table_plots::ibin.

                                                       {
    double binctr=gwt_1p->GetBinCenter(ibin);
    if( binctr<1. || binctr>3. ) continue;
    Xtrue_1p.push_back((1.+gwt_1p->GetBinContent(ibin))*hgenie->GetBinContent(ibin)); 
    Xtrue_2p.push_back((1.+gwt_2p->GetBinContent(ibin))*hgenie->GetBinContent(ibin)); 
    Xtrue_1m.push_back((1.-gwt_1m->GetBinContent(ibin))*hgenie->GetBinContent(ibin)); 
    Xtrue_2m.push_back((1.-gwt_2m->GetBinContent(ibin))*hgenie->GetBinContent(ibin)); 
    //cout<<binctr<<" "<<gwt_1p->GetBinContent(ibin)<<" "<<gwt_2p->GetBinContent(ibin)<<" "<<gwt_1m->GetBinContent(ibin)<<" "<<gwt_2m->GetBinContent(ibin)<<" "<<hgenie->GetBinContent(ibin)<<endl;
    //cout<<binctr<<" "<<Xtrue_1p[ibin-1]<<" "<<Xtrue_1m[ibin-1]<<" "<<Xtrue_2p[ibin-1]<<" "<<Xtrue_2m[ibin-1]<<endl;
  }

gStyle SetEndErrorSize

(

4

)

h1bis SetFillColor

(

3

)

Referenced by genie::masterclass::GNuMcMainFrame::BuildMCTruthTab(), compare_fits(), demoStarPlots(), canMan::drawOverlay(), drawSensitivity(), energyResolution(), extract_hadnucleus_xsec(), hA_basic_analysis(), Init(), make_starplots(), MakeNuEnergyPlots(), MakePlots(), MakeSplitCans(), PlotPreSelectionPlots(), PlotSelectionPlots(), SetHistProp(), and starPlot().

bin4_1sigma SetFillColor

(

5

)

pad3 SetFillColor

(

0

)

h1bis SetLineColor

(

3

)

bin4_1sigma SetLineColor

(

5

)

hgenie SetLineColor

(

4

)

lo_gr_750 SetLineColor

(

2

)

leg SetLineColor

(

0

)

hgenie SetLineStyle

(

2

)

grerrn SetLineWidth

(

3

)

GRdata_sys SetLineWidth

(

2

)

GRdata SetMarkerSize

(

1.

0

)

Referenced by genie::utils::style::Format(), and genie::utils::style::SetDefaultStyle().

GRdata SetMarkerStyle

(

8

)

gStyle SetOptStat

(

0

)

h2f SetXTitle

(

"E_{#nu_{e}} (GeV)"

)

h2f SetYTitle

(

"#nu_{e} CC #sigma (10^{-38} cm^{2}/nucleon)"

)

Variable Documentation

TBox* bin1_1sigma =new TBox(1.,Xtrue_1m[0],1.5,Xtrue_1p[0])

Definition at line 188 of file Xdiff_gwt.C.

TBox* bin1_2sigma =new TBox(1.,Xtrue_2m[0],1.5,Xtrue_2p[0])

Definition at line 171 of file Xdiff_gwt.C.

TBox* bin2_1sigma =new TBox(1.5,Xtrue_1m[1],2.0,Xtrue_1p[1])

Definition at line 189 of file Xdiff_gwt.C.

TBox* bin2_2sigma =new TBox(1.5,Xtrue_2m[1],2.0,Xtrue_2p[1])

Definition at line 172 of file Xdiff_gwt.C.

TBox* bin3_1sigma =new TBox(2.0,Xtrue_1m[2],2.5,Xtrue_1p[2])

Definition at line 190 of file Xdiff_gwt.C.

TBox* bin3_2sigma =new TBox(2.0,Xtrue_2m[2],2.5,Xtrue_2p[2])

Definition at line 173 of file Xdiff_gwt.C.

TBox* bin4_1sigma =new TBox(2.5,Xtrue_1m[3],3.0,Xtrue_1p[3])

Definition at line 191 of file Xdiff_gwt.C.

TBox* bin4_2sigma =new TBox(2.5,Xtrue_2m[3],3.0,Xtrue_2p[3])

Definition at line 174 of file Xdiff_gwt.C.

double bkg_sf[10] ={0,0,0.95,0.95,0.95,0.95,0,0,0,0}

Definition at line 67 of file Xdiff_gwt.C.

double bkg_sf_err[10] ={0,0,0.21,0.21,0.21,0.21,0,0,0,0}

Definition at line 68 of file Xdiff_gwt.C.

double bkg_sf_err_rel[10] ={0.}

Definition at line 69 of file Xdiff_gwt.C.

TCanvas* c1 =new TCanvas("Xdiff_bdt_gwt","",700,600)

Definition at line 162 of file Xdiff_gwt.C.

double Ebins[11] ={0.}

Definition at line 8 of file Xdiff_gwt.C.

Referenced by ana::ExPIDPlot(), and ExPIDPlot().

double ehadron_syst[10] ={0.,0.,0.099,0.069,0.054,0.015,0,0,0,0}

Definition at line 62 of file Xdiff_gwt.C.

TFile* feff =new TFile("SelEff_nueCC.root")

Definition at line 105 of file Xdiff_gwt.C.

TFile* fEreco =new TFile("Ereco.root")

Definition at line 48 of file Xdiff_gwt.C.

TFile* fflux =new TFile("newflux.root")

Definition at line 118 of file Xdiff_gwt.C.

Referenced by Make_NuMuCC_Inc_XS().

TFile* fgwt =new TFile("Etrue_gwt_total_ccnue.root")

Definition at line 20 of file Xdiff_gwt.C.

double flux[10] ={0.}

Definition at line 119 of file Xdiff_gwt.C.

double flux_syst[10] ={0.,0.,0.101,0.114,0.107,0.121,0,0,0,0}

Definition at line 61 of file Xdiff_gwt.C.

double flux_weight[10] ={0,0,0.903,0.906,0.915,0.932,0,0,0,0}

Definition at line 120 of file Xdiff_gwt.C.

TFile* funfold =new TFile("BinByBinUnfold.root")

Definition at line 96 of file Xdiff_gwt.C.

GRdata = new TGraphErrors(10,Xbins,Xdata,Xbins_err,Xdata_err)

Definition at line 220 of file Xdiff_gwt.C.

GRdata_sys = new TGraphErrors(10,Xbins,Xdata,Xbins_err,Xdata_allerr)

Definition at line 213 of file Xdiff_gwt.C.

TH1F* gwt_1m =(TH1F*)fgwt->Get("ccnue_1m_total")

Definition at line 23 of file Xdiff_gwt.C.

TH1F* gwt_1p =(TH1F*)fgwt->Get("ccnue_1p_total")

Definition at line 21 of file Xdiff_gwt.C.

TH1F* gwt_2m =(TH1F*)fgwt->Get("ccnue_2m_total")

Definition at line 24 of file Xdiff_gwt.C.

TH1F* gwt_2p =(TH1F*)fgwt->Get("ccnue_2p_total")

Definition at line 22 of file Xdiff_gwt.C.

TH2F* h2f =new TH2F("h2f","NOvA, 2.6 #times 10^{20} POT",10,1.001,2.999,10,0,5.49)

Definition at line 165 of file Xdiff_gwt.C.

TH1F* hbkg =(TH1F*)fEreco->Get("Ereco_bkg")

Definition at line 58 of file Xdiff_gwt.C.

Referenced by bin_composition_pie_chart(), fd_plot(), selection_story_plots(), and TableNuePredictionFourBins().

TH1F* hd =(TH1F*)fEreco->Get("Ereco_data")

Definition at line 57 of file Xdiff_gwt.C.

Referenced by evd::HardwarePad::Config(), om::HistoSet::CopyAndResetOne(), for(), ana::SpectrumLoaderBase::GetNextFile(), om::CaptionBox::GUIModelDataIssue(), om::ComparisonBox::GUIModelDataIssue(), om::HistoTable::LookUp(), om::HistogramBrowser::MakeCategories(), plot_nd_data_mc(), om::HistogramBrowser::Populate(), evd::RawDataDrawer::RawDigitHardware(), om::PlotViewer::Update(), om::PlotViewer::UpdateCompare(), and om::HistoSet::UTCReset().

TH1F* hflux =(TH1F*)fflux->Get("nue_05_fid")

Definition at line 121 of file Xdiff_gwt.C.

Referenced by readNtuple::Loop(), ana::CrossSectionAnalysis::Result(), ana::NumuCCIncAnalysis::Result(), ana::NumuCC2p2hAnalysis::Result1DFluxInt(), ana::SingleNucAnalysis::Result2D(), ana::NumuCC2p2hAnalysis::Result2D(), and ana::NumuCCIncAnalysis::Result2D().

TH1F* hid =(TH1F*)feff->Get("toteff")

Definition at line 106 of file Xdiff_gwt.C.

TH1F* hsig =(TH1F*)fEreco->Get("Ereco_sig")

Definition at line 59 of file Xdiff_gwt.C.

Referenced by CVNCuts(), ana::DefaultNueLegend(), MakePeriCutPlot(), PeripheralCuts(), PlotNuePredictionFourBins(), TableNuePredictionFourBins(), and test_nue2017Prediction().

TH1F* hunfold =(TH1F*)funfold->Get("ratio")

Definition at line 97 of file Xdiff_gwt.C.

double ideff[10] ={0.}

Definition at line 104 of file Xdiff_gwt.C.

Referenced by cheat::BackTracker::TracksToParticles().

TLegend* leg = new TLegend(0.6,0.6,0.89,0.89)

Definition at line 226 of file Xdiff_gwt.C.

double Nbkg[10] ={0.}

Definition at line 53 of file Xdiff_gwt.C.

Referenced by nuebar_signif().

double Nbkg_syst[10] ={0.}

Definition at line 54 of file Xdiff_gwt.C.

double Nsig[10] ={0.}

Definition at line 55 of file Xdiff_gwt.C.

double Nsub[10] ={0.}

Definition at line 49 of file Xdiff_gwt.C.

double Nsub_bkgerr[10] ={0.}

Definition at line 52 of file Xdiff_gwt.C.

double Nsub_dataerr[10] ={0.}

Definition at line 51 of file Xdiff_gwt.C.

double Nsub_err[10] ={0.}

Definition at line 50 of file Xdiff_gwt.C.

double pot_flux =1.

Definition at line 6 of file Xdiff_gwt.C.

double seleff_syst[10] ={0.,0.,0.071,0.071,0.071,0.071,0,0,0,0}

Definition at line 63 of file Xdiff_gwt.C.

double T =2.68

Definition at line 5 of file Xdiff_gwt.C.

Referenced by NdmcMetric< T >::_calculateSendValue(), dbutils::daqconfig::GeneralDataReader::_getListOfValuesForSingleConfig(), dbutils::daqconfig::GeneralDataReader::_getListOfValuesFromDB(), dbutils::daqconfig::GeneralDataReader::_getSetOfValuesForSingleConfig(), dbutils::daqconfig::GeneralDataReader::_getSetOfValuesFromDB(), NdmcMetric< T >::_getType(), dbutils::daqconfig::ConfigDataTree::_getValueFromContainer(), osc::_OscCalcDMP< T >::_OscCalcDMP(), osc::_PMNS_DMP< T >::_PMNS_DMP(), ana::PredictionAddRock::_PredictComponentSyst(), art::ServiceToken::add(), art::detail::CanBeAggregated< T, Enable >::aggregate(), art::detail::CanBeAggregated< std::array< T, N > >::aggregate(), stan::math::stack_alloc::alloc_array(), art::debugging_allocator< T >::allocate(), ana::AlmostEqual(), slicer::S4DParamCalc::analyze(), slicer::SlicerAna::analyze(), stan::math::promote_scalar_struct< T, S >::apply(), stan::math::autocorrelation(), stan::math::autocovariance(), beamlinereco::LEHitFinder< T >::BackwardFindingOfHitStart(), beamlinereco::CFDHitFinder< T >::BackwardFindingOfHitStart(), stan::math::bessel_first_kind(), stan::math::bessel_second_kind(), calib::BetheBlochTables::BetheBlochTables(), osc::_PMNSOpt< T >::BuildHlv(), stan::math::check_ldlt_factor(), cet::checked_hypot(), stan::math::cholesky_corr_constrain(), stan::math::cholesky_factor_constrain(), art::Ptr< L >::Class_Version(), cet::sqlite::detail::columns(), cet::sqlite::detail::concatenate(), fhicl::detail::concatenate_keys(), cet::sqlite::detail::constraints_info(), art::debugging_allocator< T >::construct(), art::Consumer::consumes(), art::Consumer::consumesMany(), art::Consumer::consumesView(), fhicl::OptionalTupleAs< T(ARGS...)>::conversion_comment(), stan::io::reader< T >::corr(), stan::math::corr_constrain(), HepGeom::BasicVector3D< float >::cosTheta(), chaninfo::BadChanList::CountDCMHits(), art::ServiceToken::createInvalid(), cet::cube(), fhicl::detail::decode(), calib::BetheBlochTables::dEdx(), ana::DemangledTypeName(), stan::mcmc::softabs_metric< Model, BaseRNG >::dG_dt(), cet::diff_of_squares(), cet::LibraryManager::dllExtPattern(), art::Wrapper< T >::do_getElementAddresses(), genie::utils::gsl::d2Xsec_dTCosth::DoEval(), stan::math::domain_error(), stan::math::internal::empty_broadcast_array< ViewElt, std::vector< std::vector< Op > > >::empty_broadcast_array(), calib::EnergyLossVsDistance::EnergyLossVsDistance(), fhicl::detail::ensure_callable(), osc::_PMNS< T >::EvalEqn5(), caf::EvalInstanceWrapper(), expect_same(), fhicl::detail::no_defaults::expected_types< T >::expected_types(), genie::LwlynSmithFF::F1N(), genie::LwlynSmithFF::F1P(), genie::LwlynSmithFF::F2N(), genie::LwlynSmithFF::F2P(), stan::math::F32(), stan::math::falling_factorial(), art::detail::fhicl_key(), art::View< T >::fill(), flat::Flat< int >::Fill(), flat::Flat< std::vector< T > >::Fill(), art::CannotFillView< std::vector< bool, A > >::fill(), fhicl::OptionalTuple< ARGS... >::finalize_tuple_elements(), fhicl::Tuple< TYPES >::finalize_tuple_elements(), beamlinereco::CFDHitFinder< T >::FindCFDHits(), beamlinereco::CFDHitFinder< T >::FindHitWidth(), beamlinereco::LEHitFinder< T >::FindLEHits(), caf::CAFMaker::FindManyPStrict(), beamlinereco::CFDHitFinder< T >::FindPeakValue(), beamlinereco::LEHitFinder< T >::FindPedestal(), beamlinereco::CFDHitFinder< T >::FindPedestal(), beamlinereco::LEHitFinder< T >::FindRawHitLogic(), beamlinereco::CFDHitFinder< T >::FindRawHitLogic(), flat::Flat< int >::Flat(), beamlinereco::LEHitFinder< T >::ForwardFindingOfHitFallTime(), beamlinereco::CFDHitFinder< T >::ForwardFindingOfHitFallTime(), cet::fourth(), stan::math::gamma_p(), stan::math::gamma_q(), stan::math::gaussian_dlm_obs_lpdf(), cet::sqlite::detail::constraints_info_impl<>::generate_snippet(), art::Transient< art::ProductProvenance::Transients >::get(), fhicl::DelegatedParameter::get(), dbutils::SimpleParameterSet::get(), nova::database::Column::Get(), art::ServicesManager::get(), art::detail::ServiceCacheEntry::get(), nutools::dbi::Column::Get(), fhicl::ParameterSet::get(), art::ServiceRegistry::get(), nova::dbi::Column::Get(), fhicl::OptionalDelegatedParameter::get_if_present(), fhicl::ParameterSet::get_if_present(), cet::get_underlying(), ana::GenericSystComponentScale< T >::GetCut(), art::Ptr< T >::getData_(), calib::IBetheBloch::GetdEdxGraph(), beamlinereco::CFDHitFinder< double >::GetDiscriminationThreshold(), art::getElementAddresses(), HepGeom::BasicVector3D< float >::getEta(), sim::FLSHit::getLast(), NdmcClient::getLastValueSent(), beamlinereco::LEHitFinder< double >::GetLEThreshold(), cet::LibraryManager::getLoadableLibraries(), calib::IBetheBloch::GetMPVGraph(), art::root::getObjectRequireDict(), osc::_OscCalcDMP< T >::GetParamsHash(), osc::_IOscCalcAdjustable< T >::GetParamsHashDefault(), evd::SliceNavigator::GetProductsImplFull(), evd::SliceNavigator::GetProductsImplPart(), caf::CAFMaker::GetPsetParameter(), evd::RawDataDrawer::GetRawDigits(), NdmcClient::getSendTimeInterval(), sim::FLSHit::getSum(), cet::LibraryManager::getSymbolByLibspec(), cet::LibraryManager::getSymbolByPath(), dbutils::daqconfig::ConfigDataTree::getTableName(), caf::GetTypedValueWrapper(), NdmcMetric< T >::getValue(), caf::Proxy< T >::GetValueFlatMulti(), caf::Proxy< T >::GetValueFlatSingle(), caf::Proxy< T >::GetValueNested(), genie::GetValueOrUseDefault(), ana::_HistAxis< Var >::GetVars(), stan::math::grad_2F1(), stan::math::grad_F32(), stan::math::grad_reg_inc_beta(), art::Handle< T >::Handle(), stan::math::hypot(), stan::math::inc_beta(), stan::math::inc_beta_dda(), stan::math::inc_beta_ddb(), cet::sqlite::IncompleteInsert::IncompleteInsert(), bpfit::dEdxTable::InitGlue(), bpfit::dEdxTable::InitPVC(), bpfit::dEdxTable::InitScintillator(), bpfit::dEdxTable::InitSteel(), bpfit::dEdxTable::InitVac(), art::prevent_recursion< T >::insert_if_sampled_product(), bpfit::Path::IntegrateT(), beamlinereco::LEHitFinder< T >::IntegrateWaveformInADC(), beamlinereco::CFDHitFinder< T >::IntegrateWaveformInADC(), stan::math::inv_Phi(), stan::math::inv_sqrt(), stan::math::inv_square(), stan::math::inverse_spd(), util::ipow(), fhicl::detail::is_table(), art::ServiceRegistry::isAvailable(), art::ServicesManager::isAvailable(), ana::SystShifts::IsNominal(), art::Wrapper< T >::isPresent_(), art::SummedValue< T >::isValid(), stan::math::lmgamma(), stan::math::log_modified_bessel_first_kind(), ana::StanFitter::log_prob(), stan::math::log_softmax(), stan::math::log_sum_exp(), art::detail::value_type_helper::look_past_pair(), calib::ltTPt(), stan::math::lub_constrain(), art::TFileDirectory::make(), art::WrappedTypeID::make(), hdf5::make_col(), stan::math::make_nu(), art::detail::tool_type< T, std::enable_if_t< std::is_function< T >::value > >::make_plugin(), ana::make_simple_res_plot(), art::TFileDirectory::makeAndRegister(), art::detail::makeServiceFrom(), art::detail::makeWatchFunc(), ana::MCMCSamples::MarginalizeTo(), stan::math::matrix_exp_2x2(), stan::math::max(), art::debugging_allocator< T >::max_size(), cet::sqlite::detail::maybe_quote(), art::Consumer::mayConsume(), art::Consumer::mayConsumeMany(), art::Consumer::mayConsumeView(), stan::math::mean(), art::rootNames::metaBranchRootName(), cet::sqlite::min(), stan::math::min(), calib::BetheBlochAnalytic::MIP(), stan::math::modified_bessel_first_kind(), stan::math::modified_bessel_second_kind(), art::WorkerT< OutputModule >::modifiesEvent(), art::WorkerT< OutputModule >::module(), caf::Proxy< caf::caf::SRTruthBranch >::Name(), slicer::PointManager::NeighborScore(), novaddt::DBSlicer::NeighborScore(), caf::Proxy< unsigned int >::operator unsigned int(), ana::_Cut< T >::operator!(), ana::_Cut< T >::operator&&(), stan::model::model_functional< M >::operator()(), stan::mcmc::softabs_fun< Model >::operator()(), ana::One< T >::operator()(), ana::NoCut< T >::operator()(), cet::detail::pow< 2u, T, 0u >::operator()(), fhicl::detail::it_value_get< T, typename tt::disable_if< std::is_reference< T >::value|| std::is_pointer< T >::value >::type >::operator()(), art::detail::ItemGetter< T, C >::operator()(), fhicl::YesInsert::operator()(), fhicl::detail::it_value_get< intermediate_table::table_t const & >::operator()(), art::DoNotMakePartner< T >::operator()(), art::DoNotSetPtr< T >::operator()(), ana::_Var< T >::operator*(), art::ServiceHandle< art::FileTransfer >::operator*(), art::ServiceHandle< T, art::ServiceScope::PER_SCHEDULE >::operator*(), art::Handle< T >::operator*(), caf::Proxy< T >::operator*=(), ana::_Var< T >::operator+(), caf::Proxy< T >::operator+=(), ana::_Var< T >::operator-(), ana::ThreadLocal< std::map< ana::PredictionInterp::Key_t, ana::PredictionInterp::Val_t > >::operator->(), art::Ptr< T >::operator->(), art::ValidHandle< T >::operator->(), ana::_Var< T >::operator/(), art::Transient< art::ProductProvenance::Transients >::operator=(), ana::_Cut< T >::operator||(), fhicl::OptionalTupleAs< T(ARGS...)>::OptionalTupleAs(), HepGeom::BasicVector3D< float >::orthogonal(), stan::math::owens_t(), ParseOption(), HepGeom::BasicVector3D< float >::perp(), HepGeom::BasicVector3D< float >::perp2(), stan::math::Phi(), plotHist_SpectrumCVNID(), HepGeom::Plane3D< T >::point(), stan::math::pow(), ana::SingleSampleExperiment::PredHistIncCosmics(), genie::FKR::Print(), PrintTable(), stan::io::reader< T >::prob(), stan::math::prob_constrain(), art::Path::process(), art::EndPathExecutor::process(), art::Schedule::process(), art::EventProcessor::process_(), stan::math::prod(), fuzz::FuzzyKVertex::produce(), stan::io::product(), art::Ptr< T >::productGetter(), osc::_PMNS< T >::PropMatter(), osc::_PMNSOpt< T >::PropVacuum(), art::Ptr< T >::Ptr(), art::ServicesManager::put(), stan::math::quad_form_sym(), art::rangeSetFor(), stan::math::read_corr_L(), art::detail::readMetadata(), zcl::FMMTracker::recoT(), cet::maybe_ref< T >::ref(), ana::DepMan< T >::RegisterConstruction(), ana::DepMan< T >::RegisterDestruction(), genie::RegistryItem< T >::RegistryItem(), cet::maybe_ref< T >::reseat(), genie::FKR::Reset(), resolutionplotter(), stan::math::rising_factorial(), stan::lang::function_signatures::rng_return_type(), Eigen::matrix_exp_computeUV< MatrixType >::run(), art::WorkerInPath::runWorker(), stan::model::rvalue(), art::Sampled< T >::Sampled(), stan::io::reader< T >::scalar_lb(), stan::io::reader< T >::scalar_lub(), stan::io::reader< T >::scalar_pos(), stan::io::reader< T >::scalar_ub(), cet::sqlite::select(), genie::MECGenerator::SelectNSVLeptonKinematics(), NdmcMetric< T >::send(), art::detail::ServiceCacheEntry::serviceScope(), art::detail::ServiceWrapper< T, SCOPE >::ServiceWrapper(), art::detail::ServiceWrapper< T, art::ServiceScope::PER_SCHEDULE >::ServiceWrapper(), osc::_PMNS< T >::SetDeltaMsqrs(), dq::FlasherFinder::SetInstigatorT(), genie::RegistryItem< T >::SetLocal(), HepGeom::BasicVector3D< float >::setMag(), osc::_PMNS< T >::SetMix(), osc::_PMNS_DMP< T >::SetParameters(), HepGeom::BasicVector3D< float >::setPerp(), HepGeom::BasicVector3D< float >::setPhi(), art::setPtr(), HepGeom::BasicVector3D< float >::setTheta(), osc::_PMNSOpt< T >::SetVacuumEigensystem(), ana::FitSinSq2Theta13::SetValue(), NdmcClient::setValue(), beamlinereco::LEHitFinder< T >::SetWaveform(), beamlinereco::CFDHitFinder< T >::SetWaveform(), ana::PredictionInterp::ShiftBins(), ana::Registry< T >::ShortNameToPtr(), stan::math::simplex_constrain(), stan::math::simplex_free(), stan::math::size_zero(), stan::math::softmax(), osc::analytic::SolveCubic(), osc::_PMNS< T >::SortEigenvalues(), cet::sqlite::column< std::string, Constraints... >::sqlite_type(), cet::square(), stan::io::reader< T >::std_vector(), stan::math::sum(), cet::sum_of_squares(), art::ValidHandle< T >::swap(), Tangible< T >::Tangible(), stan::math::tanh(), stan::mcmc::mock_hamiltonian< mock_model, rng_t >::tau(), stan::mcmc::divergent_hamiltonian< mock_model, rng_t >::tau(), stan::math::tgamma(), art::detail::throw_if_invalid(), stan::math::to_array_1d(), stan::math::to_array_2d(), stan::math::to_fvar(), to_stringP(), stan::math::trigamma_impl(), fhicl::TupleAs< T(ARGS...)>::TupleAs(), cet::unchecked_hypot(), HepGeom::BasicVector3D< float >::unit(), stan::math::unit_vector_constrain(), art::SummedValue< T >::update_impl(), cet::sqlite::detail::values_str_impl(), ana::VarFromNOvARwgt(), stan::math::variance(), vector_multi_test(), vector_uni_test(), cheat::BackTracker::VecToVecRawPtr(), fhicl::OptionalTuple< ARGS... >::visit_element(), fhicl::Tuple< TYPES >::visit_element(), art::Wrapper< T >::Wrapper(), zheevc3(), zheevh3(), zheevq3(), zhetrd3(), osc::_IOscCalc< T >::~_IOscCalc(), and caf::Proxy< T[N]>::~Proxy().

double total_syst[10] ={0.}

Definition at line 127 of file Xdiff_gwt.C.

double unfold_scale[10] ={0.}

Definition at line 95 of file Xdiff_gwt.C.

double Xbins[10] ={0.}

Definition at line 130 of file Xdiff_gwt.C.

double Xbins_err[10] ={0.}

Definition at line 131 of file Xdiff_gwt.C.

double Xdata[10] ={0.}

Definition at line 129 of file Xdiff_gwt.C.

double Xdata_allerr[10] ={0.}

Definition at line 133 of file Xdiff_gwt.C.

double Xdata_err[10] ={0.}

Definition at line 132 of file Xdiff_gwt.C.

vector<double> Xtrue_1m

Definition at line 28 of file Xdiff_gwt.C.

vector<double> Xtrue_1p

Definition at line 26 of file Xdiff_gwt.C.

vector<double> Xtrue_2m

Definition at line 29 of file Xdiff_gwt.C.

vector<double> Xtrue_2p

Definition at line 27 of file Xdiff_gwt.C.