#include <MuonHoughPattern.h>

Inheritance diagram for MuonHoughPattern:

Collaboration diagram for MuonHoughPattern:

Public Types
using	MuonHoughPatternCollection = DataVector< MuonHoughPattern >

Public Member Functions
	MuonHoughPattern (int id_number)
	MuonHoughPattern does not own its hits (contrary to the default) MuonHoughHitContainer! More...

virtual	~MuonHoughPattern ()=default
	destructor More...

void	resetTracksegment ()
	clear pattern More...

bool	hitInHoughPattern (const std::shared_ptr< MuonHoughHit > &hit) const
	returns if hit is in pattern More...

double	calculateEZ () const
	calculate estimated z-position of pattern More...

double	patternLength () const
	returns distance between first and last hit More...

void	printHoughPattern () const
	prints out info about hough pattern More...

int	getIdNumber () const
	returns id number of hough transform used to generate pattern More...

double	getEPhi () const
	returns phi of pattern More...

double	getERPhi () const
	returns r0/d0 of pattern More...

double	getETheta () const
	returns theta of pattern More...

double	getERTheta () const
	returns z0 of pattern More...

double	getECurvature () const
	returns curvature of pattern More...

double	getMaximumHistogram () const
	returns maximum of histogram used to generate pattern More...

void	setEPhi (double ephi)
	set phi of pattern More...

void	setERPhi (double erphi)
	set r0 of pattern More...

void	setETheta (double etheta)
	set theta of pattern More...

void	setERTheta (double ertheta)
	set z0 of pattern More...

void	setECurvature (double curvature)
	set curvature of pattern More...

void	setMaximumHistogram (double maximumhistogram)
	set maximum of histogram used to generate pattern More...

double	getEAngle () const
	returns angle in precision plane in rad More...

double	getER () const
	returns radius in precision plane in mm More...

void	setEAngle (double eangle)
	set angle in precision plane in rad More...

void	setER (double er)
	set radius in precision plane in mm More...

void	setWhichSegment (bool which_segment)
	set which segment pattern is in, not in use More...

Amg::Vector3D	getEPos () const
	calulates 3d point closest to ip More...

Amg::Vector3D	getEDir () const
	calculates direction at point closest to ip More...

void	updateParametersRPhi (bool cosmics=false)
	update parameters in rphi plane based on weighted fit More...

std::shared_ptr< MuonHoughHit >	getHit (int hitno) const
	returns Hit at position hitno More...

const std::vector< std::shared_ptr< MuonHoughHit > > &	getHits () const
	returns hit vector More...

void	addHit (const std::shared_ptr< MuonHoughHit > &hit)
	add hit to container More...

void	removeHit (unsigned int hitno)
	remove hit from container More...

int	getHitId (unsigned int hitno) const
	returns hitid of hit hitno More...

double	getHitx (unsigned int hitno) const
	returns x position of hit hitno More...

double	getHity (unsigned int hitno) const
	returns y position of hit hitno More...

double	getHitz (unsigned int hitno) const
	returns z position of hit hitno More...

double	getRadius (unsigned int hitno) const
	returns radius of hit hitno More...

double	getTheta (unsigned int hitno) const
	returns theta of hit hitno More...

double	getPhi (unsigned int hitno) const
	returns phi of hit hitno More...

double	getWeight (unsigned int hitno) const
	returns weight of hit hitno More...

double	getOrigWeight (unsigned int hitno) const
	returns the orignal weight of hit hitno More...

bool	getMeasuresPhi (unsigned int hitno) const
	returns if hit hitno measures phi More...

const Trk::PrepRawData *	getPrd (unsigned int hitno) const
	returns preprawdata pointer of hit hitno More...

std::string	getWhichDetector (unsigned int hitno) const
	returns detectortechnology in string of hit hitno More...

MuonHough::DetectorTechnology	getDetectorId (unsigned int hitno) const
	returns detectortechnology of hit hitno More...

unsigned int	size () const
	returns size of hitcontainer More...

bool	empty () const
	returns if hitcontainer is empty More...

void	reserve (int size)
	allocates memory for hitvector More...

int	getRPChitno () const
	returns number of rpc hits in container More...

int	getRPCetahitno () const
	returns number of rpc eta hits in container More...

int	getMDThitno () const
	returns number of mdt hits in container More...

int	getCSChitno () const
	returns number of csc hits in container More...

int	getTGChitno () const
	returns number of tgc hits in container More...

Protected Attributes
std::vector< std::shared_ptr< MuonHoughHit > >	m_hit
	vector of hits in container More...

Private Attributes
int	m_id_number {-1}
	id number of hough transform used to generate pattern More...

bool	m_whichsegment {false}
	which segment is pattern created in, not relevant if split search is off 0 lower segment, 1 uppersegment More...

double	m_ephi {-M_PI_2}
	phi in rad More...

double	m_erphi {0.}
	r0 in mm More...

double	m_etheta {M_PI_2}
	theta in rad More...

double	m_ertheta {0.}
	z0 in mm More...

double	m_ecurvature {1.}
	curvature of pattern in rz plane in mm More...

double	m_maximumhistogram {0.}
	maximum of histogram More...

Detailed Description

Definition at line 17 of file MuonHoughPattern.h.

Member Typedef Documentation

◆ MuonHoughPatternCollection

using MuonHoughPattern::MuonHoughPatternCollection = DataVector<MuonHoughPattern>

Definition at line 90 of file MuonHoughPattern.h.

Constructor & Destructor Documentation

◆ MuonHoughPattern()

MuonHoughPattern::MuonHoughPattern ( int id_number )

MuonHoughPattern does not own its hits (contrary to the default) MuonHoughHitContainer!

constructor

Definition at line 11 of file MuonHoughPattern.cxx.

11 : MuonHoughHitContainer(), m_id_number(id_number) {}

◆ ~MuonHoughPattern()

virtual MuonHoughPattern::~MuonHoughPattern ( )

virtualdefault

destructor

Member Function Documentation

◆ addHit()

void MuonHoughHitContainer::addHit ( const std::shared_ptr< MuonHoughHit > & hit )

inherited

add hit to container

Definition at line 8 of file MuonHoughHitContainer.cxx.

                                                                          {
     if (hit->getId() == -1) { hit->setId(size()); }
     m_hit.push_back(hit);
 }

◆ calculateEZ()

double MuonHoughPattern::calculateEZ ( ) const

calculate estimated z-position of pattern

Definition at line 65 of file MuonHoughPattern.cxx.

                                            {
     double ez = 0;
     double count = 0;
     //  double r0 = m_erphi;
  
     for (unsigned int hitno = 0; hitno < m_hit.size(); hitno++) {
         // for each hit the distance from the hit to x,y line through (0,0) in the xy plane is calcualted and then distance*cos(theta) gives
         // z_hit - z
  
         double distance = MuonHoughMathUtils::signedDistanceToLine(getHit(hitno)->getHitx(), getHit(hitno)->getHity(), 0, m_ephi - M_PI);
  
         double z = getHit(hitno)->getHitz() -
                    (distance * std::cos(m_etheta) / std::sin(m_etheta));  // distance * sin (m_etheta) = L  // - sign correct?
         // hough correction in here?
  
         // straight line fitter?
  
         ez += z;
         count++;
     }
  
     ez = ez / count;
  
     return ez;
 }

◆ empty()

bool MuonHoughHitContainer::empty ( ) const

inlineinherited

returns if hitcontainer is empty

Definition at line 105 of file MuonHoughHitContainer.h.

105 { return m_hit.empty(); }

◆ getCSChitno()

int MuonHoughHitContainer::getCSChitno ( ) const

inherited

returns number of csc hits in container

Definition at line 37 of file MuonHoughHitContainer.cxx.

                                              {
     int cschitno = 0;
     for (unsigned int i = 0; i < m_hit.size(); i++) { cschitno += getDetectorId(i) == MuonHough::CSC; }
     return cschitno;
 }

◆ getDetectorId()

MuonHough::DetectorTechnology MuonHoughHitContainer::getDetectorId ( unsigned int hitno ) const

inlineinherited

returns detectortechnology of hit hitno

Definition at line 107 of file MuonHoughHitContainer.h.

107 { return m_hit[hitno]->getDetectorId(); }

◆ getEAngle()

double MuonHoughPattern::getEAngle ( ) const

returns angle in precision plane in rad

Definition at line 105 of file MuonHoughPattern.cxx.

                                          {
     double eangle = 0;
     switch (m_id_number) {
         case MuonHough::hough_xy: eangle = m_ephi; break;
         case MuonHough::hough_rz:
         case MuonHough::hough_rzcosmics:
         case MuonHough::hough_rz_rpc:
         case MuonHough::hough_rz_mdt:
         case MuonHough::hough_curved_at_a_cylinder: eangle = m_etheta; break;
         default:
             MsgStream log(Athena::getMessageSvc(), "MuonHoughPattern::getEAngle");
             if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "MuonHoughPattern::no valid id_number" << endmsg;
     }
     return eangle;
 }

◆ getECurvature()

double MuonHoughPattern::getECurvature ( ) const

inline

returns curvature of pattern

Definition at line 120 of file MuonHoughPattern.h.

120 { return m_ecurvature; }

◆ getEDir()

Amg::Vector3D MuonHoughPattern::getEDir ( ) const

calculates direction at point closest to ip

Definition at line 185 of file MuonHoughPattern.cxx.

                                             {
     CxxUtils::sincos scphi(m_ephi);
     CxxUtils::sincos sctheta(m_etheta);
     return {scphi.cs * sctheta.sn, scphi.sn * sctheta.sn, sctheta.cs};
 }

◆ getEPhi()

double MuonHoughPattern::getEPhi ( ) const

inline

returns phi of pattern

Definition at line 116 of file MuonHoughPattern.h.

116 { return m_ephi; }

◆ getEPos()

Amg::Vector3D MuonHoughPattern::getEPos ( ) const

calulates 3d point closest to ip

Definition at line 165 of file MuonHoughPattern.cxx.

                                             {
     // similar to StandardAlgs::shortestPointOfLineToOrigin3D
  
     // should be maybe really shortest point, but
  
     // problem is that there is no starting point on the line (well this calc. pos could be available)
  
     Amg::Vector3D pos;
  
     pos[0] = m_erphi * std::sin(m_ephi);
     pos[1] = -m_erphi * std::cos(m_ephi);
  
     if (m_etheta != 0) {
         pos[2] = m_ertheta / std::sin(m_etheta);
     } else {
         pos[2] = 0;
     }
     return pos;
 }

◆ getER()

double MuonHoughPattern::getER ( ) const

returns radius in precision plane in mm

Definition at line 121 of file MuonHoughPattern.cxx.

                                      {
     double er = 0;
     switch (m_id_number) {
         case MuonHough::hough_xy: er = m_erphi; break;
         case MuonHough::hough_rz:
         case MuonHough::hough_rzcosmics:
         case MuonHough::hough_rz_rpc:
         case MuonHough::hough_rz_mdt:
         case MuonHough::hough_curved_at_a_cylinder: er = m_ertheta; break;
         default:
             MsgStream log(Athena::getMessageSvc(), "MuonHoughPattern::getER");
             if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "MuonHoughPattern::no valid id_number" << endmsg;
     }
     return er;
 }

◆ getERPhi()

double MuonHoughPattern::getERPhi ( ) const

inline

returns r0/d0 of pattern

Definition at line 117 of file MuonHoughPattern.h.

117 { return m_erphi; }

◆ getERTheta()

double MuonHoughPattern::getERTheta ( ) const

inline

returns z0 of pattern

Definition at line 119 of file MuonHoughPattern.h.

119 { return m_ertheta; }

◆ getETheta()

double MuonHoughPattern::getETheta ( ) const

inline

returns theta of pattern

Definition at line 118 of file MuonHoughPattern.h.

118 { return m_etheta; }

◆ getHit()

std::shared_ptr< MuonHoughHit > MuonHoughHitContainer::getHit ( int hitno ) const

inlineinherited

returns Hit at position hitno

Definition at line 91 of file MuonHoughHitContainer.h.

91 { return m_hit.at(hitno); }

◆ getHitId()

int MuonHoughHitContainer::getHitId ( unsigned int hitno ) const

inlineinherited

returns hitid of hit hitno

Definition at line 93 of file MuonHoughHitContainer.h.

93 { return m_hit[hitno]->getId(); }

◆ getHits()

const std::vector< std::shared_ptr< MuonHoughHit > > & MuonHoughHitContainer::getHits ( ) const

inlineinherited

returns hit vector

Definition at line 92 of file MuonHoughHitContainer.h.

92 { return m_hit; }

◆ getHitx()

double MuonHoughHitContainer::getHitx ( unsigned int hitno ) const

inlineinherited

returns x position of hit hitno

Definition at line 94 of file MuonHoughHitContainer.h.

94 { return m_hit[hitno]->getHitx(); }

◆ getHity()

double MuonHoughHitContainer::getHity ( unsigned int hitno ) const

inlineinherited

returns y position of hit hitno

Definition at line 95 of file MuonHoughHitContainer.h.

95 { return m_hit[hitno]->getHity(); }

◆ getHitz()

double MuonHoughHitContainer::getHitz ( unsigned int hitno ) const

inlineinherited

returns z position of hit hitno

Definition at line 96 of file MuonHoughHitContainer.h.

96 { return m_hit[hitno]->getHitz(); }

◆ getIdNumber()

int MuonHoughPattern::getIdNumber ( ) const

inline

returns id number of hough transform used to generate pattern

Definition at line 114 of file MuonHoughPattern.h.

114 { return m_id_number; }

◆ getMaximumHistogram()

double MuonHoughPattern::getMaximumHistogram ( ) const

inline

returns maximum of histogram used to generate pattern

Definition at line 121 of file MuonHoughPattern.h.

121 { return m_maximumhistogram; }

◆ getMDThitno()

int MuonHoughHitContainer::getMDThitno ( ) const

inherited

returns number of mdt hits in container

Definition at line 25 of file MuonHoughHitContainer.cxx.

                                              {
     int mdthitno = 0;
     for (unsigned int i = 0; i < m_hit.size(); i++) { mdthitno += getDetectorId(i) == MuonHough::MDT; }
     return mdthitno;
 }

◆ getMeasuresPhi()

bool MuonHoughHitContainer::getMeasuresPhi ( unsigned int hitno ) const

inlineinherited

returns if hit hitno measures phi

Definition at line 109 of file MuonHoughHitContainer.h.

109 { return m_hit[hitno]->getMeasuresPhi(); }

◆ getOrigWeight()

double MuonHoughHitContainer::getOrigWeight ( unsigned int hitno ) const

inlineinherited

returns the orignal weight of hit hitno

Definition at line 102 of file MuonHoughHitContainer.h.

102 { return m_hit[hitno]->getOrigWeight(); }

◆ getPhi()

double MuonHoughHitContainer::getPhi ( unsigned int hitno ) const

inlineinherited

returns phi of hit hitno

Definition at line 100 of file MuonHoughHitContainer.h.

100 { return m_hit[hitno]->getPhi(); }

◆ getPrd()

const Trk::PrepRawData * MuonHoughHitContainer::getPrd ( unsigned int hitno ) const

inlineinherited

returns preprawdata pointer of hit hitno

Definition at line 110 of file MuonHoughHitContainer.h.

110 { return m_hit[hitno]->getPrd(); }

◆ getRadius()

double MuonHoughHitContainer::getRadius ( unsigned int hitno ) const

inlineinherited

returns radius of hit hitno

Definition at line 97 of file MuonHoughHitContainer.h.

97 { return m_hit[hitno]->getRadius(); }

◆ getRPCetahitno()

int MuonHoughHitContainer::getRPCetahitno ( ) const

inherited

returns number of rpc eta hits in container

Definition at line 31 of file MuonHoughHitContainer.cxx.

                                                 {
     int rpchitno = 0;
     for (unsigned int i = 0; i < m_hit.size(); i++) { rpchitno += getDetectorId(i) == MuonHough::RPC && !getMeasuresPhi(i); }
     return rpchitno;
 }

◆ getRPChitno()

int MuonHoughHitContainer::getRPChitno ( ) const

inherited

returns number of rpc hits in container

Definition at line 19 of file MuonHoughHitContainer.cxx.

                                              {
     int rpchitno = 0;
     for (unsigned int i = 0; i < m_hit.size(); i++) { rpchitno += getDetectorId(i) == MuonHough::RPC; }
     return rpchitno;
 }

◆ getTGChitno()

int MuonHoughHitContainer::getTGChitno ( ) const

inherited

returns number of tgc hits in container

Definition at line 43 of file MuonHoughHitContainer.cxx.

                                              {
     int tgchitno = 0;
     for (unsigned int i = 0; i < m_hit.size(); i++) { tgchitno += getDetectorId(i) == MuonHough::TGC; }
     return tgchitno;
 }

◆ getTheta()

double MuonHoughHitContainer::getTheta ( unsigned int hitno ) const

inlineinherited

returns theta of hit hitno

Definition at line 99 of file MuonHoughHitContainer.h.

99 { return m_hit[hitno]->getTheta(); }

◆ getWeight()

double MuonHoughHitContainer::getWeight ( unsigned int hitno ) const

inlineinherited

returns weight of hit hitno

Definition at line 101 of file MuonHoughHitContainer.h.

101 { return m_hit[hitno]->getWeight(); }

◆ getWhichDetector()

std::string MuonHoughHitContainer::getWhichDetector ( unsigned int hitno ) const

inlineinherited

returns detectortechnology in string of hit hitno

Definition at line 108 of file MuonHoughHitContainer.h.

108 { return m_hit[hitno]->getWhichDetector(); }

◆ hitInHoughPattern()

bool MuonHoughPattern::hitInHoughPattern ( const std::shared_ptr< MuonHoughHit > & hit ) const

returns if hit is in pattern

Definition at line 17 of file MuonHoughPattern.cxx.

 {
     for (unsigned int i = 0; i < size(); i++) {
         if (m_hit[i]->getId() == hit->getId()) { return true; }
     }
  
     return false;
 }

◆ patternLength()

double MuonHoughPattern::patternLength ( ) const

returns distance between first and last hit

Definition at line 26 of file MuonHoughPattern.cxx.

                                              {
     // takes the first 2 hits and calculates distance and then takes next hit, and calculates from previous 2 hits which 2 are farthest
     // away, etc.. also possible to calculate point closest to IP and determine if left/ right to pattern, app. just as fast.
  
     double max_patternlength = 0;
  
     if (m_hit.size() >= 2) {
         double pattern_length1 = 0, pattern_length2 = 0;
         int hitno1 = 0;
         int hitno2 = 1;
  
         Amg::Vector3D diff = m_hit[hitno1]->getPosition() - m_hit[hitno2]->getPosition();
         max_patternlength = diff.mag();
  
         for (unsigned int i = 2; i < m_hit.size(); i++) {
             diff = m_hit[hitno1]->getPosition() - m_hit[i]->getPosition();
             pattern_length1 = diff.mag();
  
             diff = m_hit[hitno2]->getPosition() - m_hit[i]->getPosition();
             pattern_length2 = diff.mag();
  
             if (pattern_length1 <= max_patternlength && pattern_length2 <= max_patternlength) {
                 // nothing happens..
             } else if (pattern_length1 > max_patternlength && pattern_length1 >= pattern_length2) {
                 hitno2 = i;
                 max_patternlength = pattern_length1;
             } else if (pattern_length2 > max_patternlength && pattern_length2 > pattern_length1) {
                 hitno1 = i;
                 max_patternlength = pattern_length2;
             }
         }
     } else {
         MsgStream log(Athena::getMessageSvc(), "MuonHoughPattern::patternLength");
         if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "MuonHoughPattern::pattern_size <2" << endmsg;
     }
  
     return max_patternlength;
 }

◆ printHoughPattern()

void MuonHoughPattern::printHoughPattern ( ) const

prints out info about hough pattern

Definition at line 91 of file MuonHoughPattern.cxx.

                                                {
     MsgStream log(Athena::getMessageSvc(), "MuonHoughPattern::printHoughPattern");
     if (m_hit.empty()) {
         if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "MuonHoughPattern::No pattern_Segment" << endmsg;
     }
  
     if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "MuonHoughPattern::Size of MuonHoughPattern: " << m_hit.size() << endmsg;
     for (unsigned int i = 0; i < m_hit.size(); i++) {
         if (log.level() <= MSG::VERBOSE)
             log << MSG::VERBOSE << m_hit[i]->getHitx() << " " << m_hit[i]->getHity() << " " << m_hit[i]->getHitz() << " "
                 << m_hit[i]->getId() << endmsg;
     }
 }

◆ removeHit()

void MuonHoughHitContainer::removeHit ( unsigned int hitno )

inherited

remove hit from container

Definition at line 13 of file MuonHoughHitContainer.cxx.

                                                         {
     if (hitno >= m_hit.size()) { throw std::runtime_error("MuonHoughHitContainer::range error!"); }
     
     m_hit.erase(m_hit.begin() + hitno);
 }

◆ reserve()

void MuonHoughHitContainer::reserve ( int size )

inlineinherited

allocates memory for hitvector

Definition at line 106 of file MuonHoughHitContainer.h.

106 { m_hit.reserve(size); }

◆ resetTracksegment()

void MuonHoughPattern::resetTracksegment ( )

clear pattern

Definition at line 13 of file MuonHoughPattern.cxx.

                                          {
     m_hit.clear();
 }

◆ setEAngle()

void MuonHoughPattern::setEAngle ( double eangle )

set angle in precision plane in rad

Definition at line 137 of file MuonHoughPattern.cxx.

                                               {
     switch (m_id_number) {
         case MuonHough::hough_xy: m_ephi = eangle; break;
         case MuonHough::hough_rz:
         case MuonHough::hough_rzcosmics:
         case MuonHough::hough_rz_rpc:
         case MuonHough::hough_rz_mdt:
         case MuonHough::hough_curved_at_a_cylinder: m_etheta = eangle; break;
         default:
             MsgStream log(Athena::getMessageSvc(), "MuonHoughPattern::setEAngle");
             if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "MuonHoughPattern::no valid id_number" << endmsg;
     }
 }

◆ setECurvature()

void MuonHoughPattern::setECurvature ( double curvature )

inline

set curvature of pattern

Definition at line 127 of file MuonHoughPattern.h.

127 { m_ecurvature = ecurvature; }

◆ setEPhi()

void MuonHoughPattern::setEPhi ( double ephi )

inline

set phi of pattern

Definition at line 123 of file MuonHoughPattern.h.

123 { m_ephi = ephi; }

◆ setER()

void MuonHoughPattern::setER ( double er )

set radius in precision plane in mm

Definition at line 151 of file MuonHoughPattern.cxx.

                                       {
     switch (m_id_number) {
         case MuonHough::hough_xy: m_erphi = er; break;
         case MuonHough::hough_rz:
         case MuonHough::hough_rzcosmics:
         case MuonHough::hough_rz_rpc:
         case MuonHough::hough_rz_mdt:
         case MuonHough::hough_curved_at_a_cylinder: m_ertheta = er; break;
         default:
             MsgStream log(Athena::getMessageSvc(), "MuonHoughPattern::setER");
             if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "MuonHoughPattern::no valid id_number" << endmsg;
     }
 }

◆ setERPhi()

void MuonHoughPattern::setERPhi ( double erphi )

inline

set r0 of pattern

Definition at line 124 of file MuonHoughPattern.h.

124 { m_erphi = erphi; }

◆ setERTheta()

void MuonHoughPattern::setERTheta ( double ertheta )

inline

set z0 of pattern

Definition at line 126 of file MuonHoughPattern.h.

126 { m_ertheta = ertheta; }

◆ setETheta()

void MuonHoughPattern::setETheta ( double etheta )

inline

set theta of pattern

Definition at line 125 of file MuonHoughPattern.h.

125 { m_etheta = etheta; }

◆ setMaximumHistogram()

void MuonHoughPattern::setMaximumHistogram ( double maximumhistogram )

inline

set maximum of histogram used to generate pattern

Definition at line 128 of file MuonHoughPattern.h.

128 { m_maximumhistogram = maximumhistogram; }

◆ setWhichSegment()

void MuonHoughPattern::setWhichSegment ( bool which_segment )

inline

set which segment pattern is in, not in use

Definition at line 130 of file MuonHoughPattern.h.

130 { m_whichsegment = which_segment; }

◆ size()

unsigned int MuonHoughHitContainer::size ( ) const

inlineinherited

returns size of hitcontainer

Definition at line 104 of file MuonHoughHitContainer.h.

104 { return m_hit.size(); }

◆ updateParametersRPhi()

void MuonHoughPattern::updateParametersRPhi ( bool cosmics = false )

update parameters in rphi plane based on weighted fit

Definition at line 192 of file MuonHoughPattern.cxx.

                                                         {
     if (empty()) return;
  
     double sum_x = 0.;
     double sum_y = 0.;
  
     int hitsno = size();
  
     for (unsigned int i = 0; i < size(); i++) {
         sum_x += getHitx(i);  //*getOrigWeight(i);
         sum_y += getHity(i);  //*getOrigWeight(i);
     }
  
     // adding interaction point constraint:
     if (!cosmics || size() == 1) hitsno += size();
  
     const double av_x = sum_x / (hitsno + 0.);
     const double av_y = sum_y / (hitsno + 0.);
  
     double sumx = 0.;
     double sumy = 0.;
     for (unsigned int i = 0; i < size(); i++) {
         double hitx = getHitx(i);
         double hity = getHity(i);
         double x_offset = hitx - av_x;
         double y_offset = hity - av_y;
         double weight = (x_offset * x_offset + y_offset * y_offset);  //*getOrigWeight(i);
         int sign = 1;
         if (!cosmics) {
             if (x_offset * hitx + y_offset * hity < 0) { sign = -1; }
         } else {
             if (y_offset < 0) { sign = -1; }  // assume cosmic muon comes from above
         }
         sumx += weight * sign * x_offset;
         sumy += weight * sign * y_offset;
     }
     // adding interaction point (count size)
     if (!cosmics) {
         double weight = av_x * av_x + av_y * av_y;
         sumx += size() * weight * av_x;
         sumy += size() * weight * av_y;
     }
  
     if (std::abs(sumx) < 0.000001 || std::abs(sumy) < 0.000001) { return; }
  
     double phi = std::atan2(sumy, sumx);
     if (cosmics) {
         if (phi > 0) phi -= M_PI;  // phi between 0,-Pi for cosmics!
     }
  
     CxxUtils::sincos scphi(phi);
  
     const double r0 = scphi.apply(av_x, -av_y);
  
     setEPhi(phi);
     setERPhi(r0);
 }

Member Data Documentation

◆ m_ecurvature

double MuonHoughPattern::m_ecurvature {1.}

private

curvature of pattern in rz plane in mm

Definition at line 109 of file MuonHoughPattern.h.

◆ m_ephi

double MuonHoughPattern::m_ephi {-M_PI_2}

private

phi in rad

Definition at line 101 of file MuonHoughPattern.h.

◆ m_erphi

double MuonHoughPattern::m_erphi {0.}

private

r0 in mm

Definition at line 103 of file MuonHoughPattern.h.

◆ m_ertheta

double MuonHoughPattern::m_ertheta {0.}

private

z0 in mm

Definition at line 107 of file MuonHoughPattern.h.

◆ m_etheta

double MuonHoughPattern::m_etheta {M_PI_2}

private

theta in rad

Definition at line 105 of file MuonHoughPattern.h.

◆ m_hit

std::vector<std::shared_ptr<MuonHoughHit> > MuonHoughHitContainer::m_hit

protectedinherited

vector of hits in container

Definition at line 86 of file MuonHoughHitContainer.h.

◆ m_id_number

int MuonHoughPattern::m_id_number {-1}

private

id number of hough transform used to generate pattern

Definition at line 94 of file MuonHoughPattern.h.

◆ m_maximumhistogram

double MuonHoughPattern::m_maximumhistogram {0.}

private

maximum of histogram

Definition at line 111 of file MuonHoughPattern.h.

◆ m_whichsegment

bool MuonHoughPattern::m_whichsegment {false}

private

which segment is pattern created in, not relevant if split search is off 0 lower segment, 1 uppersegment

Definition at line 98 of file MuonHoughPattern.h.

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

Public Types

Public Member Functions

Protected Attributes

Private Attributes

Detailed Description

Member Typedef Documentation

◆ MuonHoughPatternCollection

Constructor & Destructor Documentation

◆ MuonHoughPattern()

◆ ~MuonHoughPattern()

Member Function Documentation

◆ addHit()

◆ calculateEZ()

◆ empty()

◆ getCSChitno()

◆ getDetectorId()

◆ getEAngle()

◆ getECurvature()

◆ getEDir()

◆ getEPhi()

◆ getEPos()

◆ getER()

◆ getERPhi()

◆ getERTheta()

◆ getETheta()

◆ getHit()

◆ getHitId()

◆ getHits()

◆ getHitx()

◆ getHity()

◆ getHitz()

◆ getIdNumber()

◆ getMaximumHistogram()

◆ getMDThitno()

◆ getMeasuresPhi()

◆ getOrigWeight()

◆ getPhi()

◆ getPrd()

◆ getRadius()

◆ getRPCetahitno()

◆ getRPChitno()

◆ getTGChitno()

◆ getTheta()

◆ getWeight()

◆ getWhichDetector()

◆ hitInHoughPattern()

◆ patternLength()

◆ printHoughPattern()

◆ removeHit()

◆ reserve()

◆ resetTracksegment()

◆ setEAngle()

◆ setECurvature()

◆ setEPhi()

◆ setER()

◆ setERPhi()

◆ setERTheta()

◆ setETheta()

◆ setMaximumHistogram()

◆ setWhichSegment()

◆ size()

◆ updateParametersRPhi()

Member Data Documentation

◆ m_ecurvature

◆ m_ephi

◆ m_erphi

◆ m_ertheta

◆ m_etheta

◆ m_hit

◆ m_id_number

◆ m_maximumhistogram

◆ m_whichsegment