Eric Hawker, 2/11/00 The beam monte carlo code that Heidi developed now has in it calls to LEPTO which simulates neutrino-nucleon interactions. The main code is called nuint.f, and the script that compiles and runs nuint.f is call rundis. You should edit rundis to point to the location of cernlib and lund on your local system. When running the code you need to specify 2 parameter files, a parameter file (nmu.inc) that controls what nuint.f does and a structure function file which Heidi's code uses for a separate cross section calculation. After everything is set up type rundis nuint nmu cteq4m this will use the nmu.inc and cteq4m.ini parameter files. In order for lepto to run smoothly I have had to make a 8.5 GeV cut on the total nuetrino momentum. Neutrinos with lower momentum cause lepto to freeze up. The "end" product particles of the hadron shower are output with their momentum and energy. The only output particles that I have seen are photons, muons, pions, kaons, protons, and neutrons. I have not yet tested electron neutrino interactions. I have replaced some of Heidi's original ntuple variables with lepto output. If you turn lepto off in the parameter file, then you should get Heidi's old ntuple variables back. Eric Hawker 2/16/00 I found a bug in the code, in nue mode the LEPTO events would be nue events but the incoming nuetrino momenta were that of the muon neutrino. This has been fixed now. In addintion near the end of the .log file produced there is now an entry that looks like XXXXXXXX cross section in pb estimate from MC events where the XXXXXXX is the cross section in pb estimated by LEPTO based on the events it generated. This cross section takes into account both the 1GeV Q**2 cut and the 8.5 GeV nuetrino momentum cut. Since the log files can be long the best way to see the cross section is to use tail -20 logfile.log ********************** parameter file nmu.inc ************************* * based on CJ2.0 with 600 instead of 688 straight section * HMS 11-16-99 * EAH added lepto parameters 2/1/00 parameter (straight = 60000) parameter (sin2tw = .2253) parameter (ndecays = 50000) parameter (output = .false.) ! write output record parameter (polarize = 999.) ! 1 -> spin along momentum vector parameter (qmu = -1.) ! muon charge (+1 or -1) parameter (npbeam = 50.) ! nominal beam momentum (GeV/c) parameter (dpbeam = 0.01) ! percent spread parameter (sigma_dx = 0.00073) ! Horizontal beam divergence (radians) parameter (sigma_dy = 0.00073) ! Vertical beam divergence (radians) parameter (sigma_x = 1.2333) ! Horizontal beam size (cm) parameter (sigma_y = 0.8833) ! Vertical beam size (cm) parameter (m_d = 1.0) ! detector thickness (g/cm^2) parameter (n_a = 6.0E+23) ! avagadro parameter (emass = 0.0005109990615) parameter (mumass = 0.105658389) parameter (twopi = 6.28318530717958648) parameter (n_mu = 1.E+20) c parameters for lepto parameter (leptoflag = .true.) ! generate lepto event c select type of incident neutrino parameter (lepin = 14) ! 14=numu, 12=nue, -14=numubar, -12=nuebar parameter (interaction = 2) ! 2=charged current, 3=neutral current parameter (nucleon = 1) ! nucleon hit by neutrino 1=proton, 2=neutron parameter (targa = 1.0) ! A of target nucleus parameter (targz = 1.0) ! Z of target nucleus c detector effiecencies parameter (shadcal = 0.80) !hadron calorimeter eff. parameter (semcal = 0.20) !EM calorimeter eff. parameter (smucal = 0.10) !muon energy measurement eff. parameter (sthnu = 1.0) !uncertainty in incident nuetrino angle in mrad parameter (sthmum = 3.0) !uncertainty in measured muon angle in mrad ************** ntuple list of beam properties and lepto output *********** ny(1) = float(idecay) ! event number ny(2) = pbeam ! muon beam energy ny(3) = qmu ! muon charge ny(4) = polar ! muon polarization ny(5) = ityp ! 1 for e, 2 for mu neutrino ny(6) = -qmu ! -1 for anti, 1 for neutrino c cm muon neutrino ny(7) = pnum_cms(4) ! neutrino energy cm ny(8) = pnum_cms(3)/pnum_cms(4) ! cos theta in cm c decay effects ny(9) = -straight*rndm(dummy13) ! decay position ny(10) = xprime ! beam angle in x ny(11) = yprime ! beam angle in y c lab muon neutrino ny(12) = pnum_lab(1) ! lab momentum ny(13) = pnum_lab(2) ny(14) = pnum_lab(3) ny(15) = pnum_lab(4) ny(16) = emax_num c ny(17) = nu_wt*det_wt/rescale !nue lum / gr/cm2 hydrogen ny(18) = eynum ! ybj for nu-e scatter ny(19) = ncenuml ny(20) = ncenumr ny(21) = ccenuml ny(22) = ccenumr ny(23) = rg32(dummy3)*sigma_x ny(24) = rg32(dummy3)*sigma_y ny(25) = thetaex ny(26) = thetaey ! will be theta for nu-e c ny(27) = det_wt*nu_wt/rescale ! lum for dis ny(28) = sigma_num*rescale ! lum*sigma for dis isoscalar ny(29) = i ! ith particle in lepto event ny(30) = p(i,1) ! x momentum of the ith particle ny(31) = p(i,2) ! y momentum of the ith particle ny(32) = p(i,3) ! z momentum of the ith particle ny(33) = p(i,4) ! energy of the ith particle ny(34) = p(i,5) ! rest mass of the ith particle ny(35) = x ! Bjorken-x of the lepto event ny(36) = y ! Bjorken-y of the lepto event ny(37) = w2 ! W**2 of the lepto event ny(38) = q2 ! Q**2 of the lepto event ny(39) = u ! energy transfered in the event (nu) ny(40) = smrq2 ! smeared Q**2 of the lepto event ny(41) = smrx ! smeared Bjorken-x of the lepto event ny(42) = smry ! smeared Bjorken-y of the lepto even