Functions | Variables
piN_asymmetric_decay.C File Reference
#include <TMath.h>
#include <TH1D.h>
#include <TH2D.h>
#include <TGraphAsymmErrors.h>
#include <TFile.h>
#include <TNtuple.h>
#include <TSpline.h>
#include <TRandom3.h>
#include <TGenPhaseSpace.h>
#include <TLorentzVector.h>

Go to the source code of this file.

Functions

TSpline3 splCoeffReY01 ("splCoeffReY01", Wc, CoeffReY01, kNWBins,"", 0, 0)
 
TSpline3 splCoeffReY02 ("splCoeffReY02", Wc, CoeffReY02, kNWBins,"", 1, 1)
 
TSpline3 splCoeffReY11 ("splCoeffReY11", Wc, CoeffReY11, kNWBins,"", 1, 1)
 
TSpline3 splCoeffImY11 ("splCoeffImY11", Wc, CoeffImY11, kNWBins,"", 1, 1)
 
TSpline3 splCoeffReY12 ("splCoeffReY12", Wc, CoeffReY12, kNWBins,"", 1, 1)
 
TSpline3 splCoeffImY12 ("splCoeffImY12", Wc, CoeffImY12, kNWBins,"", 1, 1)
 
TSpline3 splCoeffReY22 ("splCoeffReY22", Wc, CoeffReY22, kNWBins,"", 1, 1)
 
TSpline3 splCoeffImY22 ("splCoeffImY22", Wc, CoeffImY22, kNWBins,"", 1, 1)
 
double angular_distribution_pdf (double W, double th, double phi)
 
double coeff (double W, int m, int l, bool real)
 
double Y (double th, double phi, int m, int l, bool real)
 
double Y00 (double th, double phi, bool real)
 
double Y01 (double th, double phi, bool real)
 
double Y02 (double th, double phi, bool real)
 
double Y11 (double th, double phi, bool real)
 
double Y12 (double th, double phi, bool real)
 
double Y22 (double th, double phi, bool real)
 
void run_evgen_test (double W)
 
void plot_input_data (void)
 

Variables

const int kNWBins = 6
 
double Wmin [kNWBins] = { 1.100, 1.180, 1.230, 1.280, 1.500, 1.700 }
 
double Wmax [kNWBins] = { 1.180, 1.230, 1.280, 1.500, 1.700, 2.000 }
 
double Wc [kNWBins] = { 1.140, 1.205, 1.255, 1.390, 1.600, 1.850 }
 
double CoeffReY01 [kNWBins] = { -0.22, -0.01, 0.29, 0.32, 0.61, 0.84 }
 
double CoeffReY02 [kNWBins] = { -0.06, -0.13, -0.19, -0.08, 0.49, 0.80 }
 
double CoeffReY11 [kNWBins] = { 0.21, -0.02, -0.20, -0.05, 0.12, -0.03 }
 
double CoeffImY11 [kNWBins] = { -0.09, -0.06, -0.06, 0.01, -0.08, -0.17 }
 
double CoeffReY12 [kNWBins] = { 0.26, 0.05, 0.02, 0.00, 0.09, -0.05 }
 
double CoeffImY12 [kNWBins] = { 0.07, 0.00, 0.16, -0.01, -0.04, -0.03 }
 
double CoeffReY22 [kNWBins] = { -0.04, 0.08, 0.08, 0.04, 0.09, -0.08 }
 
double CoeffImY22 [kNWBins] = { -0.07, -0.07, -0.01, 0.06, -0.15, -0.03 }
 
double dCoeffReY01 [kNWBins] = { 0.11, 0.07, 0.09, 0.09, 0.13, 0.13 }
 
double dCoeffReY02 [kNWBins] = { 0.12, 0.07, 0.10, 0.09, 0.15, 0.14 }
 
double dCoeffReY11 [kNWBins] = { 0.09, 0.06, 0.08, 0.06, 0.08, 0.07 }
 
double dCoeffImY11 [kNWBins] = { 0.08, 0.06, 0.08, 0.06, 0.09, 0.09 }
 
double dCoeffReY12 [kNWBins] = { 0.08, 0.06, 0.08, 0.06, 0.07, 0.06 }
 
double dCoeffImY12 [kNWBins] = { 0.07, 0.05, 0.08, 0.06, 0.08, 0.07 }
 
double dCoeffReY22 [kNWBins] = { 0.07, 0.05, 0.07, 0.06, 0.09, 0.09 }
 
double dCoeffImY22 [kNWBins] = { 0.08, 0.05, 0.08, 0.06, 0.08, 0.07 }
 
const double kPi = 3.14
 

Function Documentation

double angular_distribution_pdf ( double  W,
double  th,
double  phi 
)

Definition at line 273 of file piN_asymmetric_decay.C.

References coeff(), MakeMiniprodValidationCuts::f, ana::Sqrt(), and Y().

Referenced by plot_input_data(), and run_evgen_test().

274 {
275 // return Sum{m,l} { (interpolated coeff) x (spherical harmonic)}
276 //
277  double fr = 1 +
278  coeff(W, 0,1,true ) * Y(th,phi, 0,1,true ) +
279  coeff(W, 0,2,true ) * Y(th,phi, 0,2,true ) +
280  coeff(W, 1,1,true ) * Y(th,phi, 1,1,true ) +
281  coeff(W, 1,2,true ) * Y(th,phi, 1,2,true ) +
282  coeff(W, 2,2,true ) * Y(th,phi, 2,2,true );
283 
284  double fi = coeff(W, 1,1,false) * Y(th,phi, 1,1,false) +
285  coeff(W, 1,2,false) * Y(th,phi, 1,2,false) +
286  coeff(W, 2,2,false) * Y(th,phi, 2,2,false);
287 
288  double f = TMath::Sqrt(fr*fr+fi*fi);
289  return f;
290 }
double coeff(double W, int m, int l, bool real)
double Y(double th, double phi, int m, int l, bool real)
Var Sqrt(const Var &v)
Use to take sqrt of a var.
Definition: Var.cxx:326
#define W(x)
double coeff ( double  W,
int  m,
int  l,
bool  real 
)

Definition at line 292 of file piN_asymmetric_decay.C.

References splCoeffImY11(), splCoeffImY12(), splCoeffImY22(), splCoeffReY01(), splCoeffReY02(), splCoeffReY11(), splCoeffReY12(), and splCoeffReY22().

Referenced by ana::PredictionSystJoint2018::AddNormSyst(), ana::PredictionSystJointDemo::AddNormSyst(), ana::PredictionSyst3Flavor2020::AddNormSyst(), ana::PredictionSystNue2017::AddNormSysts(), ana::PredictionSystNumu2017::AddNormSysts(), angular_distribution_pdf(), fit_thresh_corrs(), fit_thresh_corrs_fb(), and plot_input_data().

293 {
294  if (m==0 && l==0 && real) return 1;
295  else if (m==0 && l==1 && real) return splCoeffReY01.Eval(W);
296  else if (m==0 && l==2 && real) return splCoeffReY02.Eval(W);
297  else if (m==1 && l==1 && real) return splCoeffReY11.Eval(W);
298  else if (m==1 && l==1 && !real) return splCoeffImY11.Eval(W);
299  else if (m==1 && l==2 && real) return splCoeffReY12.Eval(W);
300  else if (m==1 && l==2 && !real) return splCoeffImY12.Eval(W);
301  else if (m==2 && l==2 && real) return splCoeffReY22.Eval(W);
302  else if (m==2 && l==2 && !real) return splCoeffImY22.Eval(W);
303  else
304  return 0;
305 }
TSpline3 splCoeffImY22("splCoeffImY22", Wc, CoeffImY22, kNWBins,"", 1, 1)
TSpline3 splCoeffImY12("splCoeffImY12", Wc, CoeffImY12, kNWBins,"", 1, 1)
TSpline3 splCoeffReY11("splCoeffReY11", Wc, CoeffReY11, kNWBins,"", 1, 1)
TSpline3 splCoeffReY12("splCoeffReY12", Wc, CoeffReY12, kNWBins,"", 1, 1)
TSpline3 splCoeffReY01("splCoeffReY01", Wc, CoeffReY01, kNWBins,"", 0, 0)
TSpline3 splCoeffImY11("splCoeffImY11", Wc, CoeffImY11, kNWBins,"", 1, 1)
TSpline3 splCoeffReY02("splCoeffReY02", Wc, CoeffReY02, kNWBins,"", 1, 1)
#define W(x)
TSpline3 splCoeffReY22("splCoeffReY22", Wc, CoeffReY22, kNWBins,"", 1, 1)
void plot_input_data ( void  )

Definition at line 171 of file piN_asymmetric_decay.C.

References angular_distribution_pdf(), coeff(), CoeffImY11, CoeffImY12, CoeffImY22, CoeffReY01, CoeffReY02, CoeffReY11, CoeffReY12, CoeffReY22, dCoeffImY11, dCoeffImY12, dCoeffImY22, dCoeffReY01, dCoeffReY02, dCoeffReY11, dCoeffReY12, dCoeffReY22, MakeMiniprodValidationCuts::f, MECModelEnuComparisons::i, calib::j, kNWBins, kPi, nt, sum, w, W, Wc, Wmax, and Y().

172 {
173 // save all distributions that are input to the new 2-body decay scheme
174 //
175  TFile f("./inputs.root","recreate");
176  TNtuple nt("nt","","th:phi:W:sum:r01:r02:r11:i11:r12:i12:r22:i22");
177 
178  const int nw = 11;
179  const int nth = 180;
180  const int nphi = 360;
181 
182  double th_min = 0;
183  double th_max = kPi;
184  double phi_min = 0;
185  double phi_max = 2*kPi;
186  double dth = (th_max - th_min ) / (nth -1);
187  double dphi = (phi_max - phi_min) / (nphi-1);
188 
189  double W[nw] = { 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0 };
190 
191  for(int i=0; i<nth; i++) {
192  double th = th_min + i*dth;
193  for(int j=0; j<nphi; j++) {
194  double phi = phi_min + j*dphi;
195  for(int k=0; k<nw; k++) {
196  double w = W[k];
197 
198  double sum = angular_distribution_pdf(w, th, phi);
199  double r01 = coeff(w, 0,1,true ) * Y(th,phi, 0,1,true );
200  double r02 = coeff(w, 0,2,true ) * Y(th,phi, 0,2,true );
201  double r11 = coeff(w, 1,1,true ) * Y(th,phi, 1,1,true );
202  double i11 = coeff(w, 1,1,false) * Y(th,phi, 1,1,false);
203  double r12 = coeff(w, 1,2,true ) * Y(th,phi, 1,2,true );
204  double i12 = coeff(w, 1,2,false) * Y(th,phi, 1,2,false);
205  double r22 = coeff(w, 2,2,true ) * Y(th,phi, 2,2,true );
206  double i22 = coeff(w, 2,2,false) * Y(th,phi, 2,2,false);
207 
208  nt.Fill(th,phi,w,sum,r01,r02,r11,i11,r12,i12,r22,i22);
209  }//w
210  }//phi
211  }//theta
212 
213  double wmin = 1.1;
214  double wmax = 2.0;
215  double nwb = 100;
216  TH1D cr01("cr01","coefficients of R{Y01} vs W",nwb,wmin,wmax);
217  TH1D cr02("cr02","coefficients of R{Y02} vs W",nwb,wmin,wmax);
218  TH1D cr11("cr11","coefficients of R{Y11} vs W",nwb,wmin,wmax);
219  TH1D ci11("ci11","coefficients of I{Y11} vs W",nwb,wmin,wmax);
220  TH1D cr12("cr12","coefficients of R{Y12} vs W",nwb,wmin,wmax);
221  TH1D ci12("ci12","coefficients of I{Y12} vs W",nwb,wmin,wmax);
222  TH1D cr22("cr22","coefficients of R{Y22} vs W",nwb,wmin,wmax);
223  TH1D ci22("ci22","coefficients of I{Y22} vs W",nwb,wmin,wmax);
224 
225  for(int i=1; i<cr01.GetNbinsX(); i++) {
226  double w = cr01.GetBinCenter(i);
227  cr01.Fill( w, coeff(w,0,1,true) );
228  cr02.Fill( w, coeff(w,0,2,true) );
229  cr11.Fill( w, coeff(w,1,1,true) );
230  ci11.Fill( w, coeff(w,1,1,false) );
231  cr12.Fill( w, coeff(w,1,2,true) );
232  ci12.Fill( w, coeff(w,1,2,false) );
233  cr22.Fill( w, coeff(w,2,2,true) );
234  ci22.Fill( w, coeff(w,2,2,false) );
235  }
236 
237  double dW[kNWBins];
238  for(int i=0; i<kNWBins; i++) {dW[i] = Wmax[i]-Wc[i];}
239 
240  TGraphAsymmErrors dr01(kNWBins,Wc,CoeffReY01,dW,dW,dCoeffReY01,dCoeffReY01);
241  TGraphAsymmErrors dr02(kNWBins,Wc,CoeffReY02,dW,dW,dCoeffReY02,dCoeffReY02);
242  TGraphAsymmErrors dr11(kNWBins,Wc,CoeffReY11,dW,dW,dCoeffReY11,dCoeffReY11);
243  TGraphAsymmErrors di11(kNWBins,Wc,CoeffImY11,dW,dW,dCoeffImY11,dCoeffImY11);
244  TGraphAsymmErrors dr12(kNWBins,Wc,CoeffReY12,dW,dW,dCoeffReY12,dCoeffReY12);
245  TGraphAsymmErrors di12(kNWBins,Wc,CoeffImY12,dW,dW,dCoeffImY12,dCoeffImY12);
246  TGraphAsymmErrors dr22(kNWBins,Wc,CoeffReY22,dW,dW,dCoeffReY22,dCoeffReY22);
247  TGraphAsymmErrors di22(kNWBins,Wc,CoeffImY22,dW,dW,dCoeffImY22,dCoeffImY22);
248 
249  dr01.Write("dr01");
250  dr02.Write("dr02");
251  dr11.Write("dr11");
252  di11.Write("di11");
253  dr12.Write("dr12");
254  di12.Write("di12");
255  dr22.Write("dr22");
256  di22.Write("di22");
257  cr01.Write();
258  cr02.Write();
259  cr11.Write();
260  ci11.Write();
261  cr12.Write();
262  ci12.Write();
263  cr22.Write();
264  ci22.Write();
265  nt.Write();
266  f.Close();
267 }
const double kPi
double Wc[kNWBins]
double dCoeffImY22[kNWBins]
double Wmax[kNWBins]
double CoeffImY12[kNWBins]
double CoeffReY02[kNWBins]
double dCoeffReY02[kNWBins]
double dCoeffReY22[kNWBins]
double dCoeffReY11[kNWBins]
double angular_distribution_pdf(double W, double th, double phi)
const int kNWBins
const double j
Definition: BetheBloch.cxx:29
double coeff(double W, int m, int l, bool real)
double CoeffImY11[kNWBins]
double Y(double th, double phi, int m, int l, bool real)
double CoeffReY22[kNWBins]
double dCoeffReY12[kNWBins]
double dCoeffImY11[kNWBins]
double CoeffReY12[kNWBins]
double CoeffReY01[kNWBins]
Double_t sum
Definition: plot.C:31
double CoeffImY22[kNWBins]
TNtuple * nt
Definition: drawXsec.C:2
double dCoeffImY12[kNWBins]
Float_t w
Definition: plot.C:20
#define W(x)
double CoeffReY11[kNWBins]
double dCoeffReY01[kNWBins]
void run_evgen_test ( double  W)

Definition at line 89 of file piN_asymmetric_decay.C.

References angular_distribution_pdf(), ana::assert(), MakeMiniprodValidationCuts::f, MECModelEnuComparisons::i, iev, calib::j, kPi, nt, make_associated_cosmic_defs::p4, RunSnowGlobes::pdf, confusionMatrixTree::t, and ana::weight.

90 {
91 // generate events with both the old & new scheme for comparison
92 //
93  TFile f("./events.root","recreate");
94  TNtuple nt("nt","","m:i:th:phi");
95 
96  TRandom rndgen(18928928);
97 
98  const int nev = 10000000;
99 
100  double th_min = 0;
101  double th_max = kPi;
102  double phi_min = 0;
103  double phi_max = 2*kPi;
104 
105 // scan for max
106 //
107  const int nth = 180;
108  const int nphi = 360;
109  double dth = (th_max - th_min ) / (nth -1);
110  double dphi = (phi_max - phi_min) / (nphi-1);
111 
112  double pdf_max = -1;
113  for(int i=0; i<nth; i++)
114  for(int j=0; j<nphi; j++)
115  pdf_max = TMath::Max(pdf_max,
116  angular_distribution_pdf(W,th_min+i*dth,phi_min+j*dphi));
117  pdf_max *= 1.2;
118 
119 // generate N theta,phi pairs using the new method
120 //
121  for(int iev=0; iev<nev; iev++) {
122  bool generated=false;
123  while(!generated) {
124  double th = rndgen.Rndm() * kPi;
125  double phi = rndgen.Rndm() * 2*kPi;
126  double t = rndgen.Rndm() * pdf_max;
127  double pdf = angular_distribution_pdf(W,th,phi);
128  generated = (t<pdf);
129  if(generated) { nt.Fill(1,iev,th,phi); }
130  }
131  }
132 
133 // generate N theta,phi pairs using the old method (phase space decay)
134 //
135  const int Np = 2;
136  const double MN = 0.938;
137  const double Mpi = 0.139;
138  TLorentzVector p4(0.0, 0.0, 0.0, W);
139  Double_t masses[Np] = { MN, Mpi } ;
140 
141  TGenPhaseSpace phspg;
142  Bool_t is_allowed = phspg.SetDecay(p4, Np, masses);
143  assert(is_allowed);
144 
145  double max_wght = -1;
146  for (Int_t iev=0; iev<200; iev++) {
147  double weight = phspg.Generate();
148  max_wght = TMath::Max(max_wght, weight);
149  }
150  max_wght *= 1.2;
151 
152  for(int iev=0; iev<nev; iev++) {
153  bool generated=false;
154  while(!generated) {
155  double weight = phspg.Generate();
156  double t = rndgen.Rndm() * max_wght;
157  generated = (t<weight);
158  if(generated) {
159  TLorentzVector * p4pi = phspg.GetDecay(1);
160  double th = p4pi->Theta();
161  double phi = p4pi->Phi() + kPi;
162  nt.Fill(-1,iev,th,phi);
163  }
164  }
165  }
166 
167  nt.Write();
168  f.Close();
169 }
int iev
Definition: runWimpSim.h:118
const double kPi
const Var weight
double angular_distribution_pdf(double W, double th, double phi)
const double j
Definition: BetheBloch.cxx:29
assert(nhit_max >=nhit_nbins)
TNtuple * nt
Definition: drawXsec.C:2
#define W(x)
TSpline3 splCoeffImY11 ( "splCoeffImY11"  ,
Wc  ,
CoeffImY11  ,
kNWBins  ,
""  ,
,
 
)

Referenced by coeff().

TSpline3 splCoeffImY12 ( "splCoeffImY12"  ,
Wc  ,
CoeffImY12  ,
kNWBins  ,
""  ,
,
 
)

Referenced by coeff().

TSpline3 splCoeffImY22 ( "splCoeffImY22"  ,
Wc  ,
CoeffImY22  ,
kNWBins  ,
""  ,
,
 
)

Referenced by coeff().

TSpline3 splCoeffReY01 ( "splCoeffReY01"  ,
Wc  ,
CoeffReY01  ,
kNWBins  ,
""  ,
,
 
)

Referenced by coeff().

TSpline3 splCoeffReY02 ( "splCoeffReY02"  ,
Wc  ,
CoeffReY02  ,
kNWBins  ,
""  ,
,
 
)

Referenced by coeff().

TSpline3 splCoeffReY11 ( "splCoeffReY11"  ,
Wc  ,
CoeffReY11  ,
kNWBins  ,
""  ,
,
 
)

Referenced by coeff().

TSpline3 splCoeffReY12 ( "splCoeffReY12"  ,
Wc  ,
CoeffReY12  ,
kNWBins  ,
""  ,
,
 
)

Referenced by coeff().

TSpline3 splCoeffReY22 ( "splCoeffReY22"  ,
Wc  ,
CoeffReY22  ,
kNWBins  ,
""  ,
,
 
)

Referenced by coeff().

double Y ( double  th,
double  phi,
int  m,
int  l,
bool  real 
)

Definition at line 307 of file piN_asymmetric_decay.C.

References Y00(), Y01(), Y02(), Y11(), Y12(), and Y22().

Referenced by angular_distribution_pdf(), and plot_input_data().

308 {
309  if (m==0 && l==0) return Y00(th,phi,real);
310  else if (m==0 && l==1) return Y01(th,phi,real);
311  else if (m==0 && l==2) return Y02(th,phi,real);
312  else if (m==1 && l==1) return Y11(th,phi,real);
313  else if (m==1 && l==2) return Y12(th,phi,real);
314  else if (m==2 && l==2) return Y22(th,phi,real);
315  else
316  return 0;
317 }
double Y12(double th, double phi, bool real)
double Y22(double th, double phi, bool real)
double Y00(double th, double phi, bool real)
double Y01(double th, double phi, bool real)
double Y11(double th, double phi, bool real)
double Y02(double th, double phi, bool real)
double Y00 ( double  th,
double  phi,
bool  real 
)

Definition at line 319 of file piN_asymmetric_decay.C.

References kPi, and ana::Sqrt().

Referenced by Y().

320 {
321  if(!real) return 0;
322  return 0.5/TMath::Sqrt(kPi);
323 }
const double kPi
Var Sqrt(const Var &v)
Use to take sqrt of a var.
Definition: Var.cxx:326
double Y01 ( double  th,
double  phi,
bool  real 
)

Definition at line 325 of file piN_asymmetric_decay.C.

References kPi, and ana::Sqrt().

Referenced by Y().

326 {
327  if(!real) return 0;
328  return 0.5 * TMath::Sqrt(3./kPi) * TMath::Cos(th);
329 }
const double kPi
Var Sqrt(const Var &v)
Use to take sqrt of a var.
Definition: Var.cxx:326
double Y02 ( double  th,
double  phi,
bool  real 
)

Definition at line 331 of file piN_asymmetric_decay.C.

References kPi, and ana::Sqrt().

Referenced by Y().

332 {
333  if(!real) return 0;
334  return 0.25 * TMath::Sqrt(5./kPi) * (3.* TMath::Power(TMath::Cos(th),2.) - 1);
335 }
const double kPi
Var Sqrt(const Var &v)
Use to take sqrt of a var.
Definition: Var.cxx:326
double Y11 ( double  th,
double  phi,
bool  real 
)

Definition at line 337 of file piN_asymmetric_decay.C.

References a, kPi, and ana::Sqrt().

Referenced by Y().

338 {
339  double a = -0.5 * TMath::Sqrt(1.5/kPi) * TMath::Sin(th);
340  if(!real) return a * TMath::Sin(phi);
341  else return a * TMath::Cos(phi);
342 }
const double kPi
const double a
Var Sqrt(const Var &v)
Use to take sqrt of a var.
Definition: Var.cxx:326
double Y12 ( double  th,
double  phi,
bool  real 
)

Definition at line 344 of file piN_asymmetric_decay.C.

References a, kPi, and ana::Sqrt().

Referenced by Y().

345 {
346  double a = -0.5 * TMath::Sqrt(7.5/kPi) * TMath::Cos(th) * TMath::Sin(th);
347  if(!real) return a * TMath::Sin(phi);
348  else return a * TMath::Cos(phi);
349 }
const double kPi
const double a
Var Sqrt(const Var &v)
Use to take sqrt of a var.
Definition: Var.cxx:326
double Y22 ( double  th,
double  phi,
bool  real 
)

Definition at line 351 of file piN_asymmetric_decay.C.

References a, kPi, and ana::Sqrt().

Referenced by Y().

352 {
353  double a = 0.25 * TMath::Sqrt(7.5/kPi) * TMath::Power( TMath::Sin(th), 2.);
354  if(!real) return a * TMath::Sin(2*phi);
355  else return a * TMath::Cos(2*phi);
356 }
const double kPi
const double a
Var Sqrt(const Var &v)
Use to take sqrt of a var.
Definition: Var.cxx:326

Variable Documentation

double CoeffImY11[kNWBins] = { -0.09, -0.06, -0.06, 0.01, -0.08, -0.17 }

Definition at line 46 of file piN_asymmetric_decay.C.

Referenced by plot_input_data().

double CoeffImY12[kNWBins] = { 0.07, 0.00, 0.16, -0.01, -0.04, -0.03 }

Definition at line 48 of file piN_asymmetric_decay.C.

Referenced by plot_input_data().

double CoeffImY22[kNWBins] = { -0.07, -0.07, -0.01, 0.06, -0.15, -0.03 }

Definition at line 50 of file piN_asymmetric_decay.C.

Referenced by plot_input_data().

double CoeffReY01[kNWBins] = { -0.22, -0.01, 0.29, 0.32, 0.61, 0.84 }

Definition at line 43 of file piN_asymmetric_decay.C.

Referenced by plot_input_data().

double CoeffReY02[kNWBins] = { -0.06, -0.13, -0.19, -0.08, 0.49, 0.80 }

Definition at line 44 of file piN_asymmetric_decay.C.

Referenced by plot_input_data().

double CoeffReY11[kNWBins] = { 0.21, -0.02, -0.20, -0.05, 0.12, -0.03 }

Definition at line 45 of file piN_asymmetric_decay.C.

Referenced by plot_input_data().

double CoeffReY12[kNWBins] = { 0.26, 0.05, 0.02, 0.00, 0.09, -0.05 }

Definition at line 47 of file piN_asymmetric_decay.C.

Referenced by plot_input_data().

double CoeffReY22[kNWBins] = { -0.04, 0.08, 0.08, 0.04, 0.09, -0.08 }

Definition at line 49 of file piN_asymmetric_decay.C.

Referenced by plot_input_data().

double dCoeffImY11[kNWBins] = { 0.08, 0.06, 0.08, 0.06, 0.09, 0.09 }

Definition at line 55 of file piN_asymmetric_decay.C.

Referenced by plot_input_data().

double dCoeffImY12[kNWBins] = { 0.07, 0.05, 0.08, 0.06, 0.08, 0.07 }

Definition at line 57 of file piN_asymmetric_decay.C.

Referenced by plot_input_data().

double dCoeffImY22[kNWBins] = { 0.08, 0.05, 0.08, 0.06, 0.08, 0.07 }

Definition at line 59 of file piN_asymmetric_decay.C.

Referenced by plot_input_data().

double dCoeffReY01[kNWBins] = { 0.11, 0.07, 0.09, 0.09, 0.13, 0.13 }

Definition at line 52 of file piN_asymmetric_decay.C.

Referenced by plot_input_data().

double dCoeffReY02[kNWBins] = { 0.12, 0.07, 0.10, 0.09, 0.15, 0.14 }

Definition at line 53 of file piN_asymmetric_decay.C.

Referenced by plot_input_data().

double dCoeffReY11[kNWBins] = { 0.09, 0.06, 0.08, 0.06, 0.08, 0.07 }

Definition at line 54 of file piN_asymmetric_decay.C.

Referenced by plot_input_data().

double dCoeffReY12[kNWBins] = { 0.08, 0.06, 0.08, 0.06, 0.07, 0.06 }

Definition at line 56 of file piN_asymmetric_decay.C.

Referenced by plot_input_data().

double dCoeffReY22[kNWBins] = { 0.07, 0.05, 0.07, 0.06, 0.09, 0.09 }

Definition at line 58 of file piN_asymmetric_decay.C.

Referenced by plot_input_data().

const int kNWBins = 6

Definition at line 38 of file piN_asymmetric_decay.C.

Referenced by plot_input_data().

const double kPi = 3.14

Definition at line 74 of file piN_asymmetric_decay.C.

Referenced by genie::HNIntranuke2018::AbsorbHN(), genie::MECGenerator::AddFinalStateLepton(), genie::AMNuGammaGenerator::AddPhoton(), genie::flux::GPointSourceAstroFlux::AddPointSource(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_BergerSehgalFM(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_ReinSehgal(), genie::COHKinematicsGenerator::CalculateKin_AlvarezRuso(), genie::SKKinematicsGenerator::CalculateKin_AtharSingleKaon(), CheckUnitarityLimit(), genie::NievesQELCCPXSec::CNCTCLimUcalc(), genie::SKKinematicsGenerator::ComputeMaxXSec(), genie::SmithMonizQELCCPXSec::d2sQES_dQ2dv_SM(), genie::SmithMonizQELCCPXSec::d3sQES_dQ2dvdkF_SM(), genie::mueloss::BetheBlochModel::dE_dx(), genie::KNOHadronization::DecayMethod2(), genie::NievesQELCCPXSec::deltaLindhard(), genie::alvarezruso::AlvarezRusoCOHPiPDXSec::DeltaWidthFree(), genie::alvarezruso::AlvarezRusoCOHPiPDXSec::DifferentialCrossSection(), genie::mueloss::gsl::BezrukovBugaevIntegrand::DoEval(), genie::mueloss::gsl::KokoulinPetrukhinIntegrand::DoEval(), genie::StrumiaVissaniIBDPXSec::dSigDt(), genie::SmithMonizQELCCPXSec::dsQES_dQ2_SM(), genie::flux::GFLUKAAtmoFlux::FillFluxHisto(), genie::flux::GBGLRSAtmoFlux::FillFluxHisto(), genie::StrumiaVissaniIBDPXSec::FinalStateCorr(), genie::LwlynSmithQELCCPXSec::FullDifferentialXSec(), genie::flux::GAtmoFlux::GenerateNext_1try(), genie::SpectralFunc::GenerateNucleon(), genie::EffectiveSF::GenerateNucleon(), genie::LocalFGM::GenerateNucleon(), genie::FGMBodekRitchie::GenerateNucleon(), genie::SpectralFunc1d::GenerateNucleon(), genie::flux::GCylindTH1Flux::GeneratePhi(), genie::VertexGenerator::GenerateVertex(), GetDetectorVertex(), genie::CharmHadronization::Hadronize(), genie::HAIntranuke::Inelastic(), genie::HAIntranuke2018::Inelastic(), genie::GLRESPXSec::Integral(), genie::COHXSecAR::Integrate(), genie::NewQELXSec::Integrate(), genie::COHKinematicsGenerator::MaxXSec_AlvarezRuso(), genie::utils::intranuke::PhaseSpaceDecay(), genie::utils::intranuke2018::PhaseSpaceDecay(), genie::NucDeExcitationSim::Photon4P(), plot_input_data(), genie::FGMBodekRitchie::ProbDistro(), genie::LocalFGM::ProbDistro(), genie::DFRHadronicSystemGenerator::ProcessEventRecord(), genie::NuEPrimaryLeptonGenerator::ProcessEventRecord(), genie::PrimaryLeptonGenerator::ProcessEventRecord(), genie::OutgoingDarkGenerator::ProcessEventRecord(), genie::QELEventGenerator::ProcessEventRecord(), genie::StrumiaVissaniIBDPXSec::RadiativeCorr(), genie::NievesQELCCPXSec::relLindhardIm(), run_evgen_test(), genie::MECGenerator::SelectNSVLeptonKinematics(), genie::flux::GAstroFlux::NuGenerator::SelectOrigin(), genie::flux::GFLUKAAtmoFlux::SetBinSizes(), genie::flux::GBGLRSAtmoFlux::SetBinSizes(), genie::flux::GAstroFlux::SetDetectorPosition(), genie::GHepParticle::SetPolarization(), genie::utils::intranuke::ThreeBodyKinematics(), genie::utils::intranuke2018::ThreeBodyKinematics(), genie::utils::intranuke::TwoBodyKinematics(), genie::utils::intranuke2018::TwoBodyKinematics(), genie::NievesQELCCPXSec::vcr(), genie::SlowRsclCharmDISPXSecLO::XSec(), genie::AlamSimoAtharVacasSKPXSec2014::XSec(), genie::AhrensNCELPXSec::XSec(), genie::RosenbluthPXSec::XSec(), genie::AivazisCharmPXSecLO::XSec(), genie::IMDAnnihilationPXSec::XSec(), genie::QPMDISPXSec::XSec(), genie::P33PaschosLalakulichPXSec::XSec(), genie::AhrensDMELPXSec::XSec(), genie::BergerSehgalFMCOHPiPXSec2015::XSec(), genie::LwlynSmithQELCCPXSec::XSec(), genie::QPMDMDISPXSec::XSec(), genie::BergerSehgalCOHPiPXSec2015::XSec(), genie::NuElectronPXSec::XSec(), genie::PaisQELLambdaPXSec::XSec(), genie::ReinSehgalCOHPiPXSec::XSec(), genie::KovalenkoQELCharmPXSec::XSec(), genie::BardinIMDRadCorPXSec::XSec(), genie::ReinSehgalRESPXSec::XSec(), genie::BSKLNBaseRESPXSec2014::XSec(), genie::NievesQELCCPXSec::XSec(), Y00(), Y01(), Y02(), Y11(), Y12(), and Y22().

double Wc[kNWBins] = { 1.140, 1.205, 1.255, 1.390, 1.600, 1.850 }

Definition at line 41 of file piN_asymmetric_decay.C.

Referenced by plot_input_data().

double Wmax[kNWBins] = { 1.180, 1.230, 1.280, 1.500, 1.700, 2.000 }
double Wmin[kNWBins] = { 1.100, 1.180, 1.230, 1.280, 1.500, 1.700 }