66 for(
auto anAtomPos:
GetAtomBase(anElement)->GetPos())
76 SF +=
G4complex(AFF * GFS.real(),AFF * GFS.imag());
106 for (
size_t i=0; i<6; i++) {
107 for (
size_t j=0; j<6; j++) {
124 +
" is not registered." ;
135 std::vector<G4ThreeVector>
pos;
139 vecout.insert(std::end(vecout), std::begin(pos), std::end(pos));
147 std::vector<G4ThreeVector>
pos;
152 vecout.insert(std::end(vecout), std::begin(pos), std::end(pos));
void AddAtomBase(const G4Element *anElement, G4CrystalAtomBase *aBase)
static const G4double pos
std::map< const G4Element *, G4CrystalAtomBase * > theCrystalAtomBaseMap
G4bool FillAtomicPos(G4ThreeVector &pos, std::vector< G4ThreeVector > &vecout)
const G4String & GetName() const
ReducedElasticity fElReduced
G4complex ComputeStructureFactorGeometrical(G4int h, G4int k, G4int l)
const G4ElementVector * GetElementVector() const
G4CrystalExtension(G4Material *, const G4String &name="crystal")
G4bool GetAtomPos(const G4Element *anEl, std::vector< G4ThreeVector > &vecout)
std::complex< G4double > G4complex
const XML_Char int const XML_Char * value
std::vector< G4ThreeVector > GetPos()
void SetElReduced(const ReducedElasticity &mat)
G4CrystalUnitCell * theUnitCell
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
void SetCpq(G4int p, G4int q, G4double value)
G4complex ComputeStructureFactor(G4double kScatteringVector, G4int h, G4int k, G4int l)
static constexpr double pi
G4CrystalAtomBase * GetAtomBase(const G4Element *anElement)