Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer << 3 // * DISCLAIMER * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers * 6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which * 7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: * 8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license * 9 // * include a list of copyright holders. << 10 // * 9 // * * 11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file << 14 // * use. * 16 // * for the full disclaimer and the limitatio << 17 // * 15 // * * 18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the * 19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. * 20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work * 21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this * 22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. * 23 // * acceptance of all terms of the Geant4 Sof << 24 // ******************************************* 21 // ******************************************************************** 25 // 22 // 26 // G4QuadrupoleMagField implementation << 27 // 23 // 28 // 03.02.1997, V.Grichine - Created << 24 // $Id: G4QuadrupoleMagField.cc,v 1.3 2003/10/31 14:35:55 gcosmo Exp $ 29 // 11.05.2012, B.Riese - Allow displaced orig << 25 // GEANT4 tag $Name: geant4-07-01 $ >> 26 // 30 // ------------------------------------------- 27 // ------------------------------------------------------------------- 31 28 32 #include "G4QuadrupoleMagField.hh" 29 #include "G4QuadrupoleMagField.hh" 33 #include "G4RotationMatrix.hh" << 34 << 35 namespace << 36 { << 37 G4RotationMatrix IdentityMatrix; << 38 } << 39 30 40 G4QuadrupoleMagField::G4QuadrupoleMagField(G4d 31 G4QuadrupoleMagField::G4QuadrupoleMagField(G4double pGradient) 41 { 32 { 42 fGradient = pGradient; << 43 fpMatrix = &IdentityMatrix; << 44 } << 45 << 46 // ------------------------------------------- << 47 << 48 G4QuadrupoleMagField::G4QuadrupoleMagField(G4d << 49 con << 50 G4R << 51 { << 52 fGradient = pGradient ; 33 fGradient = pGradient ; 53 fOrigin = pOrigin ; << 54 fpMatrix = pMatrix ; << 55 } 34 } 56 35 57 // ------------------------------------------- << 36 ///////////////////////////////////////////////////////////////////////// 58 37 59 G4Field* G4QuadrupoleMagField::Clone() const << 38 G4QuadrupoleMagField::~G4QuadrupoleMagField() 60 { 39 { 61 return new G4QuadrupoleMagField(fGradient, << 62 } 40 } 63 41 64 // ------------------------------------------- << 42 //////////////////////////////////////////////////////////////////////// 65 << 66 G4QuadrupoleMagField::~G4QuadrupoleMagField() << 67 43 68 // ------------------------------------------- << 44 void G4QuadrupoleMagField::GetFieldValue( const G4double y[7], 69 << 45 G4double B[3] ) const 70 void G4QuadrupoleMagField::GetFieldValue( cons << 71 << 72 { 46 { 73 // with displaced origin and rotation << 47 // G4double fGradient = 0.001 ; // Tesla/mm 74 << 48 B[0] = fGradient*y[1] ; 75 G4ThreeVector r_global = G4ThreeVector(y[0] << 49 B[1] = fGradient*y[0] ; 76 y[1] << 50 B[2] = 0 ; 77 y[2] << 78 << 79 const G4ThreeVector r_local = (*fpMatrix) * << 80 const G4ThreeVector B_local( fGradient * r_l << 81 fGradient * r_l << 82 const G4ThreeVector B_global = fpMatrix->inv << 83 << 84 B[0] = B_global.x(); << 85 B[1] = B_global.y(); << 86 B[2] = B_global.z(); << 87 } 51 } 88 52