Class to print out information on current geometry. More...

Inheritance diagram for CTPPSGeometryInfo:

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

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	printGeometry (const CTPPSGeometry &, const edm::Event &)

Static Private Member Functions
static std::string	formatDetId (const CTPPSDetId &id, bool printDetails=true)

Private Attributes
std::string	geometryType_

bool	printRPInfo_

bool	printSensorInfo_

edm::ESWatcher< IdealGeometryRecord >	watcherIdealGeometry_

edm::ESWatcher< VeryForwardMisalignedGeometryRecord >	watcherMisalignedGeometry_

edm::ESWatcher< VeryForwardRealGeometryRecord >	watcherRealGeometry_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Class to print out information on current geometry.

Definition at line 36 of file CTPPSGeometryInfo.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 59 of file CTPPSGeometryInfo.cc.

                                                                      :
   geometryType_   ( iConfig.getUntrackedParameter<std::string>( "geometryType", "real" ) ),
   printRPInfo_    ( iConfig.getUntrackedParameter<bool>( "printRPInfo", true ) ),
   printSensorInfo_( iConfig.getUntrackedParameter<bool>( "printSensorInfo", true ) )
 {}

Member Function Documentation

void CTPPSGeometryInfo::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 68 of file CTPPSGeometryInfo.cc.

References edm::ESWatcher< T >::check(), Exception, geometry, geometryType_, edm::EventSetup::get(), printGeometry(), AlCaHLTBitMon_QueryRunRegistry::string, watcherIdealGeometry_, watcherMisalignedGeometry_, and watcherRealGeometry_.

 {
   edm::ESHandle<CTPPSGeometry> geometry;
 
   if ( geometryType_ == "ideal" ) {
     if ( watcherIdealGeometry_.check( iSetup ) ) {
       iSetup.get<IdealGeometryRecord>().get( geometry );
       printGeometry( *geometry, iEvent );
     }
     return;
   }
 
   else if ( geometryType_ == "real" ) {
     if ( watcherRealGeometry_.check( iSetup ) ) {
       iSetup.get<VeryForwardRealGeometryRecord>().get( geometry );
       printGeometry( *geometry, iEvent );
     }
     return;
   }
 
   else if ( geometryType_ == "misaligned" ) {
     if ( watcherMisalignedGeometry_.check( iSetup ) ) {
       iSetup.get<VeryForwardMisalignedGeometryRecord>().get( geometry );
       printGeometry( *geometry, iEvent );
     }
     return;
   }
 
   throw cms::Exception("CTPPSGeometryInfo") << "Unknown geometry type: `" << geometryType_ << "'.";
 }

std::string CTPPSGeometryInfo::formatDetId	(	const CTPPSDetId &	id,
		bool	printDetails = `true`
	)

staticprivate

Definition at line 102 of file CTPPSGeometryInfo.cc.

References CTPPSDiamondDetId::channel(), TotemTimingDetId::channel(), CTPPSPixelDetId::plane(), CTPPSDiamondDetId::plane(), TotemRPDetId::plane(), TotemTimingDetId::plane(), CTPPSDetId::sdTimingDiamond, CTPPSDetId::sdTimingFastSilicon, CTPPSDetId::sdTrackingPixel, and CTPPSDetId::sdTrackingStrip.

Referenced by printGeometry().

 {
   std::ostringstream oss;
   oss << id.rawId();
 
   const unsigned int rpDecId = id.arm()*100 + id.station()*10 + id.rp();
 
   if ( id.subdetId() == CTPPSDetId::sdTrackingStrip ) {
     TotemRPDetId fid( id );
     oss << " (strip RP " << std::setw(3) << rpDecId;
     if ( printDetails )
       oss <<  ", plane " << fid.plane();
     oss << ")";
   }
 
   else if ( id.subdetId() == CTPPSDetId::sdTrackingPixel ) {
     CTPPSPixelDetId fid( id );
     oss << " (pixel RP " << std::setw(3) << rpDecId;
     if ( printDetails )
       oss <<  ", plane " << fid.plane();
     oss << ")";
   }
 
   else if ( id.subdetId() == CTPPSDetId::sdTimingDiamond ) {
     CTPPSDiamondDetId fid( id );
     oss << " (diamd RP " << std::setw(3) << rpDecId;
     if ( printDetails )
       oss <<  ", plane " << fid.plane() << ", channel " << std::setw(2) << fid.channel();
     oss << ")";
   }
 
   else if ( id.subdetId() == CTPPSDetId::sdTimingFastSilicon ) {
     TotemTimingDetId fid( id );
     oss << " (totim RP " << std::setw(3) << rpDecId;
     if ( printDetails )
       oss <<  ", plane " << fid.plane() << ", channel " << std::setw(2) << fid.channel();
     oss << ")";
   }
 
   return oss.str();
 }

void CTPPSGeometryInfo::printGeometry	(	const CTPPSGeometry &	geometry,
		const edm::Event &	event
	)

private

Definition at line 147 of file CTPPSGeometryInfo.cc.

References CTPPSGeometry::beginRP(), CTPPSGeometry::beginSensor(), DEFINE_FWK_MODULE, CTPPSGeometry::endRP(), CTPPSGeometry::endSensor(), alignBH_cfg::fixed, formatDetId(), geometryType_, CTPPSGeometry::localToGlobal(), printRPInfo_, printSensorInfo_, and protons_cff::t.

Referenced by analyze().

 {
   time_t unixTime = event.time().unixTime();
   char timeStr[50];
   strftime( timeStr, 50, "%F %T", localtime( &unixTime ) );
 
   std::ostringstream oss;
 
   // RP geometry
   if ( printRPInfo_ ) {
     oss << "* RPs:\n"
         << "    ce: RP center in global coordinates, in mm\n";
 
     for ( auto it = geometry.beginRP(); it != geometry.endRP(); ++it ) {
       const DetGeomDesc::Translation &t = it->second->translation();
 
       oss << formatDetId( CTPPSDetId( it->first ), false )
         << std::fixed << std::setprecision( 3 ) << std::showpos
         << " | ce=(" << t.x() << ", " << t.y() << ", " << t.z() << ")\n";
     }
 
     edm::LogVerbatim("CTPPSGeometryInfo") << oss.str();
   }
 
   // sensor geometry
   if ( printSensorInfo_ ) {
     oss << "* sensors:\n"
       << "    ce: sensor center in global coordinates, in mm\n"
       << "    a1: local axis (1, 0, 0) in global coordinates\n"
       << "    a2: local axis (0, 1, 0) in global coordinates\n"
       << "    a3: local axis (0, 0, 1) in global coordinates\n";
 
     for ( auto it = geometry.beginSensor(); it != geometry.endSensor(); ++it ) {
       CTPPSDetId detId( it->first );
 
       const CLHEP::Hep3Vector gl_o  = geometry.localToGlobal( detId, CLHEP::Hep3Vector( 0, 0, 0 ) );
       const CLHEP::Hep3Vector gl_a1 = geometry.localToGlobal( detId, CLHEP::Hep3Vector( 1, 0, 0 ) ) - gl_o;
       const CLHEP::Hep3Vector gl_a2 = geometry.localToGlobal( detId, CLHEP::Hep3Vector( 0, 1, 0 ) ) - gl_o;
       const CLHEP::Hep3Vector gl_a3 = geometry.localToGlobal( detId, CLHEP::Hep3Vector( 0, 0, 1 ) ) - gl_o;
 
       oss << formatDetId( detId )
         << std::fixed << std::setprecision( 3 ) << std::showpos
         << " | ce=(" << gl_o.x() << ", " << gl_o.y() << ", " << gl_o.z() << ")"
         << " | a1=(" << gl_a1.x() << ", " << gl_a1.y() << ", " << gl_a1.z() << ")"
         << " | a2=(" << gl_a2.x() << ", " << gl_a2.y() << ", " << gl_a2.z() << ")"
         << " | a3=(" << gl_a3.x() << ", " << gl_a3.y() << ", " << gl_a3.z() << ")\n";
     }
   }
 
   edm::LogInfo("CTPPSGeometryInfo")
     << "New " << geometryType_ << " geometry found in run="
     << event.id().run() << ", event=" << event.id().event() << ", UNIX timestamp=" << unixTime
     << " (" << timeStr << ")\n"
     << oss.str();
 }