Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/geometry/magneticfield/include/G4TQuadrupoleMagField.hh

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 ]

  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 // G4TMagFieldEquation
 27 //
 28 // Class description:
 29 //
 30 // Templated version of  motion of quadrupole magnetic field.
 31 // Enables testing of inlined code for field, equation, stepper, driver,
 32 // avoiding all virtual calls.
 33 //
 34 // --------------------------------------------------------------------
 35 // Created: Josh Xie  (Google Summer of Code 2014 )
 36 // Adapted from G4QuadrupoleMagField  ( of June 2014 .... )
 37 //
 38 
 39 #ifndef G4T_QUADRUPOLE_MAGFIELD_HH
 40 #define G4T_QUADRUPOLE_MAGFIELD_HH
 41 
 42 #include "G4ThreeVector.hh"
 43 #include "G4RotationMatrix.hh"
 44 #include "G4MagneticField.hh"
 45 
 46 static G4RotationMatrix IdentityMatrix; 
 47 
 48 class G4TQuadrupoleMagField : public G4MagneticField 
 49 {
 50   public: // with description
 51 
 52   G4TQuadrupoleMagField(G4double pGradient)
 53   {
 54     fGradient = pGradient ;
 55     fOrigin   = G4ThreeVector( 0.0, 0.0, 0.0) ;
 56     fpMatrix  = &IdentityMatrix;
 57   }
 58 
 59         G4TQuadrupoleMagField(G4double pGradient, 
 60                 G4ThreeVector pOrigin, 
 61                 G4RotationMatrix* pMatrix)
 62         {
 63             fGradient    = pGradient ;
 64             fOrigin      = pOrigin ;
 65             fpMatrix     = pMatrix ;
 66         }
 67 
 68         virtual ~G4TQuadrupoleMagField() {;}
 69 
 70         inline void GetFieldValue(const G4double y[7],
 71                 G4double B[3]     ) const
 72         {
 73             G4ThreeVector r_global = G4ThreeVector(
 74                     y[0] - fOrigin.x(), 
 75                     y[1] - fOrigin.y(), 
 76                     y[2] - fOrigin.z());
 77 
 78             G4ThreeVector r_local = G4ThreeVector(
 79                     fpMatrix->colX() * r_global,
 80                     fpMatrix->colY() * r_global,
 81                     fpMatrix->colZ() * r_global);
 82 
 83             G4ThreeVector B_local = G4ThreeVector(
 84                     fGradient * r_local.y(),
 85                     fGradient * r_local.x(),
 86                     0);
 87 
 88             G4ThreeVector B_global = G4ThreeVector(
 89                     fpMatrix->inverse().rowX() * B_local,
 90                     fpMatrix->inverse().rowY() * B_local,
 91                     fpMatrix->inverse().rowZ() * B_local);
 92 
 93             B[0] = B_global.x() ;
 94             B[1] = B_global.y() ;
 95             B[2] = B_global.z() ;
 96         }
 97 
 98         G4TQuadrupoleMagField* Clone() const
 99         {
100             //TODO: Can the fpMatrix be shared??
101             return new G4TQuadrupoleMagField(this->fGradient,
102                     this->fOrigin,
103                     this->fpMatrix);
104         }
105 
106     private:
107 
108         G4double          fGradient;
109         G4ThreeVector     fOrigin;
110         G4RotationMatrix* fpMatrix;
111 };
112 #endif
113 
114