LArGeo::HEC2WheelConstruction Class Reference

GeoModel description of LAr HEC. More...

#include <HEC2WheelConstruction.h>

Collaboration diagram for LArGeo::HEC2WheelConstruction:

Public Member Functions
	HEC2WheelConstruction ()=default
	~HEC2WheelConstruction ()=default
GeoIntrusivePtr< GeoFullPhysVol >	GetEnvelope (bool fullGeo, bool posZSide=true)

Private Attributes
GeoIntrusivePtr< GeoFullPhysVol >	m_physiHEC {}
bool	m_posZSide {false}

Detailed Description

GeoModel description of LAr HEC.

The geometry is built and placed within HEC envelope which is implemented by LArGeoEndcap.

Definition at line 25 of file HEC2WheelConstruction.h.

Constructor & Destructor Documentation

HEC2WheelConstruction()

LArGeo::HEC2WheelConstruction::HEC2WheelConstruction ( )

default

~HEC2WheelConstruction()

LArGeo::HEC2WheelConstruction::~HEC2WheelConstruction ( )

default

Member Function Documentation

GetEnvelope()

GeoIntrusivePtr< GeoFullPhysVol > LArGeo::HEC2WheelConstruction::GetEnvelope	(	bool	fullGeo,
		bool	posZSide = true )

Definition at line 74 of file HEC2WheelConstruction.cxx.

{
  m_posZSide = posZSide;
 
  if(m_physiHEC) return m_physiHEC->clone();
 
  // Get access to the material manager:
  ISvcLocator *svcLocator = Gaudi::svcLocator();
 
  MsgStream log(Athena::getMessageSvc(),"HEC2WheelConstruction " );
  
 
  log << MSG::DEBUG << "+++++++++++++++++++++++++++++++++++++++++++++++++++++" << endmsg;
  log << MSG::DEBUG << "+                                                   +" << endmsg;
  log << MSG::DEBUG << "+         Start of HEC2Wheel GeoModel definition    +" << endmsg;
  log << MSG::DEBUG << "+                                                   +" << endmsg;
  log << MSG::DEBUG << "+++++++++++++++++++++++++++++++++++++++++++++++++++++" << endmsg;
 
 
  SmartIF<StoreGateSvc> detStore{svcLocator->service("DetectorStore")};
  if(!detStore.isValid()) {
    throw std::runtime_error("Error in HEC2WheelConstruction, cannot access DetectorStore");
  }
 
  StoredMaterialManager* materialManager = nullptr;
  if (StatusCode::SUCCESS != detStore->retrieve(materialManager, std::string("MATERIALS"))) {
    throw std::runtime_error("Error in HEC2WheelConstruction, cannot access Material Manager");
  }
 
  SmartIF<IRDBAccessSvc> pAccessSvc{svcLocator->service("RDBAccessSvc")};
  if(!pAccessSvc.isValid()) {
    throw std::runtime_error ("Cannot locate RDBAccessSvc!!");
  }
 
  SmartIF<IGeoModelSvc> geoModel{svcLocator->service("GeoModelSvc")};
  if(!geoModel.isValid()) {
    throw std::runtime_error ("Cannot locate GeoModelSvc!!");
  }
 
  std::string AtlasVersion = geoModel->atlasVersion();
  std::string LArVersion = geoModel->LAr_VersionOverride();
  std::string detectorKey  = LArVersion.empty() ? AtlasVersion : LArVersion;
  std::string detectorNode = LArVersion.empty() ? "ATLAS" : "LAr";
 
  IRDBRecordset_ptr hadronicEndcap = pAccessSvc->getRecordsetPtr("HadronicEndcap",detectorKey, detectorNode); 
  if(hadronicEndcap->size()>0)
    log << MSG::ALWAYS << "Using numbers from HadronicEndcap tag:  " << hadronicEndcap->tagName() << endmsg;
  else
    throw std::runtime_error("Error in HECConstruction:  hadronicEendcap not found");
 
  IRDBRecordset_ptr hecLongitudinalBlock = pAccessSvc->getRecordsetPtr("HecLongitudinalBlock",detectorKey, detectorNode); 
  if(hecLongitudinalBlock->size()>0) 
    log << MSG::ALWAYS << "Using numbers from HecLongitudinalBlock tag:  " << hecLongitudinalBlock->tagName() << endmsg;
  else 
    throw std::runtime_error("Error in HECConstruction:  hecLongitudinalBlock not found");
 
  IRDBRecordset_ptr larPosition  =  pAccessSvc->getRecordsetPtr("LArPosition",detectorKey, detectorNode);
  if (larPosition->size()==0 ) {
    larPosition = pAccessSvc->getRecordsetPtr("LArPosition", "LArPosition-00");
    if (larPosition->size()==0 ) {
      throw std::runtime_error("Error, no lar position table in database!");
    }
  }
 
  // This is a start for letting everyhting go cold (shrinkCold = 1.0 means warm; = 0.997 means cold)
  double shrinkCold            = 1.0 ; 
 
  double rInner1               = shrinkCold * (*hecLongitudinalBlock)[0]->getDouble("BLRMN")*cm;
  double rInner2               = shrinkCold * (*hecLongitudinalBlock)[1]->getDouble("BLRMN")*cm;
  double rOuter                = shrinkCold * (*hecLongitudinalBlock)[0]->getDouble("BLRMX")*cm;
  double wheelGap              = shrinkCold * (*hadronicEndcap)[0]->getDouble("GAPWHL")*cm;
 
  
  const GeoMaterial *LAr  = materialManager->getMaterial("std::LiquidArgon");
  if (!LAr) throw std::runtime_error("Error in HEC2WheelConstruction, std::LiquidArgon is not found.");
 
  GeoPcon*           solidHEC;     //pointer to the solid HEC
  const GeoLogVol*   logicHEC;     //pointer to the logical HEC
  GeoIntrusivePtr<GeoFullPhysVol>    physiHEC;     //pointer to the physical HEC
 
  
 
  // YIKES! Hard-coded numbers!!!
 
  if (rOuter<2100.*mm) rOuter  = shrinkCold *2130*mm;  // Needs fixing in database to make room for HEC Clamping bars! 
                                           // Caution: We'll need 2 Routers - One for Wheel, one for Module!
 
  double     HEC1length = shrinkCold * 816.5*mm;
  double     HEC2length = shrinkCold * 961.0*mm;
  double     TotalWarmHEC = HEC1length + wheelGap + HEC2length ;  // For the moment, this defines the HEC rear face  
  double     TotalHECLength = TotalWarmHEC;                      // (where is that -cold- in the cryostat?
 
  double     zCoordinate[4];
  zCoordinate[0]=shrinkCold * 0.0*cm; 
  zCoordinate[1]=shrinkCold * 28.05*cm; 
  zCoordinate[2]=shrinkCold * 28.051*cm; 
  zCoordinate[3]=shrinkCold * (181.9*cm+4.0*mm); // JFB, MGF:  more room in the mother fiducal volume for misalignment.
 // 37.2, 47.5, 203 cm for R
 
  solidHEC = new GeoPcon(0., 2.*M_PI);
  for (int i=0; i< 4; i++) 
    {
      if (i<2) solidHEC->addPlane(zCoordinate[i],rInner1 - shrinkCold*1.*mm,rOuter);
      else     solidHEC->addPlane(zCoordinate[i],rInner2 - shrinkCold*1.*mm,rOuter);
   }
  
  std::string hecName = "LAr::HEC::Mother";
  logicHEC = new GeoLogVol(hecName, solidHEC , LAr);
  physiHEC = new GeoFullPhysVol(logicHEC);
  
  //--- Make the Front Wheel alignable:
 
  const IRDBRecord *posHec1 = GeoDBUtils::getTransformRecord(larPosition, m_posZSide ? "HEC1_POS":"HEC1_NEG");
  Transform3D xfPosHec1 = posHec1 ? GeoDBUtils::getTransform(posHec1) : Translate3D(0.,0.,0.);
  GeoAlignableTransform *xfHec1 = new GeoAlignableTransform(xfPosHec1);
 
  std::string tag1 = m_posZSide? std::string("HEC1_POS") : std::string("HEC1_NEG");
 
  HECWheelConstruction theFrontHEC(fullGeo,"front",false,m_posZSide) ;
  GeoIntrusivePtr<GeoFullPhysVol> EnvelopeF = theFrontHEC.GetEnvelope();
 
  StoredPhysVol *sPhysVolHec1 = new StoredPhysVol(EnvelopeF);
  StatusCode sc=detStore->record(sPhysVolHec1,tag1);
  if(!sc.isSuccess()) throw std::runtime_error ((std::string("Cannot store")+tag1).c_str());
 
  StoredAlignX *sAlignX1 = new StoredAlignX(xfHec1);
  sc=detStore->record(sAlignX1,tag1);
  if(!sc.isSuccess()) throw std::runtime_error ((std::string("Cannot store")+tag1).c_str());
 
  physiHEC->add(xfHec1);
  physiHEC->add(new GeoIdentifierTag(0));
  physiHEC->add(EnvelopeF);
 
  //--- Now the Rear Wheel:
 
  double Zpos = TotalHECLength - HEC2length  ;  // It is fixed at its end; position it from there.
  const IRDBRecord *posHec2 = GeoDBUtils::getTransformRecord(std::move(larPosition), m_posZSide ? "HEC2_POS":"HEC2_NEG");
  Transform3D xfPosHec2 = posHec2 ? GeoDBUtils::getTransform(posHec2):Translate3D(0.,0.,Zpos);
 
 
  GeoAlignableTransform *xfHec2 = new GeoAlignableTransform(xfPosHec2);
 
  std::string tag2 = m_posZSide? std::string("HEC2_POS") : std::string("HEC2_NEG");
 
  HECWheelConstruction theRearHEC(fullGeo,"rear",false,m_posZSide);
  GeoIntrusivePtr<GeoFullPhysVol> EnvelopeR = theRearHEC.GetEnvelope();
 
  StoredPhysVol *sPhysVolHec2 = new StoredPhysVol(EnvelopeR);
  sc=detStore->record(sPhysVolHec2,tag2);
  if(!sc.isSuccess()) throw std::runtime_error ((std::string("Cannot store")+tag2).c_str());
 
  StoredAlignX *sAlignX2 = new StoredAlignX(xfHec2);
  sc=detStore->record(sAlignX2,tag2);
  if(!sc.isSuccess()) throw std::runtime_error ((std::string("Cannot store")+tag2).c_str());
 
  physiHEC->add(xfHec2);
  physiHEC->add(new GeoIdentifierTag(1));
  physiHEC->add(EnvelopeR);
 
  return physiHEC;
}

Member Data Documentation

m_physiHEC

GeoIntrusivePtr<GeoFullPhysVol> LArGeo::HEC2WheelConstruction::m_physiHEC {}

private

Definition at line 35 of file HEC2WheelConstruction.h.

{};

m_posZSide

bool LArGeo::HEC2WheelConstruction::m_posZSide {false}

private

Definition at line 36 of file HEC2WheelConstruction.h.

{false};

---

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

• HEC2WheelConstruction.h
• HEC2WheelConstruction.cxx