DDEcalPreshowerAlgo.cc

Go to the documentation of this file.

#include <cmath>
#include <algorithm>

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DetectorDescription/Core/interface/DDTypes.h"
#include "DetectorDescription/Core/interface/DDLogicalPart.h"
#include "DetectorDescription/Core/interface/DDSolid.h"
#include "CLHEP/Units/GlobalSystemOfUnits.h"
#include <vector>
#include <string>

#include "DetectorDescription/Core/interface/DDAlgorithm.h"
#include "DetectorDescription/Core/interface/DDMaterial.h"

//#define EDM_ML_DEBUG

class DDEcalPreshowerAlgo : public DDAlgorithm {
public:
  DDMaterial getMaterial(unsigned int i) const { return DDMaterial(materials_[i]); }
  DDMaterial getLaddMaterial() const { return DDMaterial(LaddMaterial_); }
  std::string getLayName(unsigned int i) const { return layName_[i]; }
  std::string getLadPrefix(unsigned int i) const { return ladPfx_[i]; }

  DDEcalPreshowerAlgo();
  void initialize(const DDNumericArguments& nArgs,
                  const DDVectorArguments& vArgs,
                  const DDMapArguments& mArgs,
                  const DDStringArguments& sArgs,
                  const DDStringVectorArguments& vsArgs) override;
  void execute(DDCompactView& pos) override;

private:
  void doLayers(DDCompactView& pos);
  void doLadders(DDCompactView& pos);
  void doSens(DDCompactView& pos);

  int nmat_;                            // number of preshower layers
  double thickness_;                    // overall thickness of the preshower envelope
  std::vector<std::string> materials_;  // materials of the presh-layers
  std::vector<std::string> layName_;    // names of the presh-layers
  std::vector<std::string> ladPfx_;     // name prefix for ladders
  std::string LaddMaterial_;            // ladd material - air
  std::vector<double> thickLayers_;
  std::vector<double> abs1stx;
  std::vector<double> abs1sty;
  std::vector<double> abs2ndx;
  std::vector<double> abs2ndy;
  std::vector<double> asym_ladd_;
  std::vector<double> rminVec;
  std::vector<double> rmaxVec;
  std::vector<double> noLaddInCol_;
  std::vector<double> startOfFirstLadd_;
  std::vector<std::string> types_l5_;
  std::vector<std::string> types_l4_;
  std::vector<double> ladd_l5_map_;
  std::vector<double> ladd_l4_map_;
  std::vector<std::string> typeOfLaddRow0;
  std::vector<std::string> typeOfLaddRow1;
  std::vector<std::string> typeOfLaddRow2;
  std::vector<std::string> typeOfLaddRow3;

  double zlead1_, zlead2_, zfoam1_, zfoam2_;
  double waf_intra_col_sep, waf_inter_col_sep, waf_active, wedge_length, wedge_offset, zwedge_ceramic_diff,
      ywedge_ceramic_diff, wedge_angle, box_thick, dee_separation, In_rad_Abs_Al, In_rad_Abs_Pb;
  double ladder_thick, yladder_1stwedge_diff, ladder_width, ladder_length, micromodule_length;
  double absAlX_X_, absAlX_Y_, absAlX_subtr1_Xshift_, absAlX_subtr1_Yshift_, rMax_Abs_Al_;
  double absAlY_X_, absAlY_Y_, absAlY_subtr1_Xshift_, absAlY_subtr1_Yshift_;
  double LdrBck_Length, LdrFrnt_Length, LdrFrnt_Offset, LdrBck_Offset, ceramic_length, wedge_back_thick;
};

DDEcalPreshowerAlgo::DDEcalPreshowerAlgo() : DDAlgorithm() {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("EcalGeom") << "DDEcalPreshowerAlgo info: Creating an instance";
#endif
}

void DDEcalPreshowerAlgo::initialize(const DDNumericArguments& nArgs,
                                     const DDVectorArguments& vArgs,
                                     const DDMapArguments& mArgs,
                                     const DDStringArguments& sArgs,
                                     const DDStringVectorArguments& vsArgs) {
  edm::LogVerbatim("EcalGeom") << "DDEcalPreshowerAlgo info: Initialize";

  asym_ladd_ = vArgs["ASYMETRIC_LADDER"];
  types_l5_ = vsArgs["TYPES_OF_LADD_L5"];
  types_l4_ = vsArgs["TYPES_OF_LADD_L4"];
  ladd_l5_map_ = vArgs["LADD_L5_MAP"];
  ladd_l4_map_ = vArgs["LADD_L4_MAP"];
  noLaddInCol_ = vArgs["NUMB_OF_LADD_IN_COL"];
  startOfFirstLadd_ = vArgs["START_OF_1ST_LADD"];
  typeOfLaddRow0 = vsArgs["TYPE_OF_LADD_1"];
  typeOfLaddRow1 = vsArgs["TYPE_OF_LADD_2"];
  typeOfLaddRow2 = vsArgs["TYPE_OF_LADD_3"];
  typeOfLaddRow3 = vsArgs["TYPE_OF_LADD_4"];
  thickLayers_ = vArgs["Layers"];
  thickness_ = double(nArgs["PRESH_Z_TOTAL"]);
  materials_ = vsArgs["LayMat"];
  layName_ = vsArgs["LayName"];
  rmaxVec = vArgs["R_MAX"];  // inner radii
  rminVec = vArgs["R_MIN"];  // outer radii
  waf_intra_col_sep = double(nArgs["waf_intra_col_sep"]);
  waf_inter_col_sep = double(nArgs["waf_inter_col_sep"]);
  waf_active = double(nArgs["waf_active"]);
  wedge_length = double(nArgs["wedge_length"]);
  wedge_offset = double(nArgs["wedge_offset"]);
  zwedge_ceramic_diff = double(nArgs["zwedge_ceramic_diff"]);
  ywedge_ceramic_diff = double(nArgs["ywedge_ceramic_diff"]);
  ceramic_length = double(nArgs["ceramic_length"]);
  wedge_angle = double(nArgs["wedge_angle"]);
  wedge_back_thick = double(nArgs["wedge_back_thick"]);
  ladder_thick = double(nArgs["ladder_thick"]);
  ladder_width = double(nArgs["ladder_width"]);
  micromodule_length = double(nArgs["micromodule_length"]);
  box_thick = double(nArgs["box_thick"]);
  abs1stx = vArgs["1ST_ABSX"];
  abs1sty = vArgs["1ST_ABSY"];
  abs2ndx = vArgs["2ND_ABSX"];
  abs2ndy = vArgs["2ND_ABSY"];
  ladPfx_ = vsArgs["LadPrefix"];
  LaddMaterial_ = sArgs["LadderMaterial"];
  LdrFrnt_Length = double(nArgs["LdrFrnt_Length"]);
  LdrFrnt_Offset = double(nArgs["LdrFrnt_Offset"]);
  LdrBck_Length = double(nArgs["LdrBck_Length"]);
  LdrBck_Offset = double(nArgs["LdrBck_Offset"]);
  dee_separation = double(nArgs["dee_sep"]);
  In_rad_Abs_Al = double(nArgs["R_MIN_Abs_Al"]);
  In_rad_Abs_Pb = double(nArgs["R_MIN_Abs_Pb"]);
  rMax_Abs_Al_ = double(nArgs["R_MAX_Abs_Al"]);
  absAlX_X_ = double(nArgs["AbsAlX_X"]);
  absAlX_Y_ = double(nArgs["AbsAlX_Y"]);
  absAlX_subtr1_Xshift_ = double(nArgs["AbsAlX_subtr1_Xshift"]);
  absAlX_subtr1_Yshift_ = double(nArgs["AbsAlX_subtr1_Yshift"]);
  absAlY_X_ = double(nArgs["AbsAlY_X"]);
  absAlY_Y_ = double(nArgs["AbsAlY_Y"]);
  absAlY_subtr1_Xshift_ = double(nArgs["AbsAlY_subtr1_Xshift"]);
  absAlY_subtr1_Yshift_ = double(nArgs["AbsAlY_subtr1_Yshift"]);
}

void DDEcalPreshowerAlgo::execute(DDCompactView& cpv) {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("EcalGeom") << "******** DDEcalPreshowerAlgo execute!";
#endif
  // creates all the tube-like layers of the preshower
  doLayers(cpv);
  // creates and places the ladders
  doLadders(cpv);
  // places the slicon strips in active silicon wafers
  doSens(cpv);
}

void DDEcalPreshowerAlgo::doLayers(DDCompactView& cpv) {
  double zpos = -thickness_ / 2., sdx(0), sdy(0), bdx(0), bdy(0);

  for (size_t i = 0; i < thickLayers_.size(); ++i) {
    int I = int(i) + 1;  // FOTRAN I (offset +1)

    double rIn(0), rOut(0), zHalf(0);

    // create the name
    DDName ddname(getLayName(i), "esalgo");  // namespace:name

    // cone dimensions
    rIn = rminVec[i];
    rOut = rmaxVec[i];
    zHalf = thickLayers_[i] / 2.;

    // create a logical part representing a single layer in the preshower
    DDSolid solid = DDSolidFactory::tubs(ddname, zHalf, rIn, rOut, 0., 360. * deg);
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("SFGeomX") << ddname << " Tubs " << zHalf << ":" << rIn << ":" << rOut << ":0:360";
#endif
    DDLogicalPart layer = DDLogicalPart(ddname, getMaterial(i), solid);

    // position the logical part w.r.t. the parent volume
    zpos += zHalf;

    // create a logical part representing a single layer in the preshower
    // skip layers with detectors, front and rear window
    if (I == 2 || I == 28 || I == 13 || I == 23) {
      zpos += zHalf;
      continue;
    }

    if (I == 12) {
      zlead1_ = zpos + zHalf;
    }
    if (I == 22) {
      zlead2_ = zpos + zHalf;
    }

    if (I == 10 || I == 20) {  // New lead shape

      int absz = 0;
      double outalbx, outalby, shiftR, outalbx2, outalby2, shiftR2;

      absz = int(abs1stx.size());
      if (I == 20)
        absz = int(abs2ndx.size());
      int cutabsx = -1, cutabsy = -1;

      DDName dd_tmp_name_b(getLayName(i) + "Lcut", "esalgo");
      DDName dd_tmp_name_c(getLayName(i) + "tmpb", "esalgo");
      DDName dd_tmp_name_d(getLayName(i) + "LinPb", "esalgo");
      DDName dd_tmp_name_e(getLayName(i) + "LinAl", "esalgo");
      DDName dd_tmp_name_f(getLayName(i) + "LOutAl", "esalgo");

      std::ostringstream tmp_Alname_fin;
      tmp_Alname_fin << getLayName(i) << "LtmpAl" << absz - 1;
      DDName dd_Alname_fin(tmp_Alname_fin.str(), "esalgo");

      DDSolid Outer_Al = DDSolid(dd_Alname_fin);

      DDLogicalPart layerFinOutAl = DDLogicalPart(dd_tmp_name_f, getMaterial(i - 1), Outer_Al);

      DDName dd_Alname_f(getLayName(i) + "LOutAl", "esalgo");
      DDName dd_Alname_g(getLayName(i) + "LOutAl2", "esalgo");
      DDName dd_Alname_h(getLayName(i) + "LOutAltmp", "esalgo");
      DDName dd_Alname_i(getLayName(i) + "LOutAltmp2", "esalgo");
      DDName dd_Alname_j(getLayName(i) + "LOutAltmp3", "esalgo");
      DDName dd_Alname_k(getLayName(i) + "LOutAltmp4", "esalgo");
      DDName dd_Alname_l(getLayName(i) + "LOutAltmp5", "esalgo");
      DDName dd_Alname_m(getLayName(i) + "LOutAltmp6", "esalgo");

      DDSolid Out_Al =
          DDSolidFactory::tubs(dd_Alname_f, zHalf - 0.1 * mm, rMax_Abs_Al_ - 20 * cm, rMax_Abs_Al_, 0., 90. * deg);
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("SFGeomX") << dd_Alname_f << " Tubs " << (zHalf - 0.1 * mm) << ":" << (rMax_Abs_Al_ - 20 * cm)
                                  << ":" << rMax_Abs_Al_ << ":0:90";
#endif

      outalbx = absAlX_X_ * 0.1;
      outalby = rMax_Abs_Al_ + 0.1 * mm - absAlX_subtr1_Yshift_;
      shiftR = absAlX_subtr1_Yshift_;
      if (I == 20) {
        outalbx = absAlY_X_ * 0.1;
        outalby = rMax_Abs_Al_ + 0.1 * mm - absAlY_subtr1_Yshift_;
        shiftR = absAlY_subtr1_Xshift_;
      }
      DDSolid OutAltmp = DDSolidFactory::box(dd_Alname_h, outalbx / 2 + 0.1 * mm, outalby / 2 + 0.1 * mm, zHalf);
      DDSolid Out_Altmp3 = DDSolidFactory::subtraction(
          dd_Alname_j, Out_Al, OutAltmp, DDTranslation(outalbx / 2, outalby / 2 + shiftR, 0), DDRotation());
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("SFGeomX") << dd_Alname_h << " Box " << (outalbx / 2 + 0.1 * mm) << ":"
                                  << (outalby / 2 + 0.1 * mm) << ":" << zHalf;
      edm::LogVerbatim("SFGeomX") << dd_Alname_j << " Subtraction " << Out_Al.name() << ":" << OutAltmp.name() << " at "
                                  << DDTranslation(outalbx / 2, outalby / 2 + shiftR, 0) << " no rotation";
#endif
      outalby2 = absAlX_Y_ * 0.1;
      outalbx2 = rMax_Abs_Al_ + 0.1 * mm - absAlX_subtr1_Xshift_;
      shiftR2 = absAlX_subtr1_Xshift_;
      if (I == 20) {
        outalby2 = absAlY_Y_ * 0.1;
        outalbx2 = rMax_Abs_Al_ + 0.1 * mm - absAlY_subtr1_Xshift_;
        shiftR2 = absAlY_subtr1_Xshift_;
      }
      DDSolid OutAltmp2 = DDSolidFactory::box(dd_Alname_i, outalbx2 / 2 + 0.1 * mm, outalby2 / 2 + 0.1 * mm, zHalf);
      DDSolid Out_Altmp4 = DDSolidFactory::subtraction(
          dd_Alname_k, Out_Altmp3, OutAltmp2, DDTranslation(outalbx2 / 2 + shiftR2, outalby2 / 2, 0), DDRotation());
      DDSolid Out_Altmp5 = DDSolidFactory::unionSolid(
          dd_Alname_l, Out_Altmp4, Out_Altmp4, DDTranslation(0, 0, 0), DDRotation("esalgo:RABS90"));
      DDSolid Out_Altmp6 = DDSolidFactory::unionSolid(
          dd_Alname_m, Out_Altmp5, Out_Altmp4, DDTranslation(0, 0, 0), DDRotation("esalgo:RABS180B"));
      DDSolid Out_Al2 = DDSolidFactory::unionSolid(
          dd_Alname_g, Out_Altmp6, Out_Altmp4, DDTranslation(0, 0, 0), DDRotation("esalgo:R180"));
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("SFGeomX") << dd_Alname_i << " Box " << (outalbx2 / 2 + 0.1 * mm) << ":"
                                  << (outalby2 / 2 + 0.1 * mm) << ":" << zHalf;
      edm::LogVerbatim("SFGeomX") << dd_Alname_k << " Subtraction " << Out_Altmp3.name() << ":" << OutAltmp2.name()
                                  << " at " << DDTranslation(outalbx2 / 2 + shiftR2, outalby2 / 2, 0) << " no rotation";
      edm::LogVerbatim("SFGeomX") << dd_Alname_l << " Union " << Out_Altmp4.name() << ":" << Out_Altmp4.name() << " at "
                                  << DDTranslation(0, 0, 0) << " rotation esalgo:RABS90";
      edm::LogVerbatim("SFGeomX") << dd_Alname_m << " Union " << Out_Altmp5.name() << ":" << Out_Altmp4.name() << " at "
                                  << DDTranslation(0, 0, 0) << " rotation esalgo:RABS180B";
      edm::LogVerbatim("SFGeomX") << dd_Alname_g << " Union " << Out_Altmp6.name() << ":" << Out_Altmp4.name() << " at "
                                  << DDTranslation(0, 0, 0) << " rotation esalgo:R180";
      edm::LogVerbatim("SFGeomX") << Outer_Al;
#endif

      for (int L = 0; L < absz; ++L) {
        int K = L;
        std::ostringstream tmp_name_b, tmp_name_b2, tmp_FAl_name_c, tmp_FAl_name_d1, tmp_FAl_name_d2, tmp_FAl_name_d3,
            tmp_FAl_name_d;
        tmp_name_b << getLayName(i) << "L" << K;
        tmp_name_b2 << getLayName(i) << "Lb2" << K;

        if (L == 0)
          tmp_FAl_name_c << getLayName(i) << "LOutAl2";
        if (L > 0)
          tmp_FAl_name_c << getLayName(i) << "LtmpAl" << K - 1;

        tmp_FAl_name_d1 << getLayName(i) << "LtmpAl" << K << "_1";
        tmp_FAl_name_d2 << getLayName(i) << "LtmpAl" << K << "_2";
        tmp_FAl_name_d3 << getLayName(i) << "LtmpAl" << K << "_3";
        tmp_FAl_name_d << getLayName(i) << "LtmpAl" << K;

        DDName dd_tmp_name_b(tmp_name_b.str(), "esalgo");
        DDName dd_tmp_name_b2(tmp_name_b2.str(), "esalgo");
        DDName dd_FAl_name_c(tmp_FAl_name_c.str(), "esalgo");
        DDName dd_FAl_name_d1(tmp_FAl_name_d1.str(), "esalgo");
        DDName dd_FAl_name_d2(tmp_FAl_name_d2.str(), "esalgo");
        DDName dd_FAl_name_d3(tmp_FAl_name_d3.str(), "esalgo");
        DDName dd_FAl_name_d(tmp_FAl_name_d.str(), "esalgo");

        if (L == 0)
          bdx = std::abs(abs1stx[K]) / 2;
        if (L > 0)
          bdx = std::abs(abs1stx[K] - abs1stx[K - 1]) / 2;
        bdy = abs1sty[K];
        if (abs1stx[K] < rIn + 30 * cm) {
          bdy = abs1sty[K] / 2 - 30 * cm;
          cutabsx = K;
        }

        if (I == 20) {
          if (L == 0)
            bdx = std::abs(abs2ndx[K]) / 2;
          if (L > 0)
            bdx = std::abs(abs2ndx[K] - abs2ndx[K - 1]) / 2;
          bdy = abs2ndy[K];
        }

        if ((abs2ndx[K] < rIn + 30 * cm) && I == 20) {
          bdy = abs2ndy[K] / 2 - 30 * cm;
          cutabsy = K;
        }

        DDSolid solid_b = DDSolidFactory::box(dd_tmp_name_b, bdx, bdy, zHalf);
        DDSolid solid_b2 = DDSolidFactory::box(dd_tmp_name_b2, bdx + 0.1 * mm, bdy + 0.1 * mm, zHalf);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("SFGeomX") << dd_tmp_name_b << " Box " << bdx << ":" << bdy << ":" << zHalf;
        edm::LogVerbatim("SFGeomX") << dd_tmp_name_b2 << " Box " << (bdx + 0.1 * mm) << ":" << (bdy + 0.1 * mm) << ":"
                                    << zHalf;
#endif

        sdx = abs1stx[K] - bdx;
        sdy = 0;
        if (abs1stx[K] < rIn + 30 * cm)
          sdy = abs1sty[K] - bdy;

        if (I == 20) {
          sdx = abs2ndx[K] - bdx;
          sdy = 0;
        }
        if ((abs2ndx[K] < rIn + 30 * cm) && I == 20)
          sdy = abs2ndy[K] - bdy;

        DDLogicalPart layer = DDLogicalPart(dd_tmp_name_b, getMaterial(i), solid_b);

        cpv.position(layer, layerFinOutAl, 1, DDTranslation(sdx, sdy, 0), DDRotation());
        cpv.position(layer, layerFinOutAl, 2, DDTranslation(-sdx, sdy, 0), DDRotation());
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("SFGeom") << layer.name() << " copy 1 in " << layerFinOutAl.name() << " at "
                                   << DDTranslation(sdx, sdy, 0) << " no rotation";
        edm::LogVerbatim("SFGeom") << layer.name() << " copy 2 in " << layerFinOutAl.name() << " at "
                                   << DDTranslation(-sdx, sdy, 0) << " no rotation";
#endif
        DDSolid solid_c = DDSolid(dd_FAl_name_c);
        DDSolid solid_d1 =
            DDSolidFactory::unionSolid(dd_FAl_name_d1, solid_c, solid_b2, DDTranslation(sdx, sdy, 0), DDRotation());
        DDSolid solid_d2 =
            DDSolidFactory::unionSolid(dd_FAl_name_d, solid_d1, solid_b2, DDTranslation(-sdx, -sdy, 0), DDRotation());
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("SFGeomX") << dd_FAl_name_d1 << " Union " << solid_c.name() << ":" << solid_b2.name() << " at "
                                    << DDTranslation(sdx, sdy, 0) << " no rotation";
        edm::LogVerbatim("SFGeomX") << dd_FAl_name_d << " Union " << solid_d1.name() << ":" << solid_b2.name() << " at "
                                    << DDTranslation(-sdx, -sdy, 0) << " no rotation";
#endif
        if (((abs1stx[K] < rIn + 30 * cm) && I == 10) || ((abs2ndx[K] < rIn + 30 * cm) && I == 20)) {
          cpv.position(layer, layerFinOutAl, 3, DDTranslation(sdx, -sdy, 0), DDRotation());
          cpv.position(layer, layerFinOutAl, 4, DDTranslation(-sdx, -sdy, 0), DDRotation());
#ifdef EDM_ML_DEBUG
          edm::LogVerbatim("SFGeom") << layer.name() << " copy 3 in " << layerFinOutAl.name() << " at "
                                     << DDTranslation(sdx, -sdy, 0) << " no rotation";
          edm::LogVerbatim("SFGeom") << layer.name() << " copy 4 in " << layerFinOutAl.name() << " at "
                                     << DDTranslation(-sdx, -sdy, 0) << " no rotation";
#endif
          DDSolid solid_c = DDSolid(dd_FAl_name_c);
          DDSolid solid_d1 =
              DDSolidFactory::unionSolid(dd_FAl_name_d1, solid_c, solid_b2, DDTranslation(sdx, sdy, 0), DDRotation());
          DDSolid solid_d2 =
              DDSolidFactory::unionSolid(dd_FAl_name_d2, solid_d1, solid_b2, DDTranslation(sdx, -sdy, 0), DDRotation());
          DDSolid solid_d3 =
              DDSolidFactory::unionSolid(dd_FAl_name_d3, solid_d2, solid_b2, DDTranslation(-sdx, sdy, 0), DDRotation());
          DDSolid solid_d4 =
              DDSolidFactory::unionSolid(dd_FAl_name_d, solid_d3, solid_b2, DDTranslation(-sdx, -sdy, 0), DDRotation());
#ifdef EDM_ML_DEBUG
          edm::LogVerbatim("SFGeomX") << dd_FAl_name_d1 << " Union " << solid_c.name() << ":" << solid_b2.name()
                                      << " at " << DDTranslation(sdx, sdy, 0) << " no rotation";
          edm::LogVerbatim("SFGeomX") << dd_FAl_name_d2 << " Union " << solid_d1.name() << ":" << solid_b2.name()
                                      << " at " << DDTranslation(sdx, -sdy, 0) << " no rotation";
          edm::LogVerbatim("SFGeomX") << dd_FAl_name_d3 << " Union " << solid_d2.name() << ":" << solid_b2.name()
                                      << " at " << DDTranslation(-sdx, sdy, 0) << " no rotation";
          edm::LogVerbatim("SFGeomX") << dd_FAl_name_d << " Union " << solid_d3.name() << ":" << solid_b2.name()
                                      << " at " << DDTranslation(-sdx, -sdy, 0) << " no rotation";
#endif
        }
      }

      bdx = abs1stx[cutabsx];
      if (I == 20)
        bdx = abs2ndx[cutabsy];
      bdy = 2 * 30 * cm;

      DDSolid solidcut = DDSolidFactory::box(dd_tmp_name_b, bdx, bdy, zHalf);

      DDSolid iner = DDSolidFactory::tubs(dd_tmp_name_c, zHalf + 0.1 * mm, 0, In_rad_Abs_Pb, 0., 360. * deg);

      DDSolid final = DDSolidFactory::subtraction(dd_tmp_name_d, solidcut, iner, DDTranslation(0, 0, 0), DDRotation());
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("SFGeomX") << dd_tmp_name_b << " Box " << bdx << ":" << bdy << ":" << zHalf;
      edm::LogVerbatim("SFGeomX") << dd_tmp_name_c << " Tubs " << (zHalf + 0.1 * mm) << ":0:" << In_rad_Abs_Pb
                                  << ":0:360";
      edm::LogVerbatim("SFGeomX") << dd_tmp_name_d << " Subtraction " << solidcut.name() << ":" << iner.name()
                                  << " at (0,0,0) no rotation";
#endif

      DDLogicalPart layer = DDLogicalPart(dd_tmp_name_d, getMaterial(i), final);
      cpv.position(layer, parent(), 1, DDTranslation(0, 0, zpos), DDRotation());
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("SFGeom") << layer.name() << " copy 1 in " << parent().name() << " at "
                                 << DDTranslation(0, 0, zpos) << " no rotation";
#endif
      DDSolid iner_Al =
          DDSolidFactory::tubs(dd_tmp_name_e, zHalf, In_rad_Abs_Al, In_rad_Abs_Pb - 0.01 * mm, 0., 360. * deg);
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("SFGeomX") << dd_tmp_name_e << " Tubs " << zHalf << ":" << In_rad_Abs_Al << ":"
                                  << (In_rad_Abs_Pb - 0.01 * mm) << ":0:360";
#endif

      DDLogicalPart layerAl = DDLogicalPart(dd_tmp_name_e, getMaterial(i - 1), iner_Al);
      cpv.position(layerAl, parent(), 1, DDTranslation(0, 0, zpos), DDRotation());
      cpv.position(layerFinOutAl, parent(), 1, DDTranslation(0, 0, zpos), DDRotation());
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("SFGeom") << layerAl.name() << " copy 1 in " << parent().name() << " at "
                                 << DDTranslation(0, 0, zpos) << " no rotation";
      edm::LogVerbatim("SFGeom") << layerFinOutAl.name() << " copy 1 in " << parent().name() << " at "
                                 << DDTranslation(0, 0, zpos) << " no rotation";
#endif
    } else {
      cpv.position(layer, parent(), 1, DDTranslation(0., 0., zpos), DDRotation());
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("SFGeom") << layer.name() << " copy 1 in " << parent().name() << " at "
                                 << DDTranslation(0, 0, zpos) << " no rotation";
#endif
    }
    zpos += zHalf;
  }
}

void DDEcalPreshowerAlgo::doLadders(DDCompactView& cpv) {
  double xpos(0), ypos(0), zpos(0), sdx(0), sdy(0), sdz(0);
  double prev_length_(0), ladder_new_length_(0), ladd_shift_(0), ladder_length(0);
  int enb(0), swed_scopy_glob(0);
  double sdxe[50] = {0}, sdye[50] = {0}, sdze[50] = {0};
  double sdxe2[50] = {0}, sdye2[50] = {0}, sdze2[50] = {0}, sdxe3[50] = {0}, sdye3[50] = {0}, sdze3[50] = {0};

  for (int M = 0; M < int(types_l5_.size() + types_l4_.size()); M++) {
    int scopy(0);
    double boxax(0), boxay(0), boxaz(0);
    int ladd_not_plain(0), ladd_subtr_no(0), ladd_upper(0), ladd_side(0);

    DDSolid solid_lfront = DDSolidFactory::trap(DDName("LDRFRNT", "esalgo"),
                                                LdrFrnt_Length / 2,                                          // pDz
                                                -wedge_angle,                                                // pTheta
                                                0,                                                           // pPhi
                                                ladder_width / 2,                                            // pDy1
                                                (ladder_thick) / 2,                                          // pDx1
                                                (ladder_thick) / 2,                                          // pDx2
                                                0,                                                           // pAlp1
                                                ladder_width / 2,                                            // pDy2
                                                (ladder_thick - ceramic_length * sin(wedge_angle * 2)) / 2,  // pDx3
                                                (ladder_thick - ceramic_length * sin(wedge_angle * 2)) / 2,  // pDx4
                                                0);

    DDSolid solid_lbck = DDSolidFactory::trap(DDName("LDRBCK", "esalgo"),
                                              LdrBck_Length / 2,                                   // pDz
                                              -wedge_angle,                                        // pTheta
                                              0,                                                   // pPhi
                                              ladder_width / 2,                                    // pDy1
                                              (box_thick / cos(wedge_angle * 2) + 0.02 * mm) / 2,  // pDx1
                                              (box_thick / cos(wedge_angle * 2) + 0.02 * mm) / 2,  // pDx2
                                              0,                                                   // pAlp1
                                              ladder_width / 2,                                    // pDy2
                                              (ladder_thick - wedge_back_thick) / 2,               // pDx3
                                              (ladder_thick - wedge_back_thick) / 2,               // pDx4
                                              0);

    DDSolid solid_lfhalf = DDSolidFactory::trap(DDName("LDRFHALF", "esalgo"),
                                                LdrFrnt_Length / 2,                                          // pDz
                                                -wedge_angle,                                                // pTheta
                                                0,                                                           // pPhi
                                                (ladder_width / 2) / 2,                                      // pDy1
                                                (ladder_thick) / 2,                                          // pDx1
                                                (ladder_thick) / 2,                                          // pDx2
                                                0,                                                           // pAlp1
                                                (ladder_width / 2) / 2,                                      // pDy2
                                                (ladder_thick - ceramic_length * sin(wedge_angle * 2)) / 2,  // pDx3
                                                (ladder_thick - ceramic_length * sin(wedge_angle * 2)) / 2,  // pDx4
                                                0);

    DDSolid solid_lbhalf = DDSolidFactory::trap(DDName("LDRBHALF", "esalgo"),
                                                LdrBck_Length / 2,                                   // pDz
                                                -wedge_angle,                                        // pTheta
                                                0,                                                   // pPhi
                                                (ladder_width / 2) / 2,                              // pDy1
                                                (box_thick / cos(wedge_angle * 2) + 0.02 * mm) / 2,  // pDx1
                                                (box_thick / cos(wedge_angle * 2) + 0.02 * mm) / 2,  // pDx2
                                                0,                                                   // pAlp1
                                                (ladder_width / 2) / 2,                              // pDy2
                                                (ladder_thick - wedge_back_thick) / 2,               // pDx3
                                                (ladder_thick - wedge_back_thick) / 2,               // pDx4
                                                0);

    DDSolid solid_lfhtrunc =
        DDSolidFactory::trap(DDName("LDRFHTR", "esalgo"),
                             (LdrFrnt_Length - waf_active) / 2,                                          // pDz
                             -wedge_angle,                                                               // pTheta
                             0,                                                                          // pPhi
                             (ladder_width / 2) / 2,                                                     // pDy1
                             (ladder_thick) / 2,                                                         // pDx1
                             (ladder_thick) / 2,                                                         // pDx2
                             0,                                                                          // pAlp1
                             (ladder_width / 2) / 2,                                                     // pDy2
                             (ladder_thick - (ceramic_length - waf_active) * sin(wedge_angle * 2)) / 2,  // pDx3
                             (ladder_thick - (ceramic_length - waf_active) * sin(wedge_angle * 2)) / 2,  // pDx4
                             0);
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("SFGeomX") << DDName("LDRFRNT", "esalgo") << " Trap " << (LdrFrnt_Length / 2) << ":"
                                << -wedge_angle / CLHEP::deg << ":0:" << (ladder_width / 2) << ":" << (ladder_thick / 2)
                                << ":" << (ladder_thick / 2) << ":0:" << (ladder_width / 2) << ":"
                                << ((ladder_thick - ceramic_length * sin(wedge_angle * 2)) / 2) << ":"
                                << ((ladder_thick - ceramic_length * sin(wedge_angle * 2)) / 2) << ":0";
    edm::LogVerbatim("SFGeomX") << DDName("LDRBCK", "esalgo") << " Trap " << (LdrBck_Length / 2) << ":"
                                << -wedge_angle / CLHEP::deg << ":0:" << (ladder_width / 2) << ":"
                                << ((box_thick / cos(wedge_angle * 2) + 0.02 * mm) / 2) << ":"
                                << ((box_thick / cos(wedge_angle * 2) + 0.02 * mm) / 2) << ":0:" << (ladder_width / 2)
                                << ":" << ((ladder_thick - wedge_back_thick) / 2) << ":"
                                << ((ladder_thick - wedge_back_thick) / 2) << ":0";
    edm::LogVerbatim("SFGeomX") << DDName("LDRFHALF", "esalgo") << " Trap " << (LdrFrnt_Length / 2) << ":"
                                << -wedge_angle / CLHEP::deg << ":0:" << ((ladder_width / 2) / 2) << ":"
                                << (ladder_thick / 2) << ":" << (ladder_thick / 2) << ":0:" << ((ladder_width / 2) / 2)
                                << ":" << ((ladder_thick - ceramic_length * sin(wedge_angle * 2)) / 2) << ":"
                                << ((ladder_thick - ceramic_length * sin(wedge_angle * 2)) / 2) << ":0";
    edm::LogVerbatim("SFGeomX") << DDName("LDRBHALF", "esalgo") << " Trap " << (LdrBck_Length / 2) << ":"
                                << -wedge_angle / CLHEP::deg << ":0:" << ((ladder_width / 2) / 2) << ":"
                                << ((box_thick / cos(wedge_angle * 2) + 0.02 * mm) / 2) << ":"
                                << ((box_thick / cos(wedge_angle * 2) + 0.02 * mm) / 2)
                                << ":0:" << ((ladder_width / 2) / 2) << ":" << ((ladder_thick - wedge_back_thick) / 2)
                                << ":" << ((ladder_thick - wedge_back_thick) / 2) << ":0";
    edm::LogVerbatim("SFGeomX") << DDName("LDRFHTR", "esalgo") << " Trap " << ((LdrFrnt_Length - waf_active) / 2) << ":"
                                << -wedge_angle / CLHEP::deg << ":0:" << ((ladder_width / 2) / 2) << ":"
                                << ((ladder_thick) / 2) << ":" << ((ladder_thick) / 2)
                                << ":0:" << ((ladder_width / 2) / 2) << ":"
                                << ((ladder_thick - (ceramic_length - waf_active) * sin(wedge_angle * 2)) / 2) << ":"
                                << ((ladder_thick - (ceramic_length - waf_active) * sin(wedge_angle * 2)) / 2) << ":0";
#endif

    // Creation of ladders with 5 micromodules length

    if (M < int(types_l5_.size())) {
      for (int i = 0; i <= 1; i++) {
        for (int j = 0; j <= 3; j++) {
          if (ladd_l5_map_[(i + j * 2 + M * 10)] != 1) {
            ladd_not_plain = 1;
            ladd_subtr_no++;
            if (j > 1)
              ladd_upper = 1;
            ladd_side = i;
          }
        }
      }
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("SFGeomX") << "Side " << ladd_side << ":" << ladd_upper << ":" << ladd_subtr_no << ":"
                                  << ladd_not_plain << " Index " << M << ":" << types_l5_.size();
#endif
      DDName ddname(getLadPrefix(0) + types_l5_[M], "esalgo");
      ladder_length = micromodule_length + 4 * waf_active + 0.1 * mm;

      if (ladd_not_plain) {
        //          enb++;
        std::ostringstream tmp_name_5b, tmp_name_5c, tmp_name_5d;
        if (ladd_upper) {
        }  //upper
        else {
          enb++;
          std::ostringstream tmp_name_5b, tmp_name_5c, tmp_name_5d;
          DDName dd_tmp_name_5a(getLadPrefix(2), "esalgo");
          tmp_name_5b << getLadPrefix(3) << enb;
          DDName dd_tmp_name_5b(tmp_name_5b.str(), "esalgo");
          tmp_name_5c << getLadPrefix(4) << enb;
          DDName dd_tmp_name_5c(tmp_name_5c.str(), "esalgo");
          tmp_name_5d << getLadPrefix(5) << enb;
          DDName dd_tmp_name_5d(tmp_name_5d.str(), "esalgo");

          DDName dd_tmp_name_5e(getLadPrefix(6), "esalgo");

          boxay = ladder_length - LdrFrnt_Length - LdrBck_Length;
          boxax = ladder_width;
          boxaz = ladder_thick;

          DDSolid solid_5a = DDSolidFactory::box(dd_tmp_name_5a, boxax / 2, boxay / 2, boxaz / 2.);
#ifdef EDM_ML_DEBUG
          edm::LogVerbatim("SFGeomX") << dd_tmp_name_5a << " Box " << (boxax / 2) << ":" << (boxay / 2) << ":"
                                      << (boxaz / 2);
#endif
          if (ladd_side == 0)
            sdxe[enb] = ladder_width / 4;
          sdye[enb] = -boxay / 2 - LdrFrnt_Length / 2;
          sdze[enb] = -ladder_thick / 2. + LdrFrnt_Offset;
          if (ladd_side == 1)
            sdxe[enb] = -ladder_width / 4;

          DDSolid solid_5b = DDSolidFactory::unionSolid(dd_tmp_name_5b,
                                                        solid_5a,
                                                        solid_lfhalf,
                                                        DDTranslation(sdxe[enb], sdye[enb], sdze[enb]),
                                                        DDRotation("esalgo:RM1299"));
#ifdef EDM_ML_DEBUG
          edm::LogVerbatim("SFGeomX") << dd_tmp_name_5b << " Union " << solid_5a.name() << ":" << solid_lfhalf.name()
                                      << " at " << DDTranslation(sdxe[enb], sdye[enb], sdze[enb])
                                      << " rotation esalgo:RM1299";
#endif

          if (ladd_side == 0)
            sdxe2[enb] = -ladder_width / 4;
          sdye2[enb] = -boxay / 2 - LdrFrnt_Length / 2 + waf_active / 2;
          sdze2[enb] = -ladder_thick / 2. + LdrFrnt_Offset + (waf_active * sin(wedge_angle * 2)) / 4;
          if (ladd_side == 1)
            sdxe2[enb] = ladder_width / 4;

          DDSolid solid_5c = DDSolidFactory::unionSolid(dd_tmp_name_5c,
                                                        solid_5b,
                                                        solid_lfhtrunc,
                                                        DDTranslation(sdxe2[enb], sdye2[enb], sdze2[enb]),
                                                        DDRotation("esalgo:RM1299"));
#ifdef EDM_ML_DEBUG
          edm::LogVerbatim("SFGeomX") << dd_tmp_name_5c << " Union " << solid_5b.name() << ":" << solid_lfhtrunc.name()
                                      << " at " << DDTranslation(sdxe2[enb], sdye2[enb], sdze2[enb])
                                      << " rotation esalgo:RM1299";
#endif

          sdxe3[enb] = 0;
          sdye3[enb] = boxay / 2 + LdrBck_Length / 2;
          sdze3[enb] = -ladder_thick / 2. + LdrBck_Offset;
          DDSolid solid = DDSolidFactory::unionSolid(ddname,
                                                     solid_5c,
                                                     solid_lbck,
                                                     DDTranslation(sdxe3[enb], sdye3[enb], sdze3[enb]),
                                                     DDRotation("esalgo:RM1299"));
#ifdef EDM_ML_DEBUG
          edm::LogVerbatim("SFGeomX") << ddname << " Union " << solid_5c.name() << ":" << solid_lbck.name() << " at "
                                      << DDTranslation(sdxe3[enb], sdye3[enb], sdze3[enb]) << " rotation esalgo:RM1299";
#endif

          DDLogicalPart ladder = DDLogicalPart(ddname, getLaddMaterial(), solid);
          DDName ddname2(getLadPrefix(1) + types_l5_[M], "esalgo");
          DDLogicalPart ladder2 = DDLogicalPart(ddname2, getLaddMaterial(), solid);
        }

      }  //end of not plain ladder shape
      else {
        DDName dd_tmp_name_5pa(getLadPrefix(2) + "5p", "esalgo");
        DDName dd_tmp_name_5pb(getLadPrefix(3) + "5p", "esalgo");

        boxay = ladder_length - LdrFrnt_Length - LdrBck_Length;
        boxax = ladder_width;
        boxaz = ladder_thick;

        DDSolid solid_5pa = DDSolidFactory::box(dd_tmp_name_5pa, boxax / 2, boxay / 2, boxaz / 2.);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("SFGeomX") << dd_tmp_name_5pa << " Box " << (boxax / 2) << ":" << (boxay / 2) << ":"
                                    << (boxaz / 2);
#endif
        sdx = 0;
        sdy = -boxay / 2 - LdrFrnt_Length / 2;
        sdz = -ladder_thick / 2. + LdrFrnt_Offset;

        DDSolid solid_5pb = DDSolidFactory::unionSolid(
            dd_tmp_name_5pb, solid_5pa, solid_lfront, DDTranslation(sdx, sdy, sdz), DDRotation("esalgo:RM1299"));
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("SFGeomX") << dd_tmp_name_5pb << " Union " << solid_5pa.name() << ":" << solid_lfront.name()
                                    << " at " << DDTranslation(sdx, sdy, sdz) << " rotation esalgo:RM1299";
#endif

        sdx = 0;
        sdy = boxay / 2 + LdrBck_Length / 2;
        sdz = -ladder_thick / 2. + LdrBck_Offset;

        DDSolid solid = DDSolidFactory::unionSolid(
            ddname, solid_5pb, solid_lbck, DDTranslation(sdx, sdy, sdz), DDRotation("esalgo:RM1299"));
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("SFGeomX") << ddname << " Union " << solid_5pb.name() << ":" << solid_lbck.name() << " at "
                                    << DDTranslation(sdx, sdy, sdz) << " rotation esalgo:RM1299";
#endif

        DDLogicalPart ladder = DDLogicalPart(ddname, getLaddMaterial(), solid);
        DDName ddname2(getLadPrefix(1) + types_l5_[M], "esalgo");
        DDLogicalPart ladder2 = DDLogicalPart(ddname2, getLaddMaterial(), solid);
      }
    }

    // Creation of ladders with 4 micromodules length

    if (M >= int(types_l5_.size())) {
      int d = M - types_l5_.size();

      for (int i = 0; i <= 1; i++) {
        for (int j = 0; j <= 3; j++) {
          if (ladd_l4_map_[(i + j * 2 + (M - types_l5_.size()) * 8)] != 1) {
            ladd_not_plain = 1;
            ladd_subtr_no++;
            if (j > 1)
              ladd_upper = 1;
            ladd_side = i;
          }
        }
      }
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("SFGeomX") << "Side " << ladd_side << ":" << ladd_upper << ":" << ladd_subtr_no << ":"
                                  << ladd_not_plain << " Index " << M << ":" << types_l5_.size();
#endif
      DDName ddname(getLadPrefix(0) + types_l4_[d], "esalgo");
      ladder_length = micromodule_length + 3 * waf_active + 0.1 * mm;

      if (ladd_not_plain) {
        std::ostringstream tmp_name_b, tmp_name_c, tmp_name_d;
        if (ladd_upper) {
          enb++;

          DDName dd_tmp_name_a(getLadPrefix(7), "esalgo");
          tmp_name_b << getLadPrefix(8) << enb;
          DDName dd_tmp_name_b(tmp_name_b.str(), "esalgo");
          tmp_name_c << getLadPrefix(9) << enb;
          DDName dd_tmp_name_c(tmp_name_c.str(), "esalgo");
          tmp_name_d << getLadPrefix(10) << enb;
          DDName dd_tmp_name_d(tmp_name_d.str(), "esalgo");
          DDName dd_tmp_name_e(getLadPrefix(11), "esalgo");

          boxay = ladder_length - LdrFrnt_Length - LdrBck_Length;
          boxax = ladder_width;
          boxaz = ladder_thick;
          DDSolid solid_a = DDSolidFactory::box(dd_tmp_name_a, boxax / 2, boxay / 2, boxaz / 2.);
#ifdef EDM_ML_DEBUG
          edm::LogVerbatim("SFGeomX") << dd_tmp_name_a << " Box " << (boxax / 2) << ":" << (boxay / 2) << ":"
                                      << (boxaz / 2);
#endif

          sdxe[enb] = 0;
          sdye[enb] = -boxay / 2 - LdrFrnt_Length / 2;
          sdze[enb] = -ladder_thick / 2. + LdrFrnt_Offset;
          DDSolid solid_b = DDSolidFactory::unionSolid(dd_tmp_name_b,
                                                       solid_a,
                                                       solid_lfront,
                                                       DDTranslation(sdxe[enb], sdye[enb], sdze[enb]),
                                                       DDRotation("esalgo:RM1299"));
#ifdef EDM_ML_DEBUG
          edm::LogVerbatim("SFGeomX") << dd_tmp_name_b << " Union " << solid_a.name() << ":" << solid_lfront.name()
                                      << " at " << DDTranslation(sdxe[enb], sdye[enb], sdze[enb])
                                      << " rotation esalgo:RM1299";
#endif

          if (ladd_side == 0)
            sdxe2[enb] = ladder_width / 4;
          sdye2[enb] = boxay / 2 + LdrBck_Length / 2;
          sdze2[enb] = -ladder_thick / 2. + LdrBck_Offset;
          if (ladd_side == 1)
            sdxe2[enb] = -ladder_width / 4;
          DDSolid solid = DDSolidFactory::unionSolid(ddname,
                                                     solid_b,
                                                     solid_lbhalf,
                                                     DDTranslation(sdxe2[enb], sdye2[enb], sdze2[enb]),
                                                     DDRotation("esalgo:RM1299"));
#ifdef EDM_ML_DEBUG
          edm::LogVerbatim("SFGeomX") << ddname << " Union " << solid_b.name() << ":" << solid_lbhalf.name() << " at "
                                      << DDTranslation(sdxe2[enb], sdye2[enb], sdze2[enb]) << " rotation esalgo:RM1299";
#endif

          DDLogicalPart ladder = DDLogicalPart(ddname, getLaddMaterial(), solid);
          DDName ddname2(getLadPrefix(1) + types_l4_[d], "esalgo");
          DDLogicalPart ladder2 = DDLogicalPart(ddname2, getLaddMaterial(), solid);

        }  //upper
        else {
          if (ladd_subtr_no > 1) {
            enb++;

            DDName dd_tmp_name_a(getLadPrefix(7), "esalgo");
            tmp_name_b << getLadPrefix(8) << enb;
            DDName dd_tmp_name_b(tmp_name_b.str(), "esalgo");
            tmp_name_c << getLadPrefix(9) << enb;
            DDName dd_tmp_name_c(tmp_name_c.str(), "esalgo");
            tmp_name_d << getLadPrefix(10) << enb;
            DDName dd_tmp_name_d(tmp_name_d.str(), "esalgo");
            DDName dd_tmp_name_e(getLadPrefix(11), "esalgo");

            boxay = ladder_length - LdrFrnt_Length - LdrBck_Length;
            boxax = ladder_width;
            boxaz = ladder_thick;

            DDSolid solid_a = DDSolidFactory::box(dd_tmp_name_a, boxax / 2, boxay / 2, boxaz / 2.);
#ifdef EDM_ML_DEBUG
            edm::LogVerbatim("SFGeomX") << dd_tmp_name_a << " Box " << (boxax / 2) << ":" << (boxay / 2) << ":"
                                        << (boxaz / 2);
#endif
            if (ladd_side == 0)
              sdxe[enb] = ladder_width / 4;
            sdye[enb] = -boxay / 2 - LdrFrnt_Length / 2;
            sdze[enb] = -ladder_thick / 2. + LdrFrnt_Offset;
            if (ladd_side == 1)
              sdxe[enb] = -ladder_width / 4;

            DDSolid solid_b = DDSolidFactory::unionSolid(dd_tmp_name_b,
                                                         solid_a,
                                                         solid_lfhalf,
                                                         DDTranslation(sdxe[enb], sdye[enb], sdze[enb]),
                                                         DDRotation("esalgo:RM1299"));
#ifdef EDM_ML_DEBUG
            edm::LogVerbatim("SFGeomX") << dd_tmp_name_b << " Union " << solid_a.name() << ":" << solid_lfhalf.name()
                                        << " at " << DDTranslation(sdxe[enb], sdye[enb], sdze[enb])
                                        << " rotation esalgo:RM1299";
#endif

            sdxe2[enb] = 0;
            sdye2[enb] = boxay / 2 + LdrBck_Length / 2;
            sdze2[enb] = -ladder_thick / 2. + LdrBck_Offset;

            DDSolid solid = DDSolidFactory::unionSolid(ddname,
                                                       solid_b,
                                                       solid_lbck,
                                                       DDTranslation(sdxe2[enb], sdye2[enb], sdze2[enb]),
                                                       DDRotation("esalgo:RM1299"));
#ifdef EDM_ML_DEBUG
            edm::LogVerbatim("SFGeomX") << ddname << " Union " << solid_b.name() << ":" << solid_lbck.name() << " at "
                                        << DDTranslation(sdxe2[enb], sdye2[enb], sdze2[enb])
                                        << " rotation esalgo:RM1299";
#endif

            DDLogicalPart ladder = DDLogicalPart(ddname, getLaddMaterial(), solid);
            DDName ddname2(getLadPrefix(1) + types_l4_[d], "esalgo");
            DDLogicalPart ladder2 = DDLogicalPart(ddname2, getLaddMaterial(), solid);
          } else {
            enb++;
            DDName dd_tmp_name_a(getLadPrefix(7), "esalgo");
            tmp_name_b << getLadPrefix(8) << enb;
            DDName dd_tmp_name_b(tmp_name_b.str(), "esalgo");
            tmp_name_c << getLadPrefix(9) << enb;
            DDName dd_tmp_name_c(tmp_name_c.str(), "esalgo");
            tmp_name_d << getLadPrefix(10) << enb;
            DDName dd_tmp_name_d(tmp_name_d.str(), "esalgo");
            DDName dd_tmp_name_e(getLadPrefix(11), "esalgo");

            boxay = ladder_length - LdrFrnt_Length - LdrBck_Length;
            boxax = ladder_width;
            boxaz = ladder_thick;
            DDSolid solid_a = DDSolidFactory::box(dd_tmp_name_a, boxax / 2, boxay / 2, boxaz / 2.);
#ifdef EDM_ML_DEBUG
            edm::LogVerbatim("SFGeomX") << dd_tmp_name_a << " Box " << (boxax / 2) << ":" << (boxay / 2) << ":"
                                        << (boxaz / 2);
#endif
            if (ladd_side == 0)
              sdxe[enb] = ladder_width / 4;
            sdye[enb] = -boxay / 2 - LdrFrnt_Length / 2;
            sdze[enb] = -ladder_thick / 2. + LdrFrnt_Offset;
            if (ladd_side == 1)
              sdxe[enb] = -ladder_width / 4;

            DDSolid solid_b = DDSolidFactory::unionSolid(dd_tmp_name_b,
                                                         solid_a,
                                                         solid_lfhalf,
                                                         DDTranslation(sdxe[enb], sdye[enb], sdze[enb]),
                                                         DDRotation("esalgo:RM1299"));
#ifdef EDM_ML_DEBUG
            edm::LogVerbatim("SFGeomX") << dd_tmp_name_b << " Union " << solid_a.name() << ":" << solid_lfhalf.name()
                                        << " at " << DDTranslation(sdxe[enb], sdye[enb], sdze[enb])
                                        << " rotation esalgo:RM1299";
#endif

            if (ladd_side == 0)
              sdxe2[enb] = -ladder_width / 4;
            sdye2[enb] = -boxay / 2 - LdrFrnt_Length / 2 + waf_active / 2;
            sdze2[enb] = -ladder_thick / 2. + LdrFrnt_Offset + (waf_active * sin(wedge_angle * 2)) / 4;
            if (ladd_side == 1)
              sdxe2[enb] = ladder_width / 4;

            DDSolid solid_c = DDSolidFactory::unionSolid(dd_tmp_name_c,
                                                         solid_b,
                                                         solid_lfhtrunc,
                                                         DDTranslation(sdxe2[enb], sdye2[enb], sdze2[enb]),
                                                         DDRotation("esalgo:RM1299"));
#ifdef EDM_ML_DEBUG
            edm::LogVerbatim("SFGeomX") << dd_tmp_name_c << " Union " << solid_b.name() << ":" << solid_lfhtrunc.name()
                                        << " at " << DDTranslation(sdxe2[enb], sdye2[enb], sdze2[enb])
                                        << " rotation esalgo:RM1299";
#endif

            sdxe3[enb] = 0;
            sdye3[enb] = boxay / 2 + LdrBck_Length / 2;
            sdze3[enb] = -ladder_thick / 2. + LdrBck_Offset;
            DDSolid solid = DDSolidFactory::unionSolid(ddname,
                                                       solid_c,
                                                       solid_lbck,
                                                       DDTranslation(sdxe3[enb], sdye3[enb], sdze3[enb]),
                                                       DDRotation("esalgo:RM1299"));
#ifdef EDM_ML_DEBUG
            edm::LogVerbatim("SFGeomX") << ddname << " Union " << solid_c.name() << ":" << solid_lbck.name() << " at "
                                        << DDTranslation(sdxe3[enb], sdye3[enb], sdze3[enb])
                                        << " rotation esalgo:RM1299";
#endif

            DDLogicalPart ladder = DDLogicalPart(ddname, getLaddMaterial(), solid);
            DDName ddname2(getLadPrefix(1) + types_l4_[d], "esalgo");
            DDLogicalPart ladder2 = DDLogicalPart(ddname2, getLaddMaterial(), solid);
          }
        }

      }  //end of not plain ladder shape
      else {
        DDName dd_tmp_name_pa(getLadPrefix(2) + "p", "esalgo");
        DDName dd_tmp_name_pb(getLadPrefix(3) + "p", "esalgo");

        boxay = ladder_length - LdrFrnt_Length - LdrBck_Length;
        boxax = ladder_width;
        boxaz = ladder_thick;

        DDSolid solid_pa = DDSolidFactory::box(dd_tmp_name_pa, boxax / 2, boxay / 2, boxaz / 2.);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("SFGeomX") << dd_tmp_name_pa << " Box " << (boxax / 2) << ":" << (boxay / 2) << ":"
                                    << (boxaz / 2);
#endif
        sdx = 0;
        sdy = -boxay / 2 - LdrFrnt_Length / 2;
        sdz = -ladder_thick / 2. + LdrFrnt_Offset;

        DDSolid solid_pb = DDSolidFactory::unionSolid(
            dd_tmp_name_pb, solid_pa, solid_lfront, DDTranslation(sdx, sdy, sdz), DDRotation("esalgo:RM1299"));
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("SFGeomX") << dd_tmp_name_pb << " Union " << solid_pa.name() << ":" << solid_lfront.name()
                                    << " at " << DDTranslation(sdx, sdy, sdz) << " rotation esalgo:RM1299";
#endif

        sdx = 0;
        sdy = boxay / 2 + LdrBck_Length / 2;
        sdz = -ladder_thick / 2. + LdrBck_Offset;
        DDSolid solid = DDSolidFactory::unionSolid(
            ddname, solid_pb, solid_lbck, DDTranslation(sdx, sdy, sdz), DDRotation("esalgo:RM1299"));
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("SFGeomX") << ddname << " Union " << solid_pb.name() << ":" << solid_lbck.name() << " at "
                                    << DDTranslation(sdx, sdy, sdz) << " rotation esalgo:RM1299";
#endif
        DDLogicalPart ladder = DDLogicalPart(ddname, getLaddMaterial(), solid);
        DDName ddname2(getLadPrefix(1) + types_l4_[d], "esalgo");
        DDLogicalPart ladder2 = DDLogicalPart(ddname2, getLaddMaterial(), solid);
      }
    }

    // insert SWED, SFBX and SFBY into ladders
    swed_scopy_glob++;
    if (M < int(types_l5_.size())) {
      DDName ddname(getLadPrefix(0) + types_l5_[M], "esalgo");
      DDName ddname2(getLadPrefix(1) + types_l5_[M], "esalgo");
      for (int i = 0; i <= 1; i++) {
        for (int j = 0; j <= 4; j++) {
          xpos = (i * 2 - 1) * waf_intra_col_sep / 2.;
          ypos = -ladder_length / 2. + 0.05 * mm - (LdrFrnt_Length - LdrBck_Length) / 2 + wedge_length / 2. +
                 j * waf_active;
          zpos = -ladder_thick / 2. + 0.005 * mm + wedge_offset;
          if (ladd_l5_map_[(i + j * 2 + M * 10)] == 1) {
            scopy++;
            cpv.position(DDLogicalPart("esalgo:SWED"),
                         ddname,
                         scopy + 1000 * swed_scopy_glob,
                         DDTranslation(xpos, ypos, zpos),
                         DDRotation("esalgo:RM1299"));
            cpv.position(DDLogicalPart("esalgo:SWED"),
                         ddname2,
                         scopy + 1000 * swed_scopy_glob + 100,
                         DDTranslation(xpos, ypos, zpos),
                         DDRotation("esalgo:RM1299"));
#ifdef EDM_ML_DEBUG
            edm::LogVerbatim("SFGeom") << DDLogicalPart("esalgo:SWED").name() << " copy "
                                       << (scopy + 1000 * swed_scopy_glob) << " in " << ddname << " at "
                                       << DDTranslation(xpos, ypos, zpos) << " rotation esalgo:RM1299";
            edm::LogVerbatim("SFGeom") << DDLogicalPart("esalgo:SWED").name() << " copy "
                                       << (scopy + 1000 * swed_scopy_glob + 100) << " in " << ddname2 << " at "
                                       << DDTranslation(xpos, ypos, zpos) << " rotation esalgo:RM1299";
#endif
            ypos = ypos + ywedge_ceramic_diff;
            zpos = -ladder_thick / 2. + 0.005 * mm + zwedge_ceramic_diff;
            cpv.position(DDLogicalPart("esalgo:SFBX"),
                         ddname,
                         scopy + 1000 * swed_scopy_glob,
                         DDTranslation(xpos, ypos, zpos),
                         DDRotation("esalgo:RM1298"));
            cpv.position(DDLogicalPart("esalgo:SFBY"),
                         ddname2,
                         scopy + 1000 * swed_scopy_glob,
                         DDTranslation(xpos, ypos, zpos),
                         DDRotation("esalgo:RM1300A"));
#ifdef EDM_ML_DEBUG
            edm::LogVerbatim("SFGeom") << DDLogicalPart("esalgo:SFBX").name() << " copy "
                                       << (scopy + 1000 * swed_scopy_glob) << " in " << ddname << " at "
                                       << DDTranslation(xpos, ypos, zpos) << " rotation esalgo:RM1298";
            edm::LogVerbatim("SFGeom") << DDLogicalPart("esalgo:SFBY").name() << " copy "
                                       << (scopy + 1000 * swed_scopy_glob) << " in " << ddname2 << " at "
                                       << DDTranslation(xpos, ypos, zpos) << " rotation esalgo:RM1300A";
#endif
          }
        }
      }
    } else {
      int d = M - types_l5_.size();
      DDName ddname(getLadPrefix(0) + types_l4_[d], "esalgo");
      DDName ddname2(getLadPrefix(1) + types_l4_[d], "esalgo");
      for (int i = 0; i <= 1; i++) {
        for (int j = 0; j <= 3; j++) {
          xpos = (i * 2 - 1) * waf_intra_col_sep / 2.;
          ypos = -ladder_length / 2. + 0.05 * mm - (LdrFrnt_Length - LdrBck_Length) / 2 + wedge_length / 2. +
                 j * waf_active;
          zpos = -ladder_thick / 2. + 0.005 * mm + wedge_offset;
          if (ladd_l4_map_[(i + j * 2 + (M - types_l5_.size()) * 8)] == 1) {
            scopy++;
            cpv.position(DDLogicalPart("esalgo:SWED"),
                         ddname,
                         scopy + 1000 * swed_scopy_glob,
                         DDTranslation(xpos, ypos, zpos),
                         DDRotation("esalgo:RM1299"));
            cpv.position(DDLogicalPart("esalgo:SWED"),
                         ddname2,
                         scopy + 1000 * swed_scopy_glob + 100,
                         DDTranslation(xpos, ypos, zpos),
                         DDRotation("esalgo:RM1299"));
#ifdef EDM_ML_DEBUG
            edm::LogVerbatim("SFGeom") << DDLogicalPart("esalgo:SWED").name() << " copy "
                                       << (scopy + 1000 * swed_scopy_glob) << " in " << ddname << " at "
                                       << DDTranslation(xpos, ypos, zpos) << " rotation esalgo:RM1299";
            edm::LogVerbatim("SFGeom") << DDLogicalPart("esalgo:SWED").name() << " copy "
                                       << (scopy + 1000 * swed_scopy_glob + 100) << " in " << ddname2 << " at "
                                       << DDTranslation(xpos, ypos, zpos) << " rotation esalgo:RM1299";
#endif
            ypos = ypos + ywedge_ceramic_diff;
            zpos = -ladder_thick / 2. + 0.005 * mm + zwedge_ceramic_diff;
            cpv.position(DDLogicalPart("esalgo:SFBX"),
                         ddname,
                         scopy + 1000 * swed_scopy_glob,
                         DDTranslation(xpos, ypos, zpos),
                         DDRotation("esalgo:RM1298"));
            cpv.position(DDLogicalPart("esalgo:SFBY"),
                         ddname2,
                         scopy + 1000 * swed_scopy_glob,
                         DDTranslation(xpos, ypos, zpos),
                         DDRotation("esalgo:RM1300A"));
#ifdef EDM_ML_DEBUG
            edm::LogVerbatim("SFGeom") << DDLogicalPart("esalgo:SFBX").name() << " copy "
                                       << (scopy + 1000 * swed_scopy_glob) << " in " << ddname << " at "
                                       << DDTranslation(xpos, ypos, zpos) << " rotation esalgo:RM1298";
            edm::LogVerbatim("SFGeom") << DDLogicalPart("esalgo:SFBY").name() << " copy "
                                       << (scopy + 1000 * swed_scopy_glob) << " in " << ddname2 << " at "
                                       << DDTranslation(xpos, ypos, zpos) << " rotation esalgo:RM1300A";
#endif
          }
        }
      }
    }
  }

  // Positioning of ladders

  std::string type;
  int icopy[100] = {0};

  for (int I = -9; I <= 9; ++I) {
    prev_length_ = 0;
    int J = std::abs(I);
    for (int K = 0; K < noLaddInCol_[J]; K++) {
      std::string type;

      ladder_new_length_ = micromodule_length + 3 * waf_active;
      ladd_shift_ = 4 * waf_active;

      if (K == 0)
        type = typeOfLaddRow0[J];
      if (K == 1)
        type = typeOfLaddRow1[J];
      if (K == 2)
        type = typeOfLaddRow2[J];
      if (K == 3)
        type = typeOfLaddRow3[J];

      for (const auto& i : types_l5_)
        if (type == i) {
          ladder_new_length_ = micromodule_length + 4 * waf_active;
          ladd_shift_ = 5 * waf_active;
        }

      int j = 0;

      for (int t = 0; t < int(types_l5_.size()); t++)
        if (type == types_l5_[t]) {
          j = t;
          if (I < 0 && asym_ladd_[t] == 1) {
            j = j + 1;
            type = types_l5_[j];
          }
        }
      for (int t = 0; t < int(types_l4_.size()); t++)
        if (type == types_l4_[t]) {
          j = t + types_l5_.size();
          if (I < 0 && asym_ladd_[(t + types_l5_.size())] == 1) {
            j = j + 1;
            type = types_l4_[j - types_l5_.size()];
          }
        }

      xpos = I * (2 * waf_intra_col_sep + waf_inter_col_sep);
      if (I > 0)
        xpos = xpos + dee_separation;
      if (I < 0)
        xpos = xpos - dee_separation;

      int sz = 20;
      ypos = (sz - int(startOfFirstLadd_[J])) * waf_active - ladder_new_length_ / 2. +
             (LdrFrnt_Length - LdrBck_Length) / 2 + micromodule_length + 0.05 * cm - prev_length_;

      prev_length_ += ladd_shift_;

      zpos = zlead1_ + ladder_thick / 2. + 0.01 * mm;
      icopy[j] += 1;
      DDName ddname(getLadPrefix(0) + type, "esalgo");

      cpv.position(
          DDLogicalPart(ddname), DDName("SF", "esalgo"), icopy[j], DDTranslation(xpos, ypos, zpos), DDRotation());
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("SFGeom") << DDLogicalPart(ddname).name() << " copy " << icopy[j] << " in "
                                 << DDName("SF", "esalgo") << " at " << DDTranslation(xpos, ypos, zpos)
                                 << " no rotation";
#endif
      DDName ddname2(getLadPrefix(1) + type, "esalgo");

      xpos = I * (2 * waf_intra_col_sep + waf_inter_col_sep);

      cpv.position(DDLogicalPart(ddname2),
                   DDName("SF", "esalgo"),
                   icopy[j],
                   DDTranslation(ypos, -xpos, zpos - zlead1_ + zlead2_),
                   DDRotation("esalgo:R270"));
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("SFGeom") << DDLogicalPart(ddname2).name() << " copy " << icopy[j] << " in "
                                 << DDName("SF", "esalgo") << " at "
                                 << DDTranslation(ypos, -xpos, zpos - zlead1_ + zlead2_) << " rotation esalgo:R270";
#endif
      int changed = 0;
      for (int t = 0; t < int(types_l5_.size()); t++)
        if (type == types_l5_[t]) {
          j = t;
          if (asym_ladd_[t] == 2 && !changed) {
            j = j - 1;
            changed = 1;
          }
          if (asym_ladd_[t] == 1 && !changed) {
            j = j + 1;
            changed = 1;
          }
          type = types_l5_[j];
        }
      for (int t = 0; t < int(types_l4_.size()); t++)
        if (type == types_l4_[t]) {
          j = t + types_l5_.size();
          if (asym_ladd_[(t + types_l5_.size())] == 2 && !changed) {
            j = j - 1;
            changed = 1;
          }
          if (asym_ladd_[(t + types_l5_.size())] == 1 && !changed) {
            j = j + 1;
            changed = 1;
          }
          type = types_l4_[j - types_l5_.size()];
        }

      icopy[j] += 1;

      if (I > 0)
        xpos = xpos + dee_separation;
      if (I < 0)
        xpos = xpos - dee_separation;

      DDName ddname3(getLadPrefix(0) + type, "esalgo");
      cpv.position(DDLogicalPart(ddname3),
                   DDName("SF", "esalgo"),
                   icopy[j],
                   DDTranslation(xpos, -ypos, zpos),
                   DDRotation("esalgo:R180"));
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("SFGeom") << DDLogicalPart(ddname3).name() << " copy " << icopy[j] << " in "
                                 << DDName("SF", "esalgo") << " at " << DDTranslation(xpos, -ypos, zpos)
                                 << " rotation esalgo:R180";
#endif
      DDName ddname4(getLadPrefix(1) + type, "esalgo");

      xpos = I * (2 * waf_intra_col_sep + waf_inter_col_sep);

      cpv.position(DDLogicalPart(ddname4),
                   DDName("SF", "esalgo"),
                   icopy[j],
                   DDTranslation(-ypos, -xpos, zpos - zlead1_ + zlead2_),
                   DDRotation("esalgo:R090"));
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("SFGeom") << DDLogicalPart(ddname4).name() << " copy " << icopy[j] << " in "
                                 << DDName("SF", "esalgo") << " at "
                                 << DDTranslation(-ypos, -xpos, zpos - zlead1_ + zlead2_) << " rotation esalgo:R090";
#endif
    }
  }
}

void DDEcalPreshowerAlgo::doSens(DDCompactView& cpv) {
  double xpos(0), ypos(0);
  for (size_t i = 0; i < 32; ++i) {
    xpos = -waf_active / 2. + i * waf_active / 32. + waf_active / 64.;
    cpv.position(
        DDLogicalPart("esalgo:SFSX"), DDName("SFWX", "esalgo"), i + 1, DDTranslation(xpos, 0., 0.), DDRotation());
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("SFGeom") << DDLogicalPart("esalgo:SFSX").name() << " copy " << (i + 1) << " in "
                               << DDName("SFWX", "esalgo") << " at " << DDTranslation(xpos, 0., 0.) << " no rotation";
#endif

    ypos = -waf_active / 2. + i * waf_active / 32. + waf_active / 64.;
    cpv.position(
        DDLogicalPart("esalgo:SFSY"), DDName("SFWY", "esalgo"), i + 1, DDTranslation(0., ypos, 0.), DDRotation());
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("SFGeom") << DDLogicalPart("esalgo:SFSY").name() << " copy " << (i + 1) << " in "
                               << DDName("SFWY", "esalgo") << " at " << DDTranslation(0., ypos, 0.) << " no rotation";
#endif
  }
}

#include "DetectorDescription/Core/interface/DDAlgorithmFactory.h"

DEFINE_EDM_PLUGIN(DDAlgorithmFactory, DDEcalPreshowerAlgo, "ecal:DDEcalPreshowerAlgo");