Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/geometry/magneticfield/include/G4MagIntegratorDriver.icc

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/include/G4MagIntegratorDriver.icc (Version 11.3.0) and /geometry/magneticfield/include/G4MagIntegratorDriver.icc (Version 3.1)


                                                   >>   1 // This code implementation is the intellectual property of
                                                   >>   2 // the GEANT4 collaboration.
  1 //                                                  3 //
  2 // ******************************************* <<   4 // By copying, distributing or modifying the Program (or any work
  3 // * License and Disclaimer                    <<   5 // based on the Program) you indicate your acceptance of this statement,
  4 // *                                           <<   6 // and all its terms.
  5 // * The  Geant4 software  is  copyright of th << 
  6 // * the Geant4 Collaboration.  It is provided << 
  7 // * conditions of the Geant4 Software License << 
  8 // * LICENSE and available at  http://cern.ch/ << 
  9 // * include a list of copyright holders.      << 
 10 // *                                           << 
 11 // * Neither the authors of this software syst << 
 12 // * institutes,nor the agencies providing fin << 
 13 // * work  make  any representation or  warran << 
 14 // * regarding  this  software system or assum << 
 15 // * use.  Please see the license in the file  << 
 16 // * for the full disclaimer and the limitatio << 
 17 // *                                           << 
 18 // * This  code  implementation is the result  << 
 19 // * technical work of the GEANT4 collaboratio << 
 20 // * By using,  copying,  modifying or  distri << 
 21 // * any work based  on the software)  you  ag << 
 22 // * use  in  resulting  scientific  publicati << 
 23 // * acceptance of all terms of the Geant4 Sof << 
 24 // ******************************************* << 
 25 //                                                  7 //
 26 // G4MagInt_Driver inline methods implementati <<   8 // $Id: G4MagIntegratorDriver.icc,v 1.5 2000/11/01 15:15:50 gcosmo Exp $
                                                   >>   9 // GEANT4 tag $Name: geant4-03-01 $
 27 //                                                 10 //
 28 // V.Grichine, 07.10.1996 - Created            << 
 29 // ------------------------------------------- << 
 30                                                << 
 31 inline                                         << 
 32 G4double G4MagInt_Driver::                     << 
 33 AdvanceChordLimited(G4FieldTrack& track, G4dou << 
 34                     G4double epsStep, G4double << 
 35 {                                              << 
 36   return ChordFinderDelegate::AdvanceChordLimi << 
 37                                                << 
 38 }                                              << 
 39                                                << 
 40 inline                                         << 
 41 void G4MagInt_Driver::OnStartTracking()        << 
 42 {                                              << 
 43   ChordFinderDelegate::ResetStepEstimate();    << 
 44 }                                              << 
 45                                                    11 
 46 inline                                             12 inline
 47 G4double G4MagInt_Driver::GetHmin() const          13 G4double G4MagInt_Driver::GetHmin() const
 48 {                                                  14 {
 49   return fMinimumStep;                         <<  15       return hminimum_val;
 50 }                                                  16 } 
 51                                                    17 
 52 inline                                             18 inline
 53 G4double G4MagInt_Driver::Hmin() const             19 G4double G4MagInt_Driver::Hmin() const
 54 {                                                  20 {
 55   return fMinimumStep;                         <<  21       return hminimum_val;
 56 }                                                  22 }
 57                                                    23 
 58 inline                                             24 inline
 59 G4double G4MagInt_Driver::GetSafety() const        25 G4double G4MagInt_Driver::GetSafety() const
 60 {                                                  26 {
 61   return safety;                               <<  27       return safety;
 62 }                                                  28 }
 63                                                    29 
 64 inline                                             30 inline
 65 G4double G4MagInt_Driver::GetPshrnk() const        31 G4double G4MagInt_Driver::GetPshrnk() const
 66 {                                                  32 {
 67   return pshrnk;                               <<  33       return pshrnk;
 68 }                                                  34 } 
 69                                                    35 
 70 inline                                             36 inline
 71 G4double G4MagInt_Driver::GetPgrow() const         37 G4double G4MagInt_Driver::GetPgrow() const
 72 {                                                  38 {
 73   return pgrow;                                <<  39       return pgrow;
 74 }                                                  40 }
 75                                                    41  
 76 inline                                             42 inline
 77 G4double G4MagInt_Driver::GetErrcon() const        43 G4double G4MagInt_Driver::GetErrcon() const
 78 {                                                  44 {
 79   return errcon;                               <<  45       return errcon;
 80 }                                                  46 }
 81                                                    47 
 82 inline                                             48 inline
 83 void G4MagInt_Driver::SetHmin(G4double newval)     49 void G4MagInt_Driver::SetHmin(G4double newval)
 84 {                                                  50 {
 85   fMinimumStep = newval;                       <<  51       hminimum_val = newval;
 86 }                                                  52 } 
 87                                                    53 
 88 inline                                             54 inline
 89 G4double G4MagInt_Driver::ComputeAndSetErrcon(     55 G4double G4MagInt_Driver::ComputeAndSetErrcon()
 90 {                                                  56 {
 91   errcon = std::pow(max_stepping_increase/GetS <<  57       errcon = pow(max_stepping_increase/GetSafety(),1.0/GetPgrow());
 92   return errcon;                               <<  58       return errcon;
 93 }                                                  59 } 
 94                                                    60 
 95 inline                                             61 inline
 96 void G4MagInt_Driver::ReSetParameters(G4double     62 void G4MagInt_Driver::ReSetParameters(G4double new_safety)
 97 {                                                  63 {
 98   safety = new_safety;                         <<  64       safety = new_safety;
 99   pshrnk = -1.0 / pIntStepper->IntegratorOrder <<  65       pshrnk = -1.0 / pIntStepper->IntegratorOrder();
100   pgrow  = -1.0 / (1.0 + pIntStepper->Integrat <<  66       pgrow  = -1.0 / (1.0 + pIntStepper->IntegratorOrder());
101   ComputeAndSetErrcon();                       <<  67       ComputeAndSetErrcon();
102 }                                                  68 }
103                                                    69 
104 inline                                             70 inline
105 void G4MagInt_Driver::SetSafety(G4double val)      71 void G4MagInt_Driver::SetSafety(G4double val)
106 {                                                  72 { 
107   safety = val;                                <<  73       safety=val;
108   ComputeAndSetErrcon();                       <<  74       ComputeAndSetErrcon();
109 }                                                  75 }
110                                                    76 
111 inline                                             77 inline
112 void G4MagInt_Driver::SetPgrow(G4double  val)      78 void G4MagInt_Driver::SetPgrow(G4double  val)
113 {                                                  79 { 
114   pgrow = val;                                 <<  80       pgrow=val;
115   ComputeAndSetErrcon();                       <<  81       ComputeAndSetErrcon(); 
116 }                                                  82 }
117                                                    83 
118 inline                                             84 inline
119 void G4MagInt_Driver::SetErrcon(G4double val)      85 void G4MagInt_Driver::SetErrcon(G4double val)
120 {                                                  86 { 
121   errcon = val;                                <<  87       errcon=val;
122 }                                                  88 }
123                                                    89 
124 inline                                             90 inline
125 G4int G4MagInt_Driver::GetMaxNoSteps() const   <<  91 void G4MagInt_Driver::RenewStepperAndAdjust(G4MagIntegratorStepper *pItsStepper)
126 {                                              <<  92 {  
127   return fMaxNoSteps;                          <<  93       pIntStepper = pItsStepper; 
                                                   >>  94       ReSetParameters();
128 }                                                  95 }
129                                                    96 
130 inline                                             97 inline
131 void G4MagInt_Driver::SetMaxNoSteps(G4int val) <<  98 void G4MagInt_Driver::SetChargeMomentumMass(G4double particleCharge, // e+
132 {                                              <<  99                                 G4double MomentumXc,
133   fMaxNoSteps = val;                           << 100                                 G4double Mass )
                                                   >> 101 { 
                                                   >> 102       pIntStepper->GetEquationOfMotion()
                                                   >> 103                  ->SetChargeMomentumMass(particleCharge, MomentumXc, Mass); 
134 }                                                 104 }
135                                                   105 
136 inline                                            106 inline
137 G4int G4MagInt_Driver::GetVerboseLevel() const << 107 const G4MagIntegratorStepper* G4MagInt_Driver::GetStepper() const
138 {                                                 108 {
139   return fVerboseLevel;                        << 109   return pIntStepper;
140 }                                              << 110 }
141                                                   111 
142 inline                                         << 112 inline
143 void G4MagInt_Driver::SetVerboseLevel(G4int ne << 113 G4int G4MagInt_Driver::GetMaxNoSteps() const
144 {                                                 114 {
145   fVerboseLevel = newLevel;                    << 115   return fMaxNoSteps;
146 }                                                 116 }
147                                                   117 
148 inline                                            118 inline
149 G4double G4MagInt_Driver::GetSmallestFraction( << 119 void G4MagInt_Driver::SetMaxNoSteps(G4int val)
150 {                                                 120 {
151   return fSmallestFraction;                    << 121   fMaxNoSteps= val;
                                                   >> 122 }
                                                   >> 123 
                                                   >> 124 inline
                                                   >> 125 void G4MagInt_Driver::GetDerivatives(const G4FieldTrack y_curr,  // const, INput
                                                   >> 126                    G4double     dydx[])  // OUTput
                                                   >> 127 { 
                                                   >> 128   G4double  tmpValArr[G4FieldTrack::ncompSVEC];
                                                   >> 129   y_curr.DumpToArray( tmpValArr  );
                                                   >> 130   pIntStepper -> RightHandSide( tmpValArr , dydx );
152 }                                                 131 }
153                                                   132