BDSCrystalInfo.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 BDSCRYSTALINFO_H
#define BDSCRYSTALINFO_H
 
#include "BDSCrystalType.hh"
#include "globals.hh"         // geant4 types / globals
 
#include "CLHEP/Units/SystemOfUnits.h"
 
class G4Material;
 
class BDSCrystalInfo
{
public:
  BDSCrystalInfo() = delete;
  
  BDSCrystalInfo(G4Material*     materialIn,
                 const G4String& dataIn,
                 BDSCrystalType  crystalTypeIn,
                 G4double        lengthXIn,
                 G4double        lengthYIn,
                 G4double        lengthZIn,
                 G4double        sizeAIn,
                 G4double        sizeBIn,
                 G4double        sizeCIn,
                 G4double        alphaIn,
                 G4double        betaIn,
                 G4double        gammaIn,
                 G4int           spaceGroupIn,
                 G4double        bendingAngleYAxisIn = 0,
                 G4double        bendingAngleZAxisIn = 0,
                 G4double        miscutAngleY        = 0);
 
  G4Material*    material; 
  G4String       data;     
  BDSCrystalType shape;    
  G4double       lengthX;  
  G4double       lengthY;  
  G4double       lengthZ;  
  G4double       sizeA;
  G4double       sizeB;
  G4double       sizeC;
  G4double       alpha;
  G4double       beta;
  G4double       gamma;
  G4int          spaceGroup;
  G4double       bendingAngleYAxis; 
  G4double       bendingAngleZAxis; 
  G4double       miscutAngleY;
 
  G4double BendingRadiusHorizontal() const {return BendingRadius(lengthZ, bendingAngleYAxis);}
  G4double BendingRadiusVertical()   const {return BendingRadius(lengthY, bendingAngleZAxis);}
 
private:
  G4double BendingRadius(G4double length,
                         G4double angle) const;
};
 
#endif