00001 // 00002 // ******************************************************************** 00003 // * License and Disclaimer * 00004 // * * 00005 // * The Geant4 software is copyright of the Copyright Holders of * 00006 // * the Geant4 Collaboration. It is provided under the terms and * 00007 // * conditions of the Geant4 Software License, included in the file * 00008 // * LICENSE and available at http://cern.ch/geant4/license . These * 00009 // * include a list of copyright holders. * 00010 // * * 00011 // * Neither the authors of this software system, nor their employing * 00012 // * institutes,nor the agencies providing financial support for this * 00013 // * work make any representation or warranty, express or implied, * 00014 // * regarding this software system or assume any liability for its * 00015 // * use. Please see the license in the file LICENSE and URL above * 00016 // * for the full disclaimer and the limitation of liability. * 00017 // * * 00018 // * This code implementation is the result of the scientific and * 00019 // * technical work of the GEANT4 collaboration. * 00020 // * By using, copying, modifying or distributing the software (or * 00021 // * any work based on the software) you agree to acknowledge its * 00022 // * use in resulting scientific publications, and indicate your * 00023 // * acceptance of all terms of the Geant4 Software license. * 00024 // ******************************************************************** 00025 // 00026 // 00027 // $Id: G4tgrRotationMatrix.cc 69803 2013-05-15 15:24:50Z gcosmo $ 00028 // 00029 // 00030 // class G4tgrRotationMatrix 00031 00032 // History: 00033 // - Created. P.Arce, CIEMAT (November 2007) 00034 // ------------------------------------------------------------------------- 00035 00036 #include "G4tgrRotationMatrix.hh" 00037 00038 #include "G4SystemOfUnits.hh" 00039 #include "G4tgrRotationMatrixFactory.hh" 00040 #include "G4tgrUtils.hh" 00041 #include "G4tgrMessenger.hh" 00042 00043 // ------------------------------------------------------------------------- 00044 G4tgrRotationMatrix::G4tgrRotationMatrix() 00045 : theName("Rotation-Matrix"), theInputType(rm9) 00046 { 00047 } 00048 00049 00050 // ------------------------------------------------------------------------- 00051 G4tgrRotationMatrix::~G4tgrRotationMatrix() 00052 { 00053 } 00054 00055 00056 // ------------------------------------------------------------------------- 00057 G4tgrRotationMatrix::G4tgrRotationMatrix( const std::vector<G4String>& wl ) 00058 : theInputType(rm9) 00059 { 00060 theName = G4tgrUtils::GetString( wl[1] ); 00061 00062 switch( wl.size() ) 00063 { 00064 case 5: 00065 theInputType = rm3; 00066 break; 00067 case 8: 00068 theInputType = rm6; 00069 break; 00070 case 11: 00071 theInputType = rm9; 00072 break; 00073 default: 00074 G4Exception("G4tgrRotationMatrix::G4tgrRotationMatrix()", 00075 "InvalidMatrix", FatalException, 00076 "Input line must have 5, 8 or 11 words."); 00077 break; 00078 } 00079 00080 //-------- Fill matrix values 00081 size_t siz = wl.size() - 2; 00082 for( size_t ii = 0; ii < siz; ii++) 00083 { 00084 if( siz == 9 ) 00085 { 00086 theValues.push_back( G4tgrUtils::GetDouble( wl[ii+2] ) ); 00087 } 00088 else 00089 { 00090 theValues.push_back( G4tgrUtils::GetDouble( wl[ii+2] , deg ) ); 00091 } 00092 } 00093 #ifdef G4VERBOSE 00094 if( G4tgrMessenger::GetVerboseLevel() >= 2 ) 00095 { 00096 G4cout << " G4tgrRotationMatrix::G4tgrRotationMatrix() - Created: " 00097 << theName << G4endl; 00098 for( size_t ii = 0; ii < siz; ii++) 00099 { 00100 G4cout << " " << theValues[ii]; 00101 } 00102 G4cout << G4endl; 00103 } 00104 #endif 00105 } 00106 00107 00108 // ------------------------------------------------------------------------- 00109 std::ostream& operator<<(std::ostream& os, const G4tgrRotationMatrix& obj) 00110 { 00111 os << "G4tgrRotationMatrix= " << obj.theName 00112 << " InputTyep = " << obj.theInputType << " VALUES= "; 00113 00114 for( size_t ii = 0; ii < obj.theValues.size(); ii++ ) 00115 { 00116 os << obj.theValues[ii] << " "; 00117 } 00118 00119 os << G4endl; 00120 00121 return os; 00122 }