36 std::cout <<
" Geometry node for RPCGeom is " << &(*pDD) << std::endl;
37 cout <<
" I have " << pDD->detTypes().size() <<
" detTypes" << endl;
38 cout <<
" I have " << pDD->detUnits().size() <<
" detUnits" << endl;
39 cout <<
" I have " << pDD->dets().size() <<
" dets" << endl;
40 cout <<
" I have " << pDD->rolls().size() <<
" rolls" << endl;
41 cout <<
" I have " << pDD->chambers().size() <<
" chambers" << endl;
43 std::cout <<
myName() <<
": Begin iteration over geometry..." << std::endl;
48 for (TrackingGeometry::DetContainer::const_iterator it = pDD->dets().begin(); it != pDD->dets().end(); it++) {
51 if (dynamic_cast<const RPCChamber*>(*it) !=
nullptr) {
53 const RPCChamber* ch = dynamic_cast<const RPCChamber*>(*it);
55 std::vector<const RPCRoll*> rollsRaf = (ch->
rolls());
56 for (std::vector<const RPCRoll*>::iterator
r = rollsRaf.begin();
r != rollsRaf.end(); ++
r) {
58 std::cout <<
"Region = " << (*r)->id().region() <<
" Ring = " << (*r)->id().ring()
59 <<
" Station = " << (*r)->id().station() <<
" Sector = " << (*r)->id().sector()
60 <<
" Layer = " << (*r)->id().layer() <<
" Subsector = " << (*r)->id().subsector()
61 <<
" Roll = " << (*r)->id().roll() << std::endl;
64 std::cout <<
s.name() <<
" eta partition " <<
s.eta_partition() <<
" nroll=" << ch->
nrolls() <<
" z=" <<
g.z()
65 <<
" phi=" <<
g.phi() <<
" R=" <<
g.perp() << std::endl;
67 if ((*r)->id().region() == 0) {
70 double cmin = cic.first;
71 double cmax = cic.second;
76 std::pair<double, double> cic2(cmin, cmax);
79 std::pair<double, double> cic(
g.z(),
g.z());
82 }
else if ((*r)->id().region() == +1) {
85 double cmin = cic.first;
86 double cmax = cic.second;
91 std::pair<double, double> cic2(cmin, cmax);
94 std::pair<double, double> cic(
g.perp(),
g.perp());
97 }
else if ((*r)->id().region() == -1) {
100 double cmin = cic.first;
101 double cmax = cic.second;
106 std::pair<double, double> cic2(cmin, cmax);
109 std::pair<double, double> cic(
g.perp(),
g.perp());
118 std::map<int, std::pair<double, double> >::iterator
ieta;
121 std::cout <<
" Eta " <<
ieta->first <<
" Radii = ( " <<
ieta->second.first <<
", " <<
ieta->second.second <<
")"
126 std::cout <<
" Eta " <<
ieta->first <<
" Z = ( " <<
ieta->second.first <<
", " <<
ieta->second.second <<
")"
131 std::cout <<
" Eta " <<
ieta->first <<
" Radii = ( " <<
ieta->second.first <<
", " <<
ieta->second.second <<
")"