|Message: Times or positions don't seem correct in synchrotron production in atmosphere model||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
I have been trying to figure this out for months to no avail. It's time to consult the experts.
I have a model of the atmosphere using G4Orbs embedded like "Russian dolls" with equal depth (20 mg/cm2) between layers. I have a model of the Earth's magnetic field, and set myField = new CrestMagField(latitude, longitude, time); G4FieldManager* fieldMgr = G4TransportationManager::GetTransportationManager() ->GetFieldManager(); fieldMgr->SetDetectorField(myField); fieldMgr->CreateChordFinder(myField); after constructing all of my geometry. Field strengths are ~0.5 gauss and electron energies are ~10 TeV. A more detailed description of my geometry is in a recent posting in the Geometry section.
I record all synchrotron photon hit locations and times at a certain height. I should be able to look at two close together (within meters) photons and reconstruct their zenith angle using their time and position differences. However, the resulting zenith angle is always nearly vertical, meaning either the times are too close together or the positions are too far apart. I carefully tag photons when produced by synchrotron and not undergoing any other physics, such as Compton scattering, before getting to my recording height. So unless the Transportation process itself can change directions, these should have traveled in a straight line. I verified the photns and the primary electron travel at c_light.
The zenith angle reconstruction method works fine with an old "flat earth" model based on the old GEANT3 Runge-Kutta tracker, so that works fine.
The default values for the Chordfinder are
fieldMgr->GetDeltaIntersection(): 1 um
fieldMgr->GetDeltaOneStep(): 10 um
End of run statistics for one track for the chordfinder are G4ChordFinder statistics report: No trials: 4842 No Calls: 4701 Max-trial: 2 Parameters: fFirstFraction 0.999 fFractionLast 1 fFractionNextEstimate 0.98
Any ideas? I don't know where else to pursue this.