|
|
|
#include <RPCGeometryBuilderFromDDD.h>
Public Member Functions | |
| RPCGeometry * | build (const DDCompactView *cview, const MuonDDDConstants &muonConstants) |
| RPCGeometryBuilderFromDDD (bool comp11) | |
| ~RPCGeometryBuilderFromDDD () | |
Private Member Functions | |
| RPCGeometry * | buildGeometry (DDFilteredView &fview, const MuonDDDConstants &muonConstants) |
Private Attributes | |
| std::map< RPCDetId, std::list < RPCRoll * > > | chids |
| bool | theComp11Flag |
Build the RPCGeometry ftom the DDD description
Definition at line 23 of file RPCGeometryBuilderFromDDD.h.
| RPCGeometryBuilderFromDDD::RPCGeometryBuilderFromDDD | ( | bool | comp11 | ) |
Implementation of the RPC Geometry Builder from DDD
Definition at line 28 of file RPCGeometryBuilderFromDDD.cc.
: theComp11Flag(comp11) { }
| RPCGeometryBuilderFromDDD::~RPCGeometryBuilderFromDDD | ( | ) |
Definition at line 31 of file RPCGeometryBuilderFromDDD.cc.
{ }
| RPCGeometry * RPCGeometryBuilderFromDDD::build | ( | const DDCompactView * | cview, |
| const MuonDDDConstants & | muonConstants | ||
| ) |
Definition at line 34 of file RPCGeometryBuilderFromDDD.cc.
References DDFilteredView::addFilter(), DDSpecificsFilter::AND, buildGeometry(), align_tpl::filter, DDSpecificsFilter::matches, DDSpecificsFilter::setCriteria(), and relativeConstraints::value.
Referenced by RPCGeometryESModule::produce().
{
std::string attribute = "ReadOutName"; // could come from .orcarc
std::string value = "MuonRPCHits"; // could come from .orcarc
DDValue val(attribute, value, 0.0);
// Asking only for the MuonRPC's
DDSpecificsFilter filter;
filter.setCriteria(val, // name & value of a variable
DDSpecificsFilter::matches,
DDSpecificsFilter::AND,
true, // compare strings otherwise doubles
true // use merged-specifics or simple-specifics
);
DDFilteredView fview(*cview);
fview.addFilter(filter);
return this->buildGeometry(fview, muonConstants);
}
| RPCGeometry * RPCGeometryBuilderFromDDD::buildGeometry | ( | DDFilteredView & | fview, |
| const MuonDDDConstants & | muonConstants | ||
| ) | [private] |
Definition at line 54 of file RPCGeometryBuilderFromDDD.cc.
References RPCGeometry::add(), RPCChamber::add(), RPCNumberingScheme::baseNumberToUnitNumber(), chids, gather_cfg::cout, DDfetch(), cond::rpcobgas::detid, DDValue::doubles(), DDFilteredView::firstChild(), DDFilteredView::geoHistory(), MuonDDDNumbering::geoHistoryToBaseNumber(), geometry, RPCDetId::layer(), DDFilteredView::logicalPart(), python::multivaluedict::map(), DDName::name(), DDBase< N, C >::name(), mergeVDriftHistosByStation::name, DDFilteredView::nextSibling(), DDSolid::parameters(), pos, csvReporter::r, RPCDetId::region(), RPCDetId::ring(), TkRotation< T >::rotateAxes(), DDFilteredView::rotation(), GeomDetEnumerators::RPCBarrel, GeomDetEnumerators::RPCEndcap, RPCDetId::sector(), DDLogicalPart::solid(), DDFilteredView::specifics(), RPCDetId::station(), RPCDetId::subsector(), theComp11Flag, DDFilteredView::translation(), tablePrinter::width, x, detailsBasic3DVector::y, and z.
Referenced by build().
{
#ifdef LOCAL_DEBUG
std::cout <<"Building the geometry service"<<std::endl;
#endif
RPCGeometry* geometry = new RPCGeometry();
#ifdef LOCAL_DEBUG
std::cout << "About to run through the RPC structure" << std::endl;
std::cout <<" First logical part "
<<fview.logicalPart().name().name()<<std::endl;
#endif
bool doSubDets = fview.firstChild();
#ifdef LOCAL_DEBUG
std::cout << "doSubDets = " << doSubDets << std::endl;
#endif
while (doSubDets){
#ifdef LOCAL_DEBUG
std::cout <<"start the loop"<<std::endl;
#endif
// Get the Base Muon Number
MuonDDDNumbering mdddnum(muonConstants);
#ifdef LOCAL_DEBUG
std::cout <<"Getting the Muon base Number"<<std::endl;
#endif
MuonBaseNumber mbn=mdddnum.geoHistoryToBaseNumber(fview.geoHistory());
#ifdef LOCAL_DEBUG
std::cout <<"Start the Rpc Numbering Schema"<<std::endl;
#endif
// Get the The Rpc det Id
RPCNumberingScheme rpcnum(muonConstants);
int detid = 0;
#ifdef LOCAL_DEBUG
std::cout <<"Getting the Unit Number"<<std::endl;
#endif
detid = rpcnum.baseNumberToUnitNumber(mbn);
#ifdef LOCAL_DEBUG
std::cout <<"Getting the RPC det Id "<<detid <<std::endl;
#endif
RPCDetId rpcid(detid);
RPCDetId chid(rpcid.region(),rpcid.ring(),rpcid.station(),rpcid.sector(),rpcid.layer(),rpcid.subsector(),0);
#ifdef LOCAL_DEBUG
std::cout <<"The RPCDetid is "<<rpcid<<std::endl;
#endif
DDValue numbOfStrips("nStrips");
std::vector<const DDsvalues_type* > specs(fview.specifics());
std::vector<const DDsvalues_type* >::iterator is=specs.begin();
int nStrips=0;
for (;is!=specs.end(); is++){
if (DDfetch( *is, numbOfStrips)){
nStrips=int(numbOfStrips.doubles()[0]);
}
}
#ifdef LOCAL_DEBUG
if (nStrips == 0 )
std::cout <<"No strip found!!"<<std::endl;
#endif
std::vector<double> dpar=fview.logicalPart().solid().parameters();
std::string name=fview.logicalPart().name().name();
DDTranslation tran = fview.translation();
//removed .Inverse after comparing to DT...
DDRotationMatrix rota = fview.rotation();//.Inverse();
Surface::PositionType pos(tran.x()/cm,tran.y()/cm, tran.z()/cm);
// CLHEP way
// Surface::RotationType rot(rota.xx(),rota.xy(),rota.xz(),
// rota.yx(),rota.yy(),rota.yz(),
// rota.zx(),rota.zy(),rota.zz());
//ROOT::Math way
DD3Vector x, y, z;
rota.GetComponents(x,y,z);
// doesn't this just re-inverse???
Surface::RotationType rot (float(x.X()),float(x.Y()),float(x.Z()),
float(y.X()),float(y.Y()),float(y.Z()),
float(z.X()),float(z.Y()),float(z.Z()));
std::vector<float> pars;
RPCRollSpecs* rollspecs= 0;
Bounds* bounds = 0;
if (dpar.size()==3){
float width = dpar[0]/cm;
float length = dpar[1]/cm;
float thickness = dpar[2]/cm;
//RectangularPlaneBounds*
bounds =
new RectangularPlaneBounds(width,length,thickness);
pars.push_back(width);
pars.push_back(length);
pars.push_back(numbOfStrips.doubles()[0]); //h/2;
if (!theComp11Flag) {
//Correction of the orientation to get the REAL geometry.
//Change of axes for the +z part only.
//Including the 0 whell
if (tran.z() >-1500. ){
Basic3DVector<float> newX(-1.,0.,0.);
Basic3DVector<float> newY(0.,-1.,0.);
Basic3DVector<float> newZ(0.,0.,1.);
rot.rotateAxes (newX, newY,newZ);
}
}
rollspecs = new RPCRollSpecs(GeomDetEnumerators::RPCBarrel,name,pars);
#ifdef LOCAL_DEBUG
std::cout <<"Barrel "<<name
<<" par "<<width
<<" "<<length<<" "<<thickness;
#endif
}else{
float be = dpar[4]/cm;
float te = dpar[8]/cm;
float ap = dpar[0]/cm;
float ti = 0.4/cm;
// TrapezoidalPlaneBounds*
bounds =
new TrapezoidalPlaneBounds(be,te,ap,ti);
pars.push_back(dpar[4]/cm); //b/2;
pars.push_back(dpar[8]/cm); //B/2;
pars.push_back(dpar[0]/cm); //h/2;
pars.push_back(numbOfStrips.doubles()[0]); //h/2;
#ifdef LOCAL_DEBUG
std::cout <<"Forward "<<name
<<" par "<<dpar[4]/cm
<<" "<<dpar[8]/cm<<" "<<dpar[3]/cm<<" "
<<dpar[0];
#endif
rollspecs = new RPCRollSpecs(GeomDetEnumerators::RPCEndcap,name,pars);
//Change of axes for the forward
Basic3DVector<float> newX(1.,0.,0.);
Basic3DVector<float> newY(0.,0.,1.);
// if (tran.z() > 0. )
newY *= -1;
Basic3DVector<float> newZ(0.,1.,0.);
rot.rotateAxes (newX, newY,newZ);
}
#ifdef LOCAL_DEBUG
std::cout <<" Number of strips "<<nStrips<<std::endl;
#endif
BoundPlane* bp = new BoundPlane(pos,rot,bounds);
ReferenceCountingPointer<BoundPlane> surf(bp);
RPCRoll* r=new RPCRoll(rpcid,surf,rollspecs);
geometry->add(r);
std::list<RPCRoll *> rls;
if (chids.find(chid)!=chids.end()){
rls = chids[chid];
}
rls.push_back(r);
chids[chid]=rls;
doSubDets = fview.nextSibling(); // go to next layer
}
// Create the RPCChambers and store them on the Geometry
for( std::map<RPCDetId, std::list<RPCRoll *> >::iterator ich=chids.begin();
ich != chids.end(); ich++){
RPCDetId chid = ich->first;
std::list<RPCRoll * > rls = ich->second;
// compute the overall boundplane. At the moment we use just the last
// surface
BoundPlane* bp=0;
for(std::list<RPCRoll *>::iterator rl=rls.begin();
rl!=rls.end(); rl++){
const BoundPlane& bps = (*rl)->surface();
bp = const_cast<BoundPlane *>(&bps);
}
ReferenceCountingPointer<BoundPlane> surf(bp);
// Create the chamber
RPCChamber* ch = new RPCChamber (chid, surf);
// Add the rolls to rhe chamber
for(std::list<RPCRoll *>::iterator rl=rls.begin();
rl!=rls.end(); rl++){
ch->add(*rl);
}
// Add the chamber to the geometry
geometry->add(ch);
}
return geometry;
}
std::map<RPCDetId,std::list<RPCRoll *> > RPCGeometryBuilderFromDDD::chids [private] |
Definition at line 36 of file RPCGeometryBuilderFromDDD.h.
Referenced by buildGeometry().
bool RPCGeometryBuilderFromDDD::theComp11Flag [private] |
Definition at line 38 of file RPCGeometryBuilderFromDDD.h.
Referenced by buildGeometry().
1.7.3