Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/persistency/ascii/src/G4tgbPlaceParamCircle.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /persistency/ascii/src/G4tgbPlaceParamCircle.cc (Version 11.3.0) and /persistency/ascii/src/G4tgbPlaceParamCircle.cc (Version 10.7.p1)


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