I just found a bug in the production Fast MC clustering causing the
position of the cluster to be wrongly reported. (Theta is incorrectly
calculated from the rectangular coordinates of the cluster so that it
never has a value between pi/4 and 3*pi/4.)
This is fixed in our development version. If you are planning to use
the released version, you can patch it as follows:
1. Find the file CalSmeared.cxx. The full specification should end with
RootApps/FastMC/V1.0/src/CalSmeared.cxx
2. It should contain lines looking like this:
theta = TMath::ATan2(toBeamline, r);
if (z < 0.0) theta = TMath::Pi() - theta;
Change them to
theta = TMath::ATan2(toBeamline, z);
if (theta < 0.0) theta += TMath::Pi();
3. cd to the bld subdirectory and rebuild the Fast MC executable,
e.g., assuming you're in the directory containing CalSmeared.cxx
$ cd ../bld
$ gmake
Joanne
--
[log in to unmask] SLAC, MS 71
(650) 926-3685 Stanford University
P.O. Box 4349
Stanford CA, 94309
|