upmuana::LlrUpmu Class Reference

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-16/Eval/LlrUpmu.h"

Public Member Functions
	LlrUpmu ()
	~LlrUpmu ()
double	getErr (double PE)
double	getDist (double x1, double y1, double z1, double x2, double y2, double z2)
void	setLLR (std::set< std::pair< rb::CellHit, rb::RecoHit >, bool(*)(std::pair< rb::CellHit, rb::RecoHit >, std::pair< rb::CellHit, rb::RecoHit >)>)
double	P_up () const
double	P_dn () const
double	chi2 () const
double	slope () const
double	llr () const
std::vector< rb::RecoHit >	outliers () const

Private Member Functions
void	LLR (std::vector< double > &eT, std::vector< double > &mT, std::vector< double > &mTErr, double &slope, double &chi2, double &P_up, double &P_dn, std::vector< rb::RecoHit > &in, std::vector< rb::RecoHit > &outliers)

Private Attributes
std::vector< rb::RecoHit >	outliers_
double	P_up_
double	P_dn_
double	chi2_
double	slope_
double	llr_

Detailed Description

Definition at line 17 of file LlrUpmu.h.

Constructor & Destructor Documentation

upmuana::LlrUpmu::LlrUpmu

(

)

Definition at line 11 of file LlrUpmu.cxx.

    :
    outliers_(), P_up_(1e-30), P_dn_(1e-30), 
    chi2_(1e6), slope_(-1e6), llr_(0.)
  {}

upmuana::LlrUpmu::~LlrUpmu

(

)

Definition at line 17 of file LlrUpmu.cxx.

  {}

Member Function Documentation

double upmuana::LlrUpmu::chi2 ( ) const

inline

Definition at line 33 of file LlrUpmu.h.

References chi2_.

Referenced by setLLR().

{ return chi2_; }

double upmuana::LlrUpmu::getDist ( double  x1, double  y1, double  z1, double  x2, double  y2, double  z2  )

inline

Definition at line 25 of file LlrUpmu.cxx.

References std::sqrt().

Referenced by setLLR().

                                                         {
    return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)+(z1-z2)*(z1-z2));
  }

double upmuana::LlrUpmu::getErr ( double  PE )

inline

Definition at line 20 of file LlrUpmu.cxx.

References cet::pow().

Referenced by setLLR().

                                         {
    //return 231594.0/(2267.16+pow(PE,2.01011)) + 9.55689;
    return 165143/(1882.9+pow(PE,2.11447)) + 10.4321;
  }

double upmuana::LlrUpmu::llr ( ) const

inline

Definition at line 35 of file LlrUpmu.h.

References llr_.

{ return llr_; }

void upmuana::LlrUpmu::LLR ( std::vector< double > &  eT, std::vector< double > &  mT, std::vector< double > &  mTErr, double &  slope, double &  chi2, double &  P_up, double &  P_dn, std::vector< rb::RecoHit > &  in, std::vector< rb::RecoHit > &  outliers  )

private

Definition at line 90 of file LlrUpmu.cxx.

References a, b, e, stan::math::fabs(), stan::math::gamma_q(), MECModelEnuComparisons::i, util::LinFit(), getGoodRuns4SAM::n, cet::pow(), and std::sqrt().

Referenced by setLLR().

  {
    // eT - x, mT - y
    size_t n = mT.size();

    // y = a + bx
    double a, b;
    if (eT.size() >= 2)
      chi2 = util::LinFit(eT, mT, wgts, b, a);
    else {
      chi2 = 1e6;
      slope = -1e6;
      P_up = 1e-30;
      P_dn = 1e-30;
      return;
    }

    size_t totAllowOutlier = n/10;
    size_t nOutliers = 0;
    std::vector<double> x_filt, y_filt, w_filt;
    for (size_t i=0; i < n; i++) {
      double y_fit = a + b*eT.at(i);
      if ( fabs(mT.at(i) - y_fit) < 5*(1.0/sqrt(wgts.at(i)))
           || nOutliers>totAllowOutlier)
        {
          x_filt.push_back(eT.at(i));
          y_filt.push_back(mT.at(i));
          w_filt.push_back(wgts.at(i));
        }
      else {
        ++nOutliers;
        outliers.push_back(in[i]);
      }
    }


    if (x_filt.size() >= 2)
      chi2 = util::LinFit(x_filt, y_filt, w_filt, b, a);
    else {
      chi2 = 1e6;
      slope = -1e6;
      P_up = 1e-30;
      P_dn = 1e-30;
    }
    n = x_filt.size();
    if (n < 5) {
      slope = 0;
      chi2  = 999;
      P_up  = 1e-30;
      P_dn  = 1;
      return;
    }

    slope = b;
    chi2 /= (double)(n-2);

    double  one_over_ee = 0.0,
      x_over_ee   = 0.0,
      y_over_ee   = 0.0;

    for (size_t i=0; i<n; ++i) {
      one_over_ee += w_filt.at(i);
      x_over_ee   += x_filt.at(i)*w_filt.at(i);
      y_over_ee   += y_filt.at(i)*w_filt.at(i);
    }

    double up_inter = (y_over_ee-x_over_ee)/one_over_ee;
    double dn_inter = (y_over_ee+x_over_ee)/one_over_ee;

    double up_chi2 = 0.0, dn_chi2 = 0.0;
    for (size_t i=0; i<n; ++i) {
      double e = 1.0/sqrt(w_filt.at(i));
      double up_expected = up_inter + x_filt.at(i);
      double dn_expected = dn_inter - x_filt.at(i);
      up_chi2 += pow((y_filt.at(i)-up_expected) / e, 2.0);
      dn_chi2 += pow((y_filt.at(i)-dn_expected) / e, 2.0);
    }

    double prob_up = boost::math::gamma_q((double)(n-2)/2.0,up_chi2/2.0),
    prob_dn = boost::math::gamma_q((double)(n-2)/2.0,dn_chi2/2.0);

    if (prob_up < 1e-30) prob_up = 1e-30;
    if (prob_dn < 1e-30) prob_dn = 1e-30;

    P_up = prob_up;
    P_dn = prob_dn;

  }

std::vector<rb::RecoHit> upmuana::LlrUpmu::outliers ( ) const

inline

Definition at line 36 of file LlrUpmu.h.

References outliers_.

Referenced by setLLR().

{ return outliers_; }

double upmuana::LlrUpmu::P_dn ( ) const

inline

Definition at line 32 of file LlrUpmu.h.

References P_dn_.

Referenced by setLLR().

{ return P_dn_; }

double upmuana::LlrUpmu::P_up ( ) const

inline

Definition at line 31 of file LlrUpmu.h.

References P_up_.

Referenced by setLLR().

{ return P_up_; }

void upmuana::LlrUpmu::setLLR

(

std::set< std::pair< rb::CellHit, rb::RecoHit >, bool(*)(std::pair< rb::CellHit, rb::RecoHit >, std::pair< rb::CellHit, rb::RecoHit >)>

sortedHits

)

Definition at line 30 of file LlrUpmu.cxx.

References chi2(), chi2_, dist, e, getDist(), getErr(), in, LLR(), llr_, test_ParserArtEvents::log, outliers(), outliers_, P_dn(), P_dn_, P_up(), P_up_, slope(), and slope_.

Referenced by upmuana::UpMuProb::produce().

  {
    std::vector<rb::RecoHit> outliers;
    double P_up;
    double P_dn;
    double chi2;
    double slope;

    if (sortedHits.size() < 2) {
      P_up = 1e-30;
      P_dn = 1e-30;
      chi2 = 1e6;
      slope = -1e6;
      return;
    }

    std::vector<double> measuredTimes;
    std::vector<double> expectedTimes;
    std::vector<double> wgts;

    // first hit
    measuredTimes.push_back(0.0);
    expectedTimes.push_back(0.0);
    double err = getErr(sortedHits.begin()->first.PE());
    wgts.push_back(1.0/err/err);

    double  startX = sortedHits.begin()->second.X(),
      startY = sortedHits.begin()->second.Y(),
      startZ = sortedHits.begin()->second.Z(),
      startT = sortedHits.begin()->second.T();

    std::vector<rb::RecoHit> in;

    for (auto i_hit =  ++sortedHits.begin();
         i_hit != sortedHits.end();
         ++i_hit) {
      in.push_back(i_hit->second);
      measuredTimes.push_back(i_hit->second.T() - startT);
      double dist = getDist(i_hit->second.X(), i_hit->second.Y(),
                            i_hit->second.Z(), startX, startY, startZ);
      expectedTimes.push_back(dist/29.97);
      err = getErr(i_hit->first.PE());
      wgts.push_back(1.0/err/err);
    }

    LLR(expectedTimes, measuredTimes, wgts, slope, chi2, P_up, P_dn,
        in, outliers);

    outliers_ = outliers;
    P_up_ = P_up;
    P_dn_ = P_dn;
    chi2_ = chi2;
    slope_ = slope;
    llr_ = log(P_up/P_dn);

  }

double upmuana::LlrUpmu::slope ( ) const

inline

Definition at line 34 of file LlrUpmu.h.

References slope_.

Referenced by setLLR().

{ return slope_; }

Member Data Documentation

double upmuana::LlrUpmu::chi2_

private

Definition at line 43 of file LlrUpmu.h.

Referenced by chi2(), and setLLR().

double upmuana::LlrUpmu::llr_

private

Definition at line 45 of file LlrUpmu.h.

Referenced by llr(), and setLLR().

std::vector<rb::RecoHit> upmuana::LlrUpmu::outliers_

private

Definition at line 40 of file LlrUpmu.h.

Referenced by outliers(), and setLLR().

double upmuana::LlrUpmu::P_dn_

private

Definition at line 42 of file LlrUpmu.h.

Referenced by P_dn(), and setLLR().

double upmuana::LlrUpmu::P_up_

private

Definition at line 41 of file LlrUpmu.h.

Referenced by P_up(), and setLLR().

double upmuana::LlrUpmu::slope_

private

Definition at line 44 of file LlrUpmu.h.

Referenced by setLLR(), and slope().

---

The documentation for this class was generated from the following files:

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-16/Eval/LlrUpmu.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-16/Eval/LlrUpmu.cxx