Rick Van Kooten
Indiana University High Energy Physics
For Neal and co.
Code and examples using the amplitude method as described in Moser
et al, NIMA A384 (1997) 491.
The following code runs from within
a PAW session as kumacs and COMIS FORTRAN files. e.g. copy all the *.f
files and *.kumac files into your directory and run PAW. Start by
typing "exec bsmix" to get it going. These have been ripped out of
other code and tried to simplify to bare bones.
-
bsmix.f
Toy MC that generates distributions to test the amplitude fit on.
Uses typical dilution factors and proper time resolution factors. Fills
Histo 800 with total time distribution, 700 with mixed distribution, 600
with unmixed distribution.
-
bscalc.f
Calculates the relevant asymmetry as a function of proper time, including
proper errors. Stores answer and errors in vectors.
-
bsfun.f
Fitting function for fitting to the time-dependent asymmetry.
-
bsans.f
Finds x_s bin with largest amplitude result - that is the "answer" of the
fit.
-
bsmix.kumac
Driver kumac. Scans across different values of x_s, fitting amplitude
every time. Calls the above.
-
bsmix_loop.kumac
Start at scanning different _generated_ x_s and looking as average significance
of fitted result - not finished yet
Toy Results
(as postscript files)
-
bslife.ps
Check: total proper time distribution, smeared and unsmeared with exponential
fit.
-
Generated x_s = 20
-
Unsmeared , bs_xs20_hunt.ps,
goofy plot, scanning across the "frequency" x_s
-
Unsmeared, bs_xs20_fit.ps,
fit at best x_s giving largest amplitude
-
Unsmeared, bs_xs20_amp.ps,
fitted amplitude as function of x_s showing peak at generated x_s
-
Smeared,
bs_xs20_fit_smear.ps,
fit at best x_s giving largest amplitude
-
Smeared,
bs_xs20_amp_smear.ps,
fitted amplitude as function of x_s showing peak at generated x_s
-
Generated x_s = 16
-
Generated x_s = 12
This file is http://needmore.physics.indiana.edu/~rickv/bsosc/bsosc.html.
Last updated: 14 Sept. 1999
Rick
Van Kooten , rickv@paoli.physics.indiana.edu