MCSFitResult.h

Go to the documentation of this file.

#ifndef MCSFitResult_h
#define MCSFitResult_h

namespace recob {
  class MCSFitResult {
  public:
    MCSFitResult() = default;
    MCSFitResult(int pid,
                 float momFwd,
                 float momFwdUnc,
                 float llhdFwd,
                 float momBwd,
                 float momBwdUnc,
                 float llhdBwd,
                 const std::vector<float>& radlengths,
                 const std::vector<float>& angles)
      : pid_(pid)
      , momFwd_(momFwd)
      , momFwdUnc_(momFwdUnc)
      , llhdFwd_(llhdFwd)
      , momBwd_(momBwd)
      , momBwdUnc_(momBwdUnc)
      , llhdBwd_(llhdBwd)
      , radlengths_(radlengths)
      , angles_(angles)
    {}

    int particleIdHyp() const { return pid_; }

    float fwdMomentum() const { return momFwd_; }

    float fwdMomUncertainty() const { return momFwdUnc_; }

    float fwdLogLikelihood() const { return llhdFwd_; }

    float bwdMomentum() const { return momBwd_; }

    float bwdMomUncertainty() const { return momBwdUnc_; }

    float bwdLogLikelihood() const { return llhdBwd_; }

    const std::vector<float>& segmentRadLengths() const { return radlengths_; }

    const std::vector<float>& scatterAngles() const { return angles_; }

    bool isBestFwd() const { return llhdFwd_ < llhdBwd_; }

    float bestMomentum() const { return isBestFwd() ? momFwd_ : momBwd_; }

    float bestMomUncertainty() const { return isBestFwd() ? momFwdUnc_ : momBwdUnc_; }

    float bestLogLikelihood() const { return isBestFwd() ? llhdFwd_ : llhdBwd_; }

    float deltaLogLikelihood() const { return std::abs(llhdFwd_ - llhdBwd_); }
    //
  private:
    int pid_;         
    float momFwd_;    
    float momFwdUnc_; 
    float
      llhdFwd_; 
    float momBwd_;    
    float momBwdUnc_; 
    float
      llhdBwd_; 
    std::vector<float> radlengths_; 
    std::vector<float> angles_;     
  };
}

#endif