MuonGM::RpcLayer Class Reference

#include <RpcLayer.h>

Inheritance diagram for MuonGM::RpcLayer:

Collaboration diagram for MuonGM::RpcLayer:

Public Member Functions
	RpcLayer (const std::string &s, Rpc *t)
PVLink	build (StoredMaterialManager &matManager, const MYSQL &mysql)
PVLink	build (StoredMaterialManager &matManager, const MYSQL &mysql, int cutoutson, const std::vector< Cutout * > &)
virtual void	print () const override
void	setLogVolName (const std::string &str)

Public Attributes
double	lwidth {0.}
double	llength {0.}
double	thickness {0.}
double	llongWidth {0.}
Rpc *	m {nullptr}
std::string	name {}
std::string	logVolName {}

Detailed Description

Definition at line 22 of file RpcLayer.h.

Constructor & Destructor Documentation

RpcLayer()

MuonGM::RpcLayer::RpcLayer	(	const std::string &	s,
		Rpc *	t )

Definition at line 39 of file RpcLayer.cxx.

: DetectorElement(s) { m = t; }

Member Function Documentation

build() [1/2]

PVLink MuonGM::RpcLayer::build	(	StoredMaterialManager &	matManager,
		const MYSQL &	mysql )

Definition at line 41 of file RpcLayer.cxx.

                                                     {
        std::vector<Cutout *> vcutdef;
        int cutoutson = 0;
        return build(matManager, mysql, cutoutson, vcutdef);
    }

build() [2/2]

PVLink MuonGM::RpcLayer::build	(	StoredMaterialManager &	matManager,
		const MYSQL &	mysql,
		int	cutoutson,
		const std::vector< Cutout * > &	vcutdef )

Definition at line 48 of file RpcLayer.cxx.

                                                                              {
        MsgStream log(Athena::getMessageSvc(), "MuGM::GeoVPhysVol::build");
 
        double eps = 0.000001;
        double tol = 1.e-6;
 
        const RPC *r = dynamic_cast<const RPC*>(mysql.GetTechnology(name));
 
        if (m->nGasGaps() == 3 && r->NstripPanels_in_s != 1)
            throw std::runtime_error(Form("File: %s, Line: %d\nRpcLayer::build() - NstripPanels_in_s = %d for BI RPC, not possible", __FILE__, __LINE__, r->NstripPanels_in_s));
 
        double thickness = r->rpcLayerThickness - tol;
        if (m->nGasGaps() == 3)
            thickness = rpc3GapLayerThickness - tol; // for BI RPCs
        double length = m->length;
        double width = m->width;
 
        if (RPCprint) {
            log << MSG::VERBOSE << " RpcLayer: building an RPC layer with width = " << width << " length = " << length << " and thickness = " << thickness << endmsg;
        }
 
        const GeoShape *srpcl = new GeoTrd(thickness / 2, thickness / 2, width / 2, width / 2, length / 2);
        const GeoMaterial *mrpcl = matManager.getMaterial("std::Air");
        GeoLogVol *lrpcl = new GeoLogVol("Rpclayer", srpcl, mrpcl);
        PVLink prpcl = new GeoPhysVol(lrpcl);
 
        double newpos = -thickness / 2.;
 
        // phi strip panel(s)
        double strpanThickness = r->stripPanelThickness - tol;
        double strpanCopperThickness = r->stripPanelCopperSkinThickness;
        if (m->nGasGaps() == 3) { // for BI RPCs
            strpanThickness = rpc3GapStrPanThickness - tol;
            strpanCopperThickness = rpc3GapStrPanCopThickness;
        }
        double strpanLength = length / r->NstripPanels_in_z;
        double strpanWidth = width / r->NstripPanels_in_s - eps;
        if (RPCprint) {
            log << MSG::VERBOSE << " RpcLayer:: Building a strip panel:: div_s, div_z = " << r->NstripPanels_in_s << " " << r->NstripPanels_in_z << endmsg;
        }
        if (RPCprint) {
            log << MSG::VERBOSE << " RpcLayer:: Building a strip panel:: w,l,t = " << strpanWidth << " " << strpanLength << " " << strpanThickness << endmsg;
        }
        GeoTrd *sstrpan = new GeoTrd(strpanThickness / 2., strpanThickness / 2., strpanWidth / 2., strpanWidth / 2., strpanLength / 2.);
 
        GeoTrd *sfoamstrpan = new GeoTrd(strpanThickness / 2. - strpanCopperThickness, strpanThickness / 2. - strpanCopperThickness, strpanWidth / 2. - strpanCopperThickness,
                                         strpanWidth / 2. - strpanCopperThickness, strpanLength / 2. - strpanCopperThickness);
        const GeoShape *scustrpan = sstrpan;
 
        const auto *stripMaterial = matManager.getMaterial("muo::RpcFoam");
        if (m->nGasGaps() == 3) { // for BI RPCs
            stripMaterial = matManager.getMaterial("muo::Forex");
        }
 
        GeoLogVol *lcustrpan = new GeoLogVol("RPC_StripPanelCuSkin", scustrpan, matManager.getMaterial("std::Copper"));
        GeoLogVol *lfoamstrpan = new GeoLogVol("RPC_StripPanelFoam", sfoamstrpan, stripMaterial);
 
        newpos += strpanThickness / 2. + tol / 2.;
        GeoPhysVol *pcustrpan11 = new GeoPhysVol(lcustrpan);
        GeoPhysVol *pfoamstrpan11 = new GeoPhysVol(lfoamstrpan);
        GeoTransform *tx = new GeoTransform(GeoTrf::TranslateX3D(newpos));
        GeoTransform *ty1;
        GeoTransform *ty2;
 
        if (r->NstripPanels_in_s == 2) {
            if (RPCprint) {
                log << MSG::VERBOSE << " RpcLayer::NstripPanels_in_s == 2 " << endmsg;
            }
            GeoPhysVol *pcustrpan12 = new GeoPhysVol(lcustrpan);
            GeoPhysVol *pfoamstrpan12 = new GeoPhysVol(lfoamstrpan);
            ty1 = new GeoTransform(GeoTrf::TranslateY3D(-width / 4.));
            ty2 = new GeoTransform(GeoTrf::TranslateY3D(width / 4.));
            prpcl->add(tx);
            prpcl->add(ty1);
            if (RPCprint) {
                log << MSG::VERBOSE << "RpcLayer:: Locating the 1st phi strip panel at x, y " << newpos << " " << -width / 4. << " of width =" << strpanWidth
                    << " in a box of width =" << width << endmsg;
            }
            prpcl->add(pcustrpan11);
            pcustrpan11->add(pfoamstrpan11);
            prpcl->add(tx);
            prpcl->add(ty2);
            if (RPCprint) {
                log << MSG::VERBOSE << "RpcLayer:: Locating the 2nd phi strip panel at x, y " << newpos << " " << width / 4. << " of width =" << strpanWidth
                    << " in a box of width =" << width << endmsg;
            }
            prpcl->add(pcustrpan12);
            pcustrpan12->add(pfoamstrpan12);
 
        } else if (r->NstripPanels_in_s == 1) {
            if (RPCprint) {
                log << MSG::VERBOSE << " RpcLayer::NstripPanels_in_s == 1 " << endmsg;
            }
            prpcl->add(tx);
            prpcl->add(pcustrpan11);
            pcustrpan11->add(pfoamstrpan11);
            if (RPCprint) {
                log << MSG::VERBOSE << " RpcLayer:: Locating a single phi strip panel at x, y " << newpos << " 0 "
                    << " of width = " << strpanWidth << " in a box of width = " << width << endmsg;
            }
        }
 
        newpos += strpanThickness / 2.;
        double bakelThickness = r->bakeliteThickness;
        double gThickness = r->gasThickness;
        double totAirThickness = r->totalAirThickness;
        if (m->nGasGaps() == 3) { // for BI RPCs
            bakelThickness = rpc3GapBakelThickness;
            gThickness = rpc3GapGThickness;
            totAirThickness = rpc3GapTotAirThickness;
        }
 
        double ggThickness = 2. * bakelThickness + gThickness + totAirThickness;
        newpos += ggThickness / 2.;
        // this brings to the center of the gas gap
 
        double ggLength = length;
        double ggWidth = width / r->NGasGaps_in_s - eps;
        double gasLength = ggLength - 2. * r->bakeliteframesize;
        double gasWidth = ggWidth - 2. * r->bakeliteframesize;
 
        if (m->nGasGaps() == 3) {             // for BI RPCs
            if (name == "RPC26") {            // big RPC7
                gasLength = ggLength - 93.25; // ggLength - deadframesizeEta
                gasWidth = ggWidth - 109.52;  // ggWidth - deadframesizePhi
            } else if (name == "RPC27") {     // small RPC7
                gasLength = ggLength - 93.12; // ggLength - deadframesizeEta
                gasWidth = ggWidth - 109.52;  // ggWidth - deadframesizePhi
            } else if (name == "RPC28") {     // big RPC8
                gasLength = ggLength - 93.04; // ggLength - deadframesizeEta
                gasWidth = ggWidth - 109.52;  // ggWidth - deadframesizePhi
            } else if (name == "RPC29") {     // small RPC8
                gasLength = ggLength - 93.04; // ggLength - deadframesizeEta
                gasWidth = ggWidth - 109.2;   // ggWidth - deadframesizePhi
            }
        }
 
        if (RPCprint) {
            log << MSG::VERBOSE << "RpcLayer:: Building the gasgap:: " << r->NGasGaps_in_s << " in s direction" << endmsg;
            log << MSG::VERBOSE << "RpcLayer:: Building the gasgap:: w,l,t " << ggWidth << " " << ggLength << " " << ggThickness << endmsg;
            log << MSG::VERBOSE << "RpcLayer:: Building the gas:: w,l,t " << gasWidth << " " << gasLength << " " << gThickness << endmsg;
        }
        GeoTrd *sgg = new GeoTrd(ggThickness / 2., ggThickness / 2., ggWidth / 2., ggWidth / 2., ggLength / 2.);
        GeoTrd *sgas = new GeoTrd(gThickness / 2., gThickness / 2., gasWidth / 2., gasWidth / 2., gasLength / 2.);
        const GeoShape *sbak = sgg;
        GeoLogVol *lbak = new GeoLogVol("gas volume:" + MuonGM::buildString(r->NGasGaps_in_s, 0) + "gg_in_s_" + MuonGM::buildString(r->NstripPanels_in_s, 0) + "sp_in_s", sbak,
                                        matManager.getMaterial("std::Bakelite"));
        GeoLogVol *lgas = new GeoLogVol("gazGap", sgas, matManager.getMaterial("muo::RPCgas"));
        GeoPhysVol *pbak1 = new GeoPhysVol(lbak);
        GeoPhysVol *pgas1 = new GeoPhysVol(lgas);
        tx = new GeoTransform(GeoTrf::TranslateX3D(newpos));
 
        if (r->NGasGaps_in_s == 2) {
            GeoPhysVol *pbak2 = new GeoPhysVol(lbak);
            GeoPhysVol *pgas2 = new GeoPhysVol(lgas);
            ty1 = new GeoTransform(GeoTrf::TranslateY3D(-width / 4.));
            ty2 = new GeoTransform(GeoTrf::TranslateY3D(width / 4.));
            prpcl->add(tx);
            prpcl->add(ty1);
            if (RPCprint) {
                log << MSG::VERBOSE << "RpcLayer:: put 1st gas gap centre at depth, s " << newpos << " " << -width / 4. << endmsg;
            }
            prpcl->add(new GeoIdentifierTag(0));
            prpcl->add(pbak1);
            pbak1->add(new GeoIdentifierTag(1));
            pbak1->add(pgas1);
            prpcl->add(tx);
            prpcl->add(ty2);
            if (RPCprint) {
                log << MSG::VERBOSE << "RpcLayer:: put 2nd gas gap centre at depth, s " << newpos << " " << width / 4. << endmsg;
            }
            prpcl->add(new GeoIdentifierTag(10));
            prpcl->add(pbak2);
            pbak2->add(new GeoIdentifierTag(1));
            pbak2->add(pgas2);
        } else if (r->NGasGaps_in_s == 1) {
            prpcl->add(tx);
            if (m->nGasGaps() == 3) { // for BI RPCs
                GeoTransform *ty = new GeoTransform(GeoTrf::TranslateY3D(m->y_translation));
                GeoTransform *tz = new GeoTransform(GeoTrf::TranslateZ3D(m->z_translation));
                prpcl->add(ty);
                prpcl->add(tz);
            }
            prpcl->add(new GeoIdentifierTag(0));
            prpcl->add(pbak1);
            if (RPCprint) {
                log << MSG::VERBOSE << "RpcLayer:: put a single gas gap at depth, s " << newpos << " 0 " << endmsg;
            }
            pbak1->add(new GeoIdentifierTag(1));
            pbak1->add(pgas1);
        }
 
        newpos += ggThickness / 2.;
 
        // now eta strip panel
        newpos += strpanThickness / 2.;
 
        GeoPhysVol *pcustrpan21 = new GeoPhysVol(lcustrpan);
        GeoPhysVol *pfoamstrpan21 = new GeoPhysVol(lfoamstrpan);
        tx = new GeoTransform(GeoTrf::TranslateX3D(newpos));
 
        if (r->NstripPanels_in_s == 2) {
            GeoPhysVol *pcustrpan22 = new GeoPhysVol(lcustrpan);
            GeoPhysVol *pfoamstrpan22 = new GeoPhysVol(lfoamstrpan);
            ty1 = new GeoTransform(GeoTrf::TranslateY3D(-width / 4.));
            ty2 = new GeoTransform(GeoTrf::TranslateY3D(width / 4.));
            prpcl->add(tx);
            prpcl->add(ty1);
            if (RPCprint) {
                log << MSG::VERBOSE << "RpcLayer:: Locating the 1st eta panel at x, y " << newpos << " " << -width / 4. << " of width =" << strpanWidth
                    << " in a box of width =" << width << endmsg;
            }
            prpcl->add(pcustrpan21);
            pcustrpan21->add(pfoamstrpan21);
            prpcl->add(tx);
            prpcl->add(ty2);
            if (RPCprint) {
                log << MSG::VERBOSE << "RpcLayer:: Locating the 2nd eta panel at x, y " << newpos << " " << width / 4. << " of width =" << strpanWidth
                    << " in a box of width =" << width << endmsg;
            }
            prpcl->add(pcustrpan22);
            pcustrpan22->add(pfoamstrpan22);
        } else if (r->NstripPanels_in_s == 1) {
            prpcl->add(tx);
            if (RPCprint) {
                log << MSG::VERBOSE << "RpcLayer:: Locating a single eta panel at x, y " << newpos << " 0 " << endmsg;
            }
            prpcl->add(pcustrpan21);
            pcustrpan21->add(pfoamstrpan21);
        }
 
        // Apply cutouts
        if (cutoutson && !vcutdef.empty()) {
            PVLink tempPhys = nullptr;
            Cutout *cut = nullptr;
            GeoShape *cutoutShape = nullptr;
            GeoTrf::Transform3D cutTrans{GeoTrf::Transform3D::Identity()};
 
            for (unsigned i = 0; i < vcutdef.size(); i++) {
                cut = vcutdef[i];
                cutoutShape = new GeoTrd(thickness / 2. + 1., thickness / 2. + 1., cut->widthXs / 2. + 0.5, cut->widthXl / 2. + 0.5, cut->lengthY / 2. + tol);
                cutTrans = GeoTrf::Translate3D(0.0, cut->dx, -length / 2 + cut->dy + cut->lengthY / 2.);
 
                PVLink  volToCut{prpcl};
                GeoCutVolAction cutAction(*cutoutShape, cutTrans);
                volToCut->apply(&cutAction);
                tempPhys = cutAction.getPV();
                prpcl = tempPhys;
            }
        }
 
        return prpcl;
    }

print()

void MuonGM::RpcLayer::print ( ) const

overridevirtual

Implements MuonGM::DetectorElement.

Definition at line 303 of file RpcLayer.cxx.

                               {
        MsgStream log(Athena::getMessageSvc(), "MuGM::GeoVPhysVol");
        log << MSG::INFO << " Rpc Layer " << name + "Layer"
            << " :" << endmsg;
    }

setLogVolName()

void MuonGM::DetectorElement::setLogVolName ( const std::string & str )

inlineinherited

Definition at line 22 of file DetectorElement.h.

{ logVolName = str; }

Member Data Documentation

llength

double MuonGM::RpcLayer::llength {0.}

Definition at line 26 of file RpcLayer.h.

{0.};

llongWidth

double MuonGM::RpcLayer::llongWidth {0.}

Definition at line 28 of file RpcLayer.h.

{0.}; // for trapezoidal layers

logVolName

std::string MuonGM::DetectorElement::logVolName {}

inherited

Definition at line 18 of file DetectorElement.h.

{};

lwidth

double MuonGM::RpcLayer::lwidth {0.}

Definition at line 25 of file RpcLayer.h.

{0.};

m

Rpc* MuonGM::RpcLayer::m {nullptr}

Definition at line 30 of file RpcLayer.h.

{nullptr};

name

std::string MuonGM::DetectorElement::name {}

inherited

Definition at line 17 of file DetectorElement.h.

{};

thickness

double MuonGM::RpcLayer::thickness {0.}

Definition at line 27 of file RpcLayer.h.

{0.};

---

The documentation for this class was generated from the following files:

• RpcLayer.h
• RpcLayer.cxx