Builds ideal, real and misaligned geometries. More...

Inheritance diagram for CTPPSGeometryESModule:

Public Member Functions
	CTPPSGeometryESModule (const edm::ParameterSet &)

std::unique_ptr< DetGeomDesc >	produceIdealGD (const IdealGeometryRecord &)

std::unique_ptr< DetGeomDesc >	produceMisalignedGD (const VeryForwardMisalignedGeometryRecord &)

std::unique_ptr< CTPPSGeometry >	produceMisalignedTG (const VeryForwardMisalignedGeometryRecord &)

std::unique_ptr< DetGeomDesc >	produceRealGD (const VeryForwardRealGeometryRecord &)

std::unique_ptr< CTPPSGeometry >	produceRealTG (const VeryForwardRealGeometryRecord &)

	~CTPPSGeometryESModule () override

Public Member Functions inherited from edm::ESProducer
	ESProducer ()

	~ESProducer ()(false) override

Public Member Functions inherited from edm::ESProxyFactoryProducer
	ESProxyFactoryProducer ()

void	newInterval (const eventsetup::EventSetupRecordKey &iRecordType, const ValidityInterval &iInterval) override
	overrides DataProxyProvider method More...

	~ESProxyFactoryProducer () noexcept(false) override

Public Member Functions inherited from edm::eventsetup::DataProxyProvider
	DataProxyProvider ()

const ComponentDescription &	description () const

bool	isUsingRecord (const EventSetupRecordKey &) const

const KeyedProxies &	keyedProxies (const EventSetupRecordKey &iRecordKey) const

void	resetProxies (const EventSetupRecordKey &iRecordType)

void	resetProxiesIfTransient (const EventSetupRecordKey &iRecordType)

void	setAppendToDataLabel (const edm::ParameterSet &)

void	setDescription (const ComponentDescription &iDescription)

std::set< EventSetupRecordKey >	usingRecords () const

virtual	~DataProxyProvider () noexcept(false)

Static Protected Member Functions
static void	applyAlignments (const edm::ESHandle< DetGeomDesc > &, const edm::ESHandle< RPAlignmentCorrectionsData > &, DetGeomDesc *&)

static void	buildDetGeomDesc (DDFilteredView fv, DetGeomDesc gd)

Protected Attributes
std::string	compactViewTag_

unsigned int	verbosity_

Additional Inherited Members
Public Types inherited from edm::eventsetup::DataProxyProvider
typedef std::vector< std::pair< DataKey, edm::propagate_const< std::shared_ptr< DataProxy > > > >	KeyedProxies

typedef std::vector< EventSetupRecordKey >	Keys

typedef std::map< EventSetupRecordKey, KeyedProxies >	RecordProxies

Static Public Member Functions inherited from edm::eventsetup::DataProxyProvider
static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::ESProducer
template<typename T >
void	setWhatProduced (T *iThis, const es::Label &iLabel=es::Label())

template<typename T >
void	setWhatProduced (T iThis, const char iLabel)

template<typename T >
void	setWhatProduced (T *iThis, const std::string &iLabel)

template<typename T , typename TDecorator >
void	setWhatProduced (T *iThis, const TDecorator &iDec, const es::Label &iLabel=es::Label())

template<typename T , typename TReturn , typename TRecord >
void	setWhatProduced (T iThis, TReturn(T::iMethod)(const TRecord &), const es::Label &iLabel=es::Label())

template<typename T , typename TReturn , typename TRecord , typename TArg >
void	setWhatProduced (T iThis, TReturn(T::iMethod)(const TRecord &), const TArg &iDec, const es::Label &iLabel=es::Label())

Protected Member Functions inherited from edm::ESProxyFactoryProducer
template<class TFactory >
void	registerFactory (std::unique_ptr< TFactory > iFactory, const std::string &iLabel=std::string())

virtual void	registerFactoryWithKey (const eventsetup::EventSetupRecordKey &iRecord, std::unique_ptr< eventsetup::ProxyFactoryBase > iFactory, const std::string &iLabel=std::string())

void	registerProxies (const eventsetup::EventSetupRecordKey &iRecord, KeyedProxies &aProxyList) override
	override DataProxyProvider method More...

Protected Member Functions inherited from edm::eventsetup::DataProxyProvider
void	eraseAll (const EventSetupRecordKey &iRecordKey)
	deletes all the Proxies in aStream More...

void	invalidateProxies (const EventSetupRecordKey &iRecordKey)

template<class T >
void	usingRecord ()

void	usingRecordWithKey (const EventSetupRecordKey &)

Detailed Description

Builds ideal, real and misaligned geometries.

First, it creates a tree of DetGeomDesc from DDCompView. For real and misaligned geometries, it applies alignment corrections (RPAlignmentCorrections) found in corresponding ...GeometryRecord.

Second, it creates CTPPSGeometry from DetGeoDesc tree.

Definition at line 50 of file CTPPSGeometryESModule.cc.

Constructor & Destructor Documentation

CTPPSGeometryESModule::CTPPSGeometryESModule ( const edm::ParameterSet & iConfig )

Definition at line 76 of file CTPPSGeometryESModule.cc.

References produceIdealGD(), produceMisalignedGD(), produceMisalignedTG(), produceRealGD(), produceRealTG(), and edm::ESProducer::setWhatProduced().

                                                                            :
   verbosity_( iConfig.getUntrackedParameter<unsigned int>( "verbosity", 1 ) ),
   compactViewTag_( iConfig.getParameter<std::string>( "compactViewTag" ) )
 {
   setWhatProduced( this, &CTPPSGeometryESModule::produceIdealGD );
 
   setWhatProduced( this, &CTPPSGeometryESModule::produceRealGD );
   setWhatProduced( this, &CTPPSGeometryESModule::produceRealTG );
 
   setWhatProduced( this, &CTPPSGeometryESModule::produceMisalignedGD );
   setWhatProduced( this, &CTPPSGeometryESModule::produceMisalignedTG );
 }

CTPPSGeometryESModule::~CTPPSGeometryESModule ( )

inlineoverride

Definition at line 54 of file CTPPSGeometryESModule.cc.

References applyAlignments(), buildDetGeomDesc(), produceIdealGD(), produceMisalignedGD(), produceMisalignedTG(), produceRealGD(), and produceRealTG().

54 {}

Member Function Documentation

void CTPPSGeometryESModule::applyAlignments	(	const edm::ESHandle< DetGeomDesc > &	idealGD,
		const edm::ESHandle< RPAlignmentCorrectionsData > &	alignments,
		DetGeomDesc *&	newGD
	)

staticprotected

Definition at line 92 of file CTPPSGeometryESModule.cc.

References DetGeomDesc::addComponent(), DetGeomDesc::ApplyAlignment(), edmScanValgrind::buffer, DetGeomDesc::components(), DDD_CTPPS_DIAMONDS_RP_NAME, DDD_CTPPS_DIAMONDS_SEGMENT_NAME, DDD_CTPPS_PIXELS_RP_NAME, DDD_CTPPS_PIXELS_SENSOR_NAME, DDD_CTPPS_UFSD_SEGMENT_NAME, DDD_TOTEM_RP_RP_NAME, DDD_TOTEM_RP_SENSOR_NAME, DetGeomDesc::geographicalID(), RPAlignmentCorrectionsData::getFullSensorCorrection(), RPAlignmentCorrectionsData::getRPCorrection(), mps_fire::i, edm::ESHandleBase::isValid(), DDName::name(), dataset::name, DetGeomDesc::name(), edm::ESHandle< T >::product(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by produceMisalignedGD(), produceRealGD(), and ~CTPPSGeometryESModule().

 {
   newGD = new DetGeomDesc( *(idealGD.product()) );
   std::deque<const DetGeomDesc*> buffer;
   std::deque<DetGeomDesc*> bufferNew;
   buffer.emplace_back( idealGD.product() );
   bufferNew.emplace_back( newGD );
 
   while ( !buffer.empty() ) {
     const DetGeomDesc* sD = buffer.front();
     DetGeomDesc* pD = bufferNew.front();
     buffer.pop_front();
     bufferNew.pop_front();
 
     const std::string name = pD->name().name();
 
     // Is it sensor? If yes, apply full sensor alignments
     if ( name == DDD_TOTEM_RP_SENSOR_NAME
       || name == DDD_CTPPS_DIAMONDS_SEGMENT_NAME
       || name == DDD_CTPPS_UFSD_SEGMENT_NAME
       || name == DDD_CTPPS_PIXELS_SENSOR_NAME ) {
       unsigned int plId = pD->geographicalID();
 
       if ( alignments.isValid() ) {
         const RPAlignmentCorrectionData& ac = alignments->getFullSensorCorrection( plId );
         pD->ApplyAlignment( ac );
       }
     }
 
     // Is it RP box? If yes, apply RP alignments
     if ( name == DDD_TOTEM_RP_RP_NAME
       || name == DDD_CTPPS_DIAMONDS_RP_NAME
       || name == DDD_CTPPS_PIXELS_RP_NAME ) {
       unsigned int rpId = pD->geographicalID();
 
       if ( alignments.isValid() ) {
         const RPAlignmentCorrectionData& ac = alignments->getRPCorrection( rpId );
         pD->ApplyAlignment( ac );
       }
     }
 
     // create and add children
     for ( unsigned int i = 0; i < sD->components().size(); i++ ) {
       const DetGeomDesc* sDC = sD->components()[i];
       buffer.emplace_back( sDC );
 
       // create new node with the same information as in sDC and add it as a child of pD
       DetGeomDesc* cD = new DetGeomDesc( *sDC );
       pD->addComponent( cD );
 
       bufferNew.emplace_back( cD );
     }
   }
 }

void CTPPSGeometryESModule::buildDetGeomDesc	(	DDFilteredView *	fv,
		DetGeomDesc *	gd
	)

staticprotected

Definition at line 152 of file CTPPSGeometryESModule.cc.

References DetGeomDesc::addComponent(), DDFilteredView::copyno(), DDFilteredView::copyNumbers(), DDD_CTPPS_DIAMONDS_RP_NAME, DDD_CTPPS_DIAMONDS_SEGMENT_NAME, DDD_CTPPS_PIXELS_RP_NAME, DDD_CTPPS_PIXELS_SENSOR_NAME, DDD_CTPPS_UFSD_SEGMENT_NAME, DDD_TOTEM_RP_RP_NAME, DDD_TOTEM_RP_SENSOR_NAME, DDD_TOTEM_TIMING_RP_NAME, DDD_TOTEM_TIMING_SENSOR_TMPL, gamEcalExtractorBlocks_cff::detector, DDFilteredView::firstChild(), DDFilteredView::logicalPart(), DDName::name(), dataset::name, DDBase< N, C >::name(), DDFilteredView::nextSibling(), DDFilteredView::parent(), DetGeomDesc::setGeographicalID(), relativeConstraints::station, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by produceIdealGD(), and ~CTPPSGeometryESModule().

 {
   // try to dive into next level
   if ( !fv->firstChild() ) return;
 
   // loop over siblings in the level
   do {
     // create new DetGeomDesc node and add it to the parent's (gd) list
     DetGeomDesc* newGD = new DetGeomDesc( fv );
 
     const std::string name = fv->logicalPart().name().name();
 
     // strip sensors
     if ( name == DDD_TOTEM_RP_SENSOR_NAME ) {
       const std::vector<int>& copy_num = fv->copyNumbers();
       // check size of copy numubers array
       if ( copy_num.size() < 3 )
         throw cms::Exception("DDDTotemRPContruction") << "size of copyNumbers for strip sensor is "
           << copy_num.size() << ". It must be >= 3.";
 
       // extract information
       const unsigned int decRPId = copy_num[copy_num.size() - 3];
       const unsigned int arm = decRPId / 100;
       const unsigned int station = ( decRPId % 100 ) / 10;
       const unsigned int rp = decRPId % 10;
       const unsigned int detector = copy_num[copy_num.size() - 1];
       newGD->setGeographicalID( TotemRPDetId( arm, station, rp, detector ) );
     }
 
     // strip and pixels RPs
     else if ( name == DDD_TOTEM_RP_RP_NAME || name == DDD_CTPPS_PIXELS_RP_NAME ) {
       unsigned int decRPId = fv->copyno();
 
       // check if it is a pixel RP
       if ( decRPId >= 10000 ){
         decRPId = decRPId % 10000;
         const unsigned int armIdx = ( decRPId / 100 ) % 10;
         const unsigned int stIdx = ( decRPId / 10 ) % 10;
         const unsigned int rpIdx = decRPId % 10;
         newGD->setGeographicalID( CTPPSPixelDetId( armIdx, stIdx, rpIdx ) );
       }
       else {
         const unsigned int armIdx = ( decRPId / 100 ) % 10;
         const unsigned int stIdx = ( decRPId / 10 ) % 10;
         const unsigned int rpIdx = decRPId % 10;
         newGD->setGeographicalID( TotemRPDetId( armIdx, stIdx, rpIdx ) );
       }
     }
 
     else if ( std::regex_match( name, std::regex( DDD_TOTEM_TIMING_SENSOR_TMPL ) ) ) {
       const std::vector<int>& copy_num = fv->copyNumbers();
       // check size of copy numbers array
       if ( copy_num.size() < 4 )
         throw cms::Exception("DDDTotemRPContruction") << "size of copyNumbers for TOTEM timing sensor is "
           << copy_num.size() << ". It must be >= 4.";
 
       const unsigned int decRPId = copy_num[copy_num.size()-4];
       const unsigned int arm = decRPId / 100, station = ( decRPId % 100 )/10, rp = decRPId % 10;
       const unsigned int plane = copy_num[copy_num.size()-2], channel = copy_num[copy_num.size()-1];
       newGD->setGeographicalID( TotemTimingDetId( arm, station, rp, plane, channel ) );
     }
 
     else if ( name == DDD_TOTEM_TIMING_RP_NAME ) {
       const unsigned int arm = fv->copyno() / 100, station = ( fv->copyno() % 100 )/10, rp = fv->copyno() % 10;
       newGD->setGeographicalID( TotemTimingDetId( arm, station, rp ) );
     }
 
     // pixel sensors
     else if ( name == DDD_CTPPS_PIXELS_SENSOR_NAME ) {
       const std::vector<int>& copy_num = fv->copyNumbers();
       // check size of copy numubers array
       if ( copy_num.size() < 4 )
         throw cms::Exception("DDDTotemRPContruction") << "size of copyNumbers for pixel sensor is "
           << copy_num.size() << ". It must be >= 4.";
 
       // extract information
       const unsigned int decRPId = copy_num[copy_num.size() - 4] % 10000;
       const unsigned int arm = decRPId / 100;
       const unsigned int station = ( decRPId % 100 ) / 10;
       const unsigned int rp = decRPId % 10;
       const unsigned int detector = copy_num[copy_num.size() - 2] - 1;
       newGD->setGeographicalID( CTPPSPixelDetId( arm, station, rp, detector ) );
     }
 
     // diamond/UFSD sensors
     else if ( name == DDD_CTPPS_DIAMONDS_SEGMENT_NAME || name == DDD_CTPPS_UFSD_SEGMENT_NAME ) {
       const std::vector<int>& copy_num = fv->copyNumbers();
 
       const unsigned int id = copy_num[copy_num.size()-1];
       const unsigned int arm = copy_num[1]-1;
       const unsigned int station = 1;
       const unsigned int rp = 6;
       const unsigned int plane = ( id / 100 );
       const unsigned int channel = id % 100;
 
       newGD->setGeographicalID( CTPPSDiamondDetId( arm, station, rp, plane, channel ) );
     }
 
     // diamond/UFSD RPs
     else if ( name == DDD_CTPPS_DIAMONDS_RP_NAME ) {
       const std::vector<int>& copy_num = fv->copyNumbers();
 
       // check size of copy numubers array
       if ( copy_num.size() < 2 )
         throw cms::Exception("DDDTotemRPContruction") << "size of copyNumbers for diamond RP is "
           << copy_num.size() << ". It must be >= 2.";
 
       const unsigned int arm = copy_num[1] - 1;
       const unsigned int station = 1;
       const unsigned int rp = 6;
 
       newGD->setGeographicalID( CTPPSDiamondDetId( arm, station, rp ) );
     }
 
     // add component
     gd->addComponent( newGD );
 
     // recursion
     buildDetGeomDesc( fv, newGD );
   } while ( fv->nextSibling() );
 
   // go a level up
   fv->parent();
 }

std::unique_ptr< DetGeomDesc > CTPPSGeometryESModule::produceIdealGD ( const IdealGeometryRecord & iRecord )

Definition at line 280 of file CTPPSGeometryESModule.cc.

References buildDetGeomDesc(), compactViewTag_, ALCARECOTkAlBeamHalo_cff::filter, edm::eventsetup::EventSetupRecord::get(), edm::ESHandle< T >::product(), and pyrootRender::root.

Referenced by CTPPSGeometryESModule(), and ~CTPPSGeometryESModule().

 {
   // get the DDCompactView from EventSetup
   edm::ESHandle<DDCompactView> cpv;
   iRecord.get( compactViewTag_, cpv );
 
   // create DDFilteredView and apply the filter
   DDPassAllFilter filter;
   DDFilteredView fv( *( cpv.product() ), filter );
 
   // conversion to DetGeomDesc structure
   DetGeomDesc* root = new DetGeomDesc( &fv );
   buildDetGeomDesc( &fv, root );
 
   // construct the tree of DetGeomDesc
   return std::unique_ptr<DetGeomDesc>( const_cast<DetGeomDesc*>( root ) );
 }

std::unique_ptr< DetGeomDesc > CTPPSGeometryESModule::produceMisalignedGD ( const VeryForwardMisalignedGeometryRecord & iRecord )

Definition at line 333 of file CTPPSGeometryESModule.cc.

References applyAlignments(), edm::eventsetup::DependentRecordImplementation< RecordT, ListT >::getRecord(), RPAlignmentCorrectionsData::getRPMap(), RPAlignmentCorrectionsData::getSensorMap(), edm::ESHandleBase::isValid(), and verbosity_.

Referenced by CTPPSGeometryESModule(), and ~CTPPSGeometryESModule().

 {
   // get the input GeometricalDet
   edm::ESHandle<DetGeomDesc> idealGD;
   iRecord.getRecord<IdealGeometryRecord>().get( idealGD );
 
   // load alignments
   edm::ESHandle<RPAlignmentCorrectionsData> alignments;
   try { iRecord.getRecord<RPMisalignedAlignmentRecord>().get( alignments ); }
   catch ( cms::Exception& ) {}
 
   if ( alignments.isValid() ) {
     if ( verbosity_ )
       edm::LogVerbatim("CTPPSGeometryESModule::produceMisalignedGD")
         << ">> CTPPSGeometryESModule::produceMisalignedGD > Misaligned geometry: "
         << alignments->getRPMap().size() << " RP and "
         << alignments->getSensorMap().size() << " sensor alignments applied.";
   } else {
     if ( verbosity_ )
       edm::LogVerbatim("CTPPSGeometryESModule::produceMisalignedGD")
         << ">> CTPPSGeometryESModule::produceMisalignedGD > Misaligned geometry: No alignments applied.";
   }
 
   DetGeomDesc* newGD = nullptr;
   applyAlignments( idealGD, alignments, newGD );
   return std::unique_ptr<DetGeomDesc>( newGD );
 }

std::unique_ptr< CTPPSGeometry > CTPPSGeometryESModule::produceMisalignedTG ( const VeryForwardMisalignedGeometryRecord & iRecord )

Definition at line 375 of file CTPPSGeometryESModule.cc.

References DEFINE_FWK_EVENTSETUP_MODULE, edm::eventsetup::EventSetupRecord::get(), and edm::ESHandle< T >::product().

Referenced by CTPPSGeometryESModule(), and ~CTPPSGeometryESModule().

 {
   edm::ESHandle<DetGeomDesc> gD;
   iRecord.get( gD );
 
   return std::make_unique<CTPPSGeometry>( gD.product() );
 }

std::unique_ptr< DetGeomDesc > CTPPSGeometryESModule::produceRealGD ( const VeryForwardRealGeometryRecord & iRecord )

Definition at line 301 of file CTPPSGeometryESModule.cc.

References applyAlignments(), edm::eventsetup::DependentRecordImplementation< RecordT, ListT >::getRecord(), RPAlignmentCorrectionsData::getRPMap(), RPAlignmentCorrectionsData::getSensorMap(), edm::ESHandleBase::isValid(), and verbosity_.

Referenced by CTPPSGeometryESModule(), and ~CTPPSGeometryESModule().

 {
   // get the input GeometricalDet
   edm::ESHandle<DetGeomDesc> idealGD;
   iRecord.getRecord<IdealGeometryRecord>().get( idealGD );
 
   // load alignments
   edm::ESHandle<RPAlignmentCorrectionsData> alignments;
   try { iRecord.getRecord<RPRealAlignmentRecord>().get( alignments ); }
   catch ( cms::Exception& ) {}
 
   if ( alignments.isValid() ) {
     if ( verbosity_ )
       edm::LogVerbatim("CTPPSGeometryESModule::produceRealGD")
         << ">> CTPPSGeometryESModule::produceRealGD > Real geometry: "
         << alignments->getRPMap().size() << " RP and "
         << alignments->getSensorMap().size() << " sensor alignments applied.";
   }
   else {
     if ( verbosity_ )
       edm::LogVerbatim("CTPPSGeometryESModule::produceRealGD")
         << ">> CTPPSGeometryESModule::produceRealGD > Real geometry: No alignments applied.";
   }
 
   DetGeomDesc* newGD = nullptr;
   applyAlignments( idealGD, alignments, newGD );
   return std::unique_ptr<DetGeomDesc>( newGD );
 }

std::unique_ptr< CTPPSGeometry > CTPPSGeometryESModule::produceRealTG ( const VeryForwardRealGeometryRecord & iRecord )

Definition at line 364 of file CTPPSGeometryESModule.cc.

References edm::eventsetup::EventSetupRecord::get(), and edm::ESHandle< T >::product().

Referenced by CTPPSGeometryESModule(), and ~CTPPSGeometryESModule().

 {
   edm::ESHandle<DetGeomDesc> gD;
   iRecord.get( gD );
 
   return std::make_unique<CTPPSGeometry>( gD.product() );
 }

Member Data Documentation

std::string CTPPSGeometryESModule::compactViewTag_

protected

Definition at line 69 of file CTPPSGeometryESModule.cc.

Referenced by produceIdealGD().

unsigned int CTPPSGeometryESModule::verbosity_

protected

Definition at line 68 of file CTPPSGeometryESModule.cc.

Referenced by produceMisalignedGD(), and produceRealGD().

Public Member Functions

Static Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation