| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 // Implementation of G4LogicalCrystalVolume
27 //
28 // 21-04-16, created by E.Bagli
29 // -------------------------------------------
30
31 #include "G4LogicalCrystalVolume.hh"
32 #include "G4ExtendedMaterial.hh"
33 #include "G4CrystalExtension.hh"
34 #include "G4VMaterialExtension.hh"
35
36 std::vector<G4LogicalVolume*> G4LogicalCrystal
37
38 // -------------------------------------------
39
40 G4LogicalCrystalVolume::
41 G4LogicalCrystalVolume(G4VSolid* pSolid, G4Ext
42 const G4String& name, G
43 G4VSensitiveDetector* p
44 G4UserLimits* pULimits,
45 G4int h, G4int k, G4int
46 : G4LogicalVolume(pSolid,pMaterial,name,pField
47 {
48 SetMillerOrientation(h, k, l, rot);
49 fLCVvec.push_back(this);
50 }
51
52 // -------------------------------------------
53
54 G4LogicalCrystalVolume::~G4LogicalCrystalVolum
55 {
56 fLCVvec.erase( std::remove(fLCVvec.begin(),f
57 fLCVvec.end() );
58 }
59
60 // -------------------------------------------
61
62 G4bool G4LogicalCrystalVolume::IsLattice(G4Log
63 {
64 return std::find(fLCVvec.begin(), fLCVvec.en
65 }
66
67 // -------------------------------------------
68
69 const G4CrystalExtension* G4LogicalCrystalVolu
70 {
71 return dynamic_cast<G4CrystalExtension*>(
72 dynamic_cast<G4ExtendedMaterial*>(Get
73 ->RetrieveExte
74 }
75
76 // -------------------------------------------
77
78 const G4ThreeVector& G4LogicalCrystalVolume::G
79 {
80 return GetCrystal()->GetUnitCell()->GetBasis
81 }
82
83 // -------------------------------------------
84
85 void G4LogicalCrystalVolume::SetMillerOrientat
86
87
88
89 {
90 // Align Miller normal vector (hkl) with +Z
91
92 if (verboseLevel != 0)
93 {
94 G4cout << "G4LatticePhysical::SetMillerOri
95 << k << " " << l << ", " << rot/CLH
96 }
97
98 hMiller = h;
99 kMiller = k;
100 lMiller = l;
101 fRot = rot;
102
103 G4ThreeVector norm = (h*GetBasis(0)+k*GetBa
104
105 if (verboseLevel>1) G4cout << " norm = " <<
106
107 // Aligns geometry +Z axis with lattice (hk
108 fOrient = G4RotationMatrix::IDENTITY;
109 fOrient.rotateZ(rot).rotateY(norm.theta()).
110 fInverse = fOrient.inverse();
111
112 if (verboseLevel>1) G4cout << " fOrient = "
113
114 // FIXME: Is this equivalent to (phi,theta
115 }
116
117 // -------------------------------------------
118
119 // Rotate input vector between lattice and sol
120
121 const G4ThreeVector&
122 G4LogicalCrystalVolume::RotateToLattice(G4Thre
123 {
124 return dir.transform(fOrient);
125 }
126
127 const G4ThreeVector&
128 G4LogicalCrystalVolume::RotateToSolid(G4ThreeV
129 {
130 return dir.transform(fInverse);
131 }
132
133 // -------------------------------------------
134