BDSCurvilinearBuilder.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/>.
*/
#include "BDSBeamlineElement.hh"
#include "BDSUtilities.hh"
 
#include "globals.hh"
 
class BDSBeamline;
class BDSCurvilinearFactory;
class BDSSimpleComponent;
 
class BDSCurvilinearBuilder
{
public:
  BDSCurvilinearBuilder();
  ~BDSCurvilinearBuilder();
 
  BDSBeamline* BuildCurvilinearBeamLine1To1(BDSBeamline const* const beamline,
                                            const G4bool circular);
 
  BDSBeamline* BuildCurvilinearBridgeBeamLine(BDSBeamline const* const beamline);
 
private:
  BDSCurvilinearBuilder(const BDSCurvilinearBuilder&) = delete;
  BDSCurvilinearBuilder& operator=(const BDSCurvilinearBuilder&) = delete;
 
  G4double CurvilinearRadius(const BDSBeamlineElement* el) const;
 
  BDSBeamlineElement* CreateCurvilinearElement(const G4String&             elementName,
                                               BDSBeamline::const_iterator startElement,
                                               BDSBeamline::const_iterator finishElement,
                                               G4int                       index,
                                               G4double                    crRadius);
 
  BDSBeamlineElement* CreateElementFromComponent(BDSSimpleComponent*         component,
                                                 BDSBeamline::const_iterator startElement,
                                                 BDSBeamline::const_iterator finishElement,
                                                 G4int                       index);
 
 
  BDSBeamlineElement* CreateBonusSectionStart(BDSBeamline const* const beamline);
  BDSBeamlineElement* CreateBonusSectionEnd(BDSBeamline const* const beamline);
 
  BDSBeamlineElement* CreateBridgeSection(BDSAcceleratorComponent*    defaultBridge,
                                          BDSBeamline::const_iterator element,
                                          BDSBeamline::const_iterator nextElement,
                                          BDSBeamline::const_iterator end,
                                          G4int&                      numberOfUniqueComponents,
                                          const G4int                 beamlineIndex,
                                          G4double                    crRadius);
 
  void PreviousAndNext(BDSBeamline::const_iterator it,
                       BDSBeamline::const_iterator startIt,
                       BDSBeamline::const_iterator endIt,
                       const BDSBeamlineElement*&  previous,
                       const BDSBeamlineElement*&  next) const;
 
  BDSAcceleratorComponent* CreateDefaultBridgeComponent();
 
  BDSAcceleratorComponent* CreateStraightBridgeComponent(G4double                    width,
                                                         G4int&                      numberOfUniqueComponents);
  
  BDSAcceleratorComponent* CreateAngledBridgeComponent(BDSBeamline::const_iterator element,
                                                       G4int&                      numberOfUniqueComponents,
                                                       G4double                    suggestedRadius);
 
  BDSBeamlineElement* CreateBridgeElementFromComponent(BDSAcceleratorComponent*    component,
                                                       BDSBeamline::const_iterator element,
                                                       BDSBeamline::const_iterator nextElement,
                                                       BDSBeamline::const_iterator end,
                                                       const G4int                 beamlineIndex);
 
  inline G4bool Angled(BDSBeamlineElement const* const element) const;
 
  inline G4bool TooShort(BDSBeamlineElement const* const element) const;
 
  G4double paddingLength;
  G4double defaultBridgeLength;
  G4double curvilinearRadius; 
  G4double radiusTolerance;
  G4double bonusChordLength;  
 
  BDSCurvilinearFactory* factory;
};
 
inline G4bool BDSCurvilinearBuilder::Angled(BDSBeamlineElement const* const element) const
{return BDS::IsFinite(element->GetAngle());}