30{
31
32 std::cout << " +---------------------------------------------------+\n"
33 << " | Creating Mdc Geometry information |\n"
34 << " +---------------------------------------------------+\n";
35
36 IService* ser;
37 StatusCode sc = Gaudi::svcLocator()->getService("MdcGeomSvc",ser);
38 if (!sc.isSuccess())
39 std::cout <<" MdcDetector::Could not open Geometry Service"<<std::endl;
40 MdcGeomSvc *mdcsvc = dynamic_cast<MdcGeomSvc*> (ser);
41 if(!mdcsvc) std::cout <<"MdcDetector::Could not open Geometry Service"<<std::endl;
42
46
47 _senseWire = new MdcSWire* [_nSWire];
48 _dclayer = new MdcLayer* [_nLayer];
49 _slayList = new MdcSuperLayer* [_nSlay];
50
51
52 for( int iwire = 0; iwire < _nSWire; iwire++ ) {
53 const MdcGeoWire *geowir = mdcsvc->
Wire(iwire);
56
57 double sag=0.;
58 if(_doSag) sag = geowir->
Sag()/10.;
59
60
61 _senseWire[iwire] =
new MdcSWire(eastP, westP, sag, geowir->
Id(), geowir->
Cell());
62
63
64 }
65
66
67 for( int ilay = 0; ilay < _nLayer; ilay++ ) {
68 const MdcGeoLayer *geolay = mdcsvc->
Layer(ilay);
69 int nwir = geolay->
NCell();
70 int firstwir = geolay->
Wirst();
71 _dclayer[ilay] = new MdcLayer(ilay, nwir, &_senseWire[firstwir],*this);
72 _wires_in_layer[ilay] = nwir;
73 }
74
75
76 buildpointers();
77
78
79
80 buildSuperLayers();
81
82 double rOther;
83 for(int ilay = 0; ilay < _nLayer; ilay++) {
84 if( ilay == 0) {
86 } else {
88 }
89 double height = fabs(
Layer(ilay)->rMid() - rOther);
90 const_cast<MdcLayer*
>(
Layer(ilay))->setCellHeight(height);
91 }
92}
HepGeom::Point3D< double > HepPoint3D
const MdcLayer * Layer(unsigned id) const
const double Sag(void) const
HepPoint3D Forward(void) const
HepPoint3D Backward(void) const
const MdcGeoWire *const Wire(unsigned id)
const int getSuperLayerSize()
const MdcGeoLayer *const Layer(unsigned id)
1.13.2