#include <CmsTrackerDetIdBuilder.h>

Public Member Functions
GeometricDet *	buildId (GeometricDet *det)

	CmsTrackerDetIdBuilder (std::vector< int > detidShifts)

Protected Member Functions
void	iterate (GeometricDet *det, int level, unsigned int ID)

Private Attributes
std::array< int, nSubDet *maxLevels >	m_detidshifts

std::map< std::string, uint32_t >	m_mapNavTypeToDetId

Static Private Attributes
static const int	maxLevels =6

static const unsigned int	nSubDet =6

Detailed Description

Class to build a geographicalId.

Definition at line 15 of file CmsTrackerDetIdBuilder.h.

Constructor & Destructor Documentation

CmsTrackerDetIdBuilder::CmsTrackerDetIdBuilder ( std::vector< int > detidShifts )

Definition at line 13 of file CmsTrackerDetIdBuilder.cc.

References mps_fire::i, m_detidshifts, maxLevels, and nSubDet.

   : m_detidshifts()
 {
   if(detidShifts.size()!=nSubDet*maxLevels) 
     edm::LogError("WrongConfiguration") << "Wrong configuration of TrackerGeometricDetESModule. Vector of " 
                     << detidShifts.size() << " elements provided"; 
   else {
     for(unsigned int i=0;i<nSubDet*maxLevels;++i) { m_detidshifts[i]=detidShifts[i];}
   }
 }

Member Function Documentation

GeometricDet * CmsTrackerDetIdBuilder::buildId ( GeometricDet * det )

Definition at line 25 of file CmsTrackerDetIdBuilder.cc.

References GeometricDet::geographicalID(), iterate(), LogDebug, DetId::rawId(), GeometricDet::setGeographicalID(), lumiQTWidget::t, trackingTruthProducer_cfi::tracker, and DetId::Tracker.

Referenced by DDDCmsTrackerContruction::construct().

 {
 
   LogDebug("BuildingTrackerDetId") << "Starting to build Tracker DetIds";
 
   DetId t( DetId::Tracker, 0 );
   tracker->setGeographicalID( t );
   iterate( tracker, 0, tracker->geographicalID().rawId() );
 
   return tracker;
 }

void CmsTrackerDetIdBuilder::iterate	(	GeometricDet *	det,
		int	level,
		unsigned int	ID
	)

protected

Definition at line 38 of file CmsTrackerDetIdBuilder.cc.

References funct::abs(), GeometricDet::component(), makeMuonMisalignmentScenario::components, GeometricDet::components(), GeometricDet::geographicalID(), mps_fire::i, recoMuon::in, hcalDigis_cfi::level, LogTrace, m_detidshifts, RecoTauDiscriminantConfiguration::mask, maxLevels, nSubDet, DetId::rawId(), GeometricDet::setGeographicalID(), AlCaHLTBitMon_QueryRunRegistry::string, groupFilesInBlocks::temp, and GeometricDet::translation().

Referenced by buildId().

 {
   std::bitset<32> binary_ID(ID);
 
   // SubDetector (useful to know fron now on, valid only after level 0, where SubDetector is assigned)
   uint32_t mask = (7<<25);
   uint32_t iSubDet = ID & mask;
   iSubDet = iSubDet >> 25;
   //
 
   LogTrace("BuildingTrackerDetId") << std::string(2*level,'-') 
                    << "+" << ID << " " << iSubDet << " " << level;
   
   switch( level )
     {
       // level 0: special case because it is used to assign the proper detid bits based on the endcap-like subdetector position: +z or -z
     case 0:
       {  
     for( uint32_t i = 0; i<(in)->components().size(); i++ )
       {
         GeometricDet* component = in->component(i);
         uint32_t iSubDet = component->geographicalID().rawId();
         uint32_t temp = ID;
         temp |= (iSubDet<<25);
         component->setGeographicalID(temp); 
         
         if(iSubDet>0 && iSubDet<=nSubDet && m_detidshifts[level*nSubDet+iSubDet-1]>=0) {
           if(m_detidshifts[level*nSubDet+iSubDet-1]+2<25) temp|= (0<<(m_detidshifts[level*nSubDet+iSubDet-1]+2)); 
           bool negside = component->translation().z()<0.;
           if(std::abs(component->translation().z())<1.) negside = component->components().front()->translation().z()<0.; // needed for subdet like TID which are NOT translated
           LogTrace("BuildingTrackerDetId") << "Is negative endcap? " << negside 
                            << ", because z translation is " << component->translation().z()
                            << " and component z translation is " << component->components().front()->translation().z();
           if(negside)
         {
           temp |= (1<<m_detidshifts[level*nSubDet+iSubDet-1]);
         }
           else
         {
           temp |= (2<<m_detidshifts[level*nSubDet+iSubDet-1]);
         }
         }
         component->setGeographicalID(DetId(temp));  
         
         // next level
         iterate(component,level+1,((in)->components())[i]->geographicalID().rawId());
       } 
     break;
       }
       // level 1 to 5
     default:
       {
     for( uint32_t i = 0; i < (in)->components().size(); i++ )
       {
         auto component = in->component(i);
         uint32_t temp = ID;
         
         if(level<maxLevels) {
           if(iSubDet>0 && iSubDet <=nSubDet && m_detidshifts[level*nSubDet+iSubDet-1]>=0) {
         temp |= (component->geographicalID().rawId()<<m_detidshifts[level*nSubDet+iSubDet-1]); 
           }
           component->setGeographicalID( temp );
           // next level
           iterate(component,level+1,((in)->components())[i]->geographicalID().rawId());      
         }
       }
     
     break; 
       }    
       // level switch ends
     }
   
   return;
   
 }