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 >

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 More...

constexpr double	half_phisector = M_PI / phisectors
	angle of half a sector in rad More...

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) More...

constexpr double	radius_cylinder = 4000.
	radius of cylinder More...

constexpr double	z_cylinder = 6000.
	length of cylinder More...

constexpr double	z_end = 15000.
	z value whereafter no magnetic field / curvature More...

constexpr double	z_magnetic_range = z_end - z_cylinder
	range where hit is curved in endcap region More...

constexpr double	z_magnetic_range_squared = z_end * z_end - z_cylinder * z_cylinder
	range where hit is curved in endcap region ('squared') More...

Typedef Documentation

◆ HitVec

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

Definition at line 21 of file MuonLayerHough.h.

◆ PhiHitVec

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

Definition at line 20 of file MuonPhiLayerHough.h.

◆ RegionDescriptionVec

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

Definition at line 51 of file MuonLayerHough.h.

◆ valPair

using MuonHough::valPair = typedef 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 };

◆ 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 };

Function Documentation

◆ extrapolate()

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

Definition at line 519 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 = 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

constexpr double MuonHough::degree_rad_conversion_factor = M_PI / 180.

constexpr

Definition at line 17 of file MuonHoughMathUtils.h.

◆ half_phisector

constexpr double MuonHough::half_phisector = M_PI / phisectors

constexpr

angle of half a sector in rad

Definition at line 21 of file MuonHoughHit.h.

◆ phisectors

constexpr int MuonHough::phisectors = 16

constexpr

number of phi sectors

Definition at line 19 of file MuonHoughHit.h.

◆ rad_degree_conversion_factor

constexpr double MuonHough::rad_degree_conversion_factor = 180. / M_PI

constexpr

Definition at line 18 of file MuonHoughMathUtils.h.

◆ radius_cylinder

constexpr double MuonHough::radius_cylinder = 4000.

constexpr

radius of cylinder

Definition at line 24 of file MuonHoughMathUtils.h.

◆ tan_barrel

constexpr 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

constexpr double MuonHough::two_Pi = 2 * M_PI

constexpr

Definition at line 16 of file MuonHoughMathUtils.h.

◆ z_cylinder

constexpr double MuonHough::z_cylinder = 6000.

constexpr

length of cylinder

Definition at line 26 of file MuonHoughMathUtils.h.

◆ z_end

constexpr double MuonHough::z_end = 15000.

constexpr

z value whereafter no magnetic field / curvature

Definition at line 28 of file MuonHoughMathUtils.h.

◆ z_magnetic_range

constexpr 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

constexpr 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

◆ 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