HGCalTriggerTools Class Reference

#include <HGCalTriggerTools.h>

Public Member Functions
void	eventSetup (const edm::EventSetup &)
float	getEta (const GlobalPoint &position, const float &vertex_z=0.) const
float	getLayerZ (const unsigned &layerWithOffset) const
float	getLayerZ (const int &subdet, const unsigned &layer) const
float	getPhi (const GlobalPoint &position) const
float	getPt (const GlobalPoint &position, const float &hitEnergy, const float &vertex_z=0.) const
float	getTCEta (const DetId &id, const float &vertex_z=0.) const
float	getTCPhi (const DetId &id) const
GlobalPoint	getTCPosition (const DetId &id) const
float	getTCPt (const DetId &id, const float &hitEnergy, const float &vertex_z=0.) const
const HGCalTriggerGeometryBase *	getTriggerGeometry () const
	HGCalTriggerTools ()
unsigned	lastLayerBH () const
unsigned	lastLayerEE () const
unsigned	lastLayerFH () const
unsigned	layer (const DetId &) const
unsigned	layers (ForwardSubdetector type) const
unsigned	layerWithOffset (const DetId &) const
	~HGCalTriggerTools ()

Private Attributes
unsigned	bhLayers_
unsigned	eeLayers_
unsigned	fhLayers_
const HGCalTriggerGeometryBase *	geom_
unsigned	totalLayers_

Detailed Description

Tools for handling HGCal trigger det-ID: in the current version of trhe HGCAL simulation only HGCalDetId for the TriggerCells (TC) are used and not HcalDetId as in the offline! As a consequence the class assumes that only DetIds of the first kind are used in the getTC* methods NOTE: this uses the trigger geometry hence would give wrong results when used for offline reco!!!!

Author

G. Cerminara (CERN), heavily "inspired" by HGCalRechHitTools ;)

Definition at line 30 of file HGCalTriggerTools.h.

Constructor & Destructor Documentation

HGCalTriggerTools::HGCalTriggerTools ( )

inline

Definition at line 32 of file HGCalTriggerTools.h.

                        : geom_(nullptr),
      eeLayers_(0), fhLayers_(0), bhLayers_(0), totalLayers_(0){}

HGCalTriggerTools::~HGCalTriggerTools ( )

inline

Definition at line 34 of file HGCalTriggerTools.h.

References LayerTriplets::layers().

{}

Member Function Documentation

void HGCalTriggerTools::eventSetup

(

const edm::EventSetup &

es

)

Definition at line 59 of file HGCalTriggerTools.cc.

References bhLayers_, HGCalTriggerGeometryBase::bhTopology(), HGCalTopology::dddConstants(), HcalTopology::dddConstants(), eeLayers_, HGCalTriggerGeometryBase::eeTopology(), fhLayers_, HGCalTriggerGeometryBase::fhTopology(), geom_, edm::EventSetup::get(), HcalDDDRecConstants::getMaxDepth(), HGCalTriggerGeometryBase::hscTopology(), HGCalTriggerGeometryBase::isV9Geometry(), HGCalDDDConstants::layers(), edm::ESHandle< T >::product(), and totalLayers_.

Referenced by HGCalMulticlusteringImpl::eventSetup(), HGCalTriggerCellCalibration::eventSetup(), HGCalShowerShape::eventSetup(), HGCalTowerMap2DImpl::eventSetup(), HGCalClusteringImpl::eventSetup(), HGCalTriggerNtupleHGCClusters::fill(), HGCalTriggerNtupleHGCTriggerCells::fill(), and HGCalTriggerNtupleHGCDigis::fill().

{
  edm::ESHandle<HGCalTriggerGeometryBase> triggerGeometry_;
  es.get<CaloGeometryRecord>().get(triggerGeometry_);
  geom_ = triggerGeometry_.product();

  eeLayers_ = geom_->eeTopology().dddConstants().layers(true);
  fhLayers_ = geom_->fhTopology().dddConstants().layers(true);
  if(geom_->isV9Geometry())
  {
    bhLayers_ = geom_->hscTopology().dddConstants().layers(true);
    totalLayers_ =  eeLayers_ + fhLayers_;
  }
  else
  {
    bhLayers_ = geom_->bhTopology().dddConstants()->getMaxDepth(1);
    totalLayers_ =  eeLayers_ + fhLayers_ + bhLayers_;
  }
}

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

Definition at line 145 of file HGCalTriggerTools.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 getPt(), and getTCEta().

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

float HGCalTriggerTools::getLayerZ

(

const unsigned &

layerWithOffset

)

const

Definition at line 176 of file HGCalTriggerTools.cc.

References HcalEndcap, HGCEE, HGCHEF, lastLayerEE(), lastLayerFH(), and PFRecoTauDiscriminationByIsolation_cfi::offset.

                                                                        {
  int subdet = ForwardSubdetector::HGCEE;
  unsigned offset = 0;
  if(layerWithOffset > lastLayerEE() && layerWithOffset <= lastLayerFH()) {
    subdet = ForwardSubdetector::HGCHEF;
    offset = lastLayerEE();
  } else if(layerWithOffset > lastLayerFH()) {
    subdet = HcalSubdetector::HcalEndcap;
    offset = lastLayerFH();
  }
  return getLayerZ(subdet, layerWithOffset - offset);
}

float HGCalTriggerTools::getLayerZ	(	const int &	subdet,
		const unsigned &	layer
	)		const

Definition at line 189 of file HGCalTriggerTools.cc.

References HGCalTriggerGeometryBase::bhTopology(), HGCalTopology::dddConstants(), HcalTopology::dddConstants(), HGCalTriggerGeometryBase::eeTopology(), HGCalTriggerGeometryBase::fhTopology(), geom_, HcalDDDRecConstants::getEtaRange(), HcalDDDRecConstants::getRZ(), HcalEndcap, HGCEE, HGCHEB, HGCHEF, HGCalTriggerGeometryBase::hscTopology(), HGCalTriggerGeometryBase::isV9Geometry(), layer(), and HGCalDDDConstants::waferZ().

                                                                                 {
  float layerGlobalZ = 0.;
  if(subdet == ForwardSubdetector::HGCEE) {
    layerGlobalZ = geom_->eeTopology().dddConstants().waferZ(layer, true);
  } else if(subdet == ForwardSubdetector::HGCHEF) {
    layerGlobalZ = geom_->fhTopology().dddConstants().waferZ(layer, true);
  } else if(subdet == HcalSubdetector::HcalEndcap || subdet == ForwardSubdetector::HGCHEB) {
    if(geom_->isV9Geometry())
    {
      layerGlobalZ = geom_->hscTopology().dddConstants().waferZ(layer, true);
    }
    else
    {
      layerGlobalZ = geom_->bhTopology().dddConstants()->getRZ(
          HcalSubdetector::HcalEndcap,
          geom_->bhTopology().dddConstants()->getEtaRange(1).second,
          layer);
    }
  }
  return layerGlobalZ;
}

float HGCalTriggerTools::getPhi

(

const GlobalPoint &

position

)

const

Definition at line 155 of file HGCalTriggerTools.cc.

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

Referenced by getTCPhi().

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

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

Definition at line 165 of file HGCalTriggerTools.cc.

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

Referenced by getTCPt().

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

float HGCalTriggerTools::getTCEta	(	const DetId &	id,
		const float &	vertex_z = 0.
	)		const

Definition at line 150 of file HGCalTriggerTools.cc.

References getEta(), getTCPosition(), and position.

                                                                              {
  GlobalPoint position = getTCPosition(id);
  return getEta(position, vertex_z);
}

float HGCalTriggerTools::getTCPhi

(

const DetId &

id

)

const

Definition at line 160 of file HGCalTriggerTools.cc.

References getPhi(), getTCPosition(), and position.

                                                       {
  GlobalPoint position = getTCPosition(id);
  return getPhi(position);
}

GlobalPoint HGCalTriggerTools::getTCPosition

(

const DetId &

id

)

const

Definition at line 79 of file HGCalTriggerTools.cc.

References Exception, geom_, HGCalTriggerGeometryBase::getTriggerCellPosition(), DetId::Hcal, layers(), and position.

Referenced by getTCEta(), getTCPhi(), and getTCPt().

                                                                  {
  if(id.det() == DetId::Hcal) {
    throw cms::Exception("hgcal::HGCalTriggerTools")
      << "method getTCPosition called for DetId not belonging to a TC";
    // FIXME: this would actually need a better test...but at the moment I can not think to anything better
    // to distinguish a TC detId
  }

  GlobalPoint position = geom_->getTriggerCellPosition(id);
  return position;
}

float HGCalTriggerTools::getTCPt	(	const DetId &	id,
		const float &	hitEnergy,
		const float &	vertex_z = 0.
	)		const

Definition at line 171 of file HGCalTriggerTools.cc.

References getPt(), getTCPosition(), and position.

                                                                                                     {
  GlobalPoint position = getTCPosition(id);
  return getPt(position, hitEnergy, vertex_z);
}

const HGCalTriggerGeometryBase* HGCalTriggerTools::getTriggerGeometry ( ) const

inline

Definition at line 56 of file HGCalTriggerTools.h.

{return geom_;};

unsigned HGCalTriggerTools::lastLayerBH ( ) const

inline

Definition at line 44 of file HGCalTriggerTools.h.

References position.

Referenced by HGCalTowerMap2DImpl::newTowerMaps().

{return totalLayers_;}

unsigned HGCalTriggerTools::lastLayerEE ( ) const

inline

Definition at line 42 of file HGCalTriggerTools.h.

Referenced by HGCalTowerMap2DImpl::buildTowerMap2D(), getLayerZ(), and HGCalTowerMap2DImpl::newTowerMaps().

{return eeLayers_;}

unsigned HGCalTriggerTools::lastLayerFH ( ) const

inline

Definition at line 43 of file HGCalTriggerTools.h.

Referenced by getLayerZ().

{return eeLayers_+fhLayers_;}

unsigned HGCalTriggerTools::layer

(

const DetId &

id

)

const

Definition at line 119 of file HGCalTriggerTools.cc.

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

Referenced by getLayerZ(), geometryXMLparser.DTAlignable::index(), geometryXMLparser.CSCAlignable::index(), layers(), and layerWithOffset().

                             {
  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 HGCalTriggerTools::layers

(

ForwardSubdetector

type

)

const

Definition at line 94 of file HGCalTriggerTools.cc.

References bhLayers_, eeLayers_, fhLayers_, ForwardEmpty, HGCEE, HGCHEB, HGCHEF, layer(), and totalLayers_.

Referenced by HGCalClusteringImpl::clusterizeNN(), HGCalShowerShape::coreShowerLength(), and getTCPosition().

{
  unsigned layers = 0;
  switch(type)
  {
    case ForwardSubdetector::HGCEE:
      layers = eeLayers_;
      break;
    case ForwardSubdetector::HGCHEF:
      layers = fhLayers_;
      break;
    case ForwardSubdetector::HGCHEB:
      layers = bhLayers_;
      break;
    case ForwardSubdetector::ForwardEmpty:
      layers = totalLayers_;
      break;
    default:
      break;
  };
  return layers;
}

unsigned HGCalTriggerTools::layerWithOffset

(

const DetId &

id

)

const

Definition at line 133 of file HGCalTriggerTools.cc.

References eeLayers_, fhLayers_, DetId::Forward, geom_, DetId::Hcal, HcalEndcap, HGCHEB, HGCHEF, HGCalTriggerGeometryBase::isV9Geometry(), checklumidiff::l, and layer().

Referenced by HGCalTowerMap2DImpl::buildTowerMap2D(), HGCalTriggerCellCalibration::calibrateMipTinGeV(), HGCalClusteringImpl::calibratePt(), HGCalShowerShape::eMax(), HGCalTriggerNtupleHGCClusters::fill(), HGCalTriggerNtupleHGCTriggerCells::fill(), HGCalTriggerNtupleHGCDigis::fill(), HGCalShowerShape::firstLayer(), HGCalShowerShape::lastLayer(), layer(), HGCalShowerShape::maxLayer(), HGCalShowerShape::sigmaEtaEtaMax(), HGCalShowerShape::sigmaPhiPhiMax(), HGCalShowerShape::sigmaRRMax(), HGCalShowerShape::sigmaRRMean(), and HGCalClusteringImpl::triggerCellReshuffling().

                                       {
  unsigned int l = layer(id);
  if( id.det() == DetId::Forward && id.subdetId() == HGCHEF ) {
    l += eeLayers_;
  } else if( (id.det() == DetId::Hcal && id.subdetId() == HcalEndcap) ||
             (id.det() == DetId::Forward && id.subdetId() == HGCHEB) ) {
    if(geom_->isV9Geometry()) l += eeLayers_;
    else l += eeLayers_ + fhLayers_;
  }
  return l;
}

Member Data Documentation

unsigned HGCalTriggerTools::bhLayers_

private

Definition at line 67 of file HGCalTriggerTools.h.

Referenced by eventSetup(), and layers().

unsigned HGCalTriggerTools::eeLayers_

private

Definition at line 65 of file HGCalTriggerTools.h.

Referenced by eventSetup(), layers(), and layerWithOffset().

unsigned HGCalTriggerTools::fhLayers_

private

Definition at line 66 of file HGCalTriggerTools.h.

Referenced by eventSetup(), layers(), and layerWithOffset().

const HGCalTriggerGeometryBase* HGCalTriggerTools::geom_

private

Definition at line 64 of file HGCalTriggerTools.h.

Referenced by eventSetup(), getLayerZ(), getTCPosition(), and layerWithOffset().

unsigned HGCalTriggerTools::totalLayers_

private

Definition at line 68 of file HGCalTriggerTools.h.

Referenced by eventSetup(), and layers().