MuonHoughHit Class Reference

#include <MuonHoughHit.h>

Collaboration diagram for MuonHoughHit:

Public Member Functions
	MuonHoughHit (const Trk::PrepRawData *prd)
	default constructor (should not be used) More...
	MuonHoughHit (const Amg::Vector3D &pos_vec, bool measures_phi, MuonHough::DetectorTechnology detector_id, double prob=1, double weight=1., const Trk::PrepRawData *prd=nullptr, int id=-1)
	constructor More...
virtual	~MuonHoughHit ()=default
	destructor More...
MuonHough::DetectorTechnology	getDetectorId () const
	return DetectorTechnology More...
std::string	getWhichDetector () const
	return DetectorTechnology in string More...
const Amg::Vector3D &	getPosition () const
	return (x,y,z) vector More...
int	getId () const
	returns id More...
double	getHitx () const
	returns x position More...
double	getHity () const
	returns y position More...
double	getHitz () const
	returns z position More...
double	getWeight () const
	returns weight in histogram after rescaling More...
double	getOrigWeight () const
	returns original weight More...
double	getProbability () const
	returns probability that hit is part of pattern (true muon) More...
double	getRadius () const
	returns radius More...
double	getAbs () const
	returns radius More...
double	getTheta () const
	returns theta More...
double	getPhi () const
	returns phi More...
bool	isBarrel () const
	hit is barrel or endcap (for curved track model) More...
int	phiSector () const
	phi sector of hit More...
double	getMagneticTrackRatio () const
	ratio of the tracklength of the particle to which hit might belong would have traversed in magnetic field (for curved track model) More...
bool	getMeasuresPhi () const
	hit measures phi? More...
bool	getAssociated () const
	return if hit already associated to pattern More...
void	setId (int id)
	set id More...
void	setWeight (double weight)
	set weight More...
void	setProbability (double prob)
	set probability More...
void	setAssociated (bool associated)
	set associated More...
const Trk::PrepRawData *	getPrd () const
	return preprawdata More...

Protected Attributes
const Trk::PrepRawData *	m_prd {nullptr}
	Pointer to preprawdata. More...
Amg::Vector3D	m_pos {Amg::Vector3D::Zero()}
int	m_id {0}
	unique id More...
double	m_radius {0.}
	radius of hit More...
double	m_abs {0.}
	absolute value of hit More...
double	m_phi {0.}
	phi of hit More...
double	m_theta {0.}
	theta of hit More...
bool	m_barrel {true}
	hit is barrel / endcap More...
int	m_phi_sector {0}
	phi sector (0,15), Atlas Convention More...
double	m_magnetic_trackratio {0.}
	ratio of the tracklength of the particle to which hit might belong would have traversed in magnetic field (for curved track model) More...
double	m_weight {1.}
	weight of hit in histogram, after rescaling More...
const double	m_orig_weight {1.}
	original weight More...
double	m_probability {1.}
	probability that hit belongs to true muon (based on trigger confirmation and crude segment finding (used for mdts)) More...
MuonHough::DetectorTechnology	m_detector_id {MuonHough::MDT}
	detector technology of hit More...
bool	m_measures_phi {false}
	hit measures phi? More...
bool	m_associated {false}
	hit associated to pattern More...

Private Member Functions
int	calcPhiSector () const
	method that calculates phi sector of hit More...
double	calcMagneticTrackRatio () const
	method that calculates 'magnetic track ratio' (for curved track model) More...

Detailed Description

Definition at line 28 of file MuonHoughHit.h.

Constructor & Destructor Documentation

MuonHoughHit() [1/2]

MuonHoughHit::MuonHoughHit

(

const Trk::PrepRawData *

prd

)

default constructor (should not be used)

constructor per preprawdata (not used)

Definition at line 14 of file MuonHoughHit.cxx.

                                                    {
    const Trk::TrkDetElementBase* element = prd->detectorElement();
    const Identifier identifier = prd->identify();
 
    m_pos = element->surface(identifier).Trk::Surface::localToGlobal(prd->localPosition());
 
    m_radius = m_pos.perp();
    m_abs = m_pos.mag();
    m_theta = std::atan2(m_radius, getHitz());
    m_phi = m_pos.phi();
    m_barrel = (std::abs(m_radius / getHitz()) < MuonHough::tan_barrel) ? 0 : 1;
 
    m_phi_sector = calcPhiSector();
    m_magnetic_trackratio = calcMagneticTrackRatio();
    m_prd = prd;
}

MuonHoughHit() [2/2]

MuonHoughHit::MuonHoughHit	(	const Amg::Vector3D &	pos_vec,
		bool	measures_phi,
		MuonHough::DetectorTechnology	detector_id,
		double	prob = 1,
		double	weight = 1.,
		const Trk::PrepRawData *	prd = nullptr,
		int	id = -1
	)

constructor

Definition at line 31 of file MuonHoughHit.cxx.

                                                                             :
    m_prd(prd),
    m_pos(pos_vec),
    m_id(id),
    m_radius(m_pos.perp()),
    m_abs(m_pos.mag()),
    m_phi(m_pos.phi()),
    m_theta(std::atan2(m_radius, getHitz())),
    m_barrel ((std::abs(m_radius / getHitz()) < MuonHough::tan_barrel) ? 0 : 1),
    m_phi_sector(calcPhiSector()),
    m_magnetic_trackratio(calcMagneticTrackRatio()),
    m_weight(weight),
    m_orig_weight(weight),
    m_probability(prob),
    m_detector_id(detector_id),
    m_measures_phi(measures_phi)
{
}

~MuonHoughHit()

virtual MuonHoughHit::~MuonHoughHit ( )

virtualdefault

destructor

Member Function Documentation

calcMagneticTrackRatio()

double MuonHoughHit::calcMagneticTrackRatio ( ) const

private

method that calculates 'magnetic track ratio' (for curved track model)

Definition at line 75 of file MuonHoughHit.cxx.

                                                  {
    // for formulas see muonhoughmathutils, only here for cpu improvements
    double magnetic_trk_ratio = 0.;
 
    if (m_barrel) {
        // Barrel hypothesis
 
        if (m_radius >= MuonHough::radius_cylinder) {
            double diff_b = m_radius - MuonHough::radius_cylinder;
            magnetic_trk_ratio = diff_b * diff_b / m_abs;
        }
    } else {
        int sign = 1;
        if (getHitz() < 0) sign = -1;
        // Endcap hypothesis
        if (std::abs(getHitz()) < MuonHough::z_cylinder) {
            magnetic_trk_ratio = 0.;
        } else if (std::abs(getHitz()) < MuonHough::z_end) {
            // Inside Toroid
            double diff_e = getHitz() - sign * MuonHough::z_cylinder;
            magnetic_trk_ratio = diff_e * diff_e / m_abs;
 
        } else {
            // Outside Toroid
            magnetic_trk_ratio = (-MuonHough::z_magnetic_range_squared + sign * 2 * getHitz() * MuonHough::z_magnetic_range) / m_abs;
        }
    }
 
    return magnetic_trk_ratio;
}

calcPhiSector()

int MuonHoughHit::calcPhiSector ( ) const

private

method that calculates phi sector of hit

Definition at line 65 of file MuonHoughHit.cxx.

                                      {
    double phi = m_phi;  // [-Pi , Pi]
    phi += MuonHough::half_phisector;
 
    if (phi < 0) phi += MuonHough::two_Pi;                                    // [ 0 , 2*Pi ]
    int sectorhit = static_cast<int>(phi / (2 * MuonHough::half_phisector));  // 0..15
    if (sectorhit == MuonHough::phisectors) sectorhit += -1;                  // could happen in rare cases
    return sectorhit;
}

getAbs()

double MuonHoughHit::getAbs ( ) const

inline

returns radius

Definition at line 167 of file MuonHoughHit.h.

{ return m_abs; }

getAssociated()

bool MuonHoughHit::getAssociated ( ) const

inline

return if hit already associated to pattern

Definition at line 176 of file MuonHoughHit.h.

{ return m_associated; }

getDetectorId()

MuonHough::DetectorTechnology MuonHoughHit::getDetectorId ( ) const

inline

return DetectorTechnology

Definition at line 154 of file MuonHoughHit.h.

{ return m_detector_id; }

getHitx()

double MuonHoughHit::getHitx ( ) const

inline

returns x position

Definition at line 159 of file MuonHoughHit.h.

{ return m_pos[Amg::x]; }

getHity()

double MuonHoughHit::getHity ( ) const

inline

returns y position

Definition at line 160 of file MuonHoughHit.h.

{ return m_pos[Amg::y]; }

getHitz()

double MuonHoughHit::getHitz ( ) const

inline

returns z position

Definition at line 161 of file MuonHoughHit.h.

{ return m_pos[Amg::z]; }

getId()

int MuonHoughHit::getId ( ) const

inline

returns id

Definition at line 156 of file MuonHoughHit.h.

{ return m_id; }

getMagneticTrackRatio()

double MuonHoughHit::getMagneticTrackRatio ( ) const

inline

ratio of the tracklength of the particle to which hit might belong would have traversed in magnetic field (for curved track model)

Definition at line 173 of file MuonHoughHit.h.

{ return m_magnetic_trackratio; }

getMeasuresPhi()

bool MuonHoughHit::getMeasuresPhi ( ) const

inline

hit measures phi?

Definition at line 175 of file MuonHoughHit.h.

{ return m_measures_phi; }

getOrigWeight()

double MuonHoughHit::getOrigWeight ( ) const

inline

returns original weight

Definition at line 163 of file MuonHoughHit.h.

{ return m_orig_weight; }

getPhi()

double MuonHoughHit::getPhi ( ) const

inline

returns phi

Definition at line 169 of file MuonHoughHit.h.

{ return m_phi; }

getPosition()

const Amg::Vector3D & MuonHoughHit::getPosition ( ) const

inline

return (x,y,z) vector

Definition at line 157 of file MuonHoughHit.h.

{ return m_pos; }

getPrd()

const Trk::PrepRawData * MuonHoughHit::getPrd ( ) const

inline

return preprawdata

Definition at line 183 of file MuonHoughHit.h.

{ return m_prd; }

getProbability()

double MuonHoughHit::getProbability ( ) const

inline

returns probability that hit is part of pattern (true muon)

Definition at line 164 of file MuonHoughHit.h.

{ return m_probability; }

getRadius()

double MuonHoughHit::getRadius ( ) const

inline

returns radius

Definition at line 166 of file MuonHoughHit.h.

{ return m_radius; }

getTheta()

double MuonHoughHit::getTheta ( ) const

inline

returns theta

Definition at line 168 of file MuonHoughHit.h.

{ return m_theta; }

getWeight()

double MuonHoughHit::getWeight ( ) const

inline

returns weight in histogram after rescaling

Definition at line 162 of file MuonHoughHit.h.

{ return m_weight; }

getWhichDetector()

std::string MuonHoughHit::getWhichDetector

(

)

const

return DetectorTechnology in string

Definition at line 51 of file MuonHoughHit.cxx.

                                               {
    std::string detector_name;
    switch (m_detector_id) {
        case MuonHough::MDT: detector_name = "MDT"; break;
        case MuonHough::CSC: detector_name = "CSC"; break;
        case MuonHough::RPC: detector_name = "RPC"; break;
        case MuonHough::TGC: detector_name = "TGC"; break;
        default:
            MsgStream log(Athena::getMessageSvc(), "MuonHoughHit::getWhichDetector");    
            if (log.level() <= MSG::WARNING) log << MSG::WARNING << "MuonHoughHit:: no valid detector_id" << endmsg;
    }
    return detector_name;
}

isBarrel()

bool MuonHoughHit::isBarrel ( ) const

inline

hit is barrel or endcap (for curved track model)

Definition at line 171 of file MuonHoughHit.h.

{ return m_barrel; }

phiSector()

int MuonHoughHit::phiSector ( ) const

inline

phi sector of hit

Definition at line 172 of file MuonHoughHit.h.

{ return m_phi_sector; }

setAssociated()

void MuonHoughHit::setAssociated ( bool  associated )

inline

set associated

Definition at line 181 of file MuonHoughHit.h.

{ m_associated = associated; }

setId()

void MuonHoughHit::setId ( int  id )

inline

set id

Definition at line 178 of file MuonHoughHit.h.

{ m_id = id; }

setProbability()

void MuonHoughHit::setProbability ( double  prob )

inline

set probability

Definition at line 180 of file MuonHoughHit.h.

{ m_probability = prob; }

setWeight()

void MuonHoughHit::setWeight ( double  weight )

inline

set weight

Definition at line 179 of file MuonHoughHit.h.

{ m_weight = weight; }

Member Data Documentation

m_abs

double MuonHoughHit::m_abs {0.}

protected

absolute value of hit

Definition at line 111 of file MuonHoughHit.h.

m_associated

bool MuonHoughHit::m_associated {false}

protected

hit associated to pattern

Definition at line 137 of file MuonHoughHit.h.

m_barrel

bool MuonHoughHit::m_barrel {true}

protected

hit is barrel / endcap

Definition at line 117 of file MuonHoughHit.h.

m_detector_id

MuonHough::DetectorTechnology MuonHoughHit::m_detector_id {MuonHough::MDT}

protected

detector technology of hit

Definition at line 131 of file MuonHoughHit.h.

m_id

int MuonHoughHit::m_id {0}

protected

unique id

Definition at line 106 of file MuonHoughHit.h.

m_magnetic_trackratio

double MuonHoughHit::m_magnetic_trackratio {0.}

protected

ratio of the tracklength of the particle to which hit might belong would have traversed in magnetic field (for curved track model)

Definition at line 122 of file MuonHoughHit.h.

m_measures_phi

bool MuonHoughHit::m_measures_phi {false}

protected

hit measures phi?

Definition at line 134 of file MuonHoughHit.h.

m_orig_weight

const double MuonHoughHit::m_orig_weight {1.}

protected

original weight

Definition at line 126 of file MuonHoughHit.h.

m_phi

double MuonHoughHit::m_phi {0.}

protected

phi of hit

Definition at line 113 of file MuonHoughHit.h.

m_phi_sector

int MuonHoughHit::m_phi_sector {0}

protected

phi sector (0,15), Atlas Convention

Definition at line 119 of file MuonHoughHit.h.

m_pos

Amg::Vector3D MuonHoughHit::m_pos {Amg::Vector3D::Zero()}

protected

Definition at line 104 of file MuonHoughHit.h.

m_prd

const Trk::PrepRawData* MuonHoughHit::m_prd {nullptr}

protected

Pointer to preprawdata.

Definition at line 102 of file MuonHoughHit.h.

m_probability

double MuonHoughHit::m_probability {1.}

protected

probability that hit belongs to true muon (based on trigger confirmation and crude segment finding (used for mdts))

Definition at line 128 of file MuonHoughHit.h.

m_radius

double MuonHoughHit::m_radius {0.}

protected

radius of hit

Definition at line 109 of file MuonHoughHit.h.

m_theta

double MuonHoughHit::m_theta {0.}

protected

theta of hit

Definition at line 115 of file MuonHoughHit.h.

m_weight

double MuonHoughHit::m_weight {1.}

protected

weight of hit in histogram, after rescaling

Definition at line 124 of file MuonHoughHit.h.

---

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

• MuonHoughHit.h
• MuonHoughHit.cxx