FitterUtil.h

Go to the documentation of this file.

#ifndef FITTERUTIL_H
#define FITTERUTIL_H

#include <utl/PhysicalConstants.h>
#include <utl/AugerUnits.h>

namespace FitterUtil {

  const double gpercm2 = utl::gram / utl::cm2;

  const double tracebinning = 25. * utl::nanosecond;
  const double tracelength = 768;

  // numbering scheme for the signal components in the universality parametrisation (UnivParam)
  const int sc_mu    =  0;
  const int sc_eg    =  1;
  const int sc_egmu  =  2;
  const int sc_eghad =  3;
  const int sc_all   = -1;

  // numbers to specify reference of signal model (numbering scheme defined in UnivParam)
  const int up_qgsjet = 0;
  const int up_epos = 1;
  const int up_proton = 0;
  const int up_iron = 1;

  double getDistToFirstInteraction(double hfi, double sHeight, double zenith, double dist, double psi);
  double getFirstParticleArrivalTime(double r, double dx0);
  double getX0(int, double, double);

  double getMeanXmaxBias(double loge, double theta, bool saturated);
  double getMeanParametrizedXmax(double loge);

  double getMeanParametrizedXmaxData(double loge);

  double getMeanParametrizedNmuData(double loge, double theta);
  double getMeanParametrizedNmuData(double loge, double theta, double Xmax);

  double GetMeanMuonOffset(double r, double offset0);
  double GetMeanMuonOffsetMCCalibrated(double r, double lgE, double theta);
  double GetMeanMuonOffsetDataCalibrated(double r, double lgE, double theta);

}

#endif