hgcal::RecHitTools Class Reference

#include <RecHitTools.h>

Public Member Functions
unsigned int	getCell (const DetId &) const
float	getEta (const GlobalPoint &position, const float &vertex_z=0.) const
float	getEta (const DetId &id, const float &vertex_z=0.) const
void	getEvent (const edm::Event &)
void	getEventSetup (const edm::EventSetup &)
unsigned int	getLayer (const DetId &) const
unsigned int	getLayerWithOffset (const DetId &) const
float	getPhi (const GlobalPoint &position) const
float	getPhi (const DetId &id) const
GlobalPoint	getPosition (const DetId &id) const
float	getPt (const GlobalPoint &position, const float &hitEnergy, const float &vertex_z=0.) const
float	getPt (const DetId &id, const float &hitEnergy, const float &vertex_z=0.) const
std::float_t	getRadiusToSide (const DetId &) const
std::float_t	getSiThickness (const DetId &) const
unsigned int	getWafer (const DetId &) const
bool	isHalfCell (const DetId &) const
	RecHitTools ()
	~RecHitTools ()

Private Attributes
const CaloGeometry *	geom_

Detailed Description

Definition at line 17 of file RecHitTools.h.

Constructor & Destructor Documentation

hgcal::RecHitTools::RecHitTools ( )

inline

Definition at line 19 of file RecHitTools.h.

: geom_(nullptr) {}

hgcal::RecHitTools::~RecHitTools ( )

inline

Definition at line 20 of file RecHitTools.h.

{}

Member Function Documentation

unsigned int RecHitTools::getCell

(

const DetId &

id

)

const

Definition at line 147 of file RecHitTools.cc.

References HGCalDetId::cell(), DetId::Forward, and bookConverter::max.

                                                       {
  if( id.det() != DetId::Forward ) {
    edm::LogError("getCell::InvalidSiliconDetid")
      << "det id: " << id.rawId() << " is not HGCal silicon!";
    return std::numeric_limits<unsigned int>::max();
  }
  const HGCalDetId hid(id);
  unsigned int cell = hid.cell();
  return cell;
}

float RecHitTools::getEta	(	const GlobalPoint &	position,
		const float &	vertex_z = 0.
	)		const

Definition at line 170 of file RecHitTools.cc.

References PV3DBase< T, PVType, FrameType >::eta(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by getEta(), and getPt().

                                                                                  {
  GlobalPoint corrected_position = GlobalPoint(position.x(), position.y(), position.z()-vertex_z);
  return corrected_position.eta();
}

float RecHitTools::getEta	(	const DetId &	id,
		const float &	vertex_z = 0.
	)		const

Definition at line 175 of file RecHitTools.cc.

References eta, getEta(), getPosition(), and position.

                                                                      {
  GlobalPoint position = getPosition(id);
  float eta = getEta(position, vertex_z);
  return eta;
}

void RecHitTools::getEvent

(

const edm::Event &

ev

)

Definition at line 63 of file RecHitTools.cc.

Referenced by hgcal::ClusterTools::getEvent().

                                             {
}

void RecHitTools::getEventSetup

(

const edm::EventSetup &

es

)

Definition at line 66 of file RecHitTools.cc.

References relativeConstraints::geom, geom_, edm::EventSetup::get(), and edm::ESHandle< class >::product().

Referenced by hgcal::ClusterTools::getEventSetup(), and HGCalImagingAlgo::getEventSetup().

                                                       {
  edm::ESHandle<CaloGeometry> geom;
  es.get<CaloGeometryRecord>().get(geom);

  geom_ = geom.product();
}

unsigned int RecHitTools::getLayer

(

const DetId &

id

)

const

Definition at line 114 of file RecHitTools.cc.

References HcalDetId::depth(), DetId::Forward, DetId::Hcal, HcalEndcap, HGCalDetId::layer(), and bookConverter::max.

Referenced by getLayerWithOffset().

                                                        {
  unsigned int layer = std::numeric_limits<unsigned int>::max();
  if( id.det() == DetId::Forward) {
    const HGCalDetId hid(id);
    layer = hid.layer();
  } else if( id.det() == DetId::Hcal && id.subdetId() == HcalEndcap) {
    const HcalDetId hcid(id);
    layer = hcid.depth();
  }
  return layer;
}

unsigned int RecHitTools::getLayerWithOffset

(

const DetId &

id

)

const

Definition at line 126 of file RecHitTools.cc.

References DetId::Forward, getLayer(), DetId::Hcal, HcalEndcap, and HGCHEF.

Referenced by HGCalImagingAlgo::makeClusters().

                                                                  {
  unsigned int layer = getLayer(id);  
  if( id.det() == DetId::Forward && id.subdetId() == HGCHEF ) {
    layer += fhOffset;
  } else if( id.det() == DetId::Hcal && id.subdetId() == HcalEndcap) {
    layer += bhOffset;
  }  
  return layer;
}

float RecHitTools::getPhi

(

const GlobalPoint &

position

)

const

Definition at line 181 of file RecHitTools.cc.

References phi(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

                                                           {
  float phi = atan2(position.y(),position.x());
  return phi;
}

float RecHitTools::getPhi

(

const DetId &

id

)

const

Definition at line 186 of file RecHitTools.cc.

References getPosition(), phi(), position, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

                                               {
  GlobalPoint position = getPosition(id);
  float phi = atan2(position.y(),position.x());
  return phi;
}

GlobalPoint RecHitTools::getPosition

(

const DetId &

id

)

const

Definition at line 73 of file RecHitTools.cc.

References relativeConstraints::geom, geom_, CaloSubdetectorGeometry::getGeometry(), HGCalGeometry::getPosition(), CaloCellGeometry::getPosition(), CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, and position.

Referenced by getEta(), getPhi(), getPt(), HGCalImagingAlgo::Hexel::Hexel(), and HGCalImagingAlgo::makeClusters().

                                                          {
  auto geom = geom_->getSubdetectorGeometry(id);
  check_geom(geom);
  GlobalPoint position;
  if( id.det() == DetId::Hcal ) {
    position = geom->getGeometry(id)->getPosition();
  } else {
    const auto* hg = static_cast<const HGCalGeometry*>(geom);
    position = hg->getPosition(id);
  }
  return position;
}

float RecHitTools::getPt	(	const GlobalPoint &	position,
		const float &	hitEnergy,
		const float &	vertex_z = 0.
	)		const

Definition at line 192 of file RecHitTools.cc.

References eta, getEta(), and EnergyCorrector::pt.

                                                                                                         {
  float eta = getEta(position, vertex_z);
  float pt = hitEnergy / cosh(eta);
  return pt;
}

float RecHitTools::getPt	(	const DetId &	id,
		const float &	hitEnergy,
		const float &	vertex_z = 0.
	)		const

Definition at line 198 of file RecHitTools.cc.

References eta, getEta(), getPosition(), position, and EnergyCorrector::pt.

                                                                                             {
  GlobalPoint position = getPosition(id);
  float eta = getEta(position, vertex_z);
  float pt = hitEnergy / cosh(eta);
  return pt;
}

std::float_t RecHitTools::getRadiusToSide

(

const DetId &

id

)

const

Definition at line 100 of file RecHitTools.cc.

References HGCalDDDConstants::cellSizeHex(), DetId::Forward, geom_, CaloGeometry::getSubdetectorGeometry(), bookConverter::max, findQualityFiles::size, and HGCalDetId::waferType().

                                                             {
  auto geom = geom_->getSubdetectorGeometry(id);
  check_geom(geom);
  if( id.det() != DetId::Forward ) {
    edm::LogError("getRadiusToSide::InvalidSiliconDetid")
      << "det id: " << id.rawId() << " is not HGCal silicon!";
    return std::numeric_limits<std::float_t>::max();
  }
  const HGCalDetId hid(id);
  auto ddd = get_ddd(geom,hid);
  std::float_t size = ddd->cellSizeHex(hid.waferType());
  return size;
}

std::float_t RecHitTools::getSiThickness

(

const DetId &

id

)

const

Definition at line 86 of file RecHitTools.cc.

References DetId::Forward, geom_, CaloGeometry::getSubdetectorGeometry(), HGCalDetId::wafer(), and HGCalDDDConstants::waferTypeL().

                                                            {
  auto geom = geom_->getSubdetectorGeometry(id);
  check_geom(geom);
  if( id.det() != DetId::Forward ) {
    edm::LogError("getSiThickness::InvalidSiliconDetid")
      << "det id: " << id.rawId() << " is not HGCal silicon!";
  }  
  const HGCalDetId hid(id);
  auto ddd = get_ddd(geom,hid);
  unsigned int wafer = hid.wafer();
  int tidx = ddd->waferTypeL(wafer);
  return idx_to_thickness*tidx;
}

unsigned int RecHitTools::getWafer

(

const DetId &

id

)

const

Definition at line 136 of file RecHitTools.cc.

References DetId::Forward, bookConverter::max, and HGCalDetId::wafer().

                                                        {
  if( id.det() != DetId::Forward ) {
    edm::LogError("getWafer::InvalidSiliconDetid")
      << "det id: " << id.rawId() << " is not HGCal silicon!";
    return std::numeric_limits<unsigned int>::max();
  } 
  const HGCalDetId hid(id);
  unsigned int wafer = hid.wafer();
  return wafer;
}

bool RecHitTools::isHalfCell

(

const DetId &

id

)

const

Definition at line 158 of file RecHitTools.cc.

References HGCalDetId::cell(), DetId::Forward, geom_, CaloGeometry::getSubdetectorGeometry(), HGCalDDDConstants::isHalfCell(), HGCalDetId::waferType(), and HGCalDDDConstants::waferTypeT().

Referenced by HGCalImagingAlgo::makeClusters().

                                                  {
  if( id.det() != DetId::Forward ) {
    return false;
  }
  auto geom = geom_->getSubdetectorGeometry(id);
  check_geom(geom);  
  const HGCalDetId hid(id);
  auto ddd = get_ddd(geom,hid);
  const int waferType = ddd->waferTypeT(hid.waferType());
  return ddd->isHalfCell(waferType,hid.cell());
}

Member Data Documentation

const CaloGeometry* hgcal::RecHitTools::geom_

private

Definition at line 48 of file RecHitTools.h.

Referenced by getEventSetup(), getPosition(), getRadiusToSide(), getSiThickness(), and isHalfCell().