HGCalTriggerTools Class Reference

#include <HGCalTriggerTools.h>

Public Types
enum	SubDetectorType { hgcal_silicon_CEE, hgcal_silicon_CEH, hgcal_scintillator }

Public Member Functions
template<typename T >
std::vector< T >	bxVectorToVector (const BXVector< T > &inputBXVector)
void	eventSetup (const edm::EventSetup &)
float	getEta (const GlobalPoint &position, const float &vertex_z=0.) const
float	getLayerZ (const int &subdet, const unsigned &layer) const
float	getLayerZ (const unsigned &layerWithOffset) const
float	getPhi (const GlobalPoint &position) const
float	getPt (const GlobalPoint &position, const float &hitEnergy, const float &vertex_z=0.) const
SubDetectorType	getSubDetectorType (const DetId &id) 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 ()
bool	isEm (const DetId &) const
bool	isHad (const DetId &id) const
bool	isNose (const DetId &) const
bool	isScintillator (const DetId &id) const
bool	isSilicon (const DetId &) const
unsigned	lastLayerBH () const
unsigned	lastLayerEE () const
unsigned	lastLayerFH () const
unsigned	lastLayerNose () const
unsigned	layer (const DetId &) const
unsigned	layers (DetId::Detector type) const
unsigned	layers (ForwardSubdetector type) const
unsigned	layerWithOffset (const DetId &) const
DetId	simToReco (const DetId &, const HcalTopology &) const
DetId	simToReco (const DetId &, const HGCalTopology &) const
int	thicknessIndex (const DetId &, bool tc=false) const
unsigned	triggerLayer (const unsigned id) const
int	zside (const DetId &) const
	~HGCalTriggerTools ()

Static Public Attributes
static constexpr unsigned	kScintillatorPseudoThicknessIndex_ = 3

Private Member Functions
int	sensorCellThicknessV8 (const DetId &id) const

Private Attributes
unsigned	bhLayers_
unsigned	eeLayers_
unsigned	fhLayers_
const HGCalTriggerGeometryBase *	geom_
unsigned	noseLayers_
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 32 of file HGCalTriggerTools.h.

Member Enumeration Documentation

SubDetectorType

enum HGCalTriggerTools::SubDetectorType

Enumerator
hgcal_silicon_CEE
hgcal_silicon_CEH
hgcal_scintillator

Definition at line 88 of file HGCalTriggerTools.h.

                       {
    hgcal_silicon_CEE,
    hgcal_silicon_CEH,
    hgcal_scintillator,
  };

Constructor & Destructor Documentation

HGCalTriggerTools()

HGCalTriggerTools::HGCalTriggerTools ( )

inline

Definition at line 34 of file HGCalTriggerTools.h.

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

~HGCalTriggerTools()

HGCalTriggerTools::~HGCalTriggerTools ( )

inline

Definition at line 35 of file HGCalTriggerTools.h.

{}

Member Function Documentation

bxVectorToVector()

template<typename T >

std::vector<T> HGCalTriggerTools::bxVectorToVector ( const BXVector< T > &  inputBXVector )

inline

Definition at line 75 of file HGCalTriggerTools.h.

                                                                  {
    std::vector<T> outputVector;
    // loop over collection for a given bx and put the objects into a std::vector
    outputVector.insert(outputVector.end(), inputBXVector.begin(0), inputBXVector.end(0));
    return outputVector;
  }

References BXVector< T >::begin(), and BXVector< T >::end().

eventSetup()

void HGCalTriggerTools::eventSetup

(

const edm::EventSetup &

es

)

Definition at line 35 of file HGCalTriggerTools.cc.

                                                          {
  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_->isWithNoseGeometry())
    noseLayers_ = geom_->noseTopology().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_;
  }
}

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

Referenced by HGCalTriggerValidator::analyze(), HGCalTriggerClusterInterpretationEM::eventSetup(), HGCalTriggerCellCalibration::eventSetup(), HGCalVFESummationImpl::eventSetup(), HGCalConcentratorCoarsenerImpl::eventSetup(), HGCalClusteringDummyImpl::eventSetup(), HGCalConcentratorThresholdImpl::eventSetup(), HGCalCoarseTriggerCellMapping::eventSetup(), HGCalConcentratorTrigSumImpl::eventSetup(), HGCalMulticlusteringImpl::eventSetup(), HGCalHistoClusteringImpl::eventSetup(), HGCalConcentratorBestChoiceImpl::eventSetup(), HGCalConcentratorSuperTriggerCellImpl::eventSetup(), HGCalShowerShape::eventSetup(), HGCalClusteringImpl::eventSetup(), HGCalTriggerTowerGeometryHelper::eventSetup(), HGCalTowerMap2DImpl::eventSetup(), HGCalHistoSeedingImpl::eventSetup(), HGCalTriggerNtupleHGCClusters::fill(), HGCalTriggerNtupleHGCTriggerSums::fill(), HGCalTriggerNtupleHGCDigis::fill(), HGCalTriggerNtupleHGCTriggerCells::fill(), HGCalTriggerNtupleGen::fill(), CaloTruthCellsProducer::produce(), HGCalVFEProcessorSums::run(), and HGCalConcentratorProcessorSelection::run().

getEta()

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

Definition at line 281 of file HGCalTriggerTools.cc.

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

References PV3DBase< T, PVType, FrameType >::eta(), and position.

Referenced by getPt(), and getTCEta().

getLayerZ() [1/2]

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

Definition at line 326 of file HGCalTriggerTools.cc.

                                                                                 {
  float layerGlobalZ = 0.;
  if ((subdet == ForwardSubdetector::HGCEE) || (subdet == DetId::HGCalEE)) {
    layerGlobalZ = geom_->eeTopology().dddConstants().waferZ(layer, true);
  } else if ((subdet == ForwardSubdetector::HGCHEF) || (subdet == DetId::HGCalHSi)) {
    layerGlobalZ = geom_->fhTopology().dddConstants().waferZ(layer, true);
  } else if (subdet == ForwardSubdetector::HFNose) {
    layerGlobalZ = geom_->noseTopology().dddConstants().waferZ(layer, true);
  } else if ((subdet == HcalSubdetector::HcalEndcap) || (subdet == ForwardSubdetector::HGCHEB) ||
             (subdet == DetId::HGCalHSc)) {
    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;
}

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

getLayerZ() [2/2]

float HGCalTriggerTools::getLayerZ

(

const unsigned &

layerWithOffset

)

const

Definition at line 312 of file HGCalTriggerTools.cc.

                                                                        {
  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();
  }
  // note for HFnose offset is always zero since we have less layers than HGCEE
  return getLayerZ(subdet, layerWithOffset - offset);
}

References HcalEndcap, HGCEE, HGCHEF, lastLayerEE(), lastLayerFH(), layerWithOffset(), and hltrates_dqm_sourceclient-live_cfg::offset.

Referenced by HGCalTriggerNtupleGen::fill().

getPhi()

float HGCalTriggerTools::getPhi

(

const GlobalPoint &

position

)

const

Definition at line 291 of file HGCalTriggerTools.cc.

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

References phi, and position.

Referenced by getTCPhi().

getPt()

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

Definition at line 301 of file HGCalTriggerTools.cc.

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

References PVValHelper::eta, getEta(), position, and DiDispStaMuonMonitor_cfi::pt.

Referenced by getTCPt().

getSubDetectorType()

HGCalTriggerTools::SubDetectorType HGCalTriggerTools::getSubDetectorType

(

const DetId &

id

)

const

Definition at line 198 of file HGCalTriggerTools.cc.

                                                                                            {
  SubDetectorType subdet;
  if (!isScintillator(id)) {
    if (isEm(id))
      subdet = HGCalTriggerTools::hgcal_silicon_CEE;
    else
      subdet = HGCalTriggerTools::hgcal_silicon_CEH;
  } else
    subdet = HGCalTriggerTools::hgcal_scintillator;
  return subdet;
}

References hgcal_scintillator, hgcal_silicon_CEE, hgcal_silicon_CEH, isEm(), and isScintillator().

Referenced by HGCalConcentratorSuperTriggerCellImpl::assignSuperTriggerCellEnergyAndPosition(), HGCalConcentratorSuperTriggerCellImpl::createAllTriggerCells(), and HGCalConcentratorProcessorSelection::run().

getTCEta()

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

Definition at line 286 of file HGCalTriggerTools.cc.

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

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

getTCPhi()

float HGCalTriggerTools::getTCPhi

(

const DetId &

id

)

const

Definition at line 296 of file HGCalTriggerTools.cc.

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

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

getTCPosition()

GlobalPoint HGCalTriggerTools::getTCPosition

(

const DetId &

id

)

const

Definition at line 54 of file HGCalTriggerTools.cc.

                                                                  {
  if (id.det() == DetId::Hcal || id.det() == DetId::HGCalEE) {
    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;
}

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

Referenced by HGCalConcentratorSuperTriggerCellImpl::assignSuperTriggerCellEnergyAndPosition(), HGCalCoarseTriggerCellMapping::getCoarseTriggerCellPosition(), getTCEta(), getTCPhi(), and getTCPt().

getTCPt()

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

Definition at line 307 of file HGCalTriggerTools.cc.

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

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

getTriggerGeometry()

const HGCalTriggerGeometryBase* HGCalTriggerTools::getTriggerGeometry ( ) const

inline

Definition at line 69 of file HGCalTriggerTools.h.

{ return geom_; };

References geom_.

Referenced by HGCalConcentratorSuperTriggerCellImpl::createAllTriggerCells(), HGCalConcentratorTrigSumImpl::doSum(), HGCalCoarseTriggerCellMapping::getConstituentTriggerCells(), HGCalCoarseTriggerCellMapping::getRepresentativeDetId(), and HGCalVFESummationImpl::triggerCellSums().

isEm()

bool HGCalTriggerTools::isEm

(

const DetId &

id

)

const

Definition at line 151 of file HGCalTriggerTools.cc.

                                                  {
  bool em = false;
 
  if (id.det() == DetId::Forward && id.subdetId() != ForwardSubdetector::HFNose) {
    em = (id.subdetId() == HGCEE);
  } else if (id.det() == DetId::Forward && id.subdetId() == ForwardSubdetector::HFNose) {
    em = HFNoseDetId(id).isEE();
  } else if (id.det() == DetId::HGCalEE) {
    em = true;
  } else if (id.det() == DetId::HGCalTrigger &&
             HGCalTriggerDetId(id).subdet() == HGCalTriggerSubdetector::HGCalEETrigger) {
    em = true;
  } else if (id.det() == DetId::HGCalTrigger &&
             HGCalTriggerDetId(id).subdet() == HGCalTriggerSubdetector::HFNoseTrigger) {
    em = HFNoseTriggerDetId(id).isEE();
  }
  return em;
}

References DetId::Forward, HFNose, HFNoseTrigger, DetId::HGCalEE, HGCalEETrigger, DetId::HGCalTrigger, HGCEE, HFNoseTriggerDetId::isEE(), and HFNoseDetId::isEE().

Referenced by HGCalConcentratorSuperTriggerCellImpl::assignSuperTriggerCellEnergyAndPosition(), HGCalTriggerNtupleHGCTriggerCells::fill(), HGCalHistoSeedingImpl::fillHistoClusters(), getSubDetectorType(), isHad(), HGCalVFEProcessorSums::run(), and HGCalConcentratorProcessorSelection::run().

isHad()

bool HGCalTriggerTools::isHad ( const DetId &  id ) const

inline

Definition at line 44 of file HGCalTriggerTools.h.

{ return !isEm(id); }

References isEm().

Referenced by layerWithOffset().

isNose()

bool HGCalTriggerTools::isNose

(

const DetId &

id

)

const

Definition at line 170 of file HGCalTriggerTools.cc.

                                                    {
  bool nose = false;
  if (id.det() == DetId::Forward && id.subdetId() == ForwardSubdetector::HFNose) {
    nose = true;
  } else if (id.det() == DetId::HGCalTrigger &&
             HGCalTriggerDetId(id).subdet() == HGCalTriggerSubdetector::HFNoseTrigger) {
    nose = true;
  }
  return nose;
}

References DetId::Forward, HFNose, HFNoseTrigger, and DetId::HGCalTrigger.

Referenced by HGCalConcentratorSuperTriggerCellImpl::assignSuperTriggerCellEnergyAndPosition(), HGCalTowerMap2DImpl::buildTowerMap2D(), HGCalClusteringDummyImpl::calibratePt(), HGCalShowerShape::coreShowerLength(), layerWithOffset(), HGCalShowerShape::percentileLayer(), and HGCalVFEProcessorSums::run().

isScintillator()

bool HGCalTriggerTools::isScintillator ( const DetId &  id ) const

inline

Definition at line 46 of file HGCalTriggerTools.h.

{ return !isSilicon(id); }

References isSilicon().

Referenced by HGCalCoarseTriggerCellMapping::getCoarseTriggerCellId(), HGCalCoarseTriggerCellMapping::getConstituentTriggerCells(), getSubDetectorType(), layerWithOffset(), HGCalConcentratorThresholdImpl::select(), thicknessIndex(), and HGCalVFESummationImpl::triggerCellSums().

isSilicon()

bool HGCalTriggerTools::isSilicon

(

const DetId &

id

)

const

Definition at line 181 of file HGCalTriggerTools.cc.

                                                       {
  bool silicon = false;
  if (id.det() == DetId::Forward) {
    silicon = (id.subdetId() != HGCHEB);
  } else if (id.det() == DetId::HGCalEE || id.det() == DetId::HGCalHSi) {
    silicon = true;
  } else if (id.det() == DetId::HGCalTrigger &&
             HGCalTriggerDetId(id).subdet() == HGCalTriggerSubdetector::HFNoseTrigger) {
    silicon = true;
  } else if (id.det() == DetId::HGCalTrigger &&
             (HGCalTriggerDetId(id).subdet() == HGCalTriggerSubdetector::HGCalEETrigger ||
              HGCalTriggerDetId(id).subdet() == HGCalTriggerSubdetector::HGCalHSiTrigger)) {
    silicon = true;
  }
  return silicon;
}

References DetId::Forward, HFNoseTrigger, DetId::HGCalEE, HGCalEETrigger, DetId::HGCalHSi, HGCalHSiTrigger, DetId::HGCalTrigger, HGCHEB, and HGCalTriggerDetId::subdet().

Referenced by HGCalConcentratorSuperTriggerCellImpl::assignSuperTriggerCellEnergyAndPosition(), HGCalTriggerCellCalibration::calibrateInMipT(), HGCalTriggerCellCalibration::calibrateMipTinGeV(), HGCalTriggerNtupleHGCTriggerCells::fill(), isScintillator(), layerWithOffset(), HGCalVFEProcessorSums::run(), and HGCalVFESummationImpl::triggerCellSums().

lastLayerBH()

unsigned HGCalTriggerTools::lastLayerBH ( ) const

inline

Definition at line 56 of file HGCalTriggerTools.h.

{ return totalLayers_; }

References totalLayers_.

Referenced by HGCalHistoClusteringImpl::eventSetup(), HGCalTriggerNtupleGen::fill(), and HGCalTowerMap2DImpl::newTowerMaps().

lastLayerEE()

unsigned HGCalTriggerTools::lastLayerEE ( ) const

inline

Definition at line 54 of file HGCalTriggerTools.h.

{ return eeLayers_; }

References eeLayers_.

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

lastLayerFH()

unsigned HGCalTriggerTools::lastLayerFH ( ) const

inline

Definition at line 55 of file HGCalTriggerTools.h.

{ return eeLayers_ + fhLayers_; }

References eeLayers_, and fhLayers_.

Referenced by getLayerZ().

lastLayerNose()

unsigned HGCalTriggerTools::lastLayerNose ( ) const

inline

Definition at line 57 of file HGCalTriggerTools.h.

{ return noseLayers_; }

References noseLayers_.

layer()

unsigned HGCalTriggerTools::layer

(

const DetId &

id

)

const

Definition at line 111 of file HGCalTriggerTools.cc.

                                                       {
  unsigned int layer = std::numeric_limits<unsigned int>::max();
  if (id.det() == DetId::Forward && id.subdetId() != ForwardSubdetector::HFNose) {
    layer = HGCalDetId(id).layer();
  } else if (id.det() == DetId::Forward && id.subdetId() == ForwardSubdetector::HFNose) {
    layer = HFNoseDetId(id).layer();
  } else if (id.det() == DetId::Hcal && id.subdetId() == HcalEndcap) {
    layer = HcalDetId(id).depth();
  } else if (id.det() == DetId::HGCalEE || id.det() == DetId::HGCalHSi) {
    layer = HGCSiliconDetId(id).layer();
  } else if (id.det() == DetId::HGCalTrigger &&
             (HGCalTriggerDetId(id).subdet() == HGCalTriggerSubdetector::HGCalEETrigger ||
              HGCalTriggerDetId(id).subdet() == HGCalTriggerSubdetector::HGCalHSiTrigger)) {
    layer = HGCalTriggerDetId(id).layer();
  } else if (id.det() == DetId::HGCalTrigger &&
             HGCalTriggerDetId(id).subdet() == HGCalTriggerSubdetector::HFNoseTrigger) {
    layer = HFNoseTriggerDetId(id).layer();
  } else if (id.det() == DetId::HGCalHSc) {
    layer = HGCScintillatorDetId(id).layer();
  }
  return layer;
}

References HcalDetId::depth(), DetId::Forward, DetId::Hcal, HcalEndcap, HFNose, HFNoseTrigger, DetId::HGCalEE, HGCalEETrigger, DetId::HGCalHSc, DetId::HGCalHSi, HGCalHSiTrigger, DetId::HGCalTrigger, HGCalDetId::layer(), HGCScintillatorDetId::layer(), HFNoseTriggerDetId::layer(), HGCalTriggerDetId::layer(), HFNoseDetId::layer(), HGCSiliconDetId::layer(), and SiStripPI::max.

Referenced by getLayerZ(), geometryXMLparser.DTAlignable::index(), geometryXMLparser.CSCAlignable::index(), HGCalClusteringImpl::isPertinent(), layerWithOffset(), and simToReco().

layers() [1/2]

unsigned HGCalTriggerTools::layers

(

DetId::Detector

type

)

const

Definition at line 89 of file HGCalTriggerTools.cc.

                                                           {
  unsigned layers = 0;
  switch (type) {
    case DetId::HGCalEE:
      layers = eeLayers_;
      break;
    case DetId::HGCalHSi:
      layers = fhLayers_;
      break;
    case DetId::HGCalHSc:
      layers = bhLayers_;
      break;
    // FIXME: to do HFNose
    case DetId::Forward:
      layers = totalLayers_;
      break;
    default:
      break;
  }
  return layers;
}

References bhLayers_, eeLayers_, fhLayers_, DetId::Forward, DetId::HGCalEE, DetId::HGCalHSc, DetId::HGCalHSi, layers(), and totalLayers_.

layers() [2/2]

unsigned HGCalTriggerTools::layers

(

ForwardSubdetector

type

)

const

Definition at line 65 of file HGCalTriggerTools.cc.

                                                                {
  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::HFNose:
      layers = noseLayers_;
      break;
    case ForwardSubdetector::ForwardEmpty:
      layers = totalLayers_;
      break;
    default:
      break;
  };
  return layers;
}

References bhLayers_, eeLayers_, fhLayers_, ForwardEmpty, HFNose, HGCEE, HGCHEB, HGCHEF, noseLayers_, and totalLayers_.

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

layerWithOffset()

unsigned HGCalTriggerTools::layerWithOffset

(

const DetId &

id

)

const

Definition at line 134 of file HGCalTriggerTools.cc.

                                                                 {
  unsigned int l = layer(id);
 
  if (isNose(id)) {
    l = layer(id);  // no offset for HFnose
  } else if (isHad(id) && isSilicon(id)) {
    l += eeLayers_;
  } else if (isHad(id) && isScintillator(id)) {
    if (geom_->isV9Geometry())
      l += eeLayers_;  // mixed silicon and scintillator layers
    else
      l += eeLayers_ + fhLayers_;
  }
 
  return l;
}

References eeLayers_, fhLayers_, geom_, isHad(), isNose(), isScintillator(), isSilicon(), HGCalTriggerGeometryBase::isV9Geometry(), cmsLHEtoEOSManager::l, and layer().

Referenced by HGCalTriggerValidator::analyze(), HGCalTowerMap2DImpl::buildTowerMap2D(), HGCalTriggerCellCalibration::calibrateMipTinGeV(), HGCalClusteringDummyImpl::calibratePt(), HGCalClusteringImpl::calibratePt(), HGCalHistoClusteringImpl::clusterSeedMulticluster(), HGCalShowerShape::eMax(), HGCalTriggerNtupleHGCClusters::fill(), HGCalTriggerNtupleHGCTriggerSums::fill(), HGCalTriggerNtupleHGCDigis::fill(), HGCalTriggerNtupleHGCTriggerCells::fill(), HGCalShowerShape::firstLayer(), HGCalCoarseTriggerCellMapping::getCoarseTriggerCellId(), HGCalCoarseTriggerCellMapping::getConstituentTriggerCells(), getLayerZ(), HGCalShowerShape::lastLayer(), HGCalShowerShape::maxLayer(), HGCalShowerShape::sigmaEtaEtaMax(), HGCalShowerShape::sigmaPhiPhiMax(), HGCalShowerShape::sigmaRRMax(), HGCalShowerShape::sigmaRRMean(), and HGCalClusteringImpl::triggerCellReshuffling().

sensorCellThicknessV8()

int HGCalTriggerTools::sensorCellThicknessV8 ( const DetId &  id ) const

private

Definition at line 379 of file HGCalTriggerTools.cc.

                                                                  {
  int thickness = 0;
  switch (id.subdetId()) {
    case ForwardSubdetector::HGCEE:
      thickness = geom_->eeTopology().dddConstants().waferTypeL(HGCalDetId(id).wafer()) - 1;
      break;
    case ForwardSubdetector::HGCHEF:
      thickness = geom_->fhTopology().dddConstants().waferTypeL(HGCalDetId(id).wafer()) - 1;
      break;
    default:
      break;
  };
  return thickness;
}

References HGCalTopology::dddConstants(), HGCalTriggerGeometryBase::eeTopology(), HGCalTriggerGeometryBase::fhTopology(), geom_, HGCEE, HGCHEF, Calorimetry_cff::thickness, and HGCalDDDConstants::waferTypeL().

Referenced by thicknessIndex().

simToReco() [1/2]

DetId HGCalTriggerTools::simToReco	(	const DetId &	simid,
		const HcalTopology &	topo
	)		const

Definition at line 370 of file HGCalTriggerTools.cc.

                                                                                     {
  DetId recoid(0);
  const auto& dddConst = topo.dddConstants();
  HcalDetId id = HcalHitRelabeller::relabel(simid, dddConst);
  if (id.subdet() == int(HcalEndcap))
    recoid = id;
  return recoid;
}

References HcalTopology::dddConstants(), HcalEndcap, triggerObjects_cff::id, and HcalHitRelabeller::relabel().

simToReco() [2/2]

DetId HGCalTriggerTools::simToReco	(	const DetId &	simid,
		const HGCalTopology &	topo
	)		const

Definition at line 346 of file HGCalTriggerTools.cc.

                                                                                      {
  DetId recoid(0);
  const auto& dddConst = topo.dddConstants();
  // V9
  if (dddConst.waferHexagon8() || dddConst.tileTrapezoid()) {
    recoid = simid;
  }
  // V8
  else {
    int subdet(simid.subdetId());
    int layer = 0, cell = 0, sec = 0, subsec = 0, zp = 0;
    HGCalTestNumbering::unpackHexagonIndex(simid, subdet, zp, layer, sec, subsec, cell);
    // sec is wafer and subsec is celltype
    // skip this hit if after ganging it is not valid
    auto recoLayerCell = dddConst.simToReco(cell, layer, sec, topo.detectorType());
    cell = recoLayerCell.first;
    layer = recoLayerCell.second;
    if (layer >= 0 && cell >= 0) {
      recoid = HGCalDetId((ForwardSubdetector)subdet, zp, layer, subsec, sec, cell);
    }
  }
  return recoid;
}

References HGCalTopology::dddConstants(), HGCalTopology::detectorType(), layer(), fileinputsource_cfi::sec, DetId::subdetId(), and HGCalTestNumbering::unpackHexagonIndex().

Referenced by CaloTruthCellsProducer::makeHitMap(), HGCalTriggerNtupleHGCDigis::simhits(), and HGCalTriggerNtupleHGCTriggerCells::simhits().

thicknessIndex()

int HGCalTriggerTools::thicknessIndex	(	const DetId &	id,
		bool	tc = false
	)		const

Definition at line 233 of file HGCalTriggerTools.cc.

                                                                    {
  if (isScintillator(id)) {
    return kScintillatorPseudoThicknessIndex_;
  }
  unsigned det = id.det();
  int thickness = 0;
  // For the v8 detid scheme
  if (det == DetId::Forward && id.subdetId() != ForwardSubdetector::HFNose) {
    if (!tc)
      thickness = sensorCellThicknessV8(id);
    else {
      // For the old geometry, TCs can contain sensor cells
      // with different thicknesses.
      // Use a majority logic to find the TC thickness
      std::array<unsigned, 3> occurences = {{0, 0, 0}};
      for (const auto& c_id : geom_->getCellsFromTriggerCell(id)) {
        unsigned c_det = DetId(c_id).det();
        int c_thickness = -1;
        // For the v8 detid scheme
        if (c_det == DetId::Forward) {
          c_thickness = sensorCellThicknessV8(c_id);
        } else {
          c_thickness = HGCSiliconDetId(c_id).type();
        }
        if (c_thickness < 0 || unsigned(c_thickness) >= occurences.size()) {
          throw cms::Exception("OutOfBound") << "Found thickness index = " << c_thickness;
        }
        occurences[c_thickness]++;
      }
      thickness = std::max_element(occurences.begin(), occurences.end()) - occurences.begin();
    }
  }
  // For the v9 detid scheme
  else if (det == DetId::HGCalEE || det == DetId::HGCalHSi) {
    thickness = HGCSiliconDetId(id).type();
  } else if (det == DetId::Forward && id.subdetId() == ForwardSubdetector::HFNose) {
    thickness = HFNoseDetId(id).type();
  } else if (id.det() == DetId::HGCalTrigger &&
             (HGCalTriggerDetId(id).subdet() == HGCalTriggerSubdetector::HGCalEETrigger ||
              HGCalTriggerDetId(id).subdet() == HGCalTriggerSubdetector::HGCalHSiTrigger)) {
    thickness = HGCalTriggerDetId(id).type();
  } else if (id.det() == DetId::HGCalTrigger &&
             HGCalTriggerDetId(id).subdet() == HGCalTriggerSubdetector::HFNoseTrigger) {
    thickness = HFNoseTriggerDetId(id).type();
  }
  return thickness;
}

References DetId::det(), Exception, DetId::Forward, geom_, HGCalTriggerGeometryBase::getCellsFromTriggerCell(), HFNose, HFNoseTrigger, DetId::HGCalEE, HGCalEETrigger, DetId::HGCalHSi, HGCalHSiTrigger, DetId::HGCalTrigger, isScintillator(), kScintillatorPseudoThicknessIndex_, mps_monitormerge::occurences, sensorCellThicknessV8(), Calorimetry_cff::thickness, HGCalTriggerDetId::type(), HFNoseTriggerDetId::type(), HFNoseDetId::type(), and HGCSiliconDetId::type().

Referenced by HGCalConcentratorSuperTriggerCellImpl::assignSuperTriggerCellEnergyAndPosition(), HGCalTriggerCellCalibration::calibrateInMipT(), HGCalTriggerCellCalibration::calibrateMipTinGeV(), HGCalConcentratorCoarsenerImpl::coarsen(), HGCalConcentratorSuperTriggerCellImpl::createAllTriggerCells(), HGCalTriggerNtupleHGCTriggerCells::fill(), HGCalCoarseTriggerCellMapping::getCoarseTriggerCellId(), HGCalCoarseTriggerCellMapping::getConstituentTriggerCells(), HGCalVFEProcessorSums::run(), HGCalConcentratorProcessorSelection::run(), and HGCalVFESummationImpl::triggerCellSums().

triggerLayer()

unsigned HGCalTriggerTools::triggerLayer ( const unsigned  id ) const

inline

Definition at line 84 of file HGCalTriggerTools.h.

{ return geom_->triggerLayer(id); }

References geom_, and HGCalTriggerGeometryBase::triggerLayer().

Referenced by HGCalTriggerClusterInterpretationEM::interpret().

zside()

int HGCalTriggerTools::zside

(

const DetId &

id

)

const

Definition at line 210 of file HGCalTriggerTools.cc.

                                                  {
  int zside = 0;
  if (id.det() == DetId::Forward && id.subdetId() != ForwardSubdetector::HFNose) {
    zside = HGCalDetId(id).zside();
  } else if (id.det() == DetId::Forward && id.subdetId() == ForwardSubdetector::HFNose) {
    zside = HFNoseDetId(id).zside();
  } else if (id.det() == DetId::Hcal && id.subdetId() == HcalEndcap) {
    zside = HcalDetId(id).zside();
  } else if (id.det() == DetId::HGCalEE || id.det() == DetId::HGCalHSi) {
    zside = HGCSiliconDetId(id).zside();
  } else if (id.det() == DetId::HGCalTrigger &&
             (HGCalTriggerDetId(id).subdet() == HGCalTriggerSubdetector::HGCalEETrigger ||
              HGCalTriggerDetId(id).subdet() == HGCalTriggerSubdetector::HGCalHSiTrigger)) {
    zside = HGCalTriggerDetId(id).zside();
  } else if (id.det() == DetId::HGCalTrigger &&
             HGCalTriggerDetId(id).subdet() == HGCalTriggerSubdetector::HFNoseTrigger) {
    zside = HFNoseTriggerDetId(id).zside();
  } else if (id.det() == DetId::HGCalHSc) {
    zside = HGCScintillatorDetId(id).zside();
  }
  return zside;
}

References DetId::Forward, DetId::Hcal, HcalEndcap, HFNose, HFNoseTrigger, DetId::HGCalEE, HGCalEETrigger, DetId::HGCalHSc, DetId::HGCalHSi, HGCalHSiTrigger, DetId::HGCalTrigger, HGCScintillatorDetId::zside(), HGCalDetId::zside(), HGCalTriggerDetId::zside(), HFNoseTriggerDetId::zside(), HFNoseDetId::zside(), HGCSiliconDetId::zside(), and HcalDetId::zside().

Referenced by HGCalMulticlusteringImpl::clusterizeDBSCAN(), HGCalHistoClusteringImpl::clusterSeedMulticluster(), HGCalTriggerNtupleHGCTriggerSums::fill(), HGCalTriggerNtupleHGCDigis::fill(), HGCalTriggerNtupleHGCTriggerCells::fill(), HGCalHistoSeedingImpl::fillHistoClusters(), HGCalTriggerTowerGeometryHelper::HGCalTriggerTowerGeometryHelper(), HGCalMulticlusteringImpl::isPertinent(), HGCalClusteringImpl::isPertinent(), and HGCalClusteringImpl::triggerCellReshuffling().

Member Data Documentation

bhLayers_

unsigned HGCalTriggerTools::bhLayers_

private

Definition at line 99 of file HGCalTriggerTools.h.

Referenced by eventSetup(), and layers().

eeLayers_

unsigned HGCalTriggerTools::eeLayers_

private

Definition at line 97 of file HGCalTriggerTools.h.

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

fhLayers_

unsigned HGCalTriggerTools::fhLayers_

private

Definition at line 98 of file HGCalTriggerTools.h.

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

geom_

const HGCalTriggerGeometryBase* HGCalTriggerTools::geom_

private

Definition at line 96 of file HGCalTriggerTools.h.

Referenced by eventSetup(), getLayerZ(), getTCPosition(), getTriggerGeometry(), layerWithOffset(), sensorCellThicknessV8(), thicknessIndex(), and triggerLayer().

kScintillatorPseudoThicknessIndex_

constexpr unsigned HGCalTriggerTools::kScintillatorPseudoThicknessIndex_ = 3

staticconstexpr

Definition at line 86 of file HGCalTriggerTools.h.

Referenced by thicknessIndex().

noseLayers_

unsigned HGCalTriggerTools::noseLayers_

private

Definition at line 100 of file HGCalTriggerTools.h.

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

totalLayers_

unsigned HGCalTriggerTools::totalLayers_

private

Definition at line 101 of file HGCalTriggerTools.h.

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