plots_total.C File Reference

Go to the source code of this file.

Functions
	for (int i=1;i< nexp+1;++i)
	for (int ibin=1;ibin< hnorm->GetNbinsX()+1;++ibin)
gStyle	SetOptStat (0)
c1	cd ()
h2f	SetXTitle ("E_{#nu_{e}} (GeV)")
h2f	SetYTitle ("Flux weight")
h2f	Draw ()
hmean	SetLineColor (4)
hmean	SetLineWidth (2)
hmean	SetLineStyle (2)
hmean	Draw ("HIST same")
h2f2	SetYTitle ("Uncertainty on flux weight")

Variables
TFile *	file =new TFile("gsimple_nue_newflux.root")
TH1F *	hnorm =(TH1F*)file->Get("Enorm")
TH1F *	hweighted [nexp]
float	finalweight [40][nexp]
TH1F *	finalweight_hists [40]
float	weightMean [40] ={0.}
float	weightErr [40] ={0.}
TH1F *	hmean =new TH1F("hmean","",40,0,10)
TBox *	errband [40]
TCanvas *	c1 =new TCanvas("fluxweight_total","",700,600)
TH2F *	h2f =new TH2F("h2f","Flux weight vs. E_{#nu_{e}}",40,0.001,10,10,0.81,1.29)
TCanvas *	c2 =new TCanvas("fluxweight_total_uncertainty","",700,600)
TH2F *	h2f2 =new TH2F("h2f2","Flux weight uncertainty vs. E_{#nu_{e}}",40,0.001,10,10,-0.5,0.5)

Function Documentation

pad3 cd

(

)

h2f Draw

(

)

hmean Draw

(

"HIST same"

)

for

(

)

Definition at line 7 of file plots_total.C.

References MECModelEnuComparisons::i.

                               {
    TString tot_name=Form("Eweight_%i",i);
    hweighted[i-1]=(TH1F*)file->Get(tot_name);
  } 

for

(

int

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

)

Definition at line 14 of file plots_total.C.

References om::cout, den, allTimeWatchdog::endl, make_syst_table_plots::ibin, calib::j, num, and ana::weight.

                                                      {
    for( int j=0; j<nexp; ++j ){
      float den=hnorm->GetBinContent(ibin);
      float num=hweighted[j]->GetBinContent(ibin);
      float weight=0.;
      cout<<j<<" "<<num<<" "<<den<<endl;
      if( den>0.0 ) weight=num/den;
      finalweight[ibin-1][j]=weight;
    }//loop 100 weighted hists
  }//loop energy bins

h88 SetLineColor

(

4

)

Referenced by drawSensitivity(), energyResolution(), genie::utils::style::Format(), genie_syst_pca(), GetSpectToHist(), main(), MakeCanvas(), MakeComparisonPlots(), makeHadEFrac2DPlot(), MakeNuEnergyPlots(), MakeOverlayCanvas(), MakeSplitCans(), NuMu2020_MakePIDPlots(), plot_ehade_quantbound(), plot_quantbound_twosamples(), plot_quantile_boundaries_2020(), plot_quantile_boundaries_twosamples_2020(), plotContProf(), PlotPreSelectionPlots(), PlotSelectionPlots(), ProducePlots(), and Save1DContour().

lo_gr_750 SetLineStyle

(

2

)

Referenced by drawSensitivity(), genie::utils::style::Format(), MakeSplitCans(), plot_3NDvsFD(), plot_3NDvsFD_FHC(), plot_3NDvsFD_RHC(), plot_datamc_ND_numu(), plot_datamc_ND_numu_REW(), plot_ND_numu_NOMvsREW(), plot_NDvsFD_REW(), plot_NDvsFD_weights(), plot_NDvsFD_weights_FHC(), plot_NDvsFD_weights_RHC(), plot_nueFD_Signal_REWvsNOM(), plot_nueFD_Signal_REWvsNOM_FHC(), plot_nueFD_Signal_REWvsNOM_pTExtrap_FHC(), plot_nueFD_Signal_REWvsNOM_pTExtrap_RHC(), plot_nueFD_Signal_REWvsNOM_RHC(), plot_quantbound_twosamples(), plot_quantile_boundaries_twosamples_2020(), plotContProf(), and Save1DContour().

lo_gr_750 SetLineWidth

(

2

)

Referenced by genie::utils::style::Format(), MakeNuEnergyPlots(), MakeSplitCans(), PlotGenRatios(), PlotPreSelectionPlots(), PlotSelectionPlots(), and Save1DContour().

gStyle SetOptStat

(

0

)

h2f2 SetXTitle

(

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

)

h2f SetYTitle

(

"Flux weight"

)

h2f2 SetYTitle

(

"Uncertainty on flux weight"

)

Variable Documentation

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

Definition at line 60 of file plots_total.C.

TCanvas* c2 =new TCanvas("fluxweight_total_uncertainty","",700,600)

Definition at line 72 of file plots_total.C.

TBox* errband[40]

Definition at line 50 of file plots_total.C.

TFile* file =new TFile("gsimple_nue_newflux.root")

Definition at line 3 of file plots_total.C.

float finalweight[40][nexp]

Definition at line 13 of file plots_total.C.

TH1F* finalweight_hists[40]

Definition at line 26 of file plots_total.C.

TH2F* h2f =new TH2F("h2f","Flux weight vs. E_{#nu_{e}}",40,0.001,10,10,0.81,1.29)

Definition at line 62 of file plots_total.C.

Referenced by for().

TH2F* h2f2 =new TH2F("h2f2","Flux weight uncertainty vs. E_{#nu_{e}}",40,0.001,10,10,-0.5,0.5)

Definition at line 74 of file plots_total.C.

Referenced by for().

TH1F* hmean =new TH1F("hmean","",40,0,10)

Definition at line 40 of file plots_total.C.

Referenced by ana::MultiverseCorrelation::calculate_covariance_matrix().

TH1F* hnorm =(TH1F*)file->Get("Enorm")

Definition at line 4 of file plots_total.C.

Referenced by genie::flux::GAtmoFlux::LoadFluxData().

TH1F* hweighted[nexp]

Definition at line 6 of file plots_total.C.

float weightErr[40] ={0.}

Definition at line 39 of file plots_total.C.

float weightMean[40] ={0.}

Definition at line 38 of file plots_total.C.