CMS 3D CMS Logo

List of all members | Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Private Attributes
EcalPreshowerGeometry Class Referencefinal

#include <EcalPreshowerGeometry.h>

Inheritance diagram for EcalPreshowerGeometry:
CaloSubdetectorGeometry

Public Types

enum  { k_NumberOfCellsForCorners = ESDetId::kSizeForDenseIndexing }
 
enum  { k_NumberOfShapes = 4 }
 
enum  { k_NumberOfParametersPerShape = 4 }
 
typedef EcalPreshowerGeometryRecord AlignedRecord
 
typedef ESAlignmentRcd AlignmentRecord
 
typedef CaloCellGeometry::CCGFloat CCGFloat
 
typedef std::vector< PreshowerStripCellVec
 
typedef ESDetId DetIdType
 
typedef IdealGeometryRecord IdealRecord
 
typedef EcalPreshowerNumberingScheme NumberingScheme
 
typedef CaloSubdetectorGeometry::ParVec ParVec
 
typedef CaloSubdetectorGeometry::ParVecVec ParVecVec
 
typedef PEcalPreshowerRcd 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

 EcalPreshowerGeometry ()
 
DetId getClosestCell (const GlobalPoint &r) const override
 
virtual DetId getClosestCellInPlane (const GlobalPoint &r, int plane) const
 
void initializeParms () override
 
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
 
unsigned int numberOfTransformParms () const override
 
bool present (const DetId &id) const override
 is this detid present in the geometry? More...
 
void setzPlanes (CCGFloat z1minus, CCGFloat z2minus, CCGFloat z1plus, CCGFloat z2plus)
 
 ~EcalPreshowerGeometry () override
 The EcalPreshowerGeometry 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 DetIdSet getCells (const GlobalPoint &r, double dR) const
 Get a list of all cells within a dR of the given cell. More...
 
virtual CellSet getCellSet (const GlobalPoint &r, double dR) const
 
virtual std::shared_ptr< const CaloCellGeometrygetGeometry (const DetId &id) const
 Get the cell geometry of a given detector id. Should return false if not found. More...
 
virtual void getSummary (TrVec &trVector, IVec &iVector, DimVec &dimVector, IVec &dinsVector) 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) More...
 
CaloSubdetectorGeometryoperator= (const CaloSubdetectorGeometry &)=delete
 
ParMgrparMgr ()
 
const ParMgrparMgrConst () const
 
ParVecVecparVecVec ()
 
const ParVecVecparVecVec () const
 
virtual ~CaloSubdetectorGeometry ()
 The base class DOES assume that it owns the CaloCellGeometry objects. More...
 

Static Public Member Functions

static unsigned int alignmentTransformIndexGlobal (const DetId &id)
 
static unsigned int alignmentTransformIndexLocal (const DetId &id)
 
static std::string dbString ()
 
static DetId detIdFromLocalAlignmentIndex (unsigned int iLoc)
 
static std::string hitString ()
 
static void localCorners (Pt3DVec &lc, const CCGFloat *pv, unsigned int i, Pt3D &ref)
 
static unsigned int numberOfAlignments ()
 
static std::string producerTag ()
 

Protected Member Functions

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

Private Attributes

CellVec m_cellVec
 
const CCGFloat m_xInterLadGap
 
const CCGFloat m_xIntraLadGap
 
const CCGFloat m_xWidWaf
 
const CCGFloat m_yCtrOff
 
const CCGFloat m_yWidAct
 
CCGFloat m_zplane [4]
 

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 15 of file EcalPreshowerGeometry.h.

Member Typedef Documentation

Definition at line 24 of file EcalPreshowerGeometry.h.

Definition at line 25 of file EcalPreshowerGeometry.h.

Definition at line 19 of file EcalPreshowerGeometry.h.

Definition at line 17 of file EcalPreshowerGeometry.h.

Definition at line 31 of file EcalPreshowerGeometry.h.

Definition at line 23 of file EcalPreshowerGeometry.h.

Definition at line 28 of file EcalPreshowerGeometry.h.

Definition at line 29 of file EcalPreshowerGeometry.h.

Definition at line 30 of file EcalPreshowerGeometry.h.

Definition at line 26 of file EcalPreshowerGeometry.h.

Definition at line 20 of file EcalPreshowerGeometry.h.

Definition at line 21 of file EcalPreshowerGeometry.h.

Member Enumeration Documentation

anonymous enum
anonymous enum
Enumerator
k_NumberOfShapes 

Definition at line 35 of file EcalPreshowerGeometry.h.

anonymous enum
Enumerator
k_NumberOfParametersPerShape 

Definition at line 37 of file EcalPreshowerGeometry.h.

Constructor & Destructor Documentation

EcalPreshowerGeometry::EcalPreshowerGeometry ( )

Definition at line 13 of file EcalPreshowerGeometry.cc.

References m_zplane.

Referenced by numberOfParametersPerShape().

14  : m_xWidWaf(6.3),
15  m_xInterLadGap(0.05), // additional gap between wafers in adj ladders
16  m_xIntraLadGap(0.04), // gap between wafers in same ladder
17  m_yWidAct(6.1),
18  m_yCtrOff(0.05), // gap at center
20  m_zplane[0] = 0.;
21  m_zplane[1] = 0.;
22  m_zplane[2] = 0.;
23  m_zplane[3] = 0.;
24 }
EcalPreshowerGeometry::~EcalPreshowerGeometry ( )
override

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

Definition at line 26 of file EcalPreshowerGeometry.cc.

Referenced by numberOfParametersPerShape().

26 {}

Member Function Documentation

unsigned int EcalPreshowerGeometry::alignmentTransformIndexGlobal ( const DetId id)
static

Definition at line 69 of file EcalPreshowerGeometry.cc.

References DetId::Ecal.

Referenced by numberOfAlignments().

69  {
70  return (unsigned int)DetId::Ecal - 1;
71 }
unsigned int EcalPreshowerGeometry::alignmentTransformIndexLocal ( const DetId id)
static

Definition at line 28 of file EcalPreshowerGeometry.cc.

References CaloGenericDetId::isES(), and edm::second().

Referenced by numberOfAlignments().

28  {
29  const CaloGenericDetId gid(id);
30 
31  assert(gid.isES());
32 
33  // plane 2 is split into 2 dees along x=0 for both signs of z
34 
35  // plane 1 at zsign=-1 is split into 2 dees between six=19 and six=20 for siy<=20,
36  // and six=21 and 22 for siy>=21
37 
38  // plane 1 at zsign=+1 is split into 2 dees between six=20 and six=21 for siy<=20,
39  // and six=19 and 20 for siy>=21
40 
41  // Desired numbering
42  // LEFT RIGHT (as one faces the Dee from the IP)
43  // ES- pl=2 0 1
44  // pl=1 2 3 the reversal of pl=2 and pl=1 is intentional here (CM Kuo)
45  // ES+ pl=1 4 5
46  // pl=2 6 7
47 
48  const ESDetId esid(id);
49  const int jx(esid.six() - 1);
50  const int jy(esid.siy() - 1);
51  const int jz(esid.zside() + 1);
52  const int pl(esid.plane() - 1);
53  const bool second(1 == pl);
54  const bool top(19 < jy);
55  const bool negz(0 == jz);
56  const int lrl(19 > jx ? 0 : 1);
57  const int lrr(21 > jx ? 0 : 1);
58 
59  return (second ? jx / 20 + 3 * jz : // 2nd plane split along middle
60  (negz && !top ? lrl + 2 : // 1st plane at neg z and bottom half split at six=19&20
61  (negz && top ? lrr + 2 : // 1st plane at neg z and top half split at six=21&22
62  (!negz && !top ? lrr + 4 : lrl + 4)))); // opposite at positive z
63 }
U second(std::pair< T, U > const &p)
static std::string EcalPreshowerGeometry::dbString ( )
inlinestatic

Definition at line 39 of file EcalPreshowerGeometry.h.

39 { return "PEcalPreshowerRcd"; }
DetId EcalPreshowerGeometry::detIdFromLocalAlignmentIndex ( unsigned int  iLoc)
static

Definition at line 65 of file EcalPreshowerGeometry.cc.

Referenced by WriteESAlignments::convert(), numberOfAlignments(), FakeCaloAlignmentEP::produceESAli(), and TestCaloAlignmentEP::produceESAli().

65  {
66  return ESDetId(1, 10 + 20 * (iLoc % 2), 10, 2 > iLoc || 5 < iLoc ? 2 : 1, 2 * (iLoc / 4) - 1);
67 }
DetId EcalPreshowerGeometry::getClosestCell ( const GlobalPoint r) const
overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 128 of file EcalPreshowerGeometry.cc.

References getClosestCellInPlane().

Referenced by numberOfParametersPerShape().

128 { return getClosestCellInPlane(point, 2); }
virtual DetId getClosestCellInPlane(const GlobalPoint &r, int plane) const
*vegas h *****************************************************used in the default bin number in original ***version of VEGAS is ***a higher bin number might help to derive a more precise ***grade subtle point
Definition: invegas.h:5
DetId EcalPreshowerGeometry::getClosestCellInPlane ( const GlobalPoint r,
int  plane 
) const
virtual

Definition at line 130 of file EcalPreshowerGeometry.cc.

References cuy::col, CaloSubdetectorGeometry::getGeometry(), ESDetId::ISTRIP_MAX, ESDetId::ISTRIP_MIN, m_xInterLadGap, m_xIntraLadGap, m_xWidWaf, m_yCtrOff, m_yWidAct, m_zplane, AlCaHLTBitMon_ParallelJobs::p, present(), ESDetId::validDetId(), x, PV3DBase< T, PVType, FrameType >::x(), y, PV3DBase< T, PVType, FrameType >::y(), z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by PreshowerHitMaker::addHit(), getClosestCell(), numberOfParametersPerShape(), and ReducedESRecHitCollectionProducer::produce().

130  {
131  const CCGFloat x(point.x());
132  const CCGFloat y(point.y());
133  const CCGFloat z(point.z());
134 
135  if (0 == z || 1 > plane || 2 < plane)
136  return DetId(0);
137 
138  const unsigned int iz((0 > z ? 0 : 2) + plane - 1);
139 
140  const CCGFloat ze(m_zplane[iz]);
141  const CCGFloat xe(x * ze / z);
142  const CCGFloat ye(y * ze / z);
143 
144  const CCGFloat x0(1 == plane ? xe : ye);
145  const CCGFloat y0(1 == plane ? ye : xe);
146 
147  static const CCGFloat xWid(m_xWidWaf + m_xIntraLadGap + m_xInterLadGap);
148 
149  const int row(1 + int(y0 + 20. * m_yWidAct - m_yCtrOff) / m_yWidAct);
150  const int col(1 + int((x0 + 20. * xWid) / xWid));
151 
152  CCGFloat closest(1e9);
153 
154  DetId detId(0);
155 
156  const int jz(0 > ze ? -1 : 1);
157 
158  // std::cout<<"** p="<<point<<", ("<<xe<<", "<<ye<<", "<<ze<<"), row="<<row<<", col="<<col<<std::endl;
159 
160  for (int ix(-1); ix != 2; ++ix) // search within +-1 in row and col
161  {
162  for (int iy(-1); iy != 2; ++iy) {
163  for (int jstrip(ESDetId::ISTRIP_MIN); jstrip <= ESDetId::ISTRIP_MAX; ++jstrip) {
164  const int jx(1 == plane ? col + ix : row + iy);
165  const int jy(1 == plane ? row + iy : col + ix);
166  if (ESDetId::validDetId(jstrip, jx, jy, plane, jz)) {
167  const ESDetId esId(jstrip, jx, jy, plane, jz);
168  auto cell = getGeometry(esId);
169  if (nullptr != cell) {
170  const GlobalPoint& p(cell->getPosition());
171  const CCGFloat dist2((p.x() - xe) * (p.x() - xe) + (p.y() - ye) * (p.y() - ye));
172  if (dist2 < closest && present(esId)) {
173  closest = dist2;
174  detId = esId;
175  }
176  }
177  }
178  }
179  }
180  }
181  return detId;
182 }
static const int ISTRIP_MAX
Definition: ESDetId.h:69
CaloCellGeometry::CCGFloat CCGFloat
Definition: DetId.h:17
virtual std::shared_ptr< const CaloCellGeometry > getGeometry(const DetId &id) const
Get the cell geometry of a given detector id. Should return false if not found.
static bool validDetId(int istrip, int ixs, int iys, int iplane, int iz)
check if a valid index combination
Definition: ESDetId.cc:15
col
Definition: cuy.py:1010
bool present(const DetId &id) const override
is this detid present in the geometry?
static const int ISTRIP_MIN
Definition: ESDetId.h:68
*vegas h *****************************************************used in the default bin number in original ***version of VEGAS is ***a higher bin number might help to derive a more precise ***grade subtle point
Definition: invegas.h:5
const CaloCellGeometry * EcalPreshowerGeometry::getGeometryRawPtr ( uint32_t  index) const
overrideprotectedvirtual

Implements CaloSubdetectorGeometry.

Definition at line 198 of file EcalPreshowerGeometry.cc.

References m_cellVec, and CaloCellGeometry::param().

Referenced by present().

198  {
199  // Modify the RawPtr class
200  const CaloCellGeometry* cell(&m_cellVec[index]);
201  return (m_cellVec.size() < index || nullptr == cell->param() ? nullptr : cell);
202 }
static std::string EcalPreshowerGeometry::hitString ( )
inlinestatic

Definition at line 60 of file EcalPreshowerGeometry.h.

60 { return "EcalHitsES"; }
void EcalPreshowerGeometry::initializeParms ( )
overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 73 of file EcalPreshowerGeometry.cc.

References CaloSubdetectorGeometry::getGeometry(), CaloSubdetectorGeometry::getValidDetIds(), mps_fire::i, ESDetId::plane(), setzPlanes(), ESDetId::zside(), and geometryCSVtoXML::zz.

Referenced by numberOfParametersPerShape().

73  {
74  unsigned int n1minus(0);
75  unsigned int n2minus(0);
76  unsigned int n1plus(0);
77  unsigned int n2plus(0);
78  CCGFloat z1minus(0);
79  CCGFloat z2minus(0);
80  CCGFloat z1plus(0);
81  CCGFloat z2plus(0);
82  const std::vector<DetId>& esDetIds(getValidDetIds());
83 
84  for (unsigned int i(0); i != esDetIds.size(); ++i) {
85  const ESDetId esid(esDetIds[i]);
86  auto cell = getGeometry(esid);
87  if (nullptr != cell) {
88  const CCGFloat zz(cell->getPosition().z());
89  if (1 == esid.plane()) {
90  if (0 > esid.zside()) {
91  z1minus += zz;
92  ++n1minus;
93  } else {
94  z1plus += zz;
95  ++n1plus;
96  }
97  }
98  if (2 == esid.plane()) {
99  if (0 > esid.zside()) {
100  z2minus += zz;
101  ++n2minus;
102  } else {
103  z2plus += zz;
104  ++n2plus;
105  }
106  }
107  }
108  }
109  assert(0 != n1minus && 0 != n2minus && 0 != n1plus && 0 != n2plus);
110  z1minus /= (1. * n1minus);
111  z2minus /= (1. * n2minus);
112  z1plus /= (1. * n1plus);
113  z2plus /= (1. * n2plus);
114  assert(0 != z1minus && 0 != z2minus && 0 != z1plus && 0 != z2plus);
115  setzPlanes(z1minus, z2minus, z1plus, z2plus);
116 }
void setzPlanes(CCGFloat z1minus, CCGFloat z2minus, CCGFloat z1plus, CCGFloat z2plus)
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)
CaloCellGeometry::CCGFloat CCGFloat
virtual std::shared_ptr< const CaloCellGeometry > getGeometry(const DetId &id) const
Get the cell geometry of a given detector id. Should return false if not found.
void EcalPreshowerGeometry::localCorners ( Pt3DVec lc,
const CCGFloat pv,
unsigned int  i,
Pt3D ref 
)
static

Definition at line 184 of file EcalPreshowerGeometry.cc.

References PreshowerStrip::localCorners().

Referenced by numberOfAlignments().

184  {
186 }
static void localCorners(Pt3DVec &vec, const CCGFloat *pv, Pt3D &ref)
def pv(vc)
Definition: MetAnalyzer.py:7
void EcalPreshowerGeometry::newCell ( const GlobalPoint f1,
const GlobalPoint f2,
const GlobalPoint f3,
const CCGFloat parm,
const DetId detId 
)
overridevirtual

Implements CaloSubdetectorGeometry.

Definition at line 188 of file EcalPreshowerGeometry.cc.

References CaloSubdetectorGeometry::addValidID(), CaloSubdetectorGeometry::cornersMgr(), and m_cellVec.

Referenced by numberOfAlignments().

192  {
193  const unsigned int cellIndex(ESDetId(detId).denseIndex());
194  m_cellVec[cellIndex] = PreshowerStrip(f1, cornersMgr(), parm);
195  addValidID(detId);
196 }
A base class to handle the shape of preshower strips.
void addValidID(const DetId &id)
CaloCellGeometry::CornersMgr * cornersMgr()
static unsigned int EcalPreshowerGeometry::numberOfAlignments ( )
inlinestatic
unsigned int EcalPreshowerGeometry::numberOfParametersPerShape ( ) const
inlineoverridevirtual
unsigned int EcalPreshowerGeometry::numberOfShapes ( ) const
inlineoverridevirtual
unsigned int EcalPreshowerGeometry::numberOfTransformParms ( ) const
inlineoverridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 58 of file EcalPreshowerGeometry.h.

58 { return 3; }
bool EcalPreshowerGeometry::present ( const DetId id) const
inlineoverridevirtual

is this detid present in the geometry?

Reimplemented from CaloSubdetectorGeometry.

Definition at line 81 of file EcalPreshowerGeometry.h.

References CaloGenericDetId::denseIndex(), getGeometryRawPtr(), and m_cellVec.

Referenced by getClosestCellInPlane().

81  {
82  if (id == DetId(0))
83  return false;
84  // not needed???
85  auto index = CaloGenericDetId(id).denseIndex();
86  return index < m_cellVec.size();
87  }
uint32_t denseIndex() const
Definition: DetId.h:17
static std::string EcalPreshowerGeometry::producerTag ( )
inlinestatic

Definition at line 62 of file EcalPreshowerGeometry.h.

Referenced by PCaloGeometryBuilder::beginRun(), and CaloGeometryBuilder::CaloGeometryBuilder().

62 { return "EcalPreshower"; }
void EcalPreshowerGeometry::setzPlanes ( CCGFloat  z1minus,
CCGFloat  z2minus,
CCGFloat  z1plus,
CCGFloat  z2plus 
)

Definition at line 118 of file EcalPreshowerGeometry.cc.

References m_zplane.

Referenced by initializeParms(), and numberOfParametersPerShape().

118  {
119  assert(0 > z1minus && 0 > z2minus && 0 < z1plus && 0 < z2plus);
120 
121  m_zplane[0] = z1minus;
122  m_zplane[1] = z2minus;
123  m_zplane[2] = z1plus;
124  m_zplane[3] = z2plus;
125 }

Member Data Documentation

CellVec EcalPreshowerGeometry::m_cellVec
private

Definition at line 103 of file EcalPreshowerGeometry.h.

Referenced by getGeometryRawPtr(), newCell(), and present().

const CCGFloat EcalPreshowerGeometry::m_xInterLadGap
private

Definition at line 95 of file EcalPreshowerGeometry.h.

Referenced by getClosestCellInPlane().

const CCGFloat EcalPreshowerGeometry::m_xIntraLadGap
private

Definition at line 96 of file EcalPreshowerGeometry.h.

Referenced by getClosestCellInPlane().

const CCGFloat EcalPreshowerGeometry::m_xWidWaf
private

Definition at line 94 of file EcalPreshowerGeometry.h.

Referenced by getClosestCellInPlane().

const CCGFloat EcalPreshowerGeometry::m_yCtrOff
private

Definition at line 99 of file EcalPreshowerGeometry.h.

Referenced by getClosestCellInPlane().

const CCGFloat EcalPreshowerGeometry::m_yWidAct
private

Definition at line 98 of file EcalPreshowerGeometry.h.

Referenced by getClosestCellInPlane().

CCGFloat EcalPreshowerGeometry::m_zplane[4]
private