|
|
|
#include <DDTIBLayerAlgo.h>
Public Member Functions | |
| DDTIBLayerAlgo () | |
| void | execute (DDCompactView &cpv) |
| void | initialize (const DDNumericArguments &nArgs, const DDVectorArguments &vArgs, const DDMapArguments &mArgs, const DDStringArguments &sArgs, const DDStringVectorArguments &vsArgs) |
| virtual | ~DDTIBLayerAlgo () |
Private Attributes | |
| double | bwExtPillarDPhi |
| double | bwExtPillarDz |
| std::vector< double > | bwExtPillarPhi |
| std::vector< double > | bwExtPillarZ |
| double | bwIntPillarDPhi |
| double | bwIntPillarDz |
| std::vector< double > | bwIntPillarPhi |
| std::vector< double > | bwIntPillarZ |
| std::string | centMat |
| std::vector< double > | centRing1par |
| std::vector< double > | centRing2par |
| double | cylinderInR |
| std::string | cylinderMat |
| double | cylinderT |
| std::string | detectorLo |
| double | detectorTilt |
| std::string | detectorUp |
| double | dohmAuxL |
| std::string | dohmAuxMaterial |
| std::string | dohmAuxName |
| std::string | dohmCableMaterial |
| std::string | dohmCarrierMaterial |
| double | dohmCarrierPhiOff |
| std::vector< double > | dohmListBW |
| std::vector< double > | dohmListFW |
| double | dohmPrimL |
| std::string | dohmPrimMaterial |
| std::string | dohmPrimName |
| double | dohmtoMF |
| double | fillerDz |
| std::string | fillerMat |
| double | fwExtPillarDPhi |
| double | fwExtPillarDz |
| std::vector< double > | fwExtPillarPhi |
| std::vector< double > | fwExtPillarZ |
| double | fwIntPillarDPhi |
| double | fwIntPillarDz |
| std::vector< double > | fwIntPillarPhi |
| std::vector< double > | fwIntPillarZ |
| std::string | genMat |
| std::string | idNameSpace |
| double | layerL |
| std::string | MFExtRingMat |
| std::string | MFIntRingMat |
| double | MFRingDz |
| double | MFRingInR |
| double | MFRingOutR |
| double | MFRingT |
| double | phioffLo |
| double | phioffUp |
| std::string | pillarMaterial |
| double | radiusLo |
| double | radiusUp |
| std::string | ribMat |
| std::vector< double > | ribPhi |
| std::vector< double > | ribW |
| int | stringsLo |
| int | stringsUp |
| double | supportT |
Definition at line 10 of file DDTIBLayerAlgo.h.
| DDTIBLayerAlgo::DDTIBLayerAlgo | ( | ) |
| DDTIBLayerAlgo::~DDTIBLayerAlgo | ( | ) | [virtual] |
Definition at line 25 of file DDTIBLayerAlgo.cc.
{}
| void DDTIBLayerAlgo::execute | ( | DDCompactView & | cpv | ) |
Definition at line 176 of file DDTIBLayerAlgo.cc.
References abs, bwExtPillarDPhi, bwExtPillarDz, bwExtPillarPhi, bwExtPillarZ, bwIntPillarDPhi, bwIntPillarDz, bwIntPillarPhi, bwIntPillarZ, centMat, centRing1par, centRing2par, funct::cos(), cylinderInR, cylinderMat, cylinderT, dbl_to_string(), DDanonymousRot(), DDcreateRotationMatrix(), DDBase< N, C >::ddname(), DDrot(), DDSplit(), detectorLo, detectorTilt, detectorUp, dohmAuxL, dohmAuxName, dohmCarrierMaterial, dohmCarrierPhiOff, dohmListBW, dohmListFW, dohmPrimL, dohmPrimName, dohmtoMF, fillerDz, fillerMat, first, fwExtPillarDPhi, fwExtPillarDz, fwExtPillarPhi, fwExtPillarZ, fwIntPillarDPhi, fwIntPillarDz, fwIntPillarPhi, fwIntPillarZ, genMat, i, idNameSpace, j, layerL, LogDebug, MFExtRingMat, MFIntRingMat, MFRingDz, MFRingInR, MFRingOutR, MFRingT, n, DDName::name(), DDBase< N, C >::name(), mergeVDriftHistosByStation::name, dbtoconf::parent, phi, pillarMaterial, DDCompactView::position(), radiusLo, radiusUp, ribMat, ribPhi, ribW, idealTransformation::rotation, edm::second(), funct::sin(), stringsLo, stringsUp, supportT, theta(), DDSolidFactory::tubs(), and tablePrinter::width.
{
LogDebug("TIBGeom") << "==>> Constructing DDTIBLayerAlgo...";
DDName parentName = parent().name();
const std::string &idName = parentName.name();
double rmin = MFRingInR;
double rmax = MFRingOutR;
DDSolid solid = DDSolidFactory::tubs(DDName(idName, idNameSpace), 0.5*layerL,
rmin, rmax, 0, CLHEP::twopi);
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: "
<< DDName(idName,idNameSpace) << " Tubs made of "
<< genMat << " from 0 to " << CLHEP::twopi/CLHEP::deg
<< " with Rin " << rmin << " Rout " << rmax
<< " ZHalf " << 0.5*layerL;
DDName matname(DDSplit(genMat).first, DDSplit(genMat).second);
DDMaterial matter(matname);
DDLogicalPart layer(solid.ddname(), matter, solid);
//Internal layer first
double rin = rmin+MFRingT;
// double rout = 0.5*(radiusLo+radiusUp-cylinderT);
double rout = cylinderInR;
std::string name = idName + "Down";
solid = DDSolidFactory::tubs(DDName(name, idNameSpace), 0.5*layerL,
rin, rout, 0, CLHEP::twopi);
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: "
<< DDName(name, idNameSpace) << " Tubs made of "
<< genMat << " from 0 to " << CLHEP::twopi/CLHEP::deg
<< " with Rin " << rin << " Rout " << rout
<< " ZHalf " << 0.5*layerL;
DDLogicalPart layerIn(solid.ddname(), matter, solid);
cpv.position(layerIn, layer, 1, DDTranslation(0.0, 0.0, 0.0), DDRotation());
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: " << layerIn.name()
<< " number 1 positioned in " << layer.name()
<< " at (0,0,0) with no rotation";
double rposdet = radiusLo;
double dphi = CLHEP::twopi/stringsLo;
DDName detIn(DDSplit(detectorLo).first, DDSplit(detectorLo).second);
for (int n = 0; n < stringsLo; n++) {
double phi = (n+0.5)*dphi;
double phix = phi - detectorTilt + 90*CLHEP::deg;
double phideg = phix/CLHEP::deg;
DDRotation rotation;
if (phideg != 0) {
double theta = 90*CLHEP::deg;
double phiy = phix + 90.*CLHEP::deg;
std::string rotstr = idName + dbl_to_string(phideg*10.);
rotation = DDRotation(DDName(rotstr, idNameSpace));
if (!rotation) {
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: Creating a new "
<< "rotation: " << rotstr << "\t90., "
<< phix/CLHEP::deg << ", 90.,"
<< phiy/CLHEP::deg << ", 0, 0";
rotation = DDrot(DDName(rotstr, idNameSpace), theta, phix, theta, phiy,
0., 0.);
}
}
DDTranslation trdet(rposdet*cos(phi), rposdet*sin(phi), 0);
cpv.position(detIn, layerIn, n+1, trdet, rotation);
LogDebug("TIBGeom") << "DDTIBLayerAlgo test " << detIn.name()
<< " number " << n+1 << " positioned in "
<< layerIn.name() << " at " << trdet << " with "
<< rotation;
}
//Now the external layer
rin = cylinderInR + cylinderT;
rout = rmax-MFRingT;
name = idName + "Up";
solid = DDSolidFactory::tubs(DDName(name, idNameSpace), 0.5*layerL,
rin, rout, 0, CLHEP::twopi);
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: "
<< DDName(name, idNameSpace) << " Tubs made of "
<< genMat << " from 0 to " << CLHEP::twopi/CLHEP::deg
<< " with Rin " << rin << " Rout " << rout
<< " ZHalf " << 0.5*layerL;
DDLogicalPart layerOut(solid.ddname(), matter, solid);
cpv.position(layerOut, layer, 1, DDTranslation(0.0, 0.0, 0.0), DDRotation());
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: " << layerOut.name()
<< " number 1 positioned in " << layer.name()
<< " at (0,0,0) with no rotation";
rposdet = radiusUp;
dphi = CLHEP::twopi/stringsUp;
DDName detOut(DDSplit(detectorUp).first, DDSplit(detectorUp).second);
for (int n = 0; n < stringsUp; n++) {
double phi = (n+0.5)*dphi;
double phix = phi - detectorTilt - 90*CLHEP::deg;
double phideg = phix/CLHEP::deg;
DDRotation rotation;
if (phideg != 0) {
double theta = 90*CLHEP::deg;
double phiy = phix + 90.*CLHEP::deg;
std::string rotstr = idName + dbl_to_string(phideg*10.);
rotation = DDRotation(DDName(rotstr, idNameSpace));
if (!rotation) {
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: Creating a new "
<< "rotation: " << rotstr << "\t90., "
<< phix/CLHEP::deg << ", 90.,"
<< phiy/CLHEP::deg << ", 0, 0";
rotation = DDrot(DDName(rotstr, idNameSpace), theta, phix, theta, phiy,
0., 0.);
}
}
DDTranslation trdet(rposdet*cos(phi), rposdet*sin(phi), 0);
cpv.position(detOut, layerOut, n+1, trdet, rotation);
LogDebug("TIBGeom") << "DDTIBLayerAlgo test " << detOut.name()
<< " number " << n+1 << " positioned in "
<< layerOut.name() << " at " << trdet << " with "
<< rotation;
}
//
// Inner cylinder, support wall and ribs
//
// External skins
rin = cylinderInR;
rout = cylinderInR+cylinderT;
name = idName + "Cylinder";
solid = DDSolidFactory::tubs(DDName(name, idNameSpace), 0.5*layerL,
rin, rout, 0, CLHEP::twopi);
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: "
<< DDName(name, idNameSpace) << " Tubs made of "
<< cylinderMat << " from 0 to "
<< CLHEP::twopi/CLHEP::deg << " with Rin " << rin
<< " Rout " << rout << " ZHalf " << 0.5*layerL;
matname = DDName(DDSplit(cylinderMat).first, DDSplit(cylinderMat).second);
DDMaterial matcyl(matname);
DDLogicalPart cylinder(solid.ddname(), matcyl, solid);
cpv.position(cylinder, layer, 1, DDTranslation(0.0, 0.0, 0.0), DDRotation());
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: " << cylinder.name()
<< " number 1 positioned in " << layer.name()
<< " at (0,0,0) with no rotation";
//
// inner part of the cylinder
//
rin += supportT;
rout -= supportT;
name = idName + "CylinderIn";
solid = DDSolidFactory::tubs(DDName(name, idNameSpace), 0.5*layerL,
rin, rout, 0, CLHEP::twopi);
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: "
<< DDName(name, idNameSpace) << " Tubs made of "
<< genMat << " from 0 to " << CLHEP::twopi/CLHEP::deg
<< " with Rin " << rin << " Rout " << rout
<< " ZHalf " << 0.5*layerL;
DDLogicalPart cylinderIn(solid.ddname(), matter, solid);
cpv.position(cylinderIn, cylinder, 1, DDTranslation(0.0, 0.0, 0.0), DDRotation());
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: " << cylinderIn.name()
<< " number 1 positioned in " << cylinder.name()
<< " at (0,0,0) with no rotation";
//
// Filler Rings
//
matname = DDName(DDSplit(fillerMat).first, DDSplit(fillerMat).second);
DDMaterial matfiller(matname);
name = idName + "Filler";
solid = DDSolidFactory::tubs(DDName(name, idNameSpace), fillerDz, rin, rout,
0., CLHEP::twopi);
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: "
<< DDName(name, idNameSpace) << " Tubs made of "
<< fillerMat << " from " << 0. << " to "
<< CLHEP::twopi/CLHEP::deg << " with Rin " << rin
<< " Rout " << rout << " ZHalf " << fillerDz;
DDLogicalPart cylinderFiller(solid.ddname(), matfiller, solid);
cpv.position(cylinderFiller, cylinderIn, 1, DDTranslation(0.0, 0.0, 0.5*layerL-fillerDz), DDRotation());
cpv.position(cylinderFiller, cylinderIn, 2, DDTranslation(0.0, 0.0,-0.5*layerL+fillerDz), DDRotation());
LogDebug("TIBGeom") << "DDTIBLayerAlgo test " << cylinderFiller.name()
<< " number 1" << " positioned in "
<< cylinderIn.name() << " at " << DDTranslation(0.0, 0.0, 0.5*layerL-fillerDz)
<< " number 2" << " positioned in "
<< cylinderIn.name() << " at " << DDTranslation(0.0, 0.0,-0.5*layerL+fillerDz);
//
// Ribs
//
matname = DDName(DDSplit(ribMat).first, DDSplit(ribMat).second);
DDMaterial matrib(matname);
for (int i = 0; i < (int)(ribW.size()); i++) {
name = idName + "Rib" + dbl_to_string(i);
double width = 2.*ribW[i]/(rin+rout);
double dz = 0.5*layerL-2.*fillerDz;
solid = DDSolidFactory::tubs(DDName(name, idNameSpace), dz,
rin+0.5*CLHEP::mm, rout-0.5*CLHEP::mm,
-0.5*width, width);
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: "
<< DDName(name, idNameSpace) << " Tubs made of "
<< ribMat << " from " << -0.5*width/CLHEP::deg <<" to "
<< 0.5*width/CLHEP::deg << " with Rin "
<< rin+0.5*CLHEP::mm << " Rout "
<< rout-0.5*CLHEP::mm << " ZHalf " << dz;
DDLogicalPart cylinderRib(solid.ddname(), matrib, solid);
double phix = ribPhi[i];
double phideg = phix/CLHEP::deg;
DDRotation rotation;
if (phideg != 0) {
double theta = 90*CLHEP::deg;
double phiy = phix + 90.*CLHEP::deg;
std::string rotstr = idName + dbl_to_string(phideg*10.);
rotation = DDRotation(DDName(rotstr, idNameSpace));
if (!rotation) {
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: Creating a new "
<< "rotation: " << rotstr << "\t90., "
<< phix/CLHEP::deg << ", 90.," << phiy/CLHEP::deg
<< ", 0, 0";
rotation = DDrot(DDName(rotstr, idNameSpace), theta, phix, theta, phiy,
0., 0.);
}
}
DDTranslation tran(0, 0, 0);
cpv.position(cylinderRib, cylinderIn, 1, tran, rotation);
LogDebug("TIBGeom") << "DDTIBLayerAlgo test " << cylinderRib.name()
<< " number 1" << " positioned in "
<< cylinderIn.name() << " at " << tran << " with "
<< rotation;
}
//Manifold rings
//
// Inner ones first
matname = DDName(DDSplit(MFIntRingMat).first, DDSplit(MFIntRingMat).second);
DDMaterial matintmfr(matname);
rin = MFRingInR;
rout = rin + MFRingT;
name = idName + "InnerMFRing";
solid = DDSolidFactory::tubs(DDName(name, idNameSpace), MFRingDz,
rin, rout, 0, CLHEP::twopi);
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: "
<< DDName(name, idNameSpace) << " Tubs made of "
<< MFIntRingMat << " from 0 to "
<< CLHEP::twopi/CLHEP::deg << " with Rin " << rin
<< " Rout " << rout << " ZHalf " << MFRingDz;
DDLogicalPart inmfr(solid.ddname(), matintmfr, solid);
cpv.position(inmfr, layer, 1, DDTranslation(0.0, 0.0, -0.5*layerL+MFRingDz), DDRotation());
cpv.position(inmfr, layer, 2, DDTranslation(0.0, 0.0, +0.5*layerL-MFRingDz), DDRotation());
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: " << inmfr.name()
<< " number 1 and 2 positioned in " << layer.name()
<< " at (0,0,+-" << 0.5*layerL-MFRingDz << ") with no rotation";
// Outer ones
matname = DDName(DDSplit(MFExtRingMat).first, DDSplit(MFExtRingMat).second);
DDMaterial matextmfr(matname);
rout = MFRingOutR;
rin = rout - MFRingT;
name = idName + "OuterMFRing";
solid = DDSolidFactory::tubs(DDName(name, idNameSpace), MFRingDz,
rin, rout, 0, CLHEP::twopi);
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: "
<< DDName(name, idNameSpace) << " Tubs made of "
<< MFExtRingMat << " from 0 to "
<< CLHEP::twopi/CLHEP::deg << " with Rin " << rin
<< " Rout " << rout << " ZHalf " << MFRingDz;
DDLogicalPart outmfr(solid.ddname(), matextmfr, solid);
cpv.position(outmfr, layer, 1, DDTranslation(0.0, 0.0, -0.5*layerL+MFRingDz), DDRotation());
cpv.position(outmfr, layer, 2, DDTranslation(0.0, 0.0, +0.5*layerL-MFRingDz), DDRotation());
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: " << outmfr.name()
<< " number 1 and 2 positioned in " << layer.name()
<< " at (0,0,+-" << 0.5*layerL-MFRingDz
<< ") with no rotation";
//Central Support rings
//
matname = DDName(DDSplit(centMat).first, DDSplit(centMat).second);
DDMaterial matcent(matname);
// Ring 1
double centZ = centRing1par[0];
double centDz = 0.5*centRing1par[1];
rin = centRing1par[2];
rout = centRing1par[3];
name = idName + "CentRing1";
solid = DDSolidFactory::tubs(DDName(name, idNameSpace), centDz,
rin, rout, 0, CLHEP::twopi);
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: "
<< DDName(name, idNameSpace) << " Tubs made of "
<< centMat << " from 0 to " << CLHEP::twopi/CLHEP::deg
<< " with Rin " << rin << " Rout " << rout
<< " ZHalf " << centDz;
DDLogicalPart cent1(solid.ddname(), matcent, solid);
cpv.position(cent1, layer, 1, DDTranslation(0.0, 0.0, centZ), DDRotation());
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: " << cent1.name()
<< " positioned in " << layer.name()
<< " at (0,0," << centZ << ") with no rotation";
// Ring 2
centZ = centRing2par[0];
centDz = 0.5*centRing2par[1];
rin = centRing2par[2];
rout = centRing2par[3];
name = idName + "CentRing2";
solid = DDSolidFactory::tubs(DDName(name, idNameSpace), centDz,
rin, rout, 0, CLHEP::twopi);
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: "
<< DDName(name, idNameSpace) << " Tubs made of "
<< centMat << " from 0 to " << CLHEP::twopi/CLHEP::deg
<< " with Rin " << rin << " Rout " << rout
<< " ZHalf " << centDz;
DDLogicalPart cent2(solid.ddname(), matcent, solid);
cpv.position(cent2, layer, 1, DDTranslation(0.0, 0.0, centZ), DDRotation());
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: " << cent2.name()
<< " positioned in " << layer.name()
<< " at (0,0," << centZ << ") with no rotation";
//
// Preparing DOHM Carrier solid
name = idName + "DOHMCarrier";
double dohmCarrierRin = MFRingOutR - MFRingT;
double dohmCarrierRout = MFRingOutR;
double dohmCarrierDz = 0.5*(dohmPrimL+dohmtoMF);
double dohmCarrierZ = 0.5*layerL-2.*MFRingDz-dohmCarrierDz;
solid = DDSolidFactory::tubs(DDName(name, idNameSpace), dohmCarrierDz,
dohmCarrierRin, dohmCarrierRout,
dohmCarrierPhiOff,
180.*CLHEP::deg-2.*dohmCarrierPhiOff);
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: "
<< DDName(name, idNameSpace) << " Tubs made of "
<< dohmCarrierMaterial << " from "
<< dohmCarrierPhiOff << " to "
<< 180.*CLHEP::deg-dohmCarrierPhiOff << " with Rin "
<< dohmCarrierRin << " Rout " << MFRingOutR << " ZHalf "
<< dohmCarrierDz;
// Define FW and BW carrier logical volume and
// place DOHM Primary and auxiliary modules inside it
dphi = CLHEP::twopi/stringsUp;
DDRotation dohmRotation;
double dohmR = 0.5*(dohmCarrierRin+dohmCarrierRout);
for (int j = 0; j<4; j++) {
std::vector<double> dohmList;
DDTranslation tran;
std::string rotstr;
DDRotation rotation;
int dohmCarrierReplica=0;
int placeDohm = 0;
switch (j){
case 0:
name = idName + "DOHMCarrierFW";
dohmList = dohmListFW;
tran = DDTranslation(0., 0., dohmCarrierZ);
rotstr = idName + "FwUp";
rotation = DDRotation();
dohmCarrierReplica = 1;
placeDohm=1;
break;
case 1:
name = idName + "DOHMCarrierFW";
dohmList = dohmListFW;
tran = DDTranslation(0., 0., dohmCarrierZ);
rotstr = idName + "FwDown";
rotation = DDrot(DDName(rotstr, idNameSpace), 90.*CLHEP::deg,
180.*CLHEP::deg, 90.*CLHEP::deg,270.*CLHEP::deg, 0.,0.);
dohmCarrierReplica = 2;
placeDohm=0;
break;
case 2:
name = idName + "DOHMCarrierBW";
dohmList = dohmListBW;
tran = DDTranslation(0., 0., -dohmCarrierZ);
rotstr = idName + "BwUp";
rotation = DDrot(DDName(rotstr, idNameSpace), 90.*CLHEP::deg,
180.*CLHEP::deg, 90.*CLHEP::deg, 90.*CLHEP::deg,
180.*CLHEP::deg, 0.);
dohmCarrierReplica = 1;
placeDohm=1;
break;
case 3:
name = idName + "DOHMCarrierBW";
dohmList = dohmListBW;
tran = DDTranslation(0., 0., -dohmCarrierZ);
rotstr = idName + "BwDown";
rotation = DDrot(DDName(rotstr, idNameSpace), 90.*CLHEP::deg, 0.,
90.*CLHEP::deg, 270.*CLHEP::deg, 180.*CLHEP::deg, 0.);
dohmCarrierReplica = 2;
placeDohm=0;
break;
}
DDLogicalPart dohmCarrier(name,DDMaterial(dohmCarrierMaterial),solid);
int primReplica = 0;
int auxReplica = 0;
for (int i = 0; i < placeDohm*((int)(dohmList.size())); i++) {
double phi = (std::abs(dohmList[i])+0.5-1.)*dphi;
double phix = phi + 90*CLHEP::deg;
double phideg = phix/CLHEP::deg;
if (phideg != 0) {
double theta = 90*CLHEP::deg;
double phiy = phix + 90.*CLHEP::deg;
std::string rotstr = idName+dbl_to_string(std::abs(dohmList[i])-1.);
dohmRotation = DDRotation(DDName(rotstr, idNameSpace));
if (!dohmRotation) {
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: Creating a new "
<< "rotation: " << rotstr << "\t" << theta
<< ", " << phix/CLHEP::deg << ", " << theta
<< ", " << phiy/CLHEP::deg <<", 0, 0";
dohmRotation = DDrot(DDName(rotstr, idNameSpace), theta, phix, theta,
phiy, 0., 0.);
}
}
std::string dohmName;
int dohmReplica = 0;
double dohmZ = 0.;
if(dohmList[i]<0.) {
// Place a Auxiliary DOHM
dohmName = dohmAuxName;
dohmZ = dohmCarrierDz - 0.5*dohmAuxL - dohmtoMF;
primReplica++;
dohmReplica = primReplica;
} else {
// Place a Primary DOHM
dohmName = dohmPrimName;
dohmZ = dohmCarrierDz - 0.5*dohmPrimL - dohmtoMF;
auxReplica++;
dohmReplica = auxReplica;
}
DDName dohm(DDSplit(dohmName).first, DDSplit(dohmName).second);
DDTranslation dohmTrasl(dohmR*cos(phi), dohmR*sin(phi), dohmZ);
cpv.position(dohm, dohmCarrier, dohmReplica, dohmTrasl, dohmRotation);
LogDebug("TIBGeom") << "DDTIBLayerAlgo test " << dohm.name()
<< " replica " << dohmReplica << " positioned in "
<< dohmCarrier.name() << " at " << dohmTrasl << " with "
<< dohmRotation;
}
cpv.position(dohmCarrier, parent(), dohmCarrierReplica, tran, rotation );
LogDebug("TIBGeom") << "DDTIBLayerAlgo test "
<< dohmCarrier.name() << " positioned in " << parent().name() << " at "
<< tran << " with " << rotation;
}
for (int j = 0; j<4; j++) {
matname = DDName(DDSplit(pillarMaterial).first, DDSplit(pillarMaterial).second);
DDMaterial pillarMat(matname);
std::vector<double> pillarZ;
std::vector<double> pillarPhi;
double pillarDz=0, pillarDPhi=0, pillarRin=0, pillarRout=0;
switch (j){
case 0:
name = idName + "FWIntPillar";
pillarZ = fwIntPillarZ;
pillarPhi = fwIntPillarPhi;
pillarRin = MFRingInR;
pillarRout = MFRingInR + MFRingT;
pillarDz = fwIntPillarDz;
pillarDPhi = fwIntPillarDPhi;
break;
case 1:
name = idName + "BWIntPillar";
pillarZ = bwIntPillarZ;
pillarPhi = bwIntPillarPhi;
pillarRin = MFRingInR;
pillarRout = MFRingInR + MFRingT;
pillarDz = bwIntPillarDz;
pillarDPhi = bwIntPillarDPhi;
break;
case 2:
name = idName + "FWExtPillar";
pillarZ = fwExtPillarZ;
pillarPhi = fwExtPillarPhi;
pillarRin = MFRingOutR - MFRingT;
pillarRout = MFRingOutR;
pillarDz = fwExtPillarDz;
pillarDPhi = fwExtPillarDPhi;
break;
case 3:
name = idName + "BWExtPillar";
pillarZ = bwExtPillarZ;
pillarPhi = bwExtPillarPhi;
pillarRin = MFRingOutR - MFRingT;
pillarRout = MFRingOutR;
pillarDz = bwExtPillarDz;
pillarDPhi = bwExtPillarDPhi;
break;
}
solid = DDSolidFactory::tubs(DDName(name, idNameSpace), pillarDz,
pillarRin, pillarRout,
-pillarDPhi, 2.*pillarDPhi);
DDLogicalPart Pillar(name,DDMaterial(pillarMat),solid);
LogDebug("TIBGeom") << "DDTIBLayerAlgo test: "
<< DDName(name, idNameSpace) << " Tubs made of "
<< pillarMat << " from "
<< -pillarDPhi << " to "
<< pillarDPhi << " with Rin "
<< pillarRin << " Rout " << pillarRout << " ZHalf "
<< pillarDz;
DDTranslation pillarTran;
DDRotation pillarRota;
int pillarReplica = 0;
for (unsigned int i=0; i<pillarZ.size(); i++) {
if( pillarPhi[i]>0. ) {
pillarTran = DDTranslation(0., 0., pillarZ[i]);
pillarRota = DDanonymousRot(DDcreateRotationMatrix(90.*CLHEP::deg, pillarPhi[i], 90.*CLHEP::deg, 90.*CLHEP::deg+pillarPhi[i], 0., 0.));
cpv.position(Pillar, parent(), i, pillarTran, pillarRota);
LogDebug("TIBGeom") << "DDTIBLayerAlgo test "
<< Pillar.name() << " positioned in "
<< parent().name() << " at "
<< pillarTran << " with " << pillarRota
<< " copy number " << pillarReplica;
pillarReplica++;
}
}
}
}
| void DDTIBLayerAlgo::initialize | ( | const DDNumericArguments & | nArgs, |
| const DDVectorArguments & | vArgs, | ||
| const DDMapArguments & | mArgs, | ||
| const DDStringArguments & | sArgs, | ||
| const DDStringVectorArguments & | vsArgs | ||
| ) |
Definition at line 27 of file DDTIBLayerAlgo.cc.
References bwExtPillarDPhi, bwExtPillarDz, bwExtPillarPhi, bwExtPillarZ, bwIntPillarDPhi, bwIntPillarDz, bwIntPillarPhi, bwIntPillarZ, centMat, centRing1par, centRing2par, cylinderInR, cylinderMat, cylinderT, detectorLo, detectorTilt, detectorUp, dohmAuxL, dohmAuxMaterial, dohmAuxName, dohmCableMaterial, dohmCarrierMaterial, dohmCarrierPhiOff, dohmListBW, dohmListFW, dohmPrimL, dohmPrimMaterial, dohmPrimName, dohmtoMF, fillerDz, fillerMat, fwExtPillarDPhi, fwExtPillarDz, fwExtPillarPhi, fwExtPillarZ, fwIntPillarDPhi, fwIntPillarDz, fwIntPillarPhi, fwIntPillarZ, genMat, i, idNameSpace, layerL, LogDebug, MFExtRingMat, MFIntRingMat, MFRingDz, MFRingInR, MFRingOutR, MFRingT, DDCurrentNamespace::ns(), dbtoconf::parent, pillarMaterial, radiusLo, radiusUp, ribMat, ribPhi, ribW, stringsLo, stringsUp, and supportT.
{
idNameSpace = DDCurrentNamespace::ns();
genMat = sArgs["GeneralMaterial"];
DDName parentName = parent().name();
LogDebug("TIBGeom") << "DDTIBLayerAlgo debug: Parent " << parentName
<< " NameSpace " << idNameSpace
<< " General Material " << genMat;
detectorTilt = nArgs["DetectorTilt"];
layerL = nArgs["LayerL"];
radiusLo = nArgs["RadiusLo"];
stringsLo = int(nArgs["StringsLo"]);
detectorLo = sArgs["StringDetLoName"];
LogDebug("TIBGeom") << "DDTIBLayerAlgo debug: Lower layer Radius "
<< radiusLo << " Number " << stringsLo << " String "
<< detectorLo;
radiusUp = nArgs["RadiusUp"];
stringsUp = int(nArgs["StringsUp"]);
detectorUp = sArgs["StringDetUpName"];
LogDebug("TIBGeom") << "DDTIBLayerAlgo debug: Upper layer Radius "
<< radiusUp << " Number " << stringsUp << " String "
<< detectorUp;
cylinderT = nArgs["CylinderThickness"];
cylinderInR = nArgs["CylinderInnerRadius"];
cylinderMat = sArgs["CylinderMaterial"];
MFRingInR = nArgs["MFRingInnerRadius"];
MFRingOutR = nArgs["MFRingOuterRadius"];
MFRingT = nArgs["MFRingThickness"];
MFRingDz = nArgs["MFRingDeltaz"];
MFIntRingMat = sArgs["MFIntRingMaterial"];
MFExtRingMat = sArgs["MFExtRingMaterial"];
supportT = nArgs["SupportThickness"];
centMat = sArgs["CentRingMaterial"];
centRing1par = vArgs["CentRing1"];
centRing2par = vArgs["CentRing2"];
fillerMat = sArgs["FillerMaterial"];
fillerDz = nArgs["FillerDeltaz"];
ribMat = sArgs["RibMaterial"];
ribW = vArgs["RibWidth"];
ribPhi = vArgs["RibPhi"];
LogDebug("TIBGeom") << "DDTIBLayerAlgo debug: Cylinder Material/"
<< "thickness " << cylinderMat << " " << cylinderT
<< " Rib Material " << ribMat << " at "
<< ribW.size() << " positions with width/phi";
for (unsigned int i = 0; i < ribW.size(); i++)
LogDebug("TIBGeom") << "\tribW[" << i << "] = " << ribW[i]
<< "\tribPhi[" << i << "] = " << ribPhi[i]/CLHEP::deg;
dohmCarrierPhiOff = nArgs["DOHMCarrierPhiOffset"];
dohmtoMF = nArgs["DOHMtoMFDist"];
dohmPrimName = sArgs["StringDOHMPrimName"];
dohmAuxName = sArgs["StringDOHMAuxName"];
dohmCarrierMaterial = sArgs["DOHMCarrierMaterial"];
dohmCableMaterial = sArgs["DOHMCableMaterial"];
dohmPrimL = nArgs["DOHMPRIMLength"];
dohmPrimMaterial = sArgs["DOHMPRIMMaterial"];
dohmAuxL = nArgs["DOHMAUXLength"];
dohmAuxMaterial = sArgs["DOHMAUXMaterial"];
dohmListFW = vArgs["DOHMListFW"];
dohmListBW = vArgs["DOHMListBW"];
LogDebug("TIBGeom") << "DDTIBLayerAlgo debug: DOHM Primary "
<< " Material " << dohmPrimMaterial << " Length " << dohmPrimL;
LogDebug("TIBGeom") << "DDTIBLayerAlgo debug: DOHM Aux "
<< " Material " << dohmAuxMaterial << " Length " << dohmAuxL;
for (unsigned int i=0; i<dohmListFW.size(); i++) {
if (dohmListFW[i]>0.) LogDebug("TIBGeom") << "DOHM Primary at FW Position " << dohmListFW[i];
if (dohmListFW[i]<0.) LogDebug("TIBGeom") << "DOHM Aux at FW Position " << -dohmListFW[i];
}
for (unsigned int i=0; i<dohmListBW.size(); i++) {
if (dohmListBW[i]>0.) LogDebug("TIBGeom") << "DOHM Primary at BW Position " << dohmListBW[i];
if (dohmListBW[i]<0.) LogDebug("TIBGeom") << "DOHM Aux at BW Position " << -dohmListBW[i];
}
//Pillar Material
pillarMaterial = sArgs["PillarMaterial"];
// Internal Pillar Parameters
fwIntPillarDz = nArgs["FWIntPillarDz"];
fwIntPillarDPhi = nArgs["FWIntPillarDPhi"];
fwIntPillarZ = vArgs["FWIntPillarZ"];
fwIntPillarPhi = vArgs["FWIntPillarPhi"];
bwIntPillarDz = nArgs["BWIntPillarDz"];
bwIntPillarDPhi = nArgs["BWIntPillarDPhi"];
bwIntPillarZ = vArgs["BWIntPillarZ"];
bwIntPillarPhi = vArgs["BWIntPillarPhi"];
LogDebug("TIBGeom") << "FW Internal Pillar [Dz, DPhi] "
<< fwIntPillarDz << ", "
<< fwIntPillarDPhi;
for (unsigned int i=0; i<fwIntPillarZ.size(); i++) {
if( fwIntPillarPhi[i]>0. ) {
LogDebug("TIBGeom") << " at positions [z, phi] "
<< fwIntPillarZ[i] << " " << fwIntPillarPhi[i];
}
}
LogDebug("TIBGeom") << "BW Internal Pillar [Dz, DPhi] "
<< bwIntPillarDz << ", "
<< bwIntPillarDPhi;
for (unsigned int i=0; i<bwIntPillarZ.size(); i++) {
if( bwIntPillarPhi[i]>0. ) {
LogDebug("TIBGeom") << " at positions [z, phi] "
<< bwIntPillarZ[i] << " " << bwIntPillarPhi[i];
}
}
// External Pillar Parameters
fwExtPillarDz = nArgs["FWExtPillarDz"];
fwExtPillarDPhi = nArgs["FWExtPillarDPhi"];
fwExtPillarZ = vArgs["FWExtPillarZ"];
fwExtPillarPhi = vArgs["FWExtPillarPhi"];
bwExtPillarDz = nArgs["BWExtPillarDz"];
bwExtPillarDPhi = nArgs["BWExtPillarDPhi"];
bwExtPillarZ = vArgs["BWExtPillarZ"];
bwExtPillarPhi = vArgs["BWExtPillarPhi"];
LogDebug("TIBGeom") << "FW External Pillar [Dz, DPhi] "
<< fwExtPillarDz << ", "
<< fwExtPillarDPhi;
for (unsigned int i=0; i<fwExtPillarZ.size(); i++) {
if( fwExtPillarPhi[i]>0. ) {
LogDebug("TIBGeom") << " at positions [z, phi] "
<< fwExtPillarZ[i] << " " << fwExtPillarPhi[i];
}
}
LogDebug("TIBGeom") << "BW External Pillar [Dz, DPhi] "
<< bwExtPillarDz << ", "
<< bwExtPillarDPhi;
for (unsigned int i=0; i<bwExtPillarZ.size(); i++) {
if( bwExtPillarPhi[i]>0. ) {
LogDebug("TIBGeom") << " at positions [z, phi] "
<< bwExtPillarZ[i] << " " << bwExtPillarPhi[i];
}
}
}
double DDTIBLayerAlgo::bwExtPillarDPhi [private] |
Definition at line 95 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::bwExtPillarDz [private] |
Definition at line 94 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::vector<double> DDTIBLayerAlgo::bwExtPillarPhi [private] |
Definition at line 97 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::vector<double> DDTIBLayerAlgo::bwExtPillarZ [private] |
Definition at line 96 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::bwIntPillarDPhi [private] |
Definition at line 86 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::bwIntPillarDz [private] |
Definition at line 85 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::vector<double> DDTIBLayerAlgo::bwIntPillarPhi [private] |
Definition at line 88 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::vector<double> DDTIBLayerAlgo::bwIntPillarZ [private] |
Definition at line 87 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::string DDTIBLayerAlgo::centMat [private] |
Definition at line 53 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::vector<double> DDTIBLayerAlgo::centRing1par [private] |
Definition at line 54 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::vector<double> DDTIBLayerAlgo::centRing2par [private] |
Definition at line 55 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::cylinderInR [private] |
Definition at line 42 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::string DDTIBLayerAlgo::cylinderMat [private] |
Definition at line 43 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::cylinderT [private] |
Definition at line 41 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::string DDTIBLayerAlgo::detectorLo [private] |
Definition at line 34 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::detectorTilt [private] |
Definition at line 28 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::string DDTIBLayerAlgo::detectorUp [private] |
Definition at line 39 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::dohmAuxL [private] |
Definition at line 76 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::string DDTIBLayerAlgo::dohmAuxMaterial [private] |
Definition at line 77 of file DDTIBLayerAlgo.h.
Referenced by initialize().
std::string DDTIBLayerAlgo::dohmAuxName [private] |
Definition at line 70 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::string DDTIBLayerAlgo::dohmCableMaterial [private] |
Definition at line 73 of file DDTIBLayerAlgo.h.
Referenced by initialize().
std::string DDTIBLayerAlgo::dohmCarrierMaterial [private] |
Definition at line 72 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::dohmCarrierPhiOff [private] |
Definition at line 68 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::vector<double> DDTIBLayerAlgo::dohmListBW [private] |
Definition at line 65 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::vector<double> DDTIBLayerAlgo::dohmListFW [private] |
Definition at line 64 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::dohmPrimL [private] |
Definition at line 74 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::string DDTIBLayerAlgo::dohmPrimMaterial [private] |
Definition at line 75 of file DDTIBLayerAlgo.h.
Referenced by initialize().
std::string DDTIBLayerAlgo::dohmPrimName [private] |
Definition at line 69 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::dohmtoMF [private] |
Definition at line 67 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::fillerDz [private] |
Definition at line 58 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::string DDTIBLayerAlgo::fillerMat [private] |
Definition at line 57 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::fwExtPillarDPhi [private] |
Definition at line 91 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::fwExtPillarDz [private] |
Definition at line 90 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::vector<double> DDTIBLayerAlgo::fwExtPillarPhi [private] |
Definition at line 93 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::vector<double> DDTIBLayerAlgo::fwExtPillarZ [private] |
Definition at line 92 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::fwIntPillarDPhi [private] |
Definition at line 82 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::fwIntPillarDz [private] |
Definition at line 81 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::vector<double> DDTIBLayerAlgo::fwIntPillarPhi [private] |
Definition at line 84 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::vector<double> DDTIBLayerAlgo::fwIntPillarZ [private] |
Definition at line 83 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::string DDTIBLayerAlgo::genMat [private] |
Definition at line 27 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::string DDTIBLayerAlgo::idNameSpace [private] |
Definition at line 26 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::layerL [private] |
Definition at line 29 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::string DDTIBLayerAlgo::MFExtRingMat [private] |
Definition at line 49 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::string DDTIBLayerAlgo::MFIntRingMat [private] |
Definition at line 48 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::MFRingDz [private] |
Definition at line 47 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::MFRingInR [private] |
Definition at line 44 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::MFRingOutR [private] |
Definition at line 45 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::MFRingT [private] |
Definition at line 46 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::phioffLo [private] |
Definition at line 32 of file DDTIBLayerAlgo.h.
double DDTIBLayerAlgo::phioffUp [private] |
Definition at line 37 of file DDTIBLayerAlgo.h.
std::string DDTIBLayerAlgo::pillarMaterial [private] |
Definition at line 79 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::radiusLo [private] |
Definition at line 31 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::radiusUp [private] |
Definition at line 36 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::string DDTIBLayerAlgo::ribMat [private] |
Definition at line 60 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::vector<double> DDTIBLayerAlgo::ribPhi [private] |
Definition at line 62 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
std::vector<double> DDTIBLayerAlgo::ribW [private] |
Definition at line 61 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
int DDTIBLayerAlgo::stringsLo [private] |
Definition at line 33 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
int DDTIBLayerAlgo::stringsUp [private] |
Definition at line 38 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
double DDTIBLayerAlgo::supportT [private] |
Definition at line 51 of file DDTIBLayerAlgo.h.
Referenced by execute(), and initialize().
1.7.3