BDSTunnelFactorySquare.cc

/* 
Beam Delivery Simulation (BDSIM) Copyright (C) Royal Holloway, 
University of London 2001 - 2023.
 
This file is part of BDSIM.
 
BDSIM is free software: you can redistribute it and/or modify 
it under the terms of the GNU General Public License as published 
by the Free Software Foundation version 3 of the License.
 
BDSIM is distributed in the hope that it will be useful, but 
WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU General Public License for more details.
 
You should have received a copy of the GNU General Public License
along with BDSIM.  If not, see <http://www.gnu.org/licenses/>.
*/
#include "BDSTunnelFactoryBase.hh"
#include "BDSTunnelFactorySquare.hh"
#include "BDSTunnelFactoryRectangular.hh"
#include "BDSTunnelInfo.hh"
 
#include "globals.hh"
 
BDSTunnelFactorySquare::BDSTunnelFactorySquare()
{
  rectangular = new BDSTunnelFactoryRectangular();
}
 
BDSTunnelFactorySquare::~BDSTunnelFactorySquare()
{
  delete rectangular;
}
 
BDSTunnelSection* BDSTunnelFactorySquare::CreateTunnelSection(G4String    name,
                                                              G4double    length,
                                                              G4double    tunnelThickness,
                                                              G4double    tunnelSoilThickness,
                                                              G4Material* tunnelMaterial,
                                                              G4Material* tunnelSoilMaterial,
                                                              G4bool      tunnelFloor,
                                                              G4double    tunnelFloorOffset,
                                                              G4double    tunnel1,
                                                              G4double    /*tunnel2*/,
                                                              G4bool      visible)
{
  // test input parameters - set global options as default if not specified
  TestInputParameters(length, tunnelThickness, tunnelSoilThickness, tunnelMaterial,
                      tunnelSoilMaterial, tunnelFloorOffset, tunnel1);
  return rectangular->CreateTunnelSection(name, length, tunnelThickness,
                                          tunnelSoilThickness, tunnelMaterial,
                                          tunnelSoilMaterial, tunnelFloor,
                                          tunnelFloorOffset, tunnel1, tunnel1, visible);
}
 
BDSTunnelSection* BDSTunnelFactorySquare::CreateTunnelSectionAngled(G4String      name,
                                                                    G4double      length,
                                                                    G4ThreeVector inputFace,
                                                                    G4ThreeVector outputFace,
                                                                    G4double      tunnelThickness,
                                                                    G4double      tunnelSoilThickness,
                                                                    G4Material*   tunnelMaterial,
                                                                    G4Material*   tunnelSoilMaterial,
                                                                    G4bool        tunnelFloor,
                                                                    G4double      tunnelFloorOffset,
                                                                    G4double      tunnel1,
                                                                    G4double      /*tunnel2*/,
                                                                    G4bool        visible)
{
  // test input parameters - set global options as default if not specified
  TestInputParameters(length, tunnelThickness, tunnelSoilThickness, tunnelMaterial,
                      tunnelSoilMaterial, tunnelFloorOffset, tunnel1);
  return rectangular->CreateTunnelSectionAngled(name, length, inputFace, outputFace,
                                                tunnelThickness,
                                                tunnelSoilThickness, tunnelMaterial,
                                                tunnelSoilMaterial, tunnelFloor,
                                                tunnelFloorOffset, tunnel1, tunnel1, visible);
}
 
void BDSTunnelFactorySquare::TestInputParameters(G4double&    length,
                                                 G4double&    tunnelThickness,
                                                 G4double&    tunnelSoilThickness,
                                                 G4Material*& tunnelMaterial,
                                                 G4Material*& tunnelSoilMaterial,
                                                 G4double&    tunnelFloorOffset,
                                                 G4double&    tunnel1)
{
  CommontTestInputParameters(length, tunnelThickness, tunnelSoilThickness, tunnelMaterial, tunnelSoilMaterial);
  
  if (tunnelFloorOffset < 1e-10)
    {tunnelFloorOffset = defaultModel->floorOffset;}
  
  if (tunnel1 < 1e-10)
    {tunnel1 = defaultModel->aper1;}
}