From diwan@bnlku1.phy.bnl.gov Tue Mar 19 12:15:38 2002
Status: RO
X-VM-v5-Data: ([nil nil nil nil nil nil nil nil nil]
	["1879" "Tuesday" "19" "March" "2002" "12:15:37" "-0500" "Milind Diwan" "diwan@bnlku1.phy.bnl.gov" nil "91" "spectra" "^Subject:" nil nil "3" nil "spectra" nil nil nil nil nil nil]
	nil)
Return-path: <diwan@bnlku1.phy.bnl.gov>
Envelope-to: bviren@minos.phy.bnl.gov
Received: from smtpgw.sec.bnl.local
	([192.168.1.21] helo=smtpgw.bnl.gov ident=mail)
	by minos with esmtp (Exim 3.34 #1 (Debian))
	id 16nND0-0002Cr-00
	for <bviren@minos.phy.bnl.gov>; Tue, 19 Mar 2002 12:15:38 -0500
Received: from bnl.gov ([130.199.128.163])
	by smtpgw.bnl.gov with esmtp (Exim 3.34 #1 )
	id 16nNCz-00052V-00; Tue, 19 Mar 2002 12:15:37 -0500
Received: (from daemon@localhost)
	by bnl.gov (8.9.2/8.9.2) id MAA21891;
	Tue, 19 Mar 2002 12:15:37 -0500 (EST)
Received: from bnlku1 (bnlku1.phy.bnl.gov [130.199.36.3])
	by bnl.gov (8.9.2/8.9.2) with ESMTP id MAA21886;
	Tue, 19 Mar 2002 12:15:37 -0500 (EST)
Received: from diwan by bnlku1 with local (Exim 3.34 #1 (Debian))
	id 16nNCz-0005GN-00; Tue, 19 Mar 2002 12:15:37 -0500
Message-ID: <15511.29241.155610.992737@bnlku1.phy.bnl.gov>
Subject: spectra
From: Milind Diwan <diwan@bnlku1.phy.bnl.gov>
To: bv@bnl.gov, diwan@bnl.gov
Date: Tue, 19 Mar 2002 12:15:37 -0500
X-PH: V4.4@bnl.gov
X-Mailer: VM 6.92 under 21.4 (patch 4) "Artificial Intelligence" XEmacs Lucid
MIME-Version: 1.0
Content-Transfer-Encoding: 7bit
Content-Type: text/plain; charset=us-ascii



Each file and the spectra it contains


1) nugenwe.f  
	      
	      Old wide band beam at 0 degree. 

	      se(70,1)   is for neutrinos 

	      se(70,2)   is for anti-neutrinos 

	      dele = 0.1 GeV

	      units are arb. 


2) nugenw.f 

	    0degree wide band spectrum from 82 and 83 data

	    se(70,1) is for neutrinos	
	    secon(70,1) is contamination of anti-nu in neutrino spec

	    se(70,2) is for anti-nus 
	    secon(70,2) is contamination of nus in anti-nus. 


3) nugenne.f 

	     narrow band beam 

	     se(40,1) is for neutrinos 
	     
	     se(40,2) is for anti-neutrinos 

	     dele = 0.25 GeV 	  

4) nugennb.f  

	      reads in neutrinos events from some file. 
	      This is obsolete. I have no idea how it works 	      

5) nugenn.f  

	     a "square" narrow band beam. just input enumin and enumax

6) nugenw_test.f  

   This is the same as nugenw, except it only does neutrinos 

   se_00(100) is neutrino flux at 0 degrees with units of
    particles/m**2/POT at 1 km.    dele = 0.1 GeV 

   se_15(100) is neutrino flux at 1.5 degrees with same units 
   but with dele = 0.15 GeV. 


7) nugennb_test.f  

   reads in neutrino rays from a file and intersects them with 
   a cylinder. Produces the correct distributions, no matter where
   you put the cylindrical detector. 

   This is what I actually used for the E889 simulations. 

   I had a file containing all types of neutrinos for both 
   polarities and then I could just request a particular 
   reaction. This program fetched the correct neutrinos from
   the file.

   This neutrino file should still be soemwhere on disk. 
   


There are two makefiles 

nusim_main:  nugenn, nugennb, nugenne, nugenw, nugenwe


nusim_test: nugenn, nugennb_test, nugenne, nugenw_test, nugenwe 

	    I used the test program for the e889 simulations. 


Sorry took so long 

regards,

Milind