Re: MGDRAW Directional Cosines

From: John Clem <clem_at_bartol.udel.edu>
Date: Mon, 23 Jan 2012 14:44:00 -0500

Vasilis et al., the ultimate goal is to use FLUKA to characterize
synchrotron and ck radiation in the radio frequency regime for EM energy
showers in different materials and magnetic fields. It was my hope the
necessary kinematic parameters could be extracted from MGDRAW as input
to a RF Yield function (the Lienard–Wiechart potentials) such that each
track segment could be transform to an E-field element at given location
relative to the shower. Summing all the track segment contributions
should produce the expected E-field pulse in the time domain of course.
The list of parameters needed to do this:

1) X, Y, Z values at each track segment end-point (XTRACK, YTRACK, ZTRACK)
2) Velocity vector at each track end-point
3) Particle Time/Age at each track end-point (ATRACK for each end-point)
4) Curved path between each track segment end-point (CTRCK for each segment)
5) Electric charge of the particle "ICHRGE(JTRACK)"

Is there a way to extract these values ?? I realize this can be done
within GEANT, but for comparison purposes, I would like to perform this
calculation in FLUKA if possible. Any help would be appreciated.

Thanks, John

On 1/12/2012 10:35 AM, Vasilis Vlachoudis wrote:
> Hi John,
>
> you could do something like this (of course it would require a loop for
> j=NTRACK to 1 by -1),
> but you will not get the desired incoming direction, but the one actually
> performed due to the MCS.
> Which I would assume would be almost-equal with first sub-step direction
> DX = Xtrack(1) - Xtrack(0)
> DY = Ytrack(1) - Ytrack(0)
> DZ = Ztrack(1) - Ztrack(0)
> and then get dx/sqrt(dx^2 + dy^2 + dz^2), ...
>
> The correct incoming direction you will have it from the previous C[XYZ]TRCK variables,
> which you have to save somewhere.
>
> Vasilis Vlachoudis
> Dep EN, CERN
> CH-1211 GENEVA 23
> SWITZERLAND
>
> Phone: +41-22 767 9851
> GSM: +41-76 487 4378
> Fax: +41-22 766 9644
>

-- 
>>  John Clem
>>  Research Associate Professor
>>  Bartol Research Institute
>>  Department of Physics and Astronomy
>>  217 Sharp Lab
>>  University of Delaware
>>  Newark, DE 19716
>>  off 302-831-4354
>>  cell 302-388-1218
Received on Tue Jan 24 2012 - 00:26:02 CET

This archive was generated by hypermail 2.2.0 : Tue Jan 24 2012 - 00:26:03 CET