DDBHMAngular.cc

Go to the documentation of this file.

// File: DDBHMAngular.cc
// Description: Position inside the mother according to phi

#include <cmath>
#include <algorithm>

#include "DDBHMAngular.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DetectorDescription/Core/interface/DDutils.h"
#include "DetectorDescription/Core/interface/DDLogicalPart.h"
#include "DetectorDescription/Core/interface/DDCurrentNamespace.h"
#include "DetectorDescription/Core/interface/DDSplit.h"
#include "CLHEP/Units/GlobalSystemOfUnits.h"

DDBHMAngular::DDBHMAngular() {
  LogDebug("BHMGeom") << "DDBHMAngular test: Creating an instance";
}

DDBHMAngular::~DDBHMAngular() {}

void DDBHMAngular::initialize(const DDNumericArguments & nArgs,
                  const DDVectorArguments & ,
                  const DDMapArguments & ,
                  const DDStringArguments & sArgs,
                  const DDStringVectorArguments & ) {

  units  = int (nArgs["number"]);
  rr     = nArgs["radius"];
  dphi   = nArgs["deltaPhi"];
  LogDebug("BHMGeom") << "DDBHMAngular debug: Parameters for positioning-- "
               << units << " copies at radius " << rr/CLHEP::cm 
               << " cm with delta(phi) " << dphi/CLHEP::deg;

  rotMat      = sArgs["Rotation"];
  childName   = sArgs["ChildName"]; 
  LogDebug("BHMGeom") << "DDBHMAngular debug: Parent " << parent().name()
               << "\tChild " << childName << "\tRotation matrix " 
               << rotMat;
}

void DDBHMAngular::execute(DDCompactView& cpv) {

  DDName child(DDSplit(childName).first, DDSplit(childName).second);
  DDName parentName  = parent().name(); 
  std::string rotstr = DDSplit(rotMat).first;
  DDRotation rot;
  if (rotstr != "NULL") {
    std::string rotns  = DDSplit(rotMat).second;
    rot = DDRotation(DDName(rotstr, rotns));
  }

  for (int jj=0; jj<units; jj++) {
    double driverX(0), driverY(0), driverZ(0);
    if (jj<16){                             
      driverX = rr*cos((jj+0.5)*dphi);
      driverY = sqrt(rr*rr-driverX*driverX);
    } else if (jj==16) { 
      driverX = rr*cos(15.5*dphi);
      driverY =-sqrt(rr*rr-driverX*driverX);                      
    } else if (jj==17) { 
      driverX = rr*cos(14.5*dphi);
      driverY =-sqrt(rr*rr-driverX*driverX);                      
    } else if (jj==18) { 
      driverX = rr*cos(0.5*dphi);
      driverY =-sqrt(rr*rr-driverX*driverX);                      
    } else if (jj==19) { 
      driverX = rr*cos(1.5*dphi);
      driverY =-sqrt(rr*rr-driverX*driverX); 
    }               
    DDTranslation tran(driverX, driverY, driverZ);
  
    cpv.position(child, parentName, jj+1, tran, rot);
    LogDebug("BHMGeom") << "DDBHMAngular test: " << child << " number " <<jj+1
             << " positioned in " << parentName << " at " << tran 
             << " with " << rot;
  }
}