From Andreas:
Adding a protection for the case ref. index < 1 for which mzbook is not called
to allocate storage for the derivative.
(index < 1 is sometimes used to define the opt. sensitive materials)

Add default .gitattributes.

Convert svn:ignore properties to .gitignore.

From Andrei Gheata:
Fixed a problem in gtckov.F reported by Maria Patsyuk (PANDA collaboration):
the step length, total step length and tof where not updated when photon
gets in a geometrical boundary with a medium without optical properties

Updating version number

From A.Morsch:
Fixed the problem of time of flight for optical photons
(Geant3 was using the phase velocity instead of the group velocity);
reported by M. Patsyuk (GSI)

Removing duplicate definitions (Ramona)

Adding missing functions

Updating version number before tag

Changes as per Savannah bug #84396

 In GetSecondary(Int_t isec, Int_t& partPDG, ...): return the PDG code instead
of G3 internal particle ID
- In DefineParticle() with short argument list: print warning about
deprecated function

Increase stack to 300 particles

Updating version number

Fixing previous wrong commit

Closing Savannah #90939

Fix from Federico

In TGeant3TGeo::GetMaterial:
Replacing call to Gfmate with Gfmate2 in order to avoid a warning about
deprecating Gfmate in TVirtualMC.

Preparing for next tag

Fix for icc

In TGeant3:
Fix in printing of Random seed **** GTRIGI: IEVENT ...

Adding GCONSX

- Adding new TVirtualMC functions in TGeant3:
    Bool_t CurrentBoundaryNormal(Double_t &x, Double_t &y, Double_t &z) const;
      ... prints only warning that this function is available only in TGeant3TGeo
    Bool_t GetMaterial(Int_t imat, TString& name,
                       Double_t& a, Double_t& z, Double_t& density,
                       Double_t& radl, Double_t& inter, TArrayD& par);
     ... which will replace the current Gfmate(...) functions
- Adding a new TGeant3::Gfmate2(..) function with corrected call to g3fmate()
  and reimplemented the existing TGeant3::Gfmate(..) via a call to the new
  Gfmnate2 and with adding a warning about deprecating of this function in the
  TVirtualMC interface; also replaced the internal calls to Gfmate(..)
  via Gfmate2()
- Updated for the clean-up of G3 specific visualization functions in TVirtualMC
  (in G3Volume.cxx)

Fix a typo (gconsx instead of cconsx) in the declaration of common GCONSX

Fix in Gfmate with Double_t arguments.

From Lia Lavezzi:
     we realized that there is a little change to be done to the routine
trprfn.F, due to our last fix with IERR = 5 (for the cases when start and
end points of the step coincide): the problem is that the check for
coincidence (IERR = 5) is made before the check for 0 momentum (IERR = 2).
For this reason at low energies when the particle loses all its kinetic
energy and should stop with IERR = 2, it keeps trying to go on and hitting
the same point (IERR = 5) . This does not affect the results, but it is a
wrong behaviour.
With the fix, that just switches the two checks, the tests on the
particle propagation are made in the correct order.

replace fabs by TMath::Abs

From Lia.Lavezzi@pv.infn.it  and alberto.rotondi@pv.infn.it

in PANDA we prepared some additional bug fixes for GEANE

There are two corrections:
1) the propagation to length used to contain a bug so that if the
propagation to length was chosen but the given length was too small
(or 0, as limit case) it did not work.
We fixed it in such a way that if the track length is too small but not 0
the step size is set to PREC, if it is <= 0 it is set to 0 and the start
point is coincident with the end point.
The flag IERR has been provided of the additional value 5: in this case,
when start and end position coincide, IERR is set to 5 and this means that
the error matrix does not have to be updated. In this way if track length
is 0 the particle just does not move.
2) it has been added the possibility to track along the z axis, though
with an approximation. Up to now, due to the matrix elements with
1/cos(lambda), the propagation exactly along z axis (or below a certain
limit) was forbidden. The approximation we introduced is in TGeant3gu.gxx,
in eustep() function: if the theta angle drops below 1e-4 deg it is forced
to this limit value in order to be able to go on tracking. Obviously this
means that if you try to track exactly along z (theta = 0) you will be
forced to track at 1e-4 deg instead.

We checked the pulls and residual values for a plane @12m and @26m and
they both with and without magnetic field and they are acceptable.

From Andrei Gheata:
Implemented a new TVirtualMC function (in TGeant3TGeo):
Bool_t CurrentBoundaryNormal(Double_t &x, Double_t &y, Double_t &z) const;

Revrting previous commit (in TGeant3TGeo) as it caused problem in ALICE simulation.

Fix formats

In TGeant3TGeo:
Adding creating of TGeoManager in ctor; it is required when the VMC application builds geometry via TVirtualMC and then selects TGeant4TGeo as MC.

Updated dist target in Makefile:
Now using the new makedist.sh script (common for all VMC packages)

From Lia.Lavezzi@pv.infn.it  and alberto.rotondi@pv.infn.it :
The patch contains the following corrections:

1. fix for IERR = 3 error - trprfn.F/erprop.F
In case of IERR = 3 error, i.e. "H*ALFA/P AT X1 AND X2 DIFFER TOO MUCH",
there was a GOTO pointing back to a part of code where uninitialised
variables were used. This used to produce a crash.
Our solution is to skip the step where IERR = 3, getting rid of the wrong
GOTO, but to update the initial point position, momentum and magnetic
field value as well.

1.b In addition, related to this point, we put a new flag in the common
ERTRIO1, NOPRNT, to select whether printing or not the error messages
about IERR, that in some cases could be annoying.
By default the printings are there, but if the user sets the flag NOPRNT
to 1 they are skipped.

2. fix for crash in electron tracking - erland.F
A crash happens, though very rarely, in the part recently inserted to
include bremsstrahlung
        LX = 1.442695*STEP/RADL
        EMB=E/(2**LX)
        S2B=E*E*(1/3**LX - 1/4**LX)
        SB=SQRT(ABS(S2B))
        DEDXB = 1.2*SB
The calculation of **LX fails when LX is too big (i.e. in thick materials,
where RADL is small, or when big STEPs are performed): our solution is to
put DEDXB = 0 and recalculate it only for LX < 25.

3. fix for missing virtual detector planes - ertrch.F
In the code there is a check to understand whether in a step the plane is
reached or not: if the step is too big and the ending point is on the
other side of the plane the step size is reduced so that the ending point
of the plane is on the plane itself (or near to it within a certain
precision). The statement says IF(STEP.GE.ASCL1) and STEP is the step
size, ASCL1 the distance between the starting point and the plane.
Sometimes it happens that the ending point is right before the plane, so
that STEP is not "greater or equal" to ASCL1, but their difference is very
small (of the order of -10-7). In this case GEANE acts as there is enough
space to perform this step and then tries another one before reaching the
plane but doing this it fails in one of the checks it does at the
beginning of the step. It compares some quantities with the variable PREC
(of the order of 10-4) and there is the failure: we have too small
quantities here (our 10-7/-8)! Moreover it looks for geometrical
boundaries, but we are talking about virtual detector planes and it does
not find them and so it decides to perform too big steps and it does not
find the plane: ever!
Our solution: the idea is to consider the plane as reached if the distance
of the ending point to/from the plane is less than PREC.

4. There are, in addition:
   - a couple of corrections to our previously sent patch
   - a fix in TGeant3/TGeant3.h, where in the struct Ertrio the variable
iepred was not set up as a vector, while in the corresponding FORTRAN
COMMON ERTRIO it is IEPRED[MXPRED].

The patch has been deeply tested by us in Pavia and has been distributed
to the PANDA computing group to have a feedback from a large group of
people in case there were problems with it.

From Andrei:
The patch is replacing the compilation define COLLECT_TRACKS with a normal data member of TGeant3TGeo. Added also the interface Set/GetCollectTracks() at the level of VMC, which works only for TGeant3TGeo but leaves room for the implementation in TGeant4 (if using TGeo navigation)

From Andrei:
Typo error fixed and missing include added in TGeant3gu.cxx

From Andrei:
Mother track id searched now recursively (C.Reed)

Remove -fPIC from the makefile. This option should only be specified in
each architecture Makefile (alraedy done0

From Mohammad Al-Turany:
Moving the eustep and related commons to the C++ part of the VMC, and creating a new stepping action for Geane (GeaneStepping) in the MC application. Users can implement this method in there application and control or interact with the Geane propagation  step wise.

From Andreas Morsch:
add functions to skip neutrinos

From Andreas Morsch:
print a more informative message in case of warnings or errors.

From Lia.Lavezzi@pv.infn.it
in some occasions, when propagating particles starting from a small theta
angle, though different from 0, the code shows instabilities and a
floating point exception shows up. We identified the problematic lines in
trsdsc.F and trscsd.F and protected them with an IF statement:
1) the first problem was in the calculation of parameter lambda:
PC(2)=ASIN(TN(3))   (in trscsd.F)
in case when ABS(TN(3)).GT.1. During the EVO meeting it has been decided
to set it "by hand" to 1 (or -1) since the difference is very small, due
to roundings. We changed the code accordingly;
2) the second problem was in the filling of the error matrix:
RC(J)=S(J)   (in trscsd.F)
RD(J)=S(J)   (in trsdsc.F)
due to the fact that the left side variable is in single precision
while the right side one in double precision. The problem arises
when ABS(S(J)).GT.1E30. In this case we flag it with IERR.
When the flag reveals a problem arose in these routines we hadle it
in ertrak.F. In order to do this we added a new labelled common
(ERTRIO1) to the existing ERTRIO, filling it with 3 variables: IERTR,
IERTR1 and IERTR2.
For the moment we use only IERTR to handle the mentioned floating point
exception and stop the tracking of the particle in such cases.
We also changed TGeant3 in order to be able to retrieve the iertr value
from the interface in C++ (authors of the changes: C. Hoeppner, L. Lavezzi
and A. Rotondi).

From Federico:
Replace a very inefficient use of std::string consuming 14/% of the Alice simulation

Adding new TVirtualMC function:

virtual void SetUserParameters(Bool_t isUserParameters);

which allows to activate the parameters defined in tracking media
(DEEMAX, STMIN, STEMAX), which are, be default, ignored.

set version 1.10

From Susan Kasahara
function igiwty must return a value

From Alessandra Filippi
Fix for a systematic shift in reconstructed momentum.
see Alessandra's slides at the  VMC meeting at
http://indico.cern.ch/conferenceDisplay.py?confId=41169

s line, and those below, will be ignored--

M    erdecks/ertrch.F

In SetDecayMode():
Adding call to SetUserDecay() before redefinition
of decay mode. (G3 will not issue a warning about
overriding the decay mode.)

From Ivana:
fixes to previous patch

 From Federico & Ivana:
Fix a problem when Runge Kutta is selected in E02 example

From Federico:
Add a comment in GUFLD:
Please note that starting from svn revision 232, the  argument of this routine
have been moved from single to double precision. The corresponding changes have
been done to the calling routines. This change is backward incompatible and to
continue working with single precision  arguments, the compilation flag
-DSINGLEFIELD should be added to the compilation, both FORTRAN and C++.

From Federico:
Changes to support teh new field abstract interface

From Ivana:
Add new interface to the mag field

From Dario Berzano:
patch to the Geant3 Makefile in order to make it compile
on x86_64/amd64 arch. Basically, the patch adds -fPIC to the CXXFLAGS,
CFLAGS and FFLAGS. This patch is harmless on x86.

From Andrea Fontana:
- bug fix in ERTRCH: there was a bug in the calculation of the energy lost
  by a particle in the first step into a new medium. It is now fixed.

- change in the treatment of the truncation of the Landau distribution
  in ERLAND.F: a new approach is followed, similar to what is done by GEANT
  (see GEANT3 manual W5013), to truncate the Landau distribution so that
  the mean value of the lambda parameter corresponds to the mean energy
  loss. The urban distribution is truncated to have the same sigma
  of the Landau at the matching condition.

- new treatment of energy loss fluctuations for bremsstrahlung in ERLAND.F:
  this is due to Alberto Rotondi that implemented in GEANE the treatment of
  energy loss fluctuations for electrons. The results are in agreement with
  the theoretical treatment by Bethe and Heitler and allow the use of the
  track follower also for electrons.

In Makefile:
use root-config --etcdir.

From m.al-turany@gsi.de
"I am trying to interface Geane with EVE,  The first step in that is to be
able to record points from the Geane tracking, so I need access to the
stepping function in the VMC Main application, which is called usually from
gustep, but in Geane eustep is called instead! So what one could do is to
re-write the eustep in C++ and add it to the TGeant3gu.cxx file and call the
stepping from there  or simply call the gustep at the end of eustep
subroutine, the later has the advantage that if the track visualization is
switched on one will get the TGeoTracks immediately for Geane."

Increment version number to 1.9

TGeant3TGeo requires now to include TGeoMatrix.h explicitely.

From Federico:
Fixes for gfortran/g77/g95

From Federico:
- If the max no of steps is negative then the track is abandoned silently. This is necessary for us because we are producing huge output files with useless messages. The change is backward compatible

- If ABANDON is set to 0, then the photoelectric photon above thresholds are produced and put on the stack. This is a real bug fix, because now photons are simply not tracked if their range is smaller than safety, even if ABAN=0. This behaviour is inconsistent with the existence of the ABAN flag and may give wrong physics in gaseous detectors. Marek can explain you that only with this change he could obtain correct results for TPC simulation

From Christian Holm:
Cannot use jbit with gfortran.

Set property svn:eol-style LF on all source and Makefiles. This should avoid
problems with Win32 line endings ending up in the repository. All MS tools
support LF eols fine.

Using common makefiles from $ROOTSYS/etc/vmc instead
of those in config.

Moving file one level up, as config directory will be removed

Using common makefiles from $ROOTSYS/etc/vmc

- Updated for modified interface TVirtualMC::DefineParticle() - added more
  arguments needed for Geant4 implementation
- In Init(): added call to new MCApplication function AddIons()
From Mohammad Al-Turany:
- Implemented new function SetDecayMode()

From M.Al-Turany@gsi.de:
Update the README file explaining how to use GCalor with the VMC.
Add a new file gcalor_vmc.diff that must be applied to the Gcalor code downloaded
from the GCalor web site.

- New implementation of ProdProcess by Sue Kasahara:
  To make use of GCKING/KCASE variable for determining the production
  mechanism of the secondaries.  The old method was to pick the first
  active process from the current step's list of active processes
  that had the capability of generating secondaries.  This occasionally
  picked the wrong secondary production mechanism.

- In StepProcesses(): fill the array with kPNull if there are no
  active processes

- In G3toVMC(): use of new TMCProcess constants: kPHCElastic, kPHIElastic;
  the 'NEXT' mechanismus mapped to kPTransportation

From Ivana and Andreas:
In DefineParticles(): updated mass and lifetime definitions for
additional particles.

add svn:keywords Id property. Missed initialy due to wrong assumption in
patch up script.

From Marek Kowalski:
I have adde the ILOSS=6 parameter which
deals with the energy loss model for the Kr decay simulations.--This line, and those below, will be ignored--

M    gphys/gfluct.F

From Andrea Fontana:
 simplified the algorithm to find the PCA which is now more robust.

From Andrae Fontana:
 added method GetClose to get the last three points of extrapolation
 to the point of closest approach.

From Ivana:
Changed the fPDGCode type to TArrayI to allow
its dynamic size (thanks to Susan Kasahara).

From Ivana:
Do not add particles in TDatabasePDG if they
are were already aded;
(thanks to Susan Kasahara for this suggestion)

From Ivana:
Implemented the new pdg "standard" codes for ions
defined in pdg-2006, thanks to Susan Kasahara
for this suggestion.

Increment version number to 1.8

From Peter Hristov:
a small modification in gphys/gvaviv.F to avoid using non-initialized data.

From Andrei:
Added optional method:
 Bool_t TVirtualMCApplication::MisalignGeometry(). To be implemented by user application in case misalignment need to be applied to TGeo after ideal geometry was built.

From Ivana & Yuri Shitov:
Due to rounding errors COSTH could be greater than 1. set it to 1 in this case.

New README file extending README_geane with the new functionality
developed by Andrea Fontana and Alberto Rotondi.

A major update of GEANTE by Andrea Fontana and Alberto Rotondi

1) update of the Coulomb multiple scattering parametrization;
2) update of the straggling of energy loss for thin materials;
3) new options to extrapolate the track parameters to the point
   of closest approach to a point or to a wire (straight line).

Details on the physical motivation behind this work can be found
in our report for the Panda Collaboration, available at:

http://www.pv.infn.it/~fontana/tracking.pdf

Feel free to contact us for questions and discussions about these
topics by using the following email addresses:

alberto.rotondi@pv.infn.it
andrea.fontana@pv.infn.it

---

List of changes in the fortran and C++ routines of the geant3
VMC directory:

- gcmore.inc
  gtmore.inc
  geant3LinkDef.h
  gcomad.F

 Added a new common that contains all the new variables:
      COMMON/GCMORE/GCALPHA,ICLOSE,PFINAL(3),DSTRT,WIRE1(3),WIRE2(3),
     +              P1(3),P2(3),P3(3),CLENG(3)

     input to ERLAND:
      GCALPHA: energy cut parameter for energy loss fluctuations

     input to EUSTEP:
      ICLOSE: = 1 the use of the common is enabled for the closest
                  approach to a point PFINAL(3)
              = 2 the use of the common is enabled for the closest
                  approach to a wire of extremes WIRE1(3) and WIRE2(3)
              = 0 the common is empty and disabled
      PFINAL(3): assigned point
      DSTRT: assigned distance between initial point in ERTRAK
             and PFINAL along straight line (currently noy used)
      WIRE1(3): first point of a wire
      WIRE2(3): second point of a wire

     output from EUSTEP:
      P1(3): point previous to the point of closest approach to
             PFINAL() or wire
      P2(3): point of closest approach to PFINAL() or wire
      P3(3): point next to the point of closest approach to
             PFINAL() or wire
      CLENG(3): track length to the previous 3 points

      Important note: the calculated points of closest approach are
      depending on the GEANE steps. For calculating the true point
      of closest approach the last 3 points of the extrapolation, i.e.
      the previous to closest, the closest and the next to closest are
      returned to the user. Different algorithms can be implemented, but
      we decided to leave this to the users in the C++ interface to GEANE.

- ermcsc.F
 new expression for the variance of the Coulomb multiple scattering
 according to Fruhwirth and Regler, NIM A 456 (2001) 369

- ertrch.F
 added DESTEP in the calling string of ERLAND for calculation with
 Urban model. Added and saved previous step PRSTEP.

- erland.F
 added new calculation for sigma of straggling in energy loss
 to include in Geane the Urban/Landau approximation, as explained
 in the Geant manual and related papers.
 The model parametrization can be controlled with a user variable (GCALPHA)
 in the new GCMORE common block: 1 is for old gaussian model valid
 for dense materials, other values (see the report) are for gaseous
 materials.

- eustep.F
 added the calculation to the distance of closest approach to a point
 or to a wire.

- TGeant3.h
- TGeant3.cxx
 added the possibility to define user cuts (already present in the gccuts
 struct but not in the TGeant3::SetCUTS method) and to define the new
 variables of the GCMORE common with two new methods SetECut() and
 SetClose().
 Added new method InitGEANE() to initialize GEANE to the old behaviour
 (default) for backward compatibility. Only the multiple scattering has
 been updated to a more correct formula.
 Corrected a typo in the call to the routine Trscsd().

Make version v1-7

from Peter Hristov:
the behavior of the g77 compiler can be fixed using the
-fno-f2c option.

 From Andrea Fontana and Alberto Rotondi
Added one extra parameter (GCALPHA) to the COMMON GCCUTS to drive the
Urban/Landau models in the ERLAND routine.

 From Andrea Fontana and Alberto Rotondi
Added new calculation for sigma of straggling in energy loss to include
in Geane the Urban/Landau approximation, as explained in the Geant
manual and related papers.
The model parametrization can be controlled with a new user variable
(GCALPHA) in GCCUTS common block: 1 is for old gaussian model valid for
dense materials, other values (see the enclosed report) are for gaseous
materials.

 From Andrea Fontana and Alberto Rotondi
Added DESTEP in the calling string of ERLAND for calculation with Urban
model and call to ERLAND moved at end of routine

 From Andrea Fontana and Alberto Rotondi
New expression for the variance of the Coulomb multiple scattering
according to Fruhwirth and Regler, NIM A 456 (2001) 369

CVS-- -- ------------------------------------------------------------------
: ----------------------------------------------------------------------

Fix a problem when adding a new tracking medium to the TObjArray.

From Ivana Hrivnacova and Andrea Fontana:
Reintroduce functionality in TGeant3::SetCuts that was removed in a previous patch.

From Ivana:
Restore the function TGeant3::MediumId

From Federico:
Suppress compiler warnings coming from unused arguments

Add missing include TList.h

Add a description of the gcalor interface
=========================================
The GCALOR package contains the original CALOR simulation code (including the low energy neutron code MICAP ) developed and maintained by the Oak Ridge National Laboratory ( ORNL ) and an interface to utilize this code in the framework of GEANT .
GCALOR consists of

HETC
The High-Energy-Transport-Code is transporting charged hadrons up to an energy of 10 GeV through the materials of the setup
FLUKA
The FLUKA fragmentation model is utilized for interactions above the HETC limit
MICAP
The neutron code from the Monte-carlo-Ionization-Chamber-Analysis-Program is called for neutrons with a kinetic energy below 20 MeV. The simulation of interactions is based on measured/calculated cross sections and secondary particle energy/angular distributions

Interface
---------
The interface between CALOR and GEANT determines which code to call for the current particle. The material description is automatically extracted from GEANT The transport of electromagnetic interacting particles is performed by GEANT as well as the tracking and the calculation of the energy deposition.

Installation:
------------
 1. download the CALOR code (see below) and copy it to  geant3/gcalor directory
 2. in the geant3/Makefile  add the gcalor directory to the GDIRS
recompile geant3

Usage with vmc
--------------
 one has call TGeant3::SetHADR(5) in the initalization of Geant3.

For the initialisation, GCALOR needs the following  additional cross-section data (expected files: chetc.dat and xsneut.dat in the $G3SYS/lib directory) for all needed materials, GCALOR Calculate the macroscopic cross-section for low energy neutrons and write relevant informations for the user to file calor.out

The GCALOR package has been tested on different computer platforms:

Unix systems
DecAlpha, Silicon Graphics, Hewlett Packard, IBM-R6000, SUN Sparc and Linux on PC
DecVMS
DecALpha and VAX
MVS
IBM Mainframe 3090

The GCALOR does not work for the moment on 64 bit arcitecture. For the 64 bit, it is needed to compile with the option "-m32".

For download and  more detailed information:  http://www.staff.uni-mainz.de/zeitnitz/Gcalor/gcalor.html

Add a description of the geane interface
========================================
The Geane package allows the user to calculate the
average trajectories of particles and to calculate the
transport matrix as well as the propagated error covariance
matrix. Geane is a set of routines worked out
by the European Muon Collaboration [1] and it is integrated
to the GEANT3 system [2].
See geant3/doc/gedoc for more details.

  The user should invoke the routine Ertrak and Eufill or Eufilp or
Eufilv, to carry out the tracking.
In addition to this a series of utilities are available for the user
(e.g. to transform the track representation from one
system to another or to carry out 5 X 5 matrix multiplication
in an optimal way).

Track variables, Representations
--------------------------------
      The particle trajectory is characterized by 5 independent
variables as a function of one parameter (e.g. the pathlength).
Among the 5 variables 1 is related to the curvature (to the absolute
value of the momentum, p), 2 are related to the direction of the
particle and the other 2 are related to the spatial location.
The most usual representation of these 5 parameters are:

       I.       1/p, lambda, phi, y_perp, z_perp

where lambda and phi are the dip and azimuthal angles related
to the momentum components in the following way:

            p_x = p cos(lambda) cos(phi)
            p_y = p cos(lambda) sin(phi)
            p_z = p sin(lambda)

y_perp and z_perp are the coordinates of the trajectory in a
local orthonormal reference frame with the x_perp axis along the
particle direction, the y_perp being parallel to the x-y plane.
This representation is usually applied in the overall reference
frame. (In the EMC code this reference frame is labelled by 'SC'
since the overall system was identified with that of the Streamer
Chamber.)

       II.       1/p, y', z', y, z

where y'=dy/dx and z'=dz/dx. This representation is particularly
useful in fixed target experiments, where the trajectory is evaluated
on successive parallel planes (which are perpendicular to the x-axis).
(In the EMC code this representation is labelled by 'SP' since a
convenient mathematical description of a trajectory being approxima-
tely parallel to the x-axis is a 'spline'.)

       III.       1/p, v', w', v, w

where v'=dv/du and w'=dw/du in an orthonormal coordinate system with
axis u, v and w. This representation is paricularly useful when the
trajectory has to be evaluated on different detector planes
in a colliding beam experiment, where the planes can take a great
variety of directions.(In the EMC code this representation is
labelled by 'SD' as System of Detection.)

      Of course, all the above representations of the trajectory
are equivalent and one can go from one representation to the
other by calculating the corresponding Jacobian. These Jacobians
are provided by the following EMC routines:

     S/R Trscsp    from I   to II
     S/R Trspsc    from II  to I
     S/R Trscsd    from I   to III
     S/R Trsdsc    from III to I

User Interface
---------------------------------------------------
      Ertrak(const Float_t *x1, const Float_t *p1,
     const Float_t *x2, const Float_t *p2,
     Int_t ipa,  Option_t *chopt)
---------------------------------------------------
   Performs the tracking of the track from point x1 to  point x2
   (Before calling this routine the user should call Eufil(l/p/v)

          x1       - Starting coordinates (Cartesian)
          p1       - Starting 3-momentum  (Cartesian)
          x2       - Final coordinates    (Cartesian)
          p2       - Final 3-momentum     (Cartesian)
          ipa      - Particle code (a la GEANT) of the track

          chopt
              'B'   'Backward tracking' - i.e. energy loss
                                 added to the current energy
              'E'   'Exact' calculation of errors assuming
                                 helix (i.e. pathlength not
                                 assumed as infinitesimal)
              'L'   Tracking upto prescribed Lengths reached
              'M'   'Mixed' prediction (not yet coded)
              'O'   Tracking 'Only' without calculating errors
              'P'   Tracking upto prescribed Planes reached
              'V'   Tracking upto prescribed Volumes reached
              'X'   Tracking upto prescribed Point approached

------------------------------------------
Eufill(Int_t n,Float_t *ein,Float_t *xlf);
------------------------------------------

    User routine to set the input values  for chopt = 'L'

         n     Number of predictions where to store results
         ein   Input error matrix
         xlf   Defines the tracklengths which if passed the
                      result should be stored

----------------------------------------------------------------
Eufilp(const Int_t n, Float_t *ein, Float_t *pli, Float_t *plf);
----------------------------------------------------------------
 User routine to set the input values  for chopt = 'P'

 n     Number of predictions where to store results
         ein   Input error matrix (in the 'Plane' system )
         pli   Defines the start plane
                      PLI(3,1) - and
                      PLI(3,2) - 2 unit vectors in the plane
         plf   Defines the end plane
                      PLF(3,1,I) - and
                      PLF(3,2,I) - 2 unit vectors in the plane
                      PLF(3,3,I) - point on the plane
                                   at intermediate point I
---------------------------------------------------------------------
Eufilv(Int_t n, Float_t *ein, Char_t *namv, Int_t *numv,Int_t *iovl);
---------------------------------------------------------------------
    User routine to set the input values for chopt = 'V'

         n     Number of predictions where to store results
         ein   Input error matrix
        cnamv  Volume name of the prediction
        numv   Volume number (if 0 = all volumes)
        iovl   = 1  prediction when entering in the volume
               = 2  prediction when leaving the volume

--------------------------------------------------------------------------
void Trscsd(Float_t *pc,Float_t *rc,Float_t *pd,Float_t *rd,Float_t *h,
    Float_t *ch,Int_t *ierr,Float_t *spu,Float_t *dj,Float_t *dk);
--------------------------------------------------------------------------
transforms error matrix  from   sc   variables (1/p,lambda,phi,yt,zt)
          to   variables (1/p,v',w',v,w)
        pc[3]     1/p,lambda,phi                          input
        pd[3]     1/p,v',w'                              output
        h [3]     magnetic field                          input
        rc[15]   error matrix in   sc   variables        input     (triangle)
        rd[15]    error matrix in 1/p,v',w',v,w          output     (triangle)
        ch        charge of particle                      input
                  charge and magnetic field are needed
                  for correlation terms (v',yt),(v',zt),(w',yt),(w',zt)
                  these correlation terms appear because rc is assumed
                  to be the error matrix for fixed s (path length)
                  and rd for fixed u
        dj[3]     unit vector in v-direction
        dk[3]     unit vector in w-direction    of detector system

        ierr  =   1       particle moves perpendicular to u-axis
                         ( v',w' are not defined )
        spu       sign of u-component of particle momentum   output

--------------------------------------------------------------------------
void Trsdsc(Float_t *pd,Float_t *rd,Float_t *pc,Float_t *rc,Float_t *h,
    Float_t *ch,Int_t *ierr,Float_t *spu,Float_t *dj,Float_t *dk);
--------------------------------------------------------------------------
transforms error matrix from  variables (1/p,v',w',v,w)
 to    sc   variables (1/p,lambda,phi,yt,zt)

      pd[3]     1/p,v',w'                               input
      pc[3]     1/p,lambda,phi                         output
      h[3]      magnetic field                          input
      rd[15]    error matrix in 1/p,v',w',v,w           input      (triangle)
      rc[15]    error matrix in   sc   variables       output      (triangle)
      ch        charge of particle                      input
                charge and magnetic field are needed
                for correlation terms (lambda,v),(lambda,w),(phi,v),(phi,w)
                these correlation terms appear because rc is assumed
                to be the error matrix for fixed s (path length)
                and rd for fixed u
      dj[3]     unit vector in v-direction
      dk[3]     unit vector in w-direction    of detector system

      ierr              not used
      spu       sign of u-component of particle momentum    input

--------------------------------------------------------------------------
void Trscsp(Float_t *ps,Float_t *rs,Float_t *pc,Float_t *rc,Float_t *h,
    Float_t *ch,Int_t *ierr,Float_t *spx);
--------------------------------------------------------------------------
transforms error matrix  from   sc   variables (1/p,lambda,phi,yt,zt)
         to  spline variables (1/p,y',z',y,z)
      pc[3]     1/p,lambda,phi                          input
      ps[3]     1/p,y',z'                              output
      h[3]      magnetic field                          input
      rc[15]    error matrix in   sc   variables        input     (triangle)
      rs[15]    error matrix in spline variables       output     (triangle)
      ch        charge of particle                      input
                charge and magnetic field are needed
                for correlation terms (y',yt),(y',zt),(z',yt),(z',zt)
                these correlation terms appear because rc is assumed
                to be the error matrix for fixed s (path length)
                and rs for fixed x

      ierr  =   1       particle moves perpendicular to x-axis
                       ( y',z' are not defined )
      spx       sign of x-component of particle momentum   output

--------------------------------------------------------------------------
void Trspsc(Float_t *ps,Float_t *rs,Float_t *pc,Float_t *rc,Float_t *h,
    Float_t *ch,Int_t *ierr,Float_t *spx);
--------------------------------------------------------------------------
transforms error matrix  from spline variables (1/p,y',z',y,z)
         to    sc   variables (1/p,lambda,phi,yt,zt)
      ps[3]     1/p,y',z'                               input
      pc[3]     1/p,lambda,phi                         output
      h[3]      magnetic field                          input
      rs[15]    error matrix in spline variables        input      (triangle)
      rc[15]    error matrix in   sc   variables       output      (triangle)
      ch        charge of particle                      input
                charge and magnetic field are needed
                for correlation terms (lambda,y),(lambda,z),(phi,y),(phi,z)
                these correlation terms appear because rc is assumed
                to be the error matrix for fixed s (path length)
                and rs for fixed x

      ierr              not used
      spx       sign of x-component of particle momentum    input

Examples of Application
=======================
     1. Representing the trajectory at another point
     2. Joining track elements in different parts of the detector
     3. Prediction of the trajectory
     4. Fitting trajectory parameters
See geant3/doc/gedoc for details.

[1] W.Wittek, EMC Internal Reports: EMC/80/15,   EMCSW/80/39,
                                    EMCSW/81/13, EMCSW/81/18
    A.Haas,   The EMC Utility Package: UTIL42

[2] R.Brun, F.Bruyant, M.Maire, A.C.McPherson, P.Zanarini
    DD/EE/84-1, May 1986

from Mohammad Al-Turany & Denis Bertini

Changes in  TGeant3/TGeant3.cxx and TGeant3.h
------------------------------------
1. Geane interface functions are added:
    void  eufill(Int_t n, Float_t *ein, Float_t *xlf);
    void  eufilp(const int n,Float_t *ein,Float_t *pli,Float_t *plf);
    void  eufilv(Int_t n, Float_t *ein, Char_t *namv, Int_t *numv,Int_t *iovl);
    void  trscsd(Float_t *pc,Float_t *rc,Float_t *pd,Float_t *rd,
 Float_t *h,Float_t *ch,Int_t *ierr,Float_t *spu,Float_t *dj,Float_t *dk);
    void  trsdsc(Float_t *pd,Float_t *rd,Float_t *pc,Float_t *rc,
                           Float_t *h,Float_t *ch,Int_t *ierr,Float_t *spu,Float_t *dj,Float_t *dk);
    void  trscsp(Float_t *ps,Float_t *rs,Float_t *pc,Float_t *rc,Float_t *h,
   Float_t *ch,Int_t *ierr,Float_t *spx);
    void  trspsc(Float_t *ps,Float_t *rs,Float_t *pc,Float_t *rc,Float_t *h,
   Float_t *ch,Int_t *ierr,Float_t *spx);

2. The Gfang function wrapper is added
    void  g3fang( Float_t *, Float_t &,Float_t &, Float_t &, Float_t &,Int_t & );

changes in TGeant3/TGeant3gu.cxx
--------------------------
Adding GCalor interface
1. function calsig() and gcalor() are added
2. setting ihadr=5 will call the GCalor routine

changes in TGeant3/TGeant3.h
-----------------------
1. Structures for Geane output are setted as public so that the user can access them

    Ertrio_t *fErtrio
    Eropts_t *fEropts
    Eroptc_t *fEroptc
    Erwork_t *fErwork
    Trcom3_t *fTrcom3

2. The size of the error matrix errin is corrected to 15

From Ivana Hrivnacova:
Implemented new VMC function:
MediumId(const const Text_t* mediumName)

Change speed of light to agree with the current PDG value,
ie change from 299,792,510 to 299,792,458

Stamp version 1.6

From Peter Hristov:
Two new Makefiles