BDSBeamPipeFactoryLHCDetailed.hh

/* 
Beam Delivery Simulation (BDSIM) Copyright (C) Royal Holloway, 
University of London 2001 - 2022.
 
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/>.
*/
#ifndef BDSBEAMPIPEFACTORYLHCDETAILED_H
#define BDSBEAMPIPEFACTORYLHCDETAILED_H
 
#include "BDSBeamPipeFactoryBase.hh"
#include "BDSBeamPipe.hh"
 
class G4Material;
 
class BDSBeamPipeFactoryLHCDetailed: public BDSBeamPipeFactoryBase
{
public:
  BDSBeamPipeFactoryLHCDetailed();
  virtual ~BDSBeamPipeFactoryLHCDetailed(){;}
 
  G4double GetFullWidthOfCoolingPipe();
 
  virtual void CleanUp();
 
  virtual BDSBeamPipe* CreateBeamPipe(const G4String& nameIn,
                                      G4double    lengthIn,
                                      G4double    aper1               = 0,
                                      G4double    aper2               = 0,
                                      G4double    aper3               = 0,
                                      G4double    aper4               = 0,
                                      G4Material* vacuumMaterialIn    = nullptr,
                                      G4double    beamPipeThicknessIn = 0,
                                      G4Material* beamPipeMaterialIn  = nullptr,
                                      const G4String& pointsFileIn    = "",
                                      const G4String& pointsUnitIn    = "");
  
  virtual BDSBeamPipe* CreateBeamPipe(const G4String&      nameIn,
                                      G4double             lengthIn,
                                      const G4ThreeVector& inputFaceNormalIn,
                                      const G4ThreeVector& outputFaceNormalIn,
                                      G4double      aper1               = 0,
                                      G4double      aper2               = 0,
                                      G4double      aper3               = 0,
                                      G4double      aper4               = 0,
                                      G4Material*   vacuumMaterialIn    = nullptr,
                                      G4double      beamPipeThicknessIn = 0,
                                      G4Material*   beamPipeMaterialIn  = nullptr,
                                      const G4String& pointsFileIn      = "",
                                      const G4String& pointsUnitIn      = "");
 
private:
  void CleanUpLHCDetailed();
  
  BDSBeamPipe* CommonFinalConstruction(const G4String& name,
                                       G4Material* vacuumMaterial,
                                       G4Material* beamPipeMaterial,
                                       G4double    length,
                                       G4double    containerRadius);
  G4double CreateGeneralAngledSolids(const G4String&      name,
                                     G4double             length,
                                     const G4ThreeVector& inputface,
                                     const G4ThreeVector& outputface);
 
  virtual void BuildLogicalVolumes(const G4String& nameIn,
                                   G4Material* vacuumMaterialIn,
                                   G4Material* beamPipeMaterialIn);
  virtual void  SetVisAttributes(G4Material* beamPipeMaterialIn);
  virtual void  SetUserLimits(G4double length);
  virtual void  PlaceComponents(const G4String& nameIn);
 
  void InitialiseGeometricalParameters();
 
  void CalculateGeometricalParameters(G4double aper1,
                                      G4double aper2,
                                      G4double aper3,
                                      G4double beamPipeThickness,
                                      G4double length);
 
  G4VSolid*        copperSkinSolid;
  G4VSolid*        screenSolid;
  G4VSolid*        coolingPipeSolid;
  G4LogicalVolume* copperSkinLV;
  G4LogicalVolume* screenLV;
  G4LogicalVolume* coolingPipeLV;
  G4PVPlacement*   copperSkinPV;
  G4PVPlacement*   screenPV;
  G4PVPlacement*   coolingPipeTopPV;
  G4PVPlacement*   coolingPipeBottomPV;
 
  G4double coldBoreThickness;
  G4double coolingPipeThickness;
  G4double coolingPipeRadius;
  G4double copperSkinThickness;
 
  G4double vacRadius, vacBoxX, vacBoxY;
  
  G4double cuInnerRadius, cuInnerBoxX, cuInnerBoxY;
  G4double cuOuterRadius, cuOuterBoxX, cuOuterBoxY;
 
  G4double bsInnerRadius, bsInnerBoxX, bsInnerBoxY;
  G4double bsOuterRadius, bsOuterBoxX, bsOuterBoxY;
  
  G4double cbInnerRadius, cbOuterRadius;
 
  G4double containerRadius;
  
  G4double vacHalfLength;
  G4double halfLength;
 
  G4double coolingPipeOffset;
 
  G4bool verticalOrientation;
 
  G4bool buildCoolingPipe; 
};
  
#endif