CMS 3D CMS Logo

List of all members | Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Private Member Functions | Private Attributes | Friends
HcalGeometry Class Reference

#include <HcalGeometry.h>

Inheritance diagram for HcalGeometry:
CaloSubdetectorGeometry

Public Types

enum  { k_NumberOfParametersPerShape = 5 }
 
typedef HcalGeometryRecord AlignedRecord
 
typedef HcalAlignmentRcd AlignmentRecord
 
typedef CaloCellGeometry::CCGFloat CCGFloat
 
typedef HcalDetId DetIdType
 
typedef std::vector< IdealObliquePrismHBCellVec
 
typedef std::vector< IdealObliquePrismHECellVec
 
typedef std::vector< IdealZPrismHFCellVec
 
typedef std::vector< IdealObliquePrismHOCellVec
 
typedef PHcalRcd PGeometryRecord
 
typedef CaloCellGeometry::Pt3D Pt3D
 
typedef CaloCellGeometry::Pt3DVec Pt3DVec
 
- Public Types inherited from CaloSubdetectorGeometry
typedef CaloCellGeometry::CCGFloat CCGFloat
 
typedef std::vector< std::shared_ptr< const CaloCellGeometry > > CellSet
 
typedef std::set< DetIdDetIdSet
 
typedef std::vector< CCGFloatDimVec
 
typedef std::vector< unsigned int > IVec
 
typedef CaloCellGeometry::ParMgr ParMgr
 
typedef CaloCellGeometry::ParVec ParVec
 
typedef CaloCellGeometry::ParVecVec ParVecVec
 
typedef std::vector< CCGFloatTrVec
 

Public Member Functions

GlobalPoint getBackPosition (const DetId &id) const
 
CaloSubdetectorGeometry::DetIdSet getCells (const GlobalPoint &r, double dR) const override
 Get a list of all cells within a dR of the given cell. More...
 
DetId getClosestCell (const GlobalPoint &r) const override
 
CaloCellGeometry::CornersVec getCorners (const DetId &id) const
 
std::shared_ptr< const CaloCellGeometrygetGeometry (const DetId &id) const override
 Get the cell geometry of a given detector id. Should return false if not found. More...
 
unsigned int getHxSize (const int type) const
 
GlobalPoint getPosition (const DetId &id) const
 
void getSummary (CaloSubdetectorGeometry::TrVec &trVector, CaloSubdetectorGeometry::IVec &iVector, CaloSubdetectorGeometry::DimVec &dimVector, CaloSubdetectorGeometry::IVec &dinsVector) const override
 
const std::vector< DetId > & getValidDetIds (DetId::Detector det=DetId::Detector(0), int subdet=0) const override
 Get a list of valid detector ids (for the given subdetector) More...
 
 HcalGeometry (const HcalTopology &topology)
 
void localCorners (Pt3DVec &lc, const CCGFloat *pv, unsigned int i, Pt3D &ref)
 
void newCell (const GlobalPoint &f1, const GlobalPoint &f2, const GlobalPoint &f3, const CCGFloat *parm, const DetId &detId) override
 
unsigned int numberOfParametersPerShape () const override
 
unsigned int numberOfShapes () const override
 
const HcalTopologytopology () const
 
 ~HcalGeometry () override
 The HcalGeometry will delete all its cell geometries at destruction time. More...
 
- Public Member Functions inherited from CaloSubdetectorGeometry
void allocateCorners (CaloCellGeometry::CornersVec::size_type n)
 
void allocatePar (ParVec::size_type n, unsigned int m)
 
 CaloSubdetectorGeometry ()
 
 CaloSubdetectorGeometry (const CaloSubdetectorGeometry &)=delete
 avoid copies More...
 
CaloCellGeometry::CornersMgrcornersMgr ()
 
CCGFloat deltaEta (const DetId &detId) const
 
CCGFloat deltaPhi (const DetId &detId) const
 
virtual void fillDefaultNamedParameters () const
 
virtual CellSet getCellSet (const GlobalPoint &r, double dR) const
 
virtual void initializeParms ()
 
virtual unsigned int numberOfTransformParms () const
 
CaloSubdetectorGeometryoperator= (const CaloSubdetectorGeometry &)=delete
 
ParMgrparMgr ()
 
const ParMgrparMgrConst () const
 
ParVecVecparVecVec ()
 
const ParVecVecparVecVec () const
 
virtual bool present (const DetId &id) const
 is this detid present in the geometry? More...
 
virtual ~CaloSubdetectorGeometry ()
 The base class DOES assume that it owns the CaloCellGeometry objects. More...
 

Static Public Member Functions

static unsigned int alignmentBarEndForIndexLocal (const DetId &id, unsigned int nD)
 
static unsigned int alignmentBarrelIndexLocal (const DetId &id)
 
static unsigned int alignmentEndcapIndexLocal (const DetId &id)
 
static unsigned int alignmentForwardIndexLocal (const DetId &id)
 
static unsigned int alignmentOuterIndexLocal (const DetId &id)
 
static unsigned int alignmentTransformIndexGlobal (const DetId &id)
 
static unsigned int alignmentTransformIndexLocal (const DetId &id)
 
static std::string dbString ()
 
static DetId detIdFromBarrelAlignmentIndex (unsigned int i)
 
static DetId detIdFromEndcapAlignmentIndex (unsigned int i)
 
static DetId detIdFromForwardAlignmentIndex (unsigned int i)
 
static DetId detIdFromLocalAlignmentIndex (unsigned int i)
 
static DetId detIdFromOuterAlignmentIndex (unsigned int i)
 
static unsigned int numberOfAlignments ()
 
static unsigned int numberOfBarrelAlignments ()
 
static unsigned int numberOfEndcapAlignments ()
 
static unsigned int numberOfForwardAlignments ()
 
static unsigned int numberOfOuterAlignments ()
 
static std::string producerTag ()
 

Protected Member Functions

const CaloCellGeometrygetGeometryRawPtr (uint32_t index) const override
 
unsigned int indexFor (const DetId &id) const override
 
unsigned int sizeForDenseIndex (const DetId &id) const override
 
- Protected Member Functions inherited from CaloSubdetectorGeometry
void addValidID (const DetId &id)
 
virtual std::shared_ptr< const CaloCellGeometrycellGeomPtr (uint32_t index) const
 

Private Member Functions

DetId correctId (const DetId &id) const
 
int etaRing (HcalSubdetector bc, double abseta) const
 helper methods for getClosestCell More...
 
void fillDetIds () const
 
std::shared_ptr< const CaloCellGeometrygetGeometryBase (const DetId &id) const
 
void increaseReserve (unsigned int extra)
 
void init ()
 
void newCellFast (const GlobalPoint &f1, const GlobalPoint &f2, const GlobalPoint &f3, const CCGFloat *parm, const DetId &detId)
 
unsigned int newCellImpl (const GlobalPoint &f1, const GlobalPoint &f2, const GlobalPoint &f3, const CCGFloat *parm, const DetId &detId)
 
int phiBin (HcalSubdetector bc, int etaring, double phi) const
 
void sortValidIds ()
 

Private Attributes

CaloSubdetectorGeometry::IVec m_dins
 
edm::AtomicPtrCache< std::vector< DetId > > m_emptyIds
 
HBCellVec m_hbCellVec
 
edm::AtomicPtrCache< std::vector< DetId > > m_hbIds
 
HECellVec m_heCellVec
 
edm::AtomicPtrCache< std::vector< DetId > > m_heIds
 
HFCellVec m_hfCellVec
 
edm::AtomicPtrCache< std::vector< DetId > > m_hfIds
 
HOCellVec m_hoCellVec
 
edm::AtomicPtrCache< std::vector< DetId > > m_hoIds
 
bool m_mergePosition
 
const HcalTopologym_topology
 

Friends

class HcalFlexiHardcodeGeometryLoader
 
class HcalGeometryPlan1Tester
 
class HcalHardcodeGeometryLoader
 

Additional Inherited Members

- Static Protected Member Functions inherited from CaloSubdetectorGeometry
static CCGFloat deltaR (const GlobalPoint &p1, const GlobalPoint &p2)
 
- Protected Attributes inherited from CaloSubdetectorGeometry
ParVecVec m_parVecVec
 
std::vector< DetIdm_validIds
 

Detailed Description

Definition at line 17 of file HcalGeometry.h.

Member Typedef Documentation

Definition at line 36 of file HcalGeometry.h.

Definition at line 35 of file HcalGeometry.h.

Definition at line 31 of file HcalGeometry.h.

Definition at line 38 of file HcalGeometry.h.

Definition at line 26 of file HcalGeometry.h.

Definition at line 27 of file HcalGeometry.h.

typedef std::vector<IdealZPrism> HcalGeometry::HFCellVec

Definition at line 29 of file HcalGeometry.h.

Definition at line 28 of file HcalGeometry.h.

Definition at line 37 of file HcalGeometry.h.

Definition at line 32 of file HcalGeometry.h.

Definition at line 33 of file HcalGeometry.h.

Member Enumeration Documentation

anonymous enum
Enumerator
k_NumberOfParametersPerShape 

Definition at line 41 of file HcalGeometry.h.

Constructor & Destructor Documentation

HcalGeometry::HcalGeometry ( const HcalTopology topology)
explicit

Definition at line 16 of file HcalGeometry.cc.

References init().

Referenced by numberOfParametersPerShape().

16  :
17  m_topology(topology), m_mergePosition(topology.getMergePositionFlag()) {
18  init();
19 }
const HcalTopology & m_topology
Definition: HcalGeometry.h:159
bool m_mergePosition
Definition: HcalGeometry.h:160
bool getMergePositionFlag() const
Definition: HcalTopology.h:170
HcalGeometry::~HcalGeometry ( )
override

The HcalGeometry will delete all its cell geometries at destruction time.

Definition at line 21 of file HcalGeometry.cc.

Referenced by numberOfParametersPerShape().

21 {}

Member Function Documentation

unsigned int HcalGeometry::alignmentBarEndForIndexLocal ( const DetId id,
unsigned int  nD 
)
static

Definition at line 331 of file HcalGeometry.cc.

References HcalDetId::ieta(), and HcalDetId::iphi().

Referenced by alignmentBarrelIndexLocal(), alignmentEndcapIndexLocal(), alignmentForwardIndexLocal(), and numberOfAlignments().

332  {
333  const HcalDetId hid ( id ) ;
334  const unsigned int iphi ( hid.iphi() ) ;
335  const int ieta ( hid.ieta() ) ;
336  const unsigned int index ( ( 0 < ieta ? nD/2 : 0 ) + ( iphi + 1 )%72/4 ) ;
337  assert( index < nD ) ;
338  return index ;
339 }
unsigned int HcalGeometry::alignmentBarrelIndexLocal ( const DetId id)
static

Definition at line 341 of file HcalGeometry.cc.

References alignmentBarEndForIndexLocal(), and numberOfBarrelAlignments().

Referenced by alignmentTransformIndexLocal(), and numberOfAlignments().

341  {
343 }
static unsigned int numberOfBarrelAlignments()
Definition: HcalGeometry.h:68
static unsigned int alignmentBarEndForIndexLocal(const DetId &id, unsigned int nD)
unsigned int HcalGeometry::alignmentEndcapIndexLocal ( const DetId id)
static

Definition at line 345 of file HcalGeometry.cc.

References alignmentBarEndForIndexLocal(), and numberOfEndcapAlignments().

Referenced by alignmentTransformIndexLocal(), and numberOfAlignments().

345  {
347 }
static unsigned int numberOfEndcapAlignments()
Definition: HcalGeometry.h:70
static unsigned int alignmentBarEndForIndexLocal(const DetId &id, unsigned int nD)
unsigned int HcalGeometry::alignmentForwardIndexLocal ( const DetId id)
static

Definition at line 349 of file HcalGeometry.cc.

References alignmentBarEndForIndexLocal(), and numberOfForwardAlignments().

Referenced by alignmentTransformIndexLocal(), and numberOfAlignments().

349  {
351 }
static unsigned int numberOfForwardAlignments()
Definition: HcalGeometry.h:72
static unsigned int alignmentBarEndForIndexLocal(const DetId &id, unsigned int nD)
unsigned int HcalGeometry::alignmentOuterIndexLocal ( const DetId id)
static

Definition at line 353 of file HcalGeometry.cc.

References HcalDetId::ieta(), HcalDetId::iphi(), numberOfOuterAlignments(), and relativeConstraints::ring.

Referenced by alignmentTransformIndexLocal(), and numberOfAlignments().

353  {
354  const HcalDetId hid ( id ) ;
355  const int ieta ( hid.ieta() ) ;
356  const int iphi ( hid.iphi() ) ;
357  const int ring ( ieta < -10 ? 0 :
358  ( ieta < -4 ? 1 :
359  ( ieta < 5 ? 2 :
360  ( ieta < 11 ? 3 : 4 ) ) ) ) ;
361 
362  const unsigned int index ( 12*ring + ( iphi - 1 )/6 ) ;
363  assert( index < numberOfOuterAlignments() ) ;
364  return index ;
365 }
static unsigned int numberOfOuterAlignments()
Definition: HcalGeometry.h:74
unsigned int HcalGeometry::alignmentTransformIndexGlobal ( const DetId id)
static

Definition at line 391 of file HcalGeometry.cc.

References DetId::Hcal.

Referenced by numberOfAlignments(), and CaloGeometryDBEP< T, U >::produceAligned().

391  {
392  return (unsigned int)DetId::Hcal - 1 ;
393 }
unsigned int HcalGeometry::alignmentTransformIndexLocal ( const DetId id)
static

Definition at line 367 of file HcalGeometry.cc.

References alignmentBarrelIndexLocal(), alignmentEndcapIndexLocal(), alignmentForwardIndexLocal(), alignmentOuterIndexLocal(), DetId::Hcal, HcalBarrel, HcalEndcap, HcalForward, HcalObjRepresent::isHB(), HcalObjRepresent::isHE(), HcalObjRepresent::isHF(), numberOfAlignments(), numberOfBarrelAlignments(), numberOfEndcapAlignments(), numberOfForwardAlignments(), and HcalDetId::subdet().

Referenced by numberOfAlignments(), and CaloGeometryDBEP< T, U >::produceAligned().

367  {
368  assert(id.det() == DetId::Hcal) ;
369 
370  const HcalDetId hid ( id ) ;
371  bool isHB = (hid.subdet() == HcalBarrel);
372  bool isHE = (hid.subdet() == HcalEndcap);
373  bool isHF = (hid.subdet() == HcalForward);
374  // bool isHO = (hid.subdet() == HcalOuter);
375 
376  const unsigned int nB ( numberOfBarrelAlignments() ) ;
377  const unsigned int nE ( numberOfEndcapAlignments() ) ;
378  const unsigned int nF ( numberOfForwardAlignments() ) ;
379  // const unsigned int nO ( numberOfOuterAlignments() ) ;
380 
381  const unsigned int index (isHB ? alignmentBarrelIndexLocal(id) :
382  isHE ? alignmentEndcapIndexLocal(id) + nB :
383  isHF ? alignmentForwardIndexLocal( id ) + nB + nE :
384  alignmentOuterIndexLocal(id) + nB + nE + nF
385  );
386 
387  assert( index < numberOfAlignments() ) ;
388  return index ;
389 }
static unsigned int numberOfBarrelAlignments()
Definition: HcalGeometry.h:68
static unsigned int alignmentEndcapIndexLocal(const DetId &id)
bool isHE(int etabin, int depth)
bool isHB(int etabin, int depth)
static unsigned int alignmentBarrelIndexLocal(const DetId &id)
static unsigned int alignmentOuterIndexLocal(const DetId &id)
static unsigned int numberOfEndcapAlignments()
Definition: HcalGeometry.h:70
bool isHF(int etabin, int depth)
static unsigned int numberOfForwardAlignments()
Definition: HcalGeometry.h:72
static unsigned int alignmentForwardIndexLocal(const DetId &id)
static unsigned int numberOfAlignments()
Definition: HcalGeometry.h:78
DetId HcalGeometry::correctId ( const DetId id) const
private

Definition at line 561 of file HcalGeometry.cc.

References triggerObjects_cff::id, m_mergePosition, m_topology, and HcalTopology::mergedDepthDetId().

Referenced by getClosestCell(), and getGeometryBase().

561  {
562 
563  if (m_mergePosition) {
564  HcalDetId hid(id);
565  return ((DetId)(m_topology.mergedDepthDetId(hid)));
566  } else {
567  return id;
568  }
569 }
const HcalTopology & m_topology
Definition: HcalGeometry.h:159
bool m_mergePosition
Definition: HcalGeometry.h:160
HcalDetId mergedDepthDetId(const HcalDetId &id) const
Definition: HcalTopology.h:169
Definition: DetId.h:18
static std::string HcalGeometry::dbString ( )
inlinestatic

Definition at line 43 of file HcalGeometry.h.

Referenced by CaloGeometryDBEP< T, U >::produceAligned().

43 { return "PHcalRcd" ; }
DetId HcalGeometry::detIdFromBarrelAlignmentIndex ( unsigned int  i)
static

Definition at line 285 of file HcalGeometry.cc.

References HcalBarrel, and numberOfBarrelAlignments().

Referenced by detIdFromLocalAlignmentIndex(), numberOfAlignments(), FakeCaloAlignmentEP::produceHBAli(), and TestCaloAlignmentEP::produceHBAli().

285  {
286  assert( i < numberOfBarrelAlignments() ) ;
287  const int ieta ( i < numberOfBarrelAlignments()/2 ? -1 : 1 ) ;
288  const int iphi ( 1 + (4*i)%72 ) ;
289  return HcalDetId( HcalBarrel, ieta, iphi, 1 ) ;
290 }
static unsigned int numberOfBarrelAlignments()
Definition: HcalGeometry.h:68
DetId HcalGeometry::detIdFromEndcapAlignmentIndex ( unsigned int  i)
static

Definition at line 292 of file HcalGeometry.cc.

References HcalEndcap, and numberOfEndcapAlignments().

Referenced by detIdFromLocalAlignmentIndex(), numberOfAlignments(), FakeCaloAlignmentEP::produceHEAli(), and TestCaloAlignmentEP::produceHEAli().

292  {
293  assert( i < numberOfEndcapAlignments() ) ;
294  const int ieta ( i < numberOfEndcapAlignments()/2 ? -16 : 16 ) ;
295  const int iphi ( 1 + (4*i)%72 ) ;
296  return HcalDetId( HcalEndcap, ieta, iphi, 1 ) ;
297 }
static unsigned int numberOfEndcapAlignments()
Definition: HcalGeometry.h:70
DetId HcalGeometry::detIdFromForwardAlignmentIndex ( unsigned int  i)
static

Definition at line 299 of file HcalGeometry.cc.

References HcalForward, and numberOfForwardAlignments().

Referenced by detIdFromLocalAlignmentIndex(), numberOfAlignments(), FakeCaloAlignmentEP::produceHFAli(), and TestCaloAlignmentEP::produceHFAli().

299  {
300  assert( i < numberOfForwardAlignments() ) ;
301  const int ieta ( i < numberOfForwardAlignments()/2 ? -29 : 29 ) ;
302  const int iphi ( 1 + (4*i)%72 ) ;
303  return HcalDetId( HcalForward, ieta, iphi, 1 ) ;
304 }
static unsigned int numberOfForwardAlignments()
Definition: HcalGeometry.h:72
DetId HcalGeometry::detIdFromLocalAlignmentIndex ( unsigned int  i)
static

Definition at line 317 of file HcalGeometry.cc.

References detIdFromBarrelAlignmentIndex(), detIdFromEndcapAlignmentIndex(), detIdFromForwardAlignmentIndex(), detIdFromOuterAlignmentIndex(), numberOfAlignments(), numberOfBarrelAlignments(), numberOfEndcapAlignments(), and numberOfForwardAlignments().

Referenced by numberOfAlignments().

317  {
318  assert( i < numberOfAlignments() ) ;
319 
320  const unsigned int nB ( numberOfBarrelAlignments() ) ;
321  const unsigned int nE ( numberOfEndcapAlignments() ) ;
322  const unsigned int nF ( numberOfForwardAlignments() ) ;
323 // const unsigned int nO ( numberOfOuterAlignments() ) ;
324 
325  return ( i < nB ? detIdFromBarrelAlignmentIndex( i ) :
326  i < nB+nE ? detIdFromEndcapAlignmentIndex( i - nB ) :
327  i < nB+nE+nF ? detIdFromForwardAlignmentIndex( i - nB - nE ) :
328  detIdFromOuterAlignmentIndex( i - nB - nE - nF ) ) ;
329 }
static unsigned int numberOfBarrelAlignments()
Definition: HcalGeometry.h:68
static DetId detIdFromOuterAlignmentIndex(unsigned int i)
static DetId detIdFromForwardAlignmentIndex(unsigned int i)
static DetId detIdFromBarrelAlignmentIndex(unsigned int i)
static unsigned int numberOfEndcapAlignments()
Definition: HcalGeometry.h:70
static unsigned int numberOfForwardAlignments()
Definition: HcalGeometry.h:72
static DetId detIdFromEndcapAlignmentIndex(unsigned int i)
static unsigned int numberOfAlignments()
Definition: HcalGeometry.h:78
DetId HcalGeometry::detIdFromOuterAlignmentIndex ( unsigned int  i)
static

Definition at line 306 of file HcalGeometry.cc.

References HcalOuter, numberOfOuterAlignments(), and relativeConstraints::ring.

Referenced by detIdFromLocalAlignmentIndex(), numberOfAlignments(), FakeCaloAlignmentEP::produceHOAli(), and TestCaloAlignmentEP::produceHOAli().

306  {
307  assert( i < numberOfOuterAlignments() ) ;
308  const int ring ( i/12 ) ;
309  const int ieta ( 0 == ring ? -11 :
310  1 == ring ? -5 :
311  2 == ring ? 1 :
312  3 == ring ? 5 : 11 ) ;
313  const int iphi ( 1 + ( i - ring*12 )*6 ) ;
314  return HcalDetId( HcalOuter, ieta, iphi, 4 ) ;
315 }
static unsigned int numberOfOuterAlignments()
Definition: HcalGeometry.h:74
int HcalGeometry::etaRing ( HcalSubdetector  bc,
double  abseta 
) const
private

helper methods for getClosestCell

Definition at line 206 of file HcalGeometry.cc.

References HcalTopology::etaRing(), and m_topology.

Referenced by getCells(), getClosestCell(), and getGeometryBase().

206  {
207  return m_topology.etaRing(bc, abseta);
208 }
const HcalTopology & m_topology
Definition: HcalGeometry.h:159
int etaRing(HcalSubdetector subdet, double eta) const
eta and phi index from eta, phi values
void HcalGeometry::fillDetIds ( ) const
private

Definition at line 39 of file HcalGeometry.cc.

References CaloSubdetectorGeometry::getValidDetIds(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, mps_fire::i, triggerObjects_cff::id, edm::AtomicPtrCache< T >::isSet(), m_emptyIds, m_hbIds, m_heIds, m_hfIds, m_hoIds, eostools::move(), and edm::AtomicPtrCache< T >::set().

Referenced by getGeometryBase(), getHxSize(), and getValidDetIds().

39  {
40  // this must test the last record filled to avoid a race condition
41  if(!m_emptyIds.isSet()) {
42  std::unique_ptr<std::vector<DetId>> p_hbIds{new std::vector<DetId>};
43  std::unique_ptr<std::vector<DetId>> p_heIds{new std::vector<DetId>};
44  std::unique_ptr<std::vector<DetId>> p_hoIds{new std::vector<DetId>};
45  std::unique_ptr<std::vector<DetId>> p_hfIds{new std::vector<DetId>};
46  std::unique_ptr<std::vector<DetId>> p_emptyIds{new std::vector<DetId>};
47 
48  const std::vector<DetId>& baseIds (CaloSubdetectorGeometry::getValidDetIds());
49  for (unsigned int i ( 0 ) ; i != baseIds.size() ; ++i) {
50  const DetId id ( baseIds[i] );
51  if (id.subdetId() == HcalBarrel) {
52  p_hbIds->emplace_back( id ) ;
53  } else if (id.subdetId() == HcalEndcap) {
54  p_heIds->emplace_back( id ) ;
55  } else if (id.subdetId() == HcalOuter) {
56  p_hoIds->emplace_back( id ) ;
57  } else if (id.subdetId() == HcalForward) {
58  p_hfIds->emplace_back( id ) ;
59  }
60  }
61  std::sort( p_hbIds->begin(), p_hbIds->end() ) ;
62  std::sort( p_heIds->begin(), p_heIds->end() ) ;
63  std::sort( p_hoIds->begin(), p_hoIds->end() ) ;
64  std::sort( p_hfIds->begin(), p_hfIds->end() ) ;
65  p_emptyIds->resize( 0 ) ;
66 
67  m_hbIds.set(std::move(p_hbIds));
68  m_heIds.set(std::move(p_heIds));
69  m_hoIds.set(std::move(p_hoIds));
70  m_hfIds.set(std::move(p_hfIds));
71  m_emptyIds.set(std::move(p_emptyIds));
72  }
73 }
edm::AtomicPtrCache< std::vector< DetId > > m_heIds
Definition: HcalGeometry.h:163
edm::AtomicPtrCache< std::vector< DetId > > m_emptyIds
Definition: HcalGeometry.h:166
bool isSet() const
virtual const std::vector< DetId > & getValidDetIds(DetId::Detector det=DetId::Detector(0), int subdet=0) const
Get a list of valid detector ids (for the given subdetector)
Definition: DetId.h:18
bool set(std::unique_ptr< T > iNewValue) const
edm::AtomicPtrCache< std::vector< DetId > > m_hfIds
Definition: HcalGeometry.h:165
edm::AtomicPtrCache< std::vector< DetId > > m_hoIds
Definition: HcalGeometry.h:164
edm::AtomicPtrCache< std::vector< DetId > > m_hbIds
Definition: HcalGeometry.h:162
def move(src, dest)
Definition: eostools.py:510
GlobalPoint HcalGeometry::getBackPosition ( const DetId id) const

Definition at line 179 of file HcalGeometry.cc.

References getGeometryBase(), electrons_cff::ids, m_mergePosition, m_topology, and HcalTopology::unmergeDepthDetId().

Referenced by numberOfParametersPerShape().

179  {
180  if (!m_mergePosition) {
181  return (getGeometryBase(id)->getBackPoint());
182  } else {
183  std::vector<HcalDetId> ids;
185  return (getGeometryBase(ids.back())->getBackPoint());
186  }
187 }
const HcalTopology & m_topology
Definition: HcalGeometry.h:159
bool m_mergePosition
Definition: HcalGeometry.h:160
std::shared_ptr< const CaloCellGeometry > getGeometryBase(const DetId &id) const
Definition: HcalGeometry.h:129
void unmergeDepthDetId(const HcalDetId &id, std::vector< HcalDetId > &ids) const
Definition: HcalTopology.h:171
CaloSubdetectorGeometry::DetIdSet HcalGeometry::getCells ( const GlobalPoint r,
double  dR 
) const
overridevirtual

Get a list of all cells within a dR of the given cell.

The default implementation makes a loop over all cell geometries. Cleverer implementations are suggested to use rough conversions between eta/phi and ieta/iphi and test on the boundaries.

Reimplemented from CaloSubdetectorGeometry.

Definition at line 214 of file HcalGeometry.cc.

References reco::deltaR2(), PVValHelper::eta, PV3DBase< T, PVType, FrameType >::eta(), HcalTopology::etaMax(), etaRing(), CaloSubdetectorGeometry::getCells(), getGeometryBase(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, M_PI, m_topology, HcalTopology::maxDepth(), AlCaHLTBitMon_ParallelJobs::p, phi, PV3DBase< T, PVType, FrameType >::phi(), phiBin(), Validation_hcalonly_cfi::sign, and HcalTopology::valid().

Referenced by numberOfParametersPerShape().

215  {
216  CaloSubdetectorGeometry::DetIdSet dis; // this is the return object
217 
218  if (0.000001 < dR) {
219  if (dR > M_PI/2.) {// this version needs "small" dR
220  dis = CaloSubdetectorGeometry::getCells(r, dR); // base class version
221  } else {
222  const double dR2 ( dR*dR ) ;
223  const double reta ( r.eta() ) ;
224  const double rphi ( r.phi() ) ;
225  const double lowEta ( reta - dR ) ;
226  const double highEta ( reta + dR ) ;
227  const double lowPhi ( rphi - dR ) ;
228  const double highPhi ( rphi + dR ) ;
229 
230  const double hfEtaHi (m_topology.etaMax(HcalForward));
231 
232  if (highEta > -hfEtaHi &&
233  lowEta < hfEtaHi ) { // in hcal
235 
236  for (unsigned int is ( 0 ) ; is != 4 ; ++is ) {
237  const int sign ( reta>0 ? 1 : -1 ) ;
238  const int ieta_center ( sign*etaRing( hs[is], fabs( reta ) ) ) ;
239  const int ieta_lo ( ( 0 < lowEta*sign ? sign : -sign )*etaRing( hs[is], fabs( lowEta ) ) ) ;
240  const int ieta_hi ( ( 0 < highEta*sign ? sign : -sign )*etaRing( hs[is], fabs( highEta ) ) ) ;
241  const int iphi_lo ( phiBin( hs[is], ieta_center, lowPhi ) ) ;
242  const int iphi_hi ( phiBin( hs[is], ieta_center, highPhi ) ) ;
243  const int jphi_lo ( iphi_lo>iphi_hi ? iphi_lo - 72 : iphi_lo ) ;
244  const int jphi_hi ( iphi_hi ) ;
245 
246  const int idep_lo ( 1 == is ? 4 : 1 ) ;
247  const int idep_hi ( m_topology.maxDepth(hs[is]) );
248  for (int ieta ( ieta_lo ) ; ieta <= ieta_hi ; ++ieta) {// over eta limits
249  if (ieta != 0) {
250  for (int jphi ( jphi_lo ) ; jphi <= jphi_hi ; ++jphi) { // over phi limits
251  const int iphi ( 1 > jphi ? jphi+72 : jphi ) ;
252  for (int idep ( idep_lo ) ; idep <= idep_hi ; ++idep ) {
253  const HcalDetId did ( hs[is], ieta, iphi, idep ) ;
254  if (m_topology.valid(did)) {
255  std::shared_ptr<const CaloCellGeometry> cell ( getGeometryBase( did ) );
256  if (nullptr != cell ) {
257  const GlobalPoint& p ( cell->getPosition() ) ;
258  const double eta ( p.eta() ) ;
259  const double phi ( p.phi() ) ;
260  if (reco::deltaR2(eta, phi, reta, rphi ) < dR2) dis.insert( did ) ;
261  }
262  }
263  }
264  }
265  }
266  }
267  }
268  }
269  }
270  }
271  return dis;
272 }
bool valid(const DetId &id) const override
Geom::Phi< T > phi() const
Definition: PV3DBase.h:69
const HcalTopology & m_topology
Definition: HcalGeometry.h:159
std::shared_ptr< const CaloCellGeometry > getGeometryBase(const DetId &id) const
Definition: HcalGeometry.h:129
virtual DetIdSet getCells(const GlobalPoint &r, double dR) const
Get a list of all cells within a dR of the given cell.
int maxDepth(HcalSubdetector subdet) const
HcalSubdetector
Definition: HcalAssistant.h:31
int phiBin(HcalSubdetector bc, int etaring, double phi) const
int etaRing(HcalSubdetector bc, double abseta) const
helper methods for getClosestCell
#define M_PI
T1 deltaR2(T1 eta1, T2 phi1, T3 eta2, T4 phi2)
Definition: deltaR.h:36
T eta() const
Definition: PV3DBase.h:76
double etaMax(HcalSubdetector subdet) const
DetId HcalGeometry::getClosestCell ( const GlobalPoint r) const
overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 106 of file HcalGeometry.cc.

References funct::abs(), correctId(), HcalTopology::dddConstants(), PV3DBase< T, PVType, FrameType >::eta(), HcalTopology::etaMax(), etaRing(), getGeometry(), HcalDDDRecConstants::getMinDepth(), HcalBarrel, HcalEmpty, HcalEndcap, HcalForward, HcalTopology::incrementDepth(), m_topology, PV3DBase< T, PVType, FrameType >::mag(), PV3DBase< T, PVType, FrameType >::phi(), phiBin(), PV3DBase< T, PVType, FrameType >::z(), and ecaldqm::zside().

Referenced by ValidIsoTrkCalib::analyze(), and numberOfParametersPerShape().

106  {
107 
108  // Now find the closest eta_bin, eta value of a bin i is average
109  // of eta[i] and eta[i-1]
110  static const double z_long=1100.0;
111  double abseta = fabs(r.eta());
112  double absz = fabs(r.z());
113 
114  // figure out subdetector, giving preference to HE in HE/HF overlap region
116  if (abseta <= m_topology.etaMax(HcalBarrel) ) {
117  bc = HcalBarrel;
118  } else if (absz >= z_long) {
119  bc = HcalForward;
120  } else if (abseta <= m_topology.etaMax(HcalEndcap) ) {
121  bc = HcalEndcap;
122  } else {
123  bc = HcalForward;
124  }
125 
126  // find eta bin
127  int etaring = etaRing(bc, abseta);
128 
129  int phibin = phiBin(bc, etaring, r.phi());
130 
131  // add a sign to the etaring
132  int etabin = (r.z() > 0) ? etaring : -etaring;
133 
134  if (bc == HcalForward) {
135  static const double z_short=1137.0;
136  // Next line is premature depth 1 and 2 can coexist for large z-extent
137  // HcalDetId bestId(bc,etabin,phibin,((fabs(r.z())>=z_short)?(2):(1)));
138  // above line is no good with finite precision
139  HcalDetId bestId(bc,etabin,phibin,((fabs(r.z()) - z_short >-0.1)?(2):(1)));
140  return correctId(bestId);
141  } else {
142 
143  //Now do depth if required
144  int zside = (r.z() > 0) ? 1 : -1;
145  int dbin = m_topology.dddConstants()->getMinDepth(((int)(bc)-1),etaring,phibin,zside);
146  double pointrz(0), drz(99999.);
147  HcalDetId currentId(bc, etabin, phibin, dbin);
148  if (bc == HcalBarrel) pointrz = r.mag();
149  else pointrz = std::abs(r.z());
150  HcalDetId bestId;
151  for ( ; currentId != HcalDetId(); m_topology.incrementDepth(currentId)) {
152  std::shared_ptr<const CaloCellGeometry> cell = getGeometry(currentId);
153  if (cell == nullptr) {
154  assert (bestId != HcalDetId());
155  break;
156  } else {
157  double rz;
158  if (bc == HcalEndcap) rz = std::abs(cell->getPosition().z());
159  else rz = cell->getPosition().mag();
160  if (std::abs(pointrz-rz)<drz) {
161  bestId = currentId;
162  drz = std::abs(pointrz-rz);
163  }
164  }
165  }
166 
167  return correctId(bestId);
168  }
169 }
const HcalDDDRecConstants * dddConstants() const
Definition: HcalTopology.h:167
std::shared_ptr< const CaloCellGeometry > getGeometry(const DetId &id) const override
Get the cell geometry of a given detector id. Should return false if not found.
Definition: HcalGeometry.cc:86
Geom::Phi< T > phi() const
Definition: PV3DBase.h:69
const HcalTopology & m_topology
Definition: HcalGeometry.h:159
int zside(DetId const &)
T mag() const
Definition: PV3DBase.h:67
int getMinDepth(const int &itype, const int &ieta, const int &iphi, const int &zside) const
DetId correctId(const DetId &id) const
T z() const
Definition: PV3DBase.h:64
HcalSubdetector
Definition: HcalAssistant.h:31
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
int phiBin(HcalSubdetector bc, int etaring, double phi) const
int etaRing(HcalSubdetector bc, double abseta) const
helper methods for getClosestCell
bool incrementDepth(HcalDetId &id) const
T eta() const
Definition: PV3DBase.h:76
double etaMax(HcalSubdetector subdet) const
CaloCellGeometry::CornersVec HcalGeometry::getCorners ( const DetId id) const

Definition at line 189 of file HcalGeometry.cc.

References getGeometryBase(), electrons_cff::ids, gen::k, m_mergePosition, m_topology, and HcalTopology::unmergeDepthDetId().

Referenced by numberOfParametersPerShape().

189  {
190  if (!m_mergePosition) {
191  return (getGeometryBase(id)->getCorners());
192  } else {
193  std::vector<HcalDetId> ids;
198  for (unsigned int k=0; k<4; ++k) {
199  mcorners[k] = mcf[k];
200  mcorners[k+4] = mcb[k+4];
201  }
202  return mcorners;
203  }
204 }
const HcalTopology & m_topology
Definition: HcalGeometry.h:159
bool m_mergePosition
Definition: HcalGeometry.h:160
std::shared_ptr< const CaloCellGeometry > getGeometryBase(const DetId &id) const
Definition: HcalGeometry.h:129
void unmergeDepthDetId(const HcalDetId &id, std::vector< HcalDetId > &ids) const
Definition: HcalTopology.h:171
int k[5][pyjets_maxn]
CaloCellGeometry::CornersVec getCorners(const DetId &id) const
std::shared_ptr< const CaloCellGeometry > HcalGeometry::getGeometry ( const DetId id) const
overridevirtual

Get the cell geometry of a given detector id. Should return false if not found.

Reimplemented from CaloSubdetectorGeometry.

Definition at line 86 of file HcalGeometry.cc.

References gather_cfg::cout, HcalTopology::detId2denseId(), getGeometryBase(), HcalTopology::idFront(), m_mergePosition, and m_topology.

Referenced by getClosestCell(), and numberOfParametersPerShape().

86  {
87 #ifdef EDM_ML_DEBUG
88  std::cout << "HcalGeometry::getGeometry for " << HcalDetId(id) << " "
89  << m_mergePosition << " ";
90 #endif
91  if (!m_mergePosition) {
92 #ifdef EDM_ML_DEBUG
94  << "\n";
95 #endif
96  return getGeometryBase(id);
97  } else {
98 #ifdef EDM_ML_DEBUG
100  << getGeometryBase(m_topology.idFront(id)) << "\n";
101 #endif
102  return getGeometryBase(m_topology.idFront(id));
103  }
104 }
unsigned int detId2denseId(const DetId &id) const override
return a linear packed id
const HcalTopology & m_topology
Definition: HcalGeometry.h:159
bool m_mergePosition
Definition: HcalGeometry.h:160
std::shared_ptr< const CaloCellGeometry > getGeometryBase(const DetId &id) const
Definition: HcalGeometry.h:129
HcalDetId idFront(const HcalDetId &id) const
Definition: HcalTopology.h:176
std::shared_ptr<const CaloCellGeometry> HcalGeometry::getGeometryBase ( const DetId id) const
inlineprivate

Definition at line 129 of file HcalGeometry.h.

References CaloSubdetectorGeometry::cellGeomPtr(), correctId(), HcalTopology::detId2denseId(), etaRing(), fillDetIds(), increaseReserve(), init(), m_topology, newCellFast(), newCellImpl(), phi, phiBin(), and sortValidIds().

Referenced by getBackPosition(), getCells(), getCorners(), getGeometry(), and getPosition().

129  {
130  return cellGeomPtr( m_topology.detId2denseId( id ) ) ;
131  }
unsigned int detId2denseId(const DetId &id) const override
return a linear packed id
const HcalTopology & m_topology
Definition: HcalGeometry.h:159
virtual std::shared_ptr< const CaloCellGeometry > cellGeomPtr(uint32_t index) const
const CaloCellGeometry * HcalGeometry::getGeometryRawPtr ( uint32_t  index) const
overrideprotectedvirtual

Implements CaloSubdetectorGeometry.

Definition at line 471 of file HcalGeometry.cc.

References cropTnPTrees::din, m_hbCellVec, m_heCellVec, m_hfCellVec, m_hoCellVec, and CaloCellGeometry::param().

Referenced by sizeForDenseIndex().

471  {
472  // Modify the RawPtr class
473  const CaloCellGeometry* cell(nullptr);
474  if (m_hbCellVec.size() > din) {
475  cell = (&m_hbCellVec[din]);
476  } else if (m_hbCellVec.size()+m_heCellVec.size() > din) {
477  const unsigned int ind (din - m_hbCellVec.size() ) ;
478  cell = (&m_heCellVec[ind]);
479  } else if (m_hbCellVec.size()+m_heCellVec.size()+m_hoCellVec.size() > din) {
480  const unsigned int ind (din - m_hbCellVec.size() - m_heCellVec.size());
481  cell = (&m_hoCellVec[ind]);
482  } else if (m_hbCellVec.size()+m_heCellVec.size()+m_hoCellVec.size()+
483  m_hfCellVec.size() > din) {
484  const unsigned int ind (din - m_hbCellVec.size() - m_heCellVec.size() -
485  m_hoCellVec.size() ) ;
486  cell = (&m_hfCellVec[ind]);
487  }
488 
489  return (( nullptr == cell || nullptr == cell->param()) ? nullptr : cell ) ;
490 }
HBCellVec m_hbCellVec
Definition: HcalGeometry.h:169
HOCellVec m_hoCellVec
Definition: HcalGeometry.h:171
HFCellVec m_hfCellVec
Definition: HcalGeometry.h:172
HECellVec m_heCellVec
Definition: HcalGeometry.h:170
unsigned int HcalGeometry::getHxSize ( const int  type) const

Definition at line 274 of file HcalGeometry.cc.

References fillDetIds(), edm::AtomicPtrCache< T >::isSet(), m_emptyIds, m_hbIds, m_heIds, m_hfIds, and m_hoIds.

Referenced by HcalRecHitsAnalyzer::dqmBeginRun(), HcalDigisValidation::dqmBeginRun(), and numberOfOuterAlignments().

274  {
275  unsigned int hxsize(0);
276  if (!m_hbIds.isSet()) fillDetIds();
277  if (type == 1) hxsize = (m_hbIds.isSet()) ? m_hbIds->size() : 0;
278  else if (type == 2) hxsize = (m_heIds.isSet()) ? m_heIds->size() : 0;
279  else if (type == 3) hxsize = (m_hoIds.isSet()) ? m_hoIds->size() : 0;
280  else if (type == 4) hxsize = (m_hfIds.isSet()) ? m_hfIds->size() : 0;
281  else if (type == 0) hxsize = (m_emptyIds.isSet()) ? m_emptyIds->size() : 0;
282  return hxsize;
283 }
type
Definition: HCALResponse.h:21
edm::AtomicPtrCache< std::vector< DetId > > m_heIds
Definition: HcalGeometry.h:163
edm::AtomicPtrCache< std::vector< DetId > > m_emptyIds
Definition: HcalGeometry.h:166
bool isSet() const
edm::AtomicPtrCache< std::vector< DetId > > m_hfIds
Definition: HcalGeometry.h:165
void fillDetIds() const
Definition: HcalGeometry.cc:39
edm::AtomicPtrCache< std::vector< DetId > > m_hoIds
Definition: HcalGeometry.h:164
edm::AtomicPtrCache< std::vector< DetId > > m_hbIds
Definition: HcalGeometry.h:162
GlobalPoint HcalGeometry::getPosition ( const DetId id) const
void HcalGeometry::getSummary ( CaloSubdetectorGeometry::TrVec trVector,
CaloSubdetectorGeometry::IVec iVector,
CaloSubdetectorGeometry::DimVec dimVector,
CaloSubdetectorGeometry::IVec dinsVector 
) const
overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 492 of file HcalGeometry.cc.

References begin, CaloSubdetectorGeometry::cellGeomPtr(), runTauDisplay::gp, mps_fire::i, gen::k, m_dins, m_topology, HcalTopology::ncells(), groupFilesInBlocks::nn, numberOfParametersPerShape(), numberOfShapes(), CaloSubdetectorGeometry::numberOfTransformParms(), convertSQLiteXML::ok, CaloSubdetectorGeometry::parVecVec(), MetAnalyzer::pv(), findQualityFiles::rr, and groupFilesInBlocks::tt.

Referenced by numberOfAlignments().

495  {
496 
497  tVec.reserve( m_topology.ncells()*numberOfTransformParms());
498  iVec.reserve( numberOfShapes()==1 ? 1 : m_topology.ncells());
499  dVec.reserve( numberOfShapes()*numberOfParametersPerShape());
500  dinsVec.reserve( m_topology.ncells());
501 
502  for (const auto & pv : parVecVec()) {
503  for (float iv : pv) {
504  dVec.emplace_back( iv ) ;
505  }
506  }
507 
508  for( auto i : m_dins ) {
509  Tr3D tr ;
510  auto ptr = cellGeomPtr( i );
511 
512  if (nullptr != ptr) {
513  dinsVec.emplace_back( i );
514 
515  const CCGFloat* par ( ptr->param() ) ;
516 
517  unsigned int ishape ( 9999 ) ;
518 
519  for( unsigned int ivv ( 0 ) ; ivv != parVecVec().size() ; ++ivv ) {
520  bool ok ( true ) ;
521  const CCGFloat* pv ( &(*parVecVec()[ivv].begin() ) ) ;
522  for( unsigned int k ( 0 ) ; k != numberOfParametersPerShape() ; ++k ) {
523  ok = ok && ( fabs( par[k] - pv[k] ) < 1.e-6 ) ;
524  }
525  if( ok ) {
526  ishape = ivv;
527  break ;
528  }
529  }
530  assert( 9999 != ishape ) ;
531 
532  const unsigned int nn (( numberOfShapes()==1) ? (unsigned int)1 : m_dins.size() ) ;
533  if( iVec.size() < nn ) iVec.emplace_back( ishape ) ;
534 
535  ptr->getTransform( tr, ( Pt3DVec* ) nullptr ) ;
536 
537  if( Tr3D() == tr ) { // for preshower there is no rotation
538  const GlobalPoint& gp ( ptr->getPosition() ) ;
539  tr = HepGeom::Translate3D( gp.x(), gp.y(), gp.z() ) ;
540  }
541 
542  const CLHEP::Hep3Vector tt ( tr.getTranslation() ) ;
543  tVec.emplace_back( tt.x() ) ;
544  tVec.emplace_back( tt.y() ) ;
545  tVec.emplace_back( tt.z() ) ;
546  if (6 == numberOfTransformParms()) {
547  const CLHEP::HepRotation rr ( tr.getRotation() ) ;
548  const ROOT::Math::Transform3D rtr (rr.xx(), rr.xy(), rr.xz(), tt.x(),
549  rr.yx(), rr.yy(), rr.yz(), tt.y(),
550  rr.zx(), rr.zy(), rr.zz(), tt.z());
552  rtr.GetRotation( ea ) ;
553  tVec.emplace_back( ea.Phi() ) ;
554  tVec.emplace_back( ea.Theta() ) ;
555  tVec.emplace_back( ea.Psi() ) ;
556  }
557  }
558  }
559 }
unsigned int numberOfShapes() const override
Definition: HcalGeometry.h:45
const HcalTopology & m_topology
Definition: HcalGeometry.h:159
CaloCellGeometry::Tr3D Tr3D
Definition: HcalGeometry.cc:12
virtual std::shared_ptr< const CaloCellGeometry > cellGeomPtr(uint32_t index) const
unsigned int numberOfParametersPerShape() const override
Definition: HcalGeometry.h:46
def pv(vc)
Definition: MetAnalyzer.py:6
CaloCellGeometry::CCGFloat CCGFloat
CaloSubdetectorGeometry::IVec m_dins
Definition: HcalGeometry.h:167
int k[5][pyjets_maxn]
AlgebraicVector EulerAngles
Definition: Definitions.h:36
CaloCellGeometry::Pt3DVec Pt3DVec
CaloCellGeometry::Tr3D Tr3D
#define begin
Definition: vmac.h:32
unsigned int ncells() const override
return a count of valid cells (for dense indexing use)
virtual unsigned int numberOfTransformParms() const
const std::vector< DetId > & HcalGeometry::getValidDetIds ( DetId::Detector  det = DetId::Detector ( 0 ),
int  subdet = 0 
) const
overridevirtual

Get a list of valid detector ids (for the given subdetector)

Note
The implementation in this class is relevant for SubdetectorGeometries which handle only a single subdetector at a time. It does not look at the det and subdet arguments.

Reimplemented from CaloSubdetectorGeometry.

Definition at line 76 of file HcalGeometry.cc.

References fillDetIds(), CaloSubdetectorGeometry::getValidDetIds(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, edm::AtomicPtrCache< T >::load(), m_emptyIds, m_hbIds, m_heIds, m_hfIds, and m_hoIds.

Referenced by HcalHBHEMuonAnalyzer::beginRun(), cms::MinBias::endJob(), HCALRecHitAnalyzer::FillGeometry(), and numberOfParametersPerShape().

77  {
78  if( 0 != subdet ) fillDetIds() ;
79  return ( 0 == subdet ? CaloSubdetectorGeometry::getValidDetIds() :
80  ( HcalBarrel == subdet ? *m_hbIds.load() :
81  ( HcalEndcap == subdet ? *m_heIds.load() :
82  ( HcalOuter == subdet ? *m_hoIds.load() :
83  ( HcalForward == subdet ? *m_hfIds.load() : *m_emptyIds.load() ) ) ) ) ) ;
84 }
edm::AtomicPtrCache< std::vector< DetId > > m_heIds
Definition: HcalGeometry.h:163
edm::AtomicPtrCache< std::vector< DetId > > m_emptyIds
Definition: HcalGeometry.h:166
virtual const std::vector< DetId > & getValidDetIds(DetId::Detector det=DetId::Detector(0), int subdet=0) const
Get a list of valid detector ids (for the given subdetector)
edm::AtomicPtrCache< std::vector< DetId > > m_hfIds
Definition: HcalGeometry.h:165
void fillDetIds() const
Definition: HcalGeometry.cc:39
edm::AtomicPtrCache< std::vector< DetId > > m_hoIds
Definition: HcalGeometry.h:164
T const * load() const
edm::AtomicPtrCache< std::vector< DetId > > m_hbIds
Definition: HcalGeometry.h:162
void HcalGeometry::increaseReserve ( unsigned int  extra)
private
unsigned int HcalGeometry::indexFor ( const DetId id) const
inlineoverrideprotectedvirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 120 of file HcalGeometry.h.

References HcalTopology::detId2denseId(), and m_topology.

120 { return m_topology.detId2denseId(id); }
unsigned int detId2denseId(const DetId &id) const override
return a linear packed id
const HcalTopology & m_topology
Definition: HcalGeometry.h:159
void HcalGeometry::init ( void  )
private

Definition at line 23 of file HcalGeometry.cc.

References gather_cfg::cout, HcalTopology::getHBSize(), HcalTopology::getHESize(), HcalTopology::getHFSize(), HcalTopology::getHOSize(), HcalTopology::getMergePositionFlag(), m_hbCellVec, m_heCellVec, m_hfCellVec, m_hoCellVec, m_mergePosition, and m_topology.

Referenced by getGeometryBase(), and HcalGeometry().

23  {
25 #ifdef EDM_ML_DEBUG
26  std::cout << "HcalGeometry: " << "HcalGeometry::init() "
27  << " HBSize " << m_topology.getHBSize()
28  << " HESize " << m_topology.getHESize()
29  << " HOSize " << m_topology.getHOSize()
30  << " HFSize " << m_topology.getHFSize() << std::endl;
31 #endif
32 
37 }
unsigned int getHFSize() const
Definition: HcalTopology.h:141
unsigned int getHOSize() const
Definition: HcalTopology.h:140
const HcalTopology & m_topology
Definition: HcalGeometry.h:159
std::vector< IdealObliquePrism > HECellVec
Definition: HcalGeometry.h:27
bool m_mergePosition
Definition: HcalGeometry.h:160
bool getMergePositionFlag() const
Definition: HcalTopology.h:170
HBCellVec m_hbCellVec
Definition: HcalGeometry.h:169
HOCellVec m_hoCellVec
Definition: HcalGeometry.h:171
HFCellVec m_hfCellVec
Definition: HcalGeometry.h:172
std::vector< IdealObliquePrism > HOCellVec
Definition: HcalGeometry.h:28
std::vector< IdealZPrism > HFCellVec
Definition: HcalGeometry.h:29
std::vector< IdealObliquePrism > HBCellVec
Definition: HcalGeometry.h:26
unsigned int getHESize() const
Definition: HcalTopology.h:139
HECellVec m_heCellVec
Definition: HcalGeometry.h:170
unsigned int getHBSize() const
Definition: HcalTopology.h:138
void HcalGeometry::localCorners ( Pt3DVec lc,
const CCGFloat pv,
unsigned int  i,
Pt3D ref 
)

Definition at line 395 of file HcalGeometry.cc.

References HcalTopology::denseId2detId(), HcalForward, IdealObliquePrism::localCorners(), IdealZPrism::localCorners(), m_topology, and HcalDetId::subdet().

Referenced by numberOfAlignments().

398  {
400 
401  if (hid.subdet() == HcalForward ) {
402  IdealZPrism::localCorners( lc, pv, ref ) ;
403  } else {
404  IdealObliquePrism::localCorners( lc, pv, ref ) ;
405  }
406 }
DetId denseId2detId(unsigned int) const override
return a linear packed id
HcalSubdetector subdet() const
get the subdetector
Definition: HcalDetId.h:49
const HcalTopology & m_topology
Definition: HcalGeometry.h:159
static void localCorners(Pt3DVec &vec, const CCGFloat *pv, Pt3D &ref)
def pv(vc)
Definition: MetAnalyzer.py:6
static void localCorners(Pt3DVec &vec, const CCGFloat *pv, Pt3D &ref)
Definition: IdealZPrism.cc:122
void HcalGeometry::newCell ( const GlobalPoint f1,
const GlobalPoint f2,
const GlobalPoint f3,
const CCGFloat parm,
const DetId detId 
)
overridevirtual

Implements CaloSubdetectorGeometry.

Definition at line 447 of file HcalGeometry.cc.

References CaloSubdetectorGeometry::addValidID(), cropTnPTrees::din, m_dins, and newCellImpl().

Referenced by numberOfAlignments().

451  {
452 
453  unsigned int din = newCellImpl(f1,f2,f3,parm,detId);
454 
455  addValidID( detId ) ;
456  m_dins.emplace_back( din );
457 }
CaloSubdetectorGeometry::IVec m_dins
Definition: HcalGeometry.h:167
void addValidID(const DetId &id)
unsigned int newCellImpl(const GlobalPoint &f1, const GlobalPoint &f2, const GlobalPoint &f3, const CCGFloat *parm, const DetId &detId)
void HcalGeometry::newCellFast ( const GlobalPoint f1,
const GlobalPoint f2,
const GlobalPoint f3,
const CCGFloat parm,
const DetId detId 
)
private

Definition at line 459 of file HcalGeometry.cc.

References cropTnPTrees::din, m_dins, CaloSubdetectorGeometry::m_validIds, and newCellImpl().

Referenced by HcalHardcodeGeometryLoader::fillHBHO(), HcalFlexiHardcodeGeometryLoader::fillHBHO(), HcalHardcodeGeometryLoader::fillHE(), HcalFlexiHardcodeGeometryLoader::fillHE(), HcalHardcodeGeometryLoader::fillHF(), HcalFlexiHardcodeGeometryLoader::fillHF(), and getGeometryBase().

463  {
464 
465  unsigned int din = newCellImpl(f1,f2,f3,parm,detId);
466 
467  m_validIds.emplace_back( detId ) ;
468  m_dins.emplace_back( din );
469 }
std::vector< DetId > m_validIds
CaloSubdetectorGeometry::IVec m_dins
Definition: HcalGeometry.h:167
unsigned int newCellImpl(const GlobalPoint &f1, const GlobalPoint &f2, const GlobalPoint &f3, const CCGFloat *parm, const DetId &detId)
unsigned int HcalGeometry::newCellImpl ( const GlobalPoint f1,
const GlobalPoint f2,
const GlobalPoint f3,
const CCGFloat parm,
const DetId detId 
)
private

Definition at line 408 of file HcalGeometry.cc.

References CaloSubdetectorGeometry::cornersMgr(), gather_cfg::cout, DetId::det(), HcalTopology::detId2denseId(), cropTnPTrees::din, IdealZPrism::EM, IdealZPrism::HADR, DetId::Hcal, HcalBarrel, HcalEndcap, HcalOuter, m_hbCellVec, m_heCellVec, m_hfCellVec, m_hoCellVec, m_topology, DetId::rawId(), and DetId::subdetId().

Referenced by getGeometryBase(), newCell(), and newCellFast().

412  {
413 
414  assert (detId.det()==DetId::Hcal);
415 
416  const HcalDetId hid ( detId ) ;
417  unsigned int din=m_topology.detId2denseId(detId);
418 
419 #ifdef EDM_ML_DEBUG
420  std::cout << "HcalGeometry: " << " newCell subdet "
421  << detId.subdetId() << ", raw ID "
422  << detId.rawId() << ", hid " << hid << ", din "
423  << din << ", index " << std::endl;
424 #endif
425 
426  if (hid.subdet()==HcalBarrel) {
427  m_hbCellVec.at( din ) = IdealObliquePrism( f1, cornersMgr(), parm ) ;
428  } else if (hid.subdet()==HcalEndcap) {
429  const unsigned int index ( din - m_hbCellVec.size() ) ;
430  m_heCellVec.at( index ) = IdealObliquePrism( f1, cornersMgr(), parm ) ;
431  } else if (hid.subdet()==HcalOuter) {
432  const unsigned int index ( din
433  - m_hbCellVec.size()
434  - m_heCellVec.size() ) ;
435  m_hoCellVec.at( index ) = IdealObliquePrism( f1, cornersMgr(), parm ) ;
436  } else {
437  const unsigned int index ( din
438  - m_hbCellVec.size()
439  - m_heCellVec.size()
440  - m_hoCellVec.size() ) ;
441  m_hfCellVec.at( index ) = IdealZPrism( f1, cornersMgr(), parm, hid.depth()==1 ? IdealZPrism::EM : IdealZPrism::HADR ) ;
442  }
443 
444  return din;
445 }
unsigned int detId2denseId(const DetId &id) const override
return a linear packed id
const HcalTopology & m_topology
Definition: HcalGeometry.h:159
uint32_t rawId() const
get the raw id
Definition: DetId.h:44
HBCellVec m_hbCellVec
Definition: HcalGeometry.h:169
HOCellVec m_hoCellVec
Definition: HcalGeometry.h:171
HFCellVec m_hfCellVec
Definition: HcalGeometry.h:172
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:38
CaloCellGeometry::CornersMgr * cornersMgr()
HECellVec m_heCellVec
Definition: HcalGeometry.h:170
Detector det() const
get the detector field from this detid
Definition: DetId.h:36
static unsigned int HcalGeometry::numberOfAlignments ( )
inlinestatic
static unsigned int HcalGeometry::numberOfBarrelAlignments ( )
inlinestatic
static unsigned int HcalGeometry::numberOfEndcapAlignments ( )
inlinestatic
static unsigned int HcalGeometry::numberOfForwardAlignments ( )
inlinestatic
static unsigned int HcalGeometry::numberOfOuterAlignments ( )
inlinestatic
unsigned int HcalGeometry::numberOfParametersPerShape ( ) const
inlineoverridevirtual
unsigned int HcalGeometry::numberOfShapes ( ) const
inlineoverridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 45 of file HcalGeometry.h.

References HcalTopology::getNumberOfShapes(), and m_topology.

Referenced by getSummary(), HcalHardcodeGeometryLoader::load(), and HcalFlexiHardcodeGeometryLoader::load().

45 { return m_topology.getNumberOfShapes() ; }
const HcalTopology & m_topology
Definition: HcalGeometry.h:159
unsigned int getNumberOfShapes() const
Definition: HcalTopology.h:164
int HcalGeometry::phiBin ( HcalSubdetector  bc,
int  etaring,
double  phi 
) const
private

Definition at line 210 of file HcalGeometry.cc.

References m_topology, and HcalTopology::phiBin().

Referenced by getCells(), getClosestCell(), and getGeometryBase().

210  {
211  return m_topology.phiBin(bc, etaring, phi);
212 }
int phiBin(HcalSubdetector subdet, int etaRing, double phi) const
const HcalTopology & m_topology
Definition: HcalGeometry.h:159
static std::string HcalGeometry::producerTag ( )
inlinestatic
unsigned int HcalGeometry::sizeForDenseIndex ( const DetId id) const
inlineoverrideprotectedvirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 121 of file HcalGeometry.h.

References getGeometryRawPtr(), m_topology, and HcalTopology::ncells().

121 { return m_topology.ncells(); }
const HcalTopology & m_topology
Definition: HcalGeometry.h:159
unsigned int ncells() const override
return a count of valid cells (for dense indexing use)
void HcalGeometry::sortValidIds ( )
private

Definition at line 575 of file HcalGeometry.cc.

References CaloSubdetectorGeometry::m_validIds.

Referenced by getGeometryBase(), HcalHardcodeGeometryLoader::load(), and HcalFlexiHardcodeGeometryLoader::load().

575  {
576  std::sort(m_validIds.begin(),m_validIds.end());
577 }
std::vector< DetId > m_validIds
const HcalTopology& HcalGeometry::topology ( ) const
inline

Definition at line 116 of file HcalGeometry.h.

References m_topology.

Referenced by CaloTruthAccumulator::beginLuminosityBlock(), HGCalTriggerGeometryBase::bhTopology(), and numberOfParametersPerShape().

116 { return m_topology; }
const HcalTopology & m_topology
Definition: HcalGeometry.h:159

Friends And Related Function Documentation

friend class HcalFlexiHardcodeGeometryLoader
friend

Definition at line 23 of file HcalGeometry.h.

friend class HcalGeometryPlan1Tester
friend

Definition at line 19 of file HcalGeometry.h.

friend class HcalHardcodeGeometryLoader
friend

Definition at line 24 of file HcalGeometry.h.

Member Data Documentation

CaloSubdetectorGeometry::IVec HcalGeometry::m_dins
private

Definition at line 167 of file HcalGeometry.h.

Referenced by getSummary(), newCell(), and newCellFast().

edm::AtomicPtrCache<std::vector<DetId> > HcalGeometry::m_emptyIds
mutableprivate

Definition at line 166 of file HcalGeometry.h.

Referenced by fillDetIds(), getHxSize(), and getValidDetIds().

HBCellVec HcalGeometry::m_hbCellVec
private

Definition at line 169 of file HcalGeometry.h.

Referenced by getGeometryRawPtr(), init(), and newCellImpl().

edm::AtomicPtrCache<std::vector<DetId> > HcalGeometry::m_hbIds
mutableprivate

Definition at line 162 of file HcalGeometry.h.

Referenced by fillDetIds(), getHxSize(), and getValidDetIds().

HECellVec HcalGeometry::m_heCellVec
private

Definition at line 170 of file HcalGeometry.h.

Referenced by getGeometryRawPtr(), init(), and newCellImpl().

edm::AtomicPtrCache<std::vector<DetId> > HcalGeometry::m_heIds
mutableprivate

Definition at line 163 of file HcalGeometry.h.

Referenced by fillDetIds(), getHxSize(), and getValidDetIds().

HFCellVec HcalGeometry::m_hfCellVec
private

Definition at line 172 of file HcalGeometry.h.

Referenced by getGeometryRawPtr(), init(), and newCellImpl().

edm::AtomicPtrCache<std::vector<DetId> > HcalGeometry::m_hfIds
mutableprivate

Definition at line 165 of file HcalGeometry.h.

Referenced by fillDetIds(), getHxSize(), and getValidDetIds().

HOCellVec HcalGeometry::m_hoCellVec
private

Definition at line 171 of file HcalGeometry.h.

Referenced by getGeometryRawPtr(), init(), and newCellImpl().

edm::AtomicPtrCache<std::vector<DetId> > HcalGeometry::m_hoIds
mutableprivate

Definition at line 164 of file HcalGeometry.h.

Referenced by fillDetIds(), getHxSize(), and getValidDetIds().

bool HcalGeometry::m_mergePosition
private

Definition at line 160 of file HcalGeometry.h.

Referenced by correctId(), getBackPosition(), getCorners(), getGeometry(), getPosition(), and init().

const HcalTopology& HcalGeometry::m_topology
private