/*******************************************************************************
*         McXtrace instrument definition URL=http://www.mcxtrace.org
*
* Instrument: test (rename also the example and DEFINE lines below)
*
* %Identification
* Written by: Erik B Knudsen (erkn@fysik.dtu.dk)
* Date: Feb 1st 2013
* Origin: DTU Fysik 
* Release: McXtrace 1.1
* Version: 1.0
* %INSTRUMENT_SITE: Tests_sources
*
* Test instrument to show that the quasi-stcohastic source component works
*
* %Description
* This instrument is a unit test for the quasi-stochastic source component.
*
* %Example: Test_Source_quasi.instr SRC=0 Detector: psd_I=1.21847E-18
* %Example: Test_Source_quasi.instr SRC=1 Detector: psd_I=1.21847E-18
*
* %Parameters
* SRC: [ ] Integer parameter picks a source model. 0 normal Source_div, Quasi-stochastic 
*
* %End
*******************************************************************************/

/* Change name of instrument and input parameters with default values */
DEFINE INSTRUMENT Test_Source_quasi(SRC=1)

/* The DECLARE section allows us to declare variables or  small      */
/* functions in C syntax. These may be used in the whole instrument. */
DECLARE
%{
%}

/* The INITIALIZE section is executed when the simulation starts     */
/* (C code). You may use them as component parameter values.         */
INITIALIZE
%{
%}

/* Here comes the TRACE section, where the actual      */
/* instrument is defined as a sequence of components.  */
TRACE

/* The Arm() class component defines reference points and orientations  */
/* in 3D space. Every component instance must have a unique name. Here, */
/* Origin is used. This Arm() component is set to define the origin of  */
/* our global coordinate system (AT (0,0,0) ABSOLUTE). It may be used   */
/* for further RELATIVE reference, Other useful keywords are : ROTATED  */
/* EXTEND GROUP PREVIOUS. Also think about adding a neutron source !    */
/* Progress_bar is an Arm displaying simulation progress.               */
COMPONENT Origin = Progress_bar()
  AT (0,0,0) ABSOLUTE

COMPONENT sd = Source_div(
    xwidth=1e-3, yheight=1e-3, focus_aw=1e-6, focus_ah=1e-6, lambda0=2.5, dlambda=0.2, gauss=0
)
WHEN (SRC==0) AT (0,0,0) RELATIVE Origin

COMPONENT sdq = Source_div_quasi(
    xwidth=1e-3, yheight=1e-3, focus_aw=1e-6, focus_ah=1e-6, lambda0=2.5, dlambda=0.2, gauss=0
)
WHEN (SRC==1) AT (0,0,0) RELATIVE Origin




COMPONENT lmon= L_monitor(
   restore_neutron=1, xwidth=0.02, yheight=0.02,Lmin=2,Lmax=3.2,filename="lmon.dat")
AT(0,0,1) RELATIVE Origin

COMPONENT emon= E_monitor(
   restore_neutron=1, xwidth=0.02, yheight=0.02,Emin=5,Emax=20,filename="emon.dat")
AT(0,0,1) RELATIVE Origin

COMPONENT psd= PSD_monitor(
   restore_neutron=1, xwidth=0.0011, yheight=0.0011,nx=101, ny=101,filename="psdmon.dat")
AT(0,0,4) RELATIVE Origin



/* This section is executed when the simulation ends (C code). Other    */
/* optional sections are : SAVE                                         */
FINALLY
%{
%}
/* The END token marks the instrument definition end */
END

