#include <MMT_struct.h>

Inheritance diagram for MMT_Parameters:

Collaboration diagram for MMT_Parameters:

Public Member Functions
	MMT_Parameters (std::string layerSetup, char wedgeSize, const MuonGM::MuonDetectorManager *detManager)

std::vector< ROOT::Math::XYZVector >	MM_firststrip_positions (const MuonGM::MuonDetectorManager *detManager, const std::string &wedge, int eta) const

bool	is_x (int plane) const

bool	is_u (int plane) const

bool	is_v (int plane) const

char	getSector () const

double	getlWidth () const

double	getPitch () const

double	getLowerBoundEta1 () const

double	getLowerBoundEta2 () const

double	getMissedBottomEtaStrips () const

double	getMissedBottomStereoStrips () const

double	getYbase (int plane, int eta) const

bool	msgLvl (const MSG::Level lvl) const
	Test the output level. More...

MsgStream &	msg () const
	The standard message stream. More...

MsgStream &	msg (const MSG::Level lvl) const
	The standard message stream. More...

void	setLevel (MSG::Level lvl)
	Change the current logging level. More...

Private Member Functions
void	initMessaging () const
	Initialize our message level and MessageSvc. More...

Private Attributes
std::array< std::array< float, 2 >, 8 >	m_ybases

char	m_sector

double	m_lWidth

double	m_pitch

double	m_innerRadiusEta1

double	m_innerRadiusEta2

double	m_missedBottomEta

double	m_missedBottomStereo

std::string	m_nm
	Message source name. More...

boost::thread_specific_ptr< MsgStream >	m_msg_tls
	MsgStream instance (a std::cout like with print-out levels) More...

std::atomic< IMessageSvc * >	m_imsg { nullptr }
	MessageSvc pointer. More...

std::atomic< MSG::Level >	m_lvl { MSG::NIL }
	Current logging level. More...

std::atomic_flag m_initialized	ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
	Messaging initialized (initMessaging) More...

Detailed Description

Definition at line 24 of file MMT_struct.h.

Constructor & Destructor Documentation

◆ MMT_Parameters()

MMT_Parameters::MMT_Parameters	(	std::string	layerSetup,
		char	wedgeSize,
		const MuonGM::MuonDetectorManager *	detManager
	)

Definition at line 11 of file MMT_struct.cxx.

   : AthMessaging(Athena::getMessageSvc(), "MMT_Parameters") {
  
   // Get the modules for each multiplet, the sector and side arguement (third and fifth argument, respectively) shouldn't matter
   // since the algorithm computes locally
   char side   = 'A';
   m_sector = wedgeSize;
   MMDetectorHelper aHelper;
   MMDetectorDescription* mm_top_mult1    = aHelper.Get_MMDetector(m_sector, 2, 5, 1, side); // Parameters: sector, eta, phi, layer/multiplet
   MMDetectorDescription* mm_top_mult2    = aHelper.Get_MMDetector(m_sector, 2, 5, 2, side);
   MMDetectorDescription* mm_bottom_mult1 = aHelper.Get_MMDetector(m_sector, 1, 5, 1, side);
   MMDetectorDescription* mm_bottom_mult2 = aHelper.Get_MMDetector(m_sector, 1, 5, 2, side);
  
   MMReadoutParameters roParam_top_mult1    = mm_top_mult1->GetReadoutParameters();
   MMReadoutParameters roParam_top_mult2    = mm_top_mult2->GetReadoutParameters();
   MMReadoutParameters roParam_bottom_mult1 = mm_bottom_mult1->GetReadoutParameters();
   MMReadoutParameters roParam_bottom_mult2 = mm_bottom_mult2->GetReadoutParameters();
  
   //Needed to get the max channels for the deadZone in the x planes... a little annoying...
   //Only changes for small or large and per module. Phi, layer, etc don't matter (just set to 1)
   std::string wedgeString = (wedgeSize=='L' ? "MML" : "MMS");
  
   m_pitch = roParam_bottom_mult1.stripPitch;
   m_lWidth = mm_top_mult1->lWidth();
   m_innerRadiusEta1 = roParam_bottom_mult1.distanceFromZAxis;
   m_innerRadiusEta2 = roParam_top_mult1.distanceFromZAxis;
   m_missedBottomEta = roParam_bottom_mult1.nMissedBottomEta;
   m_missedBottomStereo = roParam_bottom_mult1.nMissedBottomStereo;
   if (m_sector == 'L') {
     ATH_MSG_DEBUG("LM1 \n" <<
                  "\t\t\t\t Total Strips: " << roParam_bottom_mult1.tStrips << " with pitch: " << roParam_bottom_mult1.stripPitch << "\n" <<
                  "\t\t\t\t KO strips TopEta: " << roParam_bottom_mult1.nMissedTopEta << " - BottomEta: " << roParam_bottom_mult1.nMissedBottomEta << "\n" <<
                  "\t\t\t\t KO strips TopStereo: " << roParam_bottom_mult1.nMissedTopStereo << " - BottomStereo: " << roParam_bottom_mult1.nMissedBottomStereo << "\n" <<
                  "\t\t\t\t (Top, Bottom, Length): (" << mm_bottom_mult1->lWidth() << ", " << mm_bottom_mult1->sWidth() << ", " << mm_bottom_mult1->Length() << ")\n" <<
                  "\t\t\t\t Zpos - Distance from ZAxis: " << roParam_bottom_mult1.zpos << " - " << roParam_bottom_mult1.distanceFromZAxis << "\n" <<
                  "\t\t\t\t dlStereoTop/Bottom: " << roParam_bottom_mult1.dlStereoTop << " " << roParam_bottom_mult1.dlStereoBottom << "\n" <<
                  "\t\t\t\t Active area --> (Bottom, Top, Height) : (" << roParam_bottom_mult1.activeBottomLength << ", " << roParam_bottom_mult1.activeTopLength << ", " << roParam_bottom_mult1.activeH << ")");
   } else if (m_sector == 'S') {
     ATH_MSG_DEBUG("SM1 \n" <<
                  "\t\t\t\t KO strips TopEta: " << roParam_bottom_mult1.nMissedTopEta << " - BottomEta: " << roParam_bottom_mult1.nMissedBottomEta << "\n" <<
                  "\t\t\t\t KO strips TopStereo: " << roParam_bottom_mult1.nMissedTopStereo << " - BottomStereo: " << roParam_bottom_mult1.nMissedBottomStereo << "\n" <<
                  "\t\t\t\t Total Strips: " << roParam_bottom_mult1.tStrips << " with pitch: " << roParam_bottom_mult1.stripPitch << "\n" <<
                  "\t\t\t\t (Top, Bottom, Length): (" << mm_bottom_mult1->lWidth() << ", " << mm_bottom_mult1->sWidth() << ", " << mm_bottom_mult1->Length() << ")\n" <<
                  "\t\t\t\t Zpos - Distance from ZAxis: " << roParam_bottom_mult1.zpos << " - " << roParam_bottom_mult1.distanceFromZAxis << "\n" <<
                  "\t\t\t\t dlStereoTop/Bottom: " << roParam_bottom_mult1.dlStereoTop << " " << roParam_bottom_mult1.dlStereoBottom << "\n" <<
                  "\t\t\t\t Active area --> (Bottom, Top, Height) : (" << roParam_bottom_mult1.activeBottomLength << ", " << roParam_bottom_mult1.activeTopLength << ", " << roParam_bottom_mult1.activeH << ")");
   }
  
   if (m_sector == 'L') {
     ATH_MSG_DEBUG("LM2 \n" <<
                  "\t\t\t\t Total Strips: " << roParam_top_mult1.tStrips << " with pitch: " << roParam_top_mult1.stripPitch << "\n" <<
                  "\t\t\t\t KO strips TopEta: " << roParam_top_mult1.nMissedTopEta << " - BottomEta: " << roParam_top_mult1.nMissedBottomEta << "\n" <<
                  "\t\t\t\t KO strips TopStereo: " << roParam_top_mult1.nMissedTopStereo << " - BottomStereo: " << roParam_top_mult1.nMissedBottomStereo << "\n" <<
                  "\t\t\t\t (Top, Bottom, Length): (" << mm_top_mult1->lWidth() << ", " << mm_top_mult1->sWidth() << ", " << mm_top_mult1->Length() << ")\n" <<
                  "\t\t\t\t Zpos - Distance from ZAxis: " << roParam_top_mult1.zpos << " - " << roParam_top_mult1.distanceFromZAxis << "\n" <<
                  "\t\t\t\t dlStereoTop/Bottom: " << roParam_top_mult1.dlStereoTop << " " << roParam_top_mult1.dlStereoBottom << "\n" <<
                  "\t\t\t\t Active area --> (Top, Bottom, Height) : (" << roParam_top_mult1.activeTopLength << ", " << roParam_top_mult1.activeBottomLength << ", " << roParam_top_mult1.activeH << ")");
   } else if (m_sector == 'S') {
     ATH_MSG_DEBUG("SM2 \n" <<
                  "\t\t\t\t KO strips TopEta: " << roParam_top_mult1.nMissedTopEta << " - BottomEta: " << roParam_top_mult1.nMissedBottomEta << "\n" <<
                  "\t\t\t\t KO strips TopStereo: " << roParam_top_mult1.nMissedTopStereo << " - BottomStereo: " << roParam_top_mult1.nMissedBottomStereo << "\n" <<
                  "\t\t\t\t (Top, Bottom, Length): (" << mm_top_mult1->lWidth() << ", " << mm_top_mult1->sWidth() << ", " << mm_top_mult1->Length() << ")\n" <<
                  "\t\t\t\t Zpos - Distance from ZAxis: " << roParam_top_mult1.zpos << " - " << roParam_top_mult1.distanceFromZAxis << "\n" <<
                  "\t\t\t\t dlStereoTop/Bottom: " << roParam_top_mult1.dlStereoTop << " " << roParam_top_mult1.dlStereoBottom << "\n" <<
                  "\t\t\t\t Active area --> (Top, Bottom, Height) : (" << roParam_top_mult1.activeTopLength << ", " << roParam_top_mult1.activeBottomLength << ", " << roParam_top_mult1.activeH << ")");
     for (const auto &angle : roParam_top_mult1.stereoAngle) ATH_MSG_DEBUG("Stereo angle: " << angle);
   }
  
   // retrieve the z-position of the planes
   std::vector<double> z_nominal;
   int eta = 1;
   for (const auto& pos: MM_firststrip_positions(detManager, wedgeString, eta)) z_nominal.push_back(pos.Z());
  
   if(z_nominal.size() != layerSetup.size()){
     ATH_MSG_WARNING("Number of planes in setup is "<< layerSetup.size() << ", but we have a nominal " << z_nominal.size() << " planes.");
     throw std::runtime_error("MMT_Parameters: Invalid number of planes");
   }
  
   float stereo_degree = Gaudi::Units::deg*roParam_top_mult1.stereoAngle[2];
  
   // MM_firststrip_positions returns the position for phi sector 1.
   // => for the small sectors, rotate this by -1/16 of a rotation to make our lives easier.
   // in this coordinate basis, x is up/down the wedge (radial), and y is left/right (phi).
   float cos_rotation = std::cos(-2*M_PI / 16.0);
   float sin_rotation = std::sin(-2*M_PI / 16.0);
   float x_rotated = 0.0;
   float y_rotated = 0.0;
  
   float st_angle = 0.0;
  
   float radial_pos      = 0;
   float radial_pos_xx_1 = 0, radial_pos_xx_2 = 0;
   float radial_pos_uv_1 = 0, radial_pos_uv_2 = 0;
  
   for (unsigned int eta = 1; eta <= 2; eta++){
     unsigned int layer = 1;
     for (const auto& pos: MM_firststrip_positions(detManager, wedgeString, eta)){
  
       if (wedgeString=="MMS"){
         x_rotated = pos.X()*cos_rotation - pos.Y()*sin_rotation;
         y_rotated = pos.X()*sin_rotation + pos.Y()*cos_rotation;
       }
       else{
         x_rotated = pos.X();
         y_rotated = pos.Y();
       }
  
       if      (is_u(layer)) st_angle = -1*std::abs(stereo_degree);
       else if (is_v(layer)) st_angle =    std::abs(stereo_degree);
       else                  st_angle = 0;
  
       // walk from the center of the strip to the position of the strip at the center of the wedge.
       // NB: for X-planes, this is simply the center of the strip: tan(0) = 0.
       radial_pos = std::abs(x_rotated - y_rotated*std::tan(st_angle * M_PI/180.0));
  
       if (is_x(layer) && eta==1) radial_pos_xx_1 = radial_pos;
       else if (is_x(layer) && eta==2) radial_pos_xx_2 = radial_pos;
       if (is_u(layer) && eta==1) radial_pos_uv_1 = radial_pos;
       else if (is_u(layer) && eta==2) radial_pos_uv_2 = radial_pos;
  
       layer++;
     }
   }
  
   // store radial positions as fixed point
   std::array<float, 2> radial_pos_xx = {radial_pos_xx_1, radial_pos_xx_2};
   std::array<float, 2> radial_pos_uv = {radial_pos_uv_1, radial_pos_uv_2};
   m_ybases[0] = radial_pos_xx;
   m_ybases[1] = radial_pos_xx;
   m_ybases[2] = radial_pos_uv;
   m_ybases[3] = radial_pos_uv;
   m_ybases[4] = radial_pos_uv;
   m_ybases[5] = radial_pos_uv;
   m_ybases[6] = radial_pos_xx;
   m_ybases[7] = radial_pos_xx;
 }

Member Function Documentation

◆ getLowerBoundEta1()

double MMT_Parameters::getLowerBoundEta1 ( ) const

inline

Definition at line 46 of file MMT_struct.h.

46 { return m_innerRadiusEta1; }

◆ getLowerBoundEta2()

double MMT_Parameters::getLowerBoundEta2 ( ) const

inline

Definition at line 47 of file MMT_struct.h.

47 { return m_innerRadiusEta2; }

◆ getlWidth()

double MMT_Parameters::getlWidth ( ) const

inline

Definition at line 44 of file MMT_struct.h.

44 { return m_lWidth; }

◆ getMissedBottomEtaStrips()

double MMT_Parameters::getMissedBottomEtaStrips ( ) const

inline

Definition at line 48 of file MMT_struct.h.

48 { return m_missedBottomEta; }

◆ getMissedBottomStereoStrips()

double MMT_Parameters::getMissedBottomStereoStrips ( ) const

inline

Definition at line 49 of file MMT_struct.h.

49 { return m_missedBottomStereo; }

◆ getPitch()

double MMT_Parameters::getPitch ( ) const

inline

Definition at line 45 of file MMT_struct.h.

45 { return m_pitch; }

◆ getSector()

char MMT_Parameters::getSector ( ) const

inline

Definition at line 43 of file MMT_struct.h.

43 { return m_sector; }

◆ getYbase()

double MMT_Parameters::getYbase	(	int	plane,
		int	eta
	)		const

inline

Definition at line 50 of file MMT_struct.h.

50 { return m_ybases[plane][eta]; }

◆ initMessaging()

void AthMessaging::initMessaging ( ) const

privateinherited

Initialize our message level and MessageSvc.

This method should only be called once.

Definition at line 39 of file AthMessaging.cxx.

 {
   m_imsg = Athena::getMessageSvc();
   m_lvl = m_imsg ?
     static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
     MSG::INFO;
 }

◆ is_u()

bool MMT_Parameters::is_u ( int plane ) const

inline

Definition at line 34 of file MMT_struct.h.

                                {
       static const std::set<unsigned int> planes_u {3, 5};
       return planes_u.count(plane);
     }

◆ is_v()

bool MMT_Parameters::is_v ( int plane ) const

inline

Definition at line 38 of file MMT_struct.h.

                                {
       static const std::set<unsigned int> planes_v{4, 6};
       return planes_v.count(plane);
     }

◆ is_x()

bool MMT_Parameters::is_x ( int plane ) const

inline

Definition at line 30 of file MMT_struct.h.

                                {
       static const std::set<unsigned int> planes_x{1, 2, 7, 8};
       return planes_x.count(plane);
     }

◆ MM_firststrip_positions()

std::vector< ROOT::Math::XYZVector > MMT_Parameters::MM_firststrip_positions	(	const MuonGM::MuonDetectorManager *	detManager,
		const std::string &	wedge,
		int	eta
	)		const

Definition at line 148 of file MMT_struct.cxx.

                                                                                                                                                        {
  
   std::vector<ROOT::Math::XYZVector> positions;
   positions.reserve(8);
  
   Identifier strip_id;
   int phi = 1;
   int strip = 200;
   unsigned int n_multiplets = 2;
   unsigned int n_layers     = 4;
  
   for (unsigned int mult = 1; mult <= n_multiplets; mult++) {
     for (unsigned int layer = 1; layer <= n_layers; layer++) {
  
       Amg::Vector3D pos(0.0, 0.0, 0.0);
       strip_id = detManager->mmIdHelper()->channelID(wedge, eta, phi, mult, layer, strip);
       const MuonGM::MMReadoutElement* readout = detManager->getMMReadoutElement(strip_id);
       ROOT::Math::XYZVector coordinates(0.,0.,0.);
       if (readout->stripGlobalPosition(strip_id, pos)) coordinates.SetXYZ(pos.x(), pos.y(), pos.z());
       else ATH_MSG_WARNING("Wedge " << wedge << " phi " << phi << " multiplet " << mult << " layer " << layer <<  " | Unable to retrieve global positions");
       positions.push_back(coordinates);
  
       ATH_MSG_DEBUG( "global z-pos. (using MMReadoutElement) for"
                            << " wedge size "       << wedge
                            << " multiplet "        << mult
                            << " layer "            << layer
                            << " | z = "            << coordinates.Z() );
     }
   }
  
    return positions;
 }

◆ msg() [1/2]

MsgStream & AthMessaging::msg ( ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 164 of file AthMessaging.h.

 {
   MsgStream* ms = m_msg_tls.get();
   if (!ms) {
     if (!m_initialized.test_and_set()) initMessaging();
     ms = new MsgStream(m_imsg,m_nm);
     m_msg_tls.reset( ms );
   }
  
   ms->setLevel (m_lvl);
   return *ms;
 }

◆ msg() [2/2]

MsgStream & AthMessaging::msg ( const MSG::Level lvl ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 179 of file AthMessaging.h.

180 { return msg() << lvl; }

◆ msgLvl()

bool AthMessaging::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Test the output level.

Parameters

lvl	The message level to test against

Returns: boolean Indicating if messages at given level will be printed

Return values

true	Messages at level "lvl" will be printed

Definition at line 151 of file AthMessaging.h.

 {
   if (!m_initialized.test_and_set()) initMessaging();
   if (m_lvl <= lvl) {
     msg() << lvl;
     return true;
   } else {
     return false;
   }
 }

◆ setLevel()

void AthMessaging::setLevel ( MSG::Level lvl )

inherited

Change the current logging level.

Use this rather than msg().setLevel() for proper operation with MT.

Definition at line 28 of file AthMessaging.cxx.

 {
   m_lvl = lvl;
 }

Member Data Documentation

◆ ATLAS_THREAD_SAFE

std::atomic_flag m_initialized AthMessaging::ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT

mutableprivateinherited

Messaging initialized (initMessaging)

Definition at line 141 of file AthMessaging.h.

◆ m_imsg

std::atomic<IMessageSvc*> AthMessaging::m_imsg { nullptr }

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

◆ m_innerRadiusEta1

double MMT_Parameters::m_innerRadiusEta1

private

Definition at line 55 of file MMT_struct.h.

◆ m_innerRadiusEta2

double MMT_Parameters::m_innerRadiusEta2

private

Definition at line 55 of file MMT_struct.h.

◆ m_lvl

std::atomic<MSG::Level> AthMessaging::m_lvl { MSG::NIL }

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

◆ m_lWidth

double MMT_Parameters::m_lWidth

private

Definition at line 55 of file MMT_struct.h.

◆ m_missedBottomEta

double MMT_Parameters::m_missedBottomEta

private

Definition at line 56 of file MMT_struct.h.

◆ m_missedBottomStereo

double MMT_Parameters::m_missedBottomStereo

private

Definition at line 56 of file MMT_struct.h.

◆ m_msg_tls

boost::thread_specific_ptr<MsgStream> AthMessaging::m_msg_tls

mutableprivateinherited

MsgStream instance (a std::cout like with print-out levels)

Definition at line 132 of file AthMessaging.h.

◆ m_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

◆ m_pitch

double MMT_Parameters::m_pitch

private

Definition at line 55 of file MMT_struct.h.

◆ m_sector

char MMT_Parameters::m_sector

private

Definition at line 54 of file MMT_struct.h.

◆ m_ybases

std::array<std::array<float, 2>, 8> MMT_Parameters::m_ybases

private

Definition at line 53 of file MMT_struct.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

◆ MMT_Parameters()

Member Function Documentation

◆ getLowerBoundEta1()

◆ getLowerBoundEta2()

◆ getlWidth()

◆ getMissedBottomEtaStrips()

◆ getMissedBottomStereoStrips()

◆ getPitch()

◆ getSector()

◆ getYbase()

◆ initMessaging()

◆ is_u()

◆ is_v()

◆ is_x()

◆ MM_firststrip_positions()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ setLevel()

Member Data Documentation

◆ ATLAS_THREAD_SAFE

◆ m_imsg

◆ m_innerRadiusEta1

◆ m_innerRadiusEta2

◆ m_lvl

◆ m_lWidth

◆ m_missedBottomEta

◆ m_missedBottomStereo

◆ m_msg_tls

◆ m_nm

◆ m_pitch

◆ m_sector

◆ m_ybases