benchmark calc. listing

From: seki@neutrino.kek.jp
Date: Mon Oct 15 2001 - 08:16:30 JST

  • Next message: Makoto Sakuda: "Re: benchmark calc. listing" 

    Date: Mon, 15 Oct 2001 08:16:30 +0900 (JST)
From: seki@neutrino.kek.jp
To: opc.nuint01@neutrino.kek.jp
Subject: benchmark calc. listing
Message-ID: <Pine.LNX.3.95LJ1.1b3.1011015080810.8685A-100000@neutrino.kek.jp>

Folks,
The following is my final version of the benchmark calculation listing, 
with the inclusion of the comments from Hayato san, Dave, and Sakuda san.

Please check it for the last time and give me your comments and
suggestions asap.  If I do not hear from you, I will give this to Ishida
san for posting in the NuInt01 HP.  Note that there will be a figure
added, for specifying \theta_\mu and \theta_p.  

Best wishes,
Seki

p.s. RIKEN has decided to cancel its members' attendance at the Hawaii
US/Japan nuclear physics meeting next week.  Is KEK going to do the same?

*****************************************

Benchmark MC and Theoretical Calculations of $\nu_\mu$ Reactions

  *WARNING:
       1)  You are not required to do all benchmark calculations listed
           below, but please try to do as many as possible.  Any number 
           of the calculations and any portion of the calculations are
           acceptable.  For example, showing the quasi-elastic contibution 
           alone is acceptable for the benchmark calculation 1, if you are 
           not set up to calculate the entire inclusive cross section.
       2)  Use the parameter values specified below, and do NOT use your
           favorite values even if you consider yours to be better.
           The objective of the benchmark calculations is a calibration 
           of the instruments, your MC codes and theoretical methods, and 
           the use of the common parameter values is thus vital. 

 *GUIDELINES:
       1)  Make the binning size as small as possible.
       2)  Draw straight plots (neither semi-log, nor log plots, 
           unless specified).
       3)  Use the neutrino-nucleon cross section as specified below.



 *BENCHMARK CALCULATION LISTING:
    1) 16O(\nu_\mu,\mu^-) Inclusive.
       (Relevant to Panel Discussions 1 and 2)

          MAIN:(Three figures) 
               {d^2\sigma}/{4\pi d(\cos\theta_\mu)d E_\mu}
               at theta_\mu = 30{}^o for E_\nu = 0.5, 1.0, and 5.0 GeV 
               as a function of $E_\nu - E_\mu$ from 0 to 500 MeV,
               also showing a curve for the quasi-elastic contribution 
               in each figure.  
               (Please draw semi-log plots for 5.0 GeV.)
          SUPPLEMENTAL: 
            a. (Three figures) 
               {d\sigma}/{d E_\mu} for E_\nu = 0.5, 1.0, and 5.0 GeV
               as a function of E_\mu from 0 to 1.0 GeV,
               also showing a curve for the quasi-elastic contribution 
               in each figure.  
            b. (Three figures) 
               {d\sigma}/{4\pi d(\cos\theta_\mu)}
               for E_\nu = 1.0 and 2.0 GeV
               as a function of cos(\theta_\mu)  from 1 to -1,
               also showing a curve for the quasi-elastic contribution 
               in each figure.   


    2) 16O(\nu_\mu,\mu^-p) Semi-inclusive.
       (Relevant to Panel Discussions 1 and 2)

            a. (Two figures)
               {d\sigma}/{d^3p_p} (\vec{k}_\mu integrated)
               at \theta_p =60^o for E_\nu = 1.0 and 2.0 GeV 
               as a function of |\vec{p}_p| from 0 to 500 MeV/c, 
               also showing a curve for the quasi-elastic contribution 
               in each figure.  
            b. (Two figures) 
               {d\sigma}/{d^3p_p} (vec{k}_\mu integrated)
               at |\vec{p}_p|= 500 MeV for E_\nu = 1.0 and 2.0 GeV 
               as a function of \theta_p from 0 to 90^o,
               also showing a curve for the quasi-elastic contribution 
               in each figure.  


    3) 16O and/or 56Fe Pion production 
       (Relevant to Panel Discussion 3) 
       (The number of the figures indicated below 
       is for each of 16O and 56Fe)
          MAIN: 
            a. [Two figures = 1 fig.(for nucleus) + 1 fig.(for n)] 
               Pion multiplicity 
               as a function of E_\nu from 0.5 GeV to 5.0 GeV 
               for \nu_\mu - nucleus and -neutron reactions.
            b. [Two figures = 1 fig.(for nucleus) + 1 fig.(for n)] 
               Average pion energy (per pion)
               as a function of E_\nu from 0.5 GeV to 5.0 GeV 
               for \nu_\mu - nucleus and -neutron reactions.
          SUPPLEMENTAL: 
            a. [Eight figures = 2 figs. (for nucleus and n)
                              x 2 figs. (for the directiuon and magnitude)
                              x 2 figs. (30 MeV/c and 50 MeV/c cut-off)] 
               Hadron emission momentm vector (direction and magnitude) 
               as a function of E_\nu from 0.5 GeV to 5.0 GeV 
               for \nu_\mu - nucleus and neutron. 
               (The vector is defined to be the weighted average of 
                all emitted hadron (nucleons and mesons) momenta 
                above the 30 MeV/c and 50 MeV/c cut-offs.)
            b. (Twelve numbers) 
               The above MAIN a. and b. and SUPPLEMENTAL a. at 10 GeV/c. 



  *\nu - Neutron cross section parameters

     1)  G_F cos(\theta_c) = 1.14 x 10^{-5} GeV^{-2}

     2)  M_V^2             = 0.71           GeV^2  

     3)  g_A               = -1.26

     4)  M_A               = 1.03 GeV.

     5)  1 + \mu_p -\mu_n  = 4.71 (nuclear magneton)

     Also set F_p(q) = 0.

     The neutrino-neutron cross section {d\sigma}/dt is given by, e.g.,
     C. H. Llewellyn Smith, Physics Report 3, no. 5 (1972) pp300-,
     Subsection 3.3 -A & B.

    This archive was generated by hypermail 2.1.0 : Mon Oct 15 2001 - 08:16:32 JST