Geant4 Cross Reference |
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 // G4tgrPlaceDivRep implementation << 27 // 26 // 28 // Author: P.Arce, CIEMAT (November 2007) << 27 // $Id: G4tgrPlaceDivRep.cc 69803 2013-05-15 15:24:50Z gcosmo $ 29 // ------------------------------------------- << 28 // >> 29 // >> 30 // class G4tgrPlaceDivRep >> 31 >> 32 // History: >> 33 // - Created. P.Arce, CIEMAT (November 2007) >> 34 // ------------------------------------------------------------------------- 30 35 31 #include "G4tgrPlaceDivRep.hh" 36 #include "G4tgrPlaceDivRep.hh" 32 37 33 #include "G4SystemOfUnits.hh" 38 #include "G4SystemOfUnits.hh" 34 #include "G4tgrUtils.hh" 39 #include "G4tgrUtils.hh" 35 #include "G4tgrMessenger.hh" 40 #include "G4tgrMessenger.hh" 36 #include "G4tgrVolume.hh" 41 #include "G4tgrVolume.hh" 37 42 38 // ------------------------------------------- << 43 //------------------------------------------------------------- 39 G4tgrPlaceDivRep::G4tgrPlaceDivRep() 44 G4tgrPlaceDivRep::G4tgrPlaceDivRep() >> 45 : theNDiv(0), theWidth(0.), theAxis(kUndefined), >> 46 theOffset(0.), theDivType(DivByNdivAndWidth) 40 { 47 { 41 } 48 } 42 49 43 // ------------------------------------------- << 50 >> 51 //------------------------------------------------------------- 44 G4tgrPlaceDivRep::~G4tgrPlaceDivRep() 52 G4tgrPlaceDivRep::~G4tgrPlaceDivRep() 45 { 53 { 46 } 54 } 47 55 48 // ------------------------------------------- << 56 //------------------------------------------------------------- 49 G4tgrPlaceDivRep::G4tgrPlaceDivRep(const std:: << 57 G4tgrPlaceDivRep::G4tgrPlaceDivRep( const std::vector<G4String>& wl ) 50 { 58 { 51 theDivType = DivByNdivAndWidth; 59 theDivType = DivByNdivAndWidth; 52 60 53 // Name parent axis nrep width offset 61 // Name parent axis nrep width offset 54 G4tgrUtils::CheckWLsize(wl, 6, WLSIZE_GE, << 62 G4tgrUtils::CheckWLsize( wl, 6, WLSIZE_GE, 55 "G4tgrPlaceDivRep::G << 63 "G4tgrPlaceDivRep::G4tgrPlaceDivRep" ); 56 G4tgrUtils::CheckWLsize(wl, 7, WLSIZE_LE, << 64 G4tgrUtils::CheckWLsize( wl, 7, WLSIZE_LE, 57 "G4tgrPlaceDivRep::G << 65 "G4tgrPlaceDivRep::G4tgrPlaceDivRep" ); 58 << 66 59 theParentName = G4tgrUtils::GetString(wl[2]) << 67 theParentName = G4tgrUtils::GetString(wl[2]); 60 theAxis = BuildAxis(G4tgrUtils::GetStr << 68 theAxis = BuildAxis( G4tgrUtils::GetString(wl[3]) ); 61 theNDiv = G4tgrUtils::GetInt(wl[4]); << 69 theNDiv = G4tgrUtils::GetInt( wl[4] ); 62 theWidth = G4tgrUtils::GetDouble(wl[5]) << 70 theWidth = G4tgrUtils::GetDouble(wl[5])*mm; // check if it is deg 63 if(wl.size() == 7) << 71 if( wl.size() == 7 ) 64 { 72 { 65 theOffset = G4tgrUtils::GetDouble(wl[6]) * << 73 theOffset = G4tgrUtils::GetDouble(wl[6])*mm; // check if it is deg 66 } 74 } 67 else 75 else 68 { 76 { 69 theOffset = 0.; 77 theOffset = 0.; 70 } 78 } 71 79 72 #ifdef G4VERBOSE 80 #ifdef G4VERBOSE 73 if(G4tgrMessenger::GetVerboseLevel() >= 1) << 81 if( G4tgrMessenger::GetVerboseLevel() >= 1 ) 74 { 82 { 75 G4cout << " Created " << *this << G4endl; << 83 G4cout << " Created " << *this << G4endl; 76 } 84 } 77 #endif 85 #endif >> 86 78 } 87 } 79 88 80 // ------------------------------------------- << 89 81 EAxis G4tgrPlaceDivRep::BuildAxis(const G4Stri << 90 //------------------------------------------------------------- >> 91 EAxis G4tgrPlaceDivRep::BuildAxis( const G4String& axisName ) 82 { 92 { 83 if(axisName == "X") << 93 if( axisName == "X" ) { 84 { << 85 return kXAxis; 94 return kXAxis; 86 } << 95 } else if( axisName == "Y" ) { 87 else if(axisName == "Y") << 88 { << 89 return kYAxis; 96 return kYAxis; 90 } << 97 } else if( axisName == "Z" ) { 91 else if(axisName == "Z") << 92 { << 93 return kZAxis; 98 return kZAxis; 94 } << 99 } else if( axisName == "R" ) { 95 else if(axisName == "R") << 96 { << 97 return kRho; 100 return kRho; 98 } << 101 } else if( axisName == "PHI" ) { 99 else if(axisName == "PHI") << 100 { << 101 return kPhi; 102 return kPhi; 102 } 103 } 103 else 104 else 104 { 105 { 105 G4String ErrMessage = "Axis type not found << 106 G4String ErrMessage = "Axis type not found: " + axisName 106 ". Only valid axis a << 107 + ". Only valid axis are: X, Y, Z, R, PHI !"; 107 G4Exception("G4tgrVolumeDivision::GetRepli << 108 G4Exception("G4tgrVolumeDivision::GetReplicaAxis()", 108 FatalException, ErrMessage); << 109 "InvalidAxis", FatalException, ErrMessage); 109 } 110 } 110 return kXAxis; // to avoid warning errors << 111 return kXAxis; // to avoid warning errors 111 } 112 } 112 113 113 // ------------------------------------------- << 114 >> 115 // ------------------------------------------------------------------------- 114 std::ostream& operator<<(std::ostream& os, con 116 std::ostream& operator<<(std::ostream& os, const G4tgrPlaceDivRep& obj) 115 { 117 { 116 os << "G4tgrPlaceDivRep= in " << obj.thePare << 118 os << "G4tgrPlaceDivRep= in " << obj.theParentName 117 << " Width= " << obj.theWidth << " Axis= << 119 << " NDiv= " << obj.theNDiv << " Width= " << obj.theWidth 118 << " Offset= " << obj.theOffset << " DivT << 120 << " Axis= " << obj.theAxis << " Offset= " << obj.theOffset 119 << G4endl; << 121 << " DivType= " << obj.theDivType << G4endl; 120 122 121 return os; 123 return os; 122 } 124 } 123 125