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