#include <DDCutTubsFromPoints.h>

Inheritance diagram for DDCutTubsFromPoints:

Classes
struct	Section

Public Member Functions
	DDCutTubsFromPoints ()

void	execute (DDCompactView &cpv) override

void	initialize (const DDNumericArguments &nArgs, const DDVectorArguments &vArgs, const DDMapArguments &mArgs, const DDStringArguments &sArgs, const DDStringVectorArguments &vsArgs) override

	~DDCutTubsFromPoints () override

Private Attributes
DDMaterial	material

double	r_max

double	r_min

std::vector< Section >	sections

DDName	solidOutput

double	z_pos

Detailed Description

Definition at line 14 of file DDCutTubsFromPoints.h.

Constructor & Destructor Documentation

DDCutTubsFromPoints::DDCutTubsFromPoints ( )

Definition at line 18 of file DDCutTubsFromPoints.cc.

References LogDebug.

                                          {
   LogDebug("TrackerGeom") <<"DDCutTubsFromPoints info: Creating an instance";
 }

DDCutTubsFromPoints::~DDCutTubsFromPoints ( )

override

Definition at line 22 of file DDCutTubsFromPoints.cc.

22 {}

Member Function Documentation

void DDCutTubsFromPoints::execute ( DDCompactView & cpv )

override

Definition at line 66 of file DDCutTubsFromPoints.cc.

References funct::cos(), DDSolidFactory::cuttubs(), PVValHelper::dz, mps_fire::i, material, DDName::name(), DDName::ns(), PFRecoTauDiscriminationByIsolation_cfi::offset, class-composition::parent, DDCutTubsFromPoints::Section::phi, DDCompactView::position(), alignCSCRings::r, r_max, r_min, indexGen::s2, sections, edm::shift, funct::sin(), solidOutput, mathSSE::sqrt(), square(), DDEnums::support, DDSolidFactory::unionSolid(), DDCutTubsFromPoints::Section::z_l, z_pos, and DDCutTubsFromPoints::Section::z_t.

                                                     {
 
   // radius for plane calculations
   // We use r_max here, since P3 later has a Z that is always more inside
   // than the extreme points. This means the cutting planes have outwards
   // slopes in r-Z, and the corner at r_max could stick out of the bounding
   // volume otherwise. 
   double r  = r_max;
 
   // min and max z for the placement in the end
   double min_z =  1e9;
   double max_z = -1e9;
 
   // counter of actually produced segments (excluding skipped ones)
   int segment = 0;
 
   // the segments and their corresponding offset (absolute, as in the input)
   std::vector<DDSolid> segments;
   std::vector<double>  offsets;
   
   Section s1 = sections[0];
   for (Section s2 : sections) {
     if (s1.phi != s2.phi) {
       segment++;
       // produce segment s1-s2.
       DDName segname  (solidOutput.name() + "_seg_" + std::to_string(segment), 
                        solidOutput.ns());
 
       double phi1 = s1.phi;
       double phi2 = s2.phi;
 
       // track the min/max to properly place&align later
       if (s2.z_l < min_z) min_z = s2.z_l;
       if (s2.z_t > max_z) max_z = s2.z_t;
 
       double P1_z_l = s1.z_l;
       double P1_z_t = s1.z_t;
       double P2_z_l = s2.z_l;
       double P2_z_t = s2.z_t;
 
       double P1_x_t = cos(phi1) * r;
       double P1_x_l = cos(phi1) * r;
       double P1_y_t = sin(phi1) * r;
       double P1_y_l = sin(phi1) * r;
 
       double P2_x_t = cos(phi2) * r;
       double P2_x_l = cos(phi2) * r;
       double P2_y_t = sin(phi2) * r;
       double P2_y_l = sin(phi2) * r;
 
       // each cutting plane is defined by P1-3. P1-2 are corners of the
       // segment, P3 is at r=0 with the "average" z to get a nice cut.
       double P3_z_l = (P1_z_l + P2_z_l) / 2;
       double P3_z_t = (P1_z_t + P2_z_t) / 2;
 
       // we only have one dz to position both planes. The anchor is implicitly
       // between the P3's, we have to use an offset later to make the segments
       // line up correctly.
       double dz = (P3_z_t - P3_z_l) / 2;
       double offset = (P3_z_t + P3_z_l) / 2;
 
       // the plane is defined by P1-P3 and P2-P3; since P3 is at r=0 we
       // only need the z.
       double D1_z_l = P1_z_l - P3_z_l;
       double D2_z_l = P2_z_l - P3_z_l;
 
       // the normal is then the cross product...
       double n_x_l = (P1_y_l * D2_z_l) - (D1_z_l * P2_y_l);
       double n_y_l = (D1_z_l * P2_x_l) - (P1_x_l * D2_z_l);
       double n_z_l = (P1_x_l * P2_y_l) - (P1_y_l * P2_x_l);
 
       // ... normalized.
       // flip the sign here (but not for t) since root wants it like that.
       double norm = -sqrt(square(n_x_l) + square(n_y_l) + square(n_z_l));
       n_x_l /= norm;
       n_y_l /= norm;
       n_z_l /= norm;
 
       // same game for the t side.
       double D1_z_t = P1_z_t - P3_z_t;
       double D2_z_t = P2_z_t - P3_z_t;
 
       double n_x_t = (P1_y_t * D2_z_t) - (D1_z_t * P2_y_t);
       double n_y_t = (D1_z_t * P2_x_t) - (P1_x_t * D2_z_t);
       double n_z_t = (P1_x_t * P2_y_t) - (P1_y_t * P2_x_t);
 
       norm = sqrt(square(n_x_t) + square(n_y_t) + square(n_z_t));
       n_x_t /= norm;
       n_y_t /= norm;
       n_z_t /= norm;
 
       // the cuttubs wants a delta phi
       double dphi = phi2 - phi1;
 
       DDSolid seg = DDSolidFactory::cuttubs(segname, dz, r_min, r_max, phi1, dphi,
                                             n_x_l, n_y_l, n_z_l,
                                             n_x_t, n_y_t, n_z_t); 
       segments.emplace_back(seg);
       offsets.emplace_back(offset);
     }
     s1 = s2;
   }
 
   assert(segments.size() >= 2); // less would be special cases
 
   DDSolid solid = segments[0];
   // placment happens relative to the first member of the union
   double shift = offsets[0];
 
   for (unsigned i = 1; i < segments.size()-1; i++) {
     // each sub-union needs a name. Well. 
     DDName unionname(solidOutput.name() + "_uni_" + std::to_string(i+1), 
                      solidOutput.ns());
     solid = DDSolidFactory::unionSolid(unionname, solid, segments[i], 
              DDTranslation(0, 0, offsets[i] - shift), 
              DDRotation());
   }
 
   solid = DDSolidFactory::unionSolid(solidOutput, solid, segments[segments.size()-1], 
             DDTranslation(0, 0, offsets[segments.size()-1] - shift),
             DDRotation());
 
   // remove the common offset from the input, to get sth. aligned at z=0.
   double offset = - shift + (min_z + (max_z-min_z)/2);
 
   DDLogicalPart logical(solidOutput, material, solid, DDEnums::support);
 
   // This is not as generic as the name promises, but it is hard to make a 
   // solid w/o placing it.
   DDRotation rot180("pixfwdCommon:Z180");
   cpv.position(logical, parent(), 1, DDTranslation(0, 0, z_pos - offset), DDRotation());
   cpv.position(logical, parent(), 2, DDTranslation(0, 0, z_pos - offset), rot180);
 
 }

void DDCutTubsFromPoints::initialize	(	const DDNumericArguments &	nArgs,
		const DDVectorArguments &	vArgs,
		const DDMapArguments &	mArgs,
		const DDStringArguments &	sArgs,
		const DDStringVectorArguments &	vsArgs
	)

override

Definition at line 24 of file DDCutTubsFromPoints.cc.

References mps_fire::i, LogDebug, material, DDCurrentNamespace::ns(), class-composition::parent, r_max, r_min, alignCSCRings::s, sections, DDVector::size(), solidOutput, AlCaHLTBitMon_QueryRunRegistry::string, and z_pos.

                                                       {
 
   r_min = nArgs["rMin"];
   r_max = nArgs["rMax"];
   z_pos = nArgs["zPos"];
   material = DDMaterial(sArgs["Material"]);
   
   auto phis_name = DDName(sArgs["Phi"]);
   auto z_ls_name = DDName(sArgs["z_l"]);
   auto z_ts_name = DDName(sArgs["z_t"]);
   DDVector phis(phis_name);
   DDVector z_ls(z_ls_name);
   DDVector z_ts(z_ts_name);
 
   assert(phis.size() == z_ls.size());
   assert(phis.size() == z_ts.size());
 
   for (unsigned i = 0; i < phis.size(); i++) {
     Section s = {phis[i], z_ls[i], z_ts[i] };
     sections.emplace_back(s);
   }
   assert(!sections.empty());
 
   solidOutput = DDName(sArgs["SolidName"]);
   
   std::string idNameSpace = DDCurrentNamespace::ns();
   DDName parentName = parent().name();
   LogDebug("TrackerGeom") << "DDCutTubsFromPoints debug: Parent " << parentName
                           << "\tSolid " << solidOutput << " NameSpace " 
                           << idNameSpace 
                           << "\tnumber of sections " << sections.size();
 }