Oct 28 2001
Fix on some undefined variables introduced. Bug could lead to crashes, but 
otherwise was not charmful. On most platforms results remain unchanged due to fix. 
                      April 26
Fix of missing parameter NNEUTR=3 introduced, if the code will be used 
to reweight just electron neutrino then NNEUTR=1 should be taken. This 
possibility was never checked so far.

                      Dec 11 2000
This directory consists of info how to install package for calculating
anomalous couplings from KORALZ to KKMC

At present only anomalous W-W-gamma couplings are included [1].
Anomalous tau-tau-gamma couplings [2] are not installed, even tough 
it should be straightforward. 

Step to proceed:

1) place directory  nunulib in the main directory of KKMC
2) take care that the makefile of your main program loads the library
   nunulib/glib.a . If you are lucky, copying Makefile.ffbench into 
   Makefile of directory ffbench will do, but it will work on linux redhat 
   6.1-7.0 only, and on package version of about Dec 2000.
3) In general  you will need to take care of compiler loader options in
   Makefile.ffbench     (to be placed in ffbench) 
                        alternativelly original Makefile for your main program
                        needs to be enriched with loading nunulib/glib.a 
library  
   and 
   makefile              of nunulib 
   by hand.
4) Initialization of the package can be done with the adopted version of 
   SUBROUTINE  kzphynew, see reference [1].
   IT must be done after complete KKMC initialization and before first
   event is generated.
   The sequence is:
      npr(15)=1
      call kzphynew(XPR,NPR)
   where:
      REAL*4   XPR(40)
      INTEGER  NPR(40)
   The remaining entries for input parameters are either completelly dummy,
   or prepared for tau-tau-gamma option (reference [2]).
5) Calculation of the weights can be invoked with the call 

          wtdum=  WTANOM(Dumm)
   after generation of the neutrino event. All variables double precision.
   Then, common block 
      common /kalinout/ wtkal(6)
      real*8 wtkal
   will be filled up as defined in ref [1], note that in KKMC case, the SM 
   WW-gamma couplings need not to be calculated by the library  thus, the 
   anomalous weights wtkal(2) to wtkal(6) need to be divided by the 
   (standard model) wtkal(1).
   This solution is in accordance to compatibility with ref [1].

6) WARNING: It is assumed that all 3 flavours of neutrinos are generated
            thus, the weights must be calculated and used for all of them,
            even though in principle effect is for electron neutrino only. 
            This is in accordance to compatibility to ref [1].

REFERENCES:

[1] LIBRARY OF SM AND ANOMALOUS W W GAMMA COUPLINGS FOR
    THE E+ E- ---> F ANTI-F N GAMMA MONTE CARLO PROGRAMS.
    By A. Jacholkowska (Orsay, LAL), J. Kalinowski (Warsaw U.), 
    Z. Was (CERN & Cracow, INP). CERN-TH-99-120, IFT-99-06, May 1999. 6pp. 
    Published in Comput.Phys.Commun.124:238-242,2000 
    e-Print Archive: hep-ph/9905225 

[2] LIBRARY OF ANOMALOUS TAU TAU GAMMA COUPLINGS FOR TAU+ TAU- (N GAMMA) 
    MONTE CARLO PROGRAMS. By T. Paul, J. Swain (Northeastern U.), 
    Z. Was (CERN & Cracow, INP). 2000. 
    Published in Comput.Phys.Commun.124:243-246,2000