LArGeo::ExcluderConstruction Class Reference

#include <ExcluderConstruction.h>

Collaboration diagram for LArGeo::ExcluderConstruction:

Public Member Functions
	ExcluderConstruction ()
virtual	~ExcluderConstruction ()
virtual GeoIntrusivePtr< GeoPhysVol >	GetEnvelope ()

Detailed Description

Definition at line 18 of file ExcluderConstruction.h.

Constructor & Destructor Documentation

ExcluderConstruction()

LArGeo::ExcluderConstruction::ExcluderConstruction ( )

default

~ExcluderConstruction()

LArGeo::ExcluderConstruction::~ExcluderConstruction ( )

virtualdefault

Member Function Documentation

GetEnvelope()

GeoIntrusivePtr< GeoPhysVol > LArGeo::ExcluderConstruction::GetEnvelope ( )

virtual

Definition at line 54 of file ExcluderConstruction.cxx.

{
 
  // Need to do the equivalent for excluder here:                   <<<================ 
  //  if (cryoMotherPhysical) return cryoMotherPhysical;
 
  // Detector Store
  ISvcLocator *svcLocator = Gaudi::svcLocator();
  StoreGateSvc *detStore;
  if (svcLocator->service("DetectorStore", detStore, false )==StatusCode::FAILURE) {
    throw std::runtime_error("Error in ExcluderConstruction, cannot access DetectorStore");
  }
 
  // Material Manager
 
  // Need to add Rohacell here!                                     <<<=============== 
 
  StoredMaterialManager* materialManager = nullptr;
  if (StatusCode::SUCCESS != detStore->retrieve(materialManager, std::string("MATERIALS"))) return nullptr;
 
 
  // (use Air for the moment....)                                  <<<================ 
  const GeoMaterial *Air  = materialManager->getMaterial("std::Air");
  if (!Air) {
    throw std::runtime_error("Error in ExcluderConstruction, std::Air is not found.");
  }
 
  // Define Rohacell Foam. 
  // Rohacell foam has density: 0.11g/cm3
  std::string name;
  double density;
  const GeoElement*  C=materialManager->getElement("Carbon");
  const GeoElement*  H=materialManager->getElement("Hydrogen");
  const GeoElement*  O=materialManager->getElement("Oxygen");
  const GeoElement*  N=materialManager->getElement("Nitrogen");
  GeoMaterial* Rohacell = new GeoMaterial(name="Rohacell", density=0.11*GeoModelKernelUnits::g/Gaudi::Units::cm3);
  Rohacell->add(C,0.6465);
  Rohacell->add(H,0.07836);
  Rohacell->add(O,0.19137);
  Rohacell->add(N,0.08377);
  Rohacell->lock();
 
  ServiceHandle<IGeoDbTagSvc> geoDbTagSvc ("GeoDbTagSvc", "WallsConstruction");
  if (geoDbTagSvc.retrieve().isFailure()) {
    throw std::runtime_error ("Cannot locate GeoDbTagSvc!!");
  }
  
  std::string AtlasVersion = geoDbTagSvc->atlasVersion();
  std::string LArVersion = geoDbTagSvc->LAr_VersionOverride();
  
  std::string detectorKey  = LArVersion.empty() ? AtlasVersion : LArVersion;
  std::string detectorNode = LArVersion.empty() ? "ATLAS" : "LAr";
 
  //--------
  // Now build the actual Excluder.
 
  // It is a Union out of a GeoBox and a GeoTubs.
  // Box Dimensions:
  double   xbox  =  300.0 *Gaudi::Units::mm;
  double   ybox  =  160.0 *Gaudi::Units::mm;
  double   zbox  =  300.7 *Gaudi::Units::mm;
  //
  // Tubs Dimensions:
  double   ztubs =  300.0 *Gaudi::Units::mm;
  double  phitubs=   76.2 *Gaudi::Units::deg;
  double  delphi =   27.6 *Gaudi::Units::deg;
  double   rcold = 1249.5 *Gaudi::Units::mm;
  double   rmin  = 1220.0 *Gaudi::Units::mm;
 
  // The radius of the cryostat cold wall is: 1250 Gaudi::Units::mm
  // Before we make the union, we have to shift the box in y (that actually along the beam axis)
  //        and there, positive y goes from the cryostat centre towards the beam window.
 
  std::string ExcluderName = "LAr::H6::Cryostat::Excluder";
 
  GeoBox* rohaBox   = new GeoBox(xbox, ybox, zbox);                      //  The rectangular part of the excluder
  const GeoShapeShift & rohaBoxShift = (*rohaBox << GeoTrf::TranslateY3D(1062.85*Gaudi::Units::mm) );
  GeoTubs* rohaTubs = new GeoTubs(rmin, rcold, ztubs, phitubs, delphi);  //  The round part of the excluder  
 
  // Combine the two parts to make one excluder of the correct shape:
  const GeoShapeUnion* Excluder = new GeoShapeUnion(&rohaBoxShift, rohaTubs);
 
  const GeoLogVol* LogExcluder = new GeoLogVol(ExcluderName, Excluder, Rohacell);  // <<<==== air should be rohacell !!!
  GeoIntrusivePtr<GeoPhysVol> PhysExcluder = new GeoPhysVol(LogExcluder);
  PhysExcluder->add( new GeoNameTag(ExcluderName) );
 
  return PhysExcluder;
 
}

---

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

• ExcluderConstruction.h
• ExcluderConstruction.cxx