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