ShowerSizeFunction.h

Go to the documentation of this file.

#ifndef _SdHorizontalReconstruction_ShowerSizeFunction_h_
#define _SdHorizontalReconstruction_ShowerSizeFunction_h_

#include "FitInterface.h"

#include <Minuit2/FCNBase.h>
#include <vector>

namespace SdHorizontalReconstructionNS {

  class SdHorizontalReconstruction;

  class ShowerSizeFunction : public ROOT::Minuit2::FCNBase {
  public:

    enum Type {
      eApprox, // ChiSquare minimization to get starting value for full -logLikelihood minimization
      eFull  // Full -logLikelihood minimization, silent/saturated stations included, less stable
    };

    ShowerSizeFunction(const SdHorizontalReconstruction& config,
                         const StationList& list,
                         const SilentStationList& slist,
                         const AxisData& ad,
                         const ExternalGeometryData& gd);

    ~ShowerSizeFunction() {}

    double
    operator()(const std::vector<double>& pars) const;

    double
    Up() const { return 0.5; }

    void
    Predict(double& meanSignal, double& sigmaSignal,
            const utl::Point& station,
            const utl::Point& origin,
            const double n19,
            const utl::Point& core,
            const double recoveryErr = 0) const;

    void
    SetType(const Type type) { fType = type; }

  protected:
    Type fType;
    const SdHorizontalReconstruction& fConfig;
    const StationList& fList;
    const SilentStationList& fSList;
    unsigned char fUseAxisCovariance;
    double fAxisPar[3];
    double fAxisInvCov[3][3];
    const ExternalGeometryData& fExt;
  };

} // NS SdHorizontalReconstructionNS

#endif