Geant4 Cross Reference |
1 # Category field History << 1 ------------------------------------------------------------------- 2 << 3 See `CONTRIBUTING.rst` for details of **requir << 4 which **must** added in reverse chronological << 5 It must **not** be used as a substitute for wr << 6 << 7 ---------------------------------------------- << 8 << 9 ## 2024-11-26 Ivana Hrivnacova (field-V11-02-0 << 10 - Better names for the G4EquationType enum con << 11 << 12 ## 2024-11-04 Philippe Canal (field-V11-02-05) << 13 - Improve stability of static initialization f << 14 StepperDescription() in G4DormandPrince745. << 15 << 16 ## 2024-09-10 Ivana Hrivnacova (field-V11-02-0 << 17 - Coverity fixes: << 18 - Use std::move, const auto& instead of auto << 19 << 20 ## 2024-06-20 Ivana Hrivnacova (field-V11-02-0 << 21 - Fix in G4FieldBuilder for local fields: << 22 apply the local fields in the order of user << 23 << 24 ## 2024-04-29 Ivana Hrivnacova (field-V11-02-0 << 25 - Added classes for automatic field constructi << 26 - G4FieldBuilder[Messenger], G4FieldParamete << 27 - No changes in other Geant4 classes << 28 << 29 ## 2024-04-17 Ben Morgan (field-V11-02-01) << 30 - Fix maybe-uninitialized warnings found in AT << 31 - Reported in internal ATLASSIM-6058 ticket. << 32 << 33 ## 2024-04-02 Gabriele Cosmo (field-V11-02-00) << 34 - Applied trivial clang-tidy fixes to classes, << 35 constructor/destructor, auto and using state << 36 - Minor cleanup in QSS classes. << 37 << 38 ## 2023-11-03 Gabriele Cosmo (field-V11-01-06) << 39 - Reinstated default DormandPrince745 stepper. << 40 << 41 ## 2023-10-27 Lucio Santi (field-V11-01-05) << 42 - Improvements on the QSS Stepper << 43 The main QSS integration loop is interrupted << 44 exceeds a predefined threshold (set to 1000 << 45 << 46 ## 2023-06-14 Gabriele Cosmo (field-V11-01-04) << 47 - Applied clang-tidy fixes (readability, moder << 48 << 49 ## 2023-06-12 John Apostolakis (field-V11-01-0 << 50 - Inclusion of QSS integration method. << 51 Most recent changes: << 52 * Fixes for issue with SetPrecision (used un << 53 * Updated values of parameters in G4Propagat << 54 * Withdrew interim changes that made QSS the << 55 << 56 ## 2023-05-26 Lucio Santi << 57 - Fixes for open issues: << 58 * default values for QSS made stricter dQRel << 59 to improve behaviour in full experiment te << 60 * implemented missing DistChord method -- ne << 61 << 62 ## 2023-01-31 John Apostolakis << 63 - Initial inclusion of QSS integration method. << 64 Renamed classes to include 'G4' prefix, and << 65 Added Geant4 disclaimer / license to files. << 66 Removed QSSmessenger from RunManager and G4P << 67 - Temporarily made QSS the default method for << 68 << 69 ## 2023-01-30 Rodrigo Castro, Lucio Santi & Le << 70 - Included Quantum State Simulation (QSS) inte << 71 Alternative method of integration which crea << 72 functions for each coordinate (x,p) - changi << 73 polynomials when the deviation from exact so << 74 a threshold value. << 75 Offers built-in interpolation capability - e << 76 of the intersection of the trajectory with s << 77 Currently restricted to using only with a pu << 78 Created by Rodrigo Castro and Lucio Santi, w << 79 by Leandro A. Gómez Vidal (and others?) << 80 << 81 ## 2023-03-16 Ben Morgan (field-V11-01-02) << 82 - Fixed compilation warnings on clang-15 for m << 83 in input parameters for steppers. << 84 << 85 ## 2023-01-09 Ben Morgan (field-V11-01-01) << 86 - Use PROJECT_ scoped _SOURCE/BINARY_DIR varia << 87 scope version to ease use of Geant4 as a CMa << 88 - Fixes [GitHub PR #52](https://github.com/G << 89 << 90 ## 2022-12-09 Laurie Nevay (field- << 91 - Reduced printout for setting any valid value << 92 G4FieldManager. << 93 << 94 ## 2022-11-28 Gabriele Cosmo (field- << 95 - Fixed restore of stream precision in G4Field << 96 << 97 ## 2022-11-14 John Apostolakis (field- << 98 - Revised G4FieldManager to ensure that epsilo << 99 are less than a 'maximum accepted' accuracy << 100 behaviour of the integration. Improved their << 101 - warnings if min > max, with corrective be << 102 - a fatal exception in case of values outsi << 103 << 104 To cope with needs of legacy applications or << 105 the value of the 'ceiling' maximum accepted << 106 the new static method << 107 G4FieldManager::SetMaxAcceptedEpsilon( << 108 but must remain under or equal to a final ce << 109 fMaxFinalEpsilon=0.03 << 110 << 111 ## 2022-11-10 Gabriele Cosmo (field-V11-00-03) << 112 - Fixed compilation warnings for implicit type << 113 << 114 ## 2022-11-05 Divyansh Tiwari, John Apostolaki << 115 << 116 - Introduced G4BorisScheme and G4BorisDriver, << 117 integration method, created as part of GSoC << 118 << 119 ## 2022-10-05 Gabriele Cosmo (field-V11-00-01) << 120 - Fixed compilation warnings on Intel/icx comp << 121 but not used. << 122 << 123 ## 2021-12-10 Ben Morgan (field-V11-00-00) << 124 - Change to new Markdown History format << 125 << 126 --- << 127 << 128 # History entries prior to 11.0 << 129 << 130 March 30, 2021 B.Morgan << 131 ----------------------- << 132 - Migrate sources.cmake to modular build API << 133 << 134 November 26, 2020 J.Apostolakis << 135 ------------------------------- << 136 - Fix to enable G4TDormandPrince45 to be used << 137 << 138 This fix was needed because G4Interpolation << 139 copies of the stepper class. Due to this it << 140 with the type of the templated stepper. << 141 << 142 For a stepper to be used with 'G4Interpolati << 143 implement the method << 144 Equation_type GetSpecificEquation() << 145 Templated steppers must ensure that the retu << 146 is the specific equation type of the templat << 147 << 148 This method was added to G4DormandPrince745 << 149 << 150 October 12, 2020 J.Apostolakis << 151 ------------------------------ << 152 - Fixes & refinements - mostly in templated cl << 153 Compilation fixes, e.g. missing include <cas << 154 Refinements, e.g. int -> unsigned int in tem << 155 << 156 October 9, 2020 J.Apostolakis << 157 ----------------------------- << 158 - Configure G4ChordFinder to use templated G4T << 159 default stepper for magnetic fields (when on << 160 explicitly by the application.) << 161 << 162 October 9, 2020 J.Apostolakis & J. Xie << 163 -------------------------------------- << 164 - Adapted classes that avoid virtual calls for << 165 created by Josh Xie (Google Summer of Code 2 << 166 S. Wenzel & J.A. << 167 << 168 The revised design and implementations inclu << 169 * equation of motion templated on the field << 170 * steppers templated on the type of equatio << 171 variables. << 172 In addition the key methods are marked as 'i << 173 they judge possible, to embed them and avoid << 174 << 175 These classes can be combined also with the << 176 (G4IntegrationStepper, G4FSALIntegrationStep << 177 to avoid virtual calls in all the levels up << 178 the integration driver. ] << 179 << 180 The types of stepper currently available are << 181 i) the originally developed intermediate ord << 182 << 183 - Name (templated) - original -- Order/ << 184 / E << 185 =========================================== << 186 G4TClassicalRK4: G4ClassicalRK4 4th / << 187 G4TCashKarpRKF45: G4CashKarpRKF45 5th / << 188 << 189 ii) the original low order steppers (for spe << 190 G4TSimpleHeum: G4SimpleHeum 3rd / << 191 G4TSimpleRunge: G4SimpleRunge 2nd / << 192 G4TExplicitEuler: G4ExplicitEuler 1st / << 193 << 194 and << 195 iii) some newly adapted stepper(s) / ( J.A. << 196 G4TDormandPrince45 G4DormandPrince745 5th/ << 197 =========================================== << 198 << 199 Note: Additional lower-order steppers are no << 200 'Vector' steppers (using blaze) are no << 201 << 202 Experimental versions of fields with inline << 203 G4TUniformMagField, G4TQuadrupoleMagFiel << 204 are included to show how the field evaluatio << 205 << 206 Note 2: << 207 Steppers which are not embedded obtain their << 208 breaking the step into two halves, then comp << 209 As a result many more field / derivative eva << 210 embedded steppers of the same order, but the << 211 estimation has different properties; so it m << 212 << 213 June 15, 2020 G.Cosmo << 214 --------------------- << 215 - Implemented move constructor and operator fo << 216 << 217 March 24, 2020 J.Apostolakis << 218 ---------------------------- << 219 - New class G4DriverReporter to print progress << 220 Used for debugging differences between Old << 221 - Revisions in testing CMake config file and s << 222 for magnetic field. << 223 << 224 March 13, 2020 J.Apostolakis << 225 ---------------------------- << 226 - G4VIntegrationDriver & dependent driver clas << 227 Added new virtual StreamInfo() method, << 228 & used it to implement operator << for << 229 Made (virtual) DoesReIntegrate() const. << 230 - G4MagInt_Driver: Fixed max iterations & clar << 231 ComputeNewStepSize does NOT (yet) respect << 232 ( Done to enable comparisons with new G4Inte << 233 - G4OldMagIntDriver maintains all old behaviou << 234 << 235 << 236 January 22, 2020 G.Cosmo << 237 ------------------------ << 238 - Fixed compilation errors and configuration f << 239 << 240 January 13, 2020 G.Cosmo << 241 ------------------------ << 242 - Turn off verbosity flags by default in G4Int << 243 G4InterpolationDriver and G4MagIntegratorDri << 244 << 245 December 10, 2019 B.Morgan << 246 -------------------------- << 247 - Cleanup CMake build, removing obsolete granu << 248 explicit include_directories. << 249 2 >> 3 ========================================================= >> 4 Geant4 - an Object-Oriented Toolkit for Simulation in HEP >> 5 ========================================================= >> 6 >> 7 Category History file >> 8 --------------------- >> 9 This file should be used by G4 developers and category coordinators >> 10 to briefly summarize all major modifications introduced in the code >> 11 and keep track of all category-tags. >> 12 It DOES NOT substitute the CVS log-message one should put at every >> 13 committal in the CVS repository ! >> 14 >> 15 ---------------------------------------------------------- >> 16 * Reverse chronological order (last date on top), please * >> 17 ---------------------------------------------------------- 250 November 14, 2019 J.Apostolakis 18 November 14, 2019 J.Apostolakis - field-V10-05-17 251 ------------------------------- 19 ------------------------------- 252 - Added method / attribute to G4VIntegrationDr 20 - Added method / attribute to G4VIntegrationDriver DoesReIntegrate() 253 to inform whether the driver recomputes seg 21 to inform whether the driver recomputes segment when 254 AccurateAdvance is called. 22 AccurateAdvance is called. 255 23 256 October 27, 2019 D.Sorokin 24 October 27, 2019 D.Sorokin - field-V10-05-16 257 ---------------------------- 25 ---------------------------- 258 - Remove statistics printouts from destructor 26 - Remove statistics printouts from destructor of G4BFieldIntegrationDriver 259 27 260 October 24, 2019 H.Burkhardt 28 October 24, 2019 H.Burkhardt - field-V10-05-15 261 ---------------------------- 29 ---------------------------- 262 - Standard sextupole magnet first implementati 30 - Standard sextupole magnet first implementation G4SextupoleMagField.cc,.hh 263 31 264 October 18, 2019 D. Sorokin 32 October 18, 2019 D. Sorokin - field-V10-05-14 265 ----------------------------------------- 33 ----------------------------------------- 266 - Fix floating point exception in G4BFieldInte 34 - Fix floating point exception in G4BFieldIntegrationDriver. 267 35 268 October 1, 2019 G. Cosmo 36 October 1, 2019 G. Cosmo - field-V10-05-13 269 ----------------------------------------- 37 ----------------------------------------- 270 - Fixed unitialised data in G4EqMagElectricFie 38 - Fixed unitialised data in G4EqMagElectricField and G4EqGravityField, 271 reported by Coverity. 39 reported by Coverity. 272 40 273 September 12, 2019 G. Cosmo 41 September 12, 2019 G. Cosmo - field-V10-05-12 274 ----------------------------------------- 42 ----------------------------------------- 275 - Minor C++11 revision. Code formatting and cl 43 - Minor C++11 revision. Code formatting and cleanup. 276 44 277 July 8, 2019 D. Sorokin 45 July 8, 2019 D. Sorokin - field-V10-05-11 278 ----------------------------------------- 46 ----------------------------------------- 279 - G4BFieldIntegrationDriver: an attempt to sol << 47 - G4BFieldIntegrationDriver: an attempt to solve looping particle issue. 280 48 281 July 3, 2019 D. Sorokin 49 July 3, 2019 D. Sorokin - field-V10-05-10 282 ----------------------------------------- 50 ----------------------------------------- 283 - Make G4InterpolationDriver default for perfo 51 - Make G4InterpolationDriver default for performance assessments. 284 52 285 June 25, 2019 D. Sorokin 53 June 25, 2019 D. Sorokin - field-V10-05-09 286 ----------------------------------------- 54 ----------------------------------------- 287 - Fix Dormand-Prince 5th order interpolation m 55 - Fix Dormand-Prince 5th order interpolation method 288 56 289 June 20, 2019 D. Sorokin 57 June 20, 2019 D. Sorokin - field-V10-05-08 290 ----------------------------------------- 58 ----------------------------------------- 291 - New version of G4InterpolationDriver with mu 59 - New version of G4InterpolationDriver with multiple steppers for each substep. 292 60 293 June 12, 2019 D. Sorokin & J. Apostolakis 61 June 12, 2019 D. Sorokin & J. Apostolakis - field-V10-05-07 294 ----------------------------------------- 62 ----------------------------------------- 295 - Revert to use G4IntegrationDriver as default << 63 - Revert to use G4IntegrationDriver as default in G4ChordFinder 296 ( c-tor for magnetic field. ). 64 ( c-tor for magnetic field. ). 297 Change to use concrete type of Stepper as te 65 Change to use concrete type of Stepper as template parameter 298 66 299 June 10, 2019 Dmitry Sorokin 67 June 10, 2019 Dmitry Sorokin - field-V10-05-06 300 ------------------------- 68 ------------------------- 301 - Remove division in G4DormandPrince745::Inter 69 - Remove division in G4DormandPrince745::Interpolate4thOrder. 302 70 303 June 7, 2019 Dmitry Sorokin 71 June 7, 2019 Dmitry Sorokin - field-V10-05-05 304 ------------------------- 72 ------------------------- 305 - Switch to original chord finding algorithm i 73 - Switch to original chord finding algorithm in G4InterpolationDriver 306 - Reset step estimate in G4MagInt_Driver only 74 - Reset step estimate in G4MagInt_Driver only at the beginning of a track 307 - Change G4VERBOSE to G4FIELD_DEBUG in G4Inter 75 - Change G4VERBOSE to G4FIELD_DEBUG in G4InterpolationDriver to speed-up 308 76 309 June 4, 2019 Dmitry Sorokin 77 June 4, 2019 Dmitry Sorokin - field-V10-05-04 310 ------------------------- 78 ------------------------- 311 - Use FSAL property of G4DormandPrince745 in G 79 - Use FSAL property of G4DormandPrince745 in G4InterpolationDriver. 312 80 313 June 3, 2019 Dmitry Sorokin 81 June 3, 2019 Dmitry Sorokin - field-V10-05-03 314 ------------------------- 82 ------------------------- 315 - Fix memory churn in G4InterpolationDriver. 83 - Fix memory churn in G4InterpolationDriver. 316 84 317 May 20, 2019 Dmitry Sorokin 85 May 20, 2019 Dmitry Sorokin - field-V10-05-02 318 ------------------------- 86 ------------------------- 319 - Make interpolation driver default for immedi 87 - Make interpolation driver default for immediate testing. 320 88 321 April 8, 2019 H.Burkhardt 89 April 8, 2019 H.Burkhardt - field-V10-05-01 322 ------------------------- 90 ------------------------- 323 - Correcting transformation to local system in 91 - Correcting transformation to local system in G4QuadrupoleMagField.cc 324 92 325 March 20, 2019 G.Cosmo 93 March 20, 2019 G.Cosmo - field-V10-05-00 326 ---------------------- 94 ---------------------- 327 - Corrected 'RegularStepperName' in G4ChordFin 95 - Corrected 'RegularStepperName' in G4ChordFinder to refer to the actual 328 default stepper G4DormandPrince745 rather th 96 default stepper G4DormandPrince745 rather than G4NystromRK4... 329 97 330 November 9, 2018 J.Apostolakis 98 November 9, 2018 J.Apostolakis - field-V10-04-17 331 ------------------------------ 99 ------------------------------ 332 - Small fixes for warnings in compilation. 100 - Small fixes for warnings in compilation. 333 - Revisions of descriptions of recent tags (fo 101 - Revisions of descriptions of recent tags (for clarification.) 334 102 335 September 26, 2018 D.Sorokin 103 September 26, 2018 D.Sorokin - field-V10-04-15, -16 336 ---------------------------- 104 ---------------------------- 337 - Implemented G4InterpolationDriver. This driv 105 - Implemented G4InterpolationDriver. This driver uses a Runge-Kutta 338 stepper with interpolation property to integ 106 stepper with interpolation property to integrate EoM with error control. 339 The field evaluations used to calculate the 107 The field evaluations used to calculate the enpoint of an 'advance' can 340 be reused to obtain one or multiple internal 108 be reused to obtain one or multiple internal points on the curve (x,p). 341 So in steps in which an intersection is foun 109 So in steps in which an intersection is found, the number of field 342 evaluations are expected to be reduced subst 110 evaluations are expected to be reduced substantially. 343 111 344 September 13, 2018 D.Sorokin 112 September 13, 2018 D.Sorokin - field-V10-04-13, -14 345 ---------------------------- 113 ---------------------------- 346 - G4NystromRK4 creates G4CachedMagneticField i 114 - G4NystromRK4 creates G4CachedMagneticField if it was not provided. 347 - Removed G4MagIntegratorStepper::ComputeRight 115 - Removed G4MagIntegratorStepper::ComputeRightHandSide() method 348 (all calls now use non-virtual RightHandSide 116 (all calls now use non-virtual RightHandSide() method instead). 349 - Reversed declaration of final for the class 117 - Reversed declaration of final for the class G4UniformMagField. 350 Returned to previous status: only GetFieldVa 118 Returned to previous status: only GetFieldValue method is final. 351 ( This avoid need to change examples/tests a << 119 ( This avoid need to change examples/tests and potentially migrate 352 users' code. ) << 120 users' code. ) 353 121 354 September 12, 2018 D.Sorokin 122 September 12, 2018 D.Sorokin - field-V10-04-11, -12 355 ---------------------------- 123 ---------------------------- 356 - Created new experimental driver G4BFieldInte 124 - Created new experimental driver G4BFieldIntegrationDriver for pure 357 magnetic fields. In QuickAdvance it uses he 125 magnetic fields. In QuickAdvance it uses helix stepper when the 358 estimated 'helix' rotation angle is larger t 126 estimated 'helix' rotation angle is larger than pi/3. 359 - G4NystromRK4 uses G4CachedMagneticField for 127 - G4NystromRK4 uses G4CachedMagneticField for caching values. 360 NOTE: It now the user's responsibility to cr << 128 NOTE: It now the user's responsibility to create G4CachedMagneticField. 361 ( The alternative is to create G4CachedMagne 129 ( The alternative is to create G4CachedMagneticField inside G4NystromRK4. 362 But this would have several disadvantages: 130 But this would have several disadvantages: 363 * It would change the field in the user's 131 * It would change the field in the user's equation of motion; 364 * G4NystromRK4 is not able to check whethe 132 * G4NystromRK4 is not able to check whether the user calls the method 365 that changes the equation using SetEquat 133 that changes the equation using SetEquationOfMotion method. ) 366 134 367 July 10, 2018 D. Sorokin - field- 135 July 10, 2018 D. Sorokin - field-V10-04-08, -09, -10 368 ------------------------ Tags by J. Apostolaki 136 ------------------------ Tags by J. Apostolakis 369 - Corrections in calculation of chord distance 137 - Corrections in calculation of chord distance in G4DormandPrince745 370 ( use 4th order interpolation from Shampine 138 ( use 4th order interpolation from Shampine 1986 for mid-point. ) 371 - Revised G4FieldUtils: moved copy() method fr 139 - Revised G4FieldUtils: moved copy() method from steppers, 372 added 2 new setValue: simple & variadic meth 140 added 2 new setValue: simple & variadic methods. 373 - Fixes to copy input time in ModifiedMidpoint 141 - Fixes to copy input time in ModifiedMidpoint, CashKarp, G4RK547FEq1/2/3 374 Bogacki-Shampine 4/5 and 2/3 (note: B-S 2/3 142 Bogacki-Shampine 4/5 and 2/3 (note: B-S 2/3 largely rewritten ) 375 Dormand Prince 4/5, 5/6 and 7/8, and Tsitou 143 Dormand Prince 4/5, 5/6 and 7/8, and Tsitouras steppers. 376 - Trial revision in G4BulirschStoer to reduce 144 - Trial revision in G4BulirschStoer to reduce number of division 377 operations (J.A.) 145 operations (J.A.) 378 - Fix for unitialised pIntStepper pointer in G 146 - Fix for unitialised pIntStepper pointer in G4FSALIntegrationDriver 379 constructor. Revision of data member variab 147 constructor. Revision of data member variable names. ( J.A. ) 380 << 148 381 June 15, 2018 G.Cosmo - fiel 149 June 15, 2018 G.Cosmo - field-V10-04-07 382 --------------------- 150 --------------------- 383 - Fixed minor Coverity defects in G4BurlischSt 151 - Fixed minor Coverity defects in G4BurlischStoer and G4ModifiedMidPoint. 384 Code cleanup. 152 Code cleanup. 385 153 386 June 13, 2018 G.Cosmo - fiel 154 June 13, 2018 G.Cosmo - field-V10-04-06 387 --------------------- 155 --------------------- 388 - Leave override final only for clone() method 156 - Leave override final only for clone() method in G4UniformMagField, 389 not in destructor... 157 not in destructor... 390 158 391 June 12, 2018 G.Cosmo - fiel 159 June 12, 2018 G.Cosmo - field-V10-04-05 392 --------------------- 160 --------------------- 393 - Use G4Exception warnings instead of streamin 161 - Use G4Exception warnings instead of streaming to cerr. 394 - Added override final for virtual methods in 162 - Added override final for virtual methods in G4UniformMagField (J.Apostolakis). 395 - Some code cleanup. 163 - Some code cleanup. 396 164 397 Apr 30, 2018 J. Apostolakis - fiel 165 Apr 30, 2018 J. Apostolakis - field-V10-04-04 398 --------------------------- 166 --------------------------- 399 - G4VIntegrationDriver and derived classes 167 - G4VIntegrationDriver and derived classes 400 168 401 Added new virtual method in G4VIntegrationDr 169 Added new virtual method in G4VIntegrationDriver : 402 void RenewStepperAndAdjust(G4MagIntegrat 170 void RenewStepperAndAdjust(G4MagIntegratorStepper *pItsStepper) 403 which used to exist in G4MagInt_Driver. << 171 which used to exist in G4MagInt_Driver. 404 ( Its absence caused a problem in some exist 172 ( Its absence caused a problem in some existing user code since 10.4. ) 405 173 406 This method must be defined in those derived 174 This method must be defined in those derived driver classes which accept 407 a Runge Kutta stepper ( i.e. some or any typ 175 a Runge Kutta stepper ( i.e. some or any type of G4MagIntegratorStepper ) 408 176 409 - G4IntegrationDriver, G4FSALIntegrationDriver 177 - G4IntegrationDriver, G4FSALIntegrationDriver 410 Defined default implementation (which throws << 178 Defined default implementation (which throws a Fatal G4Exception), and 411 implementations in the templated classes G4I 179 implementations in the templated classes G4IntegrationDriver and 412 G4FSALIntegrationDriver which check the type 180 G4FSALIntegrationDriver which check the type (via dynamic cast) and 413 change the stepper if it is appropriate. 181 change the stepper if it is appropriate. 414 182 415 Feb 13, 2018 J. Apostolakis - fiel 183 Feb 13, 2018 J. Apostolakis - field-V10-04-03 416 --------------------------- 184 --------------------------- 417 - Change in G4FieldManager::CreateChordFinder << 185 - Change in G4FieldManager::CreateChordFinder to cope 418 better with case that field argument is null 186 better with case that field argument is null. 419 187 420 Feb 13, 2018 D. Sorokin - fiel 188 Feb 13, 2018 D. Sorokin - field-V10-04-01 421 ----------------------- - fiel 189 ----------------------- - field-V10-04-02 (comp. fix) 422 - Introduced Bulirsch-Stoer method of integrat 190 - Introduced Bulirsch-Stoer method of integration, an alternative 423 to Runge-Kutta based on the mid-point method 191 to Runge-Kutta based on the mid-point method. 424 Separate driver derived G4VIntegrationDriver 192 Separate driver derived G4VIntegrationDriver and implemented 425 using specialisation of G4IntegrationDriver 193 using specialisation of G4IntegrationDriver template class. 426 Note: this method does not have a RK 'Steppe 194 Note: this method does not have a RK 'Stepper' ! 427 195 428 Jan 31, 2018 J. Apostolakis 196 Jan 31, 2018 J. Apostolakis 429 --------------------------- 197 --------------------------- 430 - G4MagneticField no longer inherits from G4El 198 - G4MagneticField no longer inherits from G4ElectroMagneticField 431 This was incorrect, and allowed a user to pr << 199 This was incorrect, and allowed a user to provide a full EM 432 field, but obtain the force only from the B- 200 field, but obtain the force only from the B-field when using 433 the equation created in the G4ChordFinder co 201 the equation created in the G4ChordFinder constructor. 434 => potential interface change. 202 => potential interface change. 435 203 436 Result: It is no longer possible to use the << 204 Result: It is no longer possible to use the equation of motion for a 437 pure magnetic field for the case of a mix 205 pure magnetic field for the case of a mixed electromagnetic field. 438 This change is desired, and is a *protect 206 This change is desired, and is a *protection* as that inadvertent 439 use led to incorrect results. << 207 use led to incorrect results. 440 Instead the full chain of classes for equati 208 Instead the full chain of classes for equation of full EM field, 441 steppers and driver must be used for the 209 steppers and driver must be used for the case of the 'full' EM field. 442 210 443 Feb 9, 2017 G.Cosmo - fiel 211 Feb 9, 2017 G.Cosmo - field-V10-04-00 444 ------------------- 212 ------------------- 445 - Removed static data from G4FSALIntegrationDr 213 - Removed static data from G4FSALIntegrationDriver, G4IntegrationDriver and 446 G4MagIntegratorDriver and moved to base clas 214 G4MagIntegratorDriver and moved to base class as normal data members. 447 Fix required for allowing proper symbol expo 215 Fix required for allowing proper symbol exporting on Windows DLLs. 448 216 449 Nov 17, 2017 J. Apostolakis - fiel 217 Nov 17, 2017 J. Apostolakis - field-V10-03-30 450 --------------------------- 218 --------------------------- 451 - G4ChordFinder: Reinstated old driver G4MagIn << 219 - G4ChordFinder: Reinstated old driver G4MagInt_Drv, for use 452 with G4DormandPrince745 (5th order) in G4C 220 with G4DormandPrince745 (5th order) in G4ChordFinder's 453 2nd constructor (with magnetic field.) 221 2nd constructor (with magnetic field.) 454 ( Reverses change made in field-V10-03-27 ) 222 ( Reverses change made in field-V10-03-27 ) 455 223 456 Nov 16, 2017 J. Apostolakis - fiel 224 Nov 16, 2017 J. Apostolakis - field-V10-03-29 457 --------------------------- 225 --------------------------- 458 - Suppressed statistics printing in G4ChordFin 226 - Suppressed statistics printing in G4ChordFinder and new drivers 459 ( G4IntegrationDriver and G4FSALIntegration 227 ( G4IntegrationDriver and G4FSALIntegrationDriver ) 460 protecting them by verbosity flags. 228 protecting them by verbosity flags. 461 << 229 462 - Minor fix to initialise scalar data member i 230 - Minor fix to initialise scalar data member in new steppers 463 G4RK547FEq1, G4RK547FEq2 and G4RK547FEq3. << 231 G4RK547FEq1, G4RK547FEq2 and G4RK547FEq3. 464 Fixes coverity issues 97886, 97888 and 97889 232 Fixes coverity issues 97886, 97888 and 97889. 465 233 466 234 467 Nov 14, 2017 J. Apostolakis - fiel 235 Nov 14, 2017 J. Apostolakis - field-V10-03-28 468 --------------------------- 236 --------------------------- 469 - Suppressed extra verbosity in G4ChordFinder. 237 - Suppressed extra verbosity in G4ChordFinder. 470 - Fix in G4FSALDormandPrince745 for unitialise 238 - Fix in G4FSALDormandPrince745 for unitialised variables. 471 Initialise ak8[], ak9[] in constructor. This 239 Initialise ak8[], ak9[] in constructor. This avoids repeated 472 memory allocation, without deletion - i.e. a 240 memory allocation, without deletion - i.e. a memory leak. 473 - Initialise fLastStepLength and/or fAuxSteppe << 241 - Initialise fLastStepLength and/or fAuxStepper data members in 474 in constructors of G4DormandPrince745 and G4 242 in constructors of G4DormandPrince745 and G4FSALDormandPrince745. 475 Coverity issues: 67394 (G4FSALDormandPrince7 243 Coverity issues: 67394 (G4FSALDormandPrince745) 476 67391 (G4DormandPrince7 244 67391 (G4DormandPrince745) 477 245 478 Nov 14, 2017 J. Apostolakis - fiel 246 Nov 14, 2017 J. Apostolakis - field-V10-03-27 (fixed) 479 --------------------------- - fiel 247 --------------------------- - field-V10-03-26 480 - G4ChordFinder: changes in default stepper & 248 - G4ChordFinder: changes in default stepper & driver in 2nd constructor. 481 in which only a magnetic field is provided ( 249 in which only a magnetic field is provided (used widely) : 482 << 250 483 * Reinstated G4DormandPrince745 (5th order) 251 * Reinstated G4DormandPrince745 (5th order) as default in G4ChordFinder. 484 (corrected -23, which missed to initialise 252 (corrected -23, which missed to initialise driver.) 485 253 486 * New templated driver G4IntegrationDriver i << 254 * New templated driver G4IntegrationDriver is chosen for the 487 the simple constuctors of G4ChordFinder an 255 the simple constuctors of G4ChordFinder and G4FieldManager 488 which only require a (pure) magnetic field 256 which only require a (pure) magnetic field. ( Details below.) 489 257 490 Additional information and details: 258 Additional information and details: 491 * New optional parameter to constructor of G 259 * New optional parameter to constructor of G4ChordFinder which 492 allow to create an FSAL stepper and corres 260 allow to create an FSAL stepper and corresponding driver instead. 493 This aims to reuse some derivatives, thus 261 This aims to reuse some derivatives, thus reducing computations 494 of derivatives and thus field evaluations. 262 of derivatives and thus field evaluations. ) 495 263 496 * The templated 'driver' classes, available << 264 * The templated 'driver' classes, available and now chosen in 497 G4ChordFinder's 2nd constructor, are used 265 G4ChordFinder's 2nd constructor, are used per type of stepper 498 class: 266 class: 499 - G4IntegrationDriver for non-FSAL 267 - G4IntegrationDriver for non-FSAL steppers 500 - G4FSALIntegrationDriver for FSAL-capab 268 - G4FSALIntegrationDriver for FSAL-capable steppers 501 These make non-virtual calls to the steppe 269 These make non-virtual calls to the stepper, as they know its type. 502 270 503 - First fix in G4FSALDormandPrince745 for memo << 271 - First fix in G4FSALDormandPrince745 for memory leak, missing 504 delete[] fInitialDyDx in destructor. Coverit 272 delete[] fInitialDyDx in destructor. Coverity issue: 96383. 505 273 506 ---------------------------------------------- 274 ---------------------------------------------------------------------- 507 Multiple tags (below) were to better test step 275 Multiple tags (below) were to better test stepper/driver combinations: 508 Date Tagger Tag name Stepper 276 Date Tagger Tag name Stepper Driver 509 ============================================== 277 ====================================================================== 510 Nov 3, D. Sorokin field-V10-03-25 G4FSALBo << 278 Nov 3, D. Sorokin field-V10-03-25 G4FSALBogackiShampine45 FSAL 511 Nov 3, J. Apostolakis field-V10-03-24 G4FSALDo 279 Nov 3, J. Apostolakis field-V10-03-24 G4FSALDormandPrince745 FSAL ok 512 Nov 2, J. Apostolakis field-V10-03-23 G4Dorman << 280 Nov 2, J. Apostolakis field-V10-03-23 G4DormandPrince745 new 513 Nov 2, J. Apostolakis field-V10-03-22 G4RK547F 281 Nov 2, J. Apostolakis field-V10-03-22 G4RK547FEq1 FSAL ok 514 ---------------------------------------------- 282 ---------------------------------------------------------------------- 515 Glossary of drivers: 283 Glossary of drivers: 516 the FSAL driver = using G4FSALIntegrationDr 284 the FSAL driver = using G4FSALIntegrationDriver (templated) 517 the new driver = using G4IntegrationDriver 285 the new driver = using G4IntegrationDriver (templated) 518 ============================================== 286 =================================================================== 519 287 520 Nov 2, 2017 D. Sorokin - fiel 288 Nov 2, 2017 D. Sorokin - field-V10-03-21 521 --------------------------------------- 289 --------------------------------------- 522 - Added new steppers G4RK547FEq{1,2,3} which i << 290 - Added new steppers G4RK547FEq{1,2,3} which implement 7 stage embedded 523 Runge-Kutta 4,5 pairs with: 291 Runge-Kutta 4,5 pairs with: 524 - FSAL property (return derivative at end 292 - FSAL property (return derivative at endpoint for use in next step), 525 - stable equilibrium states ( improved ha 293 - stable equilibrium states ( improved handling of difficulties in 526 integration: better stability in choo 294 integration: better stability in choosing step size. ) 527 Each one can be used with all existing Runge 295 Each one can be used with all existing Runge-Kutta drivers: 528 - G4FSALIntegrationDriver, using its FSAL 296 - G4FSALIntegrationDriver, using its FSAL property; 529 - G4IntegrationDriver without its FSAL 297 - G4IntegrationDriver without its FSAL property. 530 - the original G4MagInt_Driver, also igno 298 - the original G4MagInt_Driver, also ignoring FSAL. 531 299 532 - Reinstated G4DormandPrince745 ( known widel 300 - Reinstated G4DormandPrince745 ( known widely as DOPRI5 ) as 533 the default RK method in G4ChordFinder for p 301 the default RK method in G4ChordFinder for pure magnetic fields. 534 ( Withdrew G4BogackiShampine45. ) J. Apos 302 ( Withdrew G4BogackiShampine45. ) J. Apostolakis 535 << 303 536 Nov 1, 2017 J. Apostolakis - fiel 304 Nov 1, 2017 J. Apostolakis - field-V10-03-20 537 --------------------------- 305 --------------------------- 538 - Made G4BogackiShampine45 the default RK meth 306 - Made G4BogackiShampine45 the default RK method in G4ChordFinder 539 for pure magnetic fields. 307 for pure magnetic fields. 540 308 541 Oct 30, 2017 D. Sorokin/J.Apostolakis - fiel 309 Oct 30, 2017 D. Sorokin/J.Apostolakis - field-V10-03-19 542 --------------------------------------- 310 --------------------------------------- 543 - add GetStepper method to G4VIntegrationDrive 311 - add GetStepper method to G4VIntegrationDriver for backward compatibility 544 312 545 Oct 20, 2017 D. Sorokin - fiel 313 Oct 20, 2017 D. Sorokin - field-V10-03-18 546 --------------------------------------- 314 --------------------------------------- 547 Introduced base class for any driver algorithm 315 Introduced base class for any driver algorithm (G4VIntegrationDriver) 548 New concrete classes added: 316 New concrete classes added: 549 - G4IntegrationDriver: implementation of templ 317 - G4IntegrationDriver: implementation of templated driver algorithm 550 for non-fsal explicit Runge-Kutta steppers 318 for non-fsal explicit Runge-Kutta steppers 551 - G4FSALIntegrationDriver: implementation of d 319 - G4FSALIntegrationDriver: implementation of driver algorithm for 552 FSAL explicit Runge-Kutta steppers. 320 FSAL explicit Runge-Kutta steppers. 553 Existing class G4MagInt_Drv remains, deriving 321 Existing class G4MagInt_Drv remains, deriving from base class. 554 322 555 - Expanded testing in testProElectroMagField ( 323 - Expanded testing in testProElectroMagField (higher precision, new driver) 556 ( J. Apostolakis) 324 ( J. Apostolakis) 557 325 558 Oct 13, 2017 D. Sorokin - fiel 326 Oct 13, 2017 D. Sorokin - field-V10-03-17 559 --------------------------------------- 327 --------------------------------------- 560 - add G4FieldUtils.hh, G4FieldUtils.icc, G4Fie << 328 - add G4FieldUtils.hh, G4FieldUtils.icc, G4FieldUtils.cc: 561 utilities to extract data from arrays in con << 329 utilities to extract data from arrays in conventions of the 562 magneticfield integration. 330 magneticfield integration. 563 331 564 Oct 13, 2017 G. Cosmo - fiel 332 Oct 13, 2017 G. Cosmo - field-V10-03-16 565 --------------------------------------- 333 --------------------------------------- 566 - G4NystromRK4: use consistenly G4 types for d 334 - G4NystromRK4: use consistenly G4 types for double and bool in G4NystromRK4. 567 335 568 Aug 31, 2017 J. Apostolakis - fiel 336 Aug 31, 2017 J. Apostolakis - field-V10-03-15 569 --------------------------------------- 337 --------------------------------------- 570 - Made G4DormandPrince457 the default Stepper 338 - Made G4DormandPrince457 the default Stepper for magnetic fields in 571 G4ChordFinder - when a driver/stepper is not 339 G4ChordFinder - when a driver/stepper is not specified. 572 340 573 Aug 31, 2017 J. Apostolakis - fiel 341 Aug 31, 2017 J. Apostolakis - field-V10-03-14 574 --------------------------------------- 342 --------------------------------------- 575 - G4NystromRK4: used caching of momentum magni 343 - G4NystromRK4: used caching of momentum magnitude ( & its inverse) 576 to reduce overhead of fix in field-V10-03-10 344 to reduce overhead of fix in field-V10-03-10. 577 345 578 Aug 25, 2017 J. Apostolakis - fiel 346 Aug 25, 2017 J. Apostolakis - field-V10-03-12, 13 579 --------------------------------------- 347 --------------------------------------- 580 - G4MagInt_Driver: added GetEquationOfMotion m 348 - G4MagInt_Driver: added GetEquationOfMotion method for enabling common 581 code with future develoment. 349 code with future develoment. 582 350 583 Aug 24, 2017 D.Sorokin & J. Apostolakis - fiel 351 Aug 24, 2017 D.Sorokin & J. Apostolakis - field-V10-03-11 584 --------------------------------------- 352 --------------------------------------- 585 - G4NystromRK4: correction in call to getField 353 - G4NystromRK4: correction in call to getField: needs p[4], not p[3] 586 ( Spotted by D. Sorokin ) 354 ( Spotted by D. Sorokin ) 587 Also small cleanup of testing code into 2 sm 355 Also small cleanup of testing code into 2 small new methods. 588 ( J. Apostolakis ) 356 ( J. Apostolakis ) 589 357 590 Aug 23, 2017 J.Apostolakis - fiel 358 Aug 23, 2017 J.Apostolakis - field-V10-03-10 591 -------------------------- 359 -------------------------- 592 - G4NystromRK4 - Fix the caching of momentum, 360 - G4NystromRK4 - Fix the caching of momentum, field location & value. 593 Fixes issue reported by ALICE ( problem repo 361 Fixes issue reported by ALICE ( problem report #1996 ). 594 Thanks to Ruben Shahoyan for identifying and 362 Thanks to Ruben Shahoyan for identifying and reporting the issue! 595 - Temporarly made G4NystromRK4 the default Ste 363 - Temporarly made G4NystromRK4 the default Stepper in G4ChordFinder's 2nd 596 constructor (to ensure adequate testing). 364 constructor (to ensure adequate testing). 597 365 598 Jul 04, 2017 J.Apostolakis - field 366 Jul 04, 2017 J.Apostolakis - field-V10-03-09 599 -------------------------- 367 -------------------------- 600 - G4DormandPrinceRK78: initialise fLastStepLen 368 - G4DormandPrinceRK78: initialise fLastStepLength and/or fAuxStepper data 601 members in constructor. Coverity issue: 6739 369 members in constructor. Coverity issue: 67398, 602 Improved description and reference to origin 370 Improved description and reference to original 1981 DP paper. 603 - G4MagIntegratorDriver: initialised fNoCalls 371 - G4MagIntegratorDriver: initialised fNoCalls in c-tor. Coverity #15158. 604 372 605 Jul 03, 2017 J.Apostolakis - field 373 Jul 03, 2017 J.Apostolakis - field-V10-03-08 606 -------------------------- 374 -------------------------- 607 - G4FSALBogackiShampline45: initialise 'bi[][] 375 - G4FSALBogackiShampline45: initialise 'bi[][]' const coefficients only 1 time 608 in method PrepareConstants() called by const 376 in method PrepareConstants() called by constructor. 609 ( Similar to G4BogackiShampine => seek con 377 ( Similar to G4BogackiShampine => seek consolidation ) 610 Small fix: initialise fAuxStepper to nullptr 378 Small fix: initialise fAuxStepper to nullptr, if not primary. 611 ( Part of Coverity issue 67394 ) 379 ( Part of Coverity issue 67394 ) 612 - G4DormandPrinceRK56: fix for significant mem 380 - G4DormandPrinceRK56: fix for significant memory leak of ak10_low buffer 613 space in SetupInterpolate_low(), apparently 381 space in SetupInterpolate_low(), apparently not seen by Coverity 614 because it was initialised in each use, not 382 because it was initialised in each use, not in constructor! 615 Improved initialisation of ak10, 11, 12 (mov 383 Improved initialisation of ak10, 11, 12 (moved to constructor) 616 and fixed small memory leak (added deletion 384 and fixed small memory leak (added deletion in destructor). 617 - G4MagIntegratorStepper: initialise fNoRHSCal 385 - G4MagIntegratorStepper: initialise fNoRHSCalls & fIntegrationOrder in c-tor. 618 ( fIntegrationOrder is not yet used widely - 386 ( fIntegrationOrder is not yet used widely - future migration ) 619 Coverity issue 67392. 387 Coverity issue 67392. 620 - G4VFSALIntegrationStepper: initialise counte 388 - G4VFSALIntegrationStepper: initialise counter in c-tor. Coverity issue 67396. 621 - G4BogackiShampine23: initialised fIsFSAL, fI 389 - G4BogackiShampine23: initialised fIsFSAL, fIntegrationOrder in c-tor. 622 Suppressed obsolete virtual method isFSAL, 390 Suppressed obsolete virtual method isFSAL, 623 Deleted duplicate virtual method 'isFSAL()'. 391 Deleted duplicate virtual method 'isFSAL()'. Non-virtual method IsFSAL() 624 now exists in base class G4MagIntegratorStep 392 now exists in base class G4MagIntegratorStepper. 625 Made copy constructor and equal operator pri 393 Made copy constructor and equal operator private and not unimplemented. 626 Added explanatory comments in source file; c 394 Added explanatory comments in source file; can make this method FSAL easily 627 (it calculates derivative @ endpoint. ) 395 (it calculates derivative @ endpoint. ) 628 - Several steppers ( G4DoLoMcPriRK34, G4FSALBo 396 - Several steppers ( G4DoLoMcPriRK34, G4FSALBogackiShampine45, 629 G4DormandPrince745, G4BogackiShampine45 ): i 397 G4DormandPrince745, G4BogackiShampine45 ): initialise fLastStepLength and/or 630 fAuxStepper data members in constructor. 398 fAuxStepper data members in constructor. 631 Coverity issues: 67393, 67394, 67391, 67390, 399 Coverity issues: 67393, 67394, 67391, 67390, 67363 632 400 633 Jul 03, 2017 J.Apostolakis - field << 401 Jul 03, 2017 J.Apostolakis - field-V10-03-05, -06, -07 634 -------------------------- 402 -------------------------- 635 - fix for memory leak in Stepper() in two step 403 - fix for memory leak in Stepper() in two steppers: 636 G4FSALBogackiShampline45 ( Coverity issue 404 G4FSALBogackiShampline45 ( Coverity issue 67389 ) 637 G4FSALDormandPrince745 ( Coverity issue 405 G4FSALDormandPrince745 ( Coverity issue 67397 ) 638 - G4BogackiShampline23: fix for memory leak in 406 - G4BogackiShampline23: fix for memory leak in destructor 639 ( Coverity issue 67355 ) 407 ( Coverity issue 67355 ) 640 408 641 May 17, 2017 J.Apostolakis - field 409 May 17, 2017 J.Apostolakis - field-V10-03-03, -04 (resubmit) 642 -------------------------- 410 -------------------------- 643 - G4FieldManager::SetDetectorField 411 - G4FieldManager::SetDetectorField 644 Fixed default of 'fail mode' parameter in << 412 Fixed default of 'fail mode' parameter in 645 to 0 (silent) - not 1 which was used for 413 to 0 (silent) - not 1 which was used for debugging. 646 In case of existing G4ChordFinder, use 'max' 414 In case of existing G4ChordFinder, use 'max' (1, failMode) not min. 647 415 648 416 649 May 11, 2017 J.Apostolakis - field 417 May 11, 2017 J.Apostolakis - field-V10-03-02 650 -------------------------- << 418 -------------------------- 651 - Reintroduced changes of field-V10-01-05, 06 419 - Reintroduced changes of field-V10-01-05, 06 and 07 to address 652 issue recalled by ATLAS. 420 issue recalled by ATLAS. 653 - New fixes/refinements (since field-V10-01-07 421 - New fixes/refinements (since field-V10-01-07) 654 o Extra check for missing equation. (Fix) 422 o Extra check for missing equation. (Fix) 655 o Renamed G4FieldManager::SetFieldChangesEne 423 o Renamed G4FieldManager::SetFieldChangesEnergy() method to 656 InitialiseFieldChangesEnergy() to clarify 424 InitialiseFieldChangesEnergy() to clarify its utility. 657 - Original functionality change (copied from f 425 - Original functionality change (copied from field-V10-01-05, & -07) : 658 o G4FieldManager.cc : Revised behaviour of S 426 o G4FieldManager.cc : Revised behaviour of SetDetector - it now 659 passes the Field to the Equation, as is th 427 passes the Field to the Equation, as is the natural user expectation 660 The previous behaviour is possible, using 428 The previous behaviour is possible, using extra argument 'failMode' flag. 661 Values: 0 = silent ; << 429 Values: 0 = silent ; 662 1 (default) = warning; << 430 1 (default) = warning; 663 2 or other = fatal except 431 2 or other = fatal exception. 664 432 665 This change is to ensure that user expecta << 433 This change is to ensure that user expectation is met: 666 if the field is changed for an existing, w << 434 if the field is changed for an existing, working stepper it must 667 be propagated to the equation of motion. 435 be propagated to the equation of motion. 668 Note: this causes potential conflict betwe 436 Note: this causes potential conflict between the type of the equation 669 and the field, which will need to be 437 and the field, which will need to be solved in the future. 670 Old behaviour (just set our own pointer) i 438 Old behaviour (just set our own pointer) is still available in new method: 671 ProposeDetectorField(G4Field *detectorFiel << 439 ProposeDetectorField(G4Field *detectorField) 672 Note: Error reporting behaviour of SetDete << 440 Note: Error reporting behaviour of SetDetector is 'soft'. It does not 673 complain by default if ChordFinder i 441 complain by default if ChordFinder is not (yet) set. 674 ( The default values of failMode = 0 442 ( The default values of failMode = 0 i.e. quiet. ) 675 If a chord-finder exists, it will still wa 443 If a chord-finder exists, it will still warn for other errors. 676 To obtain have noisy extra checks call Set 444 To obtain have noisy extra checks call SetDetectorField with failMode > 0 677 or by using the new method for an existing 445 or by using the new method for an existing field: 678 ChangeDetectorField(G4Field *detectorField << 446 ChangeDetectorField(G4Field *detectorField) 679 OLD: Co-works and requires new tag in geomet 447 OLD: Co-works and requires new tag in geometry/navigation: geomnav-V10-01-24. 680 448 681 Mar 22, 2017 G.Cosmo - fiel 449 Mar 22, 2017 G.Cosmo - field-V10-03-01 682 --------------------- 450 --------------------- 683 - Revert change made in previous tag and add s 451 - Revert change made in previous tag and add simple accessor to return 684 singleton pointer. 452 singleton pointer. 685 453 686 Mar 17, 2017 G.Cosmo - fiel 454 Mar 17, 2017 G.Cosmo - field-V10-03-00 687 --------------------- 455 --------------------- 688 - Define G4FieldManagerStore as normal singlet 456 - Define G4FieldManagerStore as normal singleton and corrected destructor. 689 457 690 Oct 7, 2016 J.Apostolakis - fiel 458 Oct 7, 2016 J.Apostolakis - field-V10-02-24, 25 691 --------------------------- 459 --------------------------- 692 - Checked loops for termination for infinite l << 460 - Checked loops for termination for infinite loops and annotated 693 them with "Loop checking" comment 461 them with "Loop checking" comment 694 462 695 Oct 7, 2016 D.Sorokin - fiel 463 Oct 7, 2016 D.Sorokin - field-V10-02-23 696 ----------------------- Commit: J. Apostolaki 464 ----------------------- Commit: J. Apostolakis 697 - Protection for multiple inclusion of G4Class << 465 - Protection for multiple inclusion of G4ClassicalRK4 698 (thanks to Dmitrii Sorokin, GSoC 2016 stude 466 (thanks to Dmitrii Sorokin, GSoC 2016 student. ) 699 467 700 Aug 16, 2016 J.Apostolakis - fiel 468 Aug 16, 2016 J.Apostolakis - field-V10-02-22 701 --------------------------- 469 --------------------------- 702 - Added initialization of fAuxStepper in const << 470 - Added initialization of fAuxStepper in constructor of 703 G4DoLoMcPriRK34, G4DormandPrince745 & G4Dorm 471 G4DoLoMcPriRK34, G4DormandPrince745 & G4DormandPrinceRK56. 704 472 705 June 11, 2016 J.Apostolakis - fiel 473 June 11, 2016 J.Apostolakis - field-V10-02-21 706 --------------------------- 474 --------------------------- 707 - Restored G4ClassicalRK4 as default stepper f 475 - Restored G4ClassicalRK4 as default stepper for B-field, 708 (due to performance issues reported in bench 476 (due to performance issues reported in benchmarking.) 709 << 477 710 June 10, 2016 G.Cosmo - fiel 478 June 10, 2016 G.Cosmo - field-V10-02-20 711 --------------------- 479 --------------------- 712 - Fixed code indentation in G4HelixMixedSteppe 480 - Fixed code indentation in G4HelixMixedStepper; 713 fixes compilation warnings on gcc-6.1. 481 fixes compilation warnings on gcc-6.1. 714 482 715 June 1, 2016 J.Apostolakis - fiel 483 June 1, 2016 J.Apostolakis - field-V10-02-19, 18, 17 716 -------------------------- 484 -------------------------- 717 - Fix in multiple steppers for the size of two 485 - Fix in multiple steppers for the size of two C-arrays which hold 718 intermediate values; they must be at least 486 intermediate values; they must be at least 8 and at least 719 'fNumberOfStateVariables' long. << 487 'fNumberOfStateVariables' long. 720 Changed: G4CashKarp, << 488 Changed: G4CashKarp, 721 G4BogackiShampine 23 & 34 489 G4BogackiShampine 23 & 34 722 G4DormandPrince 745, RK56 & RK78 490 G4DormandPrince 745, RK56 & RK78 723 G4DoLoMcPriRK34 491 G4DoLoMcPriRK34 724 G4FSALBogackiShampine45, G4FSALDorm 492 G4FSALBogackiShampine45, G4FSALDormandPrince745 725 - Made G4CashKarp the default stepper in G4Cho 493 - Made G4CashKarp the default stepper in G4ChordFinder. 726 - Deleted trial class G4NewChordFinder. 494 - Deleted trial class G4NewChordFinder. 727 495 728 June 1, 2016 J.Apostolakis - fiel 496 June 1, 2016 J.Apostolakis - field-V10-02-16 729 -------------------------- 497 -------------------------- 730 - Deleted test code (RkTestRun) due to cmake i 498 - Deleted test code (RkTestRun) due to cmake issues 731 499 732 June 1, 2016 J.Apostolakis - fiel 500 June 1, 2016 J.Apostolakis - field-V10-02-15 (also -14) 733 -------------------------- 501 -------------------------- 734 - Made G4ClassicalRK4 the default stepper agai 502 - Made G4ClassicalRK4 the default stepper again in G4ChordFinder 735 - Fixes in test code (RkTestRun) for new names 503 - Fixes in test code (RkTestRun) for new names of classes. 736 << 504 737 May 31, 2016 J.Apostolakis - fiel 505 May 31, 2016 J.Apostolakis - field-V10-02-13 (also -12) 738 -------------------------- 506 -------------------------- 739 - New class G4NewChordFinder, that uses old G4 507 - New class G4NewChordFinder, that uses old G4MagIntegratorDriver or 740 the newer G4FSALIntegratorDriver -- first 508 the newer G4FSALIntegratorDriver -- first try, for testing FSAL. 741 Created new test testNewChordFinder.cc to ch 509 Created new test testNewChordFinder.cc to checks this. 742 510 743 - Renamed new Stepper classes to use 'G4' pref 511 - Renamed new Stepper classes to use 'G4' prefix, including 744 DormandPrince{745,RK56,RK67} TsitourasRK45 512 DormandPrince{745,RK56,RK67} TsitourasRK45 and DoLoMcPriRK34 745 513 746 - Renamed FSAL base classes to use 'G4' prefix 514 - Renamed FSAL base classes to use 'G4' prefix & 'Integrat-ion' 747 FSALMagIntegratorStepper to G4VFSALIntegr 515 FSALMagIntegratorStepper to G4VFSALIntegrationStepper 748 FSALIntegratorDriver to G4FSALIntegr 516 FSALIntegratorDriver to G4FSALIntegrationDriver 749 and the experimental FSAL versions of steppe 517 and the experimental FSAL versions of steppers, 750 FBogackiShampine45.hh to G4FSALBogackiS 518 FBogackiShampine45.hh to G4FSALBogackiShampine45.hh 751 FDormandPrince745.hh to G4FSALDormandP 519 FDormandPrince745.hh to G4FSALDormandPrince745.hh 752 ( Note: FSAL steppers remain 'behind' the impr 520 ( Note: FSAL steppers remain 'behind' the improved non-FSAL versions.) 753 521 754 - Made DormandPrince745 (DP45) the default ste 522 - Made DormandPrince745 (DP45) the default stepper (in G4ChordFinder). 755 << 523 756 May 26, 2016 J.Apostolakis & Som. Banerjee - 524 May 26, 2016 J.Apostolakis & Som. Banerjee - field-V10-02-11 757 ------------------------------------------ 525 ------------------------------------------ 758 - First version of FSAL classes (as additional 526 - First version of FSAL classes (as additional classes) 759 FSAL Integrator Driver : New *concrete*, << 527 FSAL Integrator Driver : New *concrete*, stand-alone driver 760 FSAL Integrator Stepper : New base class ( 528 FSAL Integrator Stepper : New base class (separate from G4MagIntegratorStepper) 761 FBogackiShampine45: FSAL-version of Bogac 529 FBogackiShampine45: FSAL-version of BogackiShampine45 stepper 762 FDormandPrince745: FSAL-version of FDorm 530 FDormandPrince745: FSAL-version of FDormandPrince745 stepper 763 Developed by Somnath Banerjee in GSoC 2015 531 Developed by Somnath Banerjee in GSoC 2015 764 532 765 Context: 533 Context: 766 * FSAL Stepper's key 'Stepper()' method for St 534 * FSAL Stepper's key 'Stepper()' method for Stepping has additional 767 argument to return the derivative at the e 535 argument to return the derivative at the end of a step. 768 If the step was successful, this estimate 536 If the step was successful, this estimate has the same order of 769 accuracy as the stepper (i.e. 5th order fo 537 accuracy as the stepper (i.e. 5th order for a 5th order stepper.) 770 538 771 - FSAL Driver method uses the revised Stepper 539 - FSAL Driver method uses the revised Stepper interface, using 772 the final (endpoint) derivative for succes 540 the final (endpoint) derivative for successful steps 773 as the start derivative for the next step, 541 as the start derivative for the next step, in case of 774 continuing of integration. 542 continuing of integration. 775 543 776 - New test (RkTestRun.cc) covering FSAL method 544 - New test (RkTestRun.cc) covering FSAL methods, integration 777 Developed by Somnath Banerjee in GSoC 2015 545 Developed by Somnath Banerjee in GSoC 2015 778 546 779 May 26, 2016 J.Apostolakis - 547 May 26, 2016 J.Apostolakis - field-V10-02-10 780 -------------------------- 548 -------------------------- 781 Preparation for addition of FSAL steppers (& d 549 Preparation for addition of FSAL steppers (& driver) 782 - G4MagIntegratorStepper: Added counter for ca << 550 - G4MagIntegratorStepper: Added counter for calls to equation RHS, 783 with Get/Reset methods. 551 with Get/Reset methods. 784 - Cosmetic fixes in sources.cmake (tabs to spa 552 - Cosmetic fixes in sources.cmake (tabs to spaces) 785 553 786 May 25, 2016 J.Apostolakis & Som. Banerjee - 554 May 25, 2016 J.Apostolakis & Som. Banerjee - field-V10-02-09 787 ------------------------------------------ 555 ------------------------------------------ 788 - Added new stepper classes embedded RK method 556 - Added new stepper classes embedded RK method: 789 * DoLoMcPriRK34, 6-stage 3/4 RK (interp 557 * DoLoMcPriRK34, 6-stage 3/4 RK (interpolation) 790 * DormandPrinceRK56, 9-stage 5/6 RK (interp 558 * DormandPrinceRK56, 9-stage 5/6 RK (interpolation, FSAL-able) "6(5)" 791 * DormandPrinceRK78, 13-stage 7/8 RK (interp 559 * DormandPrinceRK78, 13-stage 7/8 RK (interpolation) "8(7) 13M" 792 created by Somnath Banerjee as part of GSoC 560 created by Somnath Banerjee as part of GSoC 2015. 793 561 794 May 25, 2016 J.Apostolakis - field 562 May 25, 2016 J.Apostolakis - field-V10-02-08 795 -------------------------- 563 -------------------------- 796 - Revised testPropagateMagField, to reduce err 564 - Revised testPropagateMagField, to reduce error severity for over-long 797 steps. An error message is still written 565 steps. An error message is still written if move > (1.+1.0e-9) * step, 798 but test does not stop. 566 but test does not stop. 799 567 800 - G4HelixMixedStepper: Fixes & added new (5th 568 - G4HelixMixedStepper: Fixes & added new (5th order) stepper choices. 801 Included (G4)Bogacki-Shampine 23/45 , Dor 569 Included (G4)Bogacki-Shampine 23/45 , DormandPrince745, 802 and newest TsitourasRK45. 570 and newest TsitourasRK45. 803 571 804 May 24, 2016 J.Apostolakis - field 572 May 24, 2016 J.Apostolakis - field-V10-02-06, 07 805 -------------------------- 573 -------------------------- 806 - Added TsitourasRK45 stepper. 574 - Added TsitourasRK45 stepper. 807 - Made Cash-Karp the default stepper (in G4Cho 575 - Made Cash-Karp the default stepper (in G4ChordFinder). 808 << 576 809 May 23, 2016 J.Apostolakis - field 577 May 23, 2016 J.Apostolakis - field-V10-02-05 810 -------------------------- 578 -------------------------- 811 - Made DormandPrince745 (DP45) the default ste 579 - Made DormandPrince745 (DP45) the default stepper (in G4ChordFinder). 812 << 580 813 May 20, 2016 J.Apostolakis - field 581 May 20, 2016 J.Apostolakis - field-V10-02-04 814 ------------------------- 582 ------------------------- 815 - Added steppers created by Somnath Banerjee a 583 - Added steppers created by Somnath Banerjee as part of GSoC 2015, 816 G4BogackiShampine23 (BS23 below), G4Bogack 584 G4BogackiShampine23 (BS23 below), G4BogackiShampine45 (BS45) and 817 DormandPrince745 (DP45). 585 DormandPrince745 (DP45). 818 These implement third order (BS23) and fifth 586 These implement third order (BS23) and fifth order (BS45, DP45) 819 embedded RK tableaus. 587 embedded RK tableaus. 820 - DP45 includes interpolation method, to evalu 588 - DP45 includes interpolation method, to evaluate the interpolant 821 ('y') at intermediate points of a successf 589 ('y') at intermediate points of a successful step. 822 BS45's interpolation method exists, but is n 590 BS45's interpolation method exists, but is not validated (fails). 823 591 824 References: 592 References: 825 - BS45 is embedded pair from paper by P. Bog 593 - BS45 is embedded pair from paper by P. Bogacki and L. F. Shampine: 826 “An efficient Runge-Kutta (4,5) pair, 594 “An efficient Runge-Kutta (4,5) pair,” 827 Comput. Math. with Appl., vol. 32, no. 6 595 Comput. Math. with Appl., vol. 32, no. 6, pp. 15–28, Sep. 1996. 828 - BS23 is four stage RK method from the pape 596 - BS23 is four stage RK method from the paper of P. Bogacki and L. F. Shampine, 829 “A 3(2) pair of Runge - Kutta formulas << 597 “A 3(2) pair of Runge - Kutta formulas,” 830 Appl. Math. Lett., vol 2, no 4, p.321 (1 598 Appl. Math. Lett., vol 2, no 4, p.321 (1989). 831 - DP45 are from Table 2 of J. R. Dormand an << 599 - DP45 are from Table 2 of J. R. Dormand and P. J. Prince's 832 “A family of embedded Runge-Kutta form 600 “A family of embedded Runge-Kutta formulae,” 833 Journal of comput. and applied Math., vo << 601 Journal of comput. and applied Math., vol. 6, no. 1, pp. 19–26, 1980. 834 602 835 May 4, 2016 G.Cosmo - field- 603 May 4, 2016 G.Cosmo - field-V10-02-03 836 ------------------- 604 ------------------- 837 - Fixed Clone() signatures in derived Field cl 605 - Fixed Clone() signatures in derived Field classes. 838 More corrections on implementations of opera 606 More corrections on implementations of operator=() and copy-ctors. 839 607 840 Apr 29, 2016 G.Cosmo - field- 608 Apr 29, 2016 G.Cosmo - field-V10-02-02 841 -------------------- 609 -------------------- 842 - Corrected implementation of operator=() and 610 - Corrected implementation of operator=() and copy-ctor for G4Field, 843 G4ElectroMagneticField, G4ElectricField and 611 G4ElectroMagneticField, G4ElectricField and G4CachedMagneticField. 844 612 845 Feb 25, 2016 P.Gumplinger - field- 613 Feb 25, 2016 P.Gumplinger - field-V10-02-01 846 ------------------------- 614 ------------------------- 847 - address Problem Report: 1831 615 - address Problem Report: 1831 848 change: G4Mag_SpinEqRhs, G4EqEMFieldWithSpin 616 change: G4Mag_SpinEqRhs, G4EqEMFieldWithSpin, G4EqEMFieldWithEDM and 849 G4RepleteEofM.cc from g_BMT = (magMo 617 G4RepleteEofM.cc from g_BMT = (magMoment/muB)/spin 850 to g_BMT = (std:: 618 to g_BMT = (std::abs(magMoment)/muB)/spin 851 619 852 620 853 Dec 11, 2015 A.Dotti - field-V10-02-0 621 Dec 11, 2015 A.Dotti - field-V10-02-00 854 -------------------------- 622 -------------------------- 855 Fixing typo in tests/CMakeList.txt file: subdi << 623 Fixing typo in tests/CMakeList.txt file: subdirecotry() 856 instead of add_subdirectory() 624 instead of add_subdirectory() 857 625 858 Nov 2, 2015 J.Apostolakis - field 626 Nov 2, 2015 J.Apostolakis - field-V10-01-15 (fixed) 859 Oct 22, 2015 J.Apostolakis - field 627 Oct 22, 2015 J.Apostolakis - field-V10-01-14 860 -------------------------- 628 -------------------------- 861 - Reverted changes in G4FieldManager from tag 629 - Reverted changes in G4FieldManager from tags field-V10-01-05 to -07 862 inadvertently reintroduced in tag field-V10 630 inadvertently reintroduced in tag field-V10-01-13 863 (fixed in tag -15: committed new copy from 631 (fixed in tag -15: committed new copy from tag -04 'on top') 864 << 632 865 - Added explanation comment in G4MagHelicalSte 633 - Added explanation comment in G4MagHelicalSteper. 866 << 634 867 Oct 13, 2015 J.Apostolakis - field 635 Oct 13, 2015 J.Apostolakis - field-V10-01-13 868 -------------------------- 636 -------------------------- 869 - G4ChordFinder::FindNextChord: Added check ag << 637 - G4ChordFinder::FindNextChord: Added check against number of iterations, 870 to identify potential issue of lack of conve 638 to identify potential issue of lack of convergence. 871 - test: fixed testProElectricMagField.cc 639 - test: fixed testProElectricMagField.cc 872 640 873 - Note: tag-13 inadvertently reintroduced chan 641 - Note: tag-13 inadvertently reintroduced changes from tags field-V10-01-05 to -07 874 642 875 Jun 26, 2015 P.Gumplinger - field 643 Jun 26, 2015 P.Gumplinger - field-V10-01-12, 11 876 ------------------------- 644 ------------------------- 877 - Include 'NumberOfVariables' in G4RepleteEofM 645 - Include 'NumberOfVariables' in G4RepleteEofM constructor argument list. 878 646 879 Jun 22, 2015 G.Cosmo - field 647 Jun 22, 2015 G.Cosmo - field-V10-01-10 880 -------------------- 648 -------------------- 881 - Withdrawn also changes introduced in tag "fi 649 - Withdrawn also changes introduced in tag "field-V10-01-08", as 882 responsible for memory violation errors in V 650 responsible for memory violation errors in Valgrind. 883 Now corresponding to tag "field-V10-01-04". 651 Now corresponding to tag "field-V10-01-04". 884 652 885 Jun 18, 2015 G.Cosmo - field 653 Jun 18, 2015 G.Cosmo - field-V10-01-09 886 -------------------- 654 -------------------- 887 - Withdrawn changes introduced from tag "field 655 - Withdrawn changes introduced from tag "field-V10-01-05" to tag 888 "field-V10-01-07" included, responsible for 656 "field-V10-01-07" included, responsible for failures on SLC6 and 889 run-time errors on Valgrind checks. 657 run-time errors on Valgrind checks. 890 658 891 Jun 10, 2015 P.Gumplinger - field 659 Jun 10, 2015 P.Gumplinger - field-V10-01-08 892 ------------------------- 660 ------------------------- 893 - initialize all elements of dydx[] in G4Reple 661 - initialize all elements of dydx[] in G4RepleteEofM::EvaluateRhsGivenB 894 Thanks to Tatsuya Kikawa (Osaka Univ.) 662 Thanks to Tatsuya Kikawa (Osaka Univ.) 895 663 896 Jun 10, 2015 J.Apostolakis - field 664 Jun 10, 2015 J.Apostolakis - field-V10-01-07 897 -------------------------- 665 -------------------------- 898 - G4FieldManager.cc : Made behaviour of SetDet 666 - G4FieldManager.cc : Made behaviour of SetDetector 'softer'. No 899 longer complains by default if ChordFinder 667 longer complains by default if ChordFinder is not (yet) set. 900 ( The default values of failMode = 0 i.e. 668 ( The default values of failMode = 0 i.e. quiet. ) 901 If a chord-finder exists, it will still war 669 If a chord-finder exists, it will still warn for other errors. 902 << 670 903 It is still possible to have louder/extra c 671 It is still possible to have louder/extra checks done calling 904 SetDetectorField with failMode > 0 672 SetDetectorField with failMode > 0 905 or by using a new method, which expects an 673 or by using a new method, which expects an existing field: 906 ChangeDetectorField(G4Field *detectorFiel << 674 ChangeDetectorField(G4Field *detectorField) 907 675 908 Jun 8, 2015 J.Apostolakis - field 676 Jun 8, 2015 J.Apostolakis - field-V10-01-05, (-06 fixed) 909 -------------------------- 677 -------------------------- 910 - G4FieldManager.cc : Revised behaviour of Set 678 - G4FieldManager.cc : Revised behaviour of SetDetector - it now 911 passes the Field to the Equation, as is the 679 passes the Field to the Equation, as is the natural user expectation 912 The previous behaviour is possible, using e << 680 The previous behaviour is possible, using extra argument 'failMode' 913 flag. Values: 0 = silent ; << 681 flag. Values: 0 = silent ; 914 1 (default) = warning; << 682 1 (default) = warning; 915 2 or other = fatal except 683 2 or other = fatal exception. 916 Or can use alternative (new) suggestion met 684 Or can use alternative (new) suggestion method ProposeDetectorField(); 917 685 918 This change is to ensure that user expectat << 686 This change is to ensure that user expectation is met: 919 if the field is changed for an existing, << 687 if the field is changed for an existing, working stepper it must 920 be propagated to the equation of motion. 688 be propagated to the equation of motion. 921 Note: this causes potential conflict betwee 689 Note: this causes potential conflict between the type of the equation 922 and the field, which will need to be 690 and the field, which will need to be solved in the future. 923 691 924 Old behaviour (just set our own pointer) is 692 Old behaviour (just set our own pointer) is still available in new method: 925 ProposeDetectorField(G4Field *detectorFie << 693 ProposeDetectorField(G4Field *detectorField) 926 694 927 Co-works and requires new tag in geometry/n 695 Co-works and requires new tag in geometry/navigation: geomnav-V10-01-24. 928 696 929 Jun 4, 2015 P. Gumplinger - field 697 Jun 4, 2015 P. Gumplinger - field-V10-01-04 930 ------------------------- 698 ------------------------- 931 - change assignment to field instead of Field 699 - change assignment to field instead of Field array in G4RepleteEofM.cc 932 700 933 Jan 28, 2015 J.Apostolakis - field 701 Jan 28, 2015 J.Apostolakis - field-V10-01-01 (-02,03 fixed) 934 -------------------------- 702 -------------------------- ( History ) 935 - Revert unintended inclusion of changes in G4 703 - Revert unintended inclusion of changes in G4MagIntegratorDriver 936 704 937 Jan 15, 2015 J.Apostolakis - field << 705 Jan 15, 2015 J.Apostolakis - field-V10-01-00 938 -------------------------- 706 -------------------------- 939 - Small refinements in G4FieldTrack - access r 707 - Small refinements in G4FieldTrack - access rest mass, avoid div by 0 940 Avoid division by zero for Unit direction << 708 Avoid division by zero for Unit direction in FieldTrack 941 Added method to get Rest Mass. 709 Added method to get Rest Mass. 942 Needed for fixes related to Problem report 1 710 Needed for fixes related to Problem report 1696. 943 Thanks to Nathan S. Froemming <nfroemm@uw.ed 711 Thanks to Nathan S. Froemming <nfroemm@uw.edu> who reported the issue 944 and created dedicated test to demonstrate i 712 and created dedicated test to demonstrate it. 945 713 946 Nov 03, 2014 G.Cosmo - field 714 Nov 03, 2014 G.Cosmo - field-V10-00-04 947 -------------------- 715 -------------------- 948 - Moved constructors and simple methods to in 716 - Moved constructors and simple methods to in line in G4ChargeState and 949 G4LineSection for CPU speedup (courtesy of S 717 G4LineSection for CPU speedup (courtesy of S.Wenzel). 950 718 951 May 21, 2014 J.Apostolakis - field 719 May 21, 2014 J.Apostolakis - field-V10-00-03 952 -------------------------- 720 -------------------------- 953 - G4FieldTrack: 721 - G4FieldTrack: 954 * Fixed SetChargeAndMoments() to use the new 722 * Fixed SetChargeAndMoments() to use the new method 955 G4ChargeState::SetChargesAndMoments(); it 723 G4ChargeState::SetChargesAndMoments(); it was changing spin, because of 956 erroneous change in G4ChargeState::SetChar 724 erroneous change in G4ChargeState::SetChargesAndMoments(), that method 957 is now renamed SetChargeSpinMoments() to s 725 is now renamed SetChargeSpinMoments() to show what it does; 958 * Added pdg-spin as argument to constructors 726 * Added pdg-spin as argument to constructors (needed to track polarization); 959 * Renamed 'spin' to 'Polarization' in method 727 * Renamed 'spin' to 'Polarization' in method, arguments, data member to 960 avoid confusion; 728 avoid confusion; 961 * Revised constructors to add pdg-spin as la 729 * Revised constructors to add pdg-spin as last argument (with default -1). 962 - G4ChargeState: 730 - G4ChargeState: 963 - Addressed Coverity issue #56281 'Arguments << 731 - Addressed Coverity issue #56281 'Arguments in wrong order' 964 G4ChargeState::SetChargeMoments() when use 732 G4ChargeState::SetChargeMoments() when used in 965 G4FieldTrack::SetChargeMoments() from addi 733 G4FieldTrack::SetChargeMoments() from addition of spin to arguments, 966 without revising name of method to reflect 734 without revising name of method to reflect it; 967 - Changed order of arguments in constructor, 735 - Changed order of arguments in constructor, making spin 3rd argument 968 (was 2nd) and making all arguments mandato 736 (was 2nd) and making all arguments mandatory; 969 - Renamed SetChargeAndMoments() method with 737 - Renamed SetChargeAndMoments() method with 5 arguments including spin (2nd) 970 to SetChargeSpinMoments(); 738 to SetChargeSpinMoments(); 971 - Added new methods including SetChargesAndM 739 - Added new methods including SetChargesAndMoments(), SetChargeMdm(), 972 SetChargeDipoleMoments() (no spin), and al 740 SetChargeDipoleMoments() (no spin), and also SetChargeSpin(), 973 SetChargeMdmSpin() (which includes spin) t << 741 SetChargeMdmSpin() (which includes spin) to provide most/all potential 974 combinations. 742 combinations. 975 - G4MagIntegratorDriver::OneGoodStep(): 743 - G4MagIntegratorDriver::OneGoodStep(): 976 * Trial changes to address Valgrind issue, a 744 * Trial changes to address Valgrind issue, about value used before being set. 977 745 978 Apr 01, 2014 J.Apostolakis - fi 746 Apr 01, 2014 J.Apostolakis - field-V10-00-02 979 -------------------------- 747 -------------------------- 980 - Moved G4FieldTrack::LoadFromArray() to .cc a 748 - Moved G4FieldTrack::LoadFromArray() to .cc and removed non-necessary 981 thread-local static data local to the functi 749 thread-local static data local to the function. 982 750 983 Mar 14, 2014 J.Apostolakis - fi 751 Mar 14, 2014 J.Apostolakis - field-V10-00-01 984 -------------------------- 752 -------------------------- 985 - Fix assignment operator for G4FieldTrack. (C 753 - Fix assignment operator for G4FieldTrack. (Copied Rest Mass) 986 - Improve printing of G4FieldTrack 754 - Improve printing of G4FieldTrack 987 755 988 Jan 17, 2014 P.Gumplinger - fi 756 Jan 17, 2014 P.Gumplinger - field-V10-00-00 989 ------------------------- 757 ------------------------- 990 - add G4Mag_EqRhs::SetChargeMomentumMass( part 758 - add G4Mag_EqRhs::SetChargeMomentumMass( particleCharge, MomentumXc, mass) 991 to G4Mag_SpinEqRhs::SetChargeMomentumMass, s 759 to G4Mag_SpinEqRhs::SetChargeMomentumMass, see bug reported at: 992 http://hypernews.slac.stanford.edu/HyperNews 760 http://hypernews.slac.stanford.edu/HyperNews/geant4/get/emfields/210.html 993 761 994 Nov 18, 2013 P.Gumplinger - fi 762 Nov 18, 2013 P.Gumplinger - field-V09-06-13/12 995 ------------------------- 763 ------------------------- 996 - fix Coverity warning in G4RepleteEofM 764 - fix Coverity warning in G4RepleteEofM 997 765 998 Nov 17, 2013 J.Apostolakis << 766 Nov 17, 2013 J.Apostolakis 999 -------------------------- 767 -------------------------- 1000 - Fixed and improved HelixMixed Stepper: 768 - Fixed and improved HelixMixed Stepper: 1001 Corrected the default stepper for small s 769 Corrected the default stepper for small steps to be ClassicalRK4. 1002 Added Angle-threshold parameter. (Default 770 Added Angle-threshold parameter. (Default remains 0.33 radians.) 1003 Added more Steppers, including G4NystromR 771 Added more Steppers, including G4NystromRK4. 1004 Revised Id numbers of Steppers, to make t 772 Revised Id numbers of Steppers, to make the best more visible. 1005 More comments in Stepper(). 773 More comments in Stepper(). 1006 774 1007 Nov 05, 2013 P.Gumplinger - f 775 Nov 05, 2013 P.Gumplinger - field-V09-06-11 1008 ------------------------- 776 ------------------------- 1009 - remove writing out-of-bounds dydx[8] = 0. i 777 - remove writing out-of-bounds dydx[8] = 0. in G4RepleteEofM.cc 1010 778 1011 Nov 04, 2013 P.Gumplinger - f 779 Nov 04, 2013 P.Gumplinger - field-V09-06-10 1012 ------------------------- 780 ------------------------- 1013 - add G4RepleteEofM to sources.cmake 781 - add G4RepleteEofM to sources.cmake 1014 f 782 field-V09-06-09 - same as 1015 Nov 01st, 2013 P.Gumplinger - f 783 Nov 01st, 2013 P.Gumplinger - field-V09-06-08 1016 --------------------------- 784 --------------------------- 1017 - add G4RepleteEofM: EofM that includes force 785 - add G4RepleteEofM: EofM that includes force due to B, E, and gravity field 1018 as well as mu*grad(B) and spin tracking. 786 as well as mu*grad(B) and spin tracking. 1019 787 1020 May 10th, 2013 P.Gumplinger/J.Apostolakis - f 788 May 10th, 2013 P.Gumplinger/J.Apostolakis - field-V09-06-07 1021 ----------------------------------------- 789 ----------------------------------------- 1022 - Revised number of field components (to allo 790 - Revised number of field components (to allow up to 24). 1023 18 are required for new applications (B, E 791 18 are required for new applications (B, E, Gravity, B-gradients) 1024 - Fixes for initialisation errors found by Co 792 - Fixes for initialisation errors found by Coverity. 1025 793 1026 May 10th, 2013 P.Gumplinger - field-V09-06-06 794 May 10th, 2013 P.Gumplinger - field-V09-06-06 1027 --------------------------- 795 --------------------------- 1028 - prevent divided by zero and shortcut spin t 796 - prevent divided by zero and shortcut spin tracking for spin == 0 1029 particles: G4Mag_SpinEqRhs.cc, G4EqEMFieldW 797 particles: G4Mag_SpinEqRhs.cc, G4EqEMFieldWithSpin.cc and 1030 G4EqEMFieldWithEDM.cc 798 G4EqEMFieldWithEDM.cc 1031 f << 799 field-V09-06-05 - same 1032 May 8th, 2013 P.Gumplinger/J.Apostolakis - f 800 May 8th, 2013 P.Gumplinger/J.Apostolakis - field-V09-06-04 1033 ----------------------------------------- 801 ----------------------------------------- 1034 - Introduced new class G4ChargeState to hold 802 - Introduced new class G4ChargeState to hold charge, spin, magnetic moment 1035 (promoted from internal class of G4Field 803 (promoted from internal class of G4FieldTrack). 1036 - Changed the signature of EquationOfMotion:: 804 - Changed the signature of EquationOfMotion::SetChargeMomentumMass to 1037 take G4ChargeState in place of G4double (cha 805 take G4ChargeState in place of G4double (charge) as first argument. 1038 - Deleted SetChargeMomentumMass from G4ChordF 806 - Deleted SetChargeMomentumMass from G4ChordFinder, G4PropagatorInField, 1039 G4MagInt_Driver 807 G4MagInt_Driver 1040 808 1041 - Migrated testPropagateMagField and testProp 809 - Migrated testPropagateMagField and testPropagateSpin. 1042 810 1043 March 19th, 2013 G.Cosmo - field-V09-06-03, f 811 March 19th, 2013 G.Cosmo - field-V09-06-03, field-V09-06-02 1044 ------------------------- 812 ------------------------- 1045 - Cleanup of MT changes in G4FieldManagerStor 813 - Cleanup of MT changes in G4FieldManagerStore. 1046 814 1047 February 26th, 2013 A.Dotti - field-V09-06-01 815 February 26th, 2013 A.Dotti - field-V09-06-01 1048 --------------------------- 816 --------------------------- 1049 - First implementation of Cloning mechanism ( 817 - First implementation of Cloning mechanism (needed for MT) 1050 818 1051 February 22nd, 2013 G.Cosmo - field-V09-06-00 819 February 22nd, 2013 G.Cosmo - field-V09-06-00 1052 --------------------------- 820 --------------------------- 1053 - Cleanup of unnecessary changes for multi-th 821 - Cleanup of unnecessary changes for multi-threading in G4QuadrupoleMagField. 1054 - Removed unused private data-members in G4Ma 822 - Removed unused private data-members in G4Mag_UsualEqRhs and G4RKG3_Stepper. 1055 823 1056 November 14th, 2012 P.Mato - field-V09-05-07 824 November 14th, 2012 P.Mato - field-V09-05-07 1057 -------------------------- 825 -------------------------- 1058 - G4QuadrupoleMagField: allow displaced origi 826 - G4QuadrupoleMagField: allow displaced origin and rotation (Björn Riese) 1059 - Added G4SystemOfUnits.hh in unit tests 827 - Added G4SystemOfUnits.hh in unit tests 1060 828 1061 August 30th, 2012 P.Mato - field-V09-05-06 829 August 30th, 2012 P.Mato - field-V09-05-06 1062 ------------------------ 830 ------------------------ 1063 - Adapted CMakeLists.txt files to new unit te 831 - Adapted CMakeLists.txt files to new unit tests 1064 832 1065 June 11th, 2012 G.Cosmo - field-V09-05-05 833 June 11th, 2012 G.Cosmo - field-V09-05-05 1066 ----------------------- 834 ----------------------- 1067 - Explicitly use inclusion of headers for sys 835 - Explicitly use inclusion of headers for system of units and physical 1068 constants, in plan to remove implicit inclu 836 constants, in plan to remove implicit inclusion from globals.hh. 1069 837 1070 May 14th, 2012 G.Cosmo - field-V09-05-04 838 May 14th, 2012 G.Cosmo - field-V09-05-04 1071 ---------------------- 839 ---------------------- 1072 - Fixed spurious case of hidden variable visi 840 - Fixed spurious case of hidden variable visibility in G4NystromRK4, 1073 detected with '-Wshadow' compilation option 841 detected with '-Wshadow' compilation option on gcc compiler. 1074 842 1075 February 3rd, 2012 J.Apostolakis - field-V09 843 February 3rd, 2012 J.Apostolakis - field-V09-05-03 1076 --------------------------------- 844 --------------------------------- 1077 - G4FieldTrack: << 845 - G4FieldTrack: 1078 Added new (inline) method to set the Rest << 846 Added new (inline) method to set the Rest Mass. Needed in order to 1079 reset this - and keep reuse the same G 847 reset this - and keep reuse the same G4FieldTrack in a client 1080 (across tracks.) 848 (across tracks.) 1081 Made SetCurvePnt() method private. 849 Made SetCurvePnt() method private. 1082 Simplified implementation of constructors << 850 Simplified implementation of constructors. Gain of 1 square root and 1083 1 division per call. 851 1 division per call. 1084 852 1085 February 3rd, 2012 G.Cosmo - field-V09-05-0 853 February 3rd, 2012 G.Cosmo - field-V09-05-02 1086 ---------------------------- 854 ---------------------------- 1087 - Removed tests in CMakeLists.txt causing cat 855 - Removed tests in CMakeLists.txt causing catastrophic failure in system 1088 testing! 856 testing! 1089 857 1090 February 2nd, 2012 J.Apostolakis - field-V0 858 February 2nd, 2012 J.Apostolakis - field-V09-05-01 1091 ---------------------------------- 859 ---------------------------------- 1092 - G4Mag_UsualEqRhs: in SetChargeMomentumMass( 860 - G4Mag_UsualEqRhs: in SetChargeMomentumMass() method, fixed problem with 1093 division by zero; stopped updating 'fInvCur 861 division by zero; stopped updating 'fInvCurrentMomentumXc' (data member 1094 is obsolete). 862 is obsolete). 1095 - Added tests in CMakeLists.txt 863 - Added tests in CMakeLists.txt 1096 864 1097 December 13th, 2011 G.Cosmo - field-V09-05-0 865 December 13th, 2011 G.Cosmo - field-V09-05-00 1098 ---------------------------- 866 ---------------------------- 1099 - Fixed Coverity defects in G4SimpleRunge (re 867 - Fixed Coverity defects in G4SimpleRunge (removed spurious pointer 1100 data member) and G4FieldTrack (missing G4Ch 868 data member) and G4FieldTrack (missing G4ChargeState::operator=()). 1101 - Synchronised unit tests to trunk (new test 869 - Synchronised unit tests to trunk (new test field07 and minor updates). 1102 870 1103 November 11th, 2011 J.Apostolakis - field-V0 871 November 11th, 2011 J.Apostolakis - field-V09-04-06 1104 ---------------------------------- 872 ---------------------------------- 1105 - Added G4Field.cc to list of files in source 873 - Added G4Field.cc to list of files in sources.cmake 1106 874 1107 November 4th, 2011 J.Apostolakis - field-V0 875 November 4th, 2011 J.Apostolakis - field-V09-04-05 1108 ---------------------------------- 876 ---------------------------------- 1109 - G4Field: revised to improve ability to dete 877 - G4Field: revised to improve ability to detect gravity field: 1110 Added new data member to identify gravity f 878 Added new data member to identify gravity field; 1111 Added new method IsGravityActive() to acces 879 Added new method IsGravityActive() to access this data member; 1112 Created new source file: moving constructor 880 Created new source file: moving constructors to it. 1113 - Changed G4UniformGravityField and G4Electro << 881 - Changed G4UniformGravityField and G4ElectroMagneticField 1114 to accomodate this refinement. 882 to accomodate this refinement. 1115 883 1116 November 2nd, 2011 P.Gumplinger - field-V09- 884 November 2nd, 2011 P.Gumplinger - field-V09-04-04 1117 -------------------------------- 885 -------------------------------- 1118 - update list of files in sources.cmake 886 - update list of files in sources.cmake 1119 887 1120 October 28th, 2011 P.Gumplinger - field-V09- 888 October 28th, 2011 P.Gumplinger - field-V09-04-03 1121 -------------------------------- 889 -------------------------------- 1122 - Added (uniform) gravity field capability. << 890 - Added (uniform) gravity field capability. 1123 New classes G4EqGravityField, G4UniformGrav 891 New classes G4EqGravityField, G4UniformGravityField. 1124 892 1125 *Required* for co-working tag in Transporta 893 *Required* for co-working tag in Transportation 1126 894 1127 October 12th, 2011 P.Gumplinger - field-V09- 895 October 12th, 2011 P.Gumplinger - field-V09-04-02 1128 -------------------------------- 896 -------------------------------- 1129 - replaced: omegac = 0.105658387*GeV/mass * 2 897 - replaced: omegac = 0.105658387*GeV/mass * 2.837374841e-3*(rad/cm/kilogauss) 1130 with: omegac = (eplus/particleMass)*c_l 898 with: omegac = (eplus/particleMass)*c_light 1131 in G4EqEMFieldWithEDM.cc, G4Mag_SpinEqRhs.c 899 in G4EqEMFieldWithEDM.cc, G4Mag_SpinEqRhs.cc, G4EqEMFieldWithSpin.cc 1132 changed calculation: << 900 changed calculation: 1133 std::sqrt(kineticEnergy*kineticEnergy + 901 std::sqrt(kineticEnergy*kineticEnergy + 2.0*restMass_c2*kineticEnergy); 1134 to: 902 to: 1135 std::sqrt(kineticEnergy*(kineticEnergy+2 903 std::sqrt(kineticEnergy*(kineticEnergy+2.0*restMass_c2)) 1136 in G4FieldTrack.cc 904 in G4FieldTrack.cc 1137 905 1138 August 11th, 2011 G.Cosmo - field-V09-04-01 906 August 11th, 2011 G.Cosmo - field-V09-04-01 1139 -------------------------- 907 -------------------------- 1140 - Migrated to new G4Exception scheme for warn 908 - Migrated to new G4Exception scheme for warnings and errors. 1141 - Defined scheme for warnings/errors coding: 909 - Defined scheme for warnings/errors coding: 1142 + 1xxx for warnings, 0xxx for errors. 910 + 1xxx for warnings, 0xxx for errors. 1143 Code numbering: 911 Code numbering: 1144 + 0001: not supported or implemented featur 912 + 0001: not supported or implemented features; 1145 + 0002: invalid setup/construct; 913 + 0002: invalid setup/construct; 1146 + 0003: invalid run condition; 914 + 0003: invalid run condition; 1147 + 1001: incorrect setup; automatically adap 915 + 1001: incorrect setup; automatically adapted. 1148 + 1002: notification. 916 + 1002: notification. 1149 917 1150 May 13th, 2011 J.Apostolakis 918 May 13th, 2011 J.Apostolakis 1151 ---------------------------- 919 ---------------------------- 1152 - Added new class G4TrialsCounter for keeping 920 - Added new class G4TrialsCounter for keeping statistics of trials. 1153 - Avoid printing out pointer in G4CachedMagne 921 - Avoid printing out pointer in G4CachedMagneticField::ReportStatistics(). 1154 - Updated unit tests. 922 - Updated unit tests. 1155 923 1156 May 12th, 2011 G.Cosmo - field-V09-04-00 924 May 12th, 2011 G.Cosmo - field-V09-04-00 1157 ----------------------- 925 ----------------------- 1158 - Fixed compilation warnings from gcc-4.6.0 f 926 - Fixed compilation warnings from gcc-4.6.0 for variables set but not used. 1159 927 1160 May 5th, 2011 J.Apostolakis 928 May 5th, 2011 J.Apostolakis 1161 ---------------------------- 929 ---------------------------- 1162 - Reordered data members to group invariants, 930 - Reordered data members to group invariants, field parameters, etc... 1163 931 1164 Sep 10th, 2010 J.Apostolakis - field-V09-03- 932 Sep 10th, 2010 J.Apostolakis - field-V09-03-03 1165 ----------------------------- 933 ----------------------------- 1166 - Revised constructor of G4MagErrorStepper to 934 - Revised constructor of G4MagErrorStepper to add number of State variables 1167 - Corrected MagErrorStepper to copy State Var 935 - Corrected MagErrorStepper to copy State Variable to output 1168 - Enable G4ConstRK4 to copy remaining State V 936 - Enable G4ConstRK4 to copy remaining State Variables (must integrate 6) 1169 937 1170 Sep 10th, 2010 J.Apostolakis - field-V09-03- 938 Sep 10th, 2010 J.Apostolakis - field-V09-03-02 1171 ----------------------------- 939 ----------------------------- 1172 - New class G4MonopoleEq for rhs of d(pc)/ds 940 - New class G4MonopoleEq for rhs of d(pc)/ds propagation in field 1173 in SI units (V.Grichine). 941 in SI units (V.Grichine). 1174 - Fixed passing of time in G4NystromRK4 942 - Fixed passing of time in G4NystromRK4 1175 943 1176 Jul 21st, 2010 T.Nikitina - field-V09-03-01 944 Jul 21st, 2010 T.Nikitina - field-V09-03-01 1177 -------------------------- 945 -------------------------- 1178 - Fixed cases of memory corruption in G4RKG3_ 946 - Fixed cases of memory corruption in G4RKG3_Stepper. 1179 - Fixed case of unused array data member in G 947 - Fixed case of unused array data member in G4ExactHelixStepper. 1180 - Removed useless code never executed in G4Co 948 - Removed useless code never executed in G4ConstRK4::Stepper(). 1181 - Fixed initialization in G4NystromRK4 constr 949 - Fixed initialization in G4NystromRK4 constructor. 1182 950 1183 Jul 14th, 2010 G.Cosmo - field-V09-03-00 951 Jul 14th, 2010 G.Cosmo - field-V09-03-00 1184 ----------------------- 952 ----------------------- 1185 - Added dummy initialization of members in co 953 - Added dummy initialization of members in constructors for G4CashKarpRKF45, 1186 G4ConstRK4, G4EqEMFieldWithEDM, G4EqEMField 954 G4ConstRK4, G4EqEMFieldWithEDM, G4EqEMFieldWithSpin, G4ExactHelixStepper, 1187 G4FieldTrack, G4MagHelicalStepper, G4MagInt 955 G4FieldTrack, G4MagHelicalStepper, G4MagInt_Driver, G4Mag_EqRhs, 1188 G4Mag_SpinEqRhs, G4Mag_UsualEqRhs, G4Nystro 956 G4Mag_SpinEqRhs, G4Mag_UsualEqRhs, G4NystromRK4, G4RKG3_Stepper, 1189 G4UniformElectricField, G4UniformElectricFi 957 G4UniformElectricField, G4UniformElectricField. 1190 958 1191 Nov 12th, 2009 J.Apostolakis - field-V09-02- 959 Nov 12th, 2009 J.Apostolakis - field-V09-02-09 1192 ----------------------------- 960 ----------------------------- 1193 - G4MagIntegratorDriver: activate check on i 961 - G4MagIntegratorDriver: activate check on integration error for spin. 1194 962 1195 Nov 12th, 2009 J.Apostolakis - field-V09-02- 963 Nov 12th, 2009 J.Apostolakis - field-V09-02-08 1196 ----------------------------- (fix only) 964 ----------------------------- (fix only) 1197 - G4Nystrom: Corrected interface method getF 965 - G4Nystrom: Corrected interface method getField: array now has explicit dimension[4] 1198 (Problem found by gcc 4.3 - i 966 (Problem found by gcc 4.3 - it checked indices used in inline method! ) 1199 967 1200 Nov 6th, 2009 P.Gumplinger - field-V09-02-07 968 Nov 6th, 2009 P.Gumplinger - field-V09-02-07 1201 --------------------------- 969 --------------------------- 1202 - bug fix in G4EqEMFieldWithSpin and G4EqEMFi 970 - bug fix in G4EqEMFieldWithSpin and G4EqEMFieldWithEDM 1203 thanks to Hiromi Iinuma (KEK) see: 971 thanks to Hiromi Iinuma (KEK) see: 1204 http://hypernews.slac.stanford.edu/HyperNew 972 http://hypernews.slac.stanford.edu/HyperNews/geant4/get/emfields/161.html 1205 973 1206 Nov 5th, 2009 J.Apostolakis - field-V09-02-0 974 Nov 5th, 2009 J.Apostolakis - field-V09-02-06 1207 ---------------------------- 975 ---------------------------- 1208 - G4MagIntegratorDriver.cc : Enabled call to 976 - G4MagIntegratorDriver.cc : Enabled call to ComputeRightHandSide 1209 - G4NystromRK4.cc : Disabled auxilia 977 - G4NystromRK4.cc : Disabled auxiliary code in Stepper (needed if 1210 ComputeRightHan 978 ComputeRightHandSide is not called.) 1211 979 1212 Nov 5th, 2009 J.Apostolakis - field-V09-02-0 980 Nov 5th, 2009 J.Apostolakis - field-V09-02-05 1213 ---------------------------- 981 ---------------------------- 1214 - Added new virtual method CalculateRightHand 982 - Added new virtual method CalculateRightHandSide to G4MagIntegratorStepper for use 1215 in caching momentum (and field value) b 983 in caching momentum (and field value) by G4NystromRK4 1216 Default implementation in G4MagIntegrat 984 Default implementation in G4MagIntegratorStepper calls RightHandSide inline method. 1217 - G4Nystrom: New Set/Get method for cache di 985 - G4Nystrom: New Set/Get method for cache distance. 1218 Changed private data members in 986 Changed private data members in G4NystromRK4. 1219 - G4MagIntegratorDriver: alternative call to 987 - G4MagIntegratorDriver: alternative call to ComputeRightHandSide is not used (in comment) 1220 As a result G4NystromRK4 operat << 988 As a result G4NystromRK4 operates without reusing 1221 989 1222 Nov 5th, 2009 J.Apostolakis - field-V09-02-0 990 Nov 5th, 2009 J.Apostolakis - field-V09-02-04 1223 ---------------------------- 991 ---------------------------- 1224 - G4CachedMagneticField: New Simple class to << 992 - G4CachedMagneticField: New Simple class to cache value of Magnetic Field. 1225 Speeds up code when calculation of field 993 Speeds up code when calculation of field value is complex. 1226 - G4NystromRK4 : New Stepper with Nys 994 - G4NystromRK4 : New Stepper with Nystrom for magnetic field 1227 with analytic estimation of integration << 995 with analytic estimation of integration error. 1228 Greatly reduces number of field value pe 996 Greatly reduces number of field value per step. 1229 - Revised testPropagateMagField to use Cached 997 - Revised testPropagateMagField to use Cached Quadrupole field, 1230 and to cover G4NystromRK4 and G4ConstRK4 998 and to cover G4NystromRK4 and G4ConstRK4. 1231 999 1232 Nov 4th, 2009 P.Gumplinger - field-V09-02-03 1000 Nov 4th, 2009 P.Gumplinger - field-V09-02-03 1233 --------------------------- 1001 --------------------------- 1234 - (minor change) remove comment from G4EqEMFi 1002 - (minor change) remove comment from G4EqEMFieldWithSpin.cc 1235 add G4EqEMFieldWithEDM class: this is the R 1003 add G4EqEMFieldWithEDM class: this is the RHS of EofM in a combined 1236 electric and magnetic field, with spin trac 1004 electric and magnetic field, with spin tracking for both MDM and EDM terms. 1237 Thanks to Kevin Lynch, Phys. Dept. at Bosto 1005 Thanks to Kevin Lynch, Phys. Dept. at Boston University. 1238 1006 1239 May 18th, 2009 T.Nikitina - field-V09-02-02 1007 May 18th, 2009 T.Nikitina - field-V09-02-02 1240 ------------------------- 1008 ------------------------- 1241 - Enhanced algorithm G4ChordFinder::ApproxCur 1009 - Enhanced algorithm G4ChordFinder::ApproxCurvePointS() in order to speedup 1242 BrentLocator. 1010 BrentLocator. 1243 1011 1244 March 25th, 2009 G.Cosmo - field-V09-02-01 1012 March 25th, 2009 G.Cosmo - field-V09-02-01 1245 ------------------------ 1013 ------------------------ 1246 - Some code cleanup and formatting... 1014 - Some code cleanup and formatting... 1247 1015 1248 March 6th, 2009 P.Gumplinger - field-V09-02-0 1016 March 6th, 2009 P.Gumplinger - field-V09-02-00 1249 ---------------------------- 1017 ---------------------------- 1250 - Added 3rd term of BMT equation (Spin x Beta 1018 - Added 3rd term of BMT equation (Spin x Beta x Efield) to G4EqEMFieldWithSpin, 1251 addresses emfields forum posting #155 (bug 1019 addresses emfields forum posting #155 (bug report). Thanks to Kevin Lynch, 1252 Phys. Dept. at Boston University. 1020 Phys. Dept. at Boston University. 1253 - Moved renormalization of spin from G4EqEMFi 1021 - Moved renormalization of spin from G4EqEMFieldWithSpin and G4Mag_SpinEqRhs to 1254 G4ClassicalRK4 and G4SimpleHeum. 1022 G4ClassicalRK4 and G4SimpleHeum. 1255 - Added Spin propagation errors to the criter 1023 - Added Spin propagation errors to the criteria for 'OneGoodStep' in 1256 G4MagIntegratorDriver but actually don't ad 1024 G4MagIntegratorDriver but actually don't add it (yet) to the decision logic. 1257 1025 1258 November, 19th, 2008 P.Gumplinger - field-V0 1026 November, 19th, 2008 P.Gumplinger - field-V09-01-05 1259 ---------------------------------- 1027 ---------------------------------- 1260 - Renormalized the spin to 1 in G4EqEMFieldWi 1028 - Renormalized the spin to 1 in G4EqEMFieldWithSpin and G4Mag_SpinEqRhs. 1261 1029 1262 November, 7th, 2008 P.Gumplinger - field-V0 1030 November, 7th, 2008 P.Gumplinger - field-V09-01-04 1263 ---------------------------------- 1031 ---------------------------------- 1264 - Added SetAnomaly() and GetAnomaly() accesso 1032 - Added SetAnomaly() and GetAnomaly() accessors to G4EqEMFieldWithSpin. 1265 1033 1266 October, 28th, 2008 T.Nikitina - field-V09- 1034 October, 28th, 2008 T.Nikitina - field-V09-01-03 1267 -------------------------------- 1035 -------------------------------- 1268 - Added 'ApproxCurveV' to method G4ChordFinde << 1036 - Added 'ApproxCurveV' to method G4ChordFinder::ApproxCurvePointS() in 1269 order to better calculate the curve length. 1037 order to better calculate the curve length. 1270 - Some code cleanup. 1038 - Some code cleanup. 1271 1039 1272 September, 18th, 2008 T.Nikitina 1040 September, 18th, 2008 T.Nikitina 1273 ---------------------------------- 1041 ---------------------------------- 1274 - Added new stepper class G4ConstRK4, perform 1042 - Added new stepper class G4ConstRK4, performing integration of one step 1275 with error calculation in constant magnetic 1043 with error calculation in constant magnetic field; implementation derived 1276 from G4ClassicalRK4. 1044 from G4ClassicalRK4. 1277 1045 1278 July, 15th, 2008 J.Apostolakis - field-V09- 1046 July, 15th, 2008 J.Apostolakis - field-V09-01-02 1279 -------------------------------- 1047 -------------------------------- 1280 - Revised signature of G4ChordFinder::FinderN 1048 - Revised signature of G4ChordFinder::FinderNextChord(), making 'FieldTrack' 1281 argument passed a reference. 1049 argument passed a reference. 1282 1050 1283 May, 28th, 2008 T.Nikitina 1051 May, 28th, 2008 T.Nikitina 1284 ---------------------------- 1052 ---------------------------- 1285 - G4ChordFinder: added new method ApproxCurve 1053 - G4ChordFinder: added new method ApproxCurvePointS() using Brent second order 1286 location mechanism. 1054 location mechanism. 1287 1055 1288 April 24th, 2008 T.Nikitina - field-V09-01- 1056 April 24th, 2008 T.Nikitina - field-V09-01-01 1289 ----------------------------- 1057 ----------------------------- 1290 - Fix 'unused' non-initialized values of dydx 1058 - Fix 'unused' non-initialized values of dydx in G4EqMagElectricField 1291 and G4EqEMFieldWithSpin. 1059 and G4EqEMFieldWithSpin. 1292 1060 1293 January 17th, 2008 G.Cosmo - field-V09-01-0 1061 January 17th, 2008 G.Cosmo - field-V09-01-00 1294 ---------------------------- 1062 ---------------------------- 1295 - G4CashKarpRKF45: deleted un-needed duplicat 1063 - G4CashKarpRKF45: deleted un-needed duplicate (and unused) instance variable 1296 'fNumberOfVariables', already defined in G4M 1064 'fNumberOfVariables', already defined in G4MagIntegratorStepper. (J.Ap.) 1297 - Some code cleanup. 1065 - Some code cleanup. 1298 1066 1299 December 7th, 2007 J.Apostolakis & T.Nikiti 1067 December 7th, 2007 J.Apostolakis & T.Nikitina - field-V09-00-03 1300 --------------------------------------------- 1068 ----------------------------------------------- 1301 - New class G4FieldManagerStore. It enables r 1069 - New class G4FieldManagerStore. It enables resetting of all chord 1302 finders' state. Needed to address problem r 1070 finders' state. Needed to address problem reported by user. 1303 1071 1304 September 9th, 2007 G.Cosmo - field-V09-00-0 1072 September 9th, 2007 G.Cosmo - field-V09-00-02a 1305 ---------------------------- 1073 ---------------------------- 1306 - Some code cleanup in G4MagHelicalStepper.cc 1074 - Some code cleanup in G4MagHelicalStepper.cc. 1307 1075 1308 August 30th, 2007 P.Gumplinger - field-V09 1076 August 30th, 2007 P.Gumplinger - field-V09-00-02 1309 --------------------------------- 1077 --------------------------------- 1310 - Add G4EqEMFieldWithSpin class to allow trac 1078 - Add G4EqEMFieldWithSpin class to allow tracking of spin also in 1311 G4ElectroMagneticFields 1079 G4ElectroMagneticFields 1312 1080 1313 August 21st, 2007 T.Nikitina - field-V09-0 1081 August 21st, 2007 T.Nikitina - field-V09-00-01 1314 ------------------------------- 1082 ------------------------------- 1315 - Fix to DistChord() method for Helical Stepp 1083 - Fix to DistChord() method for Helical Steppers: adopt 'special' case 1316 also for Ang<2pi not only for Ang<pi; use c 1084 also for Ang<2pi not only for Ang<pi; use cos() in place of tan(), to 1317 improve CPU performance. 1085 improve CPU performance. 1318 - Improved AdvanceHelix() to be able to make 1086 - Improved AdvanceHelix() to be able to make half step and full step 1319 in one call. 1087 in one call. 1320 - In G4HelixMixedStepper added possibility to 1088 - In G4HelixMixedStepper added possibility to call different steppers for 1321 small steps and add some statistics how off 1089 small steps and add some statistics how offen different steppers are called. 1322 - G4RKG3_Stepper: improved internal algorithm 1090 - G4RKG3_Stepper: improved internal algorithm to use pre-computed values; 1323 the stepper now gives comparable CPU perfor << 1091 the stepper now gives comparable CPU performance as ClassicalRK4. 1324 - Updated unit tests. 1092 - Updated unit tests. 1325 1093 1326 August 17th, 2007 G.Cosmo - field-V09-00-0 1094 August 17th, 2007 G.Cosmo - field-V09-00-00 1327 ---------------------------- 1095 ---------------------------- 1328 - Changed FatalException to EventMustBeAborte 1096 - Changed FatalException to EventMustBeAborted in AccurateAdvance() 1329 if proposed integration step is negative. 1097 if proposed integration step is negative. 1330 1098 1331 June 8th, 2007 T.Nikitina - field-V08-03-0 1099 June 8th, 2007 T.Nikitina - field-V08-03-03 1332 ---------------------------- 1100 ---------------------------- 1333 - Relaxed FatalException to Warning in G4MagI << 1101 - Relaxed FatalException to Warning in G4MagIntegratorDriver::AccurateAdvance() 1334 for case of step zero and return previous v << 1102 for case of step zero and return previous value. This fixes problem 1335 observed by ATLAS in release 8.3 when requi 1103 observed by ATLAS in release 8.3 when requiring very accurate propagation 1336 -- ie using very small values for acceptabl 1104 -- ie using very small values for acceptable integration error. 1337 1105 1338 May 18th, 2007 T.Nikitina - field-V08-03-0 1106 May 18th, 2007 T.Nikitina - field-V08-03-02 1339 ---------------------------- 1107 ---------------------------- 1340 - New class G4HelixMixedStepper, mixing helic 1108 - New class G4HelixMixedStepper, mixing helical stepper and RK 1341 G4MagHelicalStepper 1109 G4MagHelicalStepper 1342 - Keep values of inverse Radius, .. 1110 - Keep values of inverse Radius, .. 1343 - New calculation of DistChord for SteppingAn 1111 - New calculation of DistChord for SteppingAngle > pi 1344 G4RKG3Stepper 1112 G4RKG3Stepper 1345 - New calculation of DistChord for SteppingAn 1113 - New calculation of DistChord for SteppingAngle > pi 1346 1114 1347 May 16th, 2007 G.Cosmo - field-V08-03-01 1115 May 16th, 2007 G.Cosmo - field-V08-03-01 1348 ------------------------- 1116 ------------------------- 1349 - Introduced new class G4ErrorMag_UsualEqRhs 1117 - Introduced new class G4ErrorMag_UsualEqRhs used to reverse the magnetic 1350 field when propagation is backwards for err 1118 field when propagation is backwards for error propagation (by P.Arce). 1351 - Coworks with tag "global-V08-03-02". 1119 - Coworks with tag "global-V08-03-02". 1352 1120 1353 May 11th, 2007 G.Cosmo - field-V08-03-00 1121 May 11th, 2007 G.Cosmo - field-V08-03-00 1354 ------------------------- 1122 ------------------------- 1355 - Use call to G4GeometryTolerance instead of 1123 - Use call to G4GeometryTolerance instead of kCarTolerance constant in 1356 G4MagIntegratorDriver. 1124 G4MagIntegratorDriver. 1357 - Requires tag "global-V08-03-00" and related 1125 - Requires tag "global-V08-03-00" and related tag set. 1358 1126 1359 May 10th, 2007 J.Apostolakis - field-V08-0 1127 May 10th, 2007 J.Apostolakis - field-V08-02-01 1360 ------------------------------- 1128 ------------------------------- 1361 - G4MagIntegratorDriver: << 1129 - G4MagIntegratorDriver: 1362 * added protections against zero initial 1130 * added protections against zero initial or internal step h=0, (exceptions) 1363 * revised (re)calculation of h, and quit << 1131 * revised (re)calculation of h, and quit loop if h/length < smallest-fraction (~ 10^-12) 1364 1132 1365 April 26th, 2007 T.Nikitina - field-V08-02-0 1133 April 26th, 2007 T.Nikitina - field-V08-02-00 1366 ---------------------------- 1134 ---------------------------- 1367 - Corrected integration constants in G4MagHel 1135 - Corrected integration constants in G4MagHelicalStepper 1368 - Added implementation of G4RKG3_Stepper, ste 1136 - Added implementation of G4RKG3_Stepper, stepper derived from Geant3. 1369 - Updated NTST unit test. New unit test "fiel 1137 - Updated NTST unit test. New unit test "field06" for comparison of steppers. 1370 1138 1371 November 14th, 2006 J.Apostolakis - field-V0 1139 November 14th, 2006 J.Apostolakis - field-V08-01-02 1372 ---------------------------------- 1140 ---------------------------------- 1373 - G4MagIntegratorDriver 1141 - G4MagIntegratorDriver 1374 * Made first argument of GetDerivatives a r 1142 * Made first argument of GetDerivatives a reference - avoid copy 1375 - G4MagIntegratorStepper 1143 - G4MagIntegratorStepper 1376 * Include the SetEquationOfMotion() method << 1144 * Include the SetEquationOfMotion() method in G4MagIntegratorStepper 1377 (for Geant4e, as in field-V08-01-00) 1145 (for Geant4e, as in field-V08-01-00) 1378 * Update of disclaimers. 1146 * Update of disclaimers. 1379 1147 1380 November 13th, 2006 G.Cosmo - field-V08-01-0 1148 November 13th, 2006 G.Cosmo - field-V08-01-01 1381 ---------------------------- 1149 ---------------------------- 1382 - G4FieldTrack: 1150 - G4FieldTrack: 1383 o Implemented missing definition for GetCha 1151 o Implemented missing definition for GetCharge(). 1384 o Corrected 'const' qualifier for accessors 1152 o Corrected 'const' qualifier for accessors in the subclass. 1385 o Cleanup ... 1153 o Cleanup ... 1386 - NOTE: this tag is NOT including the change 1154 - NOTE: this tag is NOT including the change below in G4MagIntegratorStepper. 1387 1155 1388 June 21st, 2006 J.Apostolakis - field-V08-01 1156 June 21st, 2006 J.Apostolakis - field-V08-01-00 1389 ------------------------------ 1157 ------------------------------ 1390 - Including the SetEquationOfMotion() method 1158 - Including the SetEquationOfMotion() method in G4MagIntegratorStepper 1391 1159 1392 June 22nd, 2006 J.Apostolakis - field-V08-00 1160 June 22nd, 2006 J.Apostolakis - field-V08-00-07 1393 ----------------------------- 1161 ----------------------------- 1394 - G4ExactHelixStepper: reintroduced corrected << 1162 - G4ExactHelixStepper: reintroduced corrected class 1395 --> testPropagateMagField now tests it al 1163 --> testPropagateMagField now tests it also (added) 1396 - G4MagHelicalStepper: made Stepper virtual ( << 1164 - G4MagHelicalStepper: made Stepper virtual (for now) 1397 to enable G4ExactHelixStepper to 1165 to enable G4ExactHelixStepper to override it 1398 1166 1399 June 21st, 2006 J.Apostolakis - field-V08-00 1167 June 21st, 2006 J.Apostolakis - field-V08-00-06 1400 ----------------------------- 1168 ----------------------------- 1401 - G4ExactHelixStepper: taking out class << 1169 - G4ExactHelixStepper: taking out class 1402 (provisional tag in case of difficulties w 1170 (provisional tag in case of difficulties with class, revisions) 1403 1171 1404 --- The field-V08-00-0x , x > 6, tags do not 1172 --- The field-V08-00-0x , x > 6, tags do not include the change below 1405 1173 1406 May 27th, 2006 J.Apostolakis - field-V08-00- 1174 May 27th, 2006 J.Apostolakis - field-V08-00-05 1407 ----------------------------- 1175 ----------------------------- 1408 - G4FieldTrack: Fix for simple compilation er 1176 - G4FieldTrack: Fix for simple compilation error on Sun. 1409 1177 1410 May 23rd, 2006 J.Apostolakis - field-V08-00- 1178 May 23rd, 2006 J.Apostolakis - field-V08-00-04 1411 ----------------------------- 1179 ----------------------------- 1412 - G4FieldTrack: Fix for memory leak for Charg 1180 - G4FieldTrack: Fix for memory leak for ChargeState 1413 Made G4ChargeState a data member (has an i 1181 Made G4ChargeState a data member (has an instance) 1414 1182 1415 May 18th, 2006 J.Apostolakis - field-V08-00- 1183 May 18th, 2006 J.Apostolakis - field-V08-00-03 1416 ----------------------------- 1184 ----------------------------- 1417 - G4FieldTrack: Fix for issue with gcc 4.0 ( 1185 - G4FieldTrack: Fix for issue with gcc 4.0 (mac) 1418 1186 1419 May 16th, 2006 J.Apostolakis - field-V08-00- 1187 May 16th, 2006 J.Apostolakis - field-V08-00-02 1420 ----------------------------- 1188 ----------------------------- 1421 - G4ExactHelixStepper 1189 - G4ExactHelixStepper 1422 Taken out alpha version - it was included << 1190 Taken out alpha version - it was included in error 1423 1191 1424 May 11th, 2006 J.Apostolakis - field-V08-00- 1192 May 11th, 2006 J.Apostolakis - field-V08-00-01 1425 ----------------------------- 1193 ----------------------------- 1426 - G4FieldTrack 1194 - G4FieldTrack 1427 Significant, provisional change: 1195 Significant, provisional change: 1428 * Added charge and moments to FieldTrack << 1196 * Added charge and moments to FieldTrack 1429 (SetMethods, and pointer to embedde 1197 (SetMethods, and pointer to embedded class G4ChargeState) 1430 Technical revisions: 1198 Technical revisions: 1431 * New constructor, with charge and witho 1199 * New constructor, with charge and without velocity 1432 * New almost-default constructor (taking 1200 * New almost-default constructor (taking a character) 1433 Also: << 1201 Also: 1434 * New constructor to take const referenc 1202 * New constructor to take const reference to Spin 1435 * New implementation method UpdateFourMo 1203 * New implementation method UpdateFourMomentum 1436 * Streaming operator takes const referen << 1204 * Streaming operator takes const reference, 1437 and outputs momentum magnitude & dir 1205 and outputs momentum magnitude & direction magn. 1438 1206 1439 - G4ExactHelixStepper 1207 - G4ExactHelixStepper 1440 Included in error << 1208 Included in error 1441 1209 1442 Apr 6th, 2006 J.Apostolakis - field-V08-00-0 1210 Apr 6th, 2006 J.Apostolakis - field-V08-00-00 1443 ---------------------------- 1211 ---------------------------- 1444 - In G4FieldManager 1212 - In G4FieldManager 1445 Added missing const GetChordFinder() method << 1213 Added missing const GetChordFinder() method -- Problem report #847 1446 Corrected description (at top of include fil 1214 Corrected description (at top of include file hh): 1447 - It is possible to set the field manager o 1215 - It is possible to set the field manager of a logical volume (since long) 1448 1216 1449 1217 1450 Nov 9th, 2005 G.Cosmo - field-V07-01-00 1218 Nov 9th, 2005 G.Cosmo - field-V07-01-00 1451 ---------------------- 1219 ---------------------- 1452 - Archived source code of field01 unit test, 1220 - Archived source code of field01 unit test, now moved to examples/extended. 1453 - Migrated unit test NTST to use <sstream> an 1221 - Migrated unit test NTST to use <sstream> and std:: namespace. Up-to-date 1454 with G4VParticleChange new interface. 1222 with G4VParticleChange new interface. 1455 1223 1456 Dec 2nd, 2004 G.Cosmo - field-V06-02-02 1224 Dec 2nd, 2004 G.Cosmo - field-V06-02-02 1457 ---------------------- 1225 ---------------------- 1458 - Implemented migration to <cmath> for standa 1226 - Implemented migration to <cmath> for standard mathematical functions. 1459 1227 1460 Dec 2nd, 2004 J.Apostolakis - field-V06-02-0 1228 Dec 2nd, 2004 J.Apostolakis - field-V06-02-01 1461 ---------------------------- 1229 ---------------------------- 1462 - Improvements to comments in key Field class 1230 - Improvements to comments in key Field classes 1463 - Changes to test/NTST, renaming 'main()' fil 1231 - Changes to test/NTST, renaming 'main()' file & modifying input files 1464 - Several improvements to test files 1232 - Several improvements to test files 1465 testProElectroMagField.cc 1233 testProElectroMagField.cc 1466 * updating for new Electric Field inheri 1234 * updating for new Electric Field inheritance (& need for 1467 different ChordFinder creation ) 1235 different ChordFinder creation ) 1468 * adding columns to print-outs in spin c 1236 * adding columns to print-outs in spin case for better diagnosis 1469 testPropagateMagField.cc 1237 testPropagateMagField.cc 1470 * Added verbosity level 1 to field Propa 1238 * Added verbosity level 1 to field Propagator, relabelled old method. 1471 testPropagateSpin.cc 1239 testPropagateSpin.cc 1472 * Added more printout of parameters (sta 1240 * Added more printout of parameters (state) of G4PropagatorInField 1473 and G4FieldManager before running each 1241 and G4FieldManager before running each test. 1474 test/NTST 1242 test/NTST 1475 * Renamed main() file to testNTST.cc and 1243 * Renamed main() file to testNTST.cc and added G4UItcsh terminal. 1476 1244 1477 1245 1478 Oct 20th, 2004 P.Gumplinger - field-V06-02-0 1246 Oct 20th, 2004 P.Gumplinger - field-V06-02-00 1479 ---------------------------- 1247 ---------------------------- 1480 - Added Set/GetAnomaly() methods in G4Mag_Spi 1248 - Added Set/GetAnomaly() methods in G4Mag_SpinEqRhs. 1481 1249 1482 Feb 9th, 2004 G.Cosmo - field-V06-00-00 1250 Feb 9th, 2004 G.Cosmo - field-V06-00-00 1483 ---------------------- 1251 ---------------------- 1484 - G4UniformMagField.cc: added missing impleme 1252 - G4UniformMagField.cc: added missing implementation for 1485 GetConstantFieldValue(). Fixes bug report # 1253 GetConstantFieldValue(). Fixes bug report #575. 1486 1254 1487 Nov 13th, 2003 J.Apostolakis - field-V05-02- 1255 Nov 13th, 2003 J.Apostolakis - field-V05-02-08 1488 ----------------------------- 1256 ----------------------------- 1489 - G4ChordFinder: << 1257 - G4ChordFinder: 1490 o extended the arguments of AdvanceChordLim 1258 o extended the arguments of AdvanceChordLimited and FindNextChord 1491 adding the current safety sphere (origin, 1259 adding the current safety sphere (origin, radius), for use in optimising. 1492 o fixed set functions which had side effect 1260 o fixed set functions which had side effect (on first fraction.) 1493 - Added new class G4ChordFinderSaf (using abo 1261 - Added new class G4ChordFinderSaf (using above) 1494 o which uses safety to create longer chords << 1262 o which uses safety to create longer chords, so enabling 1495 propagation in fewer steps --> performa 1263 propagation in fewer steps --> performance optimisation 1496 - G4FieldManager: << 1264 - G4FieldManager: 1497 Changed default accuracy parameters to mu 1265 Changed default accuracy parameters to much better values: 1498 Epsilon Min Default to 5.0e-5 from 1266 Epsilon Min Default to 5.0e-5 from 5.0e-7 (too small) 1499 Epsilon Max Default to 0.001 from 1267 Epsilon Max Default to 0.001 from 0.05 (far too big) 1500 1268 1501 Nov 7th, 2003 J.Apostolakis - field-V05-02- 1269 Nov 7th, 2003 J.Apostolakis - field-V05-02-07 1502 ----------------------------- 1270 ----------------------------- 1503 - G4ChordFinder << 1271 - G4ChordFinder 1504 o Made this a base class 1272 o Made this a base class 1505 * to enable specialised implementations 1273 * to enable specialised implementations 1506 * virtual methods: FindNextChord, PrintSt 1274 * virtual methods: FindNextChord, PrintStatistics) 1507 * new method AccumulateStatistics to aid 1275 * new method AccumulateStatistics to aid in this 1508 1276 1509 o New default value of 1277 o New default value of Accuracy!! 1510 Delta Chord to 0.25 mm (from 1278 Delta Chord to 0.25 mm (from 3.00 mm) <----- 1511 << 1279 1512 - G4FieldManager 1280 - G4FieldManager 1513 o Changed default values of 1281 o Changed default values of 1514 Delta One Step to 0.010 mm (from 1282 Delta One Step to 0.010 mm (from 0.25 mm) <----- 1515 Delta Intersection to 0.001 mm (from 1283 Delta Intersection to 0.001 mm (from 0.10 mm) <----- 1516 for increased default accuracy 1284 for increased default accuracy 1517 o Use field's DoesFieldChangesEnergy method 1285 o Use field's DoesFieldChangesEnergy method to set relevant member variable 1518 1286 1519 Nov 7th, 2003 J.Apostolakis - field-V05-02-0 1287 Nov 7th, 2003 J.Apostolakis - field-V05-02-06 1520 ---------------------------- 1288 ---------------------------- 1521 - G4ChordFinder 1289 - G4ChordFinder 1522 Created new class members for parameters u 1290 Created new class members for parameters used in heuristics: 1523 fFirstFraction for first st 1291 fFirstFraction for first step 1524 fFractionLast fraction of 1292 fFractionLast fraction of last step, if redoing 1525 fFractionNextEstimate confidence i 1293 fFractionNextEstimate confidence in step est from chord 1526 fMultipleRadius for small ra 1294 fMultipleRadius for small radii of curve 1527 New method to Print Statistics. 1295 New method to Print Statistics. 1528 Verbosity to enable optional printing of s 1296 Verbosity to enable optional printing of stats in destructor. 1529 1297 1530 Nov 7th, 2003 J.Apostolakis - field-V05-02-0 1298 Nov 7th, 2003 J.Apostolakis - field-V05-02-05 1531 ---------------------------- 1299 ---------------------------- 1532 - G4ChordFinder 1300 - G4ChordFinder 1533 o Added simple statistics (with access meth 1301 o Added simple statistics (with access methods) 1534 - test/NTST 1302 - test/NTST 1535 o Added new test, taken from D. Williams' e 1303 o Added new test, taken from D. Williams' exampleNTST 1536 modified and maintained by J. Apostolaki << 1304 modified and maintained by J. Apostolakis 1537 1305 1538 Nov 5th, 2003 J.Apostolakis - field-V05-02-0 1306 Nov 5th, 2003 J.Apostolakis - field-V05-02-04 1539 ---------------------------- 1307 ---------------------------- 1540 - G4FieldManager 1308 - G4FieldManager 1541 o Added Get/Set Min/Max imumEpsilonStep 1309 o Added Get/Set Min/Max imumEpsilonStep 1542 o Added fEpsilonMin fEpsilonMax & default t << 1310 o Added fEpsilonMin fEpsilonMax & default to member variables. 1543 >> Now these can be tailored for different 1311 >> Now these can be tailored for different parts of a detector or setup 1544 1312 1545 - G4Field G4ElectroMagneticField G4MagneticFi 1313 - G4Field G4ElectroMagneticField G4MagneticField 1546 o Adding the virtual method DoesFieldChan 1314 o Adding the virtual method DoesFieldChangeEnergy() [to be used later] 1547 o So that G4ElectroMagneticField derives fr 1315 o So that G4ElectroMagneticField derives from G4Field 1548 and G4MagneticField derives from G4Electr 1316 and G4MagneticField derives from G4ElectroMagneticField 1549 1317 1550 - G4MagErrorStepper and derived classes 1318 - G4MagErrorStepper and derived classes 1551 (G4ClassicalRK4 G4SimpleRunge G4ExplicitE 1319 (G4ClassicalRK4 G4SimpleRunge G4ExplicitEuler G4ImplicitEuler G4SimpleHeum) 1552 o the constructors take G4EquationOfMotion 1320 o the constructors take G4EquationOfMotion in place of a G4Mag_EqRhs 1553 as it these applicable to all types of eq << 1321 as it these applicable to all types of equations 1554 (not just ones for Magnetic field) 1322 (not just ones for Magnetic field) 1555 1323 1556 - G4EqMagElectricField 1324 - G4EqMagElectricField 1557 o derives from G4EquationOfMotion 1325 o derives from G4EquationOfMotion 1558 1326 1559 - G4ChordFinder 1327 - G4ChordFinder 1560 o Moved simple constructor to .cc file (tog 1328 o Moved simple constructor to .cc file (together with 'complex' one) 1561 o Changed fDefaultDeltaChord to a constant 1329 o Changed fDefaultDeltaChord to a constant member variable (from class var.) 1562 1330 1563 Nov 4th, 2003 G.Cosmo - field-V05-02-03 1331 Nov 4th, 2003 G.Cosmo - field-V05-02-03 1564 ---------------------- 1332 ---------------------- 1565 Undocumented. 1333 Undocumented. 1566 1334 1567 Nov 4th, 2003 G.Cosmo - field-V05-02-02 1335 Nov 4th, 2003 G.Cosmo - field-V05-02-02 1568 ---------------------- 1336 ---------------------- 1569 - G4ChordFinder.cc: 1337 - G4ChordFinder.cc: 1570 o Fixed compilation problem occuring when G 1338 o Fixed compilation problem occuring when G4DEBUG_FIELD was set. 1571 1339 1572 Oct 31st, 2003 G.Cosmo - field-V05-02-01 1340 Oct 31st, 2003 G.Cosmo - field-V05-02-01 1573 ----------------------- 1341 ----------------------- 1574 - Code cleanup and removal of redundant heade 1342 - Code cleanup and removal of redundant header includes, replaced 1575 where necessary by forward declarations. 1343 where necessary by forward declarations. 1576 1344 1577 Oct 24th, 2003 G.Cosmo - field-V05-02-00 1345 Oct 24th, 2003 G.Cosmo - field-V05-02-00 1578 ----------------------- 1346 ----------------------- 1579 - Migrated to use new G4Exception scheme. 1347 - Migrated to use new G4Exception scheme. 1580 - Updated outputs of unit tests. 1348 - Updated outputs of unit tests. 1581 1349 1582 Jun 25th, 2003 J.Apostolakis - field-V05-01- 1350 Jun 25th, 2003 J.Apostolakis - field-V05-01-06 1583 ----------------------------- 1351 ----------------------------- 1584 - Suppressed printing of statistics (by defau 1352 - Suppressed printing of statistics (by default) in MagIntegrationDriver. 1585 - Fix for unused variables warning in ChordFi 1353 - Fix for unused variables warning in ChordFinder. 1586 1354 1587 Jun 24th, 2003 J.Apostolakis - field-V05-01- 1355 Jun 24th, 2003 J.Apostolakis - field-V05-01-05 1588 ----------------------------- 1356 ----------------------------- 1589 - Fix for problems introduced in field-V05-01 1357 - Fix for problems introduced in field-V05-01-03: 1590 o G4ChordFinder: 1358 o G4ChordFinder: 1591 - Fixed unitialised variable for next ste 1359 - Fixed unitialised variable for next step (now used for optimisation). 1592 - Protected assignment of 'fLastStepEstim 1360 - Protected assignment of 'fLastStepEstimate_Unconstrained' in case of 1593 unitialised value. 1361 unitialised value. 1594 1362 1595 Jun 21st, 2003 J.Apostolakis - field-V05-01- 1363 Jun 21st, 2003 J.Apostolakis - field-V05-01-04 1596 ----------------------------- 1364 ----------------------------- 1597 - G4FieldManager: 1365 - G4FieldManager: 1598 o Revised into 'slightly' abstract class wi 1366 o Revised into 'slightly' abstract class with one (new) virtual method: 1599 ConfigureForTrack( const G4Track *) 1367 ConfigureForTrack( const G4Track *) 1600 to enable the user to configure it with d 1368 to enable the user to configure it with different accuracy parameters 1601 depending on energy, particle id, volume, 1369 depending on energy, particle id, volume, or other track state. 1602 Currently by default, no such configurati 1370 Currently by default, no such configuration action is undertaken. 1603 1371 1604 Jun 21st, 2003 J.Apostolakis - field-V05-01- 1372 Jun 21st, 2003 J.Apostolakis - field-V05-01-03 1605 ----------------------------- 1373 ----------------------------- 1606 - Revised G4ChordFinder: 1374 - Revised G4ChordFinder: 1607 * FindNextChord now provides next Step esti 1375 * FindNextChord now provides next Step estimate (for accuracy in Integr.) 1608 * AdvanceChordLimited to use 1376 * AdvanceChordLimited to use 1609 - above extra argument to FindNextChord 1377 - above extra argument to FindNextChord 1610 - the next Step argument as initial-step 1378 - the next Step argument as initial-step suggested to Accurate Advance 1611 - G4MagIntegratorDriver: 1379 - G4MagIntegratorDriver: 1612 * Added extra argument to Accurate Advance: 1380 * Added extra argument to Accurate Advance: 1613 hinitial: suggested initial step size 1381 hinitial: suggested initial step size 1614 - Note: 1382 - Note: 1615 * Purpose: performance boost with small/li << 1383 * Purpose: performance boost with small/limited changes, 1616 including minor revisions of interfaces. 1384 including minor revisions of interfaces. 1617 1385 1618 Jun 20th, 2003 J.Apostolakis - field-V05-01- 1386 Jun 20th, 2003 J.Apostolakis - field-V05-01-02 1619 ----------------------------- 1387 ----------------------------- 1620 - G4MagIntegratorDriver Quick Advance: 1388 - G4MagIntegratorDriver Quick Advance: 1621 o Took out two unneccessary square roots an 1389 o Took out two unneccessary square roots and two divisions. 1622 - Notes: 1390 - Notes: 1623 1) Pure magnetic, electric, spin field test 1391 1) Pure magnetic, electric, spin field test are completed. 1624 2) Purpose: performance boost with small c 1392 2) Purpose: performance boost with small changes. 1625 1393 1626 Jun 20th, 2003 J.Apostolakis - field-V05-01- 1394 Jun 20th, 2003 J.Apostolakis - field-V05-01-01 1627 ----------------------------- 1395 ----------------------------- 1628 - G4MagIntegratorDriver: 1396 - G4MagIntegratorDriver: 1629 o Corrected neglected conversion to <std> f 1397 o Corrected neglected conversion to <std> from "G4std". 1630 - Migrated unit tests to <std>. 1398 - Migrated unit tests to <std>. 1631 1399 1632 Jun 19th, 2003 J.Apostolakis - field-V05-01- 1400 Jun 19th, 2003 J.Apostolakis - field-V05-01-00 1633 ----------------------------- << 1401 ----------------------------- 1634 - G4ChordFinder: 1402 - G4ChordFinder: 1635 o Corrected simple constructor to initialis 1403 o Corrected simple constructor to initialise several member variables 1636 including LastStepEstimate_Unconstrained. 1404 including LastStepEstimate_Unconstrained. 1637 - G4MagIntegratorDriver: 1405 - G4MagIntegratorDriver: 1638 o Keep simple statistics on number of steps 1406 o Keep simple statistics on number of steps, total error and step size. 1639 o Took out an unneccessary square root (in 1407 o Took out an unneccessary square root (in loop). 1640 - NOTE: tests field0X are to be migrated. 1408 - NOTE: tests field0X are to be migrated. 1641 Purpose: fixes, small performance boost, si 1409 Purpose: fixes, small performance boost, simple statistics. 1642 1410 1643 Apr 2nd, 2003 G. Cosmo - field-V05-00-01 1411 Apr 2nd, 2003 G. Cosmo - field-V05-00-01 1644 --------------------------- 1412 --------------------------- 1645 - Created .cc files for G4MagneticField and G 1413 - Created .cc files for G4MagneticField and G4ElecroMagneticField. 1646 - Fixed definition of copy constructors in fi 1414 - Fixed definition of copy constructors in fields inheritance chain. 1647 - Moved obsolete method StepWithEst() to priv 1415 - Moved obsolete method StepWithEst() to private in G4RKG3_Stepper, 1648 G4CashKarpRKF45 and G4ClassicalRK4. 1416 G4CashKarpRKF45 and G4ClassicalRK4. 1649 - Cleared strict pedantic warnings. 1417 - Cleared strict pedantic warnings. 1650 1418 1651 Feb 5th, 2003 G. Cosmo - field-V05-00-00 1419 Feb 5th, 2003 G. Cosmo - field-V05-00-00 1652 --------------------------- 1420 --------------------------- 1653 - Removed redundant inclusions of PhysicalCon 1421 - Removed redundant inclusions of PhysicalConstants.h or SystemOfUnits.h. 1654 Now using directly CLHEP's headers. 1422 Now using directly CLHEP's headers. 1655 1423 1656 Dec 12th, 2002 G. Cosmo - field-V04-01-06 1424 Dec 12th, 2002 G. Cosmo - field-V04-01-06 1657 --------------------------- 1425 --------------------------- 1658 - Moved unused variable 'nStpPr' to debug sec 1426 - Moved unused variable 'nStpPr' to debug section in G4MagIntegratorDriver.cc, 1659 method AccurateAdvance(). 1427 method AccurateAdvance(). 1660 1428 1661 Dec 06th, 2002 G. Cosmo - field-V04-01-05 1429 Dec 06th, 2002 G. Cosmo - field-V04-01-05 1662 --------------------------- 1430 --------------------------- 1663 - Merged changes by Makoto for migration to n 1431 - Merged changes by Makoto for migration to new G4ApplicationState in 1664 messengers of unit tests. 1432 messengers of unit tests. 1665 1433 1666 Nov 29th, 2002 J. Apostolakis - field-V04 1434 Nov 29th, 2002 J. Apostolakis - field-V04-01-04 1667 --------------------------------- 1435 --------------------------------- 1668 Further fixes for unitialised variable - in o 1436 Further fixes for unitialised variable - in other steppers: 1669 G4ExplicitEuler, G4ImplicitEuler, G4ImplicitE 1437 G4ExplicitEuler, G4ImplicitEuler, G4ImplicitEuler 1670 1438 1671 Large memory leak fix for G4ImplicitEuler, sm 1439 Large memory leak fix for G4ImplicitEuler, small for ExplicitEuler. 1672 Test accuracy requirements for spin adjusted. 1440 Test accuracy requirements for spin adjusted. 1673 1441 1674 Nov 20th, 2002 J. Apostolakis - field-V04 1442 Nov 20th, 2002 J. Apostolakis - field-V04-01-03 1675 --------------------------------- 1443 --------------------------------- 1676 Further fixes for unitialised variables: 1444 Further fixes for unitialised variables: 1677 1445 1678 * G4MagIntegratorStepper 1446 * G4MagIntegratorStepper 1679 - Created instance variables for number of 1447 - Created instance variables for number of Integration, State variables 1680 1448 1681 * G4MagErrorStepper 1449 * G4MagErrorStepper 1682 - Utilise number of State variables to ini 1450 - Utilise number of State variables to initialise full state of Output. 1683 1451 1684 * G4Mag_SpinEqRhs 1452 * G4Mag_SpinEqRhs 1685 - Initialise the values of dydx that we do 1453 - Initialise the values of dydx that we do not update. 1686 1454 1687 Note: Undoing these changes and the ones belo << 1455 Note: Undoing these changes and the ones below would enable the user to see 1688 if undefined values are used in varia << 1456 if undefined values are used in variables that he/she expects to 1689 have integrated. 1457 have integrated. 1690 BUT the fix assists us to find truly un 1458 BUT the fix assists us to find truly undefined values elsewhere 1691 -- if they exist -- using valgrind. 1459 -- if they exist -- using valgrind. 1692 1460 1693 Nov 9th, 2002 J. Apostolakis - field-V04 1461 Nov 9th, 2002 J. Apostolakis - field-V04-01-02 1694 --------------------------------- 1462 --------------------------------- 1695 Fixes to address use of uninitialised variabl 1463 Fixes to address use of uninitialised variable in seting values of FieldTrack 1696 (this was done by design, but was neither ide 1464 (this was done by design, but was neither ideal or always well-behaved.) 1697 Changed: 1465 Changed: 1698 1466 1699 * G4MagIntegratorDriver 1467 * G4MagIntegratorDriver 1700 - Added instance variables for number of In 1468 - Added instance variables for number of Integration, State variables 1701 - Initialise all State variables of Output 1469 - Initialise all State variables of Output with values from initial variables. 1702 - Only the Integration variables will be in 1470 - Only the Integration variables will be integrated 1703 1471 1704 * G4ClassicalRK4 and G4MagErrorStepper 1472 * G4ClassicalRK4 and G4MagErrorStepper 1705 - Copy the initial time to intermediate and << 1473 - Copy the initial time to intermediate and output C-vectors, 1706 in order to provide a zero-eth estimate 1474 in order to provide a zero-eth estimate value. 1707 Note1: Should clean-up this fix, making it 1475 Note1: Should clean-up this fix, making it similar to Driver (above) 1708 Note2: Should *not* be used for time varyin << 1476 Note2: Should *not* be used for time varying field unless created 1709 with at least 8 integration variab 1477 with at least 8 integration variables 1710 --------------------------------- << 1478 --------------------------------- 1711 G4ChordFinder.cc 1479 G4ChordFinder.cc 1712 - Correction from zero dChord 1480 - Correction from zero dChord 1713 which is does not update the Unconstrai 1481 which is does not update the Unconstrained step estimate (state) 1714 - Added instance variables to store number 1482 - Added instance variables to store number of calls and chord-trials 1715 and new print method for 'debugging' 1483 and new print method for 'debugging' 1716 1484 1717 G4Mag_UsualEqRhs.cc 1485 G4Mag_UsualEqRhs.cc 1718 - Attempted approximations for speedup *ov 1486 - Attempted approximations for speedup *overwritten* by full inverse sqrt. 1719 1487 1720 1488 1721 Oct 29th, 2002 J. Apostolakis - test-V04-0 1489 Oct 29th, 2002 J. Apostolakis - test-V04-01-02 1722 -------------------------------- 1490 -------------------------------- 1723 - Changed output of tests for small output ch 1491 - Changed output of tests for small output change (GeV->MeV) 1724 1492 1725 Oct 29th, 2002 J. Apostolakis - test-V04-0 1493 Oct 29th, 2002 J. Apostolakis - test-V04-01-02 1726 -------------------------------- 1494 -------------------------------- 1727 - (Only changed and tagged test subdirectory) 1495 - (Only changed and tagged test subdirectory) 1728 Changed tests to use the new method SetChar << 1496 Changed tests to use the new method SetChargeMomentumMass of 1729 PropagatorInField in order to be ready for << 1497 PropagatorInField in order to be ready for its full use in 1730 G4PropagatorInField. 1498 G4PropagatorInField. 1731 (A key problem at end-July was this non-pre 1499 (A key problem at end-July was this non-preparedness of the unit tests.) 1732 1500 1733 Jul 24th, 2002 G. Cosmo - field-V04-01-01 1501 Jul 24th, 2002 G. Cosmo - field-V04-01-01 1734 -------------------------- 1502 -------------------------- 1735 - G4FieldManager[.hh.cc]: made fDefault_Delta 1503 - G4FieldManager[.hh.cc]: made fDefault_Delta_One_Step_Value and 1736 fDefault_Delta_Intersection_Val normal cons 1504 fDefault_Delta_Intersection_Val normal constants (not static) to 1737 resolve initialisation problems. Coworks wi 1505 resolve initialisation problems. Coworks with tag "geomvol-V04-01-01". 1738 Fixes problem report #383. 1506 Fixes problem report #383. 1739 1507 1740 Jul 23rd, 2002 G. Cosmo - field-V04-01-00 1508 Jul 23rd, 2002 G. Cosmo - field-V04-01-00 1741 -------------------------- 1509 -------------------------- 1742 - G4MagIntegratorDriver.cc: cleared warning o 1510 - G4MagIntegratorDriver.cc: cleared warning of unused variable. 1743 Protected variable 'nStpPr' with G4DEBUG_FI 1511 Protected variable 'nStpPr' with G4DEBUG_FIELD. 1744 - Updated reference outputs and some other up 1512 - Updated reference outputs and some other updates to unit tests. 1745 1513 1746 Jun 25th, 2002 J. Apostolakis - field-V04- 1514 Jun 25th, 2002 J. Apostolakis - field-V04-00-05 1747 -------------------------------- 1515 -------------------------------- 1748 - G4ChordFinder.cc 1516 - G4ChordFinder.cc 1749 In NewStep: fixed case of stepTrial >> ste 1517 In NewStep: fixed case of stepTrial >> stepTrialOld, to take large multiple - not small; 1750 (it was treated the same as stepTrial 1518 (it was treated the same as stepTrial << old) 1751 Added guard for stepTrial==0.0 in FindNextC 1519 Added guard for stepTrial==0.0 in FindNextChord 1752 - G4MagIntegratorDriver.cc 1520 - G4MagIntegratorDriver.cc 1753 Added protection against error of 0.0 in Co 1521 Added protection against error of 0.0 in ComputeNewStepSize. 1754 1522 1755 - G4EquationOfMotion 1523 - G4EquationOfMotion 1756 Made RightHandSide an inline method << 1524 Made RightHandSide an inline method 1757 In GetFieldValue signature of Point is doub 1525 In GetFieldValue signature of Point is double[4], instead of double[3] 1758 Defined G4maximum_number_of_field_component 1526 Defined G4maximum_number_of_field_components as enum 1759 1527 1760 Jun 7th, 2002 J. Apostolakis - field-V04- 1528 Jun 7th, 2002 J. Apostolakis - field-V04-00-04 1761 -------------------------------- 1529 -------------------------------- 1762 - G4Mag_SpinEqRhs.cc 1530 - G4Mag_SpinEqRhs.cc 1763 o Fix to momentum component of Rhs. 1531 o Fix to momentum component of Rhs. 1764 1532 1765 May 31st, 2002 J. Apostolakis - field-V04- 1533 May 31st, 2002 J. Apostolakis - field-V04-00-03 1766 -------------------------------- 1534 -------------------------------- 1767 - G4MagHelicalStepper.icc 1535 - G4MagHelicalStepper.icc 1768 o Fix to LinearStep (correcting previous 1536 o Fix to LinearStep (correcting previous fix/cosmetic) 1769 - G4ChordFinder.cc 1537 - G4ChordFinder.cc 1770 o Corrective behaviour for very small or 1538 o Corrective behaviour for very small or large changes of 'new step' 1771 [observed in complex 3-d fields (purgi 1539 [observed in complex 3-d fields (purging magnet), or 1772 approximate equations of 1540 approximate equations of motion ] 1773 - G4MagIntegratorDriver 1541 - G4MagIntegratorDriver 1774 o More information in debug printing (tan 1542 o More information in debug printing (tangent cosine, norm) 1775 1543 1776 May 7th, 2002 J. Apostolakis 1544 May 7th, 2002 J. Apostolakis 1777 -------------------------------- 1545 -------------------------------- 1778 - G4MagIntegratorDriver - ( on trunk ) 1546 - G4MagIntegratorDriver - ( on trunk ) 1779 * Fix to update the field track for steps s << 1547 * Fix to update the field track for steps smaller than the chord finder's 1780 minstep. << 1548 minstep. 1781 1549 1782 1550 1783 April 19th, 2002 J. Apostolakis (not inclu 1551 April 19th, 2002 J. Apostolakis (not included in tag field-V04-00-02) 1784 -------------------------------- 1552 -------------------------------- 1785 - G4MagIntegratorDriver.hh,cc (on trunk) 1553 - G4MagIntegratorDriver.hh,cc (on trunk) 1786 * Bug fix: last step can now only be trigg 1554 * Bug fix: last step can now only be triggerred if x+h > x2. 1787 * Debugging methods improved. 1555 * Debugging methods improved. 1788 1556 1789 May 7th, 2002 J. Apostolakis - field-V04- 1557 May 7th, 2002 J. Apostolakis - field-V04-00-02 1790 -------------------------------- 1558 -------------------------------- 1791 Tag including only two changes below: 1559 Tag including only two changes below: 1792 - Putting aside the other fixes/changes (as << 1560 - Putting aside the other fixes/changes (as the one above) to 1793 G4MagIntegratorDriver, tagging only the fix t << 1561 G4MagIntegratorDriver, tagging only the fix to update the field track 1794 for steps smaller than the chord finder's min 1562 for steps smaller than the chord finder's minstep. 1795 1563 1796 May 7th, 2002 J. Apostolakis 1564 May 7th, 2002 J. Apostolakis 1797 -------------------------------- 1565 -------------------------------- 1798 - G4MagIntegratorDriver - ** branch ** on 1566 - G4MagIntegratorDriver - ** branch ** on 1.20 ** 1799 * Fix to update the field track for steps s << 1567 * Fix to update the field track for steps smaller than the chord finder's 1800 minstep. << 1568 minstep. 1801 1569 1802 May 7th, 2002 J. Apostolakis 1570 May 7th, 2002 J. Apostolakis 1803 -------------------------------- 1571 -------------------------------- 1804 - G4MagIntegratorStepper 1572 - G4MagIntegratorStepper 1805 * Changed RightHandSide from a virtual meth << 1573 * Changed RightHandSide from a virtual method to an inline (non-virtual) 1806 method -- it is not and should not be c 1574 method -- it is not and should not be changed by derived classes. 1807 1575 1808 Mar 28th, 2002 J. Apostolakis - field-V04-00 1576 Mar 28th, 2002 J. Apostolakis - field-V04-00-01 1809 ----------------------- 1577 ----------------------- 1810 - Fixed LinearStep in G4MagHelicalStepper (ex << 1578 - Fixed LinearStep in G4MagHelicalStepper (expected unit direction), 1811 that had not been updated for new integra 1579 that had not been updated for new integration variable (momentum). 1812 --> found in testing helical steppers in << 1580 --> found in testing helical steppers in BaBar setup, 1813 in region of negligible field. 1581 in region of negligible field. 1814 1582 1815 - Changed maximum number of steps in driver, << 1583 - Changed maximum number of steps in driver, 1816 halving from 500/order to 250/order. 1584 halving from 500/order to 250/order. 1817 1585 1818 Mar 28th, 2002 J. Apostolakis 1586 Mar 28th, 2002 J. Apostolakis 1819 ----------------------- 1587 ----------------------- 1820 - Committed directory "OtherField" with examp 1588 - Committed directory "OtherField" with examples of fields, 1821 Delphi, Quadropole etc. 1589 Delphi, Quadropole etc. 1822 1590 1823 - Updated outputs of testPropagateMagField fo 1591 - Updated outputs of testPropagateMagField for new version of source, 1824 that has 3rd test case. 1592 that has 3rd test case. 1825 1593 1826 - Added new runs and outputs for testProElect << 1594 - Added new runs and outputs for testProElectroMagField, with 1827 several steppers, and modified test scrip 1595 several steppers, and modified test script test.sh accordingly. 1828 1596 1829 1597 1830 Feb 2nd, 2002 V. Grichine 1598 Feb 2nd, 2002 V. Grichine 1831 ------------------------ 1599 ------------------------ 1832 Update of test/field02 to utilise new feature << 1600 Update of test/field02 to utilise new features 1833 - for Electric Field, to signal that a fiel 1601 - for Electric Field, to signal that a field modifies the energy; 1834 - updating of time by integration, for time 1602 - updating of time by integration, for time dependent fields. 1835 1603 1836 Jan 10th, 2002 G. Cosmo - field-V04-00-00 1604 Jan 10th, 2002 G. Cosmo - field-V04-00-00 1837 ------------------------ 1605 ------------------------ 1838 - Activated G4DEBUG_FIELD in GNUmakefile. 1606 - Activated G4DEBUG_FIELD in GNUmakefile. 1839 1607 1840 Jan 7th, 2002 J. Apostolakis 1608 Jan 7th, 2002 J. Apostolakis 1841 ----------------------------- 1609 ----------------------------- 1842 - G4MagHelicalStepper.cc: renamed local varia 1610 - G4MagHelicalStepper.cc: renamed local variables "*Momentum*" to "*Velocity*". 1843 1611 1844 Dec 4th, 2001 G. Cosmo - field-V03-02-08 1612 Dec 4th, 2001 G. Cosmo - field-V03-02-08 1845 ----------------------- 1613 ----------------------- 1846 - Fixed interface to virtual functions in the 1614 - Fixed interface to virtual functions in the inheritance chain. 1847 Added 'virtual' qualifier to destructor and 1615 Added 'virtual' qualifier to destructor and GetFieldValue() for subclasses 1848 of G4Field which may act as base class for 1616 of G4Field which may act as base class for user-defined fields. 1849 Modified files: G4MagneticField.hh, G4Unifo 1617 Modified files: G4MagneticField.hh, G4UniformElectricField[.hh.cc], 1850 G4UniformMagField.hh. 1618 G4UniformMagField.hh. 1851 - G4ChordFinder.cc: the unphysical curve leng 1619 - G4ChordFinder.cc: the unphysical curve length exception has been moved to 1852 inside #ifdef G4DEBUG_FIELD. Other #ifdef r 1620 inside #ifdef G4DEBUG_FIELD. Other #ifdef renamed to G4DEBUG_FIELD. 1853 - G4LineSection[.hh.cc]: correction for A=B w 1621 - G4LineSection[.hh.cc]: correction for A=B when the distance from P to A 1854 is returned (V.Grichine). 1622 is returned (V.Grichine). 1855 1623 1856 Nov 28th, 2001 J. Apostolakis - field-V03-02 1624 Nov 28th, 2001 J. Apostolakis - field-V03-02-07 1857 ------------------------------ 1625 ------------------------------ 1858 G4FieldManager 1626 G4FieldManager 1859 - Added fFieldChangesEnergy data member and g 1627 - Added fFieldChangesEnergy data member and get/set methods 1860 ( required to allow transportation to know h 1628 ( required to allow transportation to know how to calculate end-time, 1861 enabling the integration of time only for 1629 enabling the integration of time only for electric fields. ) 1862 1630 1863 Nov 28th, 2001 J. Apostolakis - field-V03-02 1631 Nov 28th, 2001 J. Apostolakis - field-V03-02-06 1864 ------------------------------ 1632 ------------------------------ 1865 Small fixes only: 1633 Small fixes only: 1866 1634 1867 G4CashKarpRKF45 1635 G4CashKarpRKF45 1868 - fixed deletion of unassigned pointer/array 1636 - fixed deletion of unassigned pointer/array in destructor. (-> Volodya Velev) 1869 1637 1870 G4EqMagElectricField 1638 G4EqMagElectricField 1871 - bug fixed in the formula for inverse veloci 1639 - bug fixed in the formula for inverse velocity (-> V. Grichine) 1872 1640 1873 Nov 21st, 2001 G. Cosmo - field-V03-02-05 1641 Nov 21st, 2001 G. Cosmo - field-V03-02-05 1874 ------------------------ 1642 ------------------------ 1875 - Protected verbose warnings with G4FIELD_DEB 1643 - Protected verbose warnings with G4FIELD_DEBUG flag in G4ChordFinder.cc 1876 (V.Grichine). 1644 (V.Grichine). 1877 - Replaced G4DEBUG with G4DEBUG_FIELD in G4Ma 1645 - Replaced G4DEBUG with G4DEBUG_FIELD in G4MagIntegratorDriver.cc. 1878 - Updates to unit test field02. 1646 - Updates to unit test field02. 1879 1647 1880 Nov 13th, 2001 G. Cosmo - field-V03-02-04 1648 Nov 13th, 2001 G. Cosmo - field-V03-02-04 1881 ------------------------ 1649 ------------------------ 1882 - Removed warnings on gcc-2.95.2. 1650 - Removed warnings on gcc-2.95.2. 1883 - Included unit tests field01, field02, field 1651 - Included unit tests field01, field02, field03. 1884 1652 1885 Nov 9th, 2001 J. Apostolakis - field-V03-02- 1653 Nov 9th, 2001 J. Apostolakis - field-V03-02-03 1886 ----------------------------- 1654 ----------------------------- 1887 Motivated by the use case of an Electric Fiel 1655 Motivated by the use case of an Electric Field, 1888 - in order to achieve reasonable accuracy in 1656 - in order to achieve reasonable accuracy in time integration, 1889 - and to enable fields that vary with the gl 1657 - and to enable fields that vary with the global time) 1890 we (John Ap. and V. Grichine) changed the fol 1658 we (John Ap. and V. Grichine) changed the following: 1891 1659 1892 - Modified G4Field to allow the time in the G << 1660 - Modified G4Field to allow the time in the GetFieldValue method 1893 - as 4th component of the "position" arr 1661 - as 4th component of the "position" array. 1894 1662 1895 - Modified EquationOfMotion to provide the ti 1663 - Modified EquationOfMotion to provide the time to the GetFieldValue method 1896 of G4Field 1664 of G4Field 1897 * as the 4th component of a position / 1665 * as the 4th component of a position / time 4-array 1898 * in place of only position as a 3-arry 1666 * in place of only position as a 3-arry 1899 1667 1900 - Added the global time as variable to be int << 1668 - Added the global time as variable to be integrated over in the 1901 case of an Electric Field 1669 case of an Electric Field 1902 1670 1903 Support for these modifications is required i 1671 Support for these modifications is required in G4Transportation. 1904 1672 1905 (Better support will require either 1673 (Better support will require either 1906 - a modification in G4FieldTrack to indicate 1674 - a modification in G4FieldTrack to indicate whether the time has been 1907 integrated or not. << 1675 integrated or not. 1908 - or additionally making the time calculatio 1676 - or additionally making the time calculation a property of G4FieldTrack, 1909 bring in code to do it in from G4Transpor 1677 bring in code to do it in from G4Transportation.) 1910 1678 1911 1679 1912 Nov 8th, 2001 J. Apostolakis - field-V03-02 1680 Nov 8th, 2001 J. Apostolakis - field-V03-02-02 1913 ------------------------ 1681 ------------------------ 1914 - Introduced several changes in G4MagIntegrat 1682 - Introduced several changes in G4MagIntegratorDriver. 1915 1683 1916 It now respects its minimum driver step si 1684 It now respects its minimum driver step size in AccurateAdvance. 1917 * whereas before it did not. 1685 * whereas before it did not. 1918 1686 1919 Modified maximum number of steps, lowering 1687 Modified maximum number of steps, lowering by factor of 10. 1920 -> it should reduce the time spent integ 1688 -> it should reduce the time spent integrating low energy particles, 1921 in some cases where it is unneccessary. 1689 in some cases where it is unneccessary. 1922 1690 1923 Protected printing of (nearly?) all warnin 1691 Protected printing of (nearly?) all warnings with a G4DEBUG_FIELD #ifdef. 1924 1692 1925 1693 1926 Nov 05th, 2001 G. Cosmo - field-V03-02-01 1694 Nov 05th, 2001 G. Cosmo - field-V03-02-01 1927 ------------------------ 1695 ------------------------ 1928 - Fixed trivial warnings on Linux-g++. 1696 - Fixed trivial warnings on Linux-g++. 1929 1697 1930 Aug 27th, 2001 J. Apostolakis - field-V03-02 1698 Aug 27th, 2001 J. Apostolakis - field-V03-02-00 1931 ------------------------------ 1699 ------------------------------ 1932 Tagged last change. 1700 Tagged last change. 1933 This tag is utilised in transport-V03-02-01 ( 1701 This tag is utilised in transport-V03-02-01 (which requires it). 1934 1702 1935 Aug 16th, 2001 J. Apostolakis 1703 Aug 16th, 2001 J. Apostolakis 1936 ------------------------------ 1704 ------------------------------ 1937 G4ChordFinder: 1705 G4ChordFinder: 1938 Added method to erase/reset the state of the 1706 Added method to erase/reset the state of the StepEstimate. 1939 This is unwanted internal state that was car 1707 This is unwanted internal state that was carried between tracks and 1940 events as a side effect. 1708 events as a side effect. 1941 1709 1942 Calling this at the start of each track, it 1710 Calling this at the start of each track, it is now possible to 1943 eliminate a non-repeatability that was possi 1711 eliminate a non-repeatability that was possible due to this. 1944 ( It has been observed in particular in v 1712 ( It has been observed in particular in vaccuum. ) 1945 1713 1946 May 23rd, 2001 J. Apostolakis - field-V03-01 1714 May 23rd, 2001 J. Apostolakis - field-V03-01-02 1947 ------------------------------ 1715 ------------------------------ 1948 - Added Exception to G4RKG3Stepper constructo 1716 - Added Exception to G4RKG3Stepper constructor: the class does not work 1949 1717 1950 May 23rd, 2001 J. Apostolakis - field-V03-01 1718 May 23rd, 2001 J. Apostolakis - field-V03-01-01 1951 ------------------------------ 1719 ------------------------------ 1952 - Fixed G4FieldManager constructors to give d 1720 - Fixed G4FieldManager constructors to give default values to delta 1953 parameters. 1721 parameters. 1954 1722 1955 - Updated testPropagateSpin (it now compiles) << 1723 - Updated testPropagateSpin (it now compiles), 1956 and some outputs of testPropagateMagField 1724 and some outputs of testPropagateMagField (1 to 5) 1957 1725 1958 1726 1959 Apr 12th, 2001 P. Gumplinger - field-V03-01- 1727 Apr 12th, 2001 P. Gumplinger - field-V03-01-00 1960 ----------------------------- 1728 ----------------------------- 1961 - Updated G4Mag_SpinEqRhs.hh and G4Mag_SpinEq 1729 - Updated G4Mag_SpinEqRhs.hh and G4Mag_SpinEqRhs.cc for correction 1962 to spin tracking. 1730 to spin tracking. 1963 1731 1964 Apr 04th, 2001 G. Cosmo - field-V03-00-05 1732 Apr 04th, 2001 G. Cosmo - field-V03-00-05 1965 ------------------------ 1733 ------------------------ 1966 - G4FieldTrack.icc: fixed pedantic warnings o 1734 - G4FieldTrack.icc: fixed pedantic warnings on Linux-g++. 1967 1735 1968 Mar 23rd, 2001 J. Apostolakis - field-V03-00 1736 Mar 23rd, 2001 J. Apostolakis - field-V03-00-04 1969 ------------------------------ 1737 ------------------------------ 1970 - Electric Field unit test now works: test ru 1738 - Electric Field unit test now works: test runs until completion. 1971 Test outputs saved (on Linux). 1739 Test outputs saved (on Linux). 1972 1740 1973 - Fixed src/G4UniformElectricField.cc to set 1741 - Fixed src/G4UniformElectricField.cc to set correct field components 1974 (compat 1742 (compatible with Equation of Motion 1975 G4EqMa 1743 G4EqMagElectricField ) 1976 1744 1977 - Tidied G4FieldTrack 1745 - Tidied G4FieldTrack 1978 * deleted obsolete methods Set/Get Momentu 1746 * deleted obsolete methods Set/Get MomentumModulus() and its data member. 1979 * in constructor to set the value of the s << 1747 * in constructor to set the value of the spin to (0.0, 0.0, 0.0) 1980 if pSpin is null. (It was undefine 1748 if pSpin is null. (It was undefined.) 1981 1749 1982 Mar 23rd, 2001 J. Apostolakis - field-V03-00 1750 Mar 23rd, 2001 J. Apostolakis - field-V03-00-03 1983 ------------------------------ 1751 ------------------------------ 1984 - Fixed (hacked) CashKarp to work using Auxil << 1752 - Fixed (hacked) CashKarp to work using Auxiliary Stepper. 1985 1753 1986 - Modified G4MagIntegratorStepper << 1754 - Modified G4MagIntegratorStepper 1987 to use pointer to G4EquationOfMotion 1755 to use pointer to G4EquationOfMotion 1988 instead of G4Mag_EqRhs 1756 instead of G4Mag_EqRhs 1989 1757 1990 - Design iteration appears neccessary 1758 - Design iteration appears neccessary 1991 1759 1992 Mar 23rd, 2001 J. Apostolakis - field-V03-00 1760 Mar 23rd, 2001 J. Apostolakis - field-V03-00-02 1993 ------------------------------ 1761 ------------------------------ 1994 - Fixed G4MagHelicalStepper::AdvanceHelix for 1762 - Fixed G4MagHelicalStepper::AdvanceHelix for momentum, allowing 1995 the possibility of (correct) use of all H << 1763 the possibility of (correct) use of all Helical Steppers. 1996 1764 1997 - test/testProElectroMagField.cc updated for 1765 - test/testProElectroMagField.cc updated for momentum. 1998 (Not yet updated: Spin tests) 1766 (Not yet updated: Spin tests) 1999 1767 2000 1768 2001 Mar 19th, 2001 J. Apostolakis - field-V03-00 1769 Mar 19th, 2001 J. Apostolakis - field-V03-00-01 2002 ------------------------------ 1770 ------------------------------ 2003 - Tagging the modification of several classes 1771 - Tagging the modification of several classes to integrate in momentum*c_light 2004 instead of velocity. << 1772 instead of velocity. 2005 Reason for change: Velocity integration ca 1773 Reason for change: Velocity integration can result in superluminar 2006 velocities, due to arithmetic inaccuraci 1774 velocities, due to arithmetic inaccuracies. 2007 Classes modified: << 1775 Classes modified: 2008 G4FieldTrack, (Changed components 3,4,5 1776 G4FieldTrack, (Changed components 3,4,5 of SixVector to Momentum, 2009 G4Mag_EqRhs, G4Mag_UsualEqRhs, G4EqMagEle 1777 G4Mag_EqRhs, G4Mag_UsualEqRhs, G4EqMagElectricField 2010 1778 2011 Feb 20th, 2001 J. Apostolakis << 1779 Feb 20th, 2001 J. Apostolakis 2012 1780 2013 - Modified several classes to integrate in mo 1781 - Modified several classes to integrate in momentum instead of velocity. 2014 1782 2015 G4FieldTrack: 1783 G4FieldTrack: 2016 Changed components 3,4,5 of SixVector to 1784 Changed components 3,4,5 of SixVector to Momentum, 2017 Modified constructor 1785 Modified constructor 2018 Added energy calculation to "LoadFromArr 1786 Added energy calculation to "LoadFromArray" method 2019 GetEnergy method is now implemented. << 1787 GetEnergy method is now implemented. 2020 Deleted obsolete (old) method names for 1788 Deleted obsolete (old) method names for Curve and Position. New have "Get" 2021 1789 2022 G4ChordFinder 1790 G4ChordFinder 2023 Renamed old-method names of Field Track << 1791 Renamed old-method names of Field Track to new ones. 2024 1792 2025 Jan 29th, 2001 G. Cosmo - field-V03-00-00 1793 Jan 29th, 2001 G. Cosmo - field-V03-00-00 2026 ------------------------- 1794 ------------------------- 2027 1795 2028 - Minor fix to G4MagneticField.hh: added expl 1796 - Minor fix to G4MagneticField.hh: added explicit call to G4Field() in 2029 copy constructor. Fixes report #205. 1797 copy constructor. Fixes report #205. 2030 1798 2031 Nov 20th, 2000 G. Cosmo - field-V02-00-02 1799 Nov 20th, 2000 G. Cosmo - field-V02-00-02 2032 1800 2033 - Fixes to remove warnings from "-Wall -ansi 1801 - Fixes to remove warnings from "-Wall -ansi -pedantic" g++ compiler options: 2034 o commented out variables declared and not 1802 o commented out variables declared and not used. 2035 o fixed declaration of variables used befor 1803 o fixed declaration of variables used before being initialised. 2036 o fixed order of initialisation of member d 1804 o fixed order of initialisation of member data in constructors. 2037 o fixed usage of unsigned-int (size_t) for 1805 o fixed usage of unsigned-int (size_t) for array indeces. 2038 1806 2039 Nov 9th, 2000 G. Cosmo - field-V02-00-01 1807 Nov 9th, 2000 G. Cosmo - field-V02-00-01 2040 1808 2041 - G4FieldManager: added check on existence of 1809 - G4FieldManager: added check on existence of allocated fChordFinder 2042 in CreateChordFinder(G4MagneticField*) meth 1810 in CreateChordFinder(G4MagneticField*) method. 2043 1811 2044 Nov 1st, 2000 G. Cosmo - field-V02-00-00 1812 Nov 1st, 2000 G. Cosmo - field-V02-00-00 2045 1813 2046 - QA code revision and cleanup (+fixes from C 1814 - QA code revision and cleanup (+fixes from CodeWizard filtering): 2047 o Added (private) declarations of copy cons 1815 o Added (private) declarations of copy constructor and assignment operator 2048 where needed. 1816 where needed. 2049 o Added "const" qualifier to accessor metho 1817 o Added "const" qualifier to accessor methods, wherever needed. 2050 o Added equality check on operator= where n 1818 o Added equality check on operator= where needed. 2051 o Made destructor "virtual" for classes hav 1819 o Made destructor "virtual" for classes having virtual methods, 2052 wherever missing. 1820 wherever missing. 2053 o Changed virtual functions to non-inline w 1821 o Changed virtual functions to non-inline wherever's the case. 2054 o Changed return value to be "const" pointe 1822 o Changed return value to be "const" pointer to functions: 2055 - G4EquationOfMotion::GetFieldObj() 1823 - G4EquationOfMotion::GetFieldObj() 2056 - G4FieldManager::GetDetectorField() 1824 - G4FieldManager::GetDetectorField() 2057 - G4MagIntegratorDriver::GetStepper() 1825 - G4MagIntegratorDriver::GetStepper() 2058 o Fixed and corrected signature to G4MagInt 1826 o Fixed and corrected signature to G4MagIntegratorDriver::SetHmin(). 2059 o Added "inline" qualifier to methods' decl 1827 o Added "inline" qualifier to methods' declarations. 2060 o General cosmetics 1828 o General cosmetics 2061 1829 2062 Jun 2nd, 2000 J. Apostolakis 1830 Jun 2nd, 2000 J. Apostolakis 2063 - G4FieldManager[.hh.icc] 1831 - G4FieldManager[.hh.icc] 2064 o Added DeltaIntersection and DeltaOneStep 1832 o Added DeltaIntersection and DeltaOneStep to FieldManager: added 2065 Get/Set methods and SetAccuraciesWithDelt 1833 Get/Set methods and SetAccuraciesWithDeltaOneStep(). 2066 1834 2067 May 31st, 2000 J. Apostolakis - field-V01-01 1835 May 31st, 2000 J. Apostolakis - field-V01-01-03 2068 1836 2069 include/G4FieldManager.hh,icc 1837 include/G4FieldManager.hh,icc 2070 o Moved Delta Intersection and Delta One St 1838 o Moved Delta Intersection and Delta One Step to FieldManager 2071 (from Propagator in Field), to allow th 1839 (from Propagator in Field), to allow them to differ over the detector. 2072 1840 2073 May 11th, 2000 J. Apostolakis - field-V01-01 1841 May 11th, 2000 J. Apostolakis - field-V01-01-02 2074 1842 2075 - G4ChordFinder[.hh.cc] 1843 - G4ChordFinder[.hh.cc] 2076 o Modified method for finding new trial ste 1844 o Modified method for finding new trial step. 2077 New method assumes a quadratic relationsh 1845 New method assumes a quadratic relationship between step-length and d_chord: 2078 d_chord proportional to step-l 1846 d_chord proportional to step-lenght ^ 2 2079 o Added two small refinement in order to av 1847 o Added two small refinement in order to avoid unneccesary extra step due to 2080 small numerical inaccuracies in estimatio 1848 small numerical inaccuracies in estimation: 2081 FindNextChord 1849 FindNextChord 2082 to first step use (1-0.001) 1850 to first step use (1-0.001) * last-step-estimate 2083 New Step: 1851 New Step: 2084 to subsequent step use trial_st 1852 to subsequent step use trial_step = 0.98 * estimate 2085 1853 2086 May 10th, 2000 G. Cosmo 1854 May 10th, 2000 G. Cosmo 2087 1855 2088 - G4MagneticField.hh: minor fix to dummy impl 1856 - G4MagneticField.hh: minor fix to dummy implementation of assignment 2089 operator. 1857 operator. 2090 1858 2091 May 9th, 2000 J. Apostolakis - field-V01 1859 May 9th, 2000 J. Apostolakis - field-V01-01-01 2092 1860 2093 - G4MagErrorStepper: << 1861 - G4MagErrorStepper: 2094 Modified DistChord to check whether the c 1862 Modified DistChord to check whether the chord endpoint are the same point. 2095 1863 2096 - G4RKG3_Stepper: 1864 - G4RKG3_Stepper: 2097 Deleted error message from StepNoErr, as 1865 Deleted error message from StepNoErr, as it is used in Stepper. 2098 1866 2099 - G4MagIntegratorDriver: 1867 - G4MagIntegratorDriver: 2100 Changed the condition for the debug outpu 1868 Changed the condition for the debug output (~cosmetic) 2101 1869 2102 1870 2103 April 27th, 2000 G. Cosmo 1871 April 27th, 2000 G. Cosmo 2104 1872 2105 - Reorganised comments for the Software Refer 1873 - Reorganised comments for the Software Reference Manual. 2106 - Moved inlined definition to .icc files, whe 1874 - Moved inlined definition to .icc files, where needed. 2107 Created: G4ChordFinder.icc and G4FieldManag 1875 Created: G4ChordFinder.icc and G4FieldManager.icc. 2108 1876 2109 April 13th, 2000 J. Apostolakis - field-V01 1877 April 13th, 2000 J. Apostolakis - field-V01-01-00 2110 1878 2111 - Modified signature to function DumbStepper( 1879 - Modified signature to function DumbStepper() in classes: 2112 G4HelixExplicitEuler, G4HelixHeum, G4HelixI 1880 G4HelixExplicitEuler, G4HelixHeum, G4HelixImplicitEuler, 2113 G4HelixSimpleRunge and G4MagHelicalStepper. 1881 G4HelixSimpleRunge and G4MagHelicalStepper. 2114 - Modified signature and implementation of fu 1882 - Modified signature and implementation of function MagFieldEvaluate() 2115 in G4MagHelicalStepper, to use Field (as Th 1883 in G4MagHelicalStepper, to use Field (as ThreeVector) in place of dydx. 2116 - Made Richardson extrapolation optional (it 1884 - Made Richardson extrapolation optional (it is unclear if 2117 it offers benefit) in G4MagHelicalStepper. 1885 it offers benefit) in G4MagHelicalStepper. 2118 - Updated unit test testPropagateMagField and 1886 - Updated unit test testPropagateMagField and reference output. 2119 - G4LineSection: fixes to prevent divisions b 1887 - G4LineSection: fixes to prevent divisions by 0. and sqrt(<0) 2120 (V.Grichine). 1888 (V.Grichine). 2121 1889 2122 January 19th, 2000 G. Cosmo - field-V01-00 1890 January 19th, 2000 G. Cosmo - field-V01-00-00 2123 1891 2124 Syncronized versioning of files with HEAD aft 1892 Syncronized versioning of files with HEAD after ISO-C++ migration 2125 (essentially CVS header changes). 1893 (essentially CVS header changes). 2126 1894 2127 December 7th, 1999 J. Apostolakis - field- 1895 December 7th, 1999 J. Apostolakis - field-V00-01-01a 2128 1896 2129 For negative or very small curve lengths (h) 1897 For negative or very small curve lengths (h) we suppress writing an error 2130 (about the distance of endpoints being furthe 1898 (about the distance of endpoints being further than the curve length), 2131 as it is meaningless. 1899 as it is meaningless. 2132 1900 2133 November 23rd, 1999 J. Apostolakis - field- 1901 November 23rd, 1999 J. Apostolakis - field-V00-01-01 2134 1902 2135 Tagged this sub-category's state in geometry- 1903 Tagged this sub-category's state in geometry-V00-01-01 2136 1904 2137 1905 2138 July 6th, 1999 J. Apostolakis - field- 1906 July 6th, 1999 J. Apostolakis - field-01-00-07 2139 1907 2140 These fixes and those in field-01-00-06 and f 1908 These fixes and those in field-01-00-06 and field-01-00-06a all come 2141 from the need to cope with the imprecision of 1909 from the need to cope with the imprecision of the integration, and its 2142 interaction with paths that are rather straig << 1910 interaction with paths that are rather straight. 2143 1911 2144 G4MagIntegrationDriver: << 1912 G4MagIntegrationDriver: 2145 i) AccurateAdvance now that the chord distan 1913 i) AccurateAdvance now that the chord distance is smaller than the curve len- 2146 gth (times 1+epsilon). 1914 gth (times 1+epsilon). 2147 1915 2148 ii) OneGoodStep and QuickAdvance considers th 1916 ii) OneGoodStep and QuickAdvance considers the error of the velocity as well as 2149 the error of the position - and ensure that b 1917 the error of the position - and ensure that both relative errors are within 2150 the "eps" maximum given. 1918 the "eps" maximum given. 2151 1919 2152 G4ChordFinder: 1920 G4ChordFinder: 2153 Having found that the curve_length vs ABdis 1921 Having found that the curve_length vs ABdist disrepancy has its basis in 2154 the imprecision of the integration: 1922 the imprecision of the integration: 2155 i) Modified the exception in this case to h 1923 i) Modified the exception in this case to happen only for very large 2156 relative disrepancies ( 10 * eps ). 1924 relative disrepancies ( 10 * eps ). 2157 ii) Undertaken a default corrective action i 1925 ii) Undertaken a default corrective action in all other cases 2158 Note that this disrepancy only happens for ve 1926 Note that this disrepancy only happens for very straight sections of track - 2159 ie step lengths much smaller than the curvatu 1927 ie step lengths much smaller than the curvature of the track. 2160 1928 2161 July 1st, 1999 J. Apostolakis - field- 1929 July 1st, 1999 J. Apostolakis - field-01-00-06a 2162 1930 2163 In order to cope with the anomalous condition 1931 In order to cope with the anomalous condition created by candidate fix 2164 in G4PropagatorInField, we ensure that the f 1932 in G4PropagatorInField, we ensure that the fraction_AE is always between 2165 0 and 1.0. (giving it a default value of 0.5 1933 0 and 1.0. (giving it a default value of 0.5 whenever it is wrong). 2166 ( in G4ChordFinder::ApproxCurvePointV ) 1934 ( in G4ChordFinder::ApproxCurvePointV ) 2167 1935 2168 June 29th, 1999 J. Apostolakis - field- 1936 June 29th, 1999 J. Apostolakis - field-01-00-06 2169 1937 2170 I re-instated a check on whether the input 1938 I re-instated a check on whether the input distance along the curve between 2171 the two input curve points is truly greater t 1939 the two input curve points is truly greater than their linear distance. 2172 ( in G4ChordFinder::ApproxCurvePointV ) 1940 ( in G4ChordFinder::ApproxCurvePointV ) 2173 1941 2174 The condition that prints out a warning messa 1942 The condition that prints out a warning message has been altered to 2175 take into account the case of a very small fi << 1943 take into account the case of a very small final interval. 2176 ( in G4MagInt_Driver::AccurateAdvance ) 1944 ( in G4MagInt_Driver::AccurateAdvance ) 2177 1945 2178 June ??th, 1999 J. Apostolakis - field- 1946 June ??th, 1999 J. Apostolakis - field-01-00-05 2179 1947 2180 The condition that prints out a warning messa << 1948 The condition that prints out a warning message was changed to delete 2181 an incorrect comparison. << 1949 an incorrect comparison. 2182 ( in G4MagInt_Driver::AccurateAdvance ) 1950 ( in G4MagInt_Driver::AccurateAdvance ) 2183 1951 2184 April 19th, 1999 J. Apostolakis - field- 1952 April 19th, 1999 J. Apostolakis - field-01-00-04 2185 1953 2186 G4MagErrorStepper: Cosmetic changes to name o 1954 G4MagErrorStepper: Cosmetic changes to name of data members 2187 1955 2188 March 4th, 1999 J. Apostolakis - field- 1956 March 4th, 1999 J. Apostolakis - field-01-00-03 2189 1957 2190 Corrected delete statement of arrays to use ' 1958 Corrected delete statement of arrays to use 'delete[]' instead of simple delete 2191 in destructors of 2 steppers (G4CashKarpRKF45 1959 in destructors of 2 steppers (G4CashKarpRKF45, G4SimpleHeum) 2192 1960 2193 February 17th, 1999 J. Apostolakis - field 1961 February 17th, 1999 J. Apostolakis - field-01-00-02 2194 1962 2195 Moved G4PropagatorInField files to geometry/ << 1963 Moved G4PropagatorInField files to geometry/volumes. 2196 Updated GNUmakefile not to depend on geometr 1964 Updated GNUmakefile not to depend on geometry/volumes anymore. 2197 -> the known circular dependency is remove 1965 -> the known circular dependency is removed 2198 1966 2199 1967 2200 February 12th, 1999 J. Apostolakis - field 1968 February 12th, 1999 J. Apostolakis - field-spin1-works1 (act Feb17) 2201 1969 2202 A tag that contains the first version that 1970 A tag that contains the first version that worked with spin (in Peter G.'s 2203 directory). G4PropagateInSpin is still in t 1971 directory). G4PropagateInSpin is still in this version. A small change 2204 in a branch of the header file was needed. << 1972 in a branch of the header file was needed. 2205 1973 2206 1974 2207 February 10th, 1999 J. Apostolakis << 1975 February 10th, 1999 J. Apostolakis 2208 1976 2209 All the following changes are needed for ener 1977 All the following changes are needed for energy integration also: 2210 1978 2211 G4ChordFinder now passes to the constructed d 1979 G4ChordFinder now passes to the constructed driver the number of variables. 2212 Because of the need for the stepper to know a 1980 Because of the need for the stepper to know about it in its Abstract interface 2213 1981 2214 -> Moved theNumberOfVariables from derived c 1982 -> Moved theNumberOfVariables from derived class G4MagErrorStepper 2215 to base class 1983 to base class G4MagIntegratorStepper 2216 1984 2217 As a consequense the classes that derived fro 1985 As a consequense the classes that derived from G4MagIntegratorStepper directly 2218 were changed: G4CashKarpRKF45, G4MagHelicalSt << 1986 were changed: G4CashKarpRKF45, G4MagHelicalStepper, 2219 1987 2220 1988 2221 February 9th, 1999 J. Apostolakis 1989 February 9th, 1999 J. Apostolakis 2222 1990 2223 Equation of Motion: 1991 Equation of Motion: 2224 G4MagEqRhs.hh Made SetChargeMomentumMass vir 1992 G4MagEqRhs.hh Made SetChargeMomentumMass virtual, so that G4Mag_SpinEqRhs 2225 can initialise its constants << 1993 can initialise its constants too. 2226 1994 2227 1995 2228 February 9th, 1999 P. Gumplinger 1996 February 9th, 1999 P. Gumplinger 2229 1997 2230 Added Equation of Motion for Spin in Magnetic 1998 Added Equation of Motion for Spin in Magnetic Field. 2231 1999 2232 2000 2233 February 17th, 1999 J. Apostolakis - fiel 2001 February 17th, 1999 J. Apostolakis - field-01-00-01 2234 ---> date is correct, made sure that tag d 2002 ---> date is correct, made sure that tag does not include above fixes! 2235 << 2003 2236 Moved G4PropagatorInField files to geometry/ << 2004 Moved G4PropagatorInField files to geometry/volumes. 2237 Updated GNUmakefile not to depend on geometr 2005 Updated GNUmakefile not to depend on geometry/volumes anymore. 2238 -> the known circular dependency is remove 2006 -> the known circular dependency is removed 2239 2007 2240 2008 2241 --------------------------------------------- 2009 --------------------------------------------------------------------------- 2242 ------ Geant 4 Production Release - 2010 ------ Geant 4 Production Release ------------ 2243 --------------------------------------------- 2011 -------------------------------------------------------------------------- 2244 Nov 19th, 1998 J. Apostolakis - field-0 2012 Nov 19th, 1998 J. Apostolakis - field-00-04-01head ===> R+D version only 2245 2013 2246 G4PropagatorInField files modified: 2014 G4PropagatorInField files modified: 2247 2015 2248 - to store the last safety computed and it 2016 - to store the last safety computed and its origin. 2249 - to give them to a requester 2017 - to give them to a requester 2250 --------------------------------------------- 2018 ------------------------------------------------------------------ 2251 2019 2252 Nov 19th, 1998 J. Apostolakis - field-0 2020 Nov 19th, 1998 J. Apostolakis - field-00-03-03b 2253 2021 2254 Subtle fix to interface of G4UniformMagField << 2022 Subtle fix to interface of G4UniformMagField::GetFieldValue 2255 ( that tripped Visual C++ ) 2023 ( that tripped Visual C++ ) 2256 2024 2257 G4PropagatorInField: added printing method 2025 G4PropagatorInField: added printing method for verboseness. 2258 2026 2259 Nov 19th, 1998 J. Apostolakis - field-0 2027 Nov 19th, 1998 J. Apostolakis - field-00-03-03a 2260 2028 2261 Steppers: Fixed instance variables constr 2029 Steppers: Fixed instance variables construction (SimpleRunge) 2262 and clash of names in CashKarp: 2030 and clash of names in CashKarp::StepWithEstimate 2263 2031 2264 Added post-const to all Field classes for 2032 Added post-const to all Field classes for 2265 void GetFieldValue(const G4double 2033 void GetFieldValue(const G4double yTrack[] , 2266 G4double B[] 2034 G4double B[] ) const ; 2267 2035 2268 Added SetFieldValue methods to UniformMagF 2036 Added SetFieldValue methods to UniformMagField 2269 2037 2270 Added GetConstantFieldValue() methods to U 2038 Added GetConstantFieldValue() methods to Uniform Mag & Electric Fields 2271 2039 2272 2040 2273 Nov 19th, 1998 J. Apostolakis - field-0 2041 Nov 19th, 1998 J. Apostolakis - field-00-03-03 2274 2042 2275 Fixes to Propagate and Steppers 2043 Fixes to Propagate and Steppers 2276 2044 2277 Nov 19th, 1998 J. Apostolakis - field-0 2045 Nov 19th, 1998 J. Apostolakis - field-00-03-02 2278 2046 2279 Minor corrections (???) 2047 Minor corrections (???) 2280 2048 2281 Nov 19th, 1998 J. Apostolakis - field-0 2049 Nov 19th, 1998 J. Apostolakis - field-00-03-02 2282 << 2050 2283 Major development to allow Electric Field 2051 Major development to allow Electric Field 2284 << 2052 2285 2053 2286 2054 2287 June 17th, 1998 J. Apostolakis - magfiel 2055 June 17th, 1998 J. Apostolakis - magfield-06-03 2288 2056 2289 Modified Files: 2057 Modified Files: 2290 include/G4PropagatorInField.hh 2058 include/G4PropagatorInField.hh 2291 src/G4PropagatorInField.cc 2059 src/G4PropagatorInField.cc 2292 test/testPropagateInField.cc 2060 test/testPropagateInField.cc 2293 --------------------------------------------- 2061 ---------------------------------------------------------------------- 2294 2062 2295 1) Added an assertion after each call to Lo << 2063 1) Added an assertion after each call to LocateGlobalPointAndSetup. 2296 It checks that the located volume is equ 2064 It checks that the located volume is equal to the current volume. 2297 If this check fails, the Propagator does 2065 If this check fails, the Propagator does not know how to continue. 2298 2066 2299 2) Added an argument to ComputeStep method 2067 2) Added an argument to ComputeStep method of G4PropagatorInField. 2300 The new volume argument is used to check th 2068 The new volume argument is used to check the current volume during the step. 2301 2069 2302 3) When G4Navigator's ComputeStep returned 2070 3) When G4Navigator's ComputeStep returned stepLength==stepRequested, this 2303 used to be treated as a failure to intersect. 2071 used to be treated as a failure to intersect. This is incorrect. 2304 2072 2305 stepLength==stepRequested is an intersect 2073 stepLength==stepRequested is an intersection, 2306 so I am now treating it as such in each 2074 so I am now treating it as such in each comparison (of 3). 2307 Ignoring it causes problems, including 2075 Ignoring it causes problems, including volume being skipped. 2308 --> but treating it as an intersection 2076 --> but treating it as an intersection used to cause problems too ?? 2309 2077 2310 2078 2311 May 7th, 1998 J. Apostolakis - magfie 2079 May 7th, 1998 J. Apostolakis - magfield-06-02-helix01 2312 2080 2313 - Tag includes new version of G4ChordFinder ( 2081 - Tag includes new version of G4ChordFinder (not used by W.Wander). 2314 2082 2315 2083 2316 May 7th, 1998 W. Wander (recorded/c 2084 May 7th, 1998 W. Wander (recorded/commited by J. Apostolakis) 2317 2085 2318 - New "helical" RK steppers that move in heli 2086 - New "helical" RK steppers that move in helical segments, not linear. 2319 * Several changes to add steppers << 2087 * Several changes to add steppers 2320 2088 2321 Modified Files: 2089 Modified Files: 2322 G4MagErrorStepper.hh << 2090 G4MagErrorStepper.hh 2323 G4MagErrorStepper.cc << 2091 G4MagErrorStepper.cc 2324 added Helical stepping methods to 2092 added Helical stepping methods to this general ABC. 2325 << 2093 2326 G4MagIntegratorStepper.hh 2094 G4MagIntegratorStepper.hh 2327 made RightHandSide virtual (this 2095 made RightHandSide virtual (this will be reversed very soon) 2328 2096 2329 testPropagateMagField.cc 2097 testPropagateMagField.cc 2330 Test of Helical stepping functiona 2098 Test of Helical stepping functionality added. 2331 << 2099 2332 Added Files for new Helical RK steppers, der 2100 Added Files for new Helical RK steppers, deriving from G4MagErrorStepper: 2333 G4HelixExplicitEuler.cc ( and .hh ) 2101 G4HelixExplicitEuler.cc ( and .hh ) 2334 G4HelixHeum.cc 2102 G4HelixHeum.cc 2335 G4HelixImplicitEuler.cc 2103 G4HelixImplicitEuler.cc 2336 G4HelixSimpleRunge.cc << 2104 G4HelixSimpleRunge.cc 2337 2105 2338 Added File for implementation: 2106 Added File for implementation: 2339 G4MagErrorStepper.icc 2107 G4MagErrorStepper.icc 2340 2108 2341 2109 2342 May 7th, 1998 J. Apostolakis - magfi 2110 May 7th, 1998 J. Apostolakis - magfield-06-01 2343 2111 2344 - added small changes to << 2112 - added small changes to 2345 * G4ChordFinder: methods to get 2113 * G4ChordFinder: methods to get/set IntegratorDriver 2346 * G4MagIntegratorDriver: made OneGoodSt 2114 * G4MagIntegratorDriver: made OneGoodStep method public 2347 - created the file. 2115 - created the file.