hgcal::RecHitTools Class Reference

#include <RecHitTools.h>

Public Member Functions
std::pair< uint32_t, uint32_t >	firstAndLastLayer (DetId::Detector det, int subdet) const
unsigned int	firstLayerBH () const
std::pair< int, 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
const CaloGeometry *	getGeometry () const
int	getGeometryType () const
unsigned int	getLayer (DetId::Detector type, bool nose=false) const
unsigned int	getLayer (ForwardSubdetector type) const
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
GlobalPoint	getPositionLayer (int layer, bool nose=false) 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::pair< float, float >	getScintDEtaDPhi (const DetId &) const
int	getScintMaxIphi () const
int	getSiThickIndex (const DetId &) const
std::float_t	getSiThickness (const DetId &) const
const CaloSubdetectorGeometry *	getSubdetectorGeometry (const DetId &id) const
std::pair< int, int >	getWafer (const DetId &) const
bool	isHalfCell (const DetId &) const
bool	isOnlySilicon (const unsigned int layer) const
bool	isScintillator (const DetId &) const
bool	isSilicon (const DetId &) const
unsigned int	lastLayer (bool nose=false) const
unsigned int	lastLayerBH () const
unsigned int	lastLayerEE (bool nose=false) const
unsigned int	lastLayerFH () const
bool	maskCell (const DetId &id, int corners=3) const
unsigned int	maxNumberOfWafersPerLayer (bool nose=false) const
	RecHitTools ()
void	setGeometry (CaloGeometry const &)
int	zside (const DetId &id) const
	~RecHitTools ()

Private Attributes
unsigned int	bhFirstLayer_
unsigned int	bhLastLayer_
int	bhMaxIphi_
unsigned int	bhOffset_
unsigned int	eeOffset_
unsigned int	fhLastLayer_
unsigned int	fhOffset_
const CaloGeometry *	geom_
int	geometryType_
unsigned int	maxNumberOfWafersNose_
unsigned int	maxNumberOfWafersPerLayer_
unsigned int	noseLastLayer_

Detailed Description

Definition at line 23 of file RecHitTools.h.

Constructor & Destructor Documentation

hgcal::RecHitTools::RecHitTools ( )

inline

Definition at line 25 of file RecHitTools.h.

        : geom_(nullptr),
          eeOffset_(0),
          fhOffset_(0),
          bhFirstLayer_(0),
          bhOffset_(0),
          fhLastLayer_(0),
          noseLastLayer_(0),
          geometryType_(0) {}

hgcal::RecHitTools::~RecHitTools ( )

inline

Definition at line 34 of file RecHitTools.h.

{}

Member Function Documentation

std::pair< uint32_t, uint32_t > RecHitTools::firstAndLastLayer	(	DetId::Detector	det,
		int	subdet
	)		const

Definition at line 476 of file RecHitTools.cc.

References bhFirstLayer_, bhLastLayer_, eeOffset_, fhLastLayer_, fhOffset_, DetId::Forward, HFNose, DetId::HGCalEE, DetId::HGCalHSi, HGCEE, HGCHEF, and noseLastLayer_.

                                                                                              {
  if ((det == DetId::HGCalEE) || ((det == DetId::Forward) && (subdet == HGCEE))) {
    return std::make_pair(eeOffset_ + 1, fhOffset_);
  } else if ((det == DetId::HGCalHSi) || ((det == DetId::Forward) && (subdet == HGCHEF))) {
    return std::make_pair(fhOffset_ + 1, fhLastLayer_);
  } else if ((det == DetId::Forward) && (subdet == HFNose)) {
    return std::make_pair(1, noseLastLayer_);
  } else {
    return std::make_pair(bhFirstLayer_, bhLastLayer_);
  }
}

unsigned int hgcal::RecHitTools::firstLayerBH ( ) const

inline

Definition at line 77 of file RecHitTools.h.

References bhFirstLayer_.

Referenced by HGCalClusteringAlgoBase::getEventSetup(), HGCal3DClustering::makeClusters(), and HGCalDepthPreClusterer::makePreClusters().

{ return bhFirstLayer_; }

std::pair< int, int > RecHitTools::getCell

(

const DetId &

id

)

const

Definition at line 397 of file RecHitTools.cc.

References HGCalDetId::cell(), HFNoseDetId::cellU(), HGCSiliconDetId::cellU(), HFNoseDetId::cellV(), HGCSiliconDetId::cellV(), TauDecayModes::dec, DetId::Forward, HFNose, DetId::HGCalEE, DetId::HGCalHSi, and SiStripPI::max.

                                                            {
  int cellU = std::numeric_limits<int>::max();
  int cellV = 0;
  if ((id.det() == DetId::HGCalEE) || (id.det() == DetId::HGCalHSi)) {
    cellU = HGCSiliconDetId(id).cellU();
    cellV = HGCSiliconDetId(id).cellV();
  } else if (id.det() == DetId::Forward && id.subdetId() == static_cast<int>(HFNose)) {
    cellU = HFNoseDetId(id).cellU();
    cellV = HFNoseDetId(id).cellV();
  } else if (id.det() == DetId::Forward) {
    cellU = HGCalDetId(id).cell();
  } else {
    edm::LogError("getCell::InvalidSiliconDetid")
        << "det id: " << std::hex << id.rawId() << std::dec << ":" << id.det() << " is not HGCal silicon!";
  }
  return std::pair<int, int>(cellU, cellV);
}

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

Definition at line 441 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(), getPt(), and HGCalIsoCalculator::setRecHits().

                                                                                  {
  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 446 of file RecHitTools.cc.

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

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

const CaloGeometry* hgcal::RecHitTools::getGeometry ( ) const

inline

Definition at line 74 of file RecHitTools.h.

References geom_.

{ return geom_; };

int hgcal::RecHitTools::getGeometryType ( ) const

inline

Definition at line 85 of file RecHitTools.h.

References geometryType_.

{ return geometryType_; }

unsigned int RecHitTools::getLayer	(	DetId::Detector	type,
		bool	nose = false
	)		const

Definition at line 307 of file RecHitTools.cc.

References DetId::Forward, ForwardEmpty, geom_, CaloGeometry::getSubdetectorGeometry(), HFNose, DetId::HGCalEE, DetId::HGCalHSc, DetId::HGCalHSi, phase1PixelTopology::layer, and LayerTriplets::layers().

Referenced by getLayerWithOffset().

                                                                            {
  int layer;
  switch (type) {
    case (DetId::HGCalEE): {
      auto geomEE =
          static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(type, ForwardSubdetector::ForwardEmpty));
      layer = (geomEE->topology().dddConstants()).layers(true);
      break;
    }
    case (DetId::HGCalHSi): {
      auto geomFH =
          static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(type, ForwardSubdetector::ForwardEmpty));
      layer = (geomFH->topology().dddConstants()).layers(true);
      break;
    }
    case (DetId::HGCalHSc): {
      auto geomBH =
          static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(type, ForwardSubdetector::ForwardEmpty));
      layer = (geomBH->topology().dddConstants()).layers(true);
      break;
    }
    case (DetId::Forward): {
      if (nose) {
        auto geomHFN = static_cast<const HGCalGeometry*>(
            geom_->getSubdetectorGeometry(DetId::Forward, ForwardSubdetector::HFNose));
        layer = (geomHFN->topology().dddConstants()).layers(true);
      } else {
        auto geomEE = static_cast<const HGCalGeometry*>(
            geom_->getSubdetectorGeometry(DetId::HGCalEE, ForwardSubdetector::ForwardEmpty));
        layer = (geomEE->topology().dddConstants()).layers(true);
        auto geomFH = static_cast<const HGCalGeometry*>(
            geom_->getSubdetectorGeometry(DetId::HGCalHSi, ForwardSubdetector::ForwardEmpty));
        layer += (geomFH->topology().dddConstants()).layers(true);
      }
      break;
    }
    default:
      layer = 0;
  }
  return (unsigned int)(layer);
}

unsigned int RecHitTools::getLayer

(

ForwardSubdetector

type

)

const

Definition at line 257 of file RecHitTools.cc.

References DetId::Forward, ForwardEmpty, geom_, CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, HcalEndcap, HFNose, HGCEE, HGCHEB, HGCHEF, phase1PixelTopology::layer, and LayerTriplets::layers().

                                                                      {
  int layer(0);
  switch (type) {
    case (ForwardSubdetector::HGCEE): {
      auto geomEE =
          static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward, ForwardSubdetector::HGCEE));
      layer = (geomEE->topology().dddConstants()).layers(true);
      break;
    }
    case (ForwardSubdetector::HGCHEF): {
      auto geomFH =
          static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward, ForwardSubdetector::HGCHEF));
      layer = (geomFH->topology().dddConstants()).layers(true);
      break;
    }
    case (ForwardSubdetector::HGCHEB): {
      auto geomBH =
          static_cast<const HcalGeometry*>(geom_->getSubdetectorGeometry(DetId::Hcal, HcalSubdetector::HcalEndcap));
      layer = (geomBH->topology().dddConstants())->getMaxDepth(1);
      break;
    }
    case (ForwardSubdetector::ForwardEmpty): {
      auto geomEE =
          static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward, ForwardSubdetector::HGCEE));
      layer = (geomEE->topology().dddConstants()).layers(true);
      auto geomFH =
          static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward, ForwardSubdetector::HGCHEF));
      layer += (geomFH->topology().dddConstants()).layers(true);
      auto geomBH =
          static_cast<const HcalGeometry*>(geom_->getSubdetectorGeometry(DetId::Hcal, HcalSubdetector::HcalEndcap));
      if (geomBH)
        layer += (geomBH->topology().dddConstants())->getMaxDepth(1);
      auto geomBH2 =
          static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward, ForwardSubdetector::HGCHEB));
      if (geomBH2)
        layer += (geomBH2->topology().dddConstants()).layers(true);
      break;
    }
    case (ForwardSubdetector::HFNose): {
      auto geomHFN =
          static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward, ForwardSubdetector::HFNose));
      layer = (geomHFN->topology().dddConstants()).layers(true);
      break;
    }
    default:
      layer = 0;
  }
  return (unsigned int)(layer);
}

unsigned int RecHitTools::getLayer

(

const DetId &

id

)

const

Definition at line 349 of file RecHitTools.cc.

References HcalDetId::depth(), DetId::Forward, DetId::Hcal, HcalEndcap, HFNose, DetId::HGCalEE, DetId::HGCalHSc, DetId::HGCalHSi, HGCalDetId::layer(), HGCScintillatorDetId::layer(), HFNoseDetId::layer(), HGCSiliconDetId::layer(), phase1PixelTopology::layer, and SiStripPI::max.

                                                        {
  unsigned int layer = std::numeric_limits<unsigned int>::max();
  if (id.det() == DetId::HGCalEE || id.det() == DetId::HGCalHSi) {
    layer = HGCSiliconDetId(id).layer();
  } else if (id.det() == DetId::HGCalHSc) {
    layer = HGCScintillatorDetId(id).layer();
  } else if (id.det() == DetId::Forward && id.subdetId() == static_cast<int>(HFNose)) {
    layer = HFNoseDetId(id).layer();
  } else if (id.det() == DetId::Forward) {
    layer = HGCalDetId(id).layer();
  } else if (id.det() == DetId::Hcal && id.subdetId() == HcalEndcap) {
    layer = HcalDetId(id).depth();
  }
  return layer;
}

unsigned int RecHitTools::getLayerWithOffset

(

const DetId &

id

)

const

Definition at line 365 of file RecHitTools.cc.

References bhOffset_, fhOffset_, DetId::Forward, getLayer(), DetId::Hcal, HcalEndcap, DetId::HGCalHSc, DetId::HGCalHSi, HGCHEF, and phase1PixelTopology::layer.

Referenced by FWRecoGeometryESProducer::addCaloGeometry(), HGCalShowerSeparation::analyze(), HGCalHitCalibration::analyze(), TICLTrackstersEdgesValidation::dqmAnalyze(), TrackstersMergeProducer::energyRegressionAndID(), HGCalHitCalibration::fillWithRecHits(), ticl::ClusterFilterByAlgoAndSizeAndLayerRange::filter(), hgcal::ClusterTools::getLayer(), HGCalRecHitSimpleAlgo::makeRecHit(), HGCal3DClustering::organizeByLayer(), TICLLayerTileProducer::produce(), HGCalIsoCalculator::produceHGCalIso(), and hgcal::EGammaPCAHelper::storeRecHits().

                                                                  {
  unsigned int layer = getLayer(id);
  if (id.det() == DetId::Forward && id.subdetId() == HGCHEF) {
    layer += fhOffset_;
  } else if (id.det() == DetId::HGCalHSi || id.det() == DetId::HGCalHSc) {
    // DetId::HGCalHSc hits include the offset w.r.t. EE already
    layer += fhOffset_;
  } else if (id.det() == DetId::Hcal && id.subdetId() == HcalEndcap) {
    layer += bhOffset_;
  }
  // no need to add offset for HFnose
  return layer;
}

float RecHitTools::getPhi

(

const GlobalPoint &

position

)

const

Definition at line 452 of file RecHitTools.cc.

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

Referenced by HGCalIsoCalculator::setRecHits().

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

float RecHitTools::getPhi

(

const DetId &

id

)

const

Definition at line 457 of file RecHitTools.cc.

References getPosition(), 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 129 of file RecHitTools.cc.

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

Referenced by HGCalShowerSeparation::analyze(), getEta(), getPhi(), getPt(), hgcal::ClusterTools::getWidths(), HGCalImagingAlgo::Hexel::Hexel(), HGCalIsoCalculator::setRecHits(), and hgcal::EGammaPCAHelper::storeRecHits().

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

GlobalPoint RecHitTools::getPositionLayer	(	int	layer,
		bool	nose = false
	)		const

Definition at line 141 of file RecHitTools.cc.

References funct::abs(), fhOffset_, DetId::Forward, geom_, geometryType_, CaloGeometry::getSubdetectorGeometry(), HFNose, and HGCalDDDConstants::waferZ().

Referenced by PositionAtECalEntranceComputer::beginEvent(), hgcal::EGammaPCAHelper::findZFirstLayer(), TrackstersMergeProducer::produce(), and SimTrackstersProducer::produce().

                                                                    {
  unsigned int lay = std::abs(layer);
  double z(0);
  if (nose) {
    auto geomNose =
        static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward, ForwardSubdetector::HFNose));
    if (geomNose) {
      const HGCalDDDConstants* ddd = &(geomNose->topology().dddConstants());
      if (ddd)
        z = (layer > 0) ? ddd->waferZ(lay, true) : -ddd->waferZ(lay, true);
    }
  } else {
    const HGCalDDDConstants* ddd = get_ddd(geom_, geometryType_, fhOffset_, lay);
    if (geometryType_ == 1) {
      if (lay > fhOffset_)
        lay -= fhOffset_;
    }
    z = (layer > 0) ? ddd->waferZ(lay, true) : -ddd->waferZ(lay, true);
  }
  return GlobalPoint(0, 0, z);
}

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

Definition at line 463 of file RecHitTools.cc.

References PVValHelper::eta, getEta(), and DiDispStaMuonMonitor_cfi::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 469 of file RecHitTools.cc.

References PVValHelper::eta, getEta(), getPosition(), position, and DiDispStaMuonMonitor_cfi::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 235 of file RecHitTools.cc.

References HGCalDDDConstants::cellSizeHex(), TauDecayModes::dec, DetId::Forward, getSubdetectorGeometry(), HFNose, DetId::HGCalEE, DetId::HGCalHSi, SiStripPI::max, findQualityFiles::size, HFNoseDetId::type(), HGCSiliconDetId::type(), and HGCalDetId::waferType().

                                                             {
  auto geom = getSubdetectorGeometry(id);
  std::float_t size(std::numeric_limits<std::float_t>::max());
  if (id.det() == DetId::HGCalEE || id.det() == DetId::HGCalHSi) {
    const HGCSiliconDetId hid(id);
    auto ddd = get_ddd(geom, hid);
    size = ddd->cellSizeHex(hid.type());
  } else if (id.det() == DetId::Forward && id.subdetId() == static_cast<int>(HFNose)) {
    const HFNoseDetId hid(id);
    auto ddd = get_ddd(geom, hid);
    size = ddd->cellSizeHex(hid.type());
  } else if (id.det() == DetId::Forward) {
    const HGCalDetId hid(id);
    auto ddd = get_ddd(geom, hid);
    size = ddd->cellSizeHex(hid.waferType());
  } else {
    edm::LogError("getRadiusToSide::InvalidSiliconDetid")
        << "det id: " << std::hex << id.rawId() << std::dec << ":" << id.det() << " is not HGCal silicon!";
  }
  return size;
}

std::pair< float, float > RecHitTools::getScintDEtaDPhi

(

const DetId &

id

)

const

Definition at line 225 of file RecHitTools.cc.

References TauDecayModes::dec, CaloSubdetectorGeometry::getGeometry(), getSubdetectorGeometry(), isScintillator(), and LogDebug.

                                                                         {
  if (!isScintillator(id)) {
    LogDebug("getScintDEtaDPhi::InvalidScintDetid")
        << "det id: " << std::hex << id.rawId() << std::dec << ":" << id.det() << " is not HGCal scintillator!";
    return {0.f, 0.f};
  }
  auto cellGeom = getSubdetectorGeometry(id)->getGeometry(id);
  return {cellGeom->etaSpan(), cellGeom->phiSpan()};
}

int hgcal::RecHitTools::getScintMaxIphi ( ) const

inline

Definition at line 84 of file RecHitTools.h.

References bhMaxIphi_.

Referenced by HGCalClusteringAlgoBase::getEventSetup().

{ return bhMaxIphi_; }

int RecHitTools::getSiThickIndex

(

const DetId &

id

)

const

Definition at line 205 of file RecHitTools.cc.

References cms::cuda::assert(), DetId::Forward, getSiThickness(), HFNose, DetId::HGCalEE, DetId::HGCalHSi, TrackerMaterial_cfi::thickness, HFNoseDetId::type(), and HGCSiliconDetId::type().

Referenced by FWRecoGeometryESProducer::addCaloGeometry(), and hgcal::EGammaPCAHelper::storeRecHits().

                                                      {
  int thickIndex = -1;
  if (id.det() == DetId::HGCalEE || id.det() == DetId::HGCalHSi) {
    thickIndex = HGCSiliconDetId(id).type();
  } else if (id.det() == DetId::Forward && id.subdetId() == static_cast<int>(HFNose)) {
    thickIndex = HFNoseDetId(id).type();
  } else if (id.det() == DetId::Forward) {
    float thickness = getSiThickness(id);
    if (thickness > 99. && thickness < 121.)
      thickIndex = 0;
    else if (thickness > 199. && thickness < 201.)
      thickIndex = 1;
    else if (thickness > 299. && thickness < 301.)
      thickIndex = 2;
    else
      assert(thickIndex > 0 && "ERROR - silicon thickness has a nonsensical value");
  }
  return thickIndex;
}

std::float_t RecHitTools::getSiThickness

(

const DetId &

id

)

const

Definition at line 179 of file RecHitTools.cc.

References HGCalDDDConstants::cellThickness(), TauDecayModes::dec, DetId::Forward, getSubdetectorGeometry(), HFNose, DetId::HGCalEE, DetId::HGCalHSi, HGCalDetId::layer(), HFNoseDetId::layer(), HGCSiliconDetId::layer(), LogDebug, HGCalDetId::wafer(), HFNoseDetId::waferU(), HGCSiliconDetId::waferU(), HGCSiliconDetId::waferV(), and HFNoseDetId::waferV().

Referenced by HGCalHitCalibration::analyze(), HGCalHitCalibration::fillWithRecHits(), and getSiThickIndex().

                                                            {
  auto geom = getSubdetectorGeometry(id);
  std::float_t thick(0.37);
  if (id.det() == DetId::HGCalEE || id.det() == DetId::HGCalHSi) {
    const HGCSiliconDetId hid(id);
    auto ddd = get_ddd(geom, hid);
    thick = ddd->cellThickness(hid.layer(), hid.waferU(), hid.waferV());
  } else if (id.det() == DetId::Forward && id.subdetId() == static_cast<int>(HFNose)) {
    const HFNoseDetId hid(id);
    auto ddd = get_ddd(geom, hid);
    thick = ddd->cellThickness(hid.layer(), hid.waferU(), hid.waferV());
  } else if (id.det() == DetId::Forward) {
    const HGCalDetId hid(id);
    auto ddd = get_ddd(geom, hid);
    thick = ddd->cellThickness(hid.layer(), hid.wafer(), 0);
  } else {
    LogDebug("getSiThickness::InvalidSiliconDetid")
        << "det id: " << std::hex << id.rawId() << std::dec << ":" << id.det() << " is not HGCal silicon!";
    // It does not make any sense to return 0.37 as thickness for a detector
    // that is not Silicon(does it?). Mark it as 0. to be able to distinguish
    // the case.
    thick = 0.f;
  }
  return thick;
}

const CaloSubdetectorGeometry * RecHitTools::getSubdetectorGeometry

(

const DetId &

id

)

const

Definition at line 119 of file RecHitTools.cc.

References ForwardEmpty, relativeConstraints::geom, geom_, CaloGeometry::getSubdetectorGeometry(), DetId::HGCalEE, DetId::HGCalHSc, and DetId::HGCalHSi.

Referenced by getPosition(), getRadiusToSide(), getScintDEtaDPhi(), getSiThickness(), isHalfCell(), and maskCell().

                                                                                        {
  DetId::Detector det = id.det();
  int subdet = (det == DetId::HGCalEE || det == DetId::HGCalHSi || det == DetId::HGCalHSc)
                   ? ForwardSubdetector::ForwardEmpty
                   : id.subdetId();
  auto geom = geom_->getSubdetectorGeometry(det, subdet);
  check_geom(geom);
  return geom;
}

std::pair< int, int > RecHitTools::getWafer

(

const DetId &

id

)

const

Definition at line 379 of file RecHitTools.cc.

References TauDecayModes::dec, DetId::Forward, HFNose, DetId::HGCalEE, DetId::HGCalHSi, SiStripPI::max, HGCalDetId::wafer(), HGCalWaferIndex::waferU(), HFNoseDetId::waferU(), HGCSiliconDetId::waferU(), HGCalWaferIndex::waferV(), HGCSiliconDetId::waferV(), and HFNoseDetId::waferV().

                                                             {
  int waferU = std::numeric_limits<int>::max();
  int waferV = 0;
  if ((id.det() == DetId::HGCalEE) || (id.det() == DetId::HGCalHSi)) {
    waferU = HGCSiliconDetId(id).waferU();
    waferV = HGCSiliconDetId(id).waferV();
  } else if (id.det() == DetId::Forward && id.subdetId() == static_cast<int>(HFNose)) {
    waferU = HFNoseDetId(id).waferU();
    waferV = HFNoseDetId(id).waferV();
  } else if (id.det() == DetId::Forward) {
    waferU = HGCalDetId(id).wafer();
  } else {
    edm::LogError("getWafer::InvalidSiliconDetid")
        << "det id: " << std::hex << id.rawId() << std::dec << ":" << id.det() << " is not HGCal silicon!";
  }
  return std::pair<int, int>(waferU, waferV);
}

bool RecHitTools::isHalfCell

(

const DetId &

id

)

const

Definition at line 415 of file RecHitTools.cc.

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

                                                  {
  bool ishalf = false;
  if (id.det() == DetId::Forward) {
    HGCalDetId hid(id);
    auto geom = getSubdetectorGeometry(hid);
    auto ddd = get_ddd(geom, hid);
    const int waferType = ddd->waferTypeT(hid.waferType());
    return ddd->isHalfCell(waferType, hid.cell());
  }
  //new geometry is always false
  return ishalf;
}

bool RecHitTools::isOnlySilicon

(

const unsigned int

layer

)

const

Definition at line 435 of file RecHitTools.cc.

References bhLastLayer_, and bhOffset_.

                                                              {
  // HFnose TODO
  bool isonlysilicon = (layer % bhLastLayer_) < bhOffset_;
  return isonlysilicon;
}

bool RecHitTools::isScintillator

(

const DetId &

id

)

const

Definition at line 433 of file RecHitTools.cc.

References DetId::HGCalHSc.

Referenced by getScintDEtaDPhi().

{ return id.det() == DetId::HGCalHSc; }

bool RecHitTools::isSilicon

(

const DetId &

id

)

const

Definition at line 428 of file RecHitTools.cc.

References DetId::Forward, HFNose, DetId::HGCalEE, and DetId::HGCalHSi.

Referenced by FWRecoGeometryESProducer::addCaloGeometry(), ticl::ClusterFilterByAlgoAndSize::filter(), and ticl::ClusterFilterByAlgoAndSizeAndLayerRange::filter().

                                                 {
  return (id.det() == DetId::HGCalEE || id.det() == DetId::HGCalHSi ||
          (id.det() == DetId::Forward && id.subdetId() == static_cast<int>(HFNose)));
}

unsigned int hgcal::RecHitTools::lastLayer ( bool  nose = false ) const

inline

Definition at line 79 of file RecHitTools.h.

References bhLastLayer_, and noseLastLayer_.

Referenced by HGCalClusteringAlgoBase::getEventSetup(), TICLLayerTileProducer::produce(), and setGeometry().

{ return (nose ? noseLastLayer_ : bhLastLayer_); }

unsigned int hgcal::RecHitTools::lastLayerBH ( ) const

inline

Definition at line 78 of file RecHitTools.h.

References bhLastLayer_.

Referenced by HGCal3DClustering::getEventSetup(), and hgcal::EGammaPCAHelper::setRecHitTools().

{ return bhLastLayer_; }

unsigned int hgcal::RecHitTools::lastLayerEE ( bool  nose = false ) const

inline

Definition at line 75 of file RecHitTools.h.

References fhOffset_, and HFNoseDetId::HFNoseLayerEEmax.

Referenced by FWRecoGeometryESProducer::addCaloGeometry(), hgcal::EGammaPCAHelper::computePCA(), HGCalClusteringAlgoBase::getEventSetup(), hgcal::ClusterTools::getWidths(), HGCal3DClustering::makeClusters(), HGCalDepthPreClusterer::makePreClusters(), TrackstersMergeProducer::produce(), and SimTrackstersProducer::produce().

{ return (nose ? HFNoseDetId::HFNoseLayerEEmax : fhOffset_); }

unsigned int hgcal::RecHitTools::lastLayerFH ( ) const

inline

Definition at line 76 of file RecHitTools.h.

References fhLastLayer_.

Referenced by HGCalClusteringAlgoBase::getEventSetup().

{ return fhLastLayer_; }

bool RecHitTools::maskCell	(	const DetId &	id,
		int	corners = 3
	)		const

Definition at line 488 of file RecHitTools.cc.

References HGCalTopology::dddConstants(), getSubdetectorGeometry(), DetId::Hcal, HGCalDDDConstants::maskCell(), and HGCalGeometry::topology().

                                                             {
  if (id.det() == DetId::Hcal) {
    return false;
  } else {
    auto hg = static_cast<const HGCalGeometry*>(getSubdetectorGeometry(id));
    return hg->topology().dddConstants().maskCell(id, corners);
  }
}

unsigned int hgcal::RecHitTools::maxNumberOfWafersPerLayer ( bool  nose = false ) const

inline

Definition at line 81 of file RecHitTools.h.

References maxNumberOfWafersNose_, and maxNumberOfWafersPerLayer_.

                                                                    {
      return (nose ? maxNumberOfWafersNose_ : maxNumberOfWafersPerLayer_);
    }

void RecHitTools::setGeometry

(

CaloGeometry const &

geom

)

Definition at line 68 of file RecHitTools.cc.

References bhFirstLayer_, bhLastLayer_, bhMaxIphi_, bhOffset_, eeOffset_, fhLastLayer_, fhOffset_, DetId::Forward, ForwardEmpty, relativeConstraints::geom, geom_, geometryType_, CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, HcalEndcap, HFNose, DetId::HGCalEE, DetId::HGCalHSc, DetId::HGCalHSi, HGCEE, HGCHEF, lastLayer(), LayerTriplets::layers(), SiStripPI::max, maxNumberOfWafersNose_, maxNumberOfWafersPerLayer_, and noseLastLayer_.

Referenced by FWRecoGeometryESProducer::addCaloGeometry(), HeterogeneousHGCalRecHitsValidator::analyze(), HGCalShowerSeparation::analyze(), HGCalHitCalibration::analyze(), PositionAtECalEntranceComputer::beginEvent(), TICLLayerTileProducer::beginRun(), FilteredLayerClustersProducer::beginRun(), TICLTrackstersEdgesValidation::dqmBeginRun(), HGCalEgammaIDHelper::eventInit(), hgcal::ClusterTools::getEventSetup(), HGCal3DClustering::getEventSetup(), HGCalDepthPreClusterer::getEventSetup(), HGCalClusteringAlgoBase::getEventSetup(), PFHGCalRecHitCreator< DET, Layer, det, subdet >::importRecHits(), TrackstersMergeProducer::produce(), SimTrackstersProducer::produce(), and HGCalRecHitAbsAlgo::set().

                                                      {
  geom_ = &geom;
  unsigned int wmaxEE(0), wmaxFH(0);
  auto geomEE = static_cast<const HGCalGeometry*>(
      geom_->getSubdetectorGeometry(DetId::HGCalEE, ForwardSubdetector::ForwardEmpty));
  //check if it's the new geometry
  if (geomEE) {
    geometryType_ = 1;
    eeOffset_ = (geomEE->topology().dddConstants()).getLayerOffset();
    wmaxEE = (geomEE->topology().dddConstants()).waferCount(0);
    auto geomFH = static_cast<const HGCalGeometry*>(
        geom_->getSubdetectorGeometry(DetId::HGCalHSi, ForwardSubdetector::ForwardEmpty));
    fhOffset_ = (geomFH->topology().dddConstants()).getLayerOffset();
    wmaxFH = (geomFH->topology().dddConstants()).waferCount(0);
    fhLastLayer_ = fhOffset_ + (geomFH->topology().dddConstants()).lastLayer(true);
    auto geomBH = static_cast<const HGCalGeometry*>(
        geom_->getSubdetectorGeometry(DetId::HGCalHSc, ForwardSubdetector::ForwardEmpty));
    bhOffset_ = (geomBH->topology().dddConstants()).getLayerOffset();
    bhFirstLayer_ = bhOffset_ + (geomBH->topology().dddConstants()).firstLayer();
    bhLastLayer_ = bhOffset_ + (geomBH->topology().dddConstants()).lastLayer(true);
    bhMaxIphi_ = geomBH->topology().dddConstants().maxCells(true);
  } else {
    geometryType_ = 0;
    geomEE =
        static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward, ForwardSubdetector::HGCEE));
    eeOffset_ = (geomEE->topology().dddConstants()).getLayerOffset();
    wmaxEE = 1 + (geomEE->topology().dddConstants()).waferMax();
    auto geomFH =
        static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward, ForwardSubdetector::HGCHEF));
    fhOffset_ = (geomFH->topology().dddConstants()).getLayerOffset();
    fhLastLayer_ = fhOffset_ + (geomFH->topology().dddConstants()).layers(true);
    bhOffset_ = fhLastLayer_;
    bhFirstLayer_ = bhOffset_ + 1;
    wmaxFH = 1 + (geomFH->topology().dddConstants()).waferMax();
    auto geomBH =
        static_cast<const HcalGeometry*>(geom_->getSubdetectorGeometry(DetId::Hcal, HcalSubdetector::HcalEndcap));
    bhLastLayer_ = bhOffset_ + (geomBH->topology().dddConstants())->getMaxDepth(1);
  }
  maxNumberOfWafersPerLayer_ = std::max(wmaxEE, wmaxFH);
  // For nose geometry
  auto geomNose =
      static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward, ForwardSubdetector::HFNose));
  if (geomNose) {
    maxNumberOfWafersNose_ = (geomNose->topology().dddConstants()).waferCount(0);
    noseLastLayer_ = (geomNose->topology().dddConstants()).layers(true);
  } else {
    maxNumberOfWafersNose_ = 0;
    noseLastLayer_ = 0;
  }
}

int RecHitTools::zside

(

const DetId &

id

)

const

Definition at line 163 of file RecHitTools.cc.

References DetId::Forward, DetId::Hcal, HcalEndcap, HFNose, DetId::HGCalEE, DetId::HGCalHSc, DetId::HGCalHSi, HGCScintillatorDetId::zside(), HGCalDetId::zside(), HFNoseDetId::zside(), HGCSiliconDetId::zside(), and HcalDetId::zside().

Referenced by FWRecoGeometryESProducer::addCaloGeometry(), and TICLLayerTileProducer::produce().

                                            {
  int zside = 0;
  if (id.det() == DetId::HGCalEE || id.det() == DetId::HGCalHSi) {
    zside = HGCSiliconDetId(id).zside();
  } else if (id.det() == DetId::HGCalHSc) {
    zside = HGCScintillatorDetId(id).zside();
  } else if (id.det() == DetId::Forward && id.subdetId() == static_cast<int>(HFNose)) {
    zside = HFNoseDetId(id).zside();
  } else if (id.det() == DetId::Forward) {
    zside = HGCalDetId(id).zside();
  } else if (id.det() == DetId::Hcal && id.subdetId() == HcalEndcap) {
    zside = HcalDetId(id).zside();
  }
  return zside;
}

Member Data Documentation

unsigned int hgcal::RecHitTools::bhFirstLayer_

private

Definition at line 90 of file RecHitTools.h.

Referenced by firstAndLastLayer(), firstLayerBH(), and setGeometry().

unsigned int hgcal::RecHitTools::bhLastLayer_

private

Definition at line 90 of file RecHitTools.h.

Referenced by firstAndLastLayer(), isOnlySilicon(), lastLayer(), lastLayerBH(), and setGeometry().

int hgcal::RecHitTools::bhMaxIphi_

private

Definition at line 93 of file RecHitTools.h.

Referenced by getScintMaxIphi(), and setGeometry().

unsigned int hgcal::RecHitTools::bhOffset_

private

Definition at line 90 of file RecHitTools.h.

Referenced by getLayerWithOffset(), isOnlySilicon(), and setGeometry().

unsigned int hgcal::RecHitTools::eeOffset_

private

Definition at line 90 of file RecHitTools.h.

Referenced by firstAndLastLayer(), and setGeometry().

unsigned int hgcal::RecHitTools::fhLastLayer_

private

Definition at line 90 of file RecHitTools.h.

Referenced by firstAndLastLayer(), lastLayerFH(), and setGeometry().

unsigned int hgcal::RecHitTools::fhOffset_

private

Definition at line 90 of file RecHitTools.h.

Referenced by firstAndLastLayer(), getLayerWithOffset(), getPositionLayer(), lastLayerEE(), and setGeometry().

const CaloGeometry* hgcal::RecHitTools::geom_

private

Definition at line 89 of file RecHitTools.h.

Referenced by getGeometry(), getLayer(), getPositionLayer(), getSubdetectorGeometry(), and setGeometry().

int hgcal::RecHitTools::geometryType_

private

Definition at line 92 of file RecHitTools.h.

Referenced by getGeometryType(), getPositionLayer(), and setGeometry().

unsigned int hgcal::RecHitTools::maxNumberOfWafersNose_

private

Definition at line 91 of file RecHitTools.h.

Referenced by maxNumberOfWafersPerLayer(), and setGeometry().

unsigned int hgcal::RecHitTools::maxNumberOfWafersPerLayer_

private

Definition at line 91 of file RecHitTools.h.

Referenced by maxNumberOfWafersPerLayer(), and setGeometry().

unsigned int hgcal::RecHitTools::noseLastLayer_

private

Definition at line 90 of file RecHitTools.h.

Referenced by firstAndLastLayer(), lastLayer(), and setGeometry().