#include <MdtChamberGeometry.h>

Inheritance diagram for TrkDriftCircleMath::MdtChamberGeometry:

Collaboration diagram for TrkDriftCircleMath::MdtChamberGeometry:

Public Member Functions
	MdtChamberGeometry ()

	MdtChamberGeometry (const Identifier &id, const Muon::IMuonIdHelperSvc *idHelperSvc, unsigned int nml, unsigned int nlay, unsigned int ntubesml0, unsigned int ntubesml1, const LocVec2D &tube0ml0, const LocVec2D &tube0ml1, double tubeDist, double tubeStage, double layDist, double stationTheta)

virtual	~MdtChamberGeometry ()=default

void	init ()

void	setGeometry (unsigned int nml, unsigned int nlay, unsigned int ntubesml0, unsigned int ntubesml1, const LocVec2D &tube0ml0, const LocVec2D &tube0ml1, double tubeDist, double tubeStage, double layDist, double stationTheta)

unsigned int	nml () const

unsigned int	nlay () const override

unsigned int	ntubesml0 () const

unsigned int	ntubesml1 () const

const Identifier &	stationId () const override

const LocVec2D &	firstTubeMl0 () const

const LocVec2D &	firstTubeMl1 () const

double	tubeDist () const

double	tubeStage () const

double	layerDist () const

double	stationTheta () const override

double	tubeRadius () const override

LocVec2D	tubePosition (unsigned int ml, unsigned int lay, unsigned int tube) const override

DCVec	tubesPassedByLine (const Line &line, int ml) const

DCVec	tubesPassedByLine (const Line &line) const override
	Retrieves the vector of drift circles traveresed by the line. More...

void	tubesPassedByLine (const Line &line, int ml, DCVec &crossedTubes) const
	methods taking reference to output vector to reduce memory usage More...

void	tubesPassedByLine (const Line &line, DCVec &crossedTubes) const

void	isSecondMultiLayer (bool isSecond)
	set that this is the second multi layer instead of the first, used in case the first ML is dead More...

bool	validGeometry () const override

void	print (MsgStream &msg) const override

bool	validId (unsigned int ml, unsigned int lay, unsigned int tube) const

Private Member Functions
double	xPosTube (unsigned int ml, unsigned int lay, unsigned int tube) const

double	yPosTube (unsigned int ml, unsigned int lay) const

Private Attributes
Identifier	m_id {}

bool	m_sMdt {false}

bool	m_isBarrel {false}

unsigned int	m_nml {0}

unsigned int	m_nlay {0}

std::vector< unsigned int >	m_ntubesml {}

double	m_tubeDist {0}

double	m_tubeRad {0}

double	m_tubeStage {0}

double	m_layDist {0}

double	m_stationTheta {0}

std::vector< bool >	m_wasInit {}

std::vector< LocVec2D >	m_firstTube {}

bool	m_validGeometry {false}

bool	m_isSecondMultiLayer {false}

Detailed Description

Definition at line 18 of file MdtChamberGeometry.h.

Constructor & Destructor Documentation

◆ MdtChamberGeometry() [1/2]

TrkDriftCircleMath::MdtChamberGeometry::MdtChamberGeometry ( )

Definition at line 20 of file MdtChamberGeometry.cxx.

20 { init(); }

◆ MdtChamberGeometry() [2/2]

TrkDriftCircleMath::MdtChamberGeometry::MdtChamberGeometry	(	const Identifier &	id,
		const Muon::IMuonIdHelperSvc *	idHelperSvc,
		unsigned int	nml,
		unsigned int	nlay,
		unsigned int	ntubesml0,
		unsigned int	ntubesml1,
		const LocVec2D &	tube0ml0,
		const LocVec2D &	tube0ml1,
		double	tubeDist,
		double	tubeStage,
		double	layDist,
		double	stationTheta
	)

Definition at line 22 of file MdtChamberGeometry.cxx.

                                                                                                   :
         m_id{id} {
  
         m_sMdt = idHelperSvc->issMdt(id);
         m_isBarrel = idHelperSvc->mdtIdHelper().isBarrel(id);
         init();
         setGeometry(nml, nlay, ntubesml0, ntubesml1, tube0ml0, tube0ml1, tubeDist, tubeStage, layDist, stationTheta);
     }

◆ ~MdtChamberGeometry()

virtual TrkDriftCircleMath::MdtChamberGeometry::~MdtChamberGeometry ( )

virtualdefault

Member Function Documentation

◆ firstTubeMl0()

const LocVec2D& TrkDriftCircleMath::MdtChamberGeometry::firstTubeMl0 ( ) const

inline

Definition at line 39 of file MdtChamberGeometry.h.

39 { return m_firstTube[0]; }

◆ firstTubeMl1()

const LocVec2D& TrkDriftCircleMath::MdtChamberGeometry::firstTubeMl1 ( ) const

inline

Definition at line 40 of file MdtChamberGeometry.h.

40 { return m_firstTube[1]; }

◆ init()

void TrkDriftCircleMath::MdtChamberGeometry::init ( )

Definition at line 61 of file MdtChamberGeometry.cxx.

                                   {
         m_validGeometry = true;
         if (m_sMdt) {
             m_tubeDist = 15.10;
             m_tubeRad = 7.1;
             m_layDist = 13.085;
         } else {
             m_tubeDist = 30.035;
             m_tubeRad = 14.6;
             m_layDist = 26.011;
         }
         // initialize first tubes to zero
         m_ntubesml.push_back(0);
         m_ntubesml.push_back(0);
         m_firstTube.emplace_back(0., 0.);
         m_firstTube.emplace_back(0., 0.);
         m_wasInit.push_back(true);
         m_wasInit.push_back(true);
     }

◆ isSecondMultiLayer()

void TrkDriftCircleMath::MdtChamberGeometry::isSecondMultiLayer ( bool isSecond )

inline

set that this is the second multi layer instead of the first, used in case the first ML is dead

Definition at line 61 of file MdtChamberGeometry.h.

61 { m_isSecondMultiLayer = isSecond; }

◆ layerDist()

double TrkDriftCircleMath::MdtChamberGeometry::layerDist ( ) const

inline

Definition at line 44 of file MdtChamberGeometry.h.

44 { return m_layDist; }

◆ nlay()

unsigned int TrkDriftCircleMath::MdtChamberGeometry::nlay ( ) const

inlineoverridevirtual

Implements TrkDriftCircleMath::ChamberGeometry.

Definition at line 33 of file MdtChamberGeometry.h.

33 { return m_nlay; }

◆ nml()

unsigned int TrkDriftCircleMath::MdtChamberGeometry::nml ( ) const

inline

Definition at line 32 of file MdtChamberGeometry.h.

32 { return m_nml; }

◆ ntubesml0()

unsigned int TrkDriftCircleMath::MdtChamberGeometry::ntubesml0 ( ) const

inline

Definition at line 34 of file MdtChamberGeometry.h.

34 { return m_ntubesml[0]; }

◆ ntubesml1()

unsigned int TrkDriftCircleMath::MdtChamberGeometry::ntubesml1 ( ) const

inline

Definition at line 35 of file MdtChamberGeometry.h.

35 { return m_ntubesml[1]; }

◆ print()

void TrkDriftCircleMath::MdtChamberGeometry::print ( MsgStream & msg ) const

overridevirtual

Implements TrkDriftCircleMath::ChamberGeometry.

Definition at line 170 of file MdtChamberGeometry.cxx.

                                                        {
         msg << MSG::ALWAYS << " MdtChamberGeometry " << m_id <<std::endl
             << "  nml " << m_nml << " nlay " << m_nlay << " ntube1 " << m_ntubesml[0] << " ntube2 " << m_ntubesml[1] << std::endl
             << " pos ml1 " << m_firstTube[0] << "  ml2 " << m_firstTube[1] << std::endl
             << " tubeDist " << m_tubeDist << " tubeStage " << m_tubeStage << " layDist " << m_layDist << " tubeRad " << m_tubeRad
                   << endmsg;
     }

◆ setGeometry()

void TrkDriftCircleMath::MdtChamberGeometry::setGeometry	(	unsigned int	nml,
		unsigned int	nlay,
		unsigned int	ntubesml0,
		unsigned int	ntubesml1,
		const LocVec2D &	tube0ml0,
		const LocVec2D &	tube0ml1,
		double	tubeDist,
		double	tubeStage,
		double	layDist,
		double	stationTheta
	)

Definition at line 33 of file MdtChamberGeometry.cxx.

                                                                               {
         m_nml = nml;
         m_nlay = nlay;
         m_ntubesml.clear();
         m_ntubesml.push_back(ntubesml0);
         m_ntubesml.push_back(ntubesml1);
         m_tubeDist = tubeDist;
         m_tubeStage = tubeStage;
         m_layDist = layDist;
  
         m_firstTube[0] = tube0ml0;
         m_firstTube[1] = tube0ml1;
  
         if (m_nml < 1 || m_nml > 2)
             THROW_EXCEPTION("MdtChamberGeometry::setGeometry() - got called with nml="<<m_nml<<" which is definitely out of range.");
         if (m_nlay < 1 || m_nlay > 4)
             THROW_EXCEPTION("MdtChamberGeometry::setGeometry() - got called with nlay="<<m_nlay<<" which is definitely out of range");
         if (ntubesml0 < 1 || ntubesml0 > maxNTubesPerLayer)
             THROW_EXCEPTION("MdtChamberGeometry::setGeometry() - got called with ntubesml0="<<ntubesml0<<" which is definitely out of range");
         // there can be chambers with only 1 multilayer. Then, the second multilayer will have ntubesml1=0
         if (ntubesml1 > maxNTubesPerLayer)
             THROW_EXCEPTION("MdtChamberGeometry::setGeometry() - got called with ntubesml1="<<ntubesml1<<" which is definitely out of range");
  
         m_stationTheta = stationTheta;
     }

◆ stationId()

const Identifier& TrkDriftCircleMath::MdtChamberGeometry::stationId ( ) const

inlineoverridevirtual

Implements TrkDriftCircleMath::ChamberGeometry.

Definition at line 37 of file MdtChamberGeometry.h.

37 { return m_id; }

◆ stationTheta()

double TrkDriftCircleMath::MdtChamberGeometry::stationTheta ( ) const

inlineoverridevirtual

Implements TrkDriftCircleMath::ChamberGeometry.

Definition at line 46 of file MdtChamberGeometry.h.

46 { return m_stationTheta; }

◆ tubeDist()

double TrkDriftCircleMath::MdtChamberGeometry::tubeDist ( ) const

inline

Definition at line 42 of file MdtChamberGeometry.h.

42 { return m_tubeDist; }

◆ tubePosition()

LocVec2D TrkDriftCircleMath::MdtChamberGeometry::tubePosition	(	unsigned int	ml,
		unsigned int	lay,
		unsigned int	tube
	)		const

overridevirtual

Implements TrkDriftCircleMath::ChamberGeometry.

Definition at line 177 of file MdtChamberGeometry.cxx.

                                                                                                         {
         if (!validId(ml, lay, tube)) {
             THROW_EXCEPTION("Combination of multilayer ml: "<<ml<<", layer: "<<lay<<" and tube: "<<tube<<" given ");
         }
         LocVec2D tube_vec{xPosTube(ml, lay, tube), yPosTube(ml, lay)};
         return tube_vec;
     }

◆ tubeRadius()

double TrkDriftCircleMath::MdtChamberGeometry::tubeRadius ( ) const

inlineoverridevirtual

Implements TrkDriftCircleMath::ChamberGeometry.

Definition at line 48 of file MdtChamberGeometry.h.

48 { return m_tubeRad; }

◆ tubesPassedByLine() [1/4]

DCVec TrkDriftCircleMath::MdtChamberGeometry::tubesPassedByLine ( const Line & line ) const

inlineoverridevirtual

Retrieves the vector of drift circles traveresed by the line.

Implements TrkDriftCircleMath::ChamberGeometry.

Definition at line 53 of file MdtChamberGeometry.h.

53 { return tubesPassedByLine(line, -1); }

◆ tubesPassedByLine() [2/4]

void TrkDriftCircleMath::MdtChamberGeometry::tubesPassedByLine	(	const Line &	line,
		DCVec &	crossedTubes
	)		const

inline

Definition at line 57 of file MdtChamberGeometry.h.

57 { return tubesPassedByLine(line, -1, crossedTubes); }

◆ tubesPassedByLine() [3/4]

DCVec TrkDriftCircleMath::MdtChamberGeometry::tubesPassedByLine	(	const Line &	line,
		int	ml
	)		const

Definition at line 160 of file MdtChamberGeometry.cxx.

                                                                                         {
         DCVec crossedTubes;
  
         if (!m_validGeometry) return crossedTubes;
  
         tubesPassedByLine(line, inMultilayer, crossedTubes);
  
         std::stable_sort(crossedTubes.begin(), crossedTubes.end(), SortDcsByY());
         return crossedTubes;
     }

◆ tubesPassedByLine() [4/4]

void TrkDriftCircleMath::MdtChamberGeometry::tubesPassedByLine	(	const Line &	line,
		int	ml,
		DCVec &	crossedTubes
	)		const

methods taking reference to output vector to reduce memory usage

Definition at line 103 of file MdtChamberGeometry.cxx.

                                                                                                             {
         crossedTubes.reserve(50);
         ResidualWithLine resLine{line};
         const LocVec2D& linepos = line.position();
         const LocVec2D& linedir = line.direction();
         double dxdy = std::abs(linedir.y()) > 0.0001 ? linedir.x() / linedir.y() : linedir.x() / 0.0001;
         for (unsigned int ml = 0; ml < m_nml; ++ml) {
             // check whether geometry was initialized for given multilayer
             if (!m_wasInit[ml]) continue;
  
             // if indicated only scan single multilayer
             if (inMultilayer != -1 && inMultilayer != (int)ml) { continue; }
             for (unsigned int lay = 0; lay < m_nlay; ++lay) {
                 double ylay = yPosTube(ml, lay);
                 double xfirsttube = xPosTube(ml, lay, 0);
                 double xintersect = dxdy * (ylay - linepos.y()) + linepos.x();
                 double relpos = (xintersect - xfirsttube) / m_tubeDist;
                 int ctube = (int)relpos;
                 if (ctube < 0) ctube = 0;
                 if (ctube >= (int)m_ntubesml[ml]) ctube = m_ntubesml[ml] - 1;
  
                 if (inMultilayer != -1)
  
                     for (int i = ctube - 1; i >= 0; --i) {
                         const LocVec2D& lp = tubePosition(ml, lay, i);
                         double res = resLine.residual(lp);
                         if (std::abs(res) > m_tubeRad) {
                             if (m_tubeDist > 0) {
                                 if (res > m_tubeRad) break;
                             } else {
                                 if (res < -m_tubeRad) break;
                             }
                         } else {
                             // if this is a chamber with only the second ml, set the ml index accordingly
                             unsigned int actualMl = m_isSecondMultiLayer ? 1 : ml;
                             crossedTubes.emplace_back(lp, m_tubeRad, res, DriftCircle::EmptyTube, MdtId(m_isBarrel, actualMl, lay, i),
                                                       nullptr);
                         }
                     }
                 for (int i = ctube; i < (int)m_ntubesml[ml]; ++i) {
                     const LocVec2D& lp = tubePosition(ml, lay, i);
                     double res = resLine.residual(lp);
                     if (std::abs(res) > m_tubeRad) {
                         if (m_tubeDist > 0) {
                             if (res < -m_tubeRad) break;
                         } else {
                             if (res > m_tubeRad) break;
                         }
                     } else {
                         unsigned int actualMl = m_isSecondMultiLayer ? 1 : ml;
                         crossedTubes.emplace_back(lp, m_tubeRad, res, DriftCircle::EmptyTube, MdtId(m_isBarrel, actualMl, lay, i), nullptr);
                     }
                 }
             }
         }
     }

◆ tubeStage()

double TrkDriftCircleMath::MdtChamberGeometry::tubeStage ( ) const

inline

Definition at line 43 of file MdtChamberGeometry.h.

43 { return m_tubeStage; }

◆ validGeometry()

bool TrkDriftCircleMath::MdtChamberGeometry::validGeometry ( ) const

inlineoverridevirtual

Implements TrkDriftCircleMath::ChamberGeometry.

Definition at line 63 of file MdtChamberGeometry.h.

63 { return m_validGeometry; }

◆ validId()

bool TrkDriftCircleMath::MdtChamberGeometry::validId	(	unsigned int	ml,
		unsigned int	lay,
		unsigned int	tube
	)		const

Definition at line 81 of file MdtChamberGeometry.cxx.

                                                                                                {
         if (!m_validGeometry) return false;
         if (ml > 1) {
             MsgStream msg{Athena::getMessageSvc(), "MdtChamberGeometry::validId"};
             msg <<MSG::ERROR << " Wrong index: ml " << ml << " max " << m_nml << endmsg;
             print(msg);
            
             return false;
         } else if (lay > m_nlay) {
             MsgStream msg{Athena::getMessageSvc(), "MdtChamberGeometry::validId"};
             msg <<MSG::ERROR <<" Wrong index: lay " << lay << " max " << m_nlay << endmsg;
             print(msg);
             return false;
         } else if (tube > m_ntubesml[ml]) {
             MsgStream msg{Athena::getMessageSvc(), "MdtChamberGeometry::validId"};            
             msg << " wrong index: tube " << tube << " max " << m_ntubesml[ml] << endmsg;
             print(msg);
             return false;
         }
         return true;
     }

◆ xPosTube()

double TrkDriftCircleMath::MdtChamberGeometry::xPosTube	(	unsigned int	ml,
		unsigned int	lay,
		unsigned int	tube
	)		const

inlineprivate

Definition at line 185 of file MdtChamberGeometry.cxx.

                                                                                                          {
         double xpos = tube * m_tubeDist + m_firstTube[ml].x();
  
         if (lay % 2 == 1) {
             if (m_nlay == 4 && ml == 1)
                 xpos -= m_tubeStage;
             else
                 xpos += m_tubeStage;
         }
  
         return xpos;
     }

◆ yPosTube()

double TrkDriftCircleMath::MdtChamberGeometry::yPosTube	(	unsigned int	ml,
		unsigned int	lay
	)		const

private

Definition at line 198 of file MdtChamberGeometry.cxx.

198 { return lay * m_layDist + m_firstTube[ml].y(); }

Member Data Documentation

◆ m_firstTube

std::vector<LocVec2D> TrkDriftCircleMath::MdtChamberGeometry::m_firstTube {}

private

Definition at line 86 of file MdtChamberGeometry.h.

◆ m_id

Identifier TrkDriftCircleMath::MdtChamberGeometry::m_id {}

private

Definition at line 73 of file MdtChamberGeometry.h.

◆ m_isBarrel

bool TrkDriftCircleMath::MdtChamberGeometry::m_isBarrel {false}

private

Definition at line 75 of file MdtChamberGeometry.h.

◆ m_isSecondMultiLayer

bool TrkDriftCircleMath::MdtChamberGeometry::m_isSecondMultiLayer {false}

private

Definition at line 88 of file MdtChamberGeometry.h.

◆ m_layDist

double TrkDriftCircleMath::MdtChamberGeometry::m_layDist {0}

private

Definition at line 82 of file MdtChamberGeometry.h.

◆ m_nlay

unsigned int TrkDriftCircleMath::MdtChamberGeometry::m_nlay {0}

private

Definition at line 77 of file MdtChamberGeometry.h.

◆ m_nml

unsigned int TrkDriftCircleMath::MdtChamberGeometry::m_nml {0}

private

Definition at line 76 of file MdtChamberGeometry.h.

◆ m_ntubesml

std::vector<unsigned int> TrkDriftCircleMath::MdtChamberGeometry::m_ntubesml {}

private

Definition at line 78 of file MdtChamberGeometry.h.

◆ m_sMdt

bool TrkDriftCircleMath::MdtChamberGeometry::m_sMdt {false}

private

Definition at line 74 of file MdtChamberGeometry.h.

◆ m_stationTheta

double TrkDriftCircleMath::MdtChamberGeometry::m_stationTheta {0}

private

Definition at line 83 of file MdtChamberGeometry.h.

◆ m_tubeDist

double TrkDriftCircleMath::MdtChamberGeometry::m_tubeDist {0}

private

Definition at line 79 of file MdtChamberGeometry.h.

◆ m_tubeRad

double TrkDriftCircleMath::MdtChamberGeometry::m_tubeRad {0}

private

Definition at line 80 of file MdtChamberGeometry.h.

◆ m_tubeStage

double TrkDriftCircleMath::MdtChamberGeometry::m_tubeStage {0}

private

Definition at line 81 of file MdtChamberGeometry.h.

◆ m_validGeometry

bool TrkDriftCircleMath::MdtChamberGeometry::m_validGeometry {false}

private

Definition at line 87 of file MdtChamberGeometry.h.

◆ m_wasInit

std::vector<bool> TrkDriftCircleMath::MdtChamberGeometry::m_wasInit {}

private

Definition at line 85 of file MdtChamberGeometry.h.

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

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ MdtChamberGeometry() [1/2]

◆ MdtChamberGeometry() [2/2]

◆ ~MdtChamberGeometry()

Member Function Documentation

◆ firstTubeMl0()

◆ firstTubeMl1()

◆ init()

◆ isSecondMultiLayer()

◆ layerDist()

◆ nlay()

◆ nml()

◆ ntubesml0()

◆ ntubesml1()

◆ print()

◆ setGeometry()

◆ stationId()

◆ stationTheta()

◆ tubeDist()

◆ tubePosition()

◆ tubeRadius()

◆ tubesPassedByLine() [1/4]

◆ tubesPassedByLine() [2/4]

◆ tubesPassedByLine() [3/4]

◆ tubesPassedByLine() [4/4]

◆ tubeStage()

◆ validGeometry()

◆ validId()

◆ xPosTube()

◆ yPosTube()

Member Data Documentation

◆ m_firstTube

◆ m_id

◆ m_isBarrel

◆ m_isSecondMultiLayer

◆ m_layDist

◆ m_nlay

◆ m_nml

◆ m_ntubesml

◆ m_sMdt

◆ m_stationTheta

◆ m_tubeDist

◆ m_tubeRad

◆ m_tubeStage

◆ m_validGeometry

◆ m_wasInit