PRD2020/ResolutionFunctions.hh

Go to the documentation of this file.

double Hybrids_ICRC_2019_SD1500(const double lgE) {
  /* Preliminary version for ICRC 2019 */
  // V. Verzi 03 04 2019
  // Resolution calculated with hybrids data

  double E = pow(10, lgE);
  double res = 5.7 + 18.98 / pow(E / 1e18, 0.5);
  return res * 0.01;
}

double Hybrids_ICRC_2019_SD750(const double lgE) {
  // From Alan for ICRC2019
  // const double pars[2] = {0.084, 0.047};
  // const double pars[2] = {0.1119, 0.0377};
  const double pars[2] = {0.0523, 0.0614};

  //const double pars[2] = {0.1164, 0.0408};  //Res-Up
  //const double pars[2] = {0.1073, 0.0346};  //Res-Down

  const double resolution = pars[0] + pars[1] / sqrt(pow(10, lgE) / 1.e18);
  return sqrt(pow(resolution, 2) + pow(0.05, 2));
}

double Hybrids_ICRC_2019_Offline_11_06(const double lgE) {
  /* Preliminary version for ICRC 2019 */
  // V. Verzi 11 06 2019
  // Resolution calculated with hybrids data
  // OFFLINE
  double E = pow(10, lgE);
  double res = 4.38 + 22.36 / pow(E / 1e18, 0.5);
  return res * 0.01;
}

double Hybrids_ICRC_2019_Offline_05_07(const double lgE) {
  /* Preliminary version for ICRC 2019 */
  // V. Verzi 05 07 2019
  // Resolution calculated with hybrids data
  // OFFLINE
  double E = pow(10, lgE);
  double res = 7.76 + 15.74 * exp(-0.152 * E/1e18);
  return res * 0.01;
}

double Hybrids_ICRC_2019_Herald_11_06(const double lgE) {
  /* Preliminary version for ICRC 2019 */
  // V. Verzi 11 06 2019
  // Resolution calculated with hybrids data
  // HERALD
  double E = pow(10, lgE);
  double res = 4.57 + 21.69 / pow(E / 1e18, 0.5);
  return res * 0.01;
}


double SimulationSD1500EnergyResolution(const double lgE, const double zenith) {
  // obtained from 50-50 p/Fe simulations using QGSJetII-03
  double p0 =  1.090032150314515363e-01;
  double p1 =  4.353891452953017049e-01;
  double Esd = pow(10., lgE);

  double relErrs = p0 + p1 / sqrt(Esd / 1e+17) ;

  return relErrs;
}

double SimulationSD750EnergyResolution_ICRC2017(const double lgE, const double zenith) {
  //new infill resolution model from Daniela
  //determined together with fd energy bias correction and new invisible energy model
  //QGSJET 2.4 simulations, quality + FoV cuts, resolution determined for p and Fe
  const double pars_p[2] = {1.078638367648934449e-01 ,   1.595403269466228735e-01};
  const double pars_fe[2] = { 4.836779040905066218e-02, 1.697359371453481258e-01};

  const double pars[2] = {(pars_p[0] + pars_fe[0]) / 2., (pars_p[1] + pars_fe[1]) / 2.};

  return pars[0] + pars[1] / sqrt(pow(10, lgE) / 1e17);
}

double DataSD1500EnergyResolution(const double lgE, const double zenith, TVectorD pCal, double lgS38) {
  // obtained from data itself
  const double pars[9] = {7.964648011936218460e-04,
                          4.118659698340852438e-01,
                          -1.517612358021901198e+00,
                          9.102368706313619384e+00,
                          -9.938644262370022631e+00,
                          4.120050818744602217e-01,
                          1.244124736151181560e-03,
                          -6.535475475395227107e+00,
                          1.039836139398396853e+01
                         };
  const double lgsectheta = TMath::Log10(1 / cos(zenith));
  double res = pars[0] + (pars[1] + pars[2] * lgsectheta + pars[3] * lgsectheta * lgsectheta + pars[4] * lgsectheta * lgsectheta * lgsectheta) * pow(10., -lgS38 * 0.5)
               + (pars[5] + pars[6] * lgsectheta + pars[7] * lgsectheta * lgsectheta + pars[8] * lgsectheta * lgsectheta * lgsectheta) * pow(10., -lgS38);

  return res * pCal[1];
}

double DataSD750EnergyResolution_ICRC2017(const double lgE, const double zenith, TVectorD pCal, double lgS35) {
  const double pars[6] = {1.713893029281758462e+00,
                          -1.115258876619560580e+00,
                          1.557604707934408061e+01,
                          6.861112716697782554e-01,
                          -8.015754473911778089e-01,
                          8.852556156454199909e-01
                         };
  const double lgsectheta = TMath::Log10(1 / cos(zenith));

  return (pars[0] + pars[1] * lgsectheta + pars[2] * lgsectheta * lgsectheta)
         * (pow(10., (-pars[3] * lgS35)) + pars[4] * pow(10., (-pars[5] * lgS35)));
}



double Valerio12_11_2018(const double lgE) {
  // Comments
  // Test from Valerio
  // Resolution from Hybrids  12/11/2018

  const double pars[3] = {2.46619e1, -1.69615e1, 3.93182};
  double res = pars[0] + ( lgE - 18. ) * pars[1] + pow( lgE - 18., 2) * pars[2];

  return res * 0.01;
}

double Resolution_TestExample(const double lgE) {
  /* Test for resolution */
  //test function. Add here any function  you want to test
  // input is logE
  // output is resolution

  double res = 0.;
  return res;
}