Muon swimming code

This code is courtesy of Darin Acosta ( acosta@phys.ufl.edu).

This is a primitive track swimming program based on a very old field map and detector geometry. It was used to calculate the number of strips traversed by a muon crossing a station. But as it is a very tiny program, one could use it for many things. For example, it also calculates the phi position at each station, so you could study the track-finding parameters. Anyway, I've been meaning to update it with the latest detector positions and sizes, latest B field (haven't heard from Dick about this yet), and code for some histograms related to track finding (there are no histograms yet). So I'm not sure if it is ready for public prime time, but I attach the relevant files anyway for you and Jay. (If the attachments aren't readable, let me know and I'll send the files one at a time).

To use it, compile the Fortran. No other libraries needed. Then make the link: 'ln -s cms13hs.fieldmap fort.25' You set theta, phi, and Pt of the track in the data statements.

field_all.f
cms13hs.fieldmap

On 9 July Darin sent me an update version. Here's some information from the email he sent me:

It now uses the latest field map from Wisconsin. Also, I updated the chamber positions and strip widths, and send the output to an Ntuple. I've been using it to study the bending in phi for the track finder. Sometime soon I may have it swim tracks in R also to study the bending in the barrel, but that is on hold for the moment. Anyway, the output is fairly primitive. Station 0 is Me1/1, 1 through 4 are Me1/2(3) through ME4. The angles are in degrees, the bending is in strips.

To use the new field map, I started from Loveless's data file rather than a TZ file generated by UCDavis. It seemed easier to do that way. There is no smoothing of the field or interpolation, however. I reformatted the data into my own data file which I call 'brz100.dat'. Basically you just link fort.25 to it to run my new program 'field.f'.

field.f
brz100.dat.gz

last updated 27 July 1998 by Benn Tannenbaum (benn@physics.ucla.edu)