Message: Re: Different radii for the 'same' electron in B-field Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Different radii for the 'same' electron in B-field 

Keywords: GEATN4 magnetic field electrons curve
Forum: Fields: Magnetic and Otherwise
Re: None Different radii for the 'same' electron in B-field (Valentin)
Date: 10 May, 2005
From: John Apostolakis <John Apostolakis>

I think that you are not extracting the points on the trajectory correctly. The points at which GetFieldValue is called are not on the trajectory of the particle. Following the RungeKutta method, the field is sampled at different trial points, coming from linear and higher order approximations of the final and intermediate points of a step. Then these are used to estimate the endpoint of a step, which is accepted if its estimate is deemed accurate enough (see below).

Of course if you used a method that assumed that the field was uniform, then it would only evaluate it at the step's starting point. But since your real problem has a non-uniform field it makes sense to use the general Runge-Kutta based methods.

To get the actual points of the trajectory, you should create and use a G4Trajectory of the endpoints of the steps or another way to store or print only the endpoints of each step. To find out how to create a trajectory please see section 5.1.6 of the Users Guide for Application Developers (UGAD), which you can find at http://cern.ch/geant4/G4UsersDocuments/UsersGuides/ForApplicationDeveloper/html/TrackingAndPhysics/tracking.html

On a related topic, for more information about obtaining the accuracy you require for your application, please see the UGAD, section 4.3 and in particular the part of 4.3.2 on "How to Adjust the Accuracy of Propagation" http://cern.ch/geant4/G4UsersDocuments/UsersGuides/ForApplicationDeveloper/html/Detector/electroMagneticField.html

In particular if you require high accuracy for the endpoints, you should choose small values for the epsilon parameters. In general these should be smaller than the ratio of the error in the radius (or endpoint) that you are willing to tolerate over the total distance travelled by the track (along the helix). I would expect good results with values the order 10^-5, and do not know fields which are measured better than about one part in 10^4 or so anyway. I note that values below 10^-9 or so will be very challenging to achieve and likely not meaningful due to precision issues.

Best regards, John Apostolakis

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: Different radii for the 'same' electron in B-field   (Valentin - 16 May, 2005)
 Add Message Add Message
to: "Re: Different radii for the 'same' electron in B-field"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]