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