GeoPixelFluid.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
//
// This class builds one trapezoid
//
//
#include "GeoPixelFluid.h"
#include "GeoModelKernel/GeoTrd.h"
#include "GeoModelKernel/GeoLogVol.h"
#include "GeoModelKernel/GeoPhysVol.h"
#include "GeoModelKernel/GeoMaterial.h"
 
#include<cmath>
#include <sstream>
 
GeoPixelFluid::GeoPixelFluid(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 type)
  : GeoVPixelFactory (ddmgr, mgr, sqliteReader, mapFPV, mapAX)
{
  m_index    = m_gmt_mgr->PixelFluidIndex(type);
  double z1  = m_gmt_mgr->PixelFluidZ1(m_index);
  double z2  = m_gmt_mgr->PixelFluidZ2(m_index);
  m_posX     = m_gmt_mgr->PixelFluidX(m_index);
  m_posY     = m_gmt_mgr->PixelFluidY(m_index);
  m_posZ     = 0.5*(z1+z2);
}
 
GeoVPhysVol* GeoPixelFluid::Build()
{
  //
  // Dimensions
  //
  double z1       = m_gmt_mgr->PixelFluidZ1(m_index);
  double z2       = m_gmt_mgr->PixelFluidZ2(m_index);
  double thick1   = m_gmt_mgr->PixelFluidThick1(m_index);
  double thick2   = m_gmt_mgr->PixelFluidThick2(m_index);
  double width    = m_gmt_mgr->PixelFluidWidth(m_index);
 
  double length = std::abs(z2 - z1);
  double thickAtZmin = thick1;
  //  double thickAtZmax = thick2;
  if (z1 > z2) {
    thickAtZmin = thick2;
    //    thickAtZmax = thick1;
  }    
  
  std::ostringstream o;
  o << "Fluid" << m_gmt_mgr->PixelFluidType(m_index);
  std::string logName = o.str();
 
  const GeoMaterial* fluidMat = m_mat_mgr->getMaterial(m_gmt_mgr->PixelFluidMat(m_index));
  const GeoTrd* fluidShape  = new GeoTrd(0.5*thickAtZmin, 0.5*thickAtZmin, 0.5*width, 0.5*width, 0.5*length);
  GeoLogVol * fluidLV    = new GeoLogVol(logName, fluidShape, fluidMat);
  GeoPhysVol * fluidPhys = new GeoPhysVol(fluidLV);
 
  return fluidPhys;
}