GlobalSim::ERatio Class Reference

#include <ERatio.h>

Collaboration diagram for GlobalSim::ERatio:

Public Member Functions
	ERatio (const std::string &name, double thresh=50., double maxERCut=2.5)
virtual	~ERatio ()=default
StatusCode	run (const LArStripNeighborhood &in, bool &result)
std::string	toString () const

Private Member Functions
double	secondaryPeakEnergy_forw (const StripDataVector &, std::size_t)
double	secondaryPeakEnergy_back (const StripDataVector &, std::size_t)
bool	eratioCut (double p_peakE, double s_peakE)

Private Attributes
std::string	m_name {"ERatio_object"}
double	m_minDeltaE {50.}
double	m_maxERCut {2.5}

Detailed Description

Definition at line 25 of file ERatio.h.

Constructor & Destructor Documentation

ERatio()

GlobalSim::ERatio::ERatio	(	const std::string &	name,
		double	thresh = 50.,
		double	maxERCut = 2.5
	)

Definition at line 12 of file ERatio.cxx.

                          :
    m_name{name}, m_minDeltaE{minDeltaE}, m_maxERCut{maxERCut} {
  }

~ERatio()

virtual GlobalSim::ERatio::~ERatio ( )

virtualdefault

Member Function Documentation

eratioCut()

bool GlobalSim::ERatio::eratioCut ( double  p_peakE, double  s_peakE  )

private

Definition at line 128 of file ERatio.cxx.

                                                       {
 
    if (p_peakE <= 0 or s_peakE <= 0) {return false;}
 
    return std::max(p_peakE, s_peakE)/std::min(p_peakE, s_peakE) < m_maxERCut;
 
  }

run()

StatusCode GlobalSim::ERatio::run	(	const LArStripNeighborhood &	in,
		bool &	result
	)

Definition at line 19 of file ERatio.cxx.

                                                           {
 
    result = false;
 
    std::size_t maxCellInd = nbhd.maxCellIndex();
    auto p_peakE = nbhd.phi_center().at(maxCellInd).m_e;
 
    for (const auto& ps : {nbhd.phi_center(),
               nbhd.phi_low(),
               nbhd.phi_high()})  {
      
      double s_peakE = secondaryPeakEnergy_forw(ps, maxCellInd);
 
      result = eratioCut(p_peakE, s_peakE);
      if (result) {break;}
      
      s_peakE = secondaryPeakEnergy_back(ps, maxCellInd);
   
      result = eratioCut(p_peakE, s_peakE);
 
      if (result) {break;} 
    }
 
    return StatusCode::SUCCESS;
  }

secondaryPeakEnergy_back()

double GlobalSim::ERatio::secondaryPeakEnergy_back ( const StripDataVector &  phi_str, std::size_t  maxCellInd  )

private

Definition at line 85 of file ERatio.cxx.

                                            {
 
    
    /*
     * Look for a secondary peak in the direction of decreasing eta
     * starting at the eta location of the maximum energy strip.
     */
 
    
    constexpr double maxStripE{std::numeric_limits<double>::max()};
    double peakE{0.};
 
    // eratio needs at least three elements
    auto sz = maxCellInd;
    if (sz < 2) {return peakE;}
 
    bool valley =
      phi_str.at(maxCellInd-1).m_e < maxStripE - m_minDeltaE;
 
 
    // calculate offset to point to 2 less than phi_str[maxCellInd]
    std::size_t start = phi_str.size() - maxCellInd + 1;
    for(auto iter = std::rbegin(phi_str)+start;
    iter != std::rend(phi_str);
    ++iter){
      
      if(!valley) {
    valley = iter->m_e < ((iter-1)->m_e) - m_minDeltaE;
    continue;
      }
 
      auto c_cellE = iter->m_e;
      auto p_cellE = (iter-1)->m_e;
      if(c_cellE > p_cellE + m_minDeltaE) {
    peakE = std::max(peakE, c_cellE);
      }
      
    }
    return peakE;
  }

secondaryPeakEnergy_forw()

double GlobalSim::ERatio::secondaryPeakEnergy_forw ( const StripDataVector &  phi_str, std::size_t  maxCellInd  )

private

Definition at line 46 of file ERatio.cxx.

                                            {
 
    /*
     * Look for a secondary peak in the direction of increasing eta
     * starting at the eta location of the maximum energy strip.
     */
    
    constexpr double maxStripE{std::numeric_limits<double>::max()};
    double peakE{0.};
 
    // eratio needs at least three elements
    auto sz = phi_str.size();
    if (sz < 3 or maxCellInd > sz) {return peakE;}
 
    bool valley =
      phi_str.at(maxCellInd+1).m_e < maxStripE - m_minDeltaE;
    
    for(auto iter = std::begin(phi_str)+maxCellInd+2;
    iter != std::end(phi_str);
    ++iter){
      
      if(!valley) {
    valley = iter->m_e < ((iter-1)->m_e) - m_minDeltaE;
      }
 
      if(valley) {
    auto c_cellE = iter->m_e;
    auto p_cellE = (iter-1)->m_e;
    if(c_cellE > p_cellE + m_minDeltaE) {
      peakE = std::max(peakE, c_cellE);
    }
      }
      
    }
    return peakE;
  }

toString()

std::string GlobalSim::ERatio::toString

(

)

const

Definition at line 136 of file ERatio.cxx.

                                  {
    std::stringstream ss;
    ss << "ERatio. "
       << " name: " << m_name
       << " minimum strip E difference to be significant " << m_minDeltaE
       << " upper ERatio limit for pi0 candidate " << m_maxERCut
       << '\n';
 
    return ss.str();
  }

Member Data Documentation

m_maxERCut

double GlobalSim::ERatio::m_maxERCut {2.5}

private

Definition at line 43 of file ERatio.h.

m_minDeltaE

double GlobalSim::ERatio::m_minDeltaE {50.}

private

Definition at line 42 of file ERatio.h.

m_name

std::string GlobalSim::ERatio::m_name {"ERatio_object"}

private

Definition at line 41 of file ERatio.h.

---

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

• ERatio.h
• ERatio.cxx