Geant4  9.6.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4tgbPlaceParamCircle.cc
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
27 // $Id: G4tgbPlaceParamCircle.cc 69803 2013-05-15 15:24:50Z gcosmo $
28 //
29 //
30 // class G4tgbPlaceParamCircle
31 
32 // History:
33 // - Created. P.Arce, CIEMAT (November 2007)
34 // -------------------------------------------------------------------------
35 
36 #include "G4tgbPlaceParamCircle.hh"
38 #include "G4ThreeVector.hh"
39 #include "G4RotationMatrix.hh"
40 #include "G4VPhysicalVolume.hh"
41 #include "G4tgrUtils.hh"
42 #include "G4tgrMessenger.hh"
43 
44 // -------------------------------------------------------------------------
46 {
47 }
48 
49 
50 // -------------------------------------------------------------------------
53  : G4tgbPlaceParameterisation(tgrParam)
54 {
55  //---- Get translation and rotation
56  if( tgrParam->GetParamType() == "CIRCLE" )
57  {
58  CheckNExtraData( tgrParam, 7, WLSIZE_EQ, "G4tgbPlaceParamCircle:");
59  theCircleAxis = G4ThreeVector( tgrParam->GetExtraData()[4],
60  tgrParam->GetExtraData()[5],
61  tgrParam->GetExtraData()[6] );
62 
63  G4ThreeVector zaxis(0.,0.,-1.);
64  if( zaxis.cross(theCircleAxis).mag() > 1.E-6 )
65  {
66  theDirInPlane = zaxis.cross(theCircleAxis);
67  }
68  else
69  {
70  theDirInPlane = theCircleAxis.cross(G4ThreeVector(0.,-1.,0.));
71  }
72  theAxis = kZAxis;
73  }
74  else
75  {
76  CheckNExtraData( tgrParam, 4, WLSIZE_EQ, "G4tgbPlaceParamCircle:");
77  if( tgrParam->GetParamType() == "CIRCLE_XY" ) {
78  theCircleAxis = G4ThreeVector(0.,0.,1.);
79  theDirInPlane = G4ThreeVector(1.,0.,0.);
80  theAxis = kZAxis;
81  } else if( tgrParam->GetParamType() == "CIRCLE_XZ" ) {
82  theCircleAxis = G4ThreeVector(0.,1.,0.);
83  theDirInPlane = G4ThreeVector(1.,0.,0.);
84  theAxis = kYAxis;
85  } else if( tgrParam->GetParamType() == "CIRCLE_YZ" ) {
86  theCircleAxis = G4ThreeVector(1.,0.,0.);
87  theDirInPlane = G4ThreeVector(0.,1.,0.);
88  theAxis = kXAxis;
89  }
90  }
91 
92  if( theCircleAxis.mag() == 0. )
93  {
94  G4Exception("G4tgbPlaceParamCircle::G4tgbPlaceParamCircle()",
95  "InvalidSetup", FatalException, "Circle axis is zero !");
96  }
97  theCircleAxis /= theCircleAxis.mag();
98 
99  theAxis = kZAxis;
100 
101  theNCopies = G4int(tgrParam->GetExtraData()[0]);
102  theStep = tgrParam->GetExtraData()[1];
103  theOffset = tgrParam->GetExtraData()[2];
104  theRadius = tgrParam->GetExtraData()[3];
105 
106 #ifdef G4VERBOSE
108  {
109  G4cout << " G4tgbPlaceParamCircle::G4tgbPlaceParamCircle():" << G4endl
110  << " param type " << tgrParam->GetParamType() << G4endl
111  << " no copies - " << theNCopies << G4endl
112  << " step - " << theStep << G4endl
113  << " offset - " << theOffset << G4endl
114  << " radius - " << theRadius << G4endl
115  << " circle axis - " << theCircleAxis << G4endl
116  << " dir in plane - " << theDirInPlane << G4endl;
117  }
118 #endif
119 }
120 
121 
122 // -------------------------------------------------------------------------
124 ComputeTransformation(const G4int copyNo, G4VPhysicalVolume *physVol) const
125 {
126  G4double posi = theOffset + copyNo*theStep;
127  G4ThreeVector origin = theDirInPlane * theRadius;
128  origin.rotate( posi, theCircleAxis );
129 
130  //----- Calculate rotation matrix (so that all volumes point to the centre)
131  G4RotationMatrix rm;
132  rm.rotate( -posi, theCircleAxis );
133 
134  //----- Set translation and rotation
135  physVol->SetTranslation(origin);
136  G4RotationMatrix* pvRm = physVol->GetRotation();
137  if( pvRm == 0 )
138  {
139  pvRm = new G4RotationMatrix;
140  }
141  *pvRm = *theRotationMatrix * rm;
142  physVol->SetRotation(pvRm);
143  physVol->SetCopyNo( copyNo );
144 
145 #ifdef G4VERBOSE
147  {
148  G4cout << " G4tgbPlaceParamCircle::ComputeTransformation():"
149  << physVol->GetName() << G4endl
150  << " no copies - " << theNCopies << G4endl
151  << " centre - " << origin << G4endl
152  << " rotation-matrix - " << *pvRm << G4endl;
153  }
154 #endif
155 }