#include <MiddleBeamConstructionH62004.h>

Collaboration diagram for LArGeo::MiddleBeamConstructionH62004:

Public Member Functions
	MiddleBeamConstructionH62004 ()

virtual	~MiddleBeamConstructionH62004 ()

virtual PVLink	GetEnvelope ()

void	SetManager (LArDetDescrManager *mgr)

Private Attributes
GeoPhysVol *	m_H62004MiddleBeamPhysical

LArDetDescrManager *	m_detectorManager

Detailed Description

Definition at line 18 of file MiddleBeamConstructionH62004.h.

Constructor & Destructor Documentation

◆ MiddleBeamConstructionH62004()

LArGeo::MiddleBeamConstructionH62004::MiddleBeamConstructionH62004 ( )

Definition at line 47 of file MiddleBeamConstructionH62004.cxx.

   :m_H62004MiddleBeamPhysical(nullptr),
    m_detectorManager(nullptr)
 {
 }

◆ ~MiddleBeamConstructionH62004()

LArGeo::MiddleBeamConstructionH62004::~MiddleBeamConstructionH62004 ( )

virtualdefault

Member Function Documentation

◆ GetEnvelope()

PVLink LArGeo::MiddleBeamConstructionH62004::GetEnvelope ( )

virtual

Definition at line 59 of file MiddleBeamConstructionH62004.cxx.

 {
  
   if (m_H62004MiddleBeamPhysical) return m_H62004MiddleBeamPhysical;
  
   // Get access to the material manager:
   
   ISvcLocator *svcLocator = Gaudi::svcLocator();
   IMessageSvc * msgSvc;
   if (svcLocator->service("MessageSvc", msgSvc, true )==StatusCode::FAILURE) {
     throw std::runtime_error("Error in MiddleBeamConstructionH62004, cannot access MessageSvc");
   }
  
   MsgStream log(msgSvc, "LArGeo::MiddleBeamConstructionH62004"); 
   log << MSG::INFO << "+    HELLO from LArGeo::MiddleBeamConstructionH62004       +" << endmsg;
  
  
   StoreGateSvc *detStore;
   if (svcLocator->service("DetectorStore", detStore, false )==StatusCode::FAILURE) {
     throw std::runtime_error("Error in MiddleBeamConstructionH62004, cannot access DetectorStore");
   }
  
   // Get the materials from the 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 MiddleBeamConstructionH62004, std::Air is not found.");
    
    //-------------------------------------------------------------------------------------------------//
  
  
   DecodeVersionKey larVersion("LAr");
    
  
   // Define geometry
  
  
   //  Here we creat the envelope for the Moveable FrontBeam Instrumentation.  This code is repeated
   //  createEnvelope() method below.  There should be a way to avoid this repitition.
  
   //------ The FrontBeam 
  
   std::string baseName = "LAr::TBH62004";
   std::string H62004MiddleBeamName = baseName + "::MiddleBeam";
  
   //
   // Define dimension of Middle part & position of Front part
   //
   double bmtb_x = 12.0*Gaudi::Units::cm;
   double bmtb_y = 12.0*Gaudi::Units::cm;
   double bmtb_z = 25.0*Gaudi::Units::cm;
   //double bmtb_pos = 10.0*Gaudi::Units::cm;
   double bpco_pos[4] =  {1.*Gaudi::Units::cm, 1.*Gaudi::Units::cm, 15.3*Gaudi::Units::cm, 15.3*Gaudi::Units::cm};
   double bpco_shift[4] =  {0.*Gaudi::Units::cm, 7.8*Gaudi::Units::cm, 0.*Gaudi::Units::cm, 7.5*Gaudi::Units::cm};
   double bpc_old_z = (5.100/2)*Gaudi::Units::cm;
  
   GeoBox* H62004MiddleBeamShape = new GeoBox( bmtb_x, bmtb_y, bmtb_z );   
   const GeoLogVol* H62004MiddleBeamLogical = new GeoLogVol( H62004MiddleBeamName, H62004MiddleBeamShape, Air );
  
   m_H62004MiddleBeamPhysical = new GeoPhysVol(H62004MiddleBeamLogical);
  
  
   //----- Now create BPC
   log << MSG::INFO << " Create BPC 1&2 " << endmsg;
  
   BPCConstruction *BPC = new BPCConstruction(true);
   PVLink BPCPhysical = BPC->GetEnvelope();
   for(int i=0; i<4; ++i) {
      m_H62004MiddleBeamPhysical->add( new GeoIdentifierTag((3+i/2)*10+i) );
      switch(i) {
       case 0: case 2: {
                  m_H62004MiddleBeamPhysical->add( new GeoTransform(GeoTrf::Translate3D( 0.*Gaudi::Units::cm, 0.*Gaudi::Units::cm, bpco_pos[i]+bpco_shift[i]+bpc_old_z-bmtb_z) ) );
                             m_H62004MiddleBeamPhysical->add(BPCPhysical);
                 break;}
       case 1: case 3: { 
                             m_H62004MiddleBeamPhysical->add( new GeoTransform(GeoTrf::RotateZ3D(90.*Gaudi::Units::deg) *  GeoTrf::Translate3D( 0.*Gaudi::Units::cm, 0.*Gaudi::Units::cm, bpco_pos[i]+bpco_shift[i]+bpc_old_z-bmtb_z) ) );
                             m_H62004MiddleBeamPhysical->add(BPCPhysical);
                 break;}
     }
   }
  
   // End Middle FrontBeam detectors
  
  
   return m_H62004MiddleBeamPhysical;
 }