InDet::InDetMaterialVeto Class Reference

#include <InDetMaterialVeto.h>

Inheritance diagram for InDet::InDetMaterialVeto:

Collaboration diagram for InDet::InDetMaterialVeto:

Public Member Functions
	InDetMaterialVeto (const BeamPipeDetectorManager *beamPipeMgr, const InDetDD::PixelDetectorManager *pixelManager)
std::unique_ptr< TH2D >	ITkPixMaterialMap ()
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::unique_ptr< TH2D >	m_ITkPixMaterialMap
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 29 of file InDetMaterialVeto.h.

Constructor & Destructor Documentation

InDetMaterialVeto()

InDet::InDetMaterialVeto::InDetMaterialVeto	(	const BeamPipeDetectorManager *	beamPipeMgr,
		const InDetDD::PixelDetectorManager *	pixelManager
	)

Definition at line 22 of file InDetMaterialVeto.cxx.

                                                                        :
    AthMessaging(Gaudi::svcLocator()->service< IMessageSvc >( "MessageSvc" ),"InDetMaterialVeto")
  {
 
    ATH_MSG_DEBUG("Building material rejection map from BeamPipe and PixelDetectorManager");
    
    static constexpr int nbins_R = 76;
    double bins_R[nbins_R+1];
    for(unsigned int i=0; i<=15; i++) bins_R[i] = 2*i - 1.;          // 2 mm bin width below R=29 mm, centered on R=24 mm (beampipe)
    bins_R[16] = 30.; // 29-30 makes junction with pixel bins, veto IPT R=29.3 mm
    static constexpr double Rbinwidth_pixel = 6.;
    for(unsigned int i=1; i<=60; i++) bins_R[i+16] = 30 + Rbinwidth_pixel*i; // 6 mm bin width beyond R=30 mm
    
    static constexpr int nbins_Z = 1000;
    static constexpr double zmax = 3000.;
    static constexpr double zbinwidth = 2*zmax/nbins_Z;
    
    std::string mapName = "ITkMaterialMap";
    m_ITkPixMaterialMap = std::make_unique<TH2D>(mapName.c_str(),mapName.c_str(),nbins_Z,-zmax,zmax,nbins_R,bins_R); // x-axis = global z coordinates, -3240 mm to +3240 mm, 6 mm bin width / y-axis = global R coordinates, variable bin width
    
    
    // Retrieve the beam pipe radius from the SectionC03 volume
    
    double beamPipeRadius = 0.;
    PVConstLink beamPipeTopVolume =  beamPipeMgr->getTreeTop(0);
    const GeoLogVol* beamPipeLogVolume = beamPipeTopVolume->getLogVol();
    const GeoTube* beamPipeTube = nullptr;
    
    if (beamPipeLogVolume){
      beamPipeTube = dynamic_cast<const GeoTube*>(beamPipeLogVolume->getShape());
      if (beamPipeTube){
    
    for(unsigned int i=0;i<beamPipeTopVolume->getNChildVols();i++){
      
      if(beamPipeTopVolume->getNameOfChildVol(i)=="SectionC03"){
        PVConstLink childTopVolume =  beamPipeTopVolume->getChildVol(i);
        const GeoLogVol* childLogVolume = childTopVolume->getLogVol();
        const GeoTube* childTube = nullptr;
        
        if (childLogVolume){
          childTube = dynamic_cast<const GeoTube*>(childLogVolume->getShape());
          if (childTube){
        beamPipeRadius = 0.5 * (childTube->getRMax()+childTube->getRMin());
          }
        }
        
        break; //Exit loop after SectionC03 is found
      }
      
    } // Loop over child volumes
    
      }
    } // if(beamPipeLogVolume)
    
    ATH_MSG_DEBUG("BeamPipeRadius used for material rejection="<<beamPipeRadius);
    
    // Fill map with beam pipe radius for all z
    for(double z = -zmax + 0.5*zbinwidth; z<zmax; z+=zbinwidth) m_ITkPixMaterialMap->Fill(z,beamPipeRadius);
    
    
    // Retrieve IPT radius
    // Labelled as anonymous volume so taken as pixel volume with smallest radius
    
    double IPTRadius = -1.;
    
    PVConstLink pixelTopVolume =  pixelManager->getTreeTop(0);
    
    for(unsigned int i=0;i<pixelTopVolume->getNChildVols();i++){
      
      PVConstLink childTopVolume =  pixelTopVolume->getChildVol(i);
      const GeoLogVol* childLogVolume = childTopVolume->getLogVol();
      const GeoTube* childTube = nullptr;
      
      if (childLogVolume){
    childTube = dynamic_cast<const GeoTube*>(childLogVolume->getShape());
    if (childTube){
      double radius = 0.5*(childTube->getRMin()+childTube->getRMax());
      if(IPTRadius<0. || radius<IPTRadius) IPTRadius = radius;
    }
      }
      
    }
    
    ATH_MSG_DEBUG("IPTRadius used for material rejection="<<IPTRadius);
    
    // Fill map with IPT radius for all z
    for(double z = -zmax + 0.5*zbinwidth; z<zmax; z+=zbinwidth) m_ITkPixMaterialMap->Fill(z,IPTRadius);
    
 
    InDetDD::SiDetectorElementCollection::const_iterator iter;
    for (iter = pixelManager->getDetectorElementBegin(); iter != pixelManager->getDetectorElementEnd(); ++iter) {
      // get the ID
      Identifier Pixel_ModuleID = (*iter)->identify();
      // check the validity
      if (Pixel_ModuleID.is_valid()) {
    const InDetDD::SiDetectorElement *module = pixelManager->getDetectorElement(Pixel_ModuleID);
    
    //Take into account full module extent
    InDetDD::DetectorType type = module->design().type();
    
    if(type==InDetDD::PixelInclined){
      double zMin_mod = module->zMin();
      double zMax_mod = module->zMax();
      double rMin_mod = module->rMin();
      double rMax_mod = module->rMax();
      if(zMin_mod>0){
        double alpha = (rMin_mod-rMax_mod)/(zMax_mod-zMin_mod);
        for(double z=zMin_mod; z<=zMax_mod; z+=zbinwidth){
          double R = rMax_mod + alpha*(z-zMin_mod);
          m_ITkPixMaterialMap->Fill(z,R);
        }
      }
      else{
        double alpha = (rMax_mod-rMin_mod)/(zMax_mod-zMin_mod);
        for(double z=zMin_mod; z<=zMax_mod; z+=zbinwidth){
          double R = rMin_mod + alpha*(z-zMin_mod);
          m_ITkPixMaterialMap->Fill(z,R);
        }
      }
    }
    
    else if(type==InDetDD::PixelBarrel){
      double zMin_mod = module->zMin();
      double zMax_mod = module->zMax();
      double R = std::hypot(module->center().x(),module->center().y());
      for(double z=zMin_mod; z<=zMax_mod; z+=zbinwidth) m_ITkPixMaterialMap->Fill(z,R);
    }
    
    else if(type==InDetDD::PixelEndcap){
      double rMin_mod = module->rMin();
      double rMax_mod = module->rMax();
      double z = module->center().z();
      for(double R=rMin_mod; R<=rMax_mod; R+=Rbinwidth_pixel) m_ITkPixMaterialMap->Fill(z,R);
    }
    
      }
    } // end loop over pixel modules
 
  }

Member Function Documentation

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

ITkPixMaterialMap()

std::unique_ptr<TH2D> InDet::InDetMaterialVeto::ITkPixMaterialMap ( )

inline

Definition at line 34 of file InDetMaterialVeto.h.

{ return std::move(m_ITkPixMaterialMap); }

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.

{ 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_ITkPixMaterialMap

std::unique_ptr<TH2D> InDet::InDetMaterialVeto::m_ITkPixMaterialMap

private

Definition at line 37 of file InDetMaterialVeto.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_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.

---

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

• InDetMaterialVeto.h
• InDetMaterialVeto.cxx