Numu2017Fits.h File Reference

Go to the source code of this file.

Functions
double	correct_prod3_fd_muon_length (double trkLen, bool ismc)
double	correct_prod3_fd_calE_low_gain (double calE, bool ismc)
double	correct_prod3_fd_calE_high_gain (double calE, bool ismc)
double	correct_prod3_fd_calE (double calE, int run, bool ismc)
double	predict_prod3_fd_p1_muon_energy (double trkLen)
double	predict_prod3_fd_p2_muon_energy (double trkLen)
double	predict_prod3_fd_p3_muon_energy (double trkLen)
double	predict_prod3_fd_p4_muon_energy (double trkLen)
double	predict_prod3_fd_p5_muon_energy (double trkLen)
double	predict_prod3_fd_p1_had_energy (double hadVisE)
double	predict_prod3_fd_p2_had_energy (double hadVisE)
double	predict_prod3_fd_p3_had_energy (double hadVisE)
double	predict_prod3_fd_p4_had_energy (double hadVisE)
double	predict_prod3_fd_p5_had_energy (double hadVisE)
double	predict_prod3_fd_p1_muon_energy_corr (double trkLen, bool ismc)
double	predict_prod3_fd_p2_muon_energy_corr (double trkLen, bool ismc)
double	predict_prod3_fd_p3_muon_energy_corr (double trkLen, bool ismc)
double	predict_prod3_fd_p4_muon_energy_corr (double trkLen, bool ismc)
double	predict_prod3_fd_p5_muon_energy_corr (double trkLen, bool ismc)
double	predict_prod3_fd_p1_had_energy_corr (double hadVisE, bool ismc)
double	predict_prod3_fd_p2_had_energy_corr (double hadVisE, bool ismc)
double	predict_prod3_fd_p3_had_energy_corr (double hadVisE, bool ismc)
double	predict_prod3_fd_p4_had_energy_corr (double hadVisE, bool ismc)
double	predict_prod3_fd_p5_had_energy_corr (double hadVisE, bool ismc)
double	predict_prod3_nd_p3_act_energy (double ndtrklenact)
double	predict_prod3_nd_p3_cat_energy (double ndtrklencat)
double	predict_prod3_nd_p3_had_energy (double hadVisE)
double	predict_mprod16_nd_p4_act_energy (double ndtrklenact)
double	predict_mprod16_nd_p4_cat_energy (double ndtrklencat)
double	predict_mprod16_nd_p4_had_energy (double hadVisE)
double	predict_prod3_nd_p4_act_energy (double ndtrklenact)
double	predict_prod3_nd_p4_cat_energy (double ndtrklencat)
double	predict_prod3_nd_p4_had_energy (double hadVisE)
double	predict_prod3_fd_muon_energy (double trkLen, int run, bool ismc, bool isRHC)
double	predict_prod3_fd_had_energy (double hadVisE, int run, bool ismc, bool isRHC)
double	predict_prod3_nd_act_energy (double ndtrklenact, bool isRHC)
double	predict_prod3_nd_cat_energy (double ndtrklencat, bool isRHC)
double	predict_prod3_nd_had_energy (double hadVisE, bool isRHC)

Function Documentation

double correct_prod3_fd_calE	(	double	calE,
		int	run,
		bool	ismc
	)

Definition at line 132 of file Numu2017Fits.cxx.

References fd_period_2_run_end.

Referenced by predict_prod3_fd_had_energy(), and ana::predict_special_fd_had_energy().

{
  /* Calorimetric energy miscalibration fix: doc-23372-v2 */
  if (ismc) {
    return calE;
  }

  if (run <= fd_period_2_run_end) {
    /* low gain correction */
    return calE * 0.9844;
  }

  /* high gain correction */
  return calE * 0.9949;
}

double correct_prod3_fd_calE_high_gain	(	double	calE,
		bool	ismc
	)

Definition at line 119 of file Numu2017Fits.cxx.

Referenced by predict_prod3_fd_p3_had_energy_corr(), predict_prod3_fd_p4_had_energy_corr(), predict_prod3_fd_p5_had_energy_corr(), NumuEnergyFunc::predict_prod4_fd_had_energy_p(), ana::predict_special_fd_p3_had_energy_corr(), and ana::predict_special_fd_p5_had_energy_corr().

{
  /* Calorimetric energy miscalibration fix: doc-23372-v2 */
  if (ismc) {
    return calE;
  }

  /* high gain correction */
  return calE * 0.9949;
}

double correct_prod3_fd_calE_low_gain	(	double	calE,
		bool	ismc
	)

Definition at line 106 of file Numu2017Fits.cxx.

Referenced by predict_prod3_fd_p1_had_energy_corr(), predict_prod3_fd_p2_had_energy_corr(), NumuEnergyFunc::predict_prod4_fd_had_energy_p(), ana::predict_special_fd_p1_had_energy_corr(), and ana::predict_special_fd_p2_had_energy_corr().

{
  /* Calorimetric energy miscalibration fix: doc-23372-v2 */
  if (ismc) {
    return calE;
  }

  /* low gain correction */
  return calE * 0.9844;
}

double correct_prod3_fd_muon_length	(	double	trkLen,
		bool	ismc
	)

Definition at line 94 of file Numu2017Fits.cxx.

Referenced by predict_prod3_fd_muon_energy(), predict_prod3_fd_p1_muon_energy_corr(), predict_prod3_fd_p2_muon_energy_corr(), predict_prod3_fd_p3_muon_energy_corr(), predict_prod3_fd_p4_muon_energy_corr(), predict_prod3_fd_p5_muon_energy_corr(), NumuEnergyFunc::predict_prod4_fd_muon_energy_p(), NumuEnergyFunc::predict_prod5_fd_muon_energy_p(), ana::predict_special_fd_muon_energy(), ana::predict_special_fd_p1_muon_energy_corr(), ana::predict_special_fd_p2_muon_energy_corr(), ana::predict_special_fd_p3_muon_energy_corr(), and ana::predict_special_fd_p5_muon_energy_corr().

{
  /* muon energy scale correction doc-20816-v6 */
  if (ismc) {
    return trkLen;
  }

  return trkLen * 0.9957;
}

double predict_mprod16_nd_p4_act_energy

(

double

ndtrklenact

)

Definition at line 515 of file Numu2017Fits.cxx.

References PandAna.reco_validation.add_data::offset, and spline3().

{
  /* Note, this fit was done on the miniprod 2016 sample doc-18966-v2 */

  const double offset  = 1.560502711152023103e-01;
  const double slope0  = 1.927996549217538824e-01;
  const double stitch1 = 2.763097394311628818e+00;
  const double slope1  = 2.000138250885013236e-01;
  const double stitch2 = 4.649293917550735955e+00;
  const double slope2  = 2.075538808578699379e-01;
  const double stitch3 = 9.095060911026711992e+00;
  const double slope3  = 2.145950730831300002e-01;

  if (ndtrklenact <= 0) {
    return 0;
  }

  return spline3(ndtrklenact / 100,
                 offset,  slope0,
                 stitch1, slope1,
                 stitch2, slope2,
                 stitch3, slope3);
}

double predict_mprod16_nd_p4_cat_energy

(

double

ndtrklencat

)

Definition at line 541 of file Numu2017Fits.cxx.

References PandAna.reco_validation.add_data::offset, and spline3().

{  
  /* Note, this fit was done on the miniprod 2016 sample doc-18966-v2 */

  const double offset  = 7.059293364276031379e-02;
  const double slope0  = 1.052724154897643205e-01;
  const double stitch1 = 1.466666666666676111e-01;
  const double slope1  = 5.555668682049775819e-01;
  const double stitch2 = 2.274316786094766840e+00;
  const double slope2  = 2.015720164137814852e+00;
  const double stitch3 = 2.301657526094746054e+00;
  const double slope3  = 3.938167768948175884e-01;

  if (ndtrklencat <= 0) {
    return 0;
  }

  return spline3(ndtrklencat / 100,
                 offset,  slope0,
                 stitch1, slope1,
                 stitch2, slope2,
                 stitch3, slope3);
}

double predict_mprod16_nd_p4_had_energy

(

double

hadVisE

)

Definition at line 567 of file Numu2017Fits.cxx.

References PandAna.reco_validation.add_data::offset, and spline3().

{
  /* Note, this fit was done on the miniprod 2016 sample doc-18966-v2 */

  const double offset  = 4.370690101565655183e-02;
  const double slope0  = 2.306317746493836029e+00;
  const double stitch1 = 5.046916149809789476e-02;
  const double slope1  = 3.250145515659832895e+00;
  const double stitch2 = 1.021625084023882746e-01;
  const double slope2  = 1.624362074506924047e+00;
  const double stitch3 = 1.684459951617373330e-01;
  const double slope3  = 1.892049412037302369e+00;

  if (hadVisE <= 0) {
    return 0;
  }

  return spline3(hadVisE,
                 offset,  slope0,
                 stitch1, slope1,
                 stitch2, slope2,
                 stitch3, slope3);
}

double predict_prod3_fd_had_energy	(	double	hadVisE,
		int	run,
		bool	ismc,
		bool	isRHC
	)

Definition at line 688 of file Numu2017Fits.cxx.

References correct_prod3_fd_calE(), fd_period_1_run_end, fd_period_2_run_end, fd_period_3_run_end, fd_period_4_run_end, fd_period_5_run_end, predict_prod3_fd_p1_had_energy(), predict_prod3_fd_p2_had_energy(), predict_prod3_fd_p3_had_energy(), predict_prod3_fd_p4_had_energy(), and predict_prod3_fd_p5_had_energy().

Referenced by ana::TAHadEND().

{
  hadVisE = correct_prod3_fd_calE(hadVisE, run, ismc);

  if (run <= fd_period_1_run_end) {
    return predict_prod3_fd_p1_had_energy(hadVisE);
  }

  if (run <= fd_period_2_run_end) {
    return predict_prod3_fd_p2_had_energy(hadVisE);
  }

  if (run <= fd_period_3_run_end) {
    return predict_prod3_fd_p3_had_energy(hadVisE);
  }

  if (run <= fd_period_4_run_end) {
    return predict_prod3_fd_p4_had_energy(hadVisE);
  }

  if (run <= fd_period_5_run_end) {
    return predict_prod3_fd_p5_had_energy(hadVisE);
  }

  /* period 5 was the last period 2017 energy estimators were trained on  */

  if (isRHC) {
    return predict_prod3_fd_p4_had_energy(hadVisE);
  }

  return predict_prod3_fd_p5_had_energy(hadVisE);
}

double predict_prod3_fd_muon_energy	(	double	trkLen,
		int	run,
		bool	ismc,
		bool	isRHC
	)

Definition at line 652 of file Numu2017Fits.cxx.

References correct_prod3_fd_muon_length(), fd_period_1_run_end, fd_period_2_run_end, fd_period_3_run_end, fd_period_4_run_end, fd_period_5_run_end, predict_prod3_fd_p1_muon_energy(), predict_prod3_fd_p2_muon_energy(), predict_prod3_fd_p3_muon_energy(), predict_prod3_fd_p4_muon_energy(), and predict_prod3_fd_p5_muon_energy().

Referenced by ana::TAHadEND().

{
  trkLen = correct_prod3_fd_muon_length(trkLen, ismc);

  if (run <= fd_period_1_run_end) {
    return predict_prod3_fd_p1_muon_energy(trkLen);
  }

  if (run <= fd_period_2_run_end) {
    return predict_prod3_fd_p2_muon_energy(trkLen);
  }

  if (run <= fd_period_3_run_end) {
    return predict_prod3_fd_p3_muon_energy(trkLen);
  }

  if (run <= fd_period_4_run_end) {
    return predict_prod3_fd_p4_muon_energy(trkLen);
  }

  if (run <= fd_period_5_run_end) {
    return predict_prod3_fd_p5_muon_energy(trkLen);
  }

  /* period 5 was the last period 2017 energy estimators were trained on  */

  if (isRHC) {
    return predict_prod3_fd_p4_muon_energy(trkLen);
  }

  return predict_prod3_fd_p5_muon_energy(trkLen);
}

double predict_prod3_fd_p1_had_energy

(

double

hadVisE

)

Definition at line 258 of file Numu2017Fits.cxx.

References PandAna.reco_validation.add_data::offset, and spline3().

Referenced by predict_prod3_fd_had_energy(), and predict_prod3_fd_p1_had_energy_corr().

{
  const double offset  = 6.069740170932858003e-02;
  const double slope0  = 1.103685627830206739e+00;
  const double stitch1 = 8.221169483726585570e-02;
  const double slope1  = 2.094686049645236192e+00;
  const double stitch2 = 5.335415728146151171e-01;
  const double slope2  = 1.649382412059058556e+00;
  const double stitch3 = 9.849999979737633637e-01;
  const double slope3  = 2.153569995620331134e+00;

  if (hadVisE <= 0) {
    return 0;
  }

  return spline3(hadVisE,
                 offset,  slope0,
                 stitch1, slope1,
                 stitch2, slope2,
                 stitch3, slope3);
}

double predict_prod3_fd_p1_had_energy_corr	(	double	hadVisE,
		bool	ismc
	)

Definition at line 475 of file Numu2017Fits.cxx.

References correct_prod3_fd_calE_low_gain(), and predict_prod3_fd_p1_had_energy().

Referenced by ana::TAHadEND().

{
  hadVisE = correct_prod3_fd_calE_low_gain(hadVisE, ismc);
  return predict_prod3_fd_p1_had_energy(hadVisE);
}

double predict_prod3_fd_p1_muon_energy

(

double

trkLen

)

Definition at line 150 of file Numu2017Fits.cxx.

References PandAna.reco_validation.add_data::offset, and spline2().

Referenced by predict_prod3_fd_muon_energy(), and predict_prod3_fd_p1_muon_energy_corr().

{
  const double offset  = 1.483250632870023455e-01;
  const double slope0  = 1.955739883601782803e-01;
  const double stitch1 = 4.129999078650119060e+00;
  const double slope1  = 2.057095056740816386e-01;
  const double stitch2 = 9.414831896417680923e+00;
  const double slope2  = 2.147448462294846983e-01;

  if (trkLen <= 0) {
    return 0;
  }

  return spline2(trkLen / 100,
                 offset,  slope0,
                 stitch1, slope1,
                 stitch2, slope2);
}

double predict_prod3_fd_p1_muon_energy_corr	(	double	trkLen,
		bool	ismc
	)

Definition at line 435 of file Numu2017Fits.cxx.

References correct_prod3_fd_muon_length(), and predict_prod3_fd_p1_muon_energy().

Referenced by ana::TAHadEND().

{
  trkLen = correct_prod3_fd_muon_length(trkLen, ismc);
  return predict_prod3_fd_p1_muon_energy(trkLen);
}

double predict_prod3_fd_p2_had_energy

(

double

hadVisE

)

Definition at line 282 of file Numu2017Fits.cxx.

References PandAna.reco_validation.add_data::offset, and spline3().

Referenced by predict_prod3_fd_had_energy(), and predict_prod3_fd_p2_had_energy_corr().

{
  const double offset  = 3.487649857913088791e-02;
  const double slope0  = 2.109632444227325987e+00;
  const double stitch1 = 5.092713207810132259e-01;
  const double slope1  = 1.716494894914534886e+00;
  const double stitch2 = 9.449999998699000203e-01;
  const double slope2  = 2.338916725792882456e+00;
  const double stitch3 = 1.351608484611389382e+00;
  const double slope3  = 1.670675446405347708e+00;

  if (hadVisE <= 0) {
    return 0;
  }

  return spline3(hadVisE,
                 offset,  slope0,
                 stitch1, slope1,
                 stitch2, slope2,
                 stitch3, slope3);
}

double predict_prod3_fd_p2_had_energy_corr	(	double	hadVisE,
		bool	ismc
	)

Definition at line 483 of file Numu2017Fits.cxx.

References correct_prod3_fd_calE_low_gain(), and predict_prod3_fd_p2_had_energy().

Referenced by ana::TAHadEND().

{
  hadVisE = correct_prod3_fd_calE_low_gain(hadVisE, ismc);
  return predict_prod3_fd_p2_had_energy(hadVisE);
}

double predict_prod3_fd_p2_muon_energy

(

double

trkLen

)

Definition at line 171 of file Numu2017Fits.cxx.

References PandAna.reco_validation.add_data::offset, and spline2().

Referenced by predict_prod3_fd_muon_energy(), and predict_prod3_fd_p2_muon_energy_corr().

{
  const double offset  = 1.532385212574850497e-01;
  const double slope0  = 1.943835976098556295e-01;
  const double stitch1 = 3.854492640068833520e+00;
  const double slope1  = 2.057372525501594829e-01;
  const double stitch2 = 9.882283054843421510e+00;
  const double slope2  = 2.182658523728695510e-01;

  if (trkLen <= 0) {
    return 0;
  }

  return spline2(trkLen / 100,
                 offset,  slope0,
                 stitch1, slope1,
                 stitch2, slope2);
}

double predict_prod3_fd_p2_muon_energy_corr	(	double	trkLen,
		bool	ismc
	)

Definition at line 443 of file Numu2017Fits.cxx.

References correct_prod3_fd_muon_length(), and predict_prod3_fd_p2_muon_energy().

Referenced by ana::TAHadEND().

{
  trkLen = correct_prod3_fd_muon_length(trkLen, ismc);
  return predict_prod3_fd_p2_muon_energy(trkLen);
}

double predict_prod3_fd_p3_had_energy

(

double

hadVisE

)

Definition at line 306 of file Numu2017Fits.cxx.

References PandAna.reco_validation.add_data::offset, and spline3().

Referenced by predict_prod3_fd_had_energy(), and predict_prod3_fd_p3_had_energy_corr().

{
  const double offset  = 4.729065467246984511e-02;
  const double slope0  = 1.931526076780770929e+00;
  const double stitch1 = 2.799421333699638192e-01;
  const double slope1  = 2.266534645077610488e+00;
  const double stitch2 = 4.692746254578271659e-01;
  const double slope2  = 1.648693959937136633e+00;
  const double stitch3 = 8.921135013900192146e-01;
  const double slope3  = 2.020802581638818651e+00;

  if (hadVisE <= 0) {
    return 0;
  }

  return spline3(hadVisE,
                 offset,  slope0,
                 stitch1, slope1,
                 stitch2, slope2,
                 stitch3, slope3);
}

double predict_prod3_fd_p3_had_energy_corr	(	double	hadVisE,
		bool	ismc
	)

Definition at line 491 of file Numu2017Fits.cxx.

References correct_prod3_fd_calE_high_gain(), and predict_prod3_fd_p3_had_energy().

Referenced by ana::TAHadEND().

{
  hadVisE = correct_prod3_fd_calE_high_gain(hadVisE, ismc);
  return predict_prod3_fd_p3_had_energy(hadVisE);
}

double predict_prod3_fd_p3_muon_energy

(

double

trkLen

)

Definition at line 192 of file Numu2017Fits.cxx.

References PandAna.reco_validation.add_data::offset, and spline2().

Referenced by predict_prod3_fd_muon_energy(), and predict_prod3_fd_p3_muon_energy_corr().

{
  const double offset  = 1.476537053972520885e-01;
  const double slope0  = 1.960472718145106208e-01;
  const double stitch1 = 4.120072099634467300e+00;
  const double slope1  = 2.059084783219390091e-01;
  const double stitch2 = 9.874976681687856583e+00;
  const double slope2  = 2.178808582731997512e-01;

  if (trkLen <= 0) {
    return 0;
  }

  return spline2(trkLen / 100,
                 offset,  slope0,
                 stitch1, slope1,
                 stitch2, slope2);
}

double predict_prod3_fd_p3_muon_energy_corr	(	double	trkLen,
		bool	ismc
	)

Definition at line 451 of file Numu2017Fits.cxx.

References correct_prod3_fd_muon_length(), and predict_prod3_fd_p3_muon_energy().

Referenced by ana::TAHadEND().

{
  trkLen = correct_prod3_fd_muon_length(trkLen, ismc);
  return predict_prod3_fd_p3_muon_energy(trkLen);
}

double predict_prod3_fd_p4_had_energy

(

double

hadVisE

)

Definition at line 330 of file Numu2017Fits.cxx.

References PandAna.reco_validation.add_data::offset, and spline3().

Referenced by predict_prod3_fd_had_energy(), and predict_prod3_fd_p4_had_energy_corr().

{
  const double offset  = 4.332202531604933959e-02;
  const double slope0  = 2.366131159456818178e+00;
  const double stitch1 = 1.049999756567459169e-01;
  const double slope1  = 1.897994466055304974e+00;
  const double stitch2 = 4.598347475287963282e-01;
  const double slope2  = 1.574651945578143808e+00;
  const double stitch3 = 6.601961854211407799e-01;
  const double slope3  = 2.086566751147730159e+00;

  if (hadVisE <= 0) {
    return 0;
  }

  return spline3(hadVisE,
                 offset,  slope0,
                 stitch1, slope1,
                 stitch2, slope2,
                 stitch3, slope3);
}

double predict_prod3_fd_p4_had_energy_corr	(	double	hadVisE,
		bool	ismc
	)

Definition at line 499 of file Numu2017Fits.cxx.

References correct_prod3_fd_calE_high_gain(), and predict_prod3_fd_p4_had_energy().

Referenced by ana::TAHadEND().

{
  hadVisE = correct_prod3_fd_calE_high_gain(hadVisE, ismc);
  return predict_prod3_fd_p4_had_energy(hadVisE);
}

double predict_prod3_fd_p4_muon_energy

(

double

trkLen

)

Definition at line 213 of file Numu2017Fits.cxx.

References PandAna.reco_validation.add_data::offset, and spline3().

Referenced by predict_prod3_fd_muon_energy(), and predict_prod3_fd_p4_muon_energy_corr().

{
  const double offset  = 1.543934877851929066e-01;
  const double slope0  = 1.939411335845932782e-01;
  const double stitch1 = 3.187688625921237140e+00;
  const double slope1  = 2.020725882614244073e-01;
  const double stitch2 = 5.031281501272647283e+00;
  const double slope2  = 2.072840009844969800e-01;
  const double stitch3 = 1.031568838438018076e+01;
  const double slope3  = 2.182695636807331363e-01;

  if (trkLen <= 0) {
    return 0;
  }

  return spline3(trkLen / 100,
                 offset,  slope0,
                 stitch1, slope1,
                 stitch2, slope2,
                 stitch3, slope3);
}

double predict_prod3_fd_p4_muon_energy_corr	(	double	trkLen,
		bool	ismc
	)

Definition at line 459 of file Numu2017Fits.cxx.

References correct_prod3_fd_muon_length(), and predict_prod3_fd_p4_muon_energy().

Referenced by ana::TAHadEND().

{
  trkLen = correct_prod3_fd_muon_length(trkLen, ismc);
  return predict_prod3_fd_p4_muon_energy(trkLen);
}

double predict_prod3_fd_p5_had_energy

(

double

hadVisE

)

Definition at line 354 of file Numu2017Fits.cxx.

References PandAna.reco_validation.add_data::offset, and spline2().

Referenced by predict_prod3_fd_had_energy(), and predict_prod3_fd_p5_had_energy_corr().

{
  const double offset  = 3.442138256673876739e-02;
  const double slope0  = 2.036769676363937354e+00;
  const double stitch1 = 5.603659885239966609e-01;
  const double slope1  = 1.585650161781588352e+00;
  const double stitch2 = 8.493299384706546329e-01;
  const double slope2  = 2.050073571966713359e+00;

  if (hadVisE <= 0) {
    return 0;
  }

  return spline2(hadVisE,
                 offset,  slope0,
                 stitch1, slope1,
                 stitch2, slope2);
}

double predict_prod3_fd_p5_had_energy_corr	(	double	hadVisE,
		bool	ismc
	)

Definition at line 507 of file Numu2017Fits.cxx.

References correct_prod3_fd_calE_high_gain(), and predict_prod3_fd_p5_had_energy().

Referenced by ana::TAHadEND().

{
  hadVisE = correct_prod3_fd_calE_high_gain(hadVisE, ismc);
  return predict_prod3_fd_p5_had_energy(hadVisE);
}

double predict_prod3_fd_p5_muon_energy

(

double

trkLen

)

Definition at line 237 of file Numu2017Fits.cxx.

References PandAna.reco_validation.add_data::offset, and spline2().

Referenced by predict_prod3_fd_muon_energy(), and predict_prod3_fd_p5_muon_energy_corr().

{
  const double offset  = 1.488648665093454948e-01;
  const double slope0  = 1.955393334139079176e-01;
  const double stitch1 = 3.982132354846828814e+00;
  const double slope1  = 2.056853517430136113e-01;
  const double stitch2 = 9.610128112133583755e+00;
  const double slope2  = 2.178689785257645617e-01;

  if (trkLen <= 0) {
    return 0;
  }

  return spline2(trkLen / 100,
                 offset,  slope0,
                 stitch1, slope1,
                 stitch2, slope2);
}

double predict_prod3_fd_p5_muon_energy_corr	(	double	trkLen,
		bool	ismc
	)

Definition at line 467 of file Numu2017Fits.cxx.

References correct_prod3_fd_muon_length(), and predict_prod3_fd_p5_muon_energy().

Referenced by ana::TAHadEND().

{
  trkLen = correct_prod3_fd_muon_length(trkLen, ismc);
  return predict_prod3_fd_p5_muon_energy(trkLen);
}

double predict_prod3_nd_act_energy	(	double	ndtrklenact,
		bool	isRHC
	)

Definition at line 724 of file Numu2017Fits.cxx.

References predict_prod3_nd_p3_act_energy(), and predict_prod3_nd_p4_act_energy().

Referenced by ana::TAHadEND().

{
  if (isRHC) {
    return predict_prod3_nd_p4_act_energy(ndtrklenact);
  }

  return predict_prod3_nd_p3_act_energy(ndtrklenact);
}

double predict_prod3_nd_cat_energy	(	double	ndtrklencat,
		bool	isRHC
	)

Definition at line 735 of file Numu2017Fits.cxx.

References predict_prod3_nd_p3_cat_energy(), and predict_prod3_nd_p4_cat_energy().

Referenced by ana::TAHadEND().

{
  if (isRHC) {
    return predict_prod3_nd_p4_cat_energy(ndtrklencat);
  }

  return predict_prod3_nd_p3_cat_energy(ndtrklencat);
}

double predict_prod3_nd_had_energy	(	double	hadVisE,
		bool	isRHC
	)

Definition at line 746 of file Numu2017Fits.cxx.

References predict_prod3_nd_p3_had_energy(), and predict_prod3_nd_p4_had_energy().

Referenced by ana::TAHadEND().

{
  if (isRHC) {
    return predict_prod3_nd_p4_had_energy(hadVisE);
  }

  return predict_prod3_nd_p3_had_energy(hadVisE);
}

double predict_prod3_nd_p3_act_energy

(

double

ndtrklenact

)

Definition at line 375 of file Numu2017Fits.cxx.

References PandAna.reco_validation.add_data::offset, and spline1().

Referenced by predict_prod3_nd_act_energy(), and ana::TAHadEND().

{
  const double offset  = 1.473446267714921820e-01;
  const double slope0  = 1.958813916199254879e-01;
  const double stitch1 = 4.117512435079445865e+00;
  const double slope1  = 2.068093764675576485e-01;

  if (ndtrklenact <= 0) {
    return 0;
  }

  return spline1(ndtrklenact / 100,
                 offset,  slope0,
                 stitch1, slope1);
}

double predict_prod3_nd_p3_cat_energy

(

double

ndtrklencat

)

Definition at line 393 of file Numu2017Fits.cxx.

References PandAna.reco_validation.add_data::offset, and spline3().

Referenced by predict_prod3_nd_cat_energy(), and ana::TAHadEND().

{
  const double offset  = 7.017118321474302700e-02;
  const double slope0  = 1.009985063097527530e-01;
  const double stitch1 = 1.466666666666738839e-01;
  const double slope1  = 5.558686386480078401e-01;
  const double stitch2 = 2.269449191147328015e+00;
  const double slope2  = 1.506386606204337841e+00;
  const double stitch3 = 2.310334083126842852e+00;
  const double slope3  = 3.634575525654349382e-01;

  if (ndtrklencat <= 0) {
    return 0;
  }

  return spline3(ndtrklencat / 100,
                 offset,  slope0,
                 stitch1, slope1,
                 stitch2, slope2,
                 stitch3, slope3);
}

double predict_prod3_nd_p3_had_energy

(

double

hadVisE

)

Definition at line 417 of file Numu2017Fits.cxx.

References PandAna.reco_validation.add_data::offset, and spline1().

Referenced by predict_prod3_nd_had_energy(), and ana::TAHadEND().

{
  const double offset  = 3.766323500934987134e-02;
  const double slope0  = 2.206130121197688698e+00;
  const double stitch1 = 5.005521534140644491e-01;
  const double slope1  = 1.775841650258208393e+00;

  if (hadVisE <= 0) {
    return 0;
  }

  return spline1(hadVisE,
                 offset,  slope0,
                 stitch1, slope1);
}

double predict_prod3_nd_p4_act_energy

(

double

ndtrklenact

)

Definition at line 593 of file Numu2017Fits.cxx.

References PandAna.reco_validation.add_data::offset, and spline1().

Referenced by predict_prod3_nd_act_energy(), and ana::TAHadEND().

{
  const double offset  = 1.456667491802539072e-01;
  const double slope0  = 1.971031247457571411e-01;
  const double stitch1 = 4.359417948299743273e+00;
  const double slope1  = 2.078006187845642394e-01;

  if (ndtrklenact <= 0) {
    return 0;
  }

  return spline1(ndtrklenact / 100,
                 offset,  slope0,
                 stitch1, slope1);
}

double predict_prod3_nd_p4_cat_energy

(

double

ndtrklencat

)

Definition at line 611 of file Numu2017Fits.cxx.

References PandAna.reco_validation.add_data::offset, and spline0().

Referenced by predict_prod3_nd_cat_energy(), and ana::TAHadEND().

{  
  /* 
   * NOTE: i really do not like a negative offset here, and overall this fit 
   *       is much worse than similar fit on miniprod 2016
   */

  const double offset = -2.247455892420480161e-03;
  const double slope0 =  5.610243972312459748e-01;

  if (ndtrklencat <= 0) {
    return 0;
  }

  return spline0(ndtrklencat / 100,
                 offset,  slope0);
}

double predict_prod3_nd_p4_had_energy

(

double

hadVisE

)

Definition at line 631 of file Numu2017Fits.cxx.

References PandAna.reco_validation.add_data::offset, and spline2().

Referenced by predict_prod3_nd_had_energy(), and ana::TAHadEND().

{
  const double offset  = 4.483967184110027038e-02;
  const double slope0  = 2.965814310246884133e+00;
  const double stitch1 = 9.767960955481164165e-02;
  const double slope1  = 1.531955866917958486e+00;
  const double stitch2 = 1.708985997225620779e-01;
  const double slope2  = 1.944014774753460850e+00;

  if (hadVisE <= 0) {
    return 0;
  }

  return spline2(hadVisE,
                 offset,  slope0,
                 stitch1, slope1,
                 stitch2, slope2);
}