Syst.C
Go to the documentation of this file.
1 #ifdef __CINT__
2 
3 //void checkBdt(const std::string fGENIE, const std::string fout)
4 void Syst()
5 {
6  std::cout << "Sorry, you must run in compiled mode" << std::endl;
7 }
8 #else
9 
10 #include "CAFAna/Vars/Vars.h"
11 
12 #include "CAFAna/Core/HistAxis.h"
14 #include "CAFAna/Core/Binning.h"
15 #include "CAFAna/Core/Spectrum.h"
18 #include "CAFAna/Core/Cut.h"
19 #include "CAFAna/Cuts/Cuts.h"
20 #include "CAFAna/Cuts/SpillCuts.h"
22 #include "CAFAna/Cuts/TruthCuts.h"
23 #include "CAFAna/Core/Utilities.h"
24 #include "CAFAna/Vars/Vars.h"
25 #include "CAFAna/Systs/Systs.h"
26 
27 ///////NCPi0 Headers//////
28 #include "/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/ncpi0Bins.h"
29 #include "/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/ncpi0Cuts.h"
30 #include "/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/ncpi0Vars.h"
31 /////////////////////////
32 
33 // Names for the GENIE systs
34 #include "NuReweight/ReweightLabels.h"
35 #include "ReWeight/GSyst.h"
36 
37 #include "TCanvas.h"
38 #include "TH2.h"
39 #include "TStyle.h"
40 #include "TArrow.h"
41 #include "TH1.h"
42 #include "TH3.h"
43 #include "TFile.h"
44 #include "TLegend.h"
45 #include "TGraph.h"
46 //#include "TVector3.h"
47 #include <cmath>
48 #include <cassert>
49 #include <memory>
50 #include "TProfile.h"
51 #include <fstream>
52 #include <sstream>
53 #include <iostream>
54 #include <iomanip>
55 #include <stdarg.h>
56 #include <stdlib.h>
57 #include <stdio.h>
58 #include <string.h>
59 
60 using namespace ana;
61 const double pot = 8.08988e+20;
62 
63 void saveS(Spectrum s, string name){
64  TH1D *h = s.ToTH1(s.POT());
65  h->SetName(name.c_str());
66  h->SetLineColor(kRed);
67  h->Write();
68 }
69 
70 void saveS1(Spectrum s, string name){
71  TH2 *h1 = s.ToTH2(s.POT());
72  // TH2 *h1 = s.ToTH2(pot);
73  h1->SetName(name.c_str());
74  h1->SetLineColor(kRed);
75  h1->Write();
76 }
77 
78 //void checkBdt(const std::string fGENIE, const std::string fout)
79 void Syst()
80 {
81 
82  //For Prod3 FHC data:
83 
84  const std::string fGENIE = "dataset_def_name_newest_snapshot prod_caf_R17-03-01-prod3reco.d_nd_genie_nonswap_fhc_nova_v08_full_v1";
85 
86  SpectrumLoader loader(fGENIE);
88 
89 
90  // nominal
91  Spectrum selPNGCVN("selPNGCVN", PCIDbins, loader, kNpCVN1gamma, kIsPreSel );
92  Spectrum selSigPNGCVN("selSigPNGCVN", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel );
93  Spectrum selBkgPNGCVN("selBkgPNGCVN", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel );
94 
95  //plus 1 sigma
96 
97  Spectrum selPNGCVN_ShapeNorm_AxialMassNCEL_plus1sigma("selPNGCVN_ShapeNorm_AxialMassNCEL_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaNCEL), +1) );
98  Spectrum selSigPNGCVN_ShapeNorm_AxialMassNCEL_plus1sigma("selSigPNGCVN_ShapeNorm_AxialMassNCEL_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaNCEL), +1) );
99  Spectrum selBkgPNGCVN_ShapeNorm_AxialMassNCEL_plus1sigma("selBkgPNGCVN_ShapeNorm_AxialMassNCEL_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaNCEL), +1) );
100 
101  Spectrum selPNGCVN_ShapeNorm_EtaNCEL_plus1sigma("selPNGCVN_ShapeNorm_EtaNCEL_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightEtaNCEL), +1) );
102  Spectrum selSigPNGCVN_ShapeNorm_EtaNCEL_plus1sigma("selSigPNGCVN_ShapeNorm_EtaNCEL_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightEtaNCEL), +1) );
103  Spectrum selBkgPNGCVN_ShapeNorm_EtaNCEL_plus1sigma("selBkgPNGCVN_ShapeNorm_EtaNCEL_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightEtaNCEL), +1) );
104 
105  Spectrum selPNGCVN_ShapeNorm_AxialMassNCRES_plus1sigma("selPNGCVN_ShapeNorm_AxialMassNCRES_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaNCRES), +1));
106  Spectrum selSigPNGCVN_ShapeNorm_AxialMassNCRES_plus1sigma("selSigPNGCVN_ShapeNorm_AxialMassNCRES_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaNCRES), +1) );
107  Spectrum selBkgPNGCVN_ShapeNorm_AxialMassNCRES_plus1sigma("selBkgPNGCVN_ShapeNorm_AxialMassNCRES_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaNCRES), +1));
108 
109 
110  Spectrum selPNGCVN_ShapeNorm_VectorMassNCRES_plus1sigma("selPNGCVN_ShapeNorm_VectorMassNCRES_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMvNCRES), +1) );
111  Spectrum selSigPNGCVN_ShapeNorm_VectorMassNCRES_plus1sigma("selSigPNGCVN_ShapeNorm_VectorMassNCRES_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightMvNCRES), +1) );
112  Spectrum selBkgPNGCVN_ShapeNorm_VectorMassNCRES_plus1sigma("selBkgPNGCVN_ShapeNorm_VectorMassNCRES_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMvNCRES), +1));
113 
114 
115  Spectrum selPNGCVN_ShapeNorm_AxialMassCCRES_plus1sigma("selPNGCVN_ShapeNorm_AxialMassCCRES_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaCCRES), +1));
116  Spectrum selSigPNGCVN_ShapeNorm_AxialMassCCRES_plus1sigma("selSigPNGCVN_ShapeNorm_AxialMassCCRES_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaCCRES), +1) );
117  Spectrum selBkgPNGCVN_ShapeNorm_AxialMassCCRES_plus1sigma("selBkgPNGCVN_ShapeNorm_AxialMassCCRES_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaCCRES), +1) );
118 
119  Spectrum selPNGCVN_ShapeNorm_VectorMassCCRES_plus1sigma("selPNGCVN_ShapeNorm_VectorMassCCRES_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMvCCRES), +1) );
120  Spectrum selSigPNGCVN_ShapeNorm_VectorMassCCRES_plus1sigma("selSigPNGCVN_ShapeNorm_VectorMassCCRES_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMvCCRES), +1) );
121  Spectrum selBkgPNGCVN_ShapeNorm_VectorMassCCRES_plus1sigma("selBkgPNGCVN_ShapeNorm_VectorMassCCRES_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMvCCRES), +1) );
122 
123 
124  Spectrum selPNGCVN_All_AxialMassCOH_plus1sigma("selPNGCVN_All_AxialMassCOH_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaCOHpi), +1));
125  Spectrum selSigPNGCVN_All_AxialMassCOH_plus1sigma("selSigPNGCVN_All_AxialMassCOH_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaCOHpi), +1));
126  Spectrum selBkgPNGCVN_All_AxialMassCOH_plus1sigma("selBkgPNGCVN_All_AxialMassCOH_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaCOHpi), +1));
127 
128 
129  Spectrum selPNGCVN_All_NuclearSizeCOH_plus1sigma("selPNGCVN_All_NuclearSizeCOH_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightR0COHpi), +1) );
130  Spectrum selSigPNGCVN_All_NuclearSizeCOH_plus1sigma("selSigPNGCVN_All_NuclearSizeCOH_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightR0COHpi), +1));
131  Spectrum selBkgPNGCVN_All_NuclearSizeCOH_plus1sigma("selBkgPNGCVN_All_NuclearSizeCOH_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightR0COHpi), +1));
132 
133 
134  Spectrum selPNGCVN_ShapeNorm_DIS_Aht_plus1sigma("selPNGCVN_ShapeNorm_DIS_Aht_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightAhtBY), +1));
135  Spectrum selSigPNGCVN_ShapeNorm_DIS_Aht_plus1sigma("selSigPNGCVN_ShapeNorm_DIS_Aht_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightAhtBY), +1));
136  Spectrum selBkgPNGCVN_ShapeNorm_DIS_Aht_plus1sigma("selBkgPNGCVN_ShapeNorm_DIS_Aht_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightAhtBY), +1));
137 
138 
139  Spectrum selPNGCVN_ShapeNorm_DIS_Bht_plus1sigma("selPNGCVN_ShapeNorm_DIS_Bht_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightBhtBY), +1) );
140  Spectrum selSigPNGCVN_ShapeNorm_DIS_Bht_plus1sigma("selSigPNGCVN_ShapeNorm_DIS_Bht_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightBhtBY), +1) );
141  Spectrum selBkgPNGCVN_ShapeNorm_DIS_Bht_plus1sigma("selBkgPNGCVN_ShapeNorm_DIS_Bht_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightBhtBY), +1));
142 
143 
144  Spectrum selPNGCVN_ShapeNorm_DIS_CV1u_plus1sigma("selPNGCVN_ShapeNorm_DIS_CV1u_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightCV1uBY), +1) );
145  Spectrum selSigPNGCVN_ShapeNorm_DIS_CV1u_plus1sigma("selSigPNGCVN_ShapeNorm_DIS_CV1u_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightCV1uBY), +1) );
146  Spectrum selBkgPNGCVN_ShapeNorm_DIS_CV1u_plus1sigma("selBkgPNGCVN_ShapeNorm_DIS_CV1u_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightCV1uBY), +1));
147 
148 
149  Spectrum selPNGCVN_ShapeNorm_DIS_CV2u_plus1sigma("selPNGCVN_ShapeNorm_DIS_CV2u_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightCV2uBY), +1) );
150  Spectrum selSigPNGCVN_ShapeNorm_DIS_CV2u_plus1sigma("selSigPNGCVN_ShapeNorm_DIS_CV2u_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightCV2uBY), +1));
151  Spectrum selBkgPNGCVN_ShapeNorm_DIS_CV2u_plus1sigma("selBkgPNGCVN_ShapeNorm_DIS_CV2u_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightCV2uBY), +1));
152 
153 
154 
155  Spectrum selPNGCVN_mfp_pi_plus1sigma("selPNGCVN_mfp_pi_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMFP_pi), +1));
156  Spectrum selSigPNGCVN_mfp_pi_plus1sigma("selSigPNGCVN_mfp_pi_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMFP_pi), +1) );
157  Spectrum selBkgPNGCVN_mfp_pi_plus1sigma("selBkgPNGCVN_mfp_pi_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMFP_pi), +1));
158 
159  Spectrum selPNGCVN_mfp_N_plus1sigma("selPNGCVN_mfp_N_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMFP_N), +1));
160  Spectrum selSigPNGCVN_mfp_N_plus1sigma("selSigPNGCVN_mfp_N_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMFP_N), +1) );
161  Spectrum selBkgPNGCVN_mfp_N_plus1sigma("selBkgPNGCVN_mfp_N_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMFP_N), +1));
162 
163  Spectrum selPNGCVN_All_PiChargeExchangeProb_plus1sigma("selPNGCVN_All_PiChargeExchangeProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrCEx_pi), +1) );
164  Spectrum selSigPNGCVN_All_PiChargeExchangeProb_plus1sigma("selSigPNGCVN_All_PiChargeExchangeProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrCEx_pi), +1) );
165  Spectrum selBkgPNGCVN_All_PiChargeExchangeProb_plus1sigma("selBkgPNGCVN_All_PiChargeExchangeProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrCEx_pi), +1));
166 
167  Spectrum selPNGCVN_All_PiElasticProb_plus1sigma("selPNGCVN_All_PiElasticProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrElas_pi), +1) );
168  Spectrum selSigPNGCVN_All_PiElasticProb_plus1sigma("selSigPNGCVN_All_PiElasticProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrElas_pi), +1) );
169  Spectrum selBkgPNGCVN_All_PiElasticProb_plus1sigma("selBkgPNGCVN_All_PiElasticProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrElas_pi), +1) );
170 
171  Spectrum selPNGCVN_All_PiInElasticProb_plus1sigma("selPNGCVN_All_PiInElasticProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrInel_pi), +1) );
172  Spectrum selSigPNGCVN_All_PiInElasticProb_plus1sigma("selSigPNGCVN_All_PiInElasticProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrInel_pi), +1));
173  Spectrum selBkgPNGCVN_All_PiInElasticProb_plus1sigma("selBkgPNGCVN_All_PiInElasticProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrInel_pi), +1));
174 
175 
176 
177  Spectrum selPNGCVN_All_PiAbsorptionProb_plus1sigma("selPNGCVN_All_PiAbsorptionProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrAbs_pi), +1));
178  Spectrum selSigPNGCVN_All_PiAbsorptionProb_plus1sigma("selSigPNGCVN_All_PiAbsorptionProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrAbs_pi), +1) );
179  Spectrum selBkgPNGCVN_All_PiAbsorptionProb_plus1sigma("selBkgPNGCVN_All_PiAbsorptionProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrAbs_pi), +1));
180 
181  Spectrum selPNGCVN_All_PiProductionProb_plus1sigma("selPNGCVN_All_PiProductionProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrPiProd_pi), +1));
182  Spectrum selSigPNGCVN_All_PiProductionProb_plus1sigma("selSigPNGCVN_All_PiProductionProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrPiProd_pi), +1));
183  Spectrum selBkgPNGCVN_All_PiProductionProb_plus1sigma("selBkgPNGCVN_All_PiProductionProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrPiProd_pi), +1));
184 
185 
186  Spectrum selPNGCVN_All_NucChargeExchangeProb_plus1sigma("selPNGCVN_All_NucChargeExchangeProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrCEx_N), +1) );
187  Spectrum selSigPNGCVN_All_NucChargeExchangeProb_plus1sigma("selSigPNGCVN_All_NucChargeExchangeProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrCEx_N), +1) );
188  Spectrum selBkgPNGCVN_All_NucChargeExchangeProb_plus1sigma("selBkgPNGCVN_All_NucChargeExchangeProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrCEx_N), +1));
189 
190 
191  Spectrum selPNGCVN_All_NucElasticProb_plus1sigma("selPNGCVN_All_NucElasticProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrElas_N), +1) );
192  Spectrum selSigPNGCVN_All_NucElasticProb_plus1sigma("selSigPNGCVN_All_NucElasticProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrElas_N), +1) );
193  Spectrum selBkgPNGCVN_All_NucElasticProb_plus1sigma("selBkgPNGCVN_All_NucElasticProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrElas_N), +1) );
194 
195  Spectrum selPNGCVN_All_NucInElasticProb_plus1sigma("selPNGCVN_All_NucInElasticProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrInel_N), +1) );
196  Spectrum selSigPNGCVN_All_NucInElasticProb_plus1sigma("selSigPNGCVN_All_NucInElasticProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrInel_N), +1));
197  Spectrum selBkgPNGCVN_All_NucInElasticProb_plus1sigma("selBkgPNGCVN_All_NucInElasticProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrInel_N), +1));
198 
199 
200 
201  Spectrum selPNGCVN_All_NucAbsorptionProb_plus1sigma("selPNGCVN_All_NucAbsorptionProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrAbs_N), +1));
202  Spectrum selSigPNGCVN_All_NucAbsorptionProb_plus1sigma("selSigPNGCVN_All_NucAbsorptionProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrAbs_N), +1) );
203  Spectrum selBkgPNGCVN_All_NucAbsorptionProb_plus1sigma("selBkgPNGCVN_All_NucAbsorptionProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrAbs_N), +1));
204 
205  Spectrum selPNGCVN_All_NucProductionProb_plus1sigma("selPNGCVN_All_NucProductionProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrPiProd_N), +1));
206  Spectrum selSigPNGCVN_All_NucProductionProb_plus1sigma("selSigPNGCVN_All_NucProductionProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrPiProd_N), +1));
207  Spectrum selBkgPNGCVN_All_NucProductionProb_plus1sigma("selBkgPNGCVN_All_NucProductionProb_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrPiProd_N), +1));
208 
209 
210  Spectrum selPNGCVN_FermiMom_plus1sigma("selPNGCVN_FermiMom_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightCCQEPauliSupViaKF), +1));
211  Spectrum selSigPNGCVN_FermiMom_plus1sigma("selSigPNGCVN_FermiMom_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightCCQEPauliSupViaKF), +1));
212  Spectrum selBkgPNGCVN_FermiMom_plus1sigma("selBkgPNGCVN_FermiMom_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightCCQEPauliSupViaKF), +1));
213 
214  Spectrum selPNGCVN_SpectralFn_plus1sigma("selPNGCVN_SpectralFn_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightCCQEMomDistroFGtoSF), +1));
215  Spectrum selSigPNGCVN_SpectralFn_plus1sigma("selSigPNGCVN_SpectralFn_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightCCQEMomDistroFGtoSF), +1));
216  Spectrum selBkgPNGCVN_SpectralFn_plus1sigma("selBkgPNGCVN_SpectralFn_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightCCQEMomDistroFGtoSF), +1));
217 
218 
219  Spectrum selPNGCVN_All_GammaBR_plus1sigma("selPNGCVN_All_GammaBR_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightBR1gamma), +1));
220  Spectrum selSigPNGCVN_All_GammaBR_plus1sigma("selSigPNGCVN_All_GammaBR_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightBR1gamma), +1) );
221  Spectrum selBkgPNGCVN_All_GammaBR_plus1sigma("selBkgPNGCVN_All_GammaBR_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightBR1gamma), +1));
222 
223  Spectrum selPNGCVN_All_EtaBR_plus1sigma("selPNGCVN_All_EtaBR_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightBR1eta), +1));
224  Spectrum selSigPNGCVN_All_EtaBR_plus1sigma("selSigPNGCVN_All_EtaBR_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightBR1eta), +1) );
225  Spectrum selBkgPNGCVN_All_EtaBR_plus1sigma("selBkgPNGCVN_All_EtaBR_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightBR1eta), +1));
226 
227 
228  Spectrum selPNGCVN_All_PiAngularDist_plus1sigma("selPNGCVN_All_PiAngularDist_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightTheta_Delta2Npi), +1) );
229  Spectrum selSigPNGCVN_All_PiAngularDist_plus1sigma("selSigPNGCVN_All_PiAngularDist_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightTheta_Delta2Npi), +1));
230  Spectrum selBkgPNGCVN_All_PiAngularDist_plus1sigma("selBkgPNGCVN_All_PiAngularDist_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightTheta_Delta2Npi), +1) );
231 
232 
233  Spectrum selPNGCVN_xF_plus1sigma("selPNGCVN_xF_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightAGKY_xF1pi), +1));
234  Spectrum selSigPNGCVN_xF_plus1sigma("selSigPNGCVN_xF_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightAGKY_xF1pi), +1));
235  Spectrum selBkgPNGCVN_xF_plus1sigma("selBkgPNGCVN_xF_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightAGKY_xF1pi), +1));
236 
237 
238  Spectrum selPNGCVN_pT_plus1sigma("selPNGCVN_pT_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightAGKY_pT1pi), +1));
239  Spectrum selSigPNGCVN_pT_plus1sigma("selSigPNGCVN_pT_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightAGKY_pT1pi), +1));
240  Spectrum selBkgPNGCVN_pT_plus1sigma("selBkgPNGCVN_pT_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightAGKY_pT1pi), +1));
241 
242  Spectrum selPNGCVN_FormZone_plus1sigma("selPNGCVN_FormZone_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightFormZone), +1));
243  Spectrum selSigPNGCVN_FormZone_plus1sigma("selSigPNGCVN_FormZone_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFormZone), +1));
244  Spectrum selBkgPNGCVN_FormZone_plus1sigma("selBkgPNGCVN_FormZone_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFormZone), +1));
245 
246 
247 
248  //Nonres channels (NC)
249 
250 Spectrum selPNGCVN_All_1piNonRES_nup_plus1sigma("selPNGCVN_All_1piNonRES_nup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpNC1pi), +1) );
251  Spectrum selSigPNGCVN_All_1piNonRES_nup_plus1sigma("selSigPNGCVN_All_1piNonRES_nup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpNC1pi), +1));
252  Spectrum selBkgPNGCVN_All_1piNonRES_nup_plus1sigma("selBkgPNGCVN_All_1piNonRES_nup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpNC1pi), +1) );
253 
254 Spectrum selPNGCVN_All_2piNonRES_nup_plus1sigma("selPNGCVN_All_2piNonRES_nup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpNC2pi), +1) );
255  Spectrum selSigPNGCVN_All_2piNonRES_nup_plus1sigma("selSigPNGCVN_All_2piNonRES_nup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpNC2pi), +1));
256  Spectrum selBkgPNGCVN_All_2piNonRES_nup_plus1sigma("selBkgPNGCVN_All_2piNonRES_nup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpNC2pi), +1) );
257 
258 
259 Spectrum selPNGCVN_All_1piNonRES_nun_plus1sigma("selPNGCVN_All_1piNonRES_nun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnNC1pi), +1) );
260  Spectrum selSigPNGCVN_All_1piNonRES_nun_plus1sigma("selSigPNGCVN_All_1piNonRES_nun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnNC1pi), +1));
261  Spectrum selBkgPNGCVN_All_1piNonRES_nun_plus1sigma("selBkgPNGCVN_All_1piNonRES_nun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnNC1pi), +1) );
262 
263 Spectrum selPNGCVN_All_2piNonRES_nun_plus1sigma("selPNGCVN_All_2piNonRES_nun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnNC2pi), +1) );
264  Spectrum selSigPNGCVN_All_2piNonRES_nun_plus1sigma("selSigPNGCVN_All_2piNonRES_nun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnNC2pi), +1));
265  Spectrum selBkgPNGCVN_All_2piNonRES_nun_plus1sigma("selBkgPNGCVN_All_2piNonRES_nun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnNC2pi), +1) );
266 
267 
268 
269 
270 Spectrum selPNGCVN_All_1piNonRES_antinup_plus1sigma("selPNGCVN_All_1piNonRES_antinup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpNC1pi), +1) );
271  Spectrum selSigPNGCVN_All_1piNonRES_antinup_plus1sigma("selSigPNGCVN_All_1piNonRES_antinup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpNC1pi), +1));
272  Spectrum selBkgPNGCVN_All_1piNonRES_antinup_plus1sigma("selBkgPNGCVN_All_1piNonRES_antinup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpNC1pi), +1) );
273 
274 Spectrum selPNGCVN_All_2piNonRES_antinup_plus1sigma("selPNGCVN_All_2piNonRES_antinup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpNC2pi), +1) );
275  Spectrum selSigPNGCVN_All_2piNonRES_antinup_plus1sigma("selSigPNGCVN_All_2piNonRES_antinup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpNC2pi), +1));
276  Spectrum selBkgPNGCVN_All_2piNonRES_antinup_plus1sigma("selBkgPNGCVN_All_2piNonRES_antinup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpNC2pi), +1) );
277 
278 
279 Spectrum selPNGCVN_All_1piNonRES_antinun_plus1sigma("selPNGCVN_All_1piNonRES_antinun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnNC1pi), +1) );
280  Spectrum selSigPNGCVN_All_1piNonRES_antinun_plus1sigma("selSigPNGCVN_All_1piNonRES_antinun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnNC1pi), +1));
281  Spectrum selBkgPNGCVN_All_1piNonRES_antinun_plus1sigma("selBkgPNGCVN_All_1piNonRES_antinun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnNC1pi), +1) );
282 
283 Spectrum selPNGCVN_All_2piNonRES_antinun_plus1sigma("selPNGCVN_All_2piNonRES_antinun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnNC2pi), +1) );
284  Spectrum selSigPNGCVN_All_2piNonRES_antinun_plus1sigma("selSigPNGCVN_All_2piNonRES_antinun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnNC2pi), +1));
285  Spectrum selBkgPNGCVN_All_2piNonRES_antinun_plus1sigma("selBkgPNGCVN_All_2piNonRES_antinun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnNC2pi), +1) );
286 
287 
288  //Nonres channels (CC)
289 
290 Spectrum selPNGCVN_All_CC_1piNonRES_nup_plus1sigma("selPNGCVN_All_CC_1piNonRES_nup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpCC1pi), +1) );
291  Spectrum selSigPNGCVN_All_CC_1piNonRES_nup_plus1sigma("selSigPNGCVN_All_CC_1piNonRES_nup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpCC1pi), +1));
292  Spectrum selBkgPNGCVN_All_CC_1piNonRES_nup_plus1sigma("selBkgPNGCVN_All_CC_1piNonRES_nup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpCC1pi), +1) );
293 
294 Spectrum selPNGCVN_All_CC_2piNonRES_nup_plus1sigma("selPNGCVN_All_CC_2piNonRES_nup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpCC2pi), +1) );
295  Spectrum selSigPNGCVN_All_CC_2piNonRES_nup_plus1sigma("selSigPNGCVN_All_CC_2piNonRES_nup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpCC2pi), +1));
296  Spectrum selBkgPNGCVN_All_CC_2piNonRES_nup_plus1sigma("selBkgPNGCVN_All_CC_2piNonRES_nup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpCC2pi), +1) );
297 
298 
299 Spectrum selPNGCVN_All_CC_1piNonRES_nun_plus1sigma("selPNGCVN_All_CC_1piNonRES_nun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnCC1pi), +1) );
300  Spectrum selSigPNGCVN_All_CC_1piNonRES_nun_plus1sigma("selSigPNGCVN_All_CC_1piNonRES_nun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnCC1pi), +1));
301  Spectrum selBkgPNGCVN_All_CC_1piNonRES_nun_plus1sigma("selBkgPNGCVN_All_CC_1piNonRES_nun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnCC1pi), +1) );
302 
303 Spectrum selPNGCVN_All_CC_2piNonRES_nun_plus1sigma("selPNGCVN_All_CC_2piNonRES_nun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnCC2pi), +1) );
304  Spectrum selSigPNGCVN_All_CC_2piNonRES_nun_plus1sigma("selSigPNGCVN_All_CC_2piNonRES_nun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnCC2pi), +1));
305  Spectrum selBkgPNGCVN_All_CC_2piNonRES_nun_plus1sigma("selBkgPNGCVN_All_CC_2piNonRES_nun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnCC2pi), +1) );
306 
307 
308 
309 
310 Spectrum selPNGCVN_All_CC_1piNonRES_antinup_plus1sigma("selPNGCVN_All_CC_1piNonRES_antinup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpCC1pi), +1) );
311  Spectrum selSigPNGCVN_All_CC_1piNonRES_antinup_plus1sigma("selSigPNGCVN_All_CC_1piNonRES_antinup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpCC1pi), +1));
312  Spectrum selBkgPNGCVN_All_CC_1piNonRES_antinup_plus1sigma("selBkgPNGCVN_All_CC_1piNonRES_antinup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpCC1pi), +1) );
313 
314 Spectrum selPNGCVN_All_CC_2piNonRES_antinup_plus1sigma("selPNGCVN_All_CC_2piNonRES_antinup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpCC2pi), +1) );
315  Spectrum selSigPNGCVN_All_CC_2piNonRES_antinup_plus1sigma("selSigPNGCVN_All_CC_2piNonRES_antinup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpCC2pi), +1));
316  Spectrum selBkgPNGCVN_All_CC_2piNonRES_antinup_plus1sigma("selBkgPNGCVN_All_CC_2piNonRES_antinup_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpCC2pi), +1) );
317 
318 
319 Spectrum selPNGCVN_All_CC_1piNonRES_antinun_plus1sigma("selPNGCVN_All_CC_1piNonRES_antinun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnCC1pi), +1) );
320  Spectrum selSigPNGCVN_All_CC_1piNonRES_antinun_plus1sigma("selSigPNGCVN_All_CC_1piNonRES_antinun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnCC1pi), +1));
321  Spectrum selBkgPNGCVN_All_CC_1piNonRES_antinun_plus1sigma("selBkgPNGCVN_All_CC_1piNonRES_antinun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnCC1pi), +1) );
322 
323 Spectrum selPNGCVN_All_CC_2piNonRES_antinun_plus1sigma("selPNGCVN_All_CC_2piNonRES_antinun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnCC2pi), +1) );
324  Spectrum selSigPNGCVN_All_CC_2piNonRES_antinun_plus1sigma("selSigPNGCVN_All_CC_2piNonRES_antinun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnCC2pi), +1));
325  Spectrum selBkgPNGCVN_All_CC_2piNonRES_antinun_plus1sigma("selBkgPNGCVN_All_CC_2piNonRES_antinun_plus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnCC2pi), +1) );
326 
327 
328 
329  //Minus sigma
330 
331 
332 
333  Spectrum selPNGCVN_ShapeNorm_AxialMassNCEL_minus1sigma("selPNGCVN_ShapeNorm_AxialMassNCEL_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaNCEL), -1) );
334  Spectrum selSigPNGCVN_ShapeNorm_AxialMassNCEL_minus1sigma("selSigPNGCVN_ShapeNorm_AxialMassNCEL_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaNCEL), -1) );
335  Spectrum selBkgPNGCVN_ShapeNorm_AxialMassNCEL_minus1sigma("selBkgPNGCVN_ShapeNorm_AxialMassNCEL_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaNCEL), -1) );
336 
337 
338  Spectrum selPNGCVN_ShapeNorm_EtaNCEL_minus1sigma("selPNGCVN_ShapeNorm_EtaNCEL_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightEtaNCEL), -1) );
339  Spectrum selSigPNGCVN_ShapeNorm_EtaNCEL_minus1sigma("selSigPNGCVN_ShapeNorm_EtaNCEL_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightEtaNCEL), -1) );
340  Spectrum selBkgPNGCVN_ShapeNorm_EtaNCEL_minus1sigma("selBkgPNGCVN_ShapeNorm_EtaNCEL_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightEtaNCEL), -1) );
341 
342  Spectrum selPNGCVN_ShapeNorm_AxialMassNCRES_minus1sigma("selPNGCVN_ShapeNorm_AxialMassNCRES_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaNCRES), -1));
343  Spectrum selSigPNGCVN_ShapeNorm_AxialMassNCRES_minus1sigma("selSigPNGCVN_ShapeNorm_AxialMassNCRES_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaNCRES), -1) );
344  Spectrum selBkgPNGCVN_ShapeNorm_AxialMassNCRES_minus1sigma("selBkgPNGCVN_ShapeNorm_AxialMassNCRES_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaNCRES), -1));
345 
346 
347  Spectrum selPNGCVN_ShapeNorm_VectorMassNCRES_minus1sigma("selPNGCVN_ShapeNorm_VectorMassNCRES_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMvNCRES), -1) );
348  Spectrum selSigPNGCVN_ShapeNorm_VectorMassNCRES_minus1sigma("selSigPNGCVN_ShapeNorm_VectorMassNCRES_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightMvNCRES), -1) );
349  Spectrum selBkgPNGCVN_ShapeNorm_VectorMassNCRES_minus1sigma("selBkgPNGCVN_ShapeNorm_VectorMassNCRES_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMvNCRES), -1));
350 
351 
352  Spectrum selPNGCVN_ShapeNorm_AxialMassCCRES_minus1sigma("selPNGCVN_ShapeNorm_AxialMassCCRES_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaCCRES), -1));
353  Spectrum selSigPNGCVN_ShapeNorm_AxialMassCCRES_minus1sigma("selSigPNGCVN_ShapeNorm_AxialMassCCRES_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaCCRES), -1) );
354  Spectrum selBkgPNGCVN_ShapeNorm_AxialMassCCRES_minus1sigma("selBkgPNGCVN_ShapeNorm_AxialMassCCRES_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaCCRES), -1) );
355 
356  Spectrum selPNGCVN_ShapeNorm_VectorMassCCRES_minus1sigma("selPNGCVN_ShapeNorm_VectorMassCCRES_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMvCCRES), -1) );
357  Spectrum selSigPNGCVN_ShapeNorm_VectorMassCCRES_minus1sigma("selSigPNGCVN_ShapeNorm_VectorMassCCRES_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMvCCRES), -1) );
358  Spectrum selBkgPNGCVN_ShapeNorm_VectorMassCCRES_minus1sigma("selBkgPNGCVN_ShapeNorm_VectorMassCCRES_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMvCCRES), -1) );
359 
360 
361  Spectrum selPNGCVN_All_AxialMassCOH_minus1sigma("selPNGCVN_All_AxialMassCOH_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaCOHpi), -1));
362  Spectrum selSigPNGCVN_All_AxialMassCOH_minus1sigma("selSigPNGCVN_All_AxialMassCOH_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaCOHpi), -1));
363  Spectrum selBkgPNGCVN_All_AxialMassCOH_minus1sigma("selBkgPNGCVN_All_AxialMassCOH_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMaCOHpi), -1));
364 
365 
366  Spectrum selPNGCVN_All_NuclearSizeCOH_minus1sigma("selPNGCVN_All_NuclearSizeCOH_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightR0COHpi), -1) );
367  Spectrum selSigPNGCVN_All_NuclearSizeCOH_minus1sigma("selSigPNGCVN_All_NuclearSizeCOH_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightR0COHpi), -1));
368  Spectrum selBkgPNGCVN_All_NuclearSizeCOH_minus1sigma("selBkgPNGCVN_All_NuclearSizeCOH_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightR0COHpi), -1));
369 
370 
371  Spectrum selPNGCVN_ShapeNorm_DIS_Aht_minus1sigma("selPNGCVN_ShapeNorm_DIS_Aht_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightAhtBY), -1));
372  Spectrum selSigPNGCVN_ShapeNorm_DIS_Aht_minus1sigma("selSigPNGCVN_ShapeNorm_DIS_Aht_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightAhtBY), -1));
373  Spectrum selBkgPNGCVN_ShapeNorm_DIS_Aht_minus1sigma("selBkgPNGCVN_ShapeNorm_DIS_Aht_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightAhtBY), -1));
374 
375 
376  Spectrum selPNGCVN_ShapeNorm_DIS_Bht_minus1sigma("selPNGCVN_ShapeNorm_DIS_Bht_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightBhtBY), -1) );
377  Spectrum selSigPNGCVN_ShapeNorm_DIS_Bht_minus1sigma("selSigPNGCVN_ShapeNorm_DIS_Bht_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightBhtBY), -1) );
378  Spectrum selBkgPNGCVN_ShapeNorm_DIS_Bht_minus1sigma("selBkgPNGCVN_ShapeNorm_DIS_Bht_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightBhtBY), -1));
379 
380 
381  Spectrum selPNGCVN_ShapeNorm_DIS_CV1u_minus1sigma("selPNGCVN_ShapeNorm_DIS_CV1u_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightCV1uBY), -1) );
382  Spectrum selSigPNGCVN_ShapeNorm_DIS_CV1u_minus1sigma("selSigPNGCVN_ShapeNorm_DIS_CV1u_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightCV1uBY), -1) );
383  Spectrum selBkgPNGCVN_ShapeNorm_DIS_CV1u_minus1sigma("selBkgPNGCVN_ShapeNorm_DIS_CV1u_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightCV1uBY), -1));
384 
385 
386  Spectrum selPNGCVN_ShapeNorm_DIS_CV2u_minus1sigma("selPNGCVN_ShapeNorm_DIS_CV2u_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightCV2uBY), -1) );
387  Spectrum selSigPNGCVN_ShapeNorm_DIS_CV2u_minus1sigma("selSigPNGCVN_ShapeNorm_DIS_CV2u_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightCV2uBY), -1));
388  Spectrum selBkgPNGCVN_ShapeNorm_DIS_CV2u_minus1sigma("selBkgPNGCVN_ShapeNorm_DIS_CV2u_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightCV2uBY), -1));
389 
390 
391 
392  Spectrum selPNGCVN_mfp_pi_minus1sigma("selPNGCVN_mfp_pi_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMFP_pi), -1));
393  Spectrum selSigPNGCVN_mfp_pi_minus1sigma("selSigPNGCVN_mfp_pi_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMFP_pi), -1) );
394  Spectrum selBkgPNGCVN_mfp_pi_minus1sigma("selBkgPNGCVN_mfp_pi_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMFP_pi), -1));
395 
396  Spectrum selPNGCVN_mfp_N_minus1sigma("selPNGCVN_mfp_N_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMFP_N), -1));
397  Spectrum selSigPNGCVN_mfp_N_minus1sigma("selSigPNGCVN_mfp_N_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMFP_N), -1) );
398  Spectrum selBkgPNGCVN_mfp_N_minus1sigma("selBkgPNGCVN_mfp_N_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightMFP_N), -1));
399 
400  Spectrum selPNGCVN_All_PiChargeExchangeProb_minus1sigma("selPNGCVN_All_PiChargeExchangeProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrCEx_pi), -1) );
401  Spectrum selSigPNGCVN_All_PiChargeExchangeProb_minus1sigma("selSigPNGCVN_All_PiChargeExchangeProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrCEx_pi), -1) );
402  Spectrum selBkgPNGCVN_All_PiChargeExchangeProb_minus1sigma("selBkgPNGCVN_All_PiChargeExchangeProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrCEx_pi), -1));
403 
404  Spectrum selPNGCVN_All_PiElasticProb_minus1sigma("selPNGCVN_All_PiElasticProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrElas_pi), -1) );
405  Spectrum selSigPNGCVN_All_PiElasticProb_minus1sigma("selSigPNGCVN_All_PiElasticProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrElas_pi), -1) );
406  Spectrum selBkgPNGCVN_All_PiElasticProb_minus1sigma("selBkgPNGCVN_All_PiElasticProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrElas_pi), -1) );
407 
408  Spectrum selPNGCVN_All_PiInElasticProb_minus1sigma("selPNGCVN_All_PiInElasticProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrInel_pi), -1) );
409  Spectrum selSigPNGCVN_All_PiInElasticProb_minus1sigma("selSigPNGCVN_All_PiInElasticProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrInel_pi), -1));
410  Spectrum selBkgPNGCVN_All_PiInElasticProb_minus1sigma("selBkgPNGCVN_All_PiInElasticProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrInel_pi), -1));
411 
412 
413 
414  Spectrum selPNGCVN_All_PiAbsorptionProb_minus1sigma("selPNGCVN_All_PiAbsorptionProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrAbs_pi), -1));
415  Spectrum selSigPNGCVN_All_PiAbsorptionProb_minus1sigma("selSigPNGCVN_All_PiAbsorptionProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrAbs_pi), -1) );
416  Spectrum selBkgPNGCVN_All_PiAbsorptionProb_minus1sigma("selBkgPNGCVN_All_PiAbsorptionProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrAbs_pi), -1));
417 
418  Spectrum selPNGCVN_All_PiProductionProb_minus1sigma("selPNGCVN_All_PiProductionProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrPiProd_pi), -1));
419  Spectrum selSigPNGCVN_All_PiProductionProb_minus1sigma("selSigPNGCVN_All_PiProductionProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrPiProd_pi), -1));
420  Spectrum selBkgPNGCVN_All_PiProductionProb_minus1sigma("selBkgPNGCVN_All_PiProductionProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrPiProd_pi), -1));
421 
422 
423  Spectrum selPNGCVN_All_NucChargeExchangeProb_minus1sigma("selPNGCVN_All_NucChargeExchangeProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrCEx_N), -1) );
424  Spectrum selSigPNGCVN_All_NucChargeExchangeProb_minus1sigma("selSigPNGCVN_All_NucChargeExchangeProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrCEx_N), -1) );
425  Spectrum selBkgPNGCVN_All_NucChargeExchangeProb_minus1sigma("selBkgPNGCVN_All_NucChargeExchangeProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrCEx_N), -1));
426 
427 
428  Spectrum selPNGCVN_All_NucElasticProb_minus1sigma("selPNGCVN_All_NucElasticProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrElas_N), -1) );
429  Spectrum selSigPNGCVN_All_NucElasticProb_minus1sigma("selSigPNGCVN_All_NucElasticProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrElas_N), -1) );
430  Spectrum selBkgPNGCVN_All_NucElasticProb_minus1sigma("selBkgPNGCVN_All_NucElasticProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrElas_N), -1) );
431 
432  Spectrum selPNGCVN_All_NucInElasticProb_minus1sigma("selPNGCVN_All_NucInElasticProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrInel_N), -1) );
433  Spectrum selSigPNGCVN_All_NucInElasticProb_minus1sigma("selSigPNGCVN_All_NucInElasticProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrInel_N), -1));
434  Spectrum selBkgPNGCVN_All_NucInElasticProb_minus1sigma("selBkgPNGCVN_All_NucInElasticProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrInel_N), -1));
435 
436 
437 
438  Spectrum selPNGCVN_All_NucAbsorptionProb_minus1sigma("selPNGCVN_All_NucAbsorptionProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrAbs_N), -1));
439  Spectrum selSigPNGCVN_All_NucAbsorptionProb_minus1sigma("selSigPNGCVN_All_NucAbsorptionProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrAbs_N), -1) );
440  Spectrum selBkgPNGCVN_All_NucAbsorptionProb_minus1sigma("selBkgPNGCVN_All_NucAbsorptionProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrAbs_N), -1));
441 
442  Spectrum selPNGCVN_All_NucProductionProb_minus1sigma("selPNGCVN_All_NucProductionProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrPiProd_N), -1));
443  Spectrum selSigPNGCVN_All_NucProductionProb_minus1sigma("selSigPNGCVN_All_NucProductionProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrPiProd_N), -1));
444  Spectrum selBkgPNGCVN_All_NucProductionProb_minus1sigma("selBkgPNGCVN_All_NucProductionProb_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFrPiProd_N), -1));
445 
446 
447  Spectrum selPNGCVN_FermiMom_minus1sigma("selPNGCVN_FermiMom_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightCCQEPauliSupViaKF), -1));
448  Spectrum selSigPNGCVN_FermiMom_minus1sigma("selSigPNGCVN_FermiMom_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightCCQEPauliSupViaKF), -1));
449  Spectrum selBkgPNGCVN_FermiMom_minus1sigma("selBkgPNGCVN_FermiMom_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightCCQEPauliSupViaKF), -1));
450 
451  Spectrum selPNGCVN_SpectralFn_minus1sigma("selPNGCVN_SpectralFn_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightCCQEMomDistroFGtoSF), -1));
452  Spectrum selSigPNGCVN_SpectralFn_minus1sigma("selSigPNGCVN_SpectralFn_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightCCQEMomDistroFGtoSF), -1));
453  Spectrum selBkgPNGCVN_SpectralFn_minus1sigma("selBkgPNGCVN_SpectralFn_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightCCQEMomDistroFGtoSF), -1));
454 
455 
456  Spectrum selPNGCVN_All_GammaBR_minus1sigma("selPNGCVN_All_GammaBR_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightBR1gamma), -1));
457  Spectrum selSigPNGCVN_All_GammaBR_minus1sigma("selSigPNGCVN_All_GammaBR_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightBR1gamma), -1) );
458  Spectrum selBkgPNGCVN_All_GammaBR_minus1sigma("selBkgPNGCVN_All_GammaBR_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightBR1gamma), -1));
459 
460  Spectrum selPNGCVN_All_EtaBR_minus1sigma("selPNGCVN_All_EtaBR_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightBR1eta), -1));
461  Spectrum selSigPNGCVN_All_EtaBR_minus1sigma("selSigPNGCVN_All_EtaBR_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightBR1eta), -1) );
462  Spectrum selBkgPNGCVN_All_EtaBR_minus1sigma("selBkgPNGCVN_All_EtaBR_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightBR1eta), -1));
463 
464 
465  Spectrum selPNGCVN_All_PiAngularDist_minus1sigma("selPNGCVN_All_PiAngularDist_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightTheta_Delta2Npi), -1) );
466  Spectrum selSigPNGCVN_All_PiAngularDist_minus1sigma("selSigPNGCVN_All_PiAngularDist_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightTheta_Delta2Npi), -1));
467  Spectrum selBkgPNGCVN_All_PiAngularDist_minus1sigma("selBkgPNGCVN_All_PiAngularDist_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightTheta_Delta2Npi), -1 ) );
468 
469 
470  Spectrum selPNGCVN_xF_minus1sigma("selPNGCVN_xF_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightAGKY_xF1pi), -1));
471  Spectrum selSigPNGCVN_xF_minus1sigma("selSigPNGCVN_xF_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightAGKY_xF1pi), -1));
472  Spectrum selBkgPNGCVN_xF_minus1sigma("selBkgPNGCVN_xF_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightAGKY_xF1pi), -1));
473 
474 
475  Spectrum selPNGCVN_pT_minus1sigma("selPNGCVN_pT_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightAGKY_pT1pi), -1));
476  Spectrum selSigPNGCVN_pT_minus1sigma("selSigPNGCVN_pT_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightAGKY_pT1pi), -1));
477  Spectrum selBkgPNGCVN_pT_minus1sigma("selBkgPNGCVN_pT_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightAGKY_pT1pi), -1));
478 
479  Spectrum selPNGCVN_FormZone_minus1sigma("selPNGCVN_FormZone_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightFormZone), -1));
480  Spectrum selSigPNGCVN_FormZone_minus1sigma("selSigPNGCVN_FormZone_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFormZone), -1));
481  Spectrum selBkgPNGCVN_FormZone_minus1sigma("selBkgPNGCVN_FormZone_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightFormZone), -1));
482 
483  //non res
484 
485 
486 Spectrum selPNGCVN_All_1piNonRES_nup_minus1sigma("selPNGCVN_All_1piNonRES_nup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpNC1pi), -1) );
487  Spectrum selSigPNGCVN_All_1piNonRES_nup_minus1sigma("selSigPNGCVN_All_1piNonRES_nup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpNC1pi), -1));
488  Spectrum selBkgPNGCVN_All_1piNonRES_nup_minus1sigma("selBkgPNGCVN_All_1piNonRES_nup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpNC1pi), -1) );
489 
490 Spectrum selPNGCVN_All_2piNonRES_nup_minus1sigma("selPNGCVN_All_2piNonRES_nup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpNC2pi), -1) );
491  Spectrum selSigPNGCVN_All_2piNonRES_nup_minus1sigma("selSigPNGCVN_All_2piNonRES_nup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpNC2pi), -1));
492  Spectrum selBkgPNGCVN_All_2piNonRES_nup_minus1sigma("selBkgPNGCVN_All_2piNonRES_nup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpNC2pi), -1) );
493 
494 
495 Spectrum selPNGCVN_All_1piNonRES_nun_minus1sigma("selPNGCVN_All_1piNonRES_nun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnNC1pi), -1) );
496  Spectrum selSigPNGCVN_All_1piNonRES_nun_minus1sigma("selSigPNGCVN_All_1piNonRES_nun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnNC1pi), -1));
497  Spectrum selBkgPNGCVN_All_1piNonRES_nun_minus1sigma("selBkgPNGCVN_All_1piNonRES_nun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnNC1pi), -1) );
498 
499 Spectrum selPNGCVN_All_2piNonRES_nun_minus1sigma("selPNGCVN_All_2piNonRES_nun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnNC2pi), -1) );
500  Spectrum selSigPNGCVN_All_2piNonRES_nun_minus1sigma("selSigPNGCVN_All_2piNonRES_nun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnNC2pi), -1));
501  Spectrum selBkgPNGCVN_All_2piNonRES_nun_minus1sigma("selBkgPNGCVN_All_2piNonRES_nun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnNC2pi), -1) );
502 
503 
504 
505 
506 Spectrum selPNGCVN_All_1piNonRES_antinup_minus1sigma("selPNGCVN_All_1piNonRES_antinup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpNC1pi), -1) );
507  Spectrum selSigPNGCVN_All_1piNonRES_antinup_minus1sigma("selSigPNGCVN_All_1piNonRES_antinup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpNC1pi), -1));
508  Spectrum selBkgPNGCVN_All_1piNonRES_antinup_minus1sigma("selBkgPNGCVN_All_1piNonRES_antinup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpNC1pi), -1) );
509 
510 Spectrum selPNGCVN_All_2piNonRES_antinup_minus1sigma("selPNGCVN_All_2piNonRES_antinup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpNC2pi), -1) );
511  Spectrum selSigPNGCVN_All_2piNonRES_antinup_minus1sigma("selSigPNGCVN_All_2piNonRES_antinup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpNC2pi), -1));
512  Spectrum selBkgPNGCVN_All_2piNonRES_antinup_minus1sigma("selBkgPNGCVN_All_2piNonRES_antinup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpNC2pi), -1) );
513 
514 
515 Spectrum selPNGCVN_All_1piNonRES_antinun_minus1sigma("selPNGCVN_All_1piNonRES_antinun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnNC1pi), -1) );
516  Spectrum selSigPNGCVN_All_1piNonRES_antinun_minus1sigma("selSigPNGCVN_All_1piNonRES_antinun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnNC1pi), -1));
517  Spectrum selBkgPNGCVN_All_1piNonRES_antinun_minus1sigma("selBkgPNGCVN_All_1piNonRES_antinun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnNC1pi), -1) );
518 
519 Spectrum selPNGCVN_All_2piNonRES_antinun_minus1sigma("selPNGCVN_All_2piNonRES_antinun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnNC2pi), -1) );
520  Spectrum selSigPNGCVN_All_2piNonRES_antinun_minus1sigma("selSigPNGCVN_All_2piNonRES_antinun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnNC2pi), -1));
521  Spectrum selBkgPNGCVN_All_2piNonRES_antinun_minus1sigma("selBkgPNGCVN_All_2piNonRES_antinun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnNC2pi), -1) );
522 
523 
524  //Nonres channels (CC)
525 
526 Spectrum selPNGCVN_All_CC_1piNonRES_nup_minus1sigma("selPNGCVN_All_CC_1piNonRES_nup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpCC1pi), -1) );
527  Spectrum selSigPNGCVN_All_CC_1piNonRES_nup_minus1sigma("selSigPNGCVN_All_CC_1piNonRES_nup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpCC1pi), -1));
528  Spectrum selBkgPNGCVN_All_CC_1piNonRES_nup_minus1sigma("selBkgPNGCVN_All_CC_1piNonRES_nup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpCC1pi), -1) );
529 
530 Spectrum selPNGCVN_All_CC_2piNonRES_nup_minus1sigma("selPNGCVN_All_CC_2piNonRES_nup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpCC2pi), -1) );
531  Spectrum selSigPNGCVN_All_CC_2piNonRES_nup_minus1sigma("selSigPNGCVN_All_CC_2piNonRES_nup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpCC2pi), -1));
532  Spectrum selBkgPNGCVN_All_CC_2piNonRES_nup_minus1sigma("selBkgPNGCVN_All_CC_2piNonRES_nup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvpCC2pi), -1) );
533 
534 
535 Spectrum selPNGCVN_All_CC_1piNonRES_nun_minus1sigma("selPNGCVN_All_CC_1piNonRES_nun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnCC1pi), -1) );
536  Spectrum selSigPNGCVN_All_CC_1piNonRES_nun_minus1sigma("selSigPNGCVN_All_CC_1piNonRES_nun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnCC1pi), -1));
537  Spectrum selBkgPNGCVN_All_CC_1piNonRES_nun_minus1sigma("selBkgPNGCVN_All_CC_1piNonRES_nun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnCC1pi), -1) );
538 
539 Spectrum selPNGCVN_All_CC_2piNonRES_nun_minus1sigma("selPNGCVN_All_CC_2piNonRES_nun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnCC2pi), -1) );
540  Spectrum selSigPNGCVN_All_CC_2piNonRES_nun_minus1sigma("selSigPNGCVN_All_CC_2piNonRES_nun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnCC2pi), -1));
541  Spectrum selBkgPNGCVN_All_CC_2piNonRES_nun_minus1sigma("selBkgPNGCVN_All_CC_2piNonRES_nun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvnCC2pi), -1) );
542 
543 
544 
545 
546 Spectrum selPNGCVN_All_CC_1piNonRES_antinup_minus1sigma("selPNGCVN_All_CC_1piNonRES_antinup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpCC1pi), -1) );
547  Spectrum selSigPNGCVN_All_CC_1piNonRES_antinup_minus1sigma("selSigPNGCVN_All_CC_1piNonRES_antinup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpCC1pi), -1));
548  Spectrum selBkgPNGCVN_All_CC_1piNonRES_antinup_minus1sigma("selBkgPNGCVN_All_CC_1piNonRES_antinup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpCC1pi), -1) );
549 
550 Spectrum selPNGCVN_All_CC_2piNonRES_antinup_minus1sigma("selPNGCVN_All_CC_2piNonRES_antinup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpCC2pi), -1) );
551  Spectrum selSigPNGCVN_All_CC_2piNonRES_antinup_minus1sigma("selSigPNGCVN_All_CC_2piNonRES_antinup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpCC2pi), -1));
552  Spectrum selBkgPNGCVN_All_CC_2piNonRES_antinup_minus1sigma("selBkgPNGCVN_All_CC_2piNonRES_antinup_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarpCC2pi), -1) );
553 
554 
555 Spectrum selPNGCVN_All_CC_1piNonRES_antinun_minus1sigma("selPNGCVN_All_CC_1piNonRES_antinun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnCC1pi), -1) );
556  Spectrum selSigPNGCVN_All_CC_1piNonRES_antinun_minus1sigma("selSigPNGCVN_All_CC_1piNonRES_antinun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnCC1pi), -1));
557  Spectrum selBkgPNGCVN_All_CC_1piNonRES_antinun_minus1sigma("selBkgPNGCVN_All_CC_1piNonRES_antinun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnCC1pi), -1) );
558 
559 Spectrum selPNGCVN_All_CC_2piNonRES_antinun_minus1sigma("selPNGCVN_All_CC_2piNonRES_antinun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnCC2pi), -1) );
560  Spectrum selSigPNGCVN_All_CC_2piNonRES_antinun_minus1sigma("selSigPNGCVN_All_CC_2piNonRES_antinun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsSignal && kIsPreSel,SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnCC2pi), -1));
561  Spectrum selBkgPNGCVN_All_CC_2piNonRES_antinun_minus1sigma("selBkgPNGCVN_All_CC_2piNonRES_antinun_minus1sigma", PCIDbins, loader, kNpCVN1gamma, kIsTotBkg && kIsPreSel, SystShifts(GetGenieRwgtSyst(rwgt::fReweightRvbarnCC2pi), -1) );
562 
563 
564 
565 loader.Go();
566 
567  // const double pot = 8.08988e+20; // for the data Prod3 all files
568 
569  std::string fname = "GENIEKnobseffectCVN02.root";
570 
571  TFile fout(fname.c_str(),"RECREATE");
572 
573 
574  saveS(selPNGCVN, "selPNGCVN" );
575  saveS(selSigPNGCVN, "selSigPNGCVN" );
576  saveS(selBkgPNGCVN, "selBkgPNGCVN" );
577 
578 
579 
580  saveS(selPNGCVN_ShapeNorm_AxialMassNCEL_plus1sigma, "selPNGCVN_ShapeNorm_AxialMassNCEL_plus1sigma" );
581  saveS(selSigPNGCVN_ShapeNorm_AxialMassNCEL_plus1sigma, "selSigPNGCVN_ShapeNorm_AxialMassNCEL_plus1sigma" );
582  saveS(selBkgPNGCVN_ShapeNorm_AxialMassNCEL_plus1sigma, "selBkgPNGCVN_ShapeNorm_AxialMassNCEL_plus1sigma" );
583 
584  saveS(selPNGCVN_ShapeNorm_EtaNCEL_plus1sigma, "selPNGCVN_ShapeNorm_EtaNCEL_plus1sigma" );
585  saveS(selSigPNGCVN_ShapeNorm_EtaNCEL_plus1sigma, "selSigPNGCVN_ShapeNorm_EtaNCEL_plus1sigma" );
586  saveS(selBkgPNGCVN_ShapeNorm_EtaNCEL_plus1sigma, "selBkgPNGCVN_ShapeNorm_EtaNCEL_plus1sigma" );
587 
588  saveS(selPNGCVN_ShapeNorm_AxialMassNCRES_plus1sigma, "selPNGCVN_ShapeNorm_AxialMassNCRES_plus1sigma" );
589  saveS(selSigPNGCVN_ShapeNorm_AxialMassNCRES_plus1sigma, "selSigPNGCVN_ShapeNorm_AxialMassNCRES_plus1sigma" );
590  saveS(selBkgPNGCVN_ShapeNorm_AxialMassNCRES_plus1sigma, "selBkgPNGCVN_ShapeNorm_AxialMassNCRES_plus1sigma" );
591 
592  saveS(selPNGCVN_ShapeNorm_VectorMassNCRES_plus1sigma, "selPNGCVN_ShapeNorm_VectorMassNCRES_plus1sigma" );
593  saveS(selSigPNGCVN_ShapeNorm_VectorMassNCRES_plus1sigma, "selSigPNGCVN_ShapeNorm_VectorMassNCRES_plus1sigma" );
594  saveS(selBkgPNGCVN_ShapeNorm_VectorMassNCRES_plus1sigma, "selBkgPNGCVN_ShapeNorm_VectorMassNCRES_plus1sigma" );
595 
596 
597  saveS(selPNGCVN_ShapeNorm_AxialMassCCRES_plus1sigma, "selPNGCVN_ShapeNorm_AxialMassCCRES_plus1sigma" );
598  saveS(selSigPNGCVN_ShapeNorm_AxialMassCCRES_plus1sigma, "selSigPNGCVN_ShapeNorm_AxialMassCCRES_plus1sigma" );
599  saveS(selBkgPNGCVN_ShapeNorm_AxialMassCCRES_plus1sigma, "selBkgPNGCVN_ShapeNorm_AxialMassCCRES_plus1sigma" );
600 
601  saveS(selPNGCVN_ShapeNorm_VectorMassCCRES_plus1sigma, "selPNGCVN_ShapeNorm_VectorMassCCRES_plus1sigma" );
602  saveS(selSigPNGCVN_ShapeNorm_VectorMassCCRES_plus1sigma, "selSigPNGCVN_ShapeNorm_VectorMassCCRES_plus1sigma" );
603  saveS(selBkgPNGCVN_ShapeNorm_VectorMassCCRES_plus1sigma, "selBkgPNGCVN_ShapeNorm_VectorMassCCRES_plus1sigma" );
604 
605 
606  saveS(selPNGCVN_All_AxialMassCOH_plus1sigma, "selPNGCVN_All_AxialMassCOH_plus1sigma" );
607  saveS(selSigPNGCVN_All_AxialMassCOH_plus1sigma, "selSigPNGCVN_All_AxialMassCOH_plus1sigma" );
608  saveS(selBkgPNGCVN_All_AxialMassCOH_plus1sigma, "selBkgPNGCVN_All_AxialMassCOH_plus1sigma" );
609 
610  saveS(selPNGCVN_All_NuclearSizeCOH_plus1sigma, "selPNGCVN_All_NuclearSizeCOH_plus1sigma" );
611  saveS(selSigPNGCVN_All_NuclearSizeCOH_plus1sigma, "selSigPNGCVN_All_NuclearSizeCOH_plus1sigma" );
612  saveS(selBkgPNGCVN_All_NuclearSizeCOH_plus1sigma, "selBkgPNGCVN_All_NuclearSizeCOH_plus1sigma" );
613 
614  saveS(selPNGCVN_ShapeNorm_DIS_Aht_plus1sigma, "selPNGCVN_ShapeNorm_DIS_Aht_plus1sigma" );
615  saveS(selSigPNGCVN_ShapeNorm_DIS_Aht_plus1sigma, "selSigPNGCVN_ShapeNorm_DIS_Aht_plus1sigma" );
616  saveS(selBkgPNGCVN_ShapeNorm_DIS_Aht_plus1sigma, "selBkgPNGCVN_ShapeNorm_DIS_Aht_plus1sigma" );
617 
618  saveS(selPNGCVN_ShapeNorm_DIS_Bht_plus1sigma, "selPNGCVN_ShapeNorm_DIS_Bht_plus1sigma" );
619  saveS(selSigPNGCVN_ShapeNorm_DIS_Bht_plus1sigma, "selSigPNGCVN_ShapeNorm_DIS_Bht_plus1sigma" );
620  saveS(selBkgPNGCVN_ShapeNorm_DIS_Bht_plus1sigma, "selBkgPNGCVN_ShapeNorm_DIS_Bht_plus1sigma" );
621 
622  saveS(selPNGCVN_ShapeNorm_DIS_CV1u_plus1sigma, "selPNGCVN_ShapeNorm_DIS_CV1u_plus1sigma" );
623  saveS(selSigPNGCVN_ShapeNorm_DIS_CV1u_plus1sigma, "selSigPNGCVN_ShapeNorm_DIS_CV1u_plus1sigma" );
624  saveS(selBkgPNGCVN_ShapeNorm_DIS_CV1u_plus1sigma, "selBkgPNGCVN_ShapeNorm_DIS_CV1u_plus1sigma" );
625 
626  saveS(selPNGCVN_ShapeNorm_DIS_CV2u_plus1sigma, "selPNGCVN_ShapeNorm_DIS_CV2u_plus1sigma" );
627  saveS(selSigPNGCVN_ShapeNorm_DIS_CV2u_plus1sigma, "selSigPNGCVN_ShapeNorm_DIS_CV2u_plus1sigma" );
628  saveS(selBkgPNGCVN_ShapeNorm_DIS_CV2u_plus1sigma, "selBkgPNGCVN_ShapeNorm_DIS_CV2u_plus1sigma" );
629 
630 
631  saveS(selPNGCVN_mfp_pi_plus1sigma, "selPNGCVN_mfp_pi_plus1sigma" );
632  saveS(selSigPNGCVN_mfp_pi_plus1sigma, "selSigPNGCVN_mfp_pi_plus1sigma" );
633  saveS(selBkgPNGCVN_mfp_pi_plus1sigma, "selBkgPNGCVN_mfp_pi_plus1sigma" );
634 
635  saveS(selPNGCVN_mfp_N_plus1sigma, "selPNGCVN_mfp_N_plus1sigma" );
636  saveS(selSigPNGCVN_mfp_N_plus1sigma, "selSigPNGCVN_mfp_N_plus1sigma" );
637  saveS(selBkgPNGCVN_mfp_N_plus1sigma, "selBkgPNGCVN_mfp_N_plus1sigma" );
638 
639 
640  saveS(selPNGCVN_All_PiElasticProb_plus1sigma, "selPNGCVN_All_PiElasticProb_plus1sigma" );
641  saveS(selSigPNGCVN_All_PiElasticProb_plus1sigma, "selSigPNGCVN_All_PiElasticProb_plus1sigma" );
642  saveS(selBkgPNGCVN_All_PiElasticProb_plus1sigma, "selBkgPNGCVN_All_PiElasticProb_plus1sigma" );
643 
644  saveS(selPNGCVN_All_PiInElasticProb_plus1sigma, "selPNGCVN_All_PiInElasticProb_plus1sigma" );
645  saveS(selSigPNGCVN_All_PiInElasticProb_plus1sigma, "selSigPNGCVN_All_PiInElasticProb_plus1sigma" );
646  saveS(selBkgPNGCVN_All_PiInElasticProb_plus1sigma, "selBkgPNGCVN_All_PiInElasticProb_plus1sigma" );
647 
648 
649  saveS(selPNGCVN_All_PiChargeExchangeProb_plus1sigma, "selPNGCVN_All_PiChargeExchangeProb_plus1sigma" );
650  saveS(selSigPNGCVN_All_PiChargeExchangeProb_plus1sigma, "selSigPNGCVN_All_PiChargeExchangeProb_plus1sigma" );
651  saveS(selBkgPNGCVN_All_PiChargeExchangeProb_plus1sigma, "selBkgPNGCVN_All_PiChargeExchangeProb_plus1sigma" );
652 
653  saveS(selPNGCVN_All_PiAbsorptionProb_plus1sigma, "selPNGCVN_All_PiAbsorptionProb_plus1sigma" );
654  saveS(selSigPNGCVN_All_PiAbsorptionProb_plus1sigma, "selSigPNGCVN_All_PiAbsorptionProb_plus1sigma" );
655  saveS(selBkgPNGCVN_All_PiAbsorptionProb_plus1sigma, "selBkgPNGCVN_All_PiAbsorptionProb_plus1sigma" );
656 
657 
658  saveS(selPNGCVN_All_PiProductionProb_plus1sigma, "selPNGCVN_All_PiProductionProb_plus1sigma" );
659  saveS(selSigPNGCVN_All_PiProductionProb_plus1sigma, "selSigPNGCVN_All_PiProductionProb_plus1sigma" );
660  saveS(selBkgPNGCVN_All_PiProductionProb_plus1sigma, "selBkgPNGCVN_All_PiProductionProb_plus1sigma" );
661 
662 
663 
664  saveS(selPNGCVN_All_NucElasticProb_plus1sigma, "selPNGCVN_All_NucElasticProb_plus1sigma" );
665  saveS(selSigPNGCVN_All_NucElasticProb_plus1sigma, "selSigPNGCVN_All_NucElasticProb_plus1sigma" );
666  saveS(selBkgPNGCVN_All_NucElasticProb_plus1sigma, "selBkgPNGCVN_All_NucElasticProb_plus1sigma" );
667 
668  saveS(selPNGCVN_All_NucInElasticProb_plus1sigma, "selPNGCVN_All_NucInElasticProb_plus1sigma" );
669  saveS(selSigPNGCVN_All_NucInElasticProb_plus1sigma, "selSigPNGCVN_All_NucInElasticProb_plus1sigma" );
670  saveS(selBkgPNGCVN_All_NucInElasticProb_plus1sigma, "selBkgPNGCVN_All_NucInElasticProb_plus1sigma" );
671 
672 
673  saveS(selPNGCVN_All_NucChargeExchangeProb_plus1sigma, "selPNGCVN_All_NucChargeExchangeProb_plus1sigma" );
674  saveS(selSigPNGCVN_All_NucChargeExchangeProb_plus1sigma, "selSigPNGCVN_All_NucChargeExchangeProb_plus1sigma" );
675  saveS(selBkgPNGCVN_All_NucChargeExchangeProb_plus1sigma, "selBkgPNGCVN_All_NucChargeExchangeProb_plus1sigma" );
676 
677  saveS(selPNGCVN_All_NucAbsorptionProb_plus1sigma, "selPNGCVN_All_NucAbsorptionProb_plus1sigma" );
678  saveS(selSigPNGCVN_All_NucAbsorptionProb_plus1sigma, "selSigPNGCVN_All_NucAbsorptionProb_plus1sigma" );
679  saveS(selBkgPNGCVN_All_NucAbsorptionProb_plus1sigma, "selBkgPNGCVN_All_NucAbsorptionProb_plus1sigma" );
680 
681 
682  saveS(selPNGCVN_All_NucProductionProb_plus1sigma, "selPNGCVN_All_NucProductionProb_plus1sigma" );
683  saveS(selSigPNGCVN_All_NucProductionProb_plus1sigma, "selSigPNGCVN_All_NucProductionProb_plus1sigma" );
684  saveS(selBkgPNGCVN_All_NucProductionProb_plus1sigma, "selBkgPNGCVN_All_NucProductionProb_plus1sigma" );
685 
686  saveS(selPNGCVN_FermiMom_plus1sigma, "selPNGCVN_FermiMom_plus1sigma" );
687  saveS(selSigPNGCVN_FermiMom_plus1sigma, "selSigPNGCVN_FermiMom_plus1sigma" );
688  saveS(selBkgPNGCVN_FermiMom_plus1sigma, "selBkgPNGCVN_FermiMom_plus1sigma" );
689 
690 
691  saveS(selPNGCVN_SpectralFn_plus1sigma, "selPNGCVN_SpectralFn_plus1sigma" );
692  saveS(selSigPNGCVN_SpectralFn_plus1sigma, "selSigPNGCVN_SpectralFn_plus1sigma" );
693  saveS(selBkgPNGCVN_SpectralFn_plus1sigma, "selBkgPNGCVN_SpectralFn_plus1sigma" );
694 
695 
696 
697  saveS(selPNGCVN_All_GammaBR_plus1sigma, "selPNGCVN_All_GammaBR_plus1sigma" );
698  saveS(selSigPNGCVN_All_GammaBR_plus1sigma, "selSigPNGCVN_All_GammaBR_plus1sigma" );
699  saveS(selBkgPNGCVN_All_GammaBR_plus1sigma, "selBkgPNGCVN_All_GammaBR_plus1sigma" );
700 
701  saveS(selPNGCVN_All_EtaBR_plus1sigma, "selPNGCVN_All_EtaBR_plus1sigma" );
702  saveS(selSigPNGCVN_All_EtaBR_plus1sigma, "selSigPNGCVN_All_EtaBR_plus1sigma" );
703  saveS(selBkgPNGCVN_All_EtaBR_plus1sigma, "selBkgPNGCVN_All_EtaBR_plus1sigma" );
704 
705  saveS(selPNGCVN_All_PiAngularDist_plus1sigma, "selPNGCVN_All_PiAngularDist_plus1sigma" );
706  saveS(selSigPNGCVN_All_PiAngularDist_plus1sigma, "selSigPNGCVN_All_PiAngularDist_plus1sigma" );
707  saveS(selBkgPNGCVN_All_PiAngularDist_plus1sigma, "selBkgPNGCVN_All_PiAngularDist_plus1sigma" );
708 
709 
710  saveS(selPNGCVN_xF_plus1sigma, "selPNGCVN_xF_plus1sigma" );
711  saveS(selSigPNGCVN_xF_plus1sigma, "selSigPNGCVN_xF_plus1sigma" );
712  saveS(selBkgPNGCVN_xF_plus1sigma, "selBkgPNGCVN_xF_plus1sigma" );
713 
714  saveS(selPNGCVN_pT_plus1sigma, "selPNGCVN_pT_plus1sigma" );
715  saveS(selSigPNGCVN_pT_plus1sigma, "selSigPNGCVN_pT_plus1sigma" );
716  saveS(selBkgPNGCVN_pT_plus1sigma, "selBkgPNGCVN_pT_plus1sigma" );
717 
718  saveS(selPNGCVN_FormZone_plus1sigma, "selPNGCVN_FormZone_plus1sigma" );
719  saveS(selSigPNGCVN_FormZone_plus1sigma, "selSigPNGCVN_FormZone_plus1sigma" );
720  saveS(selBkgPNGCVN_FormZone_plus1sigma, "selBkgPNGCVN_FormZone_plus1sigma" );
721 
722 
723  saveS(selPNGCVN_All_1piNonRES_nup_plus1sigma, "selPNGCVN_All_1piNonRES_nup_plus1sigma" );
724  saveS(selSigPNGCVN_All_1piNonRES_nup_plus1sigma, "selSigPNGCVN_All_1piNonRES_nup_plus1sigma" );
725  saveS(selBkgPNGCVN_All_1piNonRES_nup_plus1sigma, "selBkgPNGCVN_All_1piNonRES_nup_plus1sigma" );
726 
727  saveS(selPNGCVN_All_2piNonRES_nup_plus1sigma, "selPNGCVN_All_2piNonRES_nup_plus1sigma" );
728  saveS(selSigPNGCVN_All_2piNonRES_nup_plus1sigma, "selSigPNGCVN_All_2piNonRES_nup_plus1sigma" );
729  saveS(selBkgPNGCVN_All_2piNonRES_nup_plus1sigma, "selBkgPNGCVN_All_2piNonRES_nup_plus1sigma" );
730 
731 saveS(selPNGCVN_All_1piNonRES_nun_plus1sigma, "selPNGCVN_All_1piNonRES_nun_plus1sigma" );
732  saveS(selSigPNGCVN_All_1piNonRES_nun_plus1sigma, "selSigPNGCVN_All_1piNonRES_nun_plus1sigma" );
733  saveS(selBkgPNGCVN_All_1piNonRES_nun_plus1sigma, "selBkgPNGCVN_All_1piNonRES_nun_plus1sigma" );
734 
735  saveS(selPNGCVN_All_2piNonRES_nun_plus1sigma, "selPNGCVN_All_2piNonRES_nun_plus1sigma" );
736  saveS(selSigPNGCVN_All_2piNonRES_nun_plus1sigma, "selSigPNGCVN_All_2piNonRES_nun_plus1sigma" );
737  saveS(selBkgPNGCVN_All_2piNonRES_nun_plus1sigma, "selBkgPNGCVN_All_2piNonRES_nun_plus1sigma" );
738 
739  saveS(selPNGCVN_All_1piNonRES_antinup_plus1sigma, "selPNGCVN_All_1piNonRES_antinup_plus1sigma" );
740  saveS(selSigPNGCVN_All_1piNonRES_antinup_plus1sigma, "selSigPNGCVN_All_1piNonRES_antinup_plus1sigma" );
741  saveS(selBkgPNGCVN_All_1piNonRES_antinup_plus1sigma, "selBkgPNGCVN_All_1piNonRES_antinup_plus1sigma" );
742 
743  saveS(selPNGCVN_All_2piNonRES_antinup_plus1sigma, "selPNGCVN_All_2piNonRES_antinup_plus1sigma" );
744  saveS(selSigPNGCVN_All_2piNonRES_antinup_plus1sigma, "selSigPNGCVN_All_2piNonRES_antinup_plus1sigma" );
745  saveS(selBkgPNGCVN_All_2piNonRES_antinup_plus1sigma, "selBkgPNGCVN_All_2piNonRES_antinup_plus1sigma" );
746 
747  saveS(selPNGCVN_All_1piNonRES_antinun_plus1sigma, "selPNGCVN_All_1piNonRES_antinun_plus1sigma" );
748  saveS(selSigPNGCVN_All_1piNonRES_antinun_plus1sigma, "selSigPNGCVN_All_1piNonRES_antinun_plus1sigma" );
749  saveS(selBkgPNGCVN_All_1piNonRES_antinun_plus1sigma, "selBkgPNGCVN_All_1piNonRES_antinun_plus1sigma" );
750 
751  saveS(selPNGCVN_All_2piNonRES_antinun_plus1sigma, "selPNGCVN_All_2piNonRES_antinun_plus1sigma" );
752  saveS(selSigPNGCVN_All_2piNonRES_antinun_plus1sigma, "selSigPNGCVN_All_2piNonRES_antinun_plus1sigma" );
753  saveS(selBkgPNGCVN_All_2piNonRES_antinun_plus1sigma, "selBkgPNGCVN_All_2piNonRES_antinun_plus1sigma" );
754 
755 
756 
757  //CC
758 
759  saveS(selPNGCVN_All_CC_1piNonRES_nup_plus1sigma, "selPNGCVN_All_CC_1piNonRES_nup_plus1sigma" );
760  saveS(selSigPNGCVN_All_CC_1piNonRES_nup_plus1sigma, "selSigPNGCVN_All_CC_1piNonRES_nup_plus1sigma" );
761  saveS(selBkgPNGCVN_All_CC_1piNonRES_nup_plus1sigma, "selBkgPNGCVN_All_CC_1piNonRES_nup_plus1sigma" );
762 
763  saveS(selPNGCVN_All_CC_2piNonRES_nup_plus1sigma, "selPNGCVN_All_CC_2piNonRES_nup_plus1sigma" );
764  saveS(selSigPNGCVN_All_CC_2piNonRES_nup_plus1sigma, "selSigPNGCVN_All_CC_2piNonRES_nup_plus1sigma" );
765  saveS(selBkgPNGCVN_All_CC_2piNonRES_nup_plus1sigma, "selBkgPNGCVN_All_CC_2piNonRES_nup_plus1sigma" );
766 
767  saveS(selPNGCVN_All_CC_1piNonRES_nun_plus1sigma, "selPNGCVN_All_CC_1piNonRES_nun_plus1sigma" );
768  saveS(selSigPNGCVN_All_CC_1piNonRES_nun_plus1sigma, "selSigPNGCVN_All_CC_1piNonRES_nun_plus1sigma" );
769  saveS(selBkgPNGCVN_All_CC_1piNonRES_nun_plus1sigma, "selBkgPNGCVN_All_CC_1piNonRES_nun_plus1sigma" );
770 
771  saveS(selPNGCVN_All_CC_2piNonRES_nun_plus1sigma, "selPNGCVN_All_CC_2piNonRES_nun_plus1sigma" );
772  saveS(selSigPNGCVN_All_CC_2piNonRES_nun_plus1sigma, "selSigPNGCVN_All_CC_2piNonRES_nun_plus1sigma" );
773  saveS(selBkgPNGCVN_All_CC_2piNonRES_nun_plus1sigma, "selBkgPNGCVN_All_CC_2piNonRES_nun_plus1sigma" );
774 
775  saveS(selPNGCVN_All_CC_1piNonRES_antinup_plus1sigma, "selPNGCVN_All_CC_1piNonRES_antinup_plus1sigma" );
776  saveS(selSigPNGCVN_All_CC_1piNonRES_antinup_plus1sigma, "selSigPNGCVN_All_CC_1piNonRES_antinup_plus1sigma" );
777  saveS(selBkgPNGCVN_All_CC_1piNonRES_antinup_plus1sigma, "selBkgPNGCVN_All_CC_1piNonRES_antinup_plus1sigma" );
778 
779  saveS(selPNGCVN_All_CC_2piNonRES_antinup_plus1sigma, "selPNGCVN_All_CC_2piNonRES_antinup_plus1sigma" );
780  saveS(selSigPNGCVN_All_CC_2piNonRES_antinup_plus1sigma, "selSigPNGCVN_All_CC_2piNonRES_antinup_plus1sigma" );
781  saveS(selBkgPNGCVN_All_CC_2piNonRES_antinup_plus1sigma, "selBkgPNGCVN_All_CC_2piNonRES_antinup_plus1sigma" );
782 
783 saveS(selPNGCVN_All_CC_1piNonRES_antinun_plus1sigma, "selPNGCVN_All_CC_1piNonRES_antinun_plus1sigma" );
784  saveS(selSigPNGCVN_All_CC_1piNonRES_antinun_plus1sigma, "selSigPNGCVN_All_CC_1piNonRES_antinun_plus1sigma" );
785  saveS(selBkgPNGCVN_All_CC_1piNonRES_antinun_plus1sigma, "selBkgPNGCVN_All_CC_1piNonRES_antinun_plus1sigma" );
786 
787  saveS(selPNGCVN_All_CC_2piNonRES_antinun_plus1sigma, "selPNGCVN_All_CC_2piNonRES_antinun_plus1sigma" );
788  saveS(selSigPNGCVN_All_CC_2piNonRES_antinun_plus1sigma, "selSigPNGCVN_All_CC_2piNonRES_antinun_plus1sigma" );
789  saveS(selBkgPNGCVN_All_CC_2piNonRES_antinun_plus1sigma, "selBkgPNGCVN_All_CC_2piNonRES_antinun_plus1sigma" );
790 
791 
792  saveS(selPNGCVN_ShapeNorm_AxialMassNCEL_minus1sigma, "selPNGCVN_ShapeNorm_AxialMassNCEL_minus1sigma" );
793  saveS(selSigPNGCVN_ShapeNorm_AxialMassNCEL_minus1sigma, "selSigPNGCVN_ShapeNorm_AxialMassNCEL_minus1sigma" );
794  saveS(selBkgPNGCVN_ShapeNorm_AxialMassNCEL_minus1sigma, "selBkgPNGCVN_ShapeNorm_AxialMassNCEL_minus1sigma" );
795 saveS(selPNGCVN_ShapeNorm_EtaNCEL_minus1sigma, "selPNGCVN_ShapeNorm_EtaNCEL_minus1sigma" );
796  saveS(selSigPNGCVN_ShapeNorm_EtaNCEL_minus1sigma, "selSigPNGCVN_ShapeNorm_EtaNCEL_minus1sigma" );
797  saveS(selBkgPNGCVN_ShapeNorm_EtaNCEL_minus1sigma, "selBkgPNGCVN_ShapeNorm_EtaNCEL_minus1sigma" );
798 
799  saveS(selPNGCVN_ShapeNorm_AxialMassNCRES_minus1sigma, "selPNGCVN_ShapeNorm_AxialMassNCRES_minus1sigma" );
800  saveS(selSigPNGCVN_ShapeNorm_AxialMassNCRES_minus1sigma, "selSigPNGCVN_ShapeNorm_AxialMassNCRES_minus1sigma" );
801  saveS(selBkgPNGCVN_ShapeNorm_AxialMassNCRES_minus1sigma, "selBkgPNGCVN_ShapeNorm_AxialMassNCRES_minus1sigma" );
802 
803  saveS(selPNGCVN_ShapeNorm_VectorMassNCRES_minus1sigma, "selPNGCVN_ShapeNorm_VectorMassNCRES_minus1sigma" );
804  saveS(selSigPNGCVN_ShapeNorm_VectorMassNCRES_minus1sigma, "selSigPNGCVN_ShapeNorm_VectorMassNCRES_minus1sigma" );
805  saveS(selBkgPNGCVN_ShapeNorm_VectorMassNCRES_minus1sigma, "selBkgPNGCVN_ShapeNorm_VectorMassNCRES_minus1sigma" );
806 
807 
808 saveS(selPNGCVN_ShapeNorm_AxialMassCCRES_minus1sigma, "selPNGCVN_ShapeNorm_AxialMassCCRES_minus1sigma" );
809  saveS(selSigPNGCVN_ShapeNorm_AxialMassCCRES_minus1sigma, "selSigPNGCVN_ShapeNorm_AxialMassCCRES_minus1sigma" );
810  saveS(selBkgPNGCVN_ShapeNorm_AxialMassCCRES_minus1sigma, "selBkgPNGCVN_ShapeNorm_AxialMassCCRES_minus1sigma" );
811 
812  saveS(selPNGCVN_ShapeNorm_VectorMassCCRES_minus1sigma, "selPNGCVN_ShapeNorm_VectorMassCCRES_minus1sigma" );
813  saveS(selSigPNGCVN_ShapeNorm_VectorMassCCRES_minus1sigma, "selSigPNGCVN_ShapeNorm_VectorMassCCRES_minus1sigma" );
814  saveS(selBkgPNGCVN_ShapeNorm_VectorMassCCRES_minus1sigma, "selBkgPNGCVN_ShapeNorm_VectorMassCCRES_minus1sigma" );
815 
816 
817 saveS(selPNGCVN_All_AxialMassCOH_minus1sigma, "selPNGCVN_All_AxialMassCOH_minus1sigma" );
818  saveS(selSigPNGCVN_All_AxialMassCOH_minus1sigma, "selSigPNGCVN_All_AxialMassCOH_minus1sigma" );
819  saveS(selBkgPNGCVN_All_AxialMassCOH_minus1sigma, "selBkgPNGCVN_All_AxialMassCOH_minus1sigma" );
820 
821 saveS(selPNGCVN_All_NuclearSizeCOH_minus1sigma, "selPNGCVN_All_NuclearSizeCOH_minus1sigma" );
822  saveS(selSigPNGCVN_All_NuclearSizeCOH_minus1sigma, "selSigPNGCVN_All_NuclearSizeCOH_minus1sigma" );
823  saveS(selBkgPNGCVN_All_NuclearSizeCOH_minus1sigma, "selBkgPNGCVN_All_NuclearSizeCOH_minus1sigma" );
824 
825 saveS(selPNGCVN_ShapeNorm_DIS_Aht_minus1sigma, "selPNGCVN_ShapeNorm_DIS_Aht_minus1sigma" );
826  saveS(selSigPNGCVN_ShapeNorm_DIS_Aht_minus1sigma, "selSigPNGCVN_ShapeNorm_DIS_Aht_minus1sigma" );
827  saveS(selBkgPNGCVN_ShapeNorm_DIS_Aht_minus1sigma, "selBkgPNGCVN_ShapeNorm_DIS_Aht_minus1sigma" );
828 
829 saveS(selPNGCVN_ShapeNorm_DIS_Bht_minus1sigma, "selPNGCVN_ShapeNorm_DIS_Bht_minus1sigma" );
830  saveS(selSigPNGCVN_ShapeNorm_DIS_Bht_minus1sigma, "selSigPNGCVN_ShapeNorm_DIS_Bht_minus1sigma" );
831  saveS(selBkgPNGCVN_ShapeNorm_DIS_Bht_minus1sigma, "selBkgPNGCVN_ShapeNorm_DIS_Bht_minus1sigma" );
832 
833 saveS(selPNGCVN_ShapeNorm_DIS_CV1u_minus1sigma, "selPNGCVN_ShapeNorm_DIS_CV1u_minus1sigma" );
834  saveS(selSigPNGCVN_ShapeNorm_DIS_CV1u_minus1sigma, "selSigPNGCVN_ShapeNorm_DIS_CV1u_minus1sigma" );
835  saveS(selBkgPNGCVN_ShapeNorm_DIS_CV1u_minus1sigma, "selBkgPNGCVN_ShapeNorm_DIS_CV1u_minus1sigma" );
836 
837 saveS(selPNGCVN_ShapeNorm_DIS_CV2u_minus1sigma, "selPNGCVN_ShapeNorm_DIS_CV2u_minus1sigma" );
838  saveS(selSigPNGCVN_ShapeNorm_DIS_CV2u_minus1sigma, "selSigPNGCVN_ShapeNorm_DIS_CV2u_minus1sigma" );
839  saveS(selBkgPNGCVN_ShapeNorm_DIS_CV2u_minus1sigma, "selBkgPNGCVN_ShapeNorm_DIS_CV2u_minus1sigma" );
840 
841 
842 saveS(selPNGCVN_mfp_pi_minus1sigma, "selPNGCVN_mfp_pi_minus1sigma" );
843  saveS(selSigPNGCVN_mfp_pi_minus1sigma, "selSigPNGCVN_mfp_pi_minus1sigma" );
844  saveS(selBkgPNGCVN_mfp_pi_minus1sigma, "selBkgPNGCVN_mfp_pi_minus1sigma" );
845 
846 saveS(selPNGCVN_mfp_N_minus1sigma, "selPNGCVN_mfp_N_minus1sigma" );
847  saveS(selSigPNGCVN_mfp_N_minus1sigma, "selSigPNGCVN_mfp_N_minus1sigma" );
848  saveS(selBkgPNGCVN_mfp_N_minus1sigma, "selBkgPNGCVN_mfp_N_minus1sigma" );
849 
850 
851 saveS(selPNGCVN_All_PiElasticProb_minus1sigma, "selPNGCVN_All_PiElasticProb_minus1sigma" );
852  saveS(selSigPNGCVN_All_PiElasticProb_minus1sigma, "selSigPNGCVN_All_PiElasticProb_minus1sigma" );
853  saveS(selBkgPNGCVN_All_PiElasticProb_minus1sigma, "selBkgPNGCVN_All_PiElasticProb_minus1sigma" );
854 
855 saveS(selPNGCVN_All_PiInElasticProb_minus1sigma, "selPNGCVN_All_PiInElasticProb_minus1sigma" );
856  saveS(selSigPNGCVN_All_PiInElasticProb_minus1sigma, "selSigPNGCVN_All_PiInElasticProb_minus1sigma" );
857  saveS(selBkgPNGCVN_All_PiInElasticProb_minus1sigma, "selBkgPNGCVN_All_PiInElasticProb_minus1sigma" );
858 
859 
860 saveS(selPNGCVN_All_PiChargeExchangeProb_minus1sigma, "selPNGCVN_All_PiChargeExchangeProb_minus1sigma" );
861  saveS(selSigPNGCVN_All_PiChargeExchangeProb_minus1sigma, "selSigPNGCVN_All_PiChargeExchangeProb_minus1sigma" );
862  saveS(selBkgPNGCVN_All_PiChargeExchangeProb_minus1sigma, "selBkgPNGCVN_All_PiChargeExchangeProb_minus1sigma" );
863 
864 saveS(selPNGCVN_All_PiAbsorptionProb_minus1sigma, "selPNGCVN_All_PiAbsorptionProb_minus1sigma" );
865  saveS(selSigPNGCVN_All_PiAbsorptionProb_minus1sigma, "selSigPNGCVN_All_PiAbsorptionProb_minus1sigma" );
866  saveS(selBkgPNGCVN_All_PiAbsorptionProb_minus1sigma, "selBkgPNGCVN_All_PiAbsorptionProb_minus1sigma" );
867 
868 
869 saveS(selPNGCVN_All_PiProductionProb_minus1sigma, "selPNGCVN_All_PiProductionProb_minus1sigma" );
870  saveS(selSigPNGCVN_All_PiProductionProb_minus1sigma, "selSigPNGCVN_All_PiProductionProb_minus1sigma" );
871  saveS(selBkgPNGCVN_All_PiProductionProb_minus1sigma, "selBkgPNGCVN_All_PiProductionProb_minus1sigma" );
872 
873 
874 
875 saveS(selPNGCVN_All_NucElasticProb_minus1sigma, "selPNGCVN_All_NucElasticProb_minus1sigma" );
876  saveS(selSigPNGCVN_All_NucElasticProb_minus1sigma, "selSigPNGCVN_All_NucElasticProb_minus1sigma" );
877  saveS(selBkgPNGCVN_All_NucElasticProb_minus1sigma, "selBkgPNGCVN_All_NucElasticProb_minus1sigma" );
878 
879 saveS(selPNGCVN_All_NucInElasticProb_minus1sigma, "selPNGCVN_All_NucInElasticProb_minus1sigma" );
880  saveS(selSigPNGCVN_All_NucInElasticProb_minus1sigma, "selSigPNGCVN_All_NucInElasticProb_minus1sigma" );
881  saveS(selBkgPNGCVN_All_NucInElasticProb_minus1sigma, "selBkgPNGCVN_All_NucInElasticProb_minus1sigma" );
882 
883 
884 saveS(selPNGCVN_All_NucChargeExchangeProb_minus1sigma, "selPNGCVN_All_NucChargeExchangeProb_minus1sigma" );
885  saveS(selSigPNGCVN_All_NucChargeExchangeProb_minus1sigma, "selSigPNGCVN_All_NucChargeExchangeProb_minus1sigma" );
886  saveS(selBkgPNGCVN_All_NucChargeExchangeProb_minus1sigma, "selBkgPNGCVN_All_NucChargeExchangeProb_minus1sigma" );
887 
888 saveS(selPNGCVN_All_NucAbsorptionProb_minus1sigma, "selPNGCVN_All_NucAbsorptionProb_minus1sigma" );
889  saveS(selSigPNGCVN_All_NucAbsorptionProb_minus1sigma, "selSigPNGCVN_All_NucAbsorptionProb_minus1sigma" );
890  saveS(selBkgPNGCVN_All_NucAbsorptionProb_minus1sigma, "selBkgPNGCVN_All_NucAbsorptionProb_minus1sigma" );
891 
892 
893 saveS(selPNGCVN_All_NucProductionProb_minus1sigma, "selPNGCVN_All_NucProductionProb_minus1sigma" );
894  saveS(selSigPNGCVN_All_NucProductionProb_minus1sigma, "selSigPNGCVN_All_NucProductionProb_minus1sigma" );
895  saveS(selBkgPNGCVN_All_NucProductionProb_minus1sigma, "selBkgPNGCVN_All_NucProductionProb_minus1sigma" );
896 
897  saveS(selPNGCVN_FermiMom_minus1sigma, "selPNGCVN_FermiMom_minus1sigma" );
898  saveS(selSigPNGCVN_FermiMom_minus1sigma, "selSigPNGCVN_FermiMom_minus1sigma" );
899  saveS(selBkgPNGCVN_FermiMom_minus1sigma, "selBkgPNGCVN_FermiMom_minus1sigma" );
900 
901 
902  saveS(selPNGCVN_SpectralFn_minus1sigma, "selPNGCVN_SpectralFn_minus1sigma" );
903  saveS(selSigPNGCVN_SpectralFn_minus1sigma, "selSigPNGCVN_SpectralFn_minus1sigma" );
904  saveS(selBkgPNGCVN_SpectralFn_minus1sigma, "selBkgPNGCVN_SpectralFn_minus1sigma" );
905 
906 
907 
908 saveS(selPNGCVN_All_GammaBR_minus1sigma, "selPNGCVN_All_GammaBR_minus1sigma" );
909  saveS(selSigPNGCVN_All_GammaBR_minus1sigma, "selSigPNGCVN_All_GammaBR_minus1sigma" );
910  saveS(selBkgPNGCVN_All_GammaBR_minus1sigma, "selBkgPNGCVN_All_GammaBR_minus1sigma" );
911 
912 saveS(selPNGCVN_All_EtaBR_minus1sigma, "selPNGCVN_All_EtaBR_minus1sigma" );
913  saveS(selSigPNGCVN_All_EtaBR_minus1sigma, "selSigPNGCVN_All_EtaBR_minus1sigma" );
914  saveS(selBkgPNGCVN_All_EtaBR_minus1sigma, "selBkgPNGCVN_All_EtaBR_minus1sigma" );
915 
916 saveS(selPNGCVN_All_PiAngularDist_minus1sigma, "selPNGCVN_All_PiAngularDist_minus1sigma" );
917  saveS(selSigPNGCVN_All_PiAngularDist_minus1sigma, "selSigPNGCVN_All_PiAngularDist_minus1sigma" );
918  saveS(selBkgPNGCVN_All_PiAngularDist_minus1sigma, "selBkgPNGCVN_All_PiAngularDist_minus1sigma" );
919 
920 
921  saveS(selPNGCVN_xF_minus1sigma, "selPNGCVN_xF_minus1sigma" );
922  saveS(selSigPNGCVN_xF_minus1sigma, "selSigPNGCVN_xF_minus1sigma" );
923  saveS(selBkgPNGCVN_xF_minus1sigma, "selBkgPNGCVN_xF_minus1sigma" );
924 
925  saveS(selPNGCVN_pT_minus1sigma, "selPNGCVN_pT_minus1sigma" );
926  saveS(selSigPNGCVN_pT_minus1sigma, "selSigPNGCVN_pT_minus1sigma" );
927  saveS(selBkgPNGCVN_pT_minus1sigma, "selBkgPNGCVN_pT_minus1sigma" );
928 
929  saveS(selPNGCVN_FormZone_minus1sigma, "selPNGCVN_FormZone_minus1sigma" );
930  saveS(selSigPNGCVN_FormZone_minus1sigma, "selSigPNGCVN_FormZone_minus1sigma" );
931  saveS(selBkgPNGCVN_FormZone_minus1sigma, "selBkgPNGCVN_FormZone_minus1sigma" );
932 
933 
934 saveS(selPNGCVN_All_1piNonRES_nup_minus1sigma, "selPNGCVN_All_1piNonRES_nup_minus1sigma" );
935  saveS(selSigPNGCVN_All_1piNonRES_nup_minus1sigma, "selSigPNGCVN_All_1piNonRES_nup_minus1sigma" );
936  saveS(selBkgPNGCVN_All_1piNonRES_nup_minus1sigma, "selBkgPNGCVN_All_1piNonRES_nup_minus1sigma" );
937 
938  saveS(selPNGCVN_All_2piNonRES_nup_minus1sigma, "selPNGCVN_All_2piNonRES_nup_minus1sigma" );
939  saveS(selSigPNGCVN_All_2piNonRES_nup_minus1sigma, "selSigPNGCVN_All_2piNonRES_nup_minus1sigma" );
940  saveS(selBkgPNGCVN_All_2piNonRES_nup_minus1sigma, "selBkgPNGCVN_All_2piNonRES_nup_minus1sigma" );
941 
942 saveS(selPNGCVN_All_1piNonRES_nun_minus1sigma, "selPNGCVN_All_1piNonRES_nun_minus1sigma" );
943  saveS(selSigPNGCVN_All_1piNonRES_nun_minus1sigma, "selSigPNGCVN_All_1piNonRES_nun_minus1sigma" );
944  saveS(selBkgPNGCVN_All_1piNonRES_nun_minus1sigma, "selBkgPNGCVN_All_1piNonRES_nun_minus1sigma" );
945 
946  saveS(selPNGCVN_All_2piNonRES_nun_minus1sigma, "selPNGCVN_All_2piNonRES_nun_minus1sigma" );
947  saveS(selSigPNGCVN_All_2piNonRES_nun_minus1sigma, "selSigPNGCVN_All_2piNonRES_nun_minus1sigma" );
948  saveS(selBkgPNGCVN_All_2piNonRES_nun_minus1sigma, "selBkgPNGCVN_All_2piNonRES_nun_minus1sigma" );
949 
950  saveS(selPNGCVN_All_1piNonRES_antinup_minus1sigma, "selPNGCVN_All_1piNonRES_antinup_minus1sigma" );
951  saveS(selSigPNGCVN_All_1piNonRES_antinup_minus1sigma, "selSigPNGCVN_All_1piNonRES_antinup_minus1sigma" );
952  saveS(selBkgPNGCVN_All_1piNonRES_antinup_minus1sigma, "selBkgPNGCVN_All_1piNonRES_antinup_minus1sigma" );
953 
954  saveS(selPNGCVN_All_2piNonRES_antinup_minus1sigma, "selPNGCVN_All_2piNonRES_antinup_minus1sigma" );
955  saveS(selSigPNGCVN_All_2piNonRES_antinup_minus1sigma, "selSigPNGCVN_All_2piNonRES_antinup_minus1sigma" );
956  saveS(selBkgPNGCVN_All_2piNonRES_antinup_minus1sigma, "selBkgPNGCVN_All_2piNonRES_antinup_minus1sigma" );
957 
958 saveS(selPNGCVN_All_1piNonRES_antinun_minus1sigma, "selPNGCVN_All_1piNonRES_antinun_minus1sigma" );
959  saveS(selSigPNGCVN_All_1piNonRES_antinun_minus1sigma, "selSigPNGCVN_All_1piNonRES_antinun_minus1sigma" );
960  saveS(selBkgPNGCVN_All_1piNonRES_antinun_minus1sigma, "selBkgPNGCVN_All_1piNonRES_antinun_minus1sigma" );
961 
962  saveS(selPNGCVN_All_2piNonRES_antinun_minus1sigma, "selPNGCVN_All_2piNonRES_antinun_minus1sigma" );
963  saveS(selSigPNGCVN_All_2piNonRES_antinun_minus1sigma, "selSigPNGCVN_All_2piNonRES_antinun_minus1sigma" );
964  saveS(selBkgPNGCVN_All_2piNonRES_antinun_minus1sigma, "selBkgPNGCVN_All_2piNonRES_antinun_minus1sigma" );
965 
966 
967  //CC
968 
969 saveS(selPNGCVN_All_CC_1piNonRES_nup_minus1sigma, "selPNGCVN_All_CC_1piNonRES_nup_minus1sigma" );
970  saveS(selSigPNGCVN_All_CC_1piNonRES_nup_minus1sigma, "selSigPNGCVN_All_CC_1piNonRES_nup_minus1sigma" );
971  saveS(selBkgPNGCVN_All_CC_1piNonRES_nup_minus1sigma, "selBkgPNGCVN_All_CC_1piNonRES_nup_minus1sigma" );
972 
973  saveS(selPNGCVN_All_CC_2piNonRES_nup_minus1sigma, "selPNGCVN_All_CC_2piNonRES_nup_minus1sigma" );
974  saveS(selSigPNGCVN_All_CC_2piNonRES_nup_minus1sigma, "selSigPNGCVN_All_CC_2piNonRES_nup_minus1sigma" );
975  saveS(selBkgPNGCVN_All_CC_2piNonRES_nup_minus1sigma, "selBkgPNGCVN_All_CC_2piNonRES_nup_minus1sigma" );
976 
977 saveS(selPNGCVN_All_CC_1piNonRES_nun_minus1sigma, "selPNGCVN_All_CC_1piNonRES_nun_minus1sigma" );
978  saveS(selSigPNGCVN_All_CC_1piNonRES_nun_minus1sigma, "selSigPNGCVN_All_CC_1piNonRES_nun_minus1sigma" );
979  saveS(selBkgPNGCVN_All_CC_1piNonRES_nun_minus1sigma, "selBkgPNGCVN_All_CC_1piNonRES_nun_minus1sigma" );
980 
981  saveS(selPNGCVN_All_CC_2piNonRES_nun_minus1sigma, "selPNGCVN_All_CC_2piNonRES_nun_minus1sigma" );
982  saveS(selSigPNGCVN_All_CC_2piNonRES_nun_minus1sigma, "selSigPNGCVN_All_CC_2piNonRES_nun_minus1sigma" );
983  saveS(selBkgPNGCVN_All_CC_2piNonRES_nun_minus1sigma, "selBkgPNGCVN_All_CC_2piNonRES_nun_minus1sigma" );
984 
985  saveS(selPNGCVN_All_CC_1piNonRES_antinup_minus1sigma, "selPNGCVN_All_CC_1piNonRES_antinup_minus1sigma" );
986  saveS(selSigPNGCVN_All_CC_1piNonRES_antinup_minus1sigma, "selSigPNGCVN_All_CC_1piNonRES_antinup_minus1sigma" );
987  saveS(selBkgPNGCVN_All_CC_1piNonRES_antinup_minus1sigma, "selBkgPNGCVN_All_CC_1piNonRES_antinup_minus1sigma" );
988 
989  saveS(selPNGCVN_All_CC_2piNonRES_antinup_minus1sigma, "selPNGCVN_All_CC_2piNonRES_antinup_minus1sigma" );
990  saveS(selSigPNGCVN_All_CC_2piNonRES_antinup_minus1sigma, "selSigPNGCVN_All_CC_2piNonRES_antinup_minus1sigma" );
991  saveS(selBkgPNGCVN_All_CC_2piNonRES_antinup_minus1sigma, "selBkgPNGCVN_All_CC_2piNonRES_antinup_minus1sigma" );
992 
993 saveS(selPNGCVN_All_CC_1piNonRES_antinun_minus1sigma, "selPNGCVN_All_CC_1piNonRES_antinun_minus1sigma" );
994  saveS(selSigPNGCVN_All_CC_1piNonRES_antinun_minus1sigma, "selSigPNGCVN_All_CC_1piNonRES_antinun_minus1sigma" );
995  saveS(selBkgPNGCVN_All_CC_1piNonRES_antinun_minus1sigma, "selBkgPNGCVN_All_CC_1piNonRES_antinun_minus1sigma" );
996 
997  saveS(selPNGCVN_All_CC_2piNonRES_antinun_minus1sigma, "selPNGCVN_All_CC_2piNonRES_antinun_minus1sigma" );
998  saveS(selSigPNGCVN_All_CC_2piNonRES_antinun_minus1sigma, "selSigPNGCVN_All_CC_2piNonRES_antinun_minus1sigma" );
999  saveS(selBkgPNGCVN_All_CC_2piNonRES_antinun_minus1sigma, "selBkgPNGCVN_All_CC_2piNonRES_antinun_minus1sigma" );
1000 
1001 
1002  fout.Write();
1003  fout.Close();
1004 
1005 
1006 
1007  std::cout<<"****Compiled Successfully*****"<<std::endl;
1008 
1009 }
1010 
1011 
1012 #endif
void saveS(Spectrum s, string name)
Definition: Syst.C:63
const XML_Char * name
Definition: expat.h:151
enum BeamMode kRed
TH2 * ToTH2(double exposure, EExposureType expotype=kPOT, EBinType bintype=kBinContent) const
Spectrum must be 2D to obtain TH2.
Definition: Spectrum.cxx:165
void saveS1(Spectrum s, string name)
Definition: Syst.C:70
Cuts and Vars for the 2020 FD DiF Study.
Definition: vars.h:6
const double pot
Definition: Syst.C:61
TH1D * ToTH1(double exposure, Color_t col=kBlack, Style_t style=kSolid, EExposureType expotype=kPOT, EBinType bintype=kBinContent) const
Histogram made from this Spectrum, scaled to some exposure.
Definition: Spectrum.cxx:148
Simple record of shifts applied to systematic parameters.
Definition: SystShifts.h:20
void SetSpillCut(const SpillCut &cut)
Representation of a spectrum in any variable, with associated POT.
Definition: Spectrum.h:40
const Var kNpCVN1gamma([](const caf::SRProxy *sr){if(sr->vtx.nelastic==0) return 0.f;if(sr->vtx.elastic[0].fuzzyk.npng< 1) return float(-1);return float(sr->vtx.elastic[0].fuzzyk.png[0].cvnpart.photonid);})
const SpillCut kTightBeamQualityCuts([](const caf::SRSpillProxy *s){if(s->ismc) return true; if(s->trigger==2) return true;if(s->spilltimesec==0 &&s->deltaspilltimensec==0 &&s->widthx==0) return bool(s->isgoodspill);if(std::abs(s->deltaspilltimensec) > 0.5e9) return false;if(s->spillpot< 2e12) return false;if(s->hornI< -202|| s->hornI >-198) return false;if(s->posx< -2.00|| s->posx >+2.00) return false;if(s->posy< -2.00|| s->posy >+2.00) return false;return kBeamWidthCut(s);})
Definition: SpillCuts.h:10
const Cut kIsTotBkg
Definition: ncpi0Cuts.h:1095
const XML_Char * s
Definition: expat.h:262
const SpillCut kStandardDQCuts([](const caf::SRSpillProxy *spill){if(spill->dcmedgematchfrac==0 &&spill->fracdcm3hits==0 &&spill->nmissingdcmslg==0) return bool(spill->isgoodspill); if(spill->det==caf::kNEARDET && (spill->fracdcm3hits > 0.45|| spill->nmissingdcms > 0)) return false; if(spill->eventincomplete) return false; if(spill->det==caf::kFARDET && spill->nmissingdcmslg > 0) return false; if(spill->det==caf::kFARDET && !spill->ismc && spill->dcmedgematchfrac<=0.2) return false;return true;})
Cut out events with a noisy detector or with parts missing.
Definition: SpillCuts.h:16
virtual void Go() override
Load all the registered spectra.
loader
Definition: demo0.py:10
const Binning PCIDbins
Definition: ncpi0Bins.h:24
const Cut kIsPreSel
Definition: ncpi0Cuts.h:1099
double POT() const
Definition: Spectrum.h:227
OStream cout
Definition: OStream.cxx:6
TH1F * h1
void Syst()
Definition: Syst.C:79
Collaborates with Spectrum and OscillatableSpectrum to fill spectra from CAF files.
const Cut kIsSignal
Definition: ncpi0Cuts.h:1094
enum BeamMode string