#include <H6CryostatConstruction.h>

Collaboration diagram for LArGeo::H6CryostatConstruction:

Public Member Functions
	H6CryostatConstruction ()

virtual	~H6CryostatConstruction ()

virtual GeoIntrusivePtr< GeoVFullPhysVol >	GetEnvelope ()

GeoIntrusivePtr< GeoPhysVol >	GetLArPhysical ()

Private Attributes
GeoFullPhysVol *	m_cryoMotherPhysical

GeoPhysVol *	m_cryoLArPhys

Detailed Description

Definition at line 22 of file H6CryostatConstruction.h.

Constructor & Destructor Documentation

◆ H6CryostatConstruction()

LArGeo::H6CryostatConstruction::H6CryostatConstruction ( )

Definition at line 46 of file H6CryostatConstruction.cxx.

                                                     :
   m_cryoMotherPhysical(nullptr),
   m_cryoLArPhys(nullptr)
 {}

◆ ~H6CryostatConstruction()

LArGeo::H6CryostatConstruction::~H6CryostatConstruction ( )

virtualdefault

Member Function Documentation

◆ GetEnvelope()

GeoIntrusivePtr< GeoVFullPhysVol > LArGeo::H6CryostatConstruction::GetEnvelope ( )

virtual

Definition at line 53 of file H6CryostatConstruction.cxx.

 {
  
   if (m_cryoMotherPhysical) return m_cryoMotherPhysical;
  
   // Detector Store
   ISvcLocator *svcLocator = Gaudi::svcLocator();
   StoreGateSvc *detStore;
   if (svcLocator->service("DetectorStore", detStore, false )==StatusCode::FAILURE) {
     throw std::runtime_error("Error in H6CryostatConstruction, cannot access DetectorStore");
   }
  
   // Material Manager
   StoredMaterialManager* materialManager = nullptr;
   if (StatusCode::SUCCESS != detStore->retrieve(materialManager, std::string("MATERIALS"))) return nullptr;
  
   const GeoMaterial *Air  = materialManager->getMaterial("std::Air");
   if (!Air) {
     throw std::runtime_error("Error in H6CryostatConstruction, std::Air is not found.");
   }
   const GeoMaterial *Iron  = materialManager->getMaterial("std::Iron");
   if (!Iron) {
     throw std::runtime_error("Error in H6CryostatConstruction, std::Iron is not found.");
   }
   const GeoMaterial *LAr  = materialManager->getMaterial("std::LiquidArgon");
   if (!LAr) {
     throw std::runtime_error("Error in H6CryostatConstruction, std::LiquidArgon is not found.");
   }
  
   // GeoModelSvc
   ServiceHandle<IGeoModelSvc> geoModelSvc ("GeoModelSvc", "H6CryostatConstruction");
   if (geoModelSvc.retrieve().isFailure()) {
     throw std::runtime_error ("Cannot locate GeoModelSvc!!");
   }
   
   std::string AtlasVersion = geoModelSvc->atlasVersion();
   std::string LArVersion = geoModelSvc->LAr_VersionOverride();
  
   std::string detectorKey  = LArVersion.empty() ? AtlasVersion : LArVersion;
   std::string detectorNode = LArVersion.empty() ? "ATLAS" : "LAr";
  
   //--------
   // Cryostat 
  
   // Here we should construct some sort of a cryostat:
   // First attempt at creating the H6 cryostat: 
   // (values taken from  HECCommonDetectorParamDef)
   //
   // A cylinder of half-height:   zcryo = 2000.0 *Gaudi::Units::mm
   // outer radius of warm-wall:   rwarm = 1295.5 *Gaudi::Units::mm
   // outer radius of vacuum:      rvac  = 1293.0 *Gaudi::Units::mm
   // outer radius of cold wall:   rcold = 1258.0 *Gaudi::Units::mm
   // outer radius of LAr:         rlar  = 1255.0 *Gaudi::Units::mm
  
   // needs to go into database: ---  4databa
   double  zcryo = 2000.0 *Gaudi::Units::mm;
   double  rwarm = 1295.5 *Gaudi::Units::mm;
   double  rvac  = 1293.0 *Gaudi::Units::mm;
   double  rcold = 1258.0 *Gaudi::Units::mm;
   double  rlar  = 1255.0 *Gaudi::Units::mm;
  
   std::string cryoMotherName = "LAr::H6::Cryostat::MotherVolume";
  
   // mother volume; cylinder of radius rwarm = outside of the cryostat warm wall 
   GeoTube* cryoMotherShape = new GeoTube(0.0 , rwarm+10.0, zcryo+10.0*Gaudi::Units::mm);  //  mother is a little bigger than warm wall   
  
   const GeoLogVol* cryoMotherLogical = new GeoLogVol(cryoMotherName, cryoMotherShape, Air);
   m_cryoMotherPhysical = new GeoFullPhysVol(cryoMotherLogical);
  
   // Now we have a physical air-filled cylindrical mother, into which we can place all we want and need.
  
   // First insert the Cryostat walls and the LAr - all simple cylinders
   //
   // 4databa ???
   std::string baseName          = "LArTB::H6::Cryostat::" ;
   std::string cryoWallName      = "LArTB::H6::Cryostat" ;
   std::string cryoWarmWallName  = baseName + "WarmWall" ;
   std::string cryoVacuumGapName = baseName + "Vacuum" ;
   std::string cryoColdWallName  = baseName + "ColdWall" ;
   std::string cryoLArName       = baseName + "LAr" ;     // <--- used to be called "Inside"
  
  
   // Warm Wall 
   GeoTube* cryoWarmWallShape = new GeoTube(0. , rwarm, zcryo);   
   const GeoLogVol* cryoWarmWallLog = new GeoLogVol(cryoWarmWallName, cryoWarmWallShape, Iron);
   //m_cryoMotherPhysical->add(new GeoNameTag(std::string("Cryostat")));    
   GeoIntrusivePtr<GeoPhysVol>  cryoWarmWallPhys = new GeoPhysVol(cryoWarmWallLog);
   m_cryoMotherPhysical->add(cryoWarmWallPhys);
  
   // "Vacuum" gap (filled with air...)
   GeoTube* cryoVacuumGapShape = new GeoTube(0. , rvac, zcryo-2.0*Gaudi::Units::mm);   // an arbitrary 2mm shorter to avoid confilct  
   const GeoLogVol* cryoVacuumGapLog = new GeoLogVol(cryoVacuumGapName, cryoVacuumGapShape, Air);
   GeoIntrusivePtr<GeoPhysVol>  cryoVacuumGapPhys = new GeoPhysVol(cryoVacuumGapLog);
   cryoWarmWallPhys->add(cryoVacuumGapPhys);
  
   // Cold Wall
   GeoTube* cryoColdWallShape = new GeoTube(0. , rcold, zcryo-4.0*Gaudi::Units::mm);  // an arbitrary 4mm shorter to avoid confilct   
   const GeoLogVol* cryoColdWallLog = new GeoLogVol(cryoColdWallName, cryoColdWallShape, Iron);
   GeoIntrusivePtr<GeoPhysVol>  cryoColdWallPhys = new GeoPhysVol(cryoColdWallLog);
   cryoVacuumGapPhys->add(cryoColdWallPhys);
  
  
  
   // At the moment the CryostatMother is LAr, so the FCal should be in there,
   // what would be better would be a mother of Air and LAr separately inside 
   // And the FCal the embedded in the LAr instead of the cryoMother!
  
   // Liquid Argon  
   GeoTube* cryoLArShape = new GeoTube(0. , rlar, zcryo-6.0*Gaudi::Units::mm);  // an arbitrary 2mm shorter to avoid confilct   
   const GeoLogVol* cryoLArLog = new GeoLogVol(cryoLArName, cryoLArShape, LAr);
   m_cryoMotherPhysical->add(new GeoNameTag(std::string("Cryostat LAr Physical")));    
   // m_cryoLArPhys is a class member so that we can place Detectors inside. 
   m_cryoLArPhys = new GeoPhysVol(cryoLArLog);
   cryoColdWallPhys->add(m_cryoLArPhys);
  
  
   // What remains needs to be inserted:
   // - Rohacell  (Needs an entry in LArCalorimeter/LArGeoModel/LArGeoAlgs/src/LArMaterialManager.cxx)
  
   return m_cryoMotherPhysical;
  
 }

◆ GetLArPhysical()

GeoIntrusivePtr< GeoPhysVol > LArGeo::H6CryostatConstruction::GetLArPhysical ( )

Definition at line 177 of file H6CryostatConstruction.cxx.

 {
   return m_cryoLArPhys;
 }

Member Data Documentation

◆ m_cryoLArPhys

GeoPhysVol* LArGeo::H6CryostatConstruction::m_cryoLArPhys

private

Definition at line 39 of file H6CryostatConstruction.h.

◆ m_cryoMotherPhysical

GeoFullPhysVol* LArGeo::H6CryostatConstruction::m_cryoMotherPhysical