PRD2020/BiasFunctions.hh

Go to the documentation of this file.

double HybridBias_2019_SD1500(double ESD, double coszenith, double maxZenith) {
  //bias as calculated by V. Verzi. Based on hybrid events
  //ICRC 2019
  // Gap2019-022
  double biasFinal = 1;

  if (acos(coszenith) < maxZenith * TMath::DegToRad() && acos(coszenith) > 0.) {
    if (ESD < 18.4)
      biasFinal = 1. + (24.3 + 77.61 * exp(-10. * (coszenith   - 0.5 ) ) ) * (18.4 - ESD) * 0.01;
  }
  return biasFinal;
}

double HybridBias_2019_SD750(double lgE, double fdB) {
  //bias as calculated by A Coleman. Based on Monte Carlo
  //ICRC 2019
  // March Auger meeting 2019
  // https://indico.nucleares.unam.mx/event/1408/session/16/contribution/118/material/slides/0.pdf

  const double pars[4] = { -0.1167, 5.022, 16.08, 21.61}; //Auger Mixed Comp (2017)
  // const double pars[4] = {-0.4843, 8.333, 16.43, 16.94}; //Proton
  // const double pars[4] = {-0.8852, 1.807, 15.75, 17.37}; //Iron
  const double biasS35 = 1 + pars[0] * exp(-pars[1] * pow(lgE - pars[2], 2.)) * (lgE - pars[3]);

  return pow(biasS35, fdB);
}

double HybridBias_2019_11_06_Offline_SD1500(double ESD, double coszenith, double maxZenith) {
  // Bias as calculated by V. Verzi. Based on hybrid events
  // ICRC 2019
  // June 11 2019
  // OFFLINE
  double biasFinal = 1;

  if (acos(coszenith) < maxZenith * TMath::DegToRad() && acos(coszenith) > 0.) {
    if (ESD < 18.4)
      biasFinal = 1. + (20.01 + 59.2 * exp(-10. * (coszenith   - 0.5 ) ) ) * (18.4 - ESD) * 0.01;
  }
  return biasFinal;
}

double HybridBias_2019_11_06_Herald_SD1500(double ESD, double coszenith, double maxZenith) {
  // Bias as calculated by V. Verzi. Based on hybrid events
  // ICRC 2019
  // June 11 2019
  // HERALD
  double biasFinal = 1;

  if (acos(coszenith) < maxZenith * TMath::DegToRad() && acos(coszenith) > 0.) {
    if (ESD < 18.4)
      biasFinal = 1. + (36.08 + 9.27 * exp(-10. * (coszenith   - 0.5 ) ) ) * (18.4 - ESD) * 0.01;
  }
  return biasFinal;
}

double Bias_TestExample(double ESD, double coszenith) {
  /* Test for Bias */
  //test function. Add here any function  you want to test
  // input is lgE and cos(zenith)
  // output is bias

  double res = 1.;
  return res;
}