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

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () 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

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

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
const std::string	geometryType_

const bool	printRPInfo_

const bool	printSensorInfo_

const edm::ESGetToken < CTPPSGeometry, IdealGeometryRecord >	tokIdeal_

const edm::ESGetToken < CTPPSGeometry, VeryForwardMisalignedGeometryRecord >	tokMis_

const edm::ESGetToken < CTPPSGeometry, VeryForwardRealGeometryRecord >	tokReal_

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)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Class to print out information on current geometry.

Definition at line 35 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)),
       tokIdeal_(esConsumes<CTPPSGeometry, IdealGeometryRecord>()),
       tokReal_(esConsumes<CTPPSGeometry, VeryForwardRealGeometryRecord>()),
       tokMis_(esConsumes<CTPPSGeometry, VeryForwardMisalignedGeometryRecord>()) {}

Member Function Documentation

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 69 of file CTPPSGeometryInfo.cc.

References edm::ESWatcher< T >::check(), Exception, relativeConstraints::geometry, geometryType_, edm::EventSetup::getData(), printGeometry(), tokIdeal_, tokMis_, tokReal_, watcherIdealGeometry_, watcherMisalignedGeometry_, and watcherRealGeometry_.

                                                                                    {
   if (geometryType_ == "ideal") {
     if (watcherIdealGeometry_.check(iSetup)) {
       const auto& geometry = iSetup.getData(tokIdeal_);
       printGeometry(geometry, iEvent);
     }
     return;
   }
 
   else if (geometryType_ == "real") {
     if (watcherRealGeometry_.check(iSetup)) {
       const auto& geometry = iSetup.getData(tokReal_);
       printGeometry(geometry, iEvent);
     }
     return;
   }
 
   else if (geometryType_ == "misaligned") {
     if (watcherMisalignedGeometry_.check(iSetup)) {
       const auto& geometry = iSetup.getData(tokMis_);
       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 99 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 142 of file CTPPSGeometryInfo.cc.

References CTPPSGeometry::beginRP(), CTPPSGeometry::beginSensor(), CTPPSGeometry::endRP(), CTPPSGeometry::endSensor(), formatDetId(), geometryType_, CTPPSGeometry::localToGlobal(), printRPInfo_, printSensorInfo_, and submitPVValidationJobs::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 auto gl_o = geometry.localToGlobal(detId, CTPPSGeometry::Vector(0, 0, 0));
       const auto gl_a1 = geometry.localToGlobal(detId, CTPPSGeometry::Vector(1, 0, 0)) - gl_o;
       const auto gl_a2 = geometry.localToGlobal(detId, CTPPSGeometry::Vector(0, 1, 0)) - gl_o;
       const auto gl_a3 = geometry.localToGlobal(detId, CTPPSGeometry::Vector(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();
 }