Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/geometry/magneticfield/src/G4DELPHIMagField.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 /geometry/magneticfield/src/G4DELPHIMagField.cc (Version 11.3.0) and /geometry/magneticfield/src/G4DELPHIMagField.cc (Version 5.2.p1)


  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 // G4DELPHIMagField implementation             << 
 27 //                                             << 
 28 // Created: V.Grichine - 03.02.1997            << 
 29 // ------------------------------------------- << 
 30                                                << 
 31 #include "G4DELPHIMagField.hh"                     23 #include "G4DELPHIMagField.hh"
 32 #include "globals.hh"                              24 #include "globals.hh"
 33 #include "G4SystemOfUnits.hh"                  <<  25 #include "geomdefs.hh"
 34                                                    26 
 35 G4DELPHIMagField::G4DELPHIMagField() = default << 
 36                                                    27 
 37 ////////////////////////////////////////////// << 
 38                                                    28 
 39 G4DELPHIMagField::~G4DELPHIMagField() = defaul <<  29 G4DELPHIMagField::G4DELPHIMagField()
                                                   >>  30 {
                                                   >>  31    ;
                                                   >>  32 }
 40                                                    33 
 41 ////////////////////////////////////////////// <<  34 ////////////////////////////////////////////////////////////////////////
 42                                                    35 
 43 G4Field* G4DELPHIMagField::Clone() const       <<  36 G4DELPHIMagField::~G4DELPHIMagField()
 44 {                                                  37 {
 45    return new G4DELPHIMagField;                <<  38    ;
 46 }                                                  39 }
 47                                                    40 
 48 //////////////////////////////////////////////     41 ///////////////////////////////////////////////////////////////////////
 49                                                    42 
 50 void G4DELPHIMagField::GetFieldValue( const G4 <<  43 
 51                                             G4 <<  44 void
                                                   >>  45    G4DELPHIMagField::GetFieldValue( const G4double yTrack[7],
                                                   >>  46                                           G4double B[3]         ) const 
 52 {                                                  47 {
 53    G4int i, n = 8 ;                                48    G4int i, n = 8 ;
 54    G4double a = 0.001 ;   // mm -> m               49    G4double a = 0.001 ;   // mm -> m 
 55    G4double x = a*yTrack[0], y = a*yTrack[1],      50    G4double x = a*yTrack[0], y = a*yTrack[1], z = a*yTrack[2] ;
 56    G4double x2 = x*x, y2 = y*y, z2 = z*z, r2 =     51    G4double x2 = x*x, y2 = y*y, z2 = z*z, r2 = x2 + y2 ;
 57    G4double r4 = r2*r2, z4 = z2*z2, r6 = r4*r2     52    G4double r4 = r2*r2, z4 = z2*z2, r6 = r4*r2, z6 = z4*z2 ;
 58    G4double r8 = r4*r4, z8 = z4*z4, r10 = r8*r     53    G4double r8 = r4*r4, z8 = z4*z4, r10 = r8*r2, z10 = z8*z2 ;
 59    G4double rz = z*std::sqrt(r2), r = std::sqr <<  54    G4double rz = z*sqrt(r2), r = sqrt(r2+a*a) ;
 60    G4double Br ;                                   55    G4double Br ;
 61    G4double P[8], Q[8] ;                           56    G4double P[8], Q[8] ; 
 62    static G4ThreadLocal G4double c[8] =        <<  57    static G4double c[8] = {
 63                           {                    << 
 64                             -9.26e-5, -3.51e-5     58                             -9.26e-5, -3.51e-5, 2.94e-6, -1.10e-6, 
 65                              6.25e-8, -1.77e-8     59                              6.25e-8, -1.77e-8, -6.88e-10, -7.52e-11 
 66                           } ;                      60                           } ;
 67    P[0] = 2*rz ;                                   61    P[0] = 2*rz ;
 68    P[1] = 4*rz*(r2 - 4.0*z2/3.0) ;                 62    P[1] = 4*rz*(r2 - 4.0*z2/3.0) ;
 69    P[2] = 6*rz*(r4 - 4*r2*z2 + 1.6*z4) ;           63    P[2] = 6*rz*(r4 - 4*r2*z2 + 1.6*z4) ;
 70    P[3] = 8*rz*(r6 - 8*r4*z2 + 9.6*r2*z4 - 64.     64    P[3] = 8*rz*(r6 - 8*r4*z2 + 9.6*r2*z4 - 64.0*z6/35.0) ;
 71    P[4] = 10*rz*(r8 - 40.0*r6*z2/3.0 + 32*r4*z <<  65    P[4] = 10*rz*(r8 - 40.0*r6*z2/3.0 + 32*r4*z4 - 128.0*r2*z6/7.0 + 128.0*z8/63.0);
 72                     - 128.0*r2*z6/7.0 + 128.0* << 
 73    P[5] = 0 ;                                      66    P[5] = 0 ;
 74    P[6] = 0 ;                                      67    P[6] = 0 ;
 75    P[7] = 0 ;                                      68    P[7] = 0 ;
 76                                                    69    
 77    Q[0] = r2 - 2*z2 ;                              70    Q[0] = r2 - 2*z2 ;
 78    Q[1] = r4 - 8*r2*z2 + 8.0*z4/3.0 ;              71    Q[1] = r4 - 8*r2*z2 + 8.0*z4/3.0 ;
 79    Q[2] = r6 - 18*r4*z2 + 24*r2*z4 - 3.2*z6 ;      72    Q[2] = r6 - 18*r4*z2 + 24*r2*z4 - 3.2*z6 ;
 80    Q[3] = r8 - 32*r6*z2 + 96*r4*z4 - 51.2*r2*z     73    Q[3] = r8 - 32*r6*z2 + 96*r4*z4 - 51.2*r2*z6 +128.0*z8/35.0 ;
 81    Q[4] = r10 - 50*r8*z2 + 800.0*r6*z4/3.0 - 3     74    Q[4] = r10 - 50*r8*z2 + 800.0*r6*z4/3.0 - 320*r4*z6 
 82           + 640.0*r2*z8/7.0 - 256.0*z10/63.0 ;     75           + 640.0*r2*z8/7.0 - 256.0*z10/63.0 ;
 83    Q[5] = 0 ;                                      76    Q[5] = 0 ;
 84    Q[6] = 0 ;                                      77    Q[6] = 0 ;
 85    Q[7] = 0 ;                                      78    Q[7] = 0 ;
 86                                                    79    
 87    Br = 0 ;                                        80    Br = 0 ;
 88    B[2] = 1.2*tesla  ;    //  the principal Bz     81    B[2] = 1.2*tesla  ;    //  the principal Bz value of DELPHI detector
 89    for(i=0; i<n; ++i)                          <<  82    for(i=0;i<n;i++)
 90    {                                               83    {
 91       Br += c[i]*P[i] ;                            84       Br += c[i]*P[i] ;
 92       B[2] += c[i]*Q[i] ;                          85       B[2] += c[i]*Q[i] ;
 93    }                                               86    }
 94    B[0] = Br*x/r ;                                 87    B[0] = Br*x/r ;
 95    B[1] = Br*y/r ;                                 88    B[1] = Br*y/r ;
 96    return ;                                        89    return ;
 97 }                                                  90 }
 98                                                    91 
 99 // -------------------------------------------     92 // -----------------------------------------------------------------
                                                   >>  93 
100                                                    94