BDSCrystalFactory.hh

/* 
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/>.
*/
#ifndef BDSCRYSTALFACTORY_H
#define BDSCRYSTALFACTORY_H
 
#include "BDSCrystal.hh"
#include "BDSCrystalInfo.hh"
#include "BDSFactoryBase.hh"
 
#include "globals.hh"         // geant4 globals / types
#include "G4RotationMatrix.hh"
 
#include <set>
 
class BDSExtent;
class G4LogicalVolume;
class G4Material;
class G4PVPlacement;
class G4UserLimits;
class G4VSolid;
 
class BDSCrystalFactory: public BDSFactoryBase
{
public:
  BDSCrystalFactory();
  virtual ~BDSCrystalFactory(){;}
 
  BDSCrystal* CreateCrystal(const G4String&       nameIn,
                            const BDSCrystalInfo* recipe);
 
private:
  void CleanUp();
  
  void CommonConstruction(const G4String&       nameIn,
                          const BDSCrystalInfo* recipe);
  
  BDSCrystal* BuildCrystalObject(const BDSCrystalInfo*  recipe,
                                 const BDSExtent& extent);
  
  void SetVisAttributes();                    
  void SetUserLimits(G4double length); 
 
  BDSCrystal* CreateCrystalBox(const G4String&       nameIn,
                               const BDSCrystalInfo* recipe);
 
  BDSCrystal* CreateCrystalCylinder(const G4String&       nameIn,
                                    const BDSCrystalInfo* recipe);
 
  BDSCrystal* CreateCrystalTorus(const G4String&       nameIn,
                                 const BDSCrystalInfo* recipe);
 
  void CalculateSolidAngles(G4double  bendingAngle,
                            G4double& startAngle,
                            G4double& sweepAngle) const;
 
  BDSExtent CalculateExtents(G4double bendingAngle,
                             G4double xBendingRadius,
                             G4double xThickness,
                             const BDSCrystalInfo* recipe) const;
  
  const G4double   maxStepFactor;      
  const G4int      nPoints;            
  
  G4VSolid*         crystalSolid;
  G4LogicalVolume*  crystalLV;
  G4ThreeVector     placementOffset;
  G4RotationMatrix* placementRotation;
};
 
#endif