mtd::MTDGeomUtil Class Reference

#include <MTDGeomUtil.h>

Public Member Functions
int	crystalInModule (const DetId &) const
float	eta (const GlobalPoint &position, const float &vertex_z=0.) const
float	eta (const DetId &id, const LocalPoint &local_point, const float &vertex_z=0.) const
const MTDGeometry *	geometry () const
GlobalPoint	globalPosition (const DetId &id, const LocalPoint &local_point) const
bool	isBTL (const DetId &) const
bool	isETL (const DetId &) const
unsigned int	layer (const DetId &) const
int	module (const DetId &) const
	MTDGeomUtil ()
float	phi (const GlobalPoint &position) const
float	phi (const DetId &id, const LocalPoint &local_point) const
std::pair< float, float >	pixelInModule (const DetId &id, const int row, const int column) const
std::pair< uint8_t, uint8_t >	pixelInModule (const DetId &id, const LocalPoint &local_point) const
std::pair< LocalPoint, GlobalPoint >	position (const DetId &id, int row=0, int column=0) const
float	pt (const GlobalPoint &position, const float &hitEnergy, const float &vertex_z=0.) const
float	pt (const DetId &id, const LocalPoint &local_point, const float &hitEnergy, const float &vertex_z=0.) const
uint32_t	sensorModuleId (const DetId &id) const
void	setGeometry (MTDGeometry const *geom)
void	setTopology (MTDTopology const *topo)
const MTDTopology *	topology () const
int	zside (const DetId &id) const
	~MTDGeomUtil ()

Private Attributes
const MTDGeometry *	geom_
const MTDTopology *	topology_

Detailed Description

Definition at line 23 of file MTDGeomUtil.h.

Constructor & Destructor Documentation

MTDGeomUtil()

mtd::MTDGeomUtil::MTDGeomUtil ( )

inline

Definition at line 25 of file MTDGeomUtil.h.

: geom_(nullptr), topology_(nullptr) {}

~MTDGeomUtil()

mtd::MTDGeomUtil::~MTDGeomUtil ( )

inline

Definition at line 26 of file MTDGeomUtil.h.

{}

Member Function Documentation

crystalInModule()

int MTDGeomUtil::crystalInModule

(

const DetId &

id

)

const

Definition at line 159 of file MTDGeomUtil.cc.

References BTLDetId::crystal().

                                                      {
  BTLDetId hid(id);
  return hid.crystal();
}

eta() [1/2]

float mtd::MTDGeomUtil::eta	(	const GlobalPoint &	position,
		const float &	vertex_z = 0.
	)		const

Referenced by Particle.Particle::__str__(), datamodel.Object::DeltaR(), Jet.Jet::jetID(), datamodel.Object::p4(), and Jet.Jet::puJetId().

eta() [2/2]

float mtd::MTDGeomUtil::eta	(	const DetId &	id,
		const LocalPoint &	local_point,
		const float &	vertex_z = 0.
	)		const

Referenced by Particle.Particle::__str__(), datamodel.Object::DeltaR(), Jet.Jet::jetID(), datamodel.Object::p4(), and Jet.Jet::puJetId().

geometry()

const MTDGeometry* mtd::MTDGeomUtil::geometry ( ) const

inline

Definition at line 58 of file MTDGeomUtil.h.

References geom_.

{ return geom_; };

globalPosition()

GlobalPoint MTDGeomUtil::globalPosition	(	const DetId &	id,
		const LocalPoint &	local_point
	)		const

Definition at line 58 of file MTDGeomUtil.cc.

References MTDTopologyMode::crysLayoutFromTopoMode(), TauDecayModes::dec, hcalRecHitTable_cff::detId, Exception, geom_, MTDTopology::getMTDTopologyMode(), MTDGeometry::idToDet(), isBTL(), RectangularMTDTopology::nrows(), RectangularMTDTopology::pixelToModuleLocalPoint(), DetId::rawId(), ProxyMTDTopology::specificTopology(), GeomDet::toGlobal(), GeomDet::topology(), and topology_.

Referenced by MtdTruthAccumulator::fillSimHits(), and phi().

                                                                                            {
  auto global_point = GlobalPoint(0., 0., 0.);
  if (isBTL(id)) {
    BTLDetId detId{id};
    DetId geoId = detId.geographicalId(MTDTopologyMode::crysLayoutFromTopoMode(topology_->getMTDTopologyMode()));
    const MTDGeomDet* thedet = geom_->idToDet(geoId);
    if (thedet == nullptr)
      throw cms::Exception("mtd::MTDGeomUtil") << "GeographicalID: " << std::hex << geoId.rawId() << " ("
                                               << detId.rawId() << ") is invalid!" << std::dec << std::endl;
    const ProxyMTDTopology& topoproxy = static_cast<const ProxyMTDTopology&>(thedet->topology());
    const RectangularMTDTopology& topo = static_cast<const RectangularMTDTopology&>(topoproxy.specificTopology());
    auto local_point_sim =
        topo.pixelToModuleLocalPoint(local_point, detId.row(topo.nrows()), detId.column(topo.nrows()));
    global_point = thedet->toGlobal(local_point_sim);
  } else {
    ETLDetId detId{id};
    DetId geoId = detId.geographicalId();
    const MTDGeomDet* thedet = geom_->idToDet(geoId);
    if (thedet == nullptr)
      throw cms::Exception("mtd::MTDGeomUtil") << "GeographicalID: " << std::hex << geoId.rawId() << " ("
                                               << detId.rawId() << ") is invalid!" << std::dec << std::endl;
    global_point = thedet->toGlobal(local_point);
  }
  return global_point;
}

isBTL()

bool MTDGeomUtil::isBTL

(

const DetId &

id

)

const

Definition at line 24 of file MTDGeomUtil.cc.

References isETL().

Referenced by MtdTruthAccumulator::finalizeEvent(), globalPosition(), pixelInModule(), position(), and sensorModuleId().

{ return !(isETL(id)); }

isETL()

bool MTDGeomUtil::isETL

(

const DetId &

id

)

const

Definition at line 14 of file MTDGeomUtil.cc.

References Exception, and MTDDetId::mtdSubDetector().

Referenced by MtdTruthAccumulator::finalizeEvent(), isBTL(), layer(), module(), and pixelInModule().

                                             {
  MTDDetId hid{id};
  const auto& subDet = hid.mtdSubDetector();
  if (subDet == 0)
    throw cms::Exception("mtd::MTDGeomUtil") << "DetId " << hid.rawId() << " not in MTD!" << std::endl;
  if (subDet == MTDDetId::MTDType::ETL)
    return true;
  return false;
}

layer()

unsigned int MTDGeomUtil::layer

(

const DetId &

id

)

const

Definition at line 89 of file MTDGeomUtil.cc.

References isETL(), and ETLDetId::nDisc().

Referenced by MtdTruthAccumulator::fillSimHits(), MtdTruthAccumulator::finalizeEvent(), geometryXMLparser.DTAlignable::index(), and geometryXMLparser.CSCAlignable::index().

                                                     {
  unsigned int layer(0);
  if (isETL(id)) {
    ETLDetId hid(id);
    layer = hid.nDisc();
  }
  return layer;
}

module()

int MTDGeomUtil::module

(

const DetId &

id

)

const

Definition at line 98 of file MTDGeomUtil.cc.

References isETL(), BTLDetId::module(), and ETLDetId::module().

Referenced by MtdTruthAccumulator::fillSimHits(), and MtdTruthAccumulator::finalizeEvent().

                                             {
  int module = -1;
  if (isETL(id)) {
    ETLDetId hid(id);
    module = hid.module();
  } else {
    BTLDetId hid(id);
    module = hid.module();
  }
  return module;
}

phi() [1/2]

float MTDGeomUtil::phi

(

const GlobalPoint &

position

)

const

Definition at line 188 of file MTDGeomUtil.cc.

References Matriplex::atan2(), and position().

Referenced by Particle.Particle::__str__(), datamodel.Object::DeltaR(), datamodel.Object::p4(), phi(), and ntupleDataFormat.Track::phiPull().

                                                        {
  float phi = (position.x() == 0 && position.y() == 0) ? 0 : atan2(position.y(), position.x());
  return phi;
}

phi() [2/2]

float MTDGeomUtil::phi	(	const DetId &	id,
		const LocalPoint &	local_point
	)		const

Definition at line 193 of file MTDGeomUtil.cc.

References Matriplex::atan2(), globalPosition(), phi(), and position().

Referenced by Particle.Particle::__str__(), datamodel.Object::DeltaR(), datamodel.Object::p4(), and ntupleDataFormat.Track::phiPull().

                                                                           {
  GlobalPoint position = globalPosition(id, local_point);
  float phi = (position.x() == 0 && position.y() == 0) ? 0 : atan2(position.y(), position.x());
  return phi;
}

pixelInModule() [1/2]

std::pair< float, float > MTDGeomUtil::pixelInModule	(	const DetId &	id,
		const int	row,
		const int	column
	)		const

Definition at line 111 of file MTDGeomUtil.cc.

References TauDecayModes::dec, hcalRecHitTable_cff::detId, Exception, geom_, MTDGeometry::idToDet(), isBTL(), RectangularMTDTopology::localX(), RectangularMTDTopology::localY(), RectangularMTDTopology::pixel(), DetId::rawId(), ProxyMTDTopology::specificTopology(), and GeomDet::topology().

Referenced by EtlSimHitsValidation::analyze(), EtlLocalRecoValidation::analyze(), and MtdTruthAccumulator::fillSimHits().

                                                                                                       {
  if (isBTL(id))
    throw cms::Exception("mtd::MTDGeomUtil")
        << "ID: " << std::hex << id.rawId() << " from BTL. This method is for ETL only." << std::endl;
  ETLDetId detId(id);
  DetId geoId = detId.geographicalId();
  const MTDGeomDet* thedet = geom_->idToDet(geoId);
  if (thedet == nullptr)
    throw cms::Exception("mtd::MTDGeomUtil") << "GeographicalID: " << std::hex << geoId.rawId() << " (" << detId.rawId()
                                             << ") is invalid!" << std::dec << std::endl;
  const ProxyMTDTopology& topoproxy = static_cast<const ProxyMTDTopology&>(thedet->topology());
  const RectangularMTDTopology& topo = static_cast<const RectangularMTDTopology&>(topoproxy.specificTopology());
  const Local3DPoint local_point(topo.localX(row), topo.localY(column), 0.);
  return topo.pixel(local_point);
}

pixelInModule() [2/2]

std::pair< uint8_t, uint8_t > MTDGeomUtil::pixelInModule	(	const DetId &	id,
		const LocalPoint &	local_point
	)		const

Definition at line 128 of file MTDGeomUtil.cc.

References cuy::col, MTDTopologyMode::crysLayoutFromTopoMode(), TauDecayModes::dec, hcalRecHitTable_cff::detId, Exception, geom_, MTDTopology::getMTDTopologyMode(), MTDGeometry::idToDet(), isETL(), RectangularMTDTopology::nrows(), RectangularMTDTopology::pixelIndex(), RectangularMTDTopology::pixelToModuleLocalPoint(), DetId::rawId(), ProxyMTDTopology::specificTopology(), GeomDet::topology(), and topology_.

                                                                                                         {
  if (isETL(id)) {
    ETLDetId detId(id);
    DetId geoId = detId.geographicalId();
    const MTDGeomDet* thedet = geom_->idToDet(geoId);
    if (thedet == nullptr)
      throw cms::Exception("mtd::MTDGeomUtil") << "GeographicalID: " << std::hex << geoId.rawId() << " ("
                                               << detId.rawId() << ") is invalid!" << std::dec << std::endl;
    const ProxyMTDTopology& topoproxy = static_cast<const ProxyMTDTopology&>(thedet->topology());
    const RectangularMTDTopology& topo = static_cast<const RectangularMTDTopology&>(topoproxy.specificTopology());
    const auto& thepixel = topo.pixelIndex(local_point);
    uint8_t row = static_cast<uint8_t>(thepixel.first);
    uint8_t col = static_cast<uint8_t>(thepixel.second);
    return std::pair<uint8_t, uint8_t>(row, col);
  } else {
    BTLDetId detId(id);
    DetId geoId = detId.geographicalId(MTDTopologyMode::crysLayoutFromTopoMode(topology_->getMTDTopologyMode()));
    const MTDGeomDet* thedet = geom_->idToDet(geoId);
    if (thedet == nullptr)
      throw cms::Exception("mtd::MTDGeomUtil") << "GeographicalID: " << std::hex << geoId.rawId() << " ("
                                               << detId.rawId() << ") is invalid!" << std::dec << std::endl;
    const ProxyMTDTopology& topoproxy = static_cast<const ProxyMTDTopology&>(thedet->topology());
    const RectangularMTDTopology& topo = static_cast<const RectangularMTDTopology&>(topoproxy.specificTopology());
    auto topo_point = topo.pixelToModuleLocalPoint(local_point, detId.row(topo.nrows()), detId.column(topo.nrows()));
    const auto& thepixel = topo.pixelIndex(topo_point);
    uint8_t row = static_cast<uint8_t>(thepixel.first);
    uint8_t col = static_cast<uint8_t>(thepixel.second);
    return std::pair<uint8_t, uint8_t>(row, col);
  }
}

position()

std::pair< LocalPoint, GlobalPoint > MTDGeomUtil::position	(	const DetId &	id,
		int	row = 0,
		int	column = 0
	)		const

Definition at line 27 of file MTDGeomUtil.cc.

References MTDTopologyMode::crysLayoutFromTopoMode(), TauDecayModes::dec, hcalRecHitTable_cff::detId, Exception, geom_, MTDTopology::getMTDTopologyMode(), MTDGeometry::idToDet(), isBTL(), RectangularMTDTopology::localX(), RectangularMTDTopology::localY(), RectangularMTDTopology::nrows(), RectangularMTDTopology::pixelToModuleLocalPoint(), DetId::rawId(), ProxyMTDTopology::specificTopology(), GeomDet::toGlobal(), GeomDet::topology(), and topology_.

Referenced by MtdTruthAccumulator::finalizeEvent(), and phi().

                                                                                                 {
  LocalPoint local_point(0., 0., 0.);
  GlobalPoint global_point(0., 0., 0.);
  if (isBTL(id)) {
    BTLDetId detId{id};
    DetId geoId = detId.geographicalId(MTDTopologyMode::crysLayoutFromTopoMode(topology_->getMTDTopologyMode()));
    const MTDGeomDet* thedet = geom_->idToDet(geoId);
    if (thedet == nullptr)
      throw cms::Exception("mtd::MTDGeomUtil") << "GeographicalID: " << std::hex << geoId.rawId() << " ("
                                               << detId.rawId() << ") is invalid!" << std::dec << std::endl;
    const ProxyMTDTopology& topoproxy = static_cast<const ProxyMTDTopology&>(thedet->topology());
    const RectangularMTDTopology& topo = static_cast<const RectangularMTDTopology&>(topoproxy.specificTopology());

    local_point = topo.pixelToModuleLocalPoint(local_point, detId.row(topo.nrows()), detId.column(topo.nrows()));
    global_point = thedet->toGlobal(local_point);
  } else {
    ETLDetId detId{id};
    DetId geoId = detId.geographicalId();
    const MTDGeomDet* thedet = geom_->idToDet(geoId);
    if (thedet == nullptr)
      throw cms::Exception("mtd::MTDGeomUtil") << "GeographicalID: " << std::hex << geoId.rawId() << " ("
                                               << detId.rawId() << ") is invalid!" << std::dec << std::endl;
    const ProxyMTDTopology& topoproxy = static_cast<const ProxyMTDTopology&>(thedet->topology());
    const RectangularMTDTopology& topo = static_cast<const RectangularMTDTopology&>(topoproxy.specificTopology());

    local_point = LocalPoint(topo.localX(row), topo.localY(column), 0.);
    global_point = thedet->toGlobal(local_point);
  }
  return {local_point, global_point};
}

pt() [1/2]

float mtd::MTDGeomUtil::pt	(	const GlobalPoint &	position,
		const float &	hitEnergy,
		const float &	vertex_z = 0.
	)		const

Referenced by Particle.Particle::__str__(), DiObject.DiMuon::__str__(), Photon.Photon::calScaledIsoValueExp(), Photon.Photon::calScaledIsoValueLin(), Photon.Photon::calScaledIsoValueQuadr(), Electron.Electron::mvaIDLoose(), Electron.Electron::mvaIDRun2(), Electron.Electron::mvaIDTight(), datamodel.Object::p4(), Electron.Electron::ptErr(), ntupleDataFormat.Track::ptPull(), Lepton.Lepton::relIso(), Lepton.Lepton::relIsoFromEA(), Lepton.Lepton::relIsoR(), and Jet.Jet::setCorrP4().

pt() [2/2]

float mtd::MTDGeomUtil::pt	(	const DetId &	id,
		const LocalPoint &	local_point,
		const float &	hitEnergy,
		const float &	vertex_z = 0.
	)		const

Referenced by Particle.Particle::__str__(), DiObject.DiMuon::__str__(), Photon.Photon::calScaledIsoValueExp(), Photon.Photon::calScaledIsoValueLin(), Photon.Photon::calScaledIsoValueQuadr(), Electron.Electron::mvaIDLoose(), Electron.Electron::mvaIDRun2(), Electron.Electron::mvaIDTight(), datamodel.Object::p4(), Electron.Electron::ptErr(), ntupleDataFormat.Track::ptPull(), Lepton.Lepton::relIso(), Lepton.Lepton::relIsoFromEA(), Lepton.Lepton::relIsoR(), and Jet.Jet::setCorrP4().

sensorModuleId()

uint32_t MTDGeomUtil::sensorModuleId

(

const DetId &

id

)

const

Definition at line 165 of file MTDGeomUtil.cc.

References MTDTopologyMode::crysLayoutFromTopoMode(), hcalRecHitTable_cff::detId, MTDTopology::getMTDTopologyMode(), isBTL(), DetId::rawId(), and topology_.

Referenced by MtdRecoClusterToSimLayerClusterAssociatorByHitsImpl::associateRecoToSim(), and MtdRecoClusterToSimLayerClusterAssociatorByHitsImpl::associateSimToReco().

                                                          {
  if (isBTL(id)) {
    BTLDetId detId(id);
    DetId geoId = detId.geographicalId(MTDTopologyMode::crysLayoutFromTopoMode(topology_->getMTDTopologyMode()));
    return (geoId.rawId());
  } else {
    ETLDetId detId(id);
    DetId geoId = detId.geographicalId();
    return (geoId.rawId());
  }
}

setGeometry()

void MTDGeomUtil::setGeometry

(

MTDGeometry const *

geom

)

Definition at line 10 of file MTDGeomUtil.cc.

References relativeConstraints::geom, and geom_.

Referenced by EtlSimHitsValidation::analyze(), EtlLocalRecoValidation::analyze(), MtdTruthAccumulator::initializeEvent(), and MtdRecoClusterToSimLayerClusterAssociatorByHitsProducer::produce().

{ geom_ = geom; }

setTopology()

void MTDGeomUtil::setTopology

(

MTDTopology const *

topo

)

Definition at line 12 of file MTDGeomUtil.cc.

References topology_.

Referenced by MtdTruthAccumulator::initializeEvent(), and MtdRecoClusterToSimLayerClusterAssociatorByHitsProducer::produce().

{ topology_ = topo; }

topology()

const MTDTopology* mtd::MTDGeomUtil::topology ( ) const

inline

Definition at line 59 of file MTDGeomUtil.h.

References topology_.

{ return topology_; };

zside()

int MTDGeomUtil::zside

(

const DetId &

id

)

const

Definition at line 84 of file MTDGeomUtil.cc.

References MTDDetId::zside().

                                            {
  const MTDDetId hid(id);
  return hid.zside();
}

Member Data Documentation

geom_

const MTDGeometry* mtd::MTDGeomUtil::geom_

private

Definition at line 59 of file MTDGeomUtil.h.

Referenced by geometry(), globalPosition(), pixelInModule(), position(), and setGeometry().

topology_

const MTDTopology* mtd::MTDGeomUtil::topology_

private

Definition at line 63 of file MTDGeomUtil.h.

Referenced by globalPosition(), pixelInModule(), position(), sensorModuleId(), setTopology(), and topology().