- 1996.05.31
- add hooks for user routines:
- USR_INIT
- Called before any standard gminos initialization
actions. This is suitable for adding FFREAD cards.
- USR_HIST
- Book histograms in this routine
- USR_BEGRUN
- Called after standard gminos initialization
but before any events are generated
- USR_EVENT
- Called after event has been generated, tracked, digitized.
- USR_FINISH
- Called after final event
- 1996.06.19
- preliminary digitization scheme for FLS hits
- 1996.06.21
- for detector planes mix up a new "Material" with the right
proportions of elements based on materials used in extr/comb/cell
etc.; create new "TMedium" and list it as the relevant "TMedium"
in PlaneSpec_TMedium.
- 1996.06.28
- GAFF card controls output gaf format (0=IE,1=FZ);
GAFN card sets file name (use single quotes around
input string). If file name has the
extension '.xx_gaf' then "xx" is substituted with "ie" or "fz"
as appropriate.
- 1996.07.02
- FATL card controls waveshifting fiber attenuation lengths
(default short=135, long=1100cm),
FATN card sets relative normalizations (default 1.0, 0.66)
- 1996.07.03
- add hooks for different kinematics methods: card MKIN
picks between usr_kine (0), hepevt_file_kine (1),
and neugen_kine (2). The neugen integration is not
yet complete, so this is a dummy routine at present. Assorted
small corrections also entered into the repository.
- 1996.07.19
- lots of changes! new flux directory does not get retrieved using
"cvs update" (cd above gminos and issue "cvs checkout gminos";
this will update what is there and retreive the rest). "arch"
routine in flux directory needs a "chmod +x arch" for the Makefile
to work. Some of the Makefiles (and Make.include) will have
CVS conflicts -- if you haven't made local modifications it might
be easier to remove them and let "cvs update" get fresh copies.
- removed conflicts in particle codes. NEUGEN now outputs
STDHEP id's for ions. LBL_PART no longer conflicts with
GPIONS in the assignment of GEANT particle codes. explicit
neutrino flavor codes are defined for GEANT and various decay
modes are updated to reflect them. Fewer heavy particle states
get defined for GEANT so event generator should decay such
states to the usual "stable" particles (ie. no charm states).
- a interface for supplying "neutrinos" has been defined
(cf. get_nu_init.F and get_nu.F). a switch allow users to
select different methods; currently only a uniformly distributed
"far" detector method is in place.
- MKIN 2 now implemented. this method pulls a neutrino from
the flux routine; translates it into the detector co-ord
system; modifies the flavor (oscillations); determines the
mass distribution along the line of sight; calculates the
cross sections for the different materials (elements);
accept/reject the neutrino based on Total cross section
(all physics processes, all mass "seen"); picks material type
by relative cross sections (sigma_AZ*mass); chooses vertex;
calls neugen to pick sub-process (QE,NC,CC,resonance) and
kinematics. NEUGEN package is called to perform the
cross section calculation via sigma_tot_za.F, which in turn
calls "rsqlib". Shadowing and other nuclear effects should
be accounted for in this routine (currently simply calculates
cross section as z*sigma_nu_p + (a-z)*sigma_nu_n.
- 1996.07.24
- bug fixes demonstrated by moving to alternative platforms.
(generates divide by zero on alpha without these fixes)
- 1996.07.24
- update init/gminos.ffr - lots of annotations and examples. users
who have made modifications for their own copy (and wish to
retain them!) should make a backup copy before any "cvs update".
- 1996.07.25
- important fix to 1996.07.24 fixes...
- 1996.07.26
- neutrino flavor changing method installed
Jeff Nelson's code for entering single particles for kinematics
- 1996.08.02
- fix bug in GFMIXT (returned unrealistic values if material wasn't
a true mixture; now returns average A,Z with WMAT=1.0). This
affected only the composition of the pseudo-material that was
meant to represent the average of an active detector plane. And
only if one of the component materials was an element rather than
itself being a mixture (Aluminum combs triggered this). This
affects the calculation of the average plane A,Z, rad_length,
abs_length (but not density), and the selection of nuclei types
for neutrino events deemed to occur within the planes (For LSTs
with Al combs this is no more than 5% (when Fe=2cm) of the events).
- 1996.08.26
- open "hepevt.dat" file only if using that method for kinematics
- 1996.09.04
- improve indexing in filwin_mdlorg and filwin_planespec (remove
reminder message). add .cvsignore files to gminos and adamoinc
directories so that CVS doesn't complain so much about "new"
files.
- 1996.09.04
- new B-field maps. The grid versions require the user to checkout
a new module from CVS: "bfield" (at the same directory level
as "gminos") -- this area holds the datafiles. Users on UNIX
machines then need the env variable BMAPPATH to point at this
directory (VMS users should create a logical of the same name).
Currently supported maps numbers:
| BMAP
| description
|
|---|
| 0
| no field
|
| 1
| old radially dependent, azimuthally symmetric field
|
| 2
| 1.5T in the +y direction
|
| 100
| 2d grid version of "1" (5cm spacing)
|
| 101
| 2d grid version of perfect octogon (5cm spacing)
|
| 102
| 2d grid, ANL calculation of reference design (5cm spacing)
|
| 103
| 2d grid, LLNL calc, reference design (8cm spacing)
|
| 104
| 2d grid, LLNL calc, reference design (5cm spacing)
|
- 1996.09.20
- Modified /NEUCOM/ interface between NEUGEN and GMINOS. This means
the next CVS update on one requires you to update (and
compile) the other as well. No warning will be issued if you
forget, but you risk getting garbage in your NeuKin ESet output.
The NeuKin ESet (ADAMO table) has also been modified. It now
includes the individual vectors for a first and second electron
and a tau. They still contribute to the P4Shw shower vector,
but are available separately for interested parties. The energy
returned in P4Shw(4) is to sum of all the P4(4) of final state
non-neutrino, non-muon particles minus any baryon and/or ion masses.
Users can re-fill NeuKin in old data files from the information that
is already available by calling the following routine for each event:
subroutine update_neukin
implicit none
#include "NeuKin.inc"
#include "partap.inc"
call filcom_stdhep ! StdHep ESet -> /HEPEVT/
NeuKin_ID = 1
call GETTAB(NeuKin) ! get the info we wish to keep in place
c re-fill the 4-vector
call fillmushw_neukin(NeuKin_P4Mu1,NeuKin_P4Mu2,NeuKin_P4Shw,
& NeuKin_EMFrac,
& NeuKin_P4El1,NeuKin_P4El2,NeuKin_P4Tau)
call REPTAB(NeuKin) ! put the updated info back out in ESet
return
end
- 1996.09.20
- A few more plane defaults to handle various near detector
configurations.
- 1996.10.04
- Alternative file names for the FFREAD input data can be used
(on unix machines) by setting the environment variable
GMINOS_FFR to the desired name.
(NEUGEN) The rsq_xsects.dat file contains the calculated
cross sections that are used for interpolation. Writing this
file once, and simply reading it back in
significantly saves time on subsequent runs, but could
cause problems used in a multiplatform filesystem (ie. FNAL's afs).
The routine mc_fill_rsqdat now includes the string
found in MINOS_ARCH before the .dat. This
environment variable is set by default when using the standard
labyrinth.setup script.
- 1996.10.10
- Added geometry and hit structures for emulsion (E) planes.
These are tiled planes with a geometry structure different
from the LST, FLS, RPC configuration.
Other additions and clean up handling of geometry:
- remove dead code for multiple sets, by removing STDG
card (selecting from palettes).
- add LENGth as a valid characteristic to RTAG. Box volumes
no longer have to be square but can be rectangular.
- re-ordered parameters to vol2geant to be
more natural.
- 1996.10.22
- (NEUGEN) Rob Edgecock has made modifications to enable neugen
to accept many more elements and isotopes. Appropriate changes
have been made to routines in gminos to allow access to those
(A,Z) combinations. Additional elements have been added to
LBL_PART. A minor re-shuffling of the geant_ids.
- 1996.10.25
- New FFREAD card: PLNL (plane limits). The two integers
represent the first and last plane for an accepted event vertex.
- 1996.10.25
- Modifications to Make.include to deal with new GEANT321
version 05 (now "current" at FNAL). This version uses imake/cvs
which changed the directory/calling structure for the include files.
The environment variable GEANTVERSION=321 (tweaks the link
library name), while GEANTINCSAUX handles the funky chain of
"gcking.inc" including "geant321/gckmax.inc". This means that we
needed a -I${GEANTINCSAUX} flag so that the compiler can
find the internally referenced include file.
The new GEANT uses #include "geant321/xxxx.inc" everywhere.
The labyrinth.setup file at FNAL tries to set GEANTINCS appropriately
for whichever GEANT was setup (setup geant [current|old]).
- 1996.11.08
- New FFREAD cards:
- MTDK
- tau decay method, not yet implemented.
- CCNC
- select CC or NC only events in NEUGEN
- DKOF
- turn off individual particle decays in NEUGEN
- 1996.11.12
- Important bug fix. Users who updated their source code
from the repository after 1996.09.25 should retrieve the new
gminos/init/uginit.F. This ensures that sufficient
space is reserved for FFREAD's internal common block. FFREAD
overwriting into other common blocks may have contributed to random
looking crashes that some users have experienced.
- 1996.11.14
- Jeff Nelson's improvement to FLS digitization (difls_1cell).
This de-emphasizes the importance of slow neutrons, by reducing
the light yield of slow protons generated by n-p scattering in
the hydrogen rich scintillator. A saturation effect comes into
play for large dE/dx losses.
Note: users should remove any instances of PEMU
and ELMU cards (or update them to use 19.557 and 0.0028)
unless they clearly understand the implications.
- 1996.11.19
- Code for displaced vertices of heavy particles (tau's, D mesons).
Tau decay method handled by MTDK card.
Note: users that rely on the misc module to supply
STDHEP and LUND code, must update that module as well (now at
STDHEP 3.02 modified for transportability and a
standard version of JETSET74).
- 1996.11.29
- Better handling of un-set parameters available via FFREAD cards
(AIX incorrectly dealt with BSLP values). Also modified
geom.ddl so that BeamSystem ESet better reflects
useful (and used) information.
- 1996.12.04
- (RECO_MINOS) The reco_minos package can now handle writing records
back out.
The strip_hit.F shows how one can use these routines to
read in files, drop the ESets containing hits and write new files.
(RECO_MINOS) Alternative input filenames can be used (on UNIX
machines) without having to recompile the main routine.
Setting the RECO_GAF_LIST env variable to a legal file
name will use instead.
If not set, the program reverts to using reco_minos.gaf_list.
- 1996.12.06
- (RECO_MINOS) Two new "packages" for the analysis of events in the
FLS detector. You must use cvs -d $REPOS update -d
to get the new directories, and you must simultaneously
update Make.include in the gminos module. You'll
also need to define the env variable RECOPATH to point
at the top directory of the reco_minos module (similar
to GMINOSPATH and NEUGENPATH). Neither is
called automatically, users must explicitly call these routines.
correct_attn converts the raw photo-electron
count in FLSDigit_ADC(1) into a value with the best
guess removal of the light attenuation in the waveshifting
fiber; the new value is placed in FLSDigit_ADC(2).
Users that have non-standard values for the FATL, FATN, RMIR
or STUB cards should fill the /scintfiber/ common
block with appropriate values (cf. datafy_scintfiber.F).
find_vtx_fls is one attempt at a vertex finder.
Results are placed in the general structure: RecoVtx.
Authors of alternative vertex finders should consider using
the same structure as a common output.
- 1996.12.09
- NEUGEN/GMINOS flags for controlling whether NEUGEN performs
nuclear rescattering, QE and/or DIS scattering, individual
tau branching fractions. (Hugh Gallagher)
- 1996.12.20
- GUFLD: this version of gufld can be called in reco_minos for any
value of z without complaining about the point being outside of a
B plane (silently returns zero). Allows for "z" dependent slop
on how close in z precision you must be within the B plane to allow
for the limitations of the floating point.
The precision was set semi-empirically by seeing how far GEANT
tended to produce a xyz_local outside the alleged size of the plane.
Still no knowledge of the transverse size of planes; relies solely
on the map to set B=0 outside the iron (ie. never for field map 1).
Also protected BFLD_MAP1 from divide by zero fault.
- 1997.03.07
- "LST" (gas detector, now "APT") hits and digitizations have
changed (Lynn Miller). These reduced-size structures require
that gminos/adamoinc be CVS updated and the library
re-made.
- 1997.03.07
- Important bug fix Fixed significant error in FLS digitization.
A typo in the satturation correction affected the events in
unpredictable ways; unless the original hits were save it is
impossible to discern which ADC (digitizations) are affected.
The major effect is anomolously large ADC values (>> 300);
an intermediate fix is to truncate such values to a lower
threshold value (eg. set them equal to 300).
- 1997.03.14
- (RECO_MINOS) New version of reco.ddl, structures now include:
| RecoVtx |
Vertex information |
| RecoShwEnergy |
Shower energy information |
| RecoShwDir |
Shower direction information |
| RecoMuon |
Muon momentum, direction, origin |
| MuonTrkHeader |
Information about MuonTrkPoint table |
| MuonTrkPoint |
Space points associated with a muon |
FSLDigiCell
FSLDigiPln
FSLDigiView |
These ESets are used in correcting for the light
attenuation along the cell's length. They also include
some cell geometry information that may be handy during
reconstruction but is too easily regenerated and
too bulky to carry around in the base data structure
FLSDigit. |
- 1997.03.17
- Increase number of RVOL,RTAG,RVAL
and IVOL,ITAG,IVAL triplets that can
be specified in an FFR file. Add more control over tracking
in active media and iron (see Sxxx/Cxxx/Fxxx xxx=LST,FLS,RPC,EMU,FE
cards). Correctly handle planes with extrusions that are thinner
than the plane (eg. extra material used as strip plane support
on LST).
- 1997.04.15
- (RECO_MINOS) Improved FLS vertex finder -- remove bias due
to beam slope.
- 1997.04.15
- (RECO_MINOS) Incorporated FFREAD into RECO_MINOS. The card
file should be specified in the environmental variable
RECO_FFR (with a backup default of "reco.ffr").
New user hook routine reco_hist added. FFREAD cards
can be defined in reco_init, then the file is
processed automatically, further intialization processing
before events are read in can be done in reco_hist
(important in cases where values depend on the data cards).
For further information see the URL:
http://mimosa.astro.indiana.edu/minos/simulation/reco_minos.html
- 1997.04.18
- (RECO_MINOS) Fixed bug in close_gaf_chn that caused
segmentation violation on DEC machines. Improved file handling
so that processing the same file twice in a row no longer causes
"attempt to read beyond EOF" error (ie. the file is now properly
closed and re-openned); this might be necessary to select a
few non-consecutive individual events from a single file.
- 1997.05.07
- (RECO_MINOS) Add RecoShwEnergy, RecoShwDir,
MuonTrkHeader, MuonTrkPoint, RecoMuon
ADAMO tables to the list of tables to be cleared before we go
into the next event.
- 1997.05.13
- (RECO_MINOS) Move call to clear_reco in
reco_minos_loop to just before reading a new event
(rather than after reading and before user processing) in
anticipation of dataflows containing reco quantities.
- 1997.05.13
- (BFIELD) additional new B-field maps. See entry dated 1996.09.04
for previous entered maps. Dimension is the face-to-face diameter.
| BMAP
| description
|
|---|
| 111
| 2d grid, near detector 4m 15kAturns, 5cm grid, LLNL calc.
|
| 112
| 2d grid, near detector 5m 15kAturns, 5cm grid, LLNL calc.
|
| 113
| 2d grid, near detector 6m 15kAturns, 5cm grid, LLNL calc.
|
| 114
| 2d grid, near detector 8m 15kAturns, 5cm grid, LLNL calc.
|
- 1997.05.20
- Fix bug in get_nu_1 that incorrectly distributed the
neutrino flux over the face of the near detector. New near flux
file generated: nuhist_near_v2.rz.
- 1997.05.22
- (RECO_MINOS) Digitization code for APT/LST (gas detectors) entered
into repository. Use cvs update -d to pick up
the new lst directory. See example usr_ files
for usage.
- 1997.05.31
- (RECO_MINOS) find_vtx_lst is the corresponding
partner to find_vtx_fls for vertex finding in the APT/LST
gas detectors. The algorithm was written along similar lines
but with provisions for using both wire and strip information.
- 1997.05.31
- (RECO_MINOS) Digitization routine for APT/LST no longer require
user to book histogram 200 (and 200 is now free for user use).
- 1997.05.31
- (RECO_MINOS) New package mhpd Minos HPlot
Drawing. This is a relatively simplistic drawing package
for generating event pictures for sites without VINES installed.
It uses standard HPLOT/HIGZ routines and currently has few
controls. The user must call mhpd_init once in the
intialization phase, and mhpd_xyuv to display an
individual event (the user is also responsible for pausing the
processing at that stage). The display consists of 2 columns
of 4 views (X,Y,U,V). The left column displays the whole event,
while the right is zoomed in around the vertex.
The user must call mhpd_init once in the intialization
phase, and mhpd_xyuv to display an individual event
(the user is also responsible for pausing the processing at that
stage). The display consists of 2 columns of 4 views (X,Y,U,V).
The left column displays the whole event, while the right is zoomed
in around the vertex. Items on the screen are
- FLS digits
- red, open circles
- LST wires
- green, open boxes
- LST strips
- blue, open diamonds
- true vertex
- magenta (purple-red) cross
- found vertex
- cyan (light-blue) cross
The size of the marker for each digitization scales with its ADC
or Charge (linearly with an overall minimum size) that is set
independently for each type and for each event.
- 1997.06.18
- (RECO_MINOS) Small modifications to various routines -- with
no substantial changes. These were introduced to make
FTNCHEK happier; they involved either eliminating
unused variables, or moving all the data statements below
all the type declarations (and saves).
- 1997.06.18
- (RECO_MINOS) Modifications to mhpd (Minos HPlot
Drawing) package so formatting (what is drawn, symbols,
colors, sizes) are controlled via /mhpd_cntrl/ common.
- 1997.06.18
- (RECO_MINOS) New package reco_muon, for muon reconstruction,
both tracking and momentum determination.