rounding errors?

[Dr. Georg Artus]

Hello,

ok, I have to confess a typing error. A 8 was missing in my Q. So the
conversion between units is ok within rounding errors.
Unfortunately this doesn't solve my problem. Your plot for the curve
intensity against wavelength shows a curve symmetrically around the
energy which should be selected by the monochromator. Mine still is way
off.
I have modified the input file slightly. It is attached. I would like to
send you my curves I produce with this file. Probrably you can solve the
problem at once. I have it in Excel. Can you use that? Or is
transmission via Fax better?

Georg


-- 
*********************************************

Dr. Georg Artus
Technische Universitaet Muenchen
Neue Forschungsneutronenquelle Garching
D-85747 Garching

Tel: +49 (0)89/289-14675
Fax: +49 (0)89/289-14666 or 12112
E-mail: gartus at ph.tum.de

-------------- next part --------------
DEFINE INSTRUMENT RESI()

DECLARE
%{
double energy = 81.804246;
/*
double d_energy = 8.2215;
*/
double d_energy = 31.804;
double l_energy = 50.000;
double h_energy = 90.000;
%}

TRACE

COMPONENT arm = Arm()
  AT (0,0,0) ABSOLUTE

COMPONENT source = Source_flat(
        radius = 0.1,
        dist = 5.0,
        xw = 0.04, yh = 0.09,
        E0 = energy,
        dE = d_energy)
  AT (0,0,0) RELATIVE arm

COMPONENT enbemon = E_monitor(
        xmin = -0.03,
        xmax =  0.03,
        ymin = -0.05,
        ymax =  0.05,
        Emin = l_energy,
        Emax = h_energy,
        nchan = 20,
        filename = "enbemon.out")
  AT (0,0,4.84) RELATIVE arm

COMPONENT arm45m = Arm()
  AT (0,0,5) RELATIVE arm ROTATED (0,45,0) RELATIVE arm

COMPONENT arm90m = Arm()
  AT (0,0,5) RELATIVE arm ROTATED (0,90,0) RELATIVE arm

COMPONENT mono = Monochromator(
        zmin = -0.03,
        zmax = 0.03,
        ymin = -0.05,
        ymax = 0.05,
        mosaich = 10,
        mosaicv = 10,
        r0 = 0.5,
        Q = 8.88577)  /* 45? at 1AA */
  AT (0,0,0) RELATIVE arm45m
/*
COMPONENT psdafmon = PSD_monitor(
        xmin = -0.03,
        xmax = 0.03,
        ymin = -0.06,
        ymax = 0.06,
        nx = 30,
        ny = 70,
        filename = "psdafmon.out")
  AT (0,0,0.1) RELATIVE arm90m
*/
COMPONENT enafmon = E_monitor(
        xmin = -0.03,
        xmax =  0.03,
        ymin = -0.06,
        ymax =  0.06,
        Emin = l_energy,
        Emax = h_energy,
        nchan = 20,
        filename = "enafmon.out")
  AT (0,0,0.1) RELATIVE arm90m

COMPONENT ensample = E_monitor(
        xmin = -0.005,
        xmax =  0.005,
        ymin = -0.005,
        ymax =  0.005,
        Emin = l_energy,
        Emax = h_energy,
        nchan = 20,
        filename = "ensample.out")
  AT (0,0,2.0) RELATIVE arm90m
/*
COMPONENT psdsamar = PSD_monitor(
        xmin = -0.03,
        xmax = 0.03,
        ymin = -0.07,
        ymax = 0.07,
        nx = 30,
        ny = 70,
        filename = "psdsamar.out")
  AT (0,0,2.0) RELATIVE arm90m
*/
END

-------------- next part --------------
A non-text attachment was scrubbed...
Name: energie.xls
Type: application/vnd.ms-excel
Size: 25600 bytes
Desc: not available
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/19990223/03f69571/attachment.xls>