Classes
class	Hit
	struct containing all hit information needed for the Hough transform More...
class	HitDebugInfo
	struct containing additional debug information on the hits that is not needed for the actual alg but very useful for debugging More...
struct	MuonDebugInfo
	struct containing truth or track information More...
class	MuonDetectorDescription
	class managing geometry of the Hough spaces More...
class	MuonDetectorHough
	class managing all Hough transforms in the detector More...
struct	MuonLayerHough
class	MuonLayerHoughSelector
struct	MuonPhiLayerHough
class	MuonSectorHough
	class managing all precision Hough transforms in a sector More...
struct	PhiHit
	struct containing all hit information needed for the Hough transform More...
struct	RegionDescriptor
	struct containing all information to build a Hough transform for a given chamber index More...
struct	SegDebugInfo
	struct containing truth or segment information More...
struct	SortHitsPerLayer
	struct to sort the hits More...

Typedefs
using	HitVec = std::vector<std::shared_ptr<MuonHough::Hit>>
using	RegionDescriptionVec = std::vector<RegionDescriptor>
using	PhiHitVec = std::vector<std::shared_ptr<MuonHough::PhiHit>>
using	valPair = std::pair<int, float>
using	ChIdx = Muon::MuonStationIndex::ChIndex

Enumerations
enum	DetectorTechnology { MDT , CSC , RPC , TGC }
	enum to identify the muondetectortechnology More...
enum	MuonHoughTransformers { hough_xy , hough_rzcosmics , hough_curved_at_a_cylinder , hough_rz , hough_rz_mdt , hough_rz_rpc , hough_yz }
	enum to identify the different houghtransformers More...

Functions
float	extrapolate (const MuonLayerHough::Maximum &ref, const MuonLayerHough::Maximum &ex, bool doparabolic=false)

Variables
constexpr int	phisectors = 16
	number of phi sectors
constexpr double	half_phisector = M_PI / phisectors
	angle of half a sector in rad
constexpr double	two_Pi = 2 * M_PI
constexpr double	degree_rad_conversion_factor = M_PI / 180.
constexpr double	rad_degree_conversion_factor = 180. / M_PI
constexpr double	tan_barrel = 11430. / 14000.
	relation for transition between endcap and barrel 11.43 m (r) / 14m (z)
constexpr double	radius_cylinder = 4000.
	radius of cylinder
constexpr double	z_cylinder = 6000.
	length of cylinder
constexpr double	z_end = 15000.
	z value whereafter no magnetic field / curvature
constexpr double	z_magnetic_range = z_end - z_cylinder
	range where hit is curved in endcap region
constexpr double	z_magnetic_range_squared = z_end * z_end - z_cylinder * z_cylinder
	range where hit is curved in endcap region ('squared')

Typedef Documentation

◆ ChIdx

using MuonHough::ChIdx = Muon::MuonStationIndex::ChIndex

Definition at line 9 of file MuonRegionHough.cxx.

◆ HitVec

using MuonHough::HitVec = std::vector<std::shared_ptr<MuonHough::Hit>>

Definition at line 21 of file MuonLayerHough.h.

◆ PhiHitVec

using MuonHough::PhiHitVec = std::vector<std::shared_ptr<MuonHough::PhiHit>>

Definition at line 20 of file MuonPhiLayerHough.h.

◆ RegionDescriptionVec

using MuonHough::RegionDescriptionVec = std::vector<RegionDescriptor>

Definition at line 56 of file MuonLayerHough.h.

◆ valPair

using MuonHough::valPair = std::pair<int, float>

Definition at line 14 of file MuonLayerHoughSelector.cxx.

Enumeration Type Documentation

◆ DetectorTechnology

enum MuonHough::DetectorTechnology

enum to identify the muondetectortechnology

Enumerator
MDT
CSC
RPC
TGC

Definition at line 17 of file MuonHoughHit.h.

17{ MDT, CSC, RPC, TGC };

MuonHough::CSC

@ CSC

Definition MuonHoughHit.h:17

MuonHough::RPC

@ RPC

Definition MuonHoughHit.h:17

MuonHough::MDT

@ MDT

Definition MuonHoughHit.h:17

TGC

Definition TgcBase.h:6

◆ MuonHoughTransformers

enum MuonHough::MuonHoughTransformers

enum to identify the different houghtransformers

Enumerator
hough_xy
hough_rzcosmics
hough_curved_at_a_cylinder
hough_rz
hough_rz_mdt
hough_rz_rpc
hough_yz

Definition at line 14 of file MuonHoughPattern.h.

14{ hough_xy, hough_rzcosmics, hough_curved_at_a_cylinder, hough_rz, hough_rz_mdt, hough_rz_rpc, hough_yz };

MuonHough::hough_xy

@ hough_xy

Definition MuonHoughPattern.h:14

MuonHough::hough_rz_mdt

@ hough_rz_mdt

Definition MuonHoughPattern.h:14

MuonHough::hough_rz_rpc

@ hough_rz_rpc

Definition MuonHoughPattern.h:14

MuonHough::hough_rz

@ hough_rz

Definition MuonHoughPattern.h:14

MuonHough::hough_curved_at_a_cylinder

@ hough_curved_at_a_cylinder

Definition MuonHoughPattern.h:14

MuonHough::hough_rzcosmics

@ hough_rzcosmics

Definition MuonHoughPattern.h:14

MuonHough::hough_yz

@ hough_yz

Definition MuonHoughPattern.h:14

Function Documentation

◆ extrapolate()

float MuonHough::extrapolate	(	const MuonLayerHough::Maximum &	ref,
		const MuonLayerHough::Maximum &	ex,
		bool	doparabolic = false )

Definition at line 521 of file MuonLayerHough.cxx.

                                                                                                           {
        // z is always the precision plane. r is the reference plane, simple and robust
        double ref_z = ref.getGlobalZ();
        double ref_r = ref.getGlobalR();
        double ex_z = ex.getGlobalZ();
        double ex_r = ex.getGlobalR();
        float theta_ref = ref.getGlobalTheta();
        if (!doparabolic || ref_z == 0 || theta_ref == 0) {  // do linear extrapolation
            if (!ex.isEndcap()) {                            // if extrapolate to barell
                return ex_z - ex_r / ref_r * ref_z;
            } else {  // if extrapolate to endcap
                return ex_r - ex_z * ref_r / ref_z;
            }
        } else {  // do parabolic
            float expected = 0;
            float extrapolated_diff = 9999;
            const float tan_theta_ref = std::tan(theta_ref);
            const float invtan_theta_ref = 1. / tan_theta_ref;
            float r_start = static_cast<int>(ref.hough->m_descriptor.chIndex) % 2 > 0
                                ? 4900.
                                : 5200.;  // start of barrel B field; values could be further optimized; 5500.:6500.
            float z_start = 8500.;        // start of endcap B field; used to be 6500; should start at 8500
            float z_end = 12500.;         // end of endcap B field; used to be 12000; should end at 12500
            float r_SL = ref_r + (ex_z - ref_z) * tan_theta_ref;
            float z_SL = ref_z + (ex_r - ref_r) * invtan_theta_ref;
            // start extrapolation
            if (!ref.isEndcap()) {  // this is starting with barrel chamber; BEE is 6
                float rgeo = ref_r * ref_r - r_start * r_start;
                float rhoInv = (invtan_theta_ref * ref_r - ref_z) / rgeo;
                if (!ex.isEndcap()) {  // this ends with barrel chamber; BEE is 6
                    expected = ref_z + (ex_r - ref_r) * invtan_theta_ref + (ex_r - ref_r) * (ex_r - ref_r) * rhoInv;
                    extrapolated_diff = ex_z - expected;
                }       // end with barrel to barrel extrapolation
                else {  // this is a barrel to endcap extrapolation, mostly in the transition region
                    // recalculate z start
                    z_start = ref_z + (r_start - ref_r) * invtan_theta_ref + (r_start - ref_r) * (r_start - ref_r) * rhoInv;
                    float zgeo = ref_z * ref_z - z_start * z_start;
                    float rho = (tan_theta_ref * ref_z - ref_r) / zgeo;
                    expected = ref_r + (ex_z - ref_z) * tan_theta_ref + (ex_z - ref_z) * (ex_z - ref_z) * rho;
                    extrapolated_diff = ex_r - expected;
                }
            } else {  // this starts with endcap chamber;
                // swap the starting position if on the other side
                if (tan_theta_ref < 0) {
                    z_start = -z_start;
                    z_end = -z_end;
                }
                if (ex.isEndcap()) {                          // extrapolate to endcap
                    if (std::abs(ref_z) < std::abs(z_end)) {  // extrapolate from EI or EE, have to use linear
                        expected = r_SL;
                    } else {                                       // from EM or EO or EE
                        if (std::abs(ex_z) > std::abs(z_start)) {  // extrapolate to EM or EO
                            // extrapolate to outer layer, just using theta of the middle measurement; only works if the theta measurement
                            // is correct can extrapolate with either the outside theta or middle theta; outside theta is better; farther
                            // away from the B field
                            expected = r_SL;
                        } else {  // to EI
                            float r_end = ref_r + (z_end - ref_z) * tan_theta_ref;
                            float zgeo = z_start * z_start - z_end * z_end;
                            float rhoInv = (r_end - z_end * tan_theta_ref) / zgeo;
                            float tantheta = tan_theta_ref - 2 * (z_end - z_start) * rhoInv;
                            expected = ex_z * tantheta + (ex_z - z_start) * (ex_z - z_start) * rhoInv;
                        }
                    }
                    extrapolated_diff = ex_r - expected;
                } else {  // exrapolate to barrel; again verly likely to be in transition region, complicated B field; just use linear
                    expected = z_SL;
                    extrapolated_diff = ex_z - expected;
                }
            }
            return extrapolated_diff;
        }  // end of parabolic extrapolation
    }

Variable Documentation

◆ degree_rad_conversion_factor

double MuonHough::degree_rad_conversion_factor = M_PI / 180.

constexpr

Definition at line 17 of file MuonHoughMathUtils.h.

◆ half_phisector

double MuonHough::half_phisector = M_PI / phisectors

constexpr

angle of half a sector in rad

Definition at line 21 of file MuonHoughHit.h.

◆ phisectors

int MuonHough::phisectors = 16

constexpr

number of phi sectors

Definition at line 19 of file MuonHoughHit.h.

◆ rad_degree_conversion_factor

double MuonHough::rad_degree_conversion_factor = 180. / M_PI

constexpr

Definition at line 18 of file MuonHoughMathUtils.h.

◆ radius_cylinder

double MuonHough::radius_cylinder = 4000.

constexpr

radius of cylinder

Definition at line 24 of file MuonHoughMathUtils.h.

◆ tan_barrel

double MuonHough::tan_barrel = 11430. / 14000.

constexpr

relation for transition between endcap and barrel 11.43 m (r) / 14m (z)

Definition at line 22 of file MuonHoughMathUtils.h.

◆ two_Pi

double MuonHough::two_Pi = 2 * M_PI

constexpr

Definition at line 16 of file MuonHoughMathUtils.h.

◆ z_cylinder

double MuonHough::z_cylinder = 6000.

constexpr

length of cylinder

Definition at line 26 of file MuonHoughMathUtils.h.

◆ z_end

double MuonHough::z_end = 15000.

constexpr

z value whereafter no magnetic field / curvature

Definition at line 28 of file MuonHoughMathUtils.h.

◆ z_magnetic_range

double MuonHough::z_magnetic_range = z_end - z_cylinder

constexpr

range where hit is curved in endcap region

Definition at line 30 of file MuonHoughMathUtils.h.

◆ z_magnetic_range_squared

double MuonHough::z_magnetic_range_squared = z_end * z_end - z_cylinder * z_cylinder

constexpr

range where hit is curved in endcap region ('squared')

Definition at line 32 of file MuonHoughMathUtils.h.

Classes

Typedefs

Enumerations

Functions

Variables

Typedef Documentation

◆ ChIdx

◆ HitVec

◆ PhiHitVec

◆ RegionDescriptionVec

◆ valPair

Enumeration Type Documentation

◆ DetectorTechnology

◆ MuonHoughTransformers

Function Documentation

◆ extrapolate()

Variable Documentation

◆ degree_rad_conversion_factor

◆ half_phisector

◆ phisectors

◆ rad_degree_conversion_factor

◆ radius_cylinder

◆ tan_barrel

◆ two_Pi

◆ z_cylinder

◆ z_end

◆ z_magnetic_range

◆ z_magnetic_range_squared