DDDividedTrdY Class Reference

#include <DDDividedTrd.h>

Inheritance diagram for DDDividedTrdY:

Public Member Functions
virtual void	checkParametersValidity () override
	DDDividedTrdY (const DDDivision &div, DDCompactView *cpv)
virtual double	getMaxParameter () const override
virtual DDLogicalPart	makeDDLogicalPart (const int copyNo) const override
virtual DDRotation	makeDDRotation (const int copyNo) const override
virtual DDTranslation	makeDDTranslation (const int copyNo) const override
Public Member Functions inherited from DDDividedGeometryObject
	DDDividedGeometryObject (const DDDivision &div, DDCompactView *cpv)
virtual void	execute (void)
virtual const std::string &	getType (void) const
virtual void	setType (const std::string &type)
int	volumeFirstCopyNo (void) const
virtual	~DDDividedGeometryObject (void)=default

Additional Inherited Members
Static Public Member Functions inherited from DDDividedGeometryObject
static const double	tolerance (void)
Protected Member Functions inherited from DDDividedGeometryObject
int	calculateNDiv (double motherDim, double width, double offset) const
double	calculateWidth (double motherDim, int nDiv, double offset) const
DDRotationMatrix *	changeRotMatrix (double rotZ=0.) const
void	checkNDivAndWidth (double maxPar)
void	checkOffset (double maxPar)
Protected Attributes inherited from DDDividedGeometryObject
int	compNDiv_
double	compWidth_
DDCompactView *	cpv_
DDDivision	div_
DivisionType	divisionType_
std::string	ftype_
int	theVoluFirstCopyNo_

Detailed Description

Definition at line 28 of file DDDividedTrd.h.

Constructor & Destructor Documentation

DDDividedTrdY::DDDividedTrdY	(	const DDDivision &	div,
		DDCompactView *	cpv
	)

Definition at line 168 of file DDDividedTrd.cc.

References DDDividedGeometryObject::calculateNDiv(), DDDividedGeometryObject::calculateWidth(), checkParametersValidity(), DDDividedGeometryObject::compNDiv_, DDDividedGeometryObject::compWidth_, DDDividedGeometryObject::div_, DDDividedGeometryObject::divisionType_, DivNDIV, DivWIDTH, DDDivision::nReplicas(), DDDivision::offset(), DDDivision::parent(), DDDividedGeometryObject::setType(), DDLogicalPart::solid(), DDDivision::width(), and DDTrap::y1().

  : DDDividedGeometryObject( div, cpv )
{
  checkParametersValidity();
  setType( "DivisionTrdY" );
  DDTrap mtrd = (DDTrap)(div_.parent().solid());

  if( divisionType_ == DivWIDTH )
  {
    compNDiv_ = calculateNDiv( 2 * mtrd.y1(), div_.width(), div_.offset() );
  }
  else if( divisionType_ == DivNDIV )
  {
    compWidth_ = calculateWidth( 2 * mtrd.y1(), div_.nReplicas(), div_.offset() );
  }
}

Member Function Documentation

void DDDividedTrdY::checkParametersValidity ( void  )

overridevirtual

Reimplemented from DDDividedGeometryObject.

Definition at line 263 of file DDDividedTrd.cc.

References DDTrap::alpha1(), DDTrap::alpha2(), DDDividedGeometryObject::checkParametersValidity(), DDDividedGeometryObject::div_, Exception, DDDivision::parent(), DDTrap::phi(), alignCSCRings::s, DDLogicalPart::solid(), AlCaHLTBitMon_QueryRunRegistry::string, DDTrap::theta(), DDDividedGeometryObject::tolerance(), DDTrap::y1(), and DDTrap::y2().

Referenced by DDDividedTrdY().

{
  DDDividedGeometryObject::checkParametersValidity();

  DDTrap mtrd = (DDTrap)(div_.parent().solid());

  double mpDy1 = mtrd.y1(); //->GetYHalfLength1();
  double mpDy2 = mtrd.y2(); //->GetYHalfLength2();
  double mpTheta = mtrd.theta();
  double mpPhi = mtrd.phi();
  double mpAlpha1 = mtrd.alpha1();
  double mpAlpha2 = mtrd.alpha2();

  if( fabs(mpDy1 - mpDy2) > tolerance() )
  {
    std::string s= "ERROR - DDDividedTrdY::checkParametersValidity()";
    s += "\n        Making a division of a TRD along axis Y while";
    s += "\n        the Y half lengths are not equal is not (yet)";
    s += "\n        supported. It will result in non-equal";
    s += "\n        division solids.";
    throw cms::Exception("DDException") << s;
  }
  // mec:  we only have traps, not trds in DDD, so I added this check
  // to make sure it is only a trd (I think! :-))
  if (mpAlpha1 != 0.*deg || mpAlpha2 != 0.*deg || mpTheta != 0.*deg || mpPhi != 0.*deg)
  {
    std::string s = "ERROR - DDDividedTrdY::checkParametersValidity()";
    s+= "\n        Making a division of a TRD along axis X,";
    s+= "\n        while the theta, phi and aplhpa2 are not zero,";
    s+= "\n        is not (yet) supported. It will result";
    s+= "\n        in non-equal division solids.";
    throw cms::Exception("DDException") << s;
  }
}

double DDDividedTrdY::getMaxParameter ( void  ) const

overridevirtual

Reimplemented from DDDividedGeometryObject.

Definition at line 186 of file DDDividedTrd.cc.

References DDDividedGeometryObject::div_, DDDivision::parent(), DDLogicalPart::solid(), and DDTrap::y1().

{
  DDTrap mtrd = (DDTrap)(div_.parent().solid());
  return 2 * mtrd.y1(); 
}

DDLogicalPart DDDividedTrdY::makeDDLogicalPart ( const int  copyNo ) const

overridevirtual

Reimplemented from DDDividedGeometryObject.

Definition at line 225 of file DDDividedTrd.cc.

References DDDividedGeometryObject::compWidth_, DDDividedGeometryObject::div_, DDTrap::halfZ(), DDLogicalPart::material(), DDBase< N, C >::name(), DDDivision::parent(), DDLogicalPart::solid(), DDSolidFactory::trap(), DDTrap::x1(), DDTrap::x2(), DDTrap::x3(), and DDTrap::x4().

{
  //---- The division along Y of a Trd will result a Trd, only 
  //--- if Y at -Z and +Z are equal, else use the G4Trap version
  DDTrap mtrd = (DDTrap)(div_.parent().solid());
  DDMaterial usemat = div_.parent().material();
  
  double pDx1 = mtrd.x1(); //->GetXHalfLength1() at Y+;
  double pDx2 = mtrd.x2(); //->GetXHalfLength2() at Y+;
  double pDx3 = mtrd.x3(); //->GetXHalfLength1() at Y-;
  double pDx4 = mtrd.x4(); //->GetXHalfLength2() at Y-;
  double pDz = mtrd.halfZ(); //->GetZHalfLength();
  double pDy = compWidth_/2.;
 
  //trd.SetAllParameters ( pDx1, pDx2, pDy, pDy, pDz );
  DDName solname(div_.name() );
  DDSolid  dsol(solname);
  DDLogicalPart ddlp(solname);
  if (!dsol.isDefined().second)
  {
    dsol = DDSolidFactory::trap(solname
                , pDz
                , 0.*deg
                , 0.*deg
                , pDy
                , pDx1
                , pDx2
                , 0*deg
                , pDy
                , pDx3
                , pDx4
                , 0.*deg);
    DDLogicalPart ddlp(solname,  usemat, dsol);
  }
  return ddlp;
}

DDRotation DDDividedTrdY::makeDDRotation ( const int  copyNo ) const

overridevirtual

Reimplemented from DDDividedGeometryObject.

Definition at line 219 of file DDDividedTrd.cc.

{
  return DDRotation();
}

DDTranslation DDDividedTrdY::makeDDTranslation ( const int  copyNo ) const

overridevirtual

Reimplemented from DDDividedGeometryObject.

Definition at line 193 of file DDDividedTrd.cc.

References DDDivision::axis(), DDDividedGeometryObject::compWidth_, DDDividedGeometryObject::div_, Exception, DDAxesNames::name(), DDDivision::offset(), DDDivision::parent(), alignCSCRings::s, DDLogicalPart::solid(), AlCaHLTBitMon_QueryRunRegistry::string, y, and DDTrap::y1().

{
  DDTrap mtrd = (DDTrap)(div_.parent().solid() );
  double mdy = mtrd.y1();

  //----- translation 
  double posi = -mdy + div_.offset() + (copyNo+0.5)*compWidth_;

  if( div_.axis() == DDAxes::y )
  {
    return DDTranslation(0.0, posi, 0.0);
  }
  else
  { 
    std::string s = "ERROR - DDDividedTrdY::makeDDTranslation()";
    s += "\n        Axis is along ";
    s += DDAxesNames::name(div_.axis());
    s += " !\n" ;
    s += "DDDividedTrdY::makeDDTranslation()";
    s += " IllegalConstruct: Only axes along y are allowed !";
    throw cms::Exception("DDException") << s;
  }
  return DDTranslation();
}