Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-09-01 02:45:08

0001 #include "DD4hep/DetFactoryHelper.h"
0002 #include "DataFormats/Math/interface/CMSUnits.h"
0003 #include "DetectorDescription/DDCMS/interface/DDPlugins.h"
0004 #include "FWCore/MessageLogger/interface/MessageLogger.h"
0005 
0006 using namespace cms_units::operators;
0007 
0008 static long algorithm(dd4hep::Detector&, cms::DDParsingContext& ctxt, xml_h e) {
0009   cms::DDNamespace ns(ctxt, e, true);
0010   cms::DDAlgoArguments args(ctxt, e);
0011   std::string parentName = args.parentName();
0012 
0013   std::string genMat = args.value<std::string>("GeneralMaterial");
0014   int number = args.value<int>("Ladders");
0015   double layerDz = args.value<double>("LayerDz");
0016   double coolDz = args.value<double>("CoolDz");
0017   double coolThick = args.value<double>("CoolThick");
0018   double coolRadius = args.value<double>("CoolRadius");
0019   double coolDist = args.value<double>("CoolDist");
0020   double cool1Offset = args.value<double>("Cool1Offset");
0021   double cool2Offset = args.value<double>("Cool2Offset");
0022   std::string coolMat = args.value<std::string>("CoolMaterial");
0023   std::string tubeMat = args.value<std::string>("CoolTubeMaterial");
0024   std::string coolMatHalf = args.value<std::string>("CoolMaterialHalf");
0025   std::string tubeMatHalf = args.value<std::string>("CoolTubeMaterialHalf");
0026   double phiFineTune = args.value<double>("PitchFineTune");
0027   double rOuterFineTune = args.value<double>("OuterOffsetFineTune");
0028   double rInnerFineTune = args.value<double>("InnerOffsetFineTune");
0029 
0030   // FixMe : Would need ns.vecStr
0031   std::string ladder = args.value<std::string>("LadderName");
0032 
0033   double ladderWidth = args.value<double>("LadderWidth");
0034   double ladderThick = args.value<double>("LadderThick");
0035   double ladderOffset = args.value<double>("LadderOffset");
0036   int outerFirst = args.value<int>("OuterFirst");
0037 
0038   edm::LogVerbatim("PixelGeom") << "DDPixBarLayerUpgradeAlgo debug: Parent " << parentName << " NameSpace "
0039                                 << ns.objName(parentName) << "\n\tLadders " << number << "\tGeneral Material " << genMat
0040                                 << "\tLength " << layerDz << "\tSpecification of Cooling Pieces:\n"
0041                                 << "\tLength " << coolDz << "\tThickness of Shell " << coolThick << "\tRadial distance "
0042                                 << coolDist << "\tMaterials " << coolMat << ", " << tubeMat << std::endl;
0043 
0044   edm::LogVerbatim("PixelGeom") << "DDPixBarLayerUpgradeAlgo debug: Full Ladder " << ladder << " width/thickness "
0045                                 << ladderWidth << ", " << ladderThick;
0046 
0047   double dphi = 2_pi / (double)number;
0048   double x2 = coolDist * sin(0.5 * dphi);
0049   double rtmi = coolDist * cos(0.5 * dphi) - (coolRadius + ladderThick) + rInnerFineTune;
0050   double rmxh = coolDist * cos(0.5 * dphi) + (coolRadius + ladderThick + ladderOffset) + rOuterFineTune;
0051   double rtmx = sqrt(rmxh * rmxh + ladderWidth * ladderWidth / 4.);
0052 
0053   std::string name = ns.objName(parentName);
0054 
0055   dd4hep::Solid solid = ns.addSolid(ns.prepend(name), dd4hep::Tube(rtmi, rtmx, 0.5 * layerDz, 0., 2._pi));
0056 
0057   edm::LogVerbatim("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: " << solid.name() << " Tubs made of " << genMat
0058                                 << " from 0. to " << convertRadToDeg(2_pi) << " with Rin " << rtmi << " Rout " << rtmx
0059                                 << " ZHalf " << 0.5 * layerDz;
0060 
0061   dd4hep::Volume layer = ns.addVolume(dd4hep::Volume(solid.name(), solid, ns.material(genMat)));
0062 
0063   // Full Tubes
0064   solid = ns.addSolid(ns.prepend(name + "CoolTube"), dd4hep::Tube(0., coolRadius, 0.5 * coolDz, 0., 2_pi));
0065 
0066   edm::LogVerbatim("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: " << solid.name() << " Tubs made of " << tubeMat
0067                                 << " from 0 to " << convertRadToDeg(2_pi) << " with Rout " << coolRadius << " ZHalf "
0068                                 << 0.5 * coolDz;
0069 
0070   dd4hep::Volume coolTube = ns.addVolume(dd4hep::Volume(solid.name(), solid, ns.material(tubeMat)));
0071 
0072   // Half Tubes
0073   solid = ns.addSolid(ns.prepend(name + "CoolTubeHalf"), dd4hep::Tube(0., coolRadius, 0.5 * coolDz, 0, 1_pi));
0074 
0075   edm::LogVerbatim("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: " << solid.name() << " Tubs made of " << tubeMatHalf
0076                                 << " from 0 to " << convertRadToDeg(2_pi) << " with Rout " << coolRadius << " ZHalf "
0077                                 << 0.5 * coolDz;
0078 
0079   dd4hep::Volume coolTubeHalf = ns.addVolume(dd4hep::Volume(solid.name(), solid, ns.material(tubeMatHalf)));
0080 
0081   // Full Coolant
0082   name = ns.objName(parentName);
0083 
0084   solid = ns.addSolid(ns.prepend(name + "Coolant"), dd4hep::Tube(0., coolRadius - coolThick, 0.5 * coolDz, 0., 2_pi));
0085 
0086   edm::LogVerbatim("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: " << solid.name() << " Tubs made of " << tubeMat
0087                                 << " from 0 to " << convertRadToDeg(2._pi) << " with Rout " << coolRadius - coolThick
0088                                 << " ZHalf " << 0.5 * coolDz;
0089 
0090   dd4hep::Volume cool = ns.addVolume(dd4hep::Volume(solid.name(), solid, ns.material(coolMat)));
0091 
0092   coolTube.placeVolume(cool, 1, dd4hep::Transform3D(dd4hep::Rotation3D(), dd4hep::Position(0., 0., 0.)));
0093 
0094   edm::LogVerbatim("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: " << cool.name() << " number 1 positioned in "
0095                                 << coolTube.name() << " at (0,0,0) with no rotation";
0096 
0097   // Half Coolant
0098 
0099   solid =
0100       ns.addSolid(ns.prepend(name + "CoolantHalf"), dd4hep::Tube(0., coolRadius - coolThick, 0.5 * coolDz, 0., 1._pi));
0101 
0102   edm::LogVerbatim("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: " << solid.name() << " Tubs made of " << tubeMatHalf
0103                                 << " from 0 to " << convertRadToDeg(2._pi) << " with Rout " << coolRadius - coolThick
0104                                 << " ZHalf " << 0.5 * coolDz;
0105 
0106   dd4hep::Volume coolHalf = ns.addVolume(dd4hep::Volume(solid.name(), solid, ns.material(coolMatHalf)));
0107 
0108   coolTubeHalf.placeVolume(coolHalf, 1, dd4hep::Transform3D(dd4hep::Rotation3D(), dd4hep::Position(0., 0., 0.)));
0109 
0110   edm::LogVerbatim("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: " << cool.name() << " number 1 positioned in "
0111                                 << coolTube.name() << " at (0,0,0) with no rotation";
0112   int copy = 1, iup = (-1) * outerFirst;
0113   int copyoffset = number + 2;
0114 
0115   for (int i = 1; i < number + 1; i++) {
0116     double phi = i * dphi + 90._deg - 0.5 * dphi + phiFineTune;  //to start with the interface ladder
0117     double phix, phiy, rrr, rrroffset;
0118     std::string rots;
0119 
0120     auto tran = dd4hep::Position();
0121     auto rot = dd4hep::Rotation3D();
0122 
0123     iup = -iup;
0124 
0125     double dr;
0126 
0127     if ((i == 1) || (i == number / 2 + 1)) {
0128       dr = coolRadius + 0.5 * ladderThick + ladderOffset;  //interface ladder offset
0129     } else {
0130       dr = coolRadius + 0.5 * ladderThick;
0131     }
0132 
0133     if (i % 2 == 1) {
0134       rrr = coolDist * cos(0.5 * dphi) + iup * dr + rOuterFineTune;
0135     } else {
0136       rrr = coolDist * cos(0.5 * dphi) + iup * dr + rInnerFineTune;
0137     }
0138 
0139     tran = dd4hep::Position(rrr * cos(phi), rrr * sin(phi), 0);
0140     rots = name + std::to_string(copy);
0141 
0142     if (iup > 0) {
0143       phix = phi - 90._deg;
0144     } else {
0145       phix = phi + 90._deg;
0146     }
0147 
0148     phiy = phix + 90._deg;
0149 
0150     edm::LogVerbatim("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: Creating a new "
0151                                   << "rotation: " << rots << "\t90., " << convertRadToDeg(phix) << ", 90.,"
0152                                   << convertRadToDeg(phiy) << ", 0, 0";
0153 
0154     rot = cms::makeRotation3D(90._deg, phix, 90._deg, phiy, 0., 0.);
0155 
0156     // FixMe : Would need ns.vecStr
0157     layer.placeVolume(ns.volume(ladder), copy, dd4hep::Transform3D(rot, tran));
0158 
0159     edm::LogVerbatim("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: " << ladder << " number " << copy
0160                                   << " positioned in " << layer.name() << " at " << tran << " with " << rot;
0161 
0162     copy++;
0163 
0164     rrr = coolDist * cos(0.5 * dphi) + coolRadius / 2.;
0165 
0166     rots = name + std::to_string(i + 100);
0167     phix = phi + 90._deg;
0168 
0169     if (iup < 0)
0170       phix += dphi;
0171 
0172     phiy = phix + 90._deg;
0173 
0174     edm::LogVerbatim("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: Creating a new "
0175                                   << "rotation: " << rots << "\t90., " << convertRadToDeg(phix) << ", 90.,"
0176                                   << convertRadToDeg(phiy) << ", 0, 0";
0177 
0178     tran = dd4hep::Position(rrr * cos(phi) - x2 * sin(phi), rrr * sin(phi) + x2 * cos(phi), 0.);
0179 
0180     rot = cms::makeRotation3D(90._deg, phix, 90._deg, phiy, 0., 0.);
0181 
0182     layer.placeVolume(coolTubeHalf, i + 1, dd4hep::Transform3D(rot, tran));
0183 
0184     if ((i == 1) || (i == number / 2 + 1)) {
0185       rrroffset = coolDist * cos(0.5 * dphi) + iup * ladderOffset + rOuterFineTune;
0186       tran = dd4hep::Position(
0187           rrroffset * cos(phi) - cool1Offset * sin(phi), rrroffset * sin(phi) + cool1Offset * cos(phi), 0.);
0188 
0189       layer.placeVolume(coolTube, copyoffset, dd4hep::Transform3D(dd4hep::Rotation3D(), tran));
0190       ;
0191 
0192       copyoffset++;
0193       tran = dd4hep::Position(
0194           rrroffset * cos(phi) - cool2Offset * sin(phi), rrroffset * sin(phi) + cool2Offset * cos(phi), 0.);
0195 
0196       layer.placeVolume(coolTube, copyoffset, dd4hep::Transform3D(dd4hep::Rotation3D(), tran));
0197       copyoffset++;
0198     }
0199 
0200     edm::LogVerbatim("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: " << coolTube.name() << " number " << i + 1
0201                                   << " positioned in " << layer.name() << " at " << tran << " with " << rot;
0202   }
0203 
0204   return cms::s_executed;
0205 }
0206 
0207 DECLARE_DDCMS_DETELEMENT(DDCMS_track_DDPixBarLayerUpgradeAlgo, algorithm);