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

Inheritance diagram for GeometryInfoModule:

Classes
struct	MeanRPData

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

	~GeometryInfoModule ()

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

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

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

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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 &	itemsToGetFromEvent () 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

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

virtual	~EDConsumerBase () noexcept(false)

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

virtual void	beginRun (edm::Run const &, edm::EventSetup const &)

virtual void	endJob ()

void	PrintGeometry (const TotemRPGeometry &, const edm::Event &)

Private Attributes
std::string	geometryType

bool	printMeanSensorInfo

bool	printRPInfo

bool	printSensorInfo

edm::ESWatcher < IdealGeometryRecord >	watcherIdealGeometry

edm::ESWatcher < VeryForwardMeasuredGeometryRecord >	watcherMeasuredGeometry

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 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.

See schema of TOTEM RP geometry classes

Definition at line 31 of file GeometryInfoModule.cc.

Constructor & Destructor Documentation

GeometryInfoModule::GeometryInfoModule ( const edm::ParameterSet & ps )

explicit

Definition at line 87 of file GeometryInfoModule.cc.

                                                                 :
   geometryType(ps.getUntrackedParameter<string>("geometryType", "real")),
   printRPInfo(ps.getUntrackedParameter<bool>("printRPInfo", true)),
   printSensorInfo(ps.getUntrackedParameter<bool>("printSensorInfo", true)),
   printMeanSensorInfo(ps.getUntrackedParameter<bool>("printMeanSensorInfo", true))
 {
 }

GeometryInfoModule::~GeometryInfoModule ( )

Definition at line 97 of file GeometryInfoModule.cc.

98 {

99 }

Member Function Documentation

void GeometryInfoModule::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	es
	)

privatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 115 of file GeometryInfoModule.cc.

References edm::ESWatcher< T >::check(), Exception, geometry, geometryType, edm::EventSetup::get(), PrintGeometry(), watcherIdealGeometry, watcherMeasuredGeometry, watcherMisalignedGeometry, and watcherRealGeometry.

 {
   ESHandle<TotemRPGeometry> geometry;
 
   if (!geometryType.compare("ideal"))
   {
     if (watcherIdealGeometry.check(es))
     {
       es.get<IdealGeometryRecord>().get(geometry);
       PrintGeometry(*geometry, event);
     }
     return;
   }
 
   if (!geometryType.compare("measured"))
   {
     if (watcherMeasuredGeometry.check(es))
     {
       es.get<VeryForwardMeasuredGeometryRecord>().get(geometry);
       PrintGeometry(*geometry, event);
     }
     return;
   }
 
   if (!geometryType.compare("real"))
   {
     if (watcherRealGeometry.check(es))
     {
       es.get<VeryForwardRealGeometryRecord>().get(geometry);
       PrintGeometry(*geometry, event);
     }
     return;
   }
 
   if (!geometryType.compare("misaligned"))
   {
     if (watcherMisalignedGeometry.check(es))
     {
       es.get<VeryForwardMisalignedGeometryRecord>().get(geometry);
       PrintGeometry(*geometry, event);
     }
     return;
   }
 
   throw cms::Exception("GeometryInfoModule") << "Unknown geometry type: `" << geometryType << "'.";
 }

void GeometryInfoModule::beginRun	(	edm::Run const &	,
		edm::EventSetup const &	iSetup
	)

privatevirtual

Definition at line 103 of file GeometryInfoModule.cc.

104 {

105 }

void GeometryInfoModule::endJob ( void )

privatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 109 of file GeometryInfoModule.cc.

110 {

111 }

void GeometryInfoModule::PrintGeometry	(	const TotemRPGeometry &	geometry,
		const edm::Event &	event
	)

private

Definition at line 164 of file GeometryInfoModule.cc.

Referenced by analyze().

 {
   time_t unixTime = event.time().unixTime();
   char timeStr[50];
   strftime(timeStr, 50, "%F %T", localtime(&unixTime));
   printf(">> GeometryInfoModule::PrintGeometry\n\tnew %s geometry found in run=%u, event=%llu, UNIX timestamp=%lu (%s)\n",
     geometryType.c_str(), event.id().run(), event.id().event(), unixTime, timeStr);
 
   // RP geometry
   if (printRPInfo)
   {
     printf("\n* RPs:\n");
     printf("  ID |   x (mm)    |   y (mm)    |   z  (m)  |\n");
     for (TotemRPGeometry::RPDeviceMapType::const_iterator it = geometry.beginRP(); it != geometry.endRP(); ++it)
     {
       const DDTranslation &t = it->second->translation();
       printf(" %3i | %+11.3f | %+11.3f | %+9.4f |\n", it->first, t.x(), t.y(), t.z() * 1E-3);
     }
   }
   
   // sensor geometry
   if (printSensorInfo)
   {
     printf("\n* silicon sensors:\n");
     printf("DetId |            detector center           |  readout direction  |\n");
     printf("      |   x (mm)   |   y (mm)   |   z  (m)   |    dx    |    dy    |\n");
     for (TotemRPGeometry::mapType::const_iterator it = geometry.beginDet(); it != geometry.endDet(); ++it)
     {
       TotemRPDetId id(it->first);
   
       double x = it->second->translation().x();
       double y = it->second->translation().y();
       double z = it->second->translation().z() * 1E-3;
   
       double dx = 0., dy = 0.;
       geometry.GetReadoutDirection(it->first, dx, dy);
   
       printf("%4i  |  %8.3f  |  %8.3f  |  %9.4f | %8.3f | %8.3f |\n", id.getPlaneDecimalId(), x, y, z, dx, dy);
     }
   }
   
   // sensor geometry averaged over 1 RP
   if (printMeanSensorInfo)
   {
     printf("\n* average over silicon sensors of 1 RP:\n");
 
     map<unsigned int, MeanRPData> data; // map: RP Id --> sums of geometrical info
 
     for (TotemRPGeometry::mapType::const_iterator it = geometry.beginDet(); it != geometry.endDet(); ++it)
     {
       TotemRPDetId plId(it->first);
       unsigned int rpDecId = plId.getRPDecimalId();
       bool uDirection = plId.isStripsCoordinateUDirection();
   
       double x = it->second->translation().x();
       double y = it->second->translation().y();
       double z = it->second->translation().z() * 1E-3;
   
       double dx = 0., dy = 0.;
       geometry.GetReadoutDirection(it->first, dx, dy);
   
       data[rpDecId].Fill(x, y, z, uDirection, dx, dy);
     }
 
     printf("RPId |                center                |      U direction    |      V direction    |\n");
     printf("     |   x (mm)   |   y (mm)   |   z  (m)   |    dx    |    dy    |    dx    |    dy    |\n");
     
     for (map<unsigned int, MeanRPData>::iterator it = data.begin(); it != data.end(); ++it)
     {
       const MeanRPData &d = it->second;
   
       double mx = (d.N > 0) ? d.Sx / d.N : 0.;
       double my = (d.N > 0) ? d.Sy / d.N : 0.;
       double mz = (d.N > 0) ? d.Sz / d.N : 0.;
 
       double mdx_u = (d.N_u > 0) ? d.Sdx_u / d.N_u : 0.;
       double mdy_u = (d.N_u > 0) ? d.Sdy_u / d.N_u : 0.;
 
       double mdx_v = (d.N_v > 0) ? d.Sdx_v / d.N_v : 0.;
       double mdy_v = (d.N_v > 0) ? d.Sdy_v / d.N_v : 0.;
 
       
       printf(" %3i |  %8.3f  |  %8.3f  |  %9.4f | %8.3f | %8.3f | %8.3f | %8.3f |\n", it->first, mx, my, mz, mdx_u, mdy_u, mdx_v, mdy_v);
     }
   }
 }