|
|
|
#include <GEMGeometryBuilderFromDDD.h>
Public Member Functions | |
| GEMGeometry * | build (const DDCompactView *cview, const MuonDDDConstants &muonConstants) |
| GEMGeometryBuilderFromDDD (bool comp11) | |
| ~GEMGeometryBuilderFromDDD () | |
Private Member Functions | |
| GEMGeometry * | buildGeometry (DDFilteredView &fview, const MuonDDDConstants &muonConstants) |
Private Attributes | |
| bool | theComp11Flag |
Build the GEMGeometry ftom the DDD description
Definition at line 22 of file GEMGeometryBuilderFromDDD.h.
| GEMGeometryBuilderFromDDD::GEMGeometryBuilderFromDDD | ( | bool | comp11 | ) |
Implementation of the GEM Geometry Builder from DDD
Definition at line 28 of file GEMGeometryBuilderFromDDD.cc.
: theComp11Flag(comp11) { }
| GEMGeometryBuilderFromDDD::~GEMGeometryBuilderFromDDD | ( | ) |
Definition at line 31 of file GEMGeometryBuilderFromDDD.cc.
{ }
| GEMGeometry * GEMGeometryBuilderFromDDD::build | ( | const DDCompactView * | cview, |
| const MuonDDDConstants & | muonConstants | ||
| ) |
Definition at line 34 of file GEMGeometryBuilderFromDDD.cc.
References DDFilteredView::addFilter(), DDSpecificsFilter::AND, buildGeometry(), alcazmumu_cfi::filter, DDSpecificsFilter::matches, DDSpecificsFilter::setCriteria(), AlCaHLTBitMon_QueryRunRegistry::string, and relativeConstraints::value.
Referenced by GEMGeometryESModule::produce().
{
std::string attribute = "ReadOutName"; // could come from .orcarc
std::string value = "MuonGEMHits"; // could come from .orcarc
DDValue val(attribute, value, 0.0);
// Asking only for the MuonGEM'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);
}
| GEMGeometry * GEMGeometryBuilderFromDDD::buildGeometry | ( | DDFilteredView & | fview, |
| const MuonDDDConstants & | muonConstants | ||
| ) | [private] |
Definition at line 54 of file GEMGeometryBuilderFromDDD.cc.
References GEMGeometry::add(), GEMNumberingScheme::baseNumberToUnitNumber(), DDfetch(), cond::rpcobgas::detid, DDValue::doubles(), DDFilteredView::firstChild(), GeomDetEnumerators::GEM, DDFilteredView::geoHistory(), MuonDDDNumbering::geoHistoryToBaseNumber(), geometry, LogDebug, DDFilteredView::logicalPart(), DDName::name(), DDBase< N, C >::name(), mergeVDriftHistosByStation::name, DDFilteredView::nextSibling(), DDSolid::parameters(), pos, makeMuonMisalignmentScenario::rot, TkRotation< T >::rotateAxes(), DDFilteredView::rotation(), DDLogicalPart::solid(), DDFilteredView::specifics(), AlCaHLTBitMon_QueryRunRegistry::string, DDFilteredView::translation(), x, detailsBasic3DVector::y, and z.
Referenced by build().
{
LogDebug("GEMGeometryBuilderFromDDD") <<"Building the geometry service";
GEMGeometry* geometry = new GEMGeometry();
LogDebug("GEMGeometryBuilderFromDDD") << "About to run through the GEM structure\n"
<<" First logical part "
<<fview.logicalPart().name().name();
bool doSubDets = fview.firstChild();
LogDebug("GEMGeometryBuilderFromDDD") << "doSubDets = " << doSubDets;
while (doSubDets)
{
LogDebug("GEMGeometryBuilderFromDDD") <<"start the loop";
// Get the Base Muon Number
MuonDDDNumbering mdddnum(muonConstants);
LogDebug("GEMGeometryBuilderFromDDD") <<"Getting the Muon base Number";
MuonBaseNumber mbn = mdddnum.geoHistoryToBaseNumber(fview.geoHistory());
LogDebug("GEMGeometryBuilderFromDDD") <<"Start the GEM Numbering Schema";
// Get the The GEM det Id
GEMNumberingScheme gemnum(muonConstants);
int detid = 0;
LogDebug("GEMGeometryBuilderFromDDD") <<"Getting the Unit Number";
detid = gemnum.baseNumberToUnitNumber(mbn);
LogDebug("GEMGeometryBuilderFromDDD") <<"Getting the GEM det Id "<<detid;
GEMDetId gemmid(detid);
// GEMDetId chid(gemmid.region(),gemmid.ring(),gemmid.station(),gemmid.sector(),gemmid.layer(),gemmid.subsector(),0);
LogDebug("GEMGeometryBuilderFromDDD") <<"The GEMDetid is "<<gemmid;
DDValue numbOfStrips("nStrips");
DDValue numbOfPads("nPads");
std::vector<const DDsvalues_type* > specs(fview.specifics());
std::vector<const DDsvalues_type* >::iterator is = specs.begin();
double nStrips = 0., nPads = 0.;
for (;is != specs.end(); is++)
{
if (DDfetch( *is, numbOfStrips)) nStrips = numbOfStrips.doubles()[0];
if (DDfetch( *is, numbOfPads)) nPads = numbOfPads.doubles()[0];
}
LogDebug("GEMGeometryBuilderFromDDD") << ((nStrips == 0. ) ? ("No nStrips found!!") : (""));
LogDebug("GEMGeometryBuilderFromDDD") << ((nPads == 0. ) ? ("No nPads found!!") : (""));
//std::cout <<"# strips, pads: "<< nStrips<<", "<<nPads<<std::endl;
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()));
Bounds* bounds = 0;
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);
std::vector<float> pars;
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(nStrips);
pars.push_back(nPads);
LogDebug("GEMGeometryBuilderFromDDD") <<"GEM "<<name
<<" par "<<dpar[4]/cm
<<" "<<dpar[8]/cm<<" "<<dpar[3]/cm<<" "
<<dpar[0];
GEMEtaPartitionSpecs* e_p_specs = new GEMEtaPartitionSpecs(GeomDetEnumerators::GEM, 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);
LogDebug("GEMGeometryBuilderFromDDD") <<" Number of strips "<<nStrips<<" Number of pads "<<nPads;
BoundPlane* bp = new BoundPlane(pos, rot, bounds);
ReferenceCountingPointer<BoundPlane> surf(bp);
GEMEtaPartition* gep = new GEMEtaPartition(gemmid, surf, e_p_specs);
geometry->add(gep);
std::list<GEMEtaPartition *> gepls;
/*
if (chids.find(chid)!=chids.end()){
gepls = chids[chid];
}
*/
gepls.push_back(gep);
//chids[chid]=gepls;
doSubDets = fview.nextSibling(); // go to next layer
}
/*
// Create the GEMChambers and store them on the Geometry
for( std::map<GEMDetId, std::list<GEMEtaPartition *> >::iterator ich=chids.begin();
ich != chids.end(); ich++){
GEMDetId chid = ich->first;
std::list<GEMEtaPartition * > gepls = ich->second;5D
// compute the overall boundplane. At the moment we use just the last
// surface
BoundPlane* bp=0;
for(std::list<GEMEtaPartition *>::iterator gepl=gepls.begin();
gepl!=gepls.end(); gepl++){
const BoundPlane& bps = (*gepl)->surface();
bp = const_cast<BoundPlane *>(&bps);
}
ReferenceCountingPointer<BoundPlane> surf(bp);
// Create the chamber
GEMChamber* ch = new GEMChamber (chid, surf);
// Add the etapartitions to the chamber
for(std::list<GEMEtaPartition *>::iterator gepl=gepls.begin();
gepl!=gepls.end(); rl++){
ch->add(*gepl);
}
// Add the chamber to the geometry
geometry->add(ch);
}
*/
return geometry;
}
bool GEMGeometryBuilderFromDDD::theComp11Flag [private] |
Definition at line 37 of file GEMGeometryBuilderFromDDD.h.
1.7.3