22 if (ieta==0)
continue;
23 for (
int iphi=1; iphi<=72; iphi++) {
35 if (ieta==0)
continue;
36 for (
int iphi=1; iphi<=72; iphi++) {
43 edm::LogInfo(
"Geometry") <<
"CaloTowersHardcodeGeometry made " << n <<
" towers.";
44 return std::auto_ptr<CaloSubdetectorGeometry>(
geom);
51 const double EBradius = 143.0;
52 const double HOradius = 406.0+1.0;
53 const double EEz = 320.0;
54 const double HEz = 568.0;
55 const double HFz = 1100.0;
56 const double HFthick = 165;
59 int etaRing=
abs(ieta);
60 int sign=(ieta>0)?(1):(-1);
69 double eta = 0.5*(eta1+eta2);
70 double deta = (eta2-eta1);
76 double phi_low = dphi_nominal*(iphi-1);
77 double phi = phi_low+dphi_half;
79 double x,
y,
z,thickness;
86 thickness=HFthick/tanh(eta);
87 }
else if (etaRing>17) {
92 thickness=(HEz-EEz)/tanh(eta);
94 x=EBradius *
cos(phi);
95 y=EBradius *
sin(phi);
97 z=EBradius * sinh(eta);
98 thickness=(HOradius-EBradius) * cosh(eta);
104 const double mysign ( !alongZ ? 1 : -1 ) ;
105 std::vector<CCGFloat> hh ;
107 hh.push_back( deta/2 ) ;
108 hh.push_back( dphi_half ) ;
109 hh.push_back( mysign*thickness/2. ) ;
111 hh.push_back( fabs( eta ) ) ;
112 hh.push_back( fabs( z ) ) ;
114 geom->
newCell( point, point, point,
void makeCell(int ieta, int iphi, CaloSubdetectorGeometry *geom) const
virtual void newCell(const GlobalPoint &f1, const GlobalPoint &f2, const GlobalPoint &f3, const CCGFloat *parm, const DetId &detId)=0
Sin< T >::type sin(const T &t)
int nPhiBins(int etaRing) const
how many phi segments in this ring
void allocatePar(ParVec::size_type n, unsigned int m)
std::auto_ptr< CaloSubdetectorGeometry > load()
Cos< T >::type cos(const T &t)
CaloCellGeometry::CCGFloat CCGFloat
static const CCGFloat * getParmPtr(const std::vector< CCGFloat > &vd, ParMgr *mgr, ParVecVec &pvv)
int firstHEDoublePhiRing() const
static const double theHFEtaBounds[]
CaloCellGeometry::CornersMgr * cornersMgr()
int firstHFQuadPhiRing() const
void allocateCorners(CaloCellGeometry::CornersVec::size_type n)
static const double theHBHEEtaBounds[]
*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