GeoPixelIBLFwdSvcCADModel Class Reference

#include <GeoPixelIBLFwdSvcCADModel.h>

Inheritance diagram for GeoPixelIBLFwdSvcCADModel:

Collaboration diagram for GeoPixelIBLFwdSvcCADModel:

Public Member Functions
	GeoPixelIBLFwdSvcCADModel (InDetDD::PixelDetectorManager *ddmgr, PixelGeometryManager *mgr, GeoModelIO::ReadGeoModel *sqliteReader, std::shared_ptr< std::map< std::string, GeoFullPhysVol * > > mapFPV, std::shared_ptr< std::map< std::string, GeoAlignableTransform * > > mapAX, int)
virtual GeoVPhysVol *	Build () override
GeoPhysVol *	getSupportA ()
GeoPhysVol *	getSupportC ()
GeoTransform *	getSupportTrfA ()
GeoTransform *	getSupportTrfC ()

Protected Attributes
PixelGeometryManager *	m_gmt_mgr
InDetMaterialManager *	m_mat_mgr
InDetDD::PixelDetectorManager *	m_DDmgr
GeoModelIO::ReadGeoModel *	m_sqliteReader
std::shared_ptr< std::map< std::string, GeoFullPhysVol * > >	m_mapFPV
std::shared_ptr< std::map< std::string, GeoAlignableTransform * > >	m_mapAX
const double	m_epsilon

Private Member Functions
const GeoShape *	addShape (const GeoShape *lastShape, const GeoShape *nextShape, const GeoTrf::Transform3D &trans)

Private Attributes
GeoPhysVol *	m_supportPhysA
GeoPhysVol *	m_supportPhysC
GeoTransform *	m_xformSupportA
GeoTransform *	m_xformSupportC

Detailed Description

Definition at line 19 of file GeoPixelIBLFwdSvcCADModel.h.

Constructor & Destructor Documentation

GeoPixelIBLFwdSvcCADModel()

GeoPixelIBLFwdSvcCADModel::GeoPixelIBLFwdSvcCADModel	(	InDetDD::PixelDetectorManager *	ddmgr,
		PixelGeometryManager *	mgr,
		GeoModelIO::ReadGeoModel *	sqliteReader,
		std::shared_ptr< std::map< std::string, GeoFullPhysVol * > >	mapFPV,
		std::shared_ptr< std::map< std::string, GeoAlignableTransform * > >	mapAX,
		int	)

Definition at line 36 of file GeoPixelIBLFwdSvcCADModel.cxx.

  : GeoVPixelFactory(ddmgr, mgr, sqliteReader, std::move(mapFPV), std::move(mapAX)), 
  m_supportPhysA(nullptr),
  m_supportPhysC(nullptr),
  m_xformSupportA(nullptr),
  m_xformSupportC(nullptr)
{
}

Member Function Documentation

addShape()

const GeoShape * GeoPixelIBLFwdSvcCADModel::addShape ( const GeoShape * lastShape, const GeoShape * nextShape, const GeoTrf::Transform3D & trans )

private

Definition at line 278 of file GeoPixelIBLFwdSvcCADModel.cxx.

{
  const GeoShape * shiftedShape = &(*nextShape << trans);
  if (lastShape) {
    lastShape = &(lastShape->add(*shiftedShape));
  } else {
    lastShape = shiftedShape;
  }
  return lastShape;
}

Build()

GeoVPhysVol * GeoPixelIBLFwdSvcCADModel::Build ( )

overridevirtual

Implements GeoVPixelFactory.

Definition at line 50 of file GeoPixelIBLFwdSvcCADModel.cxx.

{
 
  m_gmt_mgr->msg(MSG::INFO) <<"Build IBL fwd services - CAD tool design - Build()"<<endmsg;
  
#if !defined(GEOTORUS) || GEOTORUS==0
 
  m_gmt_mgr->msg(MSG::INFO) <<"Build IBL fwd services - CAD tool design - Torus object is not defined "<<GEOTORUS<<endmsg;
  return 0;
 
#else
 
  m_gmt_mgr->msg(MSG::INFO) <<"Build IBL fwd services - CAD tool design - Torus object is defined"<<endmsg;
 
  //  double safety = 0.01*Gaudi::Units::mm;
 
  // IBL layer shift ( 2mm shift issue )
  double layerZshift = m_gmt_mgr->PixelLayerGlobalShift();
  int nSectors = m_gmt_mgr->NPixelSectors();
  double phiOfModuleZero =  m_gmt_mgr->PhiOfModuleZero();  
  double layerRadius = m_gmt_mgr->PixelLayerRadius();
 
  // check if sectors are properly defined
  if(nSectors==0) return nullptr;
  double angle=360./(double)nSectors*Gaudi::Units::deg;
  
  // Defines the IBL_Fwd02 section in the IBL services area
  double innerRadius = 33.;
  double outerRadius = 42.499;
  double zMin = 834.607;
  double zMax = 3290.795;
 
  double zLength = zMax-zMin;
  double halfLength = zLength*.5;
  double zMiddle = (zMin+zMax)*.5;
 
  // Define IBL fwd svc section for side A and C
  std::ostringstream lnameA;
  lnameA<<"Brl0A_FwdSvc";
  std::ostringstream lnameC;
  lnameC<<"Brl0C_FwdSvc";
 
  // Build encompassing volume for both A and C sides (assemblies)
  const GeoTube* supportShapeA = new GeoTube(innerRadius,outerRadius,halfLength);
  const GeoTube* supportShapeC = new GeoTube(innerRadius,outerRadius,halfLength);
  const GeoMaterial* ether = m_mat_mgr->getMaterial("special::Ether");  
  GeoLogVol* supportLogVol_A = new GeoLogVol(lnameA.str(),supportShapeA,ether);
  GeoLogVol* supportLogVol_C = new GeoLogVol(lnameA.str(),supportShapeC,ether);
 
  // Data taken fron CATIA desgin file
  double delta = 18.;
  double deltaJunction = 0.;
  double lgFwdSvc[4]={511., 561., 560., 706. };
  //double devLgFwdSvc[4]={512., 562., 562., 707. };
  double devTotalLength = 2460.188;
  double pi = Gaudi::Units::pi;
 
  // Cable bundle sizes
  double rminCable = 0.;
  double rmaxCable = 8.*.5;
  double surfCable = rmaxCable*rmaxCable*pi;
 
  // Cooling tube sizes
  double rminCooling = 0.;
  double rmaxCooling = 1.75*.5;
  double surfCooling = rmaxCooling*rmaxCooling*pi;
 
  // IPT clip starting position
  double zpos0 = 903.;
  double zpos = 0.;
  double hermJunction = .4;
  double breakAngle = 11.*Gaudi::Units::deg;
 
  //  Loop over the wavy shape sections to build a cable and a cooling pipe
  const GeoShape * gblShapeCableA = nullptr;
  const GeoShape * gblShapeCoolingA = nullptr;
  const GeoShape * gblShapeCableC = nullptr;
  const GeoShape * gblShapeCoolingC = nullptr;
  for(int i=0; i<4; i++)
    {
      double initHalfLength = lgFwdSvc[i]*.5;
      zpos += initHalfLength;
      double corHalfLength = initHalfLength-deltaJunction;
      double radius = .5*(corHalfLength*corHalfLength+delta*delta)/delta;
      double theta = asin(corHalfLength/radius);
 
      if(i==0){
        double tubeLength=deltaJunction+zpos0-zMin+hermJunction;
        double zposTrans=zMin+tubeLength*.5-zpos0;
    const GeoTube* tubeShape = new GeoTube(0., rmaxCable, tubeLength*.5);
    GeoTrf::Transform3D trfA =  GeoTrf::TranslateZ3D(zpos0+zposTrans-zMiddle);
    gblShapeCableA = addShape(gblShapeCableA, tubeShape, trfA );
    GeoTrf::Transform3D trfC =  GeoTrf::RotateY3D(pi)*GeoTrf::TranslateZ3D(zpos0+zposTrans-zMiddle);
    gblShapeCableC = addShape(gblShapeCableC, tubeShape, trfC );
    
    //  std::cout<<"IBL_fwd : junction "<<std::setprecision(15)<<i<<" : "<<zpos0+zposTrans-tubeLength*.5<<" "<<zpos0+zposTrans+tubeLength*.5<<" // "<<zMin<<std::endl;
    
    const GeoTube* tubeShape_cooling = new GeoTube(0., rmaxCooling, tubeLength*.5);
    gblShapeCoolingA = addShape(gblShapeCoolingA, tubeShape_cooling, trfA );
    gblShapeCoolingC = addShape(gblShapeCoolingC, tubeShape_cooling, trfC );
      }
 
      //      std::cout<<"Torus radius : "<<radius<<"   "<<perimeter*(2.*theta)/(2.*pi)+2*deltaJunction<<" / "<<devLgFwdSvc[i]<<std::endl; 
 
      double rtor = radius;
      const GeoTorus* torusShape = new GeoTorus(rminCable,rmaxCable,rtor,pi*.5-theta,2.*theta);
      double angle=(i%2==0)?breakAngle:pi-breakAngle;
      GeoTrf::Transform3D trfA = GeoTrf::RotateZ3D(angle)*GeoTrf::RotateY3D(-pi*.5)*GeoTrf::TranslateY3D(-radius+delta)* GeoTrf::TranslateX3D(zpos0+zpos-zMiddle);
      gblShapeCableA = addShape(gblShapeCableA, torusShape, trfA );
      GeoTrf::Transform3D trfC = GeoTrf::RotateZ3D(angle)*GeoTrf::RotateY3D(pi*.5)*GeoTrf::TranslateY3D(-radius+delta)* GeoTrf::TranslateX3D(zpos0+zpos-zMiddle);
      gblShapeCableC = addShape(gblShapeCableC, torusShape, trfC );
 
      const GeoTorus* torusShape_cooling = new GeoTorus(rminCooling,rmaxCooling,rtor,pi*.5-theta,2.*theta);
      gblShapeCoolingA = addShape(gblShapeCoolingA, torusShape_cooling, trfA );
      gblShapeCoolingC = addShape(gblShapeCoolingC, torusShape_cooling, trfC );
      
      zpos += initHalfLength;
 
      if(i==3){
    zpos=zpos0+zpos-deltaJunction;
    double tubeLength=zMax-zpos+hermJunction;
    double zposTrans=zMax-tubeLength*.5;
    const GeoTube* tubeShape = new GeoTube(0., rmaxCable, tubeLength*.5);
    GeoTrf::Transform3D trfA =  GeoTrf::TranslateZ3D(zposTrans-zMiddle);
    gblShapeCableA = addShape(gblShapeCableA, tubeShape, trfA );
    GeoTrf::Transform3D trfC =  GeoTrf::RotateY3D(pi)*GeoTrf::TranslateZ3D(zposTrans-zMiddle);
    gblShapeCableC = addShape(gblShapeCableC, tubeShape, trfC );
 
    const GeoTube* tubeShape_cooling = new GeoTube(0., rmaxCooling, tubeLength*.5);
    gblShapeCoolingA = addShape(gblShapeCoolingA, tubeShape_cooling, trfA );
    gblShapeCoolingC = addShape(gblShapeCoolingC, tubeShape_cooling, trfC );
      }
 
    }
 
 
  // Material (material budget computed based on the detailed description of a cable bundle and cooling pipe
  const GeoMaterial* cableMat = m_mat_mgr->getMaterialForVolume("pix::IBL_Fwd02_Cable_Wvy",surfCable*devTotalLength);
  GeoLogVol * cable_logA = new GeoLogVol("IBL_Fwd02_Cable_A",gblShapeCableA,cableMat);
  GeoLogVol * cable_logC = new GeoLogVol("IBL_Fwd02_Cable_C",gblShapeCableC,cableMat);
 
  const GeoMaterial* coolingMat = m_mat_mgr->getMaterialForVolume("pix::IBL_Fwd02_Cooling_Wvy",surfCooling*devTotalLength);
  GeoLogVol * cooling_logA = new GeoLogVol("IBL_Fwd02_Cooling_A",gblShapeCoolingA,coolingMat);
  GeoLogVol * cooling_logC = new GeoLogVol("IBL_Fwd02_Cooling_C",gblShapeCoolingC,coolingMat);
 
  // ----------- Create the PhysVol object
  m_supportPhysA = new GeoPhysVol(supportLogVol_A);
  m_supportPhysC = new GeoPhysVol(supportLogVol_C);
 
  // -------------- Alignement of the services with the previous services (Fwd01 area)
  double cooling_radius = 35.1;
  double cooling_angle = -2.154*Gaudi::Units::deg;
 
  if(m_gmt_mgr->PixelStaveAxe()==1)   
    {
      cooling_radius = 34.7 + layerRadius-33.25;
      cooling_angle = -.1*Gaudi::Units::deg;
    }
 
  double cable_radius = 36.501;
  GeoPhysVol* cablePhysVolA = new GeoPhysVol(cable_logA);
  GeoPhysVol* cablePhysVolC = new GeoPhysVol(cable_logC);
  GeoPhysVol* coolingPhysVolA = new GeoPhysVol(cooling_logA);
  GeoPhysVol* coolingPhysVolC = new GeoPhysVol(cooling_logC);
 
  //    nSectors=1;
 
  // Loop over the sectors to place the 14 bundles and cooling pipes
  for(int ii=0; ii<nSectors; ii++)
    {
      m_gmt_mgr->SetPhi(ii);    
      
      // cooling transform 
      double phiOfCooling = phiOfModuleZero+ cooling_angle + ii*angle;
      
      std::ostringstream tmp1; 
      tmp1 << "IBL_Fwd02_Cooling_A" << ii;
      GeoNameTag * tag1 = new GeoNameTag(tmp1.str());
      GeoTransform* xformA1 = new GeoTransform(GeoTrf::RotateZ3D(phiOfCooling)*GeoTrf::TranslateX3D(cooling_radius));
      m_supportPhysA->add(tag1);
      m_supportPhysA->add(xformA1);
      m_supportPhysA->add(coolingPhysVolA);
      
      std::ostringstream tmp1C; 
      tmp1C << "IBL_Fwd02_Cooling_C" << ii;
      GeoNameTag * tag1C = new GeoNameTag(tmp1C.str());
      GeoTransform* xformC1 = new GeoTransform(GeoTrf::RotateZ3D(phiOfCooling)*GeoTrf::TranslateX3D(cooling_radius));
      m_supportPhysC->add(tag1C);
      m_supportPhysC->add(xformC1);
      m_supportPhysC->add(coolingPhysVolC);
      
      double phiOfCable = phiOfCooling + angle*.5;
      
      std::ostringstream tmp2; 
      tmp2 << "IBL_Fwd02_Cable_A" << ii;
      GeoNameTag * tag2 = new GeoNameTag(tmp2.str());
      GeoTransform* xformA2 = new GeoTransform(GeoTrf::RotateZ3D(phiOfCable)*GeoTrf::TranslateX3D(cable_radius));
      m_supportPhysA->add(tag2);
      m_supportPhysA->add(xformA2);
      m_supportPhysA->add(cablePhysVolA);
      
      std::ostringstream tmp2C; 
      tmp2C << "IBL_Fwd02_Cable_C" << ii;
      GeoNameTag * tag2C = new GeoNameTag(tmp2C.str());
      GeoTransform* xformC2 = new GeoTransform(GeoTrf::RotateZ3D(phiOfCable)*GeoTrf::TranslateX3D(cable_radius));
      m_supportPhysC->add(tag2C);
      m_supportPhysC->add(xformC2);
      m_supportPhysC->add(cablePhysVolC);
    
    }
 
  // Global transforms that include the IBL shift in Z
  double middleA = zMiddle+layerZshift;
  double middleC = -zMiddle+layerZshift;
 
  GeoTrf::Transform3D supportTrfA = GeoTrf::TranslateZ3D(middleA);
  m_xformSupportA = new GeoTransform(supportTrfA);
  
  GeoTrf::Transform3D supportTrfC = GeoTrf::TranslateZ3D(middleC);
  m_xformSupportC = new GeoTransform(supportTrfC);
 
  return nullptr;
 
#endif
 
}

getSupportA()

GeoPhysVol * GeoPixelIBLFwdSvcCADModel::getSupportA ( )

inline

Definition at line 31 of file GeoPixelIBLFwdSvcCADModel.h.

{ return m_supportPhysA; }

getSupportC()

GeoPhysVol * GeoPixelIBLFwdSvcCADModel::getSupportC ( )

inline

Definition at line 32 of file GeoPixelIBLFwdSvcCADModel.h.

{ return m_supportPhysC; }

getSupportTrfA()

GeoTransform * GeoPixelIBLFwdSvcCADModel::getSupportTrfA ( )

inline

Definition at line 33 of file GeoPixelIBLFwdSvcCADModel.h.

{ return m_xformSupportA; }

getSupportTrfC()

GeoTransform * GeoPixelIBLFwdSvcCADModel::getSupportTrfC ( )

inline

Definition at line 34 of file GeoPixelIBLFwdSvcCADModel.h.

{ return m_xformSupportC; }

Member Data Documentation

m_DDmgr

InDetDD::PixelDetectorManager* GeoVPixelFactory::m_DDmgr

protectedinherited

Definition at line 45 of file GeoVPixelFactory.h.

m_epsilon

const double GeoVPixelFactory::m_epsilon

protectedinherited

Definition at line 49 of file GeoVPixelFactory.h.

m_gmt_mgr

PixelGeometryManager* GeoVPixelFactory::m_gmt_mgr

protectedinherited

Definition at line 43 of file GeoVPixelFactory.h.

m_mapAX

std::shared_ptr<std::map<std::string, GeoAlignableTransform*> > GeoVPixelFactory::m_mapAX

protectedinherited

Definition at line 48 of file GeoVPixelFactory.h.

m_mapFPV

std::shared_ptr<std::map<std::string, GeoFullPhysVol*> > GeoVPixelFactory::m_mapFPV

protectedinherited

Definition at line 47 of file GeoVPixelFactory.h.

m_mat_mgr

InDetMaterialManager* GeoVPixelFactory::m_mat_mgr

protectedinherited

Definition at line 44 of file GeoVPixelFactory.h.

m_sqliteReader

GeoModelIO::ReadGeoModel* GeoVPixelFactory::m_sqliteReader

protectedinherited

Definition at line 46 of file GeoVPixelFactory.h.

m_supportPhysA

GeoPhysVol* GeoPixelIBLFwdSvcCADModel::m_supportPhysA

private

Definition at line 38 of file GeoPixelIBLFwdSvcCADModel.h.

m_supportPhysC

GeoPhysVol* GeoPixelIBLFwdSvcCADModel::m_supportPhysC

private

Definition at line 39 of file GeoPixelIBLFwdSvcCADModel.h.

m_xformSupportA

GeoTransform* GeoPixelIBLFwdSvcCADModel::m_xformSupportA

private

Definition at line 40 of file GeoPixelIBLFwdSvcCADModel.h.

m_xformSupportC

GeoTransform* GeoPixelIBLFwdSvcCADModel::m_xformSupportC

private

Definition at line 41 of file GeoPixelIBLFwdSvcCADModel.h.

---

The documentation for this class was generated from the following files:

• GeoPixelIBLFwdSvcCADModel.h
• GeoPixelIBLFwdSvcCADModel.cxx