Compare commits
1 Commits
| Author | SHA1 | Date | |
|---|---|---|---|
 glavic_a
|
1415a39ff0 |
@@ -8,7 +8,6 @@ old
|
||||
*.backup
|
||||
.project
|
||||
.settings
|
||||
.idea
|
||||
*.pyc
|
||||
*.pyo
|
||||
hosts
|
||||
|
||||
@@ -1,9 +0,0 @@
|
||||
#!/bin/bash
|
||||
# compress McStas hdf5 files to increase performance and save disk space
|
||||
|
||||
for fi in $*
|
||||
do
|
||||
echo $fi
|
||||
h5repack -i $fi -o $fi.c -f /entry1/data/detector_list_p_x_y_t_L_sx_sy_sz/events:GZIP=5 -l /entry1/data/detector_list_p_x_y_t_L_sx_sy_sz/events:CHUNK=3072x8
|
||||
mv $fi.c $fi
|
||||
done
|
||||
@@ -8,7 +8,7 @@ import sys, os
|
||||
from numpy import *
|
||||
seterr(all='ignore')
|
||||
import estia_help as eh
|
||||
from aglib import mcstas_reader as mr
|
||||
import mcstas_reader as mr
|
||||
|
||||
|
||||
# Scaling factor of source monitor after normalizing by area to get to brilliance.
|
||||
|
||||
@@ -0,0 +1,422 @@
|
||||
#-*- coding: utf-8 -*-
|
||||
'''
|
||||
Simple support library for reading, analyzing and plotting of McStas results
|
||||
from the Estia instrument simulations.
|
||||
Meant to be used from IPython Notebook or QtConsole, but can also be run stand alone.
|
||||
'''
|
||||
|
||||
import os, sys
|
||||
from numpy import *
|
||||
|
||||
try:
|
||||
import h5py
|
||||
except ImportError:
|
||||
print "h5py not found, modern NeXuS format will not be readable."
|
||||
|
||||
try:
|
||||
from IPython import display #@UnusedImport
|
||||
from IPython.core.pylabtools import print_figure
|
||||
from matplotlib.figure import Figure
|
||||
from matplotlib.backends.backend_agg import FigureCanvasAgg
|
||||
from matplotlib.colors import LogNorm
|
||||
except ImportError:
|
||||
# IPython and/or matplotlib not available, no interactive functionalities
|
||||
display=None
|
||||
|
||||
MAX_EVTS_BATCH=50000
|
||||
|
||||
class McSim(object):
|
||||
'''
|
||||
Object representing complete simulation from McStas.
|
||||
Supports either the old style McStas output with separate ASCII files
|
||||
or HDF5 single file output.
|
||||
|
||||
Different monitors can be accessed as keys like in a dictionary. The data is only loaded when the monitor is first accessed.
|
||||
'''
|
||||
|
||||
def __init__(self, path):
|
||||
'''
|
||||
Initialize the simulation. path should be the file name to either the NeXuS or the mccode.sim file.
|
||||
'''
|
||||
self._data={}
|
||||
if path.endswith('.h5'):
|
||||
self._init_hdf(path)
|
||||
elif path.endswith('.sim'):
|
||||
self._init_old(path)
|
||||
else:
|
||||
if os.path.exists(os.path.join(path,'mccode.h5')):
|
||||
self._init_hdf(os.path.join(path,'mccode.h5'))
|
||||
elif os.path.exists(os.path.join(path, 'mccode.sim')):
|
||||
self._init_old(os.path.join(path, 'mccode.sim'))
|
||||
else:
|
||||
raise IOError, "Can't locate mccode.h5 of mccode.sim file in %s"%path
|
||||
|
||||
def _init_hdf(self, path):
|
||||
self.hdf=h5py.File(path, 'r')
|
||||
self.data_loader=DataLoaderHDF(self.hdf)
|
||||
self.info=self.data_loader.info
|
||||
|
||||
def _init_old(self, path):
|
||||
self.info=HeaderFile(path)
|
||||
self.data_loader=DataLoaderOld(self.info, os.path.dirname(path))
|
||||
|
||||
def keys(self):
|
||||
return self.info['data'].keys()
|
||||
|
||||
def monitors(self):
|
||||
output=[]
|
||||
for val in self.info['data'].values():
|
||||
xy=(val['xvar'], val['yvar'])
|
||||
if not xy in output:
|
||||
output.append(xy)
|
||||
output.sort()
|
||||
return output
|
||||
|
||||
def plot(self, monitors=None):
|
||||
graphs={}
|
||||
for key, val in self.info['data'].items():
|
||||
xy=(val['xvar'], val['yvar'])
|
||||
if monitors is not None and xy!=monitors:
|
||||
continue
|
||||
data=self[key]
|
||||
if xy in graphs:
|
||||
graphs[xy].append(data)
|
||||
else:
|
||||
graphs[xy]=[data]
|
||||
for xy, datasets in sorted(graphs.items()):
|
||||
cols=min(len(datasets), 3)
|
||||
rows=len(datasets)/3+1
|
||||
|
||||
fig=Figure(figsize=(12, 5*rows), dpi=300, facecolor='#FFFFFF')
|
||||
FigureCanvasAgg(fig)
|
||||
|
||||
for i, data in enumerate(datasets):
|
||||
ax=fig.add_subplot(rows, cols, i+1)
|
||||
data.plot(ax=ax)
|
||||
graphs[xy]=fig
|
||||
if monitors is None:
|
||||
return graphs
|
||||
else:
|
||||
return fig
|
||||
|
||||
def __getitem__(self, item):
|
||||
if item in self._data:
|
||||
return self._data[item]
|
||||
elif item in self.keys():
|
||||
data=self.data_loader.load_item(item)
|
||||
self._data[item]=data
|
||||
return data
|
||||
else:
|
||||
raise KeyError, "Can't find dataset %s"%item
|
||||
|
||||
class HeaderFile(object):
|
||||
'''
|
||||
Analyze McStas mccode.sim header files.
|
||||
'''
|
||||
_data=None
|
||||
|
||||
@property
|
||||
def data(self):
|
||||
return self._data['data']
|
||||
|
||||
def __init__(self, path):
|
||||
self._data={}
|
||||
data=open(path, 'r').read()
|
||||
data_lines=data.splitlines()
|
||||
if not data.startswith('McStas simulation description file'):
|
||||
raise IOError, 'Not a valid McStas description file.'
|
||||
self._data['start_time']=data_lines[1].split(':', 1)[1].strip()
|
||||
self._data['program_name']=data_lines[2].split(':', 1)[1].strip()
|
||||
self._data['data']=self.get_data(data)
|
||||
|
||||
def get_data(self, data):
|
||||
output={}
|
||||
start_idx=0
|
||||
while start_idx<len(data):
|
||||
start_idx=data.find('begin data', start_idx)
|
||||
if start_idx==-1:
|
||||
break
|
||||
end_idx=data.find('end data', start_idx)
|
||||
block=data[start_idx:end_idx].splitlines()[1:]
|
||||
block_info={}
|
||||
for bline in block:
|
||||
item, value=bline.strip().split(':', 1)
|
||||
block_info[item]=value.strip()
|
||||
if item=='component':
|
||||
output[value.strip()]=block_info
|
||||
|
||||
start_idx=end_idx
|
||||
return output
|
||||
|
||||
def __getitem__(self, item):
|
||||
return self._data[item]
|
||||
|
||||
class DataLoaderOld(object):
|
||||
'''
|
||||
Load and analyze a old style McStas format with a simulation and a set of data text files.
|
||||
'''
|
||||
def __init__(self, info, root):
|
||||
self.info=info
|
||||
self.root=root
|
||||
|
||||
def load_item(self, item):
|
||||
item_info=self.info['data'][item]
|
||||
if not item_info['type'].startswith('array_2d'):
|
||||
return self.load_item_1d(item)
|
||||
|
||||
fname=os.path.join(self.root, item_info['filename'])
|
||||
x_col=item_info['xvar']
|
||||
y_col=item_info['yvar']
|
||||
if x_col=='Li' and y_col=='p': # Detector_nD
|
||||
cols=item_info['variables'].split()
|
||||
data=loadtxt(fname, dtype={'names': cols, 'formats': ['f4']*len(cols)})
|
||||
return TofData(data, item_info)
|
||||
else:
|
||||
raw=loadtxt(fname)
|
||||
data=raw[:len(raw)/3]
|
||||
return Dataset(data, item_info)
|
||||
|
||||
def load_item_1d(self, item):
|
||||
item_info=self.info['data'][item]
|
||||
fname=os.path.join(self.root, item_info['filename'])
|
||||
raw=loadtxt(fname).T
|
||||
data=raw[1]
|
||||
errors=raw[2]
|
||||
return Dataset1D(data, errors, item_info)
|
||||
|
||||
class DataLoaderHDF(object):
|
||||
'''
|
||||
Load and analyze a new style NeXuS file format.
|
||||
'''
|
||||
def __init__(self, hdf):
|
||||
self.hdf=hdf['entry1']
|
||||
self.info={}
|
||||
self.info['data']={}
|
||||
for item in self.hdf['data'].keys():
|
||||
node=self.hdf['data/'+item]
|
||||
info={}
|
||||
for key, value in node.attrs.items():
|
||||
info[key.strip()]=value.strip()
|
||||
if info['filename'][-4:]=='.dat' or '_list.' in info['filename']:
|
||||
self.info['data'][info['component']]=info
|
||||
else:
|
||||
self.info['data'][info['filename']]=info
|
||||
info['datapath']='data/'+item
|
||||
|
||||
def load_item(self, item):
|
||||
item_info=self.info['data'][item]
|
||||
if not item_info['type'].startswith('array_2d'):
|
||||
return self.load_item_1d(item)
|
||||
|
||||
node=self.hdf[item_info['datapath']]
|
||||
x_col=item_info['xvar']
|
||||
y_col=item_info['yvar']
|
||||
if x_col=='Li' and y_col=='p': # Detector_nD
|
||||
cols=item_info['variables'].split()
|
||||
evds=node['events']
|
||||
if len(evds)<=MAX_EVTS_BATCH:
|
||||
data=evds.value.astype(float32).view(
|
||||
dtype={'names': cols, 'formats': ['f4']*len(cols)}).flatten()
|
||||
else:
|
||||
ds=[]
|
||||
sys.stdout.write('Reading large dataset:\n')
|
||||
for i in range(len(evds)//MAX_EVTS_BATCH+1):
|
||||
sys.stdout.write('\r%i/%i'%(i*MAX_EVTS_BATCH, len(evds)))
|
||||
sys.stdout.flush()
|
||||
ds.append(evds[i*MAX_EVTS_BATCH:(i+1)*MAX_EVTS_BATCH])
|
||||
sys.stdout.write('\r%i/%i\n'%(len(evds), len(evds)))
|
||||
data=vstack(ds).astype(float32).view(
|
||||
dtype={'names': cols, 'formats': ['f4']*len(cols)}).flatten()
|
||||
return TofData(data, item_info)
|
||||
else:
|
||||
data=node['data'].value.T
|
||||
return Dataset(data, item_info)
|
||||
|
||||
def load_item_1d(self, item):
|
||||
item_info=self.info['data'][item]
|
||||
node=self.hdf[item_info['datapath']]
|
||||
data=node['data'].value
|
||||
errors=node['errors'].value
|
||||
return Dataset1D(data, errors, item_info)
|
||||
|
||||
class Dataset1D(object):
|
||||
'''
|
||||
Representation of a standard McStas dataset of one variable.
|
||||
'''
|
||||
def __init__(self, data, errors, info):
|
||||
self.data=data
|
||||
self.errors=errors
|
||||
self.info=info
|
||||
|
||||
def plot(self, log=False, ax=None):
|
||||
if ax is None:
|
||||
import pylab
|
||||
ax=pylab.gca()
|
||||
|
||||
limits=map(float, self.info['xlimits'].split())
|
||||
x=linspace(limits[0], limits[1], len(self.data))
|
||||
|
||||
ax.errorbar(x, self.data, yerr=self.errors)
|
||||
if log:
|
||||
ax.set_yscale('log')
|
||||
else:
|
||||
ax.set_yscale('linear')
|
||||
ax.set_xlabel(self.info['xlabel'])
|
||||
ax.set_ylabel(self.info['ylabel'])
|
||||
ax.set_title(self.info['component'])
|
||||
|
||||
def _repr_png_(self):
|
||||
'''
|
||||
Image representation form ipython console/notebook.
|
||||
'''
|
||||
fig=Figure(figsize=(8, 5), dpi=300, facecolor='#FFFFFF')
|
||||
FigureCanvasAgg(fig)
|
||||
ax=fig.add_subplot(111)
|
||||
|
||||
self.plot(ax=ax)
|
||||
|
||||
return print_figure(fig, dpi=72)
|
||||
|
||||
|
||||
class Dataset(object):
|
||||
'''
|
||||
Representation of a standard McStas dataset with 2D view.
|
||||
'''
|
||||
def __init__(self, data, info):
|
||||
self.data=data
|
||||
self.info=info
|
||||
|
||||
def plot(self, log=False, ax=None):
|
||||
if ax is None:
|
||||
import pylab
|
||||
ax=pylab.gca()
|
||||
|
||||
limits=map(float, self.info['xylimits'].split())
|
||||
|
||||
if log:
|
||||
ax.imshow(self.data, origin='lower', extent=limits, aspect='auto', norm=LogNorm())
|
||||
else:
|
||||
ax.imshow(self.data, origin='lower', extent=limits, aspect='auto')
|
||||
ax.set_xlabel(self.info['xlabel'])
|
||||
ax.set_ylabel(self.info['ylabel'])
|
||||
ax.set_title(self.info['component'])
|
||||
|
||||
def _repr_png_(self):
|
||||
'''
|
||||
Image representation form ipython console/notebook.
|
||||
'''
|
||||
fig=Figure(figsize=(8, 5), dpi=300, facecolor='#FFFFFF')
|
||||
FigureCanvasAgg(fig)
|
||||
ax=fig.add_subplot(111)
|
||||
|
||||
self.plot(ax=ax)
|
||||
|
||||
return print_figure(fig, dpi=72)
|
||||
|
||||
class TofData(Dataset):
|
||||
'''
|
||||
Representation of a dataset collected with Monitor_nD
|
||||
'''
|
||||
|
||||
def project1d(self, col, bins=50, fltr=None, newcols=None, norm=None):
|
||||
'''
|
||||
Generate binned data for arbitrary binning and columns.
|
||||
|
||||
The fltr argument can be a string with a filtering condition on the available columns.
|
||||
Syntax for this filtering follows numpy array convetions like (x>5)&(abs(L)<0.3) would
|
||||
be a valid statement.
|
||||
|
||||
newcols can be a list of (name, code) tuples, that calculate new columns from
|
||||
existing ones, like newcols=[('q', '4*pi/L*sin(theta)')].
|
||||
'''
|
||||
columns=dict([(coli, self.data[coli]) for coli in self.data.dtype.names])
|
||||
if newcols is not None:
|
||||
for name, code in newcols:
|
||||
columns[name]=eval(code, globals(), columns)
|
||||
|
||||
if norm is None:
|
||||
w=self.data['p']
|
||||
else:
|
||||
w=eval(norm+'*p', globals(), columns)
|
||||
|
||||
if fltr is None:
|
||||
I, x=histogram(columns[col], bins=bins, weights=w)
|
||||
else:
|
||||
if isinstance(fltr, basestring):
|
||||
fltr=eval(fltr, globals(), columns)
|
||||
I, x=histogram(columns[col][fltr], bins=bins, weights=w[fltr])
|
||||
return x, I
|
||||
|
||||
def project2d(self, xcol, ycol, bins=50, fltr=None, newcols=None):
|
||||
'''
|
||||
Generate 2D binned data for arbitrary binning and columns.
|
||||
|
||||
The fltr argument can be a string with a filtering condition on the available columns.
|
||||
Syntax for this filtering follows numpy array convetions like (x>5)&(abs(L)<0.3) would
|
||||
be a valid statement.
|
||||
|
||||
newcols can be a list of (name, code) tuples, that calculate new columns from
|
||||
existing ones, like newcols=[('q', '4*pi/L*sin(theta)')].
|
||||
'''
|
||||
columns=dict([(coli, self.data[coli]) for coli in self.data.dtype.names])
|
||||
if newcols is not None:
|
||||
for name, code in newcols:
|
||||
columns[name]=eval(code, globals(), columns)
|
||||
|
||||
if fltr is None:
|
||||
I, y, x=histogram2d(columns[ycol], columns[xcol],
|
||||
bins=bins, weights=self.data['p'])
|
||||
else:
|
||||
if isinstance(fltr, basestring):
|
||||
fltr=eval(fltr, globals(), columns)
|
||||
I, y, x=histogram2d(columns[ycol][fltr], columns[xcol][fltr],
|
||||
bins=bins, weights=columns['p'][fltr])
|
||||
return x, y, I
|
||||
|
||||
def plot(self, xcol='x', ycol='y', log=False, ax=None, bins=50, fltr=None, newcols=None):
|
||||
if ax is None:
|
||||
import pylab
|
||||
ax=pylab.gca()
|
||||
|
||||
x, y, I=self.project2d(xcol, ycol, bins=bins, fltr=fltr, newcols=newcols)
|
||||
|
||||
if log:
|
||||
ax.pcolormesh(x, y, I, norm=LogNorm())
|
||||
else:
|
||||
ax.pcolormesh(x, y, I)
|
||||
|
||||
col_names=self.info['title'].split()
|
||||
cols=self.info['variables'].split()
|
||||
|
||||
try:
|
||||
ax.set_xlabel(xcol+'-'+col_names[cols.index(xcol)])
|
||||
except ValueError:
|
||||
ax.set_xlabel(xcol)
|
||||
try:
|
||||
ax.set_ylabel(ycol+'-'+col_names[cols.index(ycol)])
|
||||
except ValueError:
|
||||
ax.set_ylabel(ycol)
|
||||
ax.set_title(self.info['component'])
|
||||
|
||||
|
||||
def plot1d(self, col='x', log=False, ax=None, bins=50, fltr=None, newcols=None):
|
||||
if ax is None:
|
||||
import pylab
|
||||
ax=pylab.gca()
|
||||
|
||||
x, I=self.project1d(col, bins=bins, fltr=fltr, newcols=newcols)
|
||||
|
||||
if log:
|
||||
ax.semilogy((x[:-1]+x[1:])/2., I)
|
||||
else:
|
||||
ax.plot((x[:-1]+x[1:])/2., I)
|
||||
col_names=self.info['title'].split()
|
||||
cols=self.info['variables'].split()
|
||||
|
||||
try:
|
||||
ax.set_xlabel(col+'-'+col_names[cols.index(col)])
|
||||
except ValueError:
|
||||
ax.set_xlabel(col)
|
||||
ax.set_ylabel('Intensity')
|
||||
ax.set_title(self.info['component'])
|
||||
|
||||
@@ -53,8 +53,7 @@ DEFINE INSTRUMENT ESS_reflectometer_Estia
|
||||
int operationmode = 0, double over_illumination = 0.0, double theta_resolution = 0.04,
|
||||
double lambda_min = 3.75, double lambda_start = 3.0, double lambda_end = 12.0,
|
||||
int enable_chopper = 0, int enable_gravity=0, int enable_windows=1,
|
||||
int enable_polarizer = 0, int enable_analyzer = 0,
|
||||
double pol1_start=0.3, double pol1_angle=1.66, double pol2_start=0.3, double pol2_angle=1.66,
|
||||
int enable_polarizer = 0, int enable_analyzer = 0,
|
||||
double source_power = 2,
|
||||
double selene1_foot1y =0.0, double selene1_foot2y = 0.0
|
||||
)
|
||||
@@ -105,20 +104,12 @@ double velocity_max ; // m/s neutron velocity of lambda_min
|
||||
|
||||
|
||||
double analyzer_max_height = 0.01; // Maximum sample height to be covered by te polarization analyzers
|
||||
double analyzer_flipper_start = 0.4; // Distance from sample to analyzer flipper
|
||||
|
||||
|
||||
double analyzer1_start = 0.65; // Distance to sample to start the first analyzer
|
||||
double analyzer1_length = 1.3; // length of first analyzer
|
||||
double analyzer2_start = 0.7; // Distance to sample to start the first analyzer
|
||||
double analyzer2_length = 1.6; // length of first analyzer
|
||||
|
||||
double Theta1_analyzer1, Theta1_analyzer2, Theta2_analyzer1, Theta2_analyzer2, dist_ana_vfocus; // quantities calculated out of values above and 1.5 degree covered divergence
|
||||
int p_int;
|
||||
%}
|
||||
|
||||
USERVARS %{
|
||||
int p_int; // a flag that gets incremented if a polarizer mirror scatters
|
||||
%}
|
||||
|
||||
INITIALIZE
|
||||
@@ -127,11 +118,11 @@ INITIALIZE
|
||||
velocity_max= 3.956034E3 / lambda_min; // h/m_n over lambda - ((3.9..e-7m^2/s))/(1e-10m)
|
||||
chopper_phase=360.0*(chopper_pos*chopper_freq/velocity_max+(pulse_zero-opening_time)*chopper_freq);
|
||||
chopper_open = 98.0;
|
||||
p_int=0;
|
||||
|
||||
|
||||
/* print out some calculated parameter for checking purposes */
|
||||
printf(" Chopper phase = %.1f deg\n", chopper_phase);
|
||||
//printf(" Chopper open = %.1f deg\n", chopper_open);
|
||||
printf(" Chopper open = %.1f deg\n", chopper_open);
|
||||
|
||||
|
||||
dist_ana_vfocus=analyzer_max_height/2.0/tan(1.5*PI/360.0); // the virtual focus point infront of the actual sample focus where the beams furthest out meet
|
||||
@@ -144,8 +135,8 @@ Theta2_analyzer2=atan((dist_ana_vfocus+analyzer2_start+analyzer2_length)/dist_an
|
||||
selene1_center=(selene1_foot1+selene1_foot2)/2.;
|
||||
selene1_shift=(selene1_foot1y+selene1_foot2y)/1000.0;
|
||||
selene1_rot=atan((selene1_foot2y-selene1_foot1y)/(selene1_foot2-selene1_foot1)/1000.0)*180.0/PI;
|
||||
//printf(" Selene 1 rotation = %.4f deg\n", selene1_rot);
|
||||
//printf(" Selene 1 shift = %.1f mm\n", selene1_shift*1e3);
|
||||
printf(" Selene 1 rotation = %.4f deg\n", selene1_rot);
|
||||
printf(" Selene 1 shift = %.1f mm\n", selene1_shift*1e3);
|
||||
|
||||
%}
|
||||
TRACE
|
||||
@@ -211,7 +202,7 @@ COMPONENT moderator = ESS_butterfly(
|
||||
sector = "E", beamline = 2, yheight = 0.03, cold_frac = 0.9,
|
||||
focus_xw = E02_01_01_Cu_in_xmax-E02_01_01_Cu_in_xmin+0.002,
|
||||
focus_yh = E02_01_01_Cu_in_yheight+0.002,
|
||||
target_index=3, Lmin = lambda_start, Lmax = lambda_end,
|
||||
target_index=4, Lmin = lambda_start, Lmax = lambda_end,
|
||||
n_pulses = 1+enable_chopper, acc_power=source_power)
|
||||
AT (0, 0, 0) RELATIVE origin
|
||||
|
||||
@@ -243,15 +234,6 @@ COMPONENT chopper = DiskChopper(radius=chopper_diameter/2.0, yheight=0.02,
|
||||
AT (0, 0, chopper_pos-2*NBOA_c) RELATIVE arm_virtual_source_beam
|
||||
|
||||
|
||||
COMPONENT vs_divergence_h = DivPos_monitor(nb=21, ndiv=41, filename="vs_hordiv",
|
||||
xmin=-0.02, xmax=0.02, ymin=-0.02, ymax=0.02, maxdiv=2.0,
|
||||
restore_neutron=1)
|
||||
AT (0, 0, -0.5*sample_length-0.001) RELATIVE arm_selene1
|
||||
|
||||
COMPONENT vs_divergence_v = DivPos_monitor(nb=21, ndiv=41, filename="vs_verdiv",
|
||||
xmin=-0.02, xmax=0.02, ymin=-0.02, ymax=0.02, maxdiv=2.0,
|
||||
restore_neutron=1, vertical=1)
|
||||
AT (0, 0, -0.5*sample_length-0.001) RELATIVE arm_selene1
|
||||
|
||||
/* The actual virtual source mask, two L-shaped absorbers (first top-right) */
|
||||
COMPONENT virtual_source_TR = Slit(
|
||||
@@ -300,8 +282,8 @@ COMPONENT ac_slit = Slit(
|
||||
* Sample area *
|
||||
***************/
|
||||
COMPONENT tof_sample = Monitor_nD(
|
||||
filename = "tof_sample", user1=p_int,
|
||||
options = "x limits=[-0.025 0.025] bins=1000 y limits=[-0.025 0.025] bins=1000 xdiv limits=[-0.75 0.75] bins=150 ydiv limits=[-2.0 2.0] bins=400 time limits=[0 0.6] bins=6000 lambda limits=[0 35] bins=3500 sy limits=[-1 1] bins=2000 user1 list all",
|
||||
filename = "tof_sample",
|
||||
options = "x limits=[-0.025 0.025] bins=1000 y limits=[-0.025 0.025] bins=1000 xdiv limits=[-0.75 0.75] bins=150 ydiv limits=[-2.0 2.0] bins=400 time limits=[0 0.6] bins=6000 lambda limits=[0 35] bins=3500 list all",
|
||||
xwidth=0.05, yheight = 0.05)
|
||||
WHEN sample==4
|
||||
AT (0, 0, 0) RELATIVE arm_sample_beam
|
||||
@@ -348,68 +330,63 @@ COMPONENT si_sample = Mirror(
|
||||
AT (0, 0, 0) RELATIVE arm_sample
|
||||
ROTATED (0, 90, 0) RELATIVE arm_sample
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
/* Rotate spin 90 degrees*/
|
||||
COMPONENT arm_hack = Arm()
|
||||
AT (0, 0, 0) RELATIVE arm_sample
|
||||
EXTEND %{
|
||||
double tmp = sy;
|
||||
sy = sx;
|
||||
sx = -tmp;
|
||||
%}
|
||||
|
||||
|
||||
|
||||
|
||||
/* flipped arm detector to position analyzer correctly for beam path 1*/
|
||||
COMPONENT arm_detector2 = Arm()
|
||||
AT (0, 0, 0) RELATIVE arm_detector
|
||||
ROTATED (0,0,180) RELATIVE arm_detector
|
||||
|
||||
/* Asymmetric lipid bilayer */
|
||||
COMPONENT lipid_h20_sample = Mirror(
|
||||
xwidth = sample_length, yheight = sample_height,
|
||||
center = 1, transmit = 0,
|
||||
reflect = "lipid_h2o.ref"
|
||||
)
|
||||
WHEN sample==4
|
||||
AT (0, 0, 0) RELATIVE arm_sample
|
||||
ROTATED (0, 90, 0) RELATIVE arm_sample
|
||||
COMPONENT lipid_d20_sample = Mirror(
|
||||
xwidth = sample_length, yheight = sample_height,
|
||||
center = 1, transmit = 0,
|
||||
reflect = "lipid_d2o.ref"
|
||||
)
|
||||
WHEN sample==5
|
||||
AT (0, 0, 0) RELATIVE arm_sample
|
||||
ROTATED (0, 90, 0) RELATIVE arm_sample
|
||||
COMPONENT lipid_smw_sample = Mirror(
|
||||
xwidth = sample_length, yheight = sample_height,
|
||||
center = 1, transmit = 0,
|
||||
reflect = "lipid_smw.ref"
|
||||
)
|
||||
WHEN sample==6
|
||||
AT (0, 0, 0) RELATIVE arm_sample
|
||||
ROTATED (0, 90, 0) RELATIVE arm_sample
|
||||
|
||||
|
||||
COMPONENT arm_analyzer = Arm()
|
||||
AT (0, 0, 0) RELATIVE arm_detector2
|
||||
ROTATED (-selene_theta+(Theta1_analyzer1-Theta2_analyzer1)/2.0, 0, 0) RELATIVE arm_detector2
|
||||
|
||||
COMPONENT virtual_analyzer_flipper = Arm() // Gone -> Pol_SF_ideal(ny=1, xwidth=1, yheight=1, zdepth=0.001)
|
||||
WHEN (enable_analyzer>2)
|
||||
AT (0, 0, analyzer_flipper_start) RELATIVE arm_analyzer
|
||||
ROTATED (0,0,0.0) RELATIVE arm_analyzer
|
||||
EXTEND %{
|
||||
if(INSTRUMENT_GETPAR(enable_analyzer)>2) sx *=-1;
|
||||
%}
|
||||
AT (0, 0, 0) RELATIVE arm_detector
|
||||
ROTATED (-selene_theta+(Theta1_analyzer1-Theta2_analyzer1)/2.0, 0, 0) RELATIVE arm_detector
|
||||
|
||||
COMPONENT arm_analyzer2 = Arm()
|
||||
AT (0, 0, 0) RELATIVE arm_detector2
|
||||
ROTATED (-selene_theta+(Theta1_analyzer2-Theta2_analyzer2)/2.0, 0, 0) RELATIVE arm_detector2
|
||||
AT (0, 0, 0) RELATIVE arm_detector
|
||||
ROTATED (-selene_theta+(Theta1_analyzer2-Theta2_analyzer2)/2.0, 0, 0) RELATIVE arm_detector
|
||||
|
||||
/* polarization analyser */
|
||||
COMPONENT analyzer1 = Polariser(enable_ref=1, d_substrate = 5e-4, reflect_d=0, reflect_u=0, lin=analyzer1_start, length=analyzer1_length,
|
||||
COMPONENT analyzer1 = Polariser(nIncRefr=1, d_substrate = 5e-4, reflect_d=0, reflect_u=0, lin=analyzer1_start, length=analyzer1_length,
|
||||
delta_theta=(Theta1_analyzer1+Theta2_analyzer1+0.05)*PI/180.0, h2=0.14, h1=0.05, abs_ref=0,
|
||||
m_u=5.5, m_d=0.45, both_coated=0, alpha=2.3, W = 0.0014)
|
||||
WHEN enable_analyzer
|
||||
AT (0, 0.0, analyzer1_start) RELATIVE arm_analyzer
|
||||
ROTATED (0,0,0.0) RELATIVE arm_analyzer
|
||||
|
||||
COMPONENT analyzer2 = Polariser(enable_ref=1, d_substrate = 5e-4, reflect_d=0, reflect_u=0, lin=analyzer2_start, length=analyzer2_length,
|
||||
COMPONENT analyzer2 = Polariser(nIncRefr=1, d_substrate = 5e-4, reflect_d=0, reflect_u=0, lin=analyzer2_start, length=analyzer2_length,
|
||||
delta_theta=(Theta1_analyzer2+Theta2_analyzer2+0.05)*PI/180.0, h2=0.2, h1=0.05, abs_out=0,
|
||||
m_u=5.0, m_d=0.65, both_coated=1, alpha=2.3, W = 0.0014)
|
||||
WHEN enable_analyzer==2
|
||||
AT (0, 0.0, analyzer2_start) RELATIVE arm_analyzer2
|
||||
ROTATED (0,0,0.0) RELATIVE arm_analyzer2
|
||||
AT (0, 0.0, analyzer2_start) RELATIVE arm_analyzer
|
||||
ROTATED (0,0,0.0) RELATIVE arm_analyzer
|
||||
|
||||
/* detector */
|
||||
COMPONENT tof_detector = Monitor_nD(
|
||||
filename = "tof_detector",
|
||||
options = "x limits=[-0.25 0.25] bins=1000 y limits=[-0.5 0.5] bins=1000 time limits=[0 0.6] bins=6000 lambda limits=[0 35] bins=3500 sx sy sz user1, list all",
|
||||
options = "x limits=[-0.25 0.25] bins=1000 y limits=[-0.5 0.5] bins=1000 time limits=[0 0.6] bins=6000 lambda limits=[0 35] bins=3500 sx limits=[-1 1] bins=1000 sy limits=[-1 1] bins=1000, list all",
|
||||
xwidth = 0.5, yheight = 1.0)
|
||||
WHEN sample!=4
|
||||
AT (0, 0, detector_arm+0.00001) RELATIVE arm_detector
|
||||
ROTATED (0, 0, 0) RELATIVE arm_detector
|
||||
|
||||
/***********************************************************************/
|
||||
|
||||
|
||||
@@ -144,7 +144,7 @@ COMPONENT NBOA_window=Al_window(thickness=NBOA_Al_entrance_length)
|
||||
AT (0,0,NBOA_Al_entrance_start) RELATIVE ISCS
|
||||
|
||||
|
||||
COMPONENT NBOA_side = Absorber(xmin=NBOA_side_x, xmax=NBOA_side_x+0.012,
|
||||
COMPONENT NBOA_side = AbsorberFixed(xmin=NBOA_side_x, xmax=NBOA_side_x+0.012,
|
||||
ymin=-0.2, ymax=0.2,
|
||||
zmin=0.0, zmax=E02_01_01_Cu_length+E02_01_01_length+E02_01_02_length+E02_01_03_length)
|
||||
AT (0, 0, E02_01_01_Cu_start) RELATIVE ISCS
|
||||
@@ -181,7 +181,7 @@ COMPONENT NBOA_Cu_collimator=Guide_gravity(
|
||||
* feeder neutron guide *
|
||||
************************/
|
||||
|
||||
COMPONENT E02_01_011_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
COMPONENT E02_01_011_wedge = AbsorberFixed(xmin=-0.1, xmax=0.0,
|
||||
ymin=-E02_01_011_wh/2.0, ymax=E02_01_011_wh/2.0,
|
||||
zmin=0.0, zmax=E02_01_01_length-0.5)
|
||||
AT (0, 0, E02_01_01_start) RELATIVE ISCS
|
||||
@@ -190,7 +190,7 @@ COMPONENT E02_01_011_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
ALLOW_BACKPROP;
|
||||
PROP_Z0;
|
||||
%}
|
||||
COMPONENT E02_01_012_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
COMPONENT E02_01_012_wedge = AbsorberFixed(xmin=-0.1, xmax=0.0,
|
||||
ymin=-E02_01_012_wh/2.0, ymax=E02_01_012_wh/2.0,
|
||||
zmin=0.5-E02_01_01_Cu_length, zmax=E02_01_01_length)
|
||||
AT (0, 0, E02_01_01_start) RELATIVE ISCS
|
||||
@@ -210,7 +210,7 @@ COMPONENT E02_01_01 = Elliptic_guide_gravity(
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, E02_01_01_start) RELATIVE ISCS
|
||||
|
||||
COMPONENT E02_01_021_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
COMPONENT E02_01_021_wedge = AbsorberFixed(xmin=-0.1, xmax=0.0,
|
||||
ymin=-E02_01_021_wh/2.0, ymax=E02_01_021_wh/2.0,
|
||||
zmin=0.0, zmax=E02_01_02_length/2.0)
|
||||
AT (0, 0, E02_01_02_start) RELATIVE ISCS
|
||||
@@ -219,7 +219,7 @@ COMPONENT E02_01_021_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
ALLOW_BACKPROP;
|
||||
PROP_Z0;
|
||||
%}
|
||||
COMPONENT E02_01_022_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
COMPONENT E02_01_022_wedge = AbsorberFixed(xmin=-0.1, xmax=0.0,
|
||||
ymin=-E02_01_022_wh/2.0, ymax=E02_01_022_wh/2.0,
|
||||
zmin=E02_01_02_length/2.0, zmax=E02_01_02_length)
|
||||
AT (0, 0, E02_01_02_start) RELATIVE ISCS
|
||||
@@ -238,7 +238,7 @@ COMPONENT E02_01_02 = Elliptic_guide_gravity(
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, E02_01_02_start) RELATIVE ISCS
|
||||
|
||||
COMPONENT E02_01_031_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
COMPONENT E02_01_031_wedge = AbsorberFixed(xmin=-0.1, xmax=0.0,
|
||||
ymin=-E02_01_031_wh/2.0, ymax=E02_01_031_wh/2.0,
|
||||
zmin=0.0, zmax=0.5)
|
||||
AT (0, 0, E02_01_03_start) RELATIVE ISCS
|
||||
@@ -247,7 +247,7 @@ COMPONENT E02_01_031_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
ALLOW_BACKPROP;
|
||||
PROP_Z0;
|
||||
%}
|
||||
COMPONENT E02_01_032_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
COMPONENT E02_01_032_wedge = AbsorberFixed(xmin=-0.1, xmax=0.0,
|
||||
ymin=-E02_01_032_wh/2.0, ymax=E02_01_032_wh/2.0,
|
||||
zmin=0.5, zmax=1.0)
|
||||
AT (0, 0, E02_01_03_start) RELATIVE ISCS
|
||||
@@ -256,7 +256,7 @@ COMPONENT E02_01_032_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
ALLOW_BACKPROP;
|
||||
PROP_Z0;
|
||||
%}
|
||||
COMPONENT E02_01_033_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
COMPONENT E02_01_033_wedge = AbsorberFixed(xmin=-0.1, xmax=0.0,
|
||||
ymin=-E02_01_033_wh/2.0, ymax=E02_01_033_wh/2.0,
|
||||
zmin=1.0, zmax=E02_01_03_length)
|
||||
AT (0, 0, E02_01_03_start) RELATIVE ISCS
|
||||
@@ -299,7 +299,7 @@ COMPONENT NFGA_Al_window_in=Al_window(
|
||||
WHEN enable_windows
|
||||
AT (0,0,NFGA_Al_start) RELATIVE ISCS
|
||||
|
||||
COMPONENT NFGA_side1 = Absorber(xmin=NFGA_side_x1, xmax=NFGA_side_x1+0.008,
|
||||
COMPONENT NFGA_side1 = AbsorberFixed(xmin=NFGA_side_x1, xmax=NFGA_side_x1+0.008,
|
||||
ymin=-0.2, ymax=0.2,
|
||||
zmin=0.002, zmax=0.502)
|
||||
AT (0, 0, E02_02_01_start-0.002) RELATIVE ISCS
|
||||
@@ -309,7 +309,7 @@ COMPONENT NFGA_side1 = Absorber(xmin=NFGA_side_x1, xmax=NFGA_side_x1+0.008,
|
||||
ALLOW_BACKPROP;
|
||||
PROP_Z0;
|
||||
%}
|
||||
COMPONENT NFGA_side2 = Absorber(xmin=NFGA_side_x2, xmax=NFGA_side_x2+0.008,
|
||||
COMPONENT NFGA_side2 = AbsorberFixed(xmin=NFGA_side_x2, xmax=NFGA_side_x2+0.008,
|
||||
ymin=-0.2, ymax=0.2,
|
||||
zmin=0.502, zmax=E02_02_01_length+E02_02_02_length+E02_02_03_length+0.002)
|
||||
AT (0, 0, E02_02_01_start-0.002) RELATIVE ISCS
|
||||
@@ -321,7 +321,7 @@ COMPONENT NFGA_side2 = Absorber(xmin=NFGA_side_x2, xmax=NFGA_side_x2+0.008,
|
||||
%}
|
||||
|
||||
// in-bunker feeder segments
|
||||
COMPONENT E02_02_011_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
COMPONENT E02_02_011_wedge = AbsorberFixed(xmin=-0.1, xmax=0.0,
|
||||
ymin=-E02_02_011_wh/2.0, ymax=E02_02_011_wh/2.0,
|
||||
zmin=0.0, zmax=E02_02_01_length/2.0)
|
||||
AT (0, 0, E02_02_01_start) RELATIVE ISCS
|
||||
@@ -330,7 +330,7 @@ COMPONENT E02_02_011_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
ALLOW_BACKPROP;
|
||||
PROP_Z0;
|
||||
%}
|
||||
COMPONENT E02_02_012_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
COMPONENT E02_02_012_wedge = AbsorberFixed(xmin=-0.1, xmax=0.0,
|
||||
ymin=-E02_02_012_wh/2.0, ymax=E02_02_012_wh/2.0,
|
||||
zmin=E02_02_01_length/2.0, zmax=E02_02_01_length)
|
||||
AT (0, 0, E02_02_01_start) RELATIVE ISCS
|
||||
@@ -349,7 +349,7 @@ COMPONENT E02_02_01 = Elliptic_guide_gravity(
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, E02_02_01_start) RELATIVE ISCS
|
||||
|
||||
COMPONENT E02_02_021_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
COMPONENT E02_02_021_wedge = AbsorberFixed(xmin=-0.1, xmax=0.0,
|
||||
ymin=-E02_02_021_wh/2.0, ymax=E02_02_021_wh/2.0,
|
||||
zmin=0.0, zmax=E02_02_02_length/2.0)
|
||||
AT (0, 0, E02_02_02_start) RELATIVE ISCS
|
||||
@@ -358,7 +358,7 @@ COMPONENT E02_02_021_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
ALLOW_BACKPROP;
|
||||
PROP_Z0;
|
||||
%}
|
||||
COMPONENT E02_02_022_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
COMPONENT E02_02_022_wedge = AbsorberFixed(xmin=-0.1, xmax=0.0,
|
||||
ymin=-E02_02_022_wh/2.0, ymax=E02_02_022_wh/2.0,
|
||||
zmin=E02_02_02_length/2.0, zmax=E02_02_02_length)
|
||||
AT (0, 0, E02_02_02_start) RELATIVE ISCS
|
||||
@@ -377,7 +377,7 @@ COMPONENT E02_02_02 = Elliptic_guide_gravity(
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, E02_02_02_start) RELATIVE ISCS
|
||||
|
||||
COMPONENT E02_02_031_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
COMPONENT E02_02_031_wedge = AbsorberFixed(xmin=-0.1, xmax=0.0,
|
||||
ymin=-E02_02_031_wh/2.0, ymax=E02_02_031_wh/2.0,
|
||||
zmin=0.0, zmax=0.5)
|
||||
AT (0, 0, E02_02_03_start) RELATIVE ISCS
|
||||
@@ -386,7 +386,7 @@ COMPONENT E02_02_031_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
ALLOW_BACKPROP;
|
||||
PROP_Z0;
|
||||
%}
|
||||
COMPONENT E02_02_032_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
COMPONENT E02_02_032_wedge = AbsorberFixed(xmin=-0.1, xmax=0.0,
|
||||
ymin=-E02_02_032_wh/2.0, ymax=E02_02_032_wh/2.0,
|
||||
zmin=0.5, zmax=1.0)
|
||||
AT (0, 0, E02_02_03_start) RELATIVE ISCS
|
||||
@@ -395,7 +395,7 @@ COMPONENT E02_02_032_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
ALLOW_BACKPROP;
|
||||
PROP_Z0;
|
||||
%}
|
||||
COMPONENT E02_02_033_wedge = Absorber(xmin=-0.1, xmax=0.0,
|
||||
COMPONENT E02_02_033_wedge = AbsorberFixed(xmin=-0.1, xmax=0.0,
|
||||
ymin=-E02_02_033_wh/2.0, ymax=E02_02_033_wh/2.0,
|
||||
zmin=1.0, zmax=E02_02_03_length)
|
||||
AT (0, 0, E02_02_03_start) RELATIVE ISCS
|
||||
|
||||
+21
-167
@@ -28,80 +28,20 @@ DEFINE INSTRUMENT Estia_selene(int enable_gravity=0)
|
||||
|
||||
DECLARE
|
||||
%{
|
||||
|
||||
// Polarizer parameters used for the detailed design
|
||||
double pol1_xstart = 1.007; // [m] up-stream of focus point
|
||||
double pol1_ystart = 0.0064; // [m] (relative to c-axis)
|
||||
double pol1_xend = 0.365; // [m] up-stream of focus point
|
||||
double pol1_yend = 0.013; // [m]
|
||||
double pol1_gamma = 1.65; // [deg] optimal incident angle for beams passing through focus point
|
||||
double pol1_m = 5.08;
|
||||
double pol1_lin_xstart = 1.007; //[m]
|
||||
double pol1_lin_ystart = 0.0064; //[m]
|
||||
double pol1_lin_xend = 0.812; //[m]
|
||||
double pol1_lin_yend = -0.0014; //[m]
|
||||
double pol1_lin_m = 5.5;
|
||||
|
||||
|
||||
double pol2_xstart = 1.286; // [m] up-stream of focus point
|
||||
double pol2_ystart = 0.009; // [m] (relative to c-axis)
|
||||
double pol2_xend = 0.500; // [m] up-stream of focus point
|
||||
double pol2_yend = 0.0175; // [m]
|
||||
double pol2_gamma = 1.70; // [deg] optimal incident angle for beams passing through focus point
|
||||
double pol2_m = 5.08;
|
||||
double pol2_lin_xstart = 1.286; //[m]
|
||||
double pol2_lin_ystart = 0.009; //[m]
|
||||
double pol2_lin_xend = 1.093; //[m]
|
||||
double pol2_lin_yend = 0.0014; //[m]
|
||||
double pol2_lin_m = 5.08;
|
||||
|
||||
double pol_hclose = 0.12; // [m] height of entrance
|
||||
double pol_hfar = 0.12; // [m] height of exit
|
||||
|
||||
// parameters for polarizer component calculated from the above
|
||||
double pol1_length, pol2_length, pol1_lin_length, pol2_lin_length;
|
||||
double Theta_pol1, Theta_pol2, Theta_rot1, Theta_rot2;
|
||||
double Theta1_pol1, Theta1_pol2, Theta2_pol1, Theta2_pol2;
|
||||
double pol1_lin_rot, pol2_lin_rot;
|
||||
|
||||
double polarizer_max_width = 0.0015; // Maximum sample height to be covered by te polarization analyzers
|
||||
double polarizer_start = 0.3; // Distance to sample to start the first analyzer
|
||||
double polarizer_length = 0.6; // length of first analyzer
|
||||
double Theta1_polarizer, Theta2_polarizer, dist_pol_vfocus; // quantities calculated out of values above and 1.5 degree covered divergence
|
||||
|
||||
%}
|
||||
|
||||
INITIALIZE
|
||||
%{
|
||||
dist_pol_vfocus=polarizer_max_width/2.0/tan(1.5*PI/360.0); // the virtual focus point infront of the actual sample focus where the beams furthest out meet
|
||||
|
||||
pol1_length = pol1_xstart-pol1_xend;
|
||||
pol2_length = pol2_xstart-pol2_xend;
|
||||
pol1_lin_length = pol1_lin_xstart-pol1_lin_xend;
|
||||
pol2_lin_length = pol2_lin_xstart-pol2_lin_xend;
|
||||
|
||||
Theta1_pol1=atan2(pol1_ystart, pol1_xstart); // [rad] angle between c-axis and entrance point
|
||||
Theta2_pol1=atan2(pol1_yend, pol1_xend); // [rad] angle between c-axis and exit point
|
||||
Theta1_pol2=atan2(pol2_ystart, pol2_xstart); // [rad] angle between c-axis and entrance point
|
||||
Theta2_pol2=atan2(pol2_yend, pol2_xend); // [rad] angle between c-axis and exit point
|
||||
|
||||
Theta_pol1=(Theta2_pol1-Theta1_pol1)*180.0/PI; // [deg] full covered divergence angle
|
||||
Theta_pol2=(Theta2_pol2-Theta1_pol2)*180.0/PI; // [deg]
|
||||
Theta_rot1=(Theta1_pol1+Theta2_pol1)/2.0*180.0/PI; // [deg] rotation of center of polarizer
|
||||
Theta_rot2=(Theta1_pol2+Theta2_pol2)/2.0*180.0/PI; // [deg]
|
||||
|
||||
pol1_lin_rot=atan2(pol1_lin_ystart-pol1_lin_yend,
|
||||
pol1_lin_xstart-pol1_lin_xend)*180.0/PI; // [deg] rotation angle of polarizer
|
||||
pol2_lin_rot=atan2(pol2_lin_ystart-pol2_lin_yend,
|
||||
pol2_lin_xstart-pol2_lin_xend)*180.0/PI; // [deg] rotation angle of polarizer
|
||||
|
||||
|
||||
|
||||
printf(" Polarizer 1 angle = %.2f deg\n", Theta_rot1);
|
||||
printf(" Polarizer 1 divergence = %.2f deg\n", Theta_pol1);
|
||||
printf(" Polarizer 1 distance = %.1f mm\n", pol1_xstart*1000.0);
|
||||
printf(" Polarizer 1 length = %.1f mm\n", pol1_length*1000.0);
|
||||
printf(" PolStraight 1 angle = %.2f deg\n", pol1_lin_rot);
|
||||
|
||||
printf(" Polarizer 2 angle = %.2f deg\n", Theta_rot1);
|
||||
printf(" Polarizer 2 divergence = %.2f deg\n", Theta_pol2);
|
||||
printf(" Polarizer 2 distance = %.1f mm\n", pol2_xstart*1000.0);
|
||||
printf(" Polarizer 2 length = %.1f mm\n", pol2_length*1000.0);
|
||||
printf(" PolStraight 2 angle = %.2f deg\n", pol2_lin_rot);
|
||||
Theta1_polarizer=atan((dist_pol_vfocus+polarizer_start)/dist_pol_vfocus*polarizer_max_width/2.0/polarizer_start)*180.0/PI;
|
||||
Theta2_polarizer=atan((dist_pol_vfocus+polarizer_start+polarizer_length)/dist_pol_vfocus*polarizer_max_width/2.0/(polarizer_start+polarizer_length))*180.0/PI;
|
||||
|
||||
%}
|
||||
TRACE
|
||||
@@ -127,112 +67,26 @@ REMOVABLE COMPONENT source = Moderator(radius = 0.001, focus_xw = selene_entry+0
|
||||
/**************************************
|
||||
* Middle focus between Selene guides *
|
||||
**************************************/
|
||||
COMPONENT arm_polref = Arm()
|
||||
COMPONENT arm_polarizer = Arm()
|
||||
AT (0, 0, 2*selene_c) RELATIVE arm_selene1
|
||||
ROTATED (selene_theta, -selene_theta, 0) RELATIVE ISCS
|
||||
|
||||
COMPONENT arm_polarizer = Arm()
|
||||
AT (0, 0, 2*selene_c) RELATIVE arm_selene1
|
||||
ROTATED (0, 0, 0) RELATIVE arm_selene1
|
||||
|
||||
COMPONENT arm_pol1 = Arm()
|
||||
AT (0, 0, 0) RELATIVE arm_polarizer
|
||||
ROTATED (0, -Theta_rot1, 0) RELATIVE arm_polarizer
|
||||
|
||||
COMPONENT arm_pol2 = Arm()
|
||||
AT (0, 0, 0) RELATIVE arm_polarizer
|
||||
ROTATED (0, -Theta_rot2, 0) RELATIVE arm_polarizer
|
||||
COMPONENT polarizer1 = Polariser(nIncRefr=1, d_substrate = 5e-4, reflect_d=0, reflect_u=0, lin=-(polarizer_start+polarizer_length), length=polarizer_length,
|
||||
delta_theta=(Theta1_polarizer+Theta2_polarizer)*PI/180.0, h2=0.1, h1=0.05, abs_ref=1, m_u=4.0, m_d=0.65, both_coated=1, alpha=2.3, W = 0.0014)
|
||||
WHEN enable_polarizer
|
||||
AT (0, 0, -(polarizer_start+polarizer_length)) RELATIVE arm_polarizer
|
||||
ROTATED (0,0,90.0) RELATIVE arm_polarizer
|
||||
|
||||
|
||||
COMPONENT polarizer2_lin = Pol_mirror(zwidth = pol2_lin_length, yheight = pol_hfar, p_reflect=0,
|
||||
rUpData="SNAG_m5p5_Tup.ref", rDownData="SNAG_m5p5_Tdown.ref")
|
||||
WHEN (enable_polarizer>0)
|
||||
AT ((pol2_lin_ystart+pol2_lin_yend)/2.0, 0, -(pol2_lin_xstart+pol2_lin_xend)/2.0) RELATIVE arm_polarizer
|
||||
ROTATED (0,-pol2_lin_rot,0.0) RELATIVE arm_polarizer
|
||||
GROUP polarizer2_set
|
||||
EXTEND
|
||||
%{
|
||||
if (SCATTERED==2) {
|
||||
p_int +=2;
|
||||
}
|
||||
%}
|
||||
|
||||
COMPONENT polarizer2 = Polariser(lin=-pol2_xstart, length=pol2_length,
|
||||
enable_ref=1, abs_ref=1, abs_out=1, both_coated=1,
|
||||
d_substrate = 5e-4, T_loss=4.0e3,
|
||||
m_u=pol2_m, m_d=0.6, m_residual=0.55,
|
||||
alpha=2.3, W = 0.0014, reflect_d=0, reflect_u=0,
|
||||
delta_theta=Theta_pol2*PI/180.0,
|
||||
h2=pol_hfar, h1=pol_hclose)
|
||||
WHEN (enable_polarizer>0)
|
||||
AT (0, 0, -pol2_xstart) RELATIVE arm_pol2
|
||||
ROTATED (0,0,90.0) RELATIVE arm_pol2
|
||||
GROUP polarizer2_set
|
||||
EXTEND
|
||||
%{
|
||||
if (SCATTERED) {
|
||||
ALLOW_BACKPROP;
|
||||
PROP_Z0;
|
||||
p_int +=1;
|
||||
} else ABSORB;
|
||||
%}
|
||||
|
||||
COMPONENT replacement_pol2_set = Slit(xwidth=0.2, yheight=0.2)
|
||||
WHEN (enable_polarizer==0)
|
||||
AT (0, 0, -0.2) RELATIVE arm_polarizer
|
||||
GROUP polarizer2_set
|
||||
|
||||
COMPONENT polarizer1_lin = Pol_mirror(zwidth = pol1_lin_length, yheight = pol_hfar, p_reflect=0,
|
||||
rUpData="SNAG_m5p5_Tup.ref", rDownData="SNAG_m5p5_Tdown.ref")
|
||||
WHEN (enable_polarizer>0)
|
||||
AT ((pol1_lin_ystart+pol1_lin_yend)/2.0, 0, -(pol1_lin_xstart+pol1_lin_xend)/2.0) RELATIVE arm_polarizer
|
||||
ROTATED (0,-pol1_lin_rot,0.0) RELATIVE arm_polarizer
|
||||
GROUP polarizer1_set
|
||||
EXTEND
|
||||
%{
|
||||
if (SCATTERED==2) {
|
||||
p_int +=8;
|
||||
}
|
||||
%}
|
||||
|
||||
COMPONENT polarizer1 = Polariser(lin=-pol1_xstart, length=pol1_length,
|
||||
enable_ref=1, abs_ref=1, abs_out=1, both_coated=1,
|
||||
d_substrate = 5e-4, T_loss=4.0e3,
|
||||
m_u=pol1_m, m_d=0.6, m_residual=0.55,
|
||||
alpha=2.3, W = 0.0014, reflect_d=0, reflect_u=0,
|
||||
delta_theta=Theta_pol1*PI/180.0,
|
||||
h2=pol_hfar, h1=pol_hclose)
|
||||
WHEN (enable_polarizer>0)
|
||||
AT (0, 0, -pol1_xstart) RELATIVE arm_pol1
|
||||
ROTATED (0,0,90.0) RELATIVE arm_pol1
|
||||
GROUP polarizer1_set
|
||||
EXTEND
|
||||
%{
|
||||
if (SCATTERED) {
|
||||
p_int +=4;
|
||||
}
|
||||
%}
|
||||
|
||||
COMPONENT replacement_pol1_set = Slit(xwidth=0.2, yheight=0.2)
|
||||
WHEN (enable_polarizer==0)
|
||||
AT (0, 0, -0.1) RELATIVE arm_polarizer
|
||||
GROUP polarizer1_set
|
||||
/* The proximal polariser comes next */
|
||||
COMPONENT polarizer2 = Polariser(nIncRefr=1, d_substrate = 5e-4, reflect_d=0, reflect_u=0, lin=polarizer_start, length=polarizer_length,
|
||||
delta_theta=(Theta1_polarizer+Theta2_polarizer)*PI/180.0, h2=0.1, h1=0.05, abs_ref=1, m_u=4.0, m_d=0.65, both_coated=1, alpha=2.3, W = 0.0014)
|
||||
WHEN (enable_polarizer>1)
|
||||
AT (0, 0, polarizer_start) RELATIVE arm_polarizer
|
||||
ROTATED (0.0,0,90.0) RELATIVE arm_polarizer
|
||||
|
||||
|
||||
|
||||
|
||||
COMPONENT mf_divergence_h = DivPos_monitor(nb=21, ndiv=41, filename="mf_hordiv",
|
||||
xmin=-0.02, xmax=0.02, ymin=-0.02, ymax=0.02, maxdiv=2.0)
|
||||
AT (0, 0, 0) RELATIVE arm_polarizer
|
||||
|
||||
COMPONENT mf_divergence_v = DivPos_monitor(nb=21, ndiv=41, filename="mf_verdiv",
|
||||
xmin=-0.02, xmax=0.02, ymin=-0.02, ymax=0.02, maxdiv=2.0, vertical=1)
|
||||
AT (0, 0, 0) RELATIVE arm_polarizer
|
||||
|
||||
|
||||
|
||||
COMPONENT virtual_flipper = Pol_SF_ideal(ny=1, xwidth=0.1, yheight=0.1, zdepth=0.001)
|
||||
WHEN (enable_polarizer>2)
|
||||
AT (0, 0, 0) RELATIVE arm_polarizer
|
||||
|
||||
|
||||
FINALLY
|
||||
|
||||
@@ -102,7 +102,7 @@ REMOVABLE COMPONENT source = Moderator(radius = 0.001, focus_xw = selene_entry+0
|
||||
/* Absorber to cut direct view beam (Copper in CAD model) */
|
||||
COMPONENT slit_before_selene_guide_1 = Slit(
|
||||
xmin = 0, xmax = selene_entry+0.001,
|
||||
ymin = -selene_entry-0.005, ymax = 0)
|
||||
ymin = 0, ymax = selene_entry+0.005)
|
||||
AT (0, 0, selene_distance-0.01) RELATIVE arm_selene1
|
||||
|
||||
COMPONENT block_before_selene_guide_1 = Absorber(
|
||||
@@ -118,7 +118,7 @@ COMPONENT E02_03_01 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance, loutyh= selene_distance+14*selene_segment,
|
||||
linxw = selene_distance, loutxw= selene_distance+14*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=0, mleft=scoating_01, mbottom=scoating_01, mtop=0,
|
||||
mright=0, mleft=scoating_01, mtop=scoating_01, mbottom=0,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance) RELATIVE arm_selene1
|
||||
|
||||
@@ -127,7 +127,7 @@ COMPONENT E02_03_02 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+1*selene_segment, loutyh= selene_distance+13*selene_segment,
|
||||
linxw = selene_distance+1*selene_segment, loutxw= selene_distance+13*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=0, mleft=scoating_02, mbottom=scoating_02, mtop=0,
|
||||
mright=0, mleft=scoating_02, mtop=scoating_02, mbottom=0,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+1*selene_segment) RELATIVE arm_selene1
|
||||
|
||||
@@ -136,7 +136,7 @@ COMPONENT E02_03_03 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+2*selene_segment, loutyh= selene_distance+12*selene_segment,
|
||||
linxw = selene_distance+2*selene_segment, loutxw= selene_distance+12*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=0, mleft=scoating_03, mbottom=scoating_03, mtop=0,
|
||||
mright=0, mleft=scoating_03, mtop=scoating_03, mbottom=0,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+2*selene_segment) RELATIVE arm_selene1
|
||||
|
||||
@@ -145,7 +145,7 @@ COMPONENT E02_03_04 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+3*selene_segment, loutyh= selene_distance+11*selene_segment,
|
||||
linxw = selene_distance+3*selene_segment, loutxw= selene_distance+11*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=0, mleft=scoating_04, mbottom=scoating_04, mtop=0,
|
||||
mright=0, mleft=scoating_04, mtop=scoating_04, mbottom=0,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+3*selene_segment) RELATIVE arm_selene1
|
||||
|
||||
@@ -154,7 +154,7 @@ COMPONENT E02_03_05 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+4*selene_segment, loutyh= selene_distance+10*selene_segment,
|
||||
linxw = selene_distance+4*selene_segment, loutxw= selene_distance+10*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=0, mleft=scoating_05, mbottom=scoating_05, mtop=0,
|
||||
mright=0, mleft=scoating_05, mtop=scoating_05, mbottom=0,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+4*selene_segment) RELATIVE arm_selene1
|
||||
|
||||
@@ -163,7 +163,7 @@ COMPONENT E02_03_06 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+5*selene_segment, loutyh= selene_distance+9*selene_segment,
|
||||
linxw = selene_distance+5*selene_segment, loutxw= selene_distance+9*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=0, mleft=scoating_06, mbottom=scoating_06, mtop=0,
|
||||
mright=0, mleft=scoating_06, mtop=scoating_06, mbottom=0,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+5*selene_segment) RELATIVE arm_selene1
|
||||
|
||||
@@ -172,7 +172,7 @@ COMPONENT E02_03_07 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+6*selene_segment, loutyh= selene_distance+8*selene_segment,
|
||||
linxw = selene_distance+6*selene_segment, loutxw= selene_distance+8*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=0, mleft=scoating_07, mbottom=scoating_07, mtop=0,
|
||||
mright=0, mleft=scoating_07, mtop=scoating_07, mbottom=0,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+6*selene_segment) RELATIVE arm_selene1
|
||||
|
||||
@@ -181,7 +181,7 @@ COMPONENT E02_03_08 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+7*selene_segment, loutyh= selene_distance+7*selene_segment,
|
||||
linxw = selene_distance+7*selene_segment, loutxw= selene_distance+7*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=0, mleft=scoating_08, mbottom=scoating_08, mtop=0,
|
||||
mright=0, mleft=scoating_08, mtop=scoating_08, mbottom=0,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+7*selene_segment) RELATIVE arm_selene1
|
||||
|
||||
@@ -190,7 +190,7 @@ COMPONENT E02_03_09 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+8*selene_segment, loutyh= selene_distance+6*selene_segment,
|
||||
linxw = selene_distance+8*selene_segment, loutxw= selene_distance+6*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=0, mleft=scoating_09, mbottom=scoating_09, mtop=0,
|
||||
mright=0, mleft=scoating_09, mtop=scoating_09, mbottom=0,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+8*selene_segment) RELATIVE arm_selene1
|
||||
|
||||
@@ -199,7 +199,7 @@ COMPONENT E02_03_10 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+9*selene_segment, loutyh= selene_distance+5*selene_segment,
|
||||
linxw = selene_distance+9*selene_segment, loutxw= selene_distance+5*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=0, mleft=scoating_10, mbottom=scoating_10, mtop=0,
|
||||
mright=0, mleft=scoating_10, mtop=scoating_10, mbottom=0,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+9*selene_segment) RELATIVE arm_selene1
|
||||
|
||||
@@ -208,7 +208,7 @@ COMPONENT E02_03_11 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+10*selene_segment, loutyh= selene_distance+4*selene_segment,
|
||||
linxw = selene_distance+10*selene_segment, loutxw= selene_distance+4*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=0, mleft=scoating_11, mbottom=scoating_11, mtop=0,
|
||||
mright=0, mleft=scoating_11, mtop=scoating_11, mbottom=0,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+10*selene_segment) RELATIVE arm_selene1
|
||||
|
||||
@@ -217,7 +217,7 @@ COMPONENT E02_03_12 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+11*selene_segment, loutyh= selene_distance+3*selene_segment,
|
||||
linxw = selene_distance+11*selene_segment, loutxw= selene_distance+3*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=0, mleft=scoating_12, mbottom=scoating_12, mtop=0,
|
||||
mright=0, mleft=scoating_12, mtop=scoating_12, mbottom=0,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+11*selene_segment) RELATIVE arm_selene1
|
||||
|
||||
@@ -226,7 +226,7 @@ COMPONENT E02_03_13 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+12*selene_segment, loutyh= selene_distance+2*selene_segment,
|
||||
linxw = selene_distance+12*selene_segment, loutxw= selene_distance+2*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=0, mleft=scoating_13, mbottom=scoating_13, mtop=0,
|
||||
mright=0, mleft=scoating_13, mtop=scoating_13, mbottom=0,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+12*selene_segment) RELATIVE arm_selene1
|
||||
|
||||
@@ -235,7 +235,7 @@ COMPONENT E02_03_14 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+13*selene_segment, loutyh= selene_distance+1*selene_segment,
|
||||
linxw = selene_distance+13*selene_segment, loutxw= selene_distance+1*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=0, mleft=scoating_14, mbottom=scoating_14, mtop=0,
|
||||
mright=0, mleft=scoating_14, mtop=scoating_14, mbottom=0,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+13*selene_segment) RELATIVE arm_selene1
|
||||
|
||||
@@ -244,14 +244,14 @@ COMPONENT E02_03_15 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+14*selene_segment, loutyh= selene_distance,
|
||||
linxw = selene_distance+14*selene_segment, loutxw= selene_distance,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=0, mleft=scoating_15, mbottom=scoating_15, mtop=0,
|
||||
mright=0, mleft=scoating_15, mtop=scoating_15, mbottom=0,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+14*selene_segment) RELATIVE arm_selene1
|
||||
|
||||
/* Absorber to cut direct view beam (Copper in CAD model) */
|
||||
COMPONENT slit_after_selene_guide_1 = Slit(
|
||||
xmin = 0, xmax = selene_entry+0.005,
|
||||
ymin = -selene_entry-0.005, ymax = 0)
|
||||
ymin = 0, ymax = selene_entry+0.01)
|
||||
AT (0, 0, 2*selene_c-selene_distance+0.001) RELATIVE arm_selene1
|
||||
|
||||
COMPONENT block_after_selene_guide_1 = Absorber(
|
||||
|
||||
@@ -42,7 +42,7 @@ TRACE
|
||||
/* Absorber to cut direct view beam (Bor-Al in CAD model) */
|
||||
COMPONENT slit_before_selene_guide_2 = Slit(
|
||||
xmin = -selene_entry-0.005, xmax=0,
|
||||
ymin = 0, ymax = selene_entry+0.01)
|
||||
ymin = -selene_entry-0.01, ymax = 0)
|
||||
AT (0, 0, selene_distance-0.002) RELATIVE arm_selene2
|
||||
|
||||
COMPONENT block_before_selene_guide_2 = Absorber(
|
||||
@@ -59,7 +59,7 @@ COMPONENT E02_04_01 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance, loutyh= selene_distance+14*selene_segment,
|
||||
linxw = selene_distance, loutxw= selene_distance+14*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=scoating_15, mleft=0, mbottom=0, mtop=scoating_15,
|
||||
mright=scoating_15, mleft=0, mtop=0, mbottom=scoating_15,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance) RELATIVE arm_selene2
|
||||
|
||||
@@ -68,7 +68,7 @@ COMPONENT E02_04_02 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+1*selene_segment, loutyh= selene_distance+13*selene_segment,
|
||||
linxw = selene_distance+1*selene_segment, loutxw= selene_distance+13*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=scoating_14, mleft=0, mbottom=0, mtop=scoating_14,
|
||||
mright=scoating_14, mleft=0, mtop=0, mbottom=scoating_14,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+1*selene_segment) RELATIVE arm_selene2
|
||||
|
||||
@@ -77,7 +77,7 @@ COMPONENT E02_04_03 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+2*selene_segment, loutyh= selene_distance+12*selene_segment,
|
||||
linxw = selene_distance+2*selene_segment, loutxw= selene_distance+12*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=scoating_13, mleft=0, mbottom=0, mtop=scoating_13,
|
||||
mright=scoating_13, mleft=0, mtop=0, mbottom=scoating_13,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+2*selene_segment) RELATIVE arm_selene2
|
||||
|
||||
@@ -86,7 +86,7 @@ COMPONENT E02_04_04 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+3*selene_segment, loutyh= selene_distance+11*selene_segment,
|
||||
linxw = selene_distance+3*selene_segment, loutxw= selene_distance+11*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=scoating_12, mleft=0, mbottom=0, mtop=scoating_12,
|
||||
mright=scoating_12, mleft=0, mtop=0, mbottom=scoating_12,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+3*selene_segment) RELATIVE arm_selene2
|
||||
|
||||
@@ -95,7 +95,7 @@ COMPONENT E02_04_05 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+4*selene_segment, loutyh= selene_distance+10*selene_segment,
|
||||
linxw = selene_distance+4*selene_segment, loutxw= selene_distance+10*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=scoating_11, mleft=0, mbottom=0, mtop=scoating_11,
|
||||
mright=scoating_11, mleft=0, mtop=0, mbottom=scoating_11,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+4*selene_segment) RELATIVE arm_selene2
|
||||
|
||||
@@ -104,7 +104,7 @@ COMPONENT E02_04_06 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+5*selene_segment, loutyh= selene_distance+9*selene_segment,
|
||||
linxw = selene_distance+5*selene_segment, loutxw= selene_distance+9*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=scoating_10, mleft=0, mbottom=0, mtop=scoating_10,
|
||||
mright=scoating_10, mleft=0, mtop=0, mbottom=scoating_10,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+5*selene_segment) RELATIVE arm_selene2
|
||||
|
||||
@@ -113,7 +113,7 @@ COMPONENT E02_04_07 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+6*selene_segment, loutyh= selene_distance+8*selene_segment,
|
||||
linxw = selene_distance+6*selene_segment, loutxw= selene_distance+8*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=scoating_09, mleft=0, mbottom=0, mtop=scoating_09,
|
||||
mright=scoating_09, mleft=0, mtop=0, mbottom=scoating_09,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+6*selene_segment) RELATIVE arm_selene2
|
||||
|
||||
@@ -122,7 +122,7 @@ COMPONENT E02_04_08 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+7*selene_segment, loutyh= selene_distance+7*selene_segment,
|
||||
linxw = selene_distance+7*selene_segment, loutxw= selene_distance+7*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=scoating_08, mleft=0, mbottom=0, mtop=scoating_08,
|
||||
mright=scoating_08, mleft=0, mtop=0, mbottom=scoating_08,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+7*selene_segment) RELATIVE arm_selene2
|
||||
|
||||
@@ -131,7 +131,7 @@ COMPONENT E02_04_09 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+8*selene_segment, loutyh= selene_distance+6*selene_segment,
|
||||
linxw = selene_distance+8*selene_segment, loutxw= selene_distance+6*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=scoating_07, mleft=0, mbottom=0, mtop=scoating_07,
|
||||
mright=scoating_07, mleft=0, mtop=0, mbottom=scoating_07,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+8*selene_segment) RELATIVE arm_selene2
|
||||
|
||||
@@ -140,7 +140,7 @@ COMPONENT E02_04_10 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+9*selene_segment, loutyh= selene_distance+5*selene_segment,
|
||||
linxw = selene_distance+9*selene_segment, loutxw= selene_distance+5*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=scoating_06, mleft=0, mbottom=0, mtop=scoating_06,
|
||||
mright=scoating_06, mleft=0, mtop=0, mbottom=scoating_06,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+9*selene_segment) RELATIVE arm_selene2
|
||||
|
||||
@@ -149,7 +149,7 @@ COMPONENT E02_04_11 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+10*selene_segment, loutyh= selene_distance+4*selene_segment,
|
||||
linxw = selene_distance+10*selene_segment, loutxw= selene_distance+4*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=scoating_05, mleft=0, mbottom=0, mtop=scoating_05,
|
||||
mright=scoating_05, mleft=0, mtop=0, mbottom=scoating_05,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+10*selene_segment) RELATIVE arm_selene2
|
||||
|
||||
@@ -158,7 +158,7 @@ COMPONENT E02_04_12 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+11*selene_segment, loutyh= selene_distance+3*selene_segment,
|
||||
linxw = selene_distance+11*selene_segment, loutxw= selene_distance+3*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=scoating_04, mleft=0, mbottom=0, mtop=scoating_04,
|
||||
mright=scoating_04, mleft=0, mtop=0, mbottom=scoating_04,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+11*selene_segment) RELATIVE arm_selene2
|
||||
|
||||
@@ -167,7 +167,7 @@ COMPONENT E02_04_13 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+12*selene_segment, loutyh= selene_distance+2*selene_segment,
|
||||
linxw = selene_distance+12*selene_segment, loutxw= selene_distance+2*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=scoating_03, mleft=0, mbottom=0, mtop=scoating_03,
|
||||
mright=scoating_03, mleft=0, mtop=0, mbottom=scoating_03,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+12*selene_segment) RELATIVE arm_selene2
|
||||
|
||||
@@ -176,7 +176,7 @@ COMPONENT E02_04_14 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+13*selene_segment, loutyh= selene_distance+1*selene_segment,
|
||||
linxw = selene_distance+13*selene_segment, loutxw= selene_distance+1*selene_segment,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=scoating_02, mleft=0, mbottom=0, mtop=scoating_02,
|
||||
mright=scoating_02, mleft=0, mtop=0, mbottom=scoating_02,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+13*selene_segment) RELATIVE arm_selene2
|
||||
|
||||
@@ -185,14 +185,14 @@ COMPONENT E02_04_15 = Elliptic_guide_gravity(
|
||||
linyh = selene_distance+14*selene_segment, loutyh= selene_distance,
|
||||
linxw = selene_distance+14*selene_segment, loutxw= selene_distance,
|
||||
xwidth=selene_b*2, yheight=selene_b*2,
|
||||
mright=scoating_01, mleft=0, mbottom=0, mtop=scoating_01,
|
||||
mright=scoating_01, mleft=0, mtop=0, mbottom=scoating_01,
|
||||
enableGravity=enable_gravity)
|
||||
AT (0, 0, selene_distance+14*selene_segment) RELATIVE arm_selene2
|
||||
|
||||
/* Absorber to cut direct view beam (Bor-Al in CAD model) */
|
||||
COMPONENT slit_after_selene_guide_2 = Slit(
|
||||
xmin = -selene_entry-0.001, xmax=0, xmax = 0.0,
|
||||
ymin = 0, ymax = selene_entry+0.01)
|
||||
ymin = -selene_entry-0.005, ymax = 0)
|
||||
AT (0, 0, 2*selene_c-selene_distance+0.001) RELATIVE arm_selene2
|
||||
|
||||
COMPONENT block_after_selene_guide_2 = Absorber(
|
||||
|
||||
+515
-697
File diff suppressed because it is too large
Load Diff
@@ -1,58 +0,0 @@
|
||||
4.37173E-05 0.253838498
|
||||
0.001967279 0.548126562
|
||||
0.004240578 0.781912929
|
||||
0.006513878 0.831781609
|
||||
0.008656026 0.935568791
|
||||
0.010841891 0.996046699
|
||||
0.012984039 0.994254859
|
||||
0.015038752 0.780745054
|
||||
0.017443204 0.571166135
|
||||
0.019497917 0.421992388
|
||||
0.021771217 0.347878292
|
||||
0.023913365 0.314863976
|
||||
0.026317816 0.273648917
|
||||
0.028241378 0.248910367
|
||||
0.030514678 0.229613667
|
||||
0.032700543 0.210315136
|
||||
0.034711539 0.201079179
|
||||
0.03711599 0.179153501
|
||||
0.039301855 0.174572119
|
||||
0.041444003 0.161570874
|
||||
0.043629868 0.159330731
|
||||
0.045772016 0.158627627
|
||||
0.047957881 0.133760075
|
||||
0.050231181 0.128230852
|
||||
0.052373329 0.135460425
|
||||
0.054559194 0.125901964
|
||||
0.056701342 0.120064418
|
||||
0.058974642 0.122922579
|
||||
0.061160507 0.127084935
|
||||
0.063171503 0.12641248
|
||||
0.06548852 0.117637177
|
||||
0.067630668 0.105706827
|
||||
0.069903968 0.112415362
|
||||
0.072089833 0.119429932
|
||||
0.074231981 0.105676178
|
||||
0.076417846 0.107716872
|
||||
0.078428842 0.107129046
|
||||
0.080833293 0.10262909
|
||||
0.082888007 0.094164854
|
||||
0.085161306 0.092239896
|
||||
0.087303454 0.085655788
|
||||
0.089620471 0.081615114
|
||||
0.091631467 0.076984785
|
||||
0.093817332 0.076157254
|
||||
0.096090632 0.070972035
|
||||
0.09823278 0.060551277
|
||||
0.100287493 0.071366358
|
||||
0.102691945 0.065426801
|
||||
0.104834093 0.059045344
|
||||
0.107019958 0.056229269
|
||||
0.109162106 0.061290519
|
||||
0.111347971 0.051160701
|
||||
0.113490119 0.053836338
|
||||
0.115763418 0.05839439
|
||||
0.117949284 0.057514252
|
||||
0.120091431 0.047536765
|
||||
0.122364731 0.04877829
|
||||
0.300000000 0.000000000
|
||||
@@ -1,58 +0,0 @@
|
||||
4.37173E-05 0.224586133
|
||||
0.001967279 0.518057252
|
||||
0.004240578 0.775672976
|
||||
0.006513878 0.825879049
|
||||
0.008656026 0.899985882
|
||||
0.010841891 0.967789786
|
||||
0.012984039 0.99196296
|
||||
0.015038752 0.992909754
|
||||
0.017443204 0.994177124
|
||||
0.019497917 0.991424846
|
||||
0.021771217 0.993887344
|
||||
0.023913365 0.998648926
|
||||
0.026317816 0.999361915
|
||||
0.028241378 0.998678679
|
||||
0.030514678 0.998109885
|
||||
0.032700543 0.997772638
|
||||
0.034711539 0.997392554
|
||||
0.03711599 0.999600524
|
||||
0.039301855 0.994779551
|
||||
0.041444003 0.996274647
|
||||
0.043629868 0.996149207
|
||||
0.045772016 0.998305774
|
||||
0.047957881 0.99658353
|
||||
0.050231181 0.995950077
|
||||
0.052373329 0.995563046
|
||||
0.054559194 0.995479284
|
||||
0.056701342 0.993347667
|
||||
0.058974642 0.990508095
|
||||
0.061160507 0.990476431
|
||||
0.063171503 0.988754593
|
||||
0.06548852 0.987959462
|
||||
0.067630668 0.987058572
|
||||
0.069903968 0.988939946
|
||||
0.072089833 0.986821236
|
||||
0.074231981 0.985149296
|
||||
0.076417846 0.987362592
|
||||
0.078428842 0.986507086
|
||||
0.080833293 0.98295765
|
||||
0.082888007 0.983313364
|
||||
0.085161306 0.983817788
|
||||
0.087303454 0.983554288
|
||||
0.089620471 0.979763465
|
||||
0.091631467 0.975348682
|
||||
0.093817332 0.975887375
|
||||
0.096090632 0.972236927
|
||||
0.09823278 0.965515882
|
||||
0.100287493 0.962852402
|
||||
0.102691945 0.956266464
|
||||
0.104834093 0.953509902
|
||||
0.107019958 0.945236222
|
||||
0.109162106 0.808306574
|
||||
0.111347971 0.186001557
|
||||
0.113490119 0.078929784
|
||||
0.115763418 0.051976893
|
||||
0.117949284 0.05338487
|
||||
0.120091431 0.04844847
|
||||
0.122364731 0.053399168
|
||||
0.300000000 0.000000000
|
||||
Binary file not shown.
@@ -0,0 +1,93 @@
|
||||
/*******************************************************************************
|
||||
*
|
||||
* McStas, neutron ray-tracing package
|
||||
* Copyright (C) 1997-2011, All rights reserved
|
||||
* Risoe National Laboratory, Roskilde, Denmark
|
||||
* Institut Laue Langevin, Grenoble, France
|
||||
*
|
||||
* Component: ScanningSlit
|
||||
*
|
||||
* %I
|
||||
* Written by: Jochen Stahn, based on the Slip component by
|
||||
* Kim Lefmann and Henrik M. Roennow
|
||||
* Date: 15. 10. 2013
|
||||
* Version: $Revision: 0.00 $
|
||||
* Origin: PSI
|
||||
* Release: McStas 1.12c
|
||||
*
|
||||
* moving rectangular slit with optional insignificance cut
|
||||
*
|
||||
* %D
|
||||
* A simple rectangular slit, where the blades in x-direction move with time.
|
||||
* No transmission around the slit is allowed.
|
||||
* If cutting option is used, low-weight neutron rays are ABSORBED
|
||||
*
|
||||
* Example: FastSlit(height=0.1)
|
||||
*
|
||||
* %P
|
||||
* INPUT PARAMETERS
|
||||
*
|
||||
* height: height of the opening in y, centered.
|
||||
*
|
||||
* Optional parameters:
|
||||
* cut: Lower limit for allowed weight (1)
|
||||
*
|
||||
* %E
|
||||
*******************************************************************************/
|
||||
|
||||
|
||||
DEFINE COMPONENT ScanningSlit
|
||||
DEFINITION PARAMETERS ()
|
||||
SETTING PARAMETERS (xmin=0, xmax=0, ymin=0, ymax=0, height=0.1, dist=2.4, lmin=4, lmax=10, omeg=0, eps=1.25, reso=0.01, dthet=1.5, posi=33.7)
|
||||
OUTPUT PARAMETERS ()
|
||||
/* STATE PARAMETERS (x,y,z,vx,vy,vz,t,s1,s2,p) */
|
||||
|
||||
DECLARE
|
||||
%{
|
||||
double lact;
|
||||
%}
|
||||
|
||||
|
||||
INITIALIZE
|
||||
%{
|
||||
if (height > 0) { ymax=height/2; ymin=-ymax; }
|
||||
lmax *= (1+reso/2) ;
|
||||
lmin *= (1-reso/2) ;
|
||||
eps *= PI/180 ;
|
||||
omeg *= PI/180 ;
|
||||
dthet *= PI/180;
|
||||
%}
|
||||
|
||||
TRACE
|
||||
%{
|
||||
PROP_Z0;
|
||||
lact = 3.956e3*t/posi ;
|
||||
xmin = -dist * (1+reso/2)*( eps + 0.5*dthet - dthet*(lact-lmin)/(lmax-lmin) ) ;
|
||||
xmax = -dist * (1-reso/2)*( eps + 0.5*dthet - dthet*(lact-lmin)/(lmax-lmin) ) ;
|
||||
if ((x<xmin || x>xmax || y<ymin || y>ymax))
|
||||
ABSORB;
|
||||
else
|
||||
SCATTER;
|
||||
%}
|
||||
|
||||
MCDISPLAY
|
||||
%{
|
||||
magnify("xy");
|
||||
double xw, yh;
|
||||
xw = (xmax - xmin)/2.0;
|
||||
yh = (ymax - ymin)/2.0;
|
||||
multiline(3, xmin-xw, (double)ymax, 0.0,
|
||||
(double)xmin, (double)ymax, 0.0,
|
||||
(double)xmin, ymax+yh, 0.0);
|
||||
multiline(3, xmax+xw, (double)ymax, 0.0,
|
||||
(double)xmax, (double)ymax, 0.0,
|
||||
(double)xmax, ymax+yh, 0.0);
|
||||
multiline(3, xmin-xw, (double)ymin, 0.0,
|
||||
(double)xmin, (double)ymin, 0.0,
|
||||
(double)xmin, ymin-yh, 0.0);
|
||||
multiline(3, xmax+xw, (double)ymin, 0.0,
|
||||
(double)xmax, (double)ymin, 0.0,
|
||||
(double)xmax, ymin-yh, 0.0);
|
||||
%}
|
||||
|
||||
END
|
||||
@@ -0,0 +1,145 @@
|
||||
/*******************************************************************************
|
||||
* McStas instrument definition URL=http://www.mcstas.org
|
||||
*
|
||||
* Instrument: Selene_geometry
|
||||
*
|
||||
* %Identification
|
||||
* Written by: Artur Glavic (artur.glavic@psi.ch); Jochen Stahn (jochen.stahn@psi.ch); Christine Klauser (christine.klauser@psi.ch)
|
||||
* Date: 01. 11. 2017
|
||||
* Origin: PSI
|
||||
* Release: McStas 2.4.1
|
||||
* Version: 1.0
|
||||
* %INSTRUMENT_SITE: Estia (E02)
|
||||
*
|
||||
* Estia is a vertical sample, focusing reflectometer for small sample
|
||||
*
|
||||
* %Description
|
||||
* This is a test model to analyze the impact of the displacement of the two Selene guide
|
||||
* mirrors to each other and the virtual source. With this it is possible to evaluate
|
||||
* necessary engineering requirements for the guide carriers.
|
||||
*
|
||||
* %Parameters
|
||||
* sample_length: [m] Size of sample in beam direction, also controls virtual source opening
|
||||
* sample_height: [m] Size of sample in vertical direction, also controls virtual source opening
|
||||
* omegaa: [deg] sample rotation omega
|
||||
*
|
||||
* %End
|
||||
*******************************************************************************/
|
||||
|
||||
DEFINE INSTRUMENT Selene_geometry(double sample_length=0.01, double sample_height=0.01, omegaa=1.2,
|
||||
double tx_1=0.0, double tz_1=0.0, double rz_1=0.0, double ry_1=0.0,
|
||||
double tx_2=0.0, double tz_2=0.0, double rz_2=0.0, double ry_2=0.0
|
||||
)
|
||||
|
||||
DECLARE
|
||||
%{
|
||||
|
||||
%}
|
||||
|
||||
INITIALIZE
|
||||
%{
|
||||
|
||||
%}
|
||||
TRACE
|
||||
|
||||
COMPONENT origin = Progress_bar()
|
||||
AT (0,0,0) ABSOLUTE
|
||||
|
||||
COMPONENT ISCS = Arm()
|
||||
AT (0, 0, 0) RELATIVE origin
|
||||
|
||||
COMPONENT arm_selene1 = Arm()
|
||||
AT (tx_1, 0, tz_1) RELATIVE ISCS
|
||||
ROTATED (0, ry_1, rz_1) RELATIVE ISCS
|
||||
|
||||
COMPONENT arm_selene2 = Arm()
|
||||
AT (tx_2, 0, 2*selene_c+tz_2) RELATIVE ISCS
|
||||
ROTATED (0, ry_2, rz_2) RELATIVE ISCS
|
||||
|
||||
COMPONENT arm_virtual_source_beam = Arm()
|
||||
AT (0, 0, 0) RELATIVE arm_selene1
|
||||
ROTATED (0, 1.25, 0) RELATIVE ISCS
|
||||
// Axes used for the virtual source masks, this depends on the sample angle to achieve the correct footprint
|
||||
COMPONENT arm_virtual_source = Arm()
|
||||
AT (0, 0, 0) RELATIVE arm_virtual_source_beam
|
||||
ROTATED (0, omegaa, 0) RELATIVE arm_virtual_source_beam
|
||||
|
||||
|
||||
COMPONENT source = Moderator(radius = 0.01, focus_xw = selene_entry+0.002, focus_yh = selene_entry+0.002,
|
||||
target_index=3, Emin=0.75, Emax=6)
|
||||
AT (0,0,-0.1) RELATIVE ISCS
|
||||
ROTATED (-1.25, 1.25, 0) RELATIVE ISCS
|
||||
|
||||
/* The actual virtual source mask, two L-shaped absorbers (first top-right) */
|
||||
COMPONENT virtual_source_TR = Slit(
|
||||
xmin = 0.0, xmax = 1.0, ymin = -1.0, ymax = sample_height/2)
|
||||
AT (0, 0, -0.5*sample_length) RELATIVE arm_virtual_source
|
||||
|
||||
/* The actual virtual source mask, two L-shaped absorbers (second bottom-left) */
|
||||
COMPONENT virtual_source_BL = Slit(
|
||||
xmin = -1.0, xmax = 0.0, ymin = -sample_height/2, ymax = 1.0)
|
||||
AT (0, 0, 0.5*sample_length) RELATIVE arm_virtual_source
|
||||
|
||||
%include "Estia_selene.instr"
|
||||
|
||||
COMPONENT M_n3 = Monitor_nD(
|
||||
filename = "M_n3",
|
||||
options = "x limits=[-0.05 0.05] bins=1000 y limits=[-0.05 0.05] bins=1000 lambda limits=[0 35] bins=350",
|
||||
xwidth = 0.5, yheight = 0.5)
|
||||
AT (0, 0, 4*selene_c-0.01) RELATIVE ISCS
|
||||
ROTATED (0, 1.25+omegaa, 0) RELATIVE ISCS
|
||||
|
||||
COMPONENT M_n2 = Monitor_nD(
|
||||
filename = "M_n2",
|
||||
options = "x limits=[-0.05 0.05] bins=1000 y limits=[-0.05 0.05] bins=1000 lambda limits=[0 35] bins=350",
|
||||
xwidth = 0.5, yheight = 0.5)
|
||||
AT (0, 0, 4*selene_c-0.005) RELATIVE ISCS
|
||||
ROTATED (0, 1.25+omegaa, 0) RELATIVE ISCS
|
||||
|
||||
COMPONENT M_n1 = Monitor_nD(
|
||||
filename = "M_n1",
|
||||
options = "x limits=[-0.05 0.05] bins=1000 y limits=[-0.05 0.05] bins=1000 lambda limits=[0 35] bins=350",
|
||||
xwidth = 0.5, yheight = 0.5)
|
||||
AT (0, 0, 4*selene_c-0.0025) RELATIVE ISCS
|
||||
ROTATED (0, 1.25+omegaa, 0) RELATIVE ISCS
|
||||
|
||||
COMPONENT M_0 = Monitor_nD(
|
||||
filename = "M_0",
|
||||
options = "x limits=[-0.05 0.05] bins=1000 y limits=[-0.05 0.05] bins=1000 lambda limits=[0 35] bins=350",
|
||||
xwidth = 0.5, yheight = 0.5)
|
||||
AT (0, 0, 4*selene_c) RELATIVE ISCS
|
||||
ROTATED (0, 1.25+omegaa, 0) RELATIVE ISCS
|
||||
|
||||
COMPONENT MFP = Monitor_nD(
|
||||
filename = "MFP",
|
||||
options = "x limits=[-0.05 0.05] bins=1000 y limits=[-0.05 0.05] bins=1000 lambda limits=[0 35] bins=350",
|
||||
xwidth = 0.5, yheight = 0.5)
|
||||
AT (0, 0, 4*selene_c+0.0001) RELATIVE ISCS
|
||||
ROTATED (0, 1.25+omegaa-90, 0) RELATIVE ISCS
|
||||
|
||||
COMPONENT M_p1 = Monitor_nD(
|
||||
filename = "M_p1",
|
||||
options = "x limits=[-0.05 0.05] bins=1000 y limits=[-0.05 0.05] bins=1000 lambda limits=[0 35] bins=350",
|
||||
xwidth = 0.5, yheight = 0.5)
|
||||
AT (0, 0, 4*selene_c+0.0025) RELATIVE ISCS
|
||||
ROTATED (0, 1.25+omegaa, 0) RELATIVE ISCS
|
||||
|
||||
COMPONENT M_p2 = Monitor_nD(
|
||||
filename = "M_p2",
|
||||
options = "x limits=[-0.05 0.05] bins=1000 y limits=[-0.05 0.05] bins=1000 lambda limits=[0 35] bins=350",
|
||||
xwidth = 0.5, yheight = 0.5)
|
||||
AT (0, 0, 4*selene_c+0.005) RELATIVE ISCS
|
||||
ROTATED (0, 1.25+omegaa, 0) RELATIVE ISCS
|
||||
|
||||
COMPONENT M_p3 = Monitor_nD(
|
||||
filename = "M_p3",
|
||||
options = "x limits=[-0.05 0.05] bins=1000 y limits=[-0.05 0.05] bins=1000 lambda limits=[0 35] bins=350",
|
||||
xwidth = 0.5, yheight = 0.5)
|
||||
AT (0, 0, 4*selene_c+0.01) RELATIVE ISCS
|
||||
ROTATED (0, 1.25+omegaa, 0) RELATIVE ISCS
|
||||
|
||||
FINALLY
|
||||
%{
|
||||
%}
|
||||
|
||||
END
|
||||
@@ -1,6 +1,4 @@
|
||||
#!/bin/bash
|
||||
export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:/afs/psi.ch/project/sinq/sl6-64/mcstas2.4/mcstas/2.4/libs
|
||||
|
||||
|
||||
if [ Estia_baseline.instr -nt Estia_baseline.out ] || [ ! -f Estia_baseline.out ] \
|
||||
|| [ Estia_feeder.instr -nt Estia_baseline.out ] \
|
||||
@@ -8,6 +6,6 @@ if [ Estia_baseline.instr -nt Estia_baseline.out ] || [ ! -f Estia_baseline.out
|
||||
|| [ Estia_mf.instr -nt Estia_baseline.out ] \
|
||||
|| [ Estia_selene2.instr -nt Estia_baseline.out ]; then
|
||||
rm Estia_baseline.c Estia_baseline.out
|
||||
mcstas -o Estia_baseline.c Estia_baseline.instr --trace
|
||||
mpicc -O2 -o Estia_baseline.out Estia_baseline.c -lm -DUSE_MPI -DUSE_NEXUS -lNeXus
|
||||
mcstas -o Estia_baseline.c Estia_baseline.instr
|
||||
mpicc -O3 -o Estia_baseline.out Estia_baseline.c -lm -DUSE_MPI
|
||||
fi
|
||||
|
||||
@@ -1,3 +0,0 @@
|
||||
mcc05:72
|
||||
mcc06:72
|
||||
|
||||
Executable
+1002
File diff suppressed because it is too large
Load Diff
Executable
+1002
File diff suppressed because it is too large
Load Diff
Executable
+102
@@ -0,0 +1,102 @@
|
||||
0.00 1.000
|
||||
0.01 0.0099682
|
||||
0.0149 0.0042479
|
||||
0.0198 0.0022557
|
||||
0.0247 0.0013447
|
||||
0.0296 0.00086176
|
||||
0.0345 0.000582
|
||||
0.0394 0.0004107
|
||||
0.0443 0.00030205
|
||||
0.0492 0.00023151
|
||||
0.0541 0.0001849
|
||||
0.059 0.00015356
|
||||
0.0639 0.00013197
|
||||
0.0688 0.00011657
|
||||
0.0737 0.00010502
|
||||
0.0786 9.58E-05
|
||||
0.0835 8.80E-05
|
||||
0.0884 8.11E-05
|
||||
0.0933 7.47E-05
|
||||
0.0982 6.88E-05
|
||||
0.1031 6.33E-05
|
||||
0.108 5.82E-05
|
||||
0.1129 5.35E-05
|
||||
0.1178 4.94E-05
|
||||
0.1227 4.57E-05
|
||||
0.1276 4.24E-05
|
||||
0.1325 3.96E-05
|
||||
0.1374 3.71E-05
|
||||
0.1423 3.48E-05
|
||||
0.1472 3.27E-05
|
||||
0.1521 3.08E-05
|
||||
0.157 2.89E-05
|
||||
0.1619 2.71E-05
|
||||
0.1668 2.53E-05
|
||||
0.1717 2.35E-05
|
||||
0.1766 2.17E-05
|
||||
0.1815 2.00E-05
|
||||
0.1864 1.83E-05
|
||||
0.1913 1.66E-05
|
||||
0.1962 1.51E-05
|
||||
0.2011 1.36E-05
|
||||
0.206 1.22E-05
|
||||
0.2109 1.08E-05
|
||||
0.2158 9.61E-06
|
||||
0.2207 8.46E-06
|
||||
0.2256 7.39E-06
|
||||
0.2305 6.39E-06
|
||||
0.2354 5.46E-06
|
||||
0.2403 4.61E-06
|
||||
0.2452 3.84E-06
|
||||
0.2501 3.16E-06
|
||||
0.255 2.55E-06
|
||||
0.2599 2.04E-06
|
||||
0.2648 1.61E-06
|
||||
0.2697 1.26E-06
|
||||
0.2746 9.93E-07
|
||||
0.2795 7.98E-07
|
||||
0.2844 6.64E-07
|
||||
0.2893 5.82E-07
|
||||
0.2942 5.39E-07
|
||||
0.2991 5.24E-07
|
||||
0.304 5.28E-07
|
||||
0.3089 5.41E-07
|
||||
0.3138 5.58E-07
|
||||
0.3187 5.73E-07
|
||||
0.3236 5.83E-07
|
||||
0.3285 5.88E-07
|
||||
0.3334 5.88E-07
|
||||
0.3383 5.84E-07
|
||||
0.3432 5.77E-07
|
||||
0.3481 5.69E-07
|
||||
0.353 5.62E-07
|
||||
0.3579 5.57E-07
|
||||
0.3628 5.53E-07
|
||||
0.3677 5.52E-07
|
||||
0.3726 5.53E-07
|
||||
0.3775 5.56E-07
|
||||
0.3824 5.61E-07
|
||||
0.3873 5.67E-07
|
||||
0.3922 5.74E-07
|
||||
0.3971 5.81E-07
|
||||
0.402 5.89E-07
|
||||
0.4069 5.96E-07
|
||||
0.4118 6.04E-07
|
||||
0.4167 6.12E-07
|
||||
0.4216 6.19E-07
|
||||
0.4265 6.27E-07
|
||||
0.4314 6.35E-07
|
||||
0.4363 6.43E-07
|
||||
0.4412 6.50E-07
|
||||
0.4461 6.58E-07
|
||||
0.451 6.66E-07
|
||||
0.4559 6.74E-07
|
||||
0.4608 6.82E-07
|
||||
0.4657 6.89E-07
|
||||
0.4706 6.97E-07
|
||||
0.4755 7.04E-07
|
||||
0.4804 7.11E-07
|
||||
0.4853 7.16E-07
|
||||
0.4902 7.20E-07
|
||||
0.4951 7.21E-07
|
||||
0.5 7.23E-07
|
||||
Binary file not shown.
@@ -0,0 +1,12 @@
|
||||
#!/bin/bash
|
||||
|
||||
mcstas -o Estia_monitor.c Estia_monitor.instr
|
||||
mpicc -O3 -o Estia_monitor.out Estia_monitor.c -lm -DUSE_MPI
|
||||
|
||||
|
||||
mpirun -np 6 Estia_monitor.out --ncount=1e8 --dir=../results/monitor --gravitation enable_windows=1 lambda_start=1.0 lambda_end=35.0 direct_beam=0
|
||||
mpirun -np 6 Estia_monitor.out --ncount=1e8 --dir=../results/monitor_ref --gravitation enable_windows=1 lambda_start=1.0 lambda_end=35.0 direct_beam=2
|
||||
mpirun -np 6 Estia_monitor.out --ncount=1e8 --dir=../results/monitor_trans --gravitation enable_windows=1 lambda_start=1.0 lambda_end=35.0 direct_beam=1
|
||||
|
||||
mpirun -np 6 Estia_monitor.out --ncount=1e8 --dir=../results/monitor100 --gravitation enable_windows=1 lambda_start=1.0 lambda_end=35.0 direct_beam=0 foil_thickness=0.0001
|
||||
mpirun -np 6 Estia_monitor.out --ncount=1e8 --dir=../results/monitor001 --gravitation enable_windows=1 lambda_start=1.0 lambda_end=35.0 direct_beam=0 foil_thickness=0.000001
|
||||
Executable
+213
@@ -0,0 +1,213 @@
|
||||
#!/bin/bash
|
||||
|
||||
DEST=../results
|
||||
ncount=2e7
|
||||
use_cores=4
|
||||
sample=4
|
||||
omega=0.8
|
||||
sample_length=0.01
|
||||
sample_height=0.01
|
||||
lambda_start=3.0
|
||||
lambda_end=32.0
|
||||
|
||||
###################### Brilliance Transfer 10x10mm² and 1x1mm² VS ####################
|
||||
if [ Estia_baseline.instr -nt Estia_baseline.out ]; then
|
||||
rm Estia_baseline.c Estia_baseline.out
|
||||
mcstas -o Estia_baseline.c Estia_baseline.instr
|
||||
mpicc -O3 -o Estia_baseline.out Estia_baseline.c -lm -DUSE_MPI -DUSE_NEXUS -lNeXus
|
||||
fi
|
||||
#
|
||||
# if [ Estia_baseline_ana1.instr -nt Estia_baseline_ana1.out ]; then
|
||||
# rm Estia_baseline_ana1.c Estia_baseline_ana1.out
|
||||
# mcstas -o Estia_baseline_ana1.c Estia_baseline_ana1.instr
|
||||
# mpicc -O3 -o Estia_baseline_ana1.out Estia_baseline_ana1.c -lm -DUSE_MPI -DUSE_NEXUS -lNeXus
|
||||
# fi
|
||||
#
|
||||
# if [ Estia_baseline_ana2.instr -nt Estia_baseline_ana2.out ]; then
|
||||
# rm Estia_baseline_ana2.c Estia_baseline_ana2.out
|
||||
# mcstas -o Estia_baseline_ana2.c Estia_baseline_ana2.instr
|
||||
# mpicc -O3 -o Estia_baseline_ana2.out Estia_baseline_ana2.c -lm -DUSE_MPI -DUSE_NEXUS -lNeXus
|
||||
# fi
|
||||
|
||||
|
||||
###################### Reference and Ni-layer measurement 10x10mm² sample ####################
|
||||
# ncount=1e10
|
||||
sample_length=0.01
|
||||
sample_height=0.01
|
||||
lambda_start=3.5
|
||||
lambda_end=30.0
|
||||
sample=1
|
||||
|
||||
# omega=1.0
|
||||
# DESTi=$DEST/pol_ref_10x10_10
|
||||
# if [ -e "$DESTi" ]; then
|
||||
# rm -r "$DESTi"
|
||||
# fi
|
||||
#
|
||||
# mpirun -np $use_cores Estia_baseline.out \
|
||||
# --dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
# omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.000 \
|
||||
# sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
# lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=0 enable_polarizer=0 enable_analyzer=0
|
||||
#
|
||||
# DESTi=$DEST/pol1_10x10_10
|
||||
# if [ -e "$DESTi" ]; then
|
||||
# rm -r "$DESTi"
|
||||
# fi
|
||||
#
|
||||
# mpirun -np $use_cores Estia_baseline.out \
|
||||
# --dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
# omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.000 \
|
||||
# sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
# lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=0 enable_polarizer=1 enable_analyzer=0
|
||||
#
|
||||
# DESTi=$DEST/pol2_10x10_10
|
||||
# if [ -e "$DESTi" ]; then
|
||||
# rm -r "$DESTi"
|
||||
# fi
|
||||
#
|
||||
# mpirun -np $use_cores Estia_baseline.out \
|
||||
# --dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
# omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.000 \
|
||||
# sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
# lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=0 enable_polarizer=2 enable_analyzer=0
|
||||
# #
|
||||
# omega=7.0
|
||||
# DESTi=$DEST/pol1_10x10_70
|
||||
# if [ -e "$DESTi" ]; then
|
||||
# rm -r "$DESTi"
|
||||
# fi
|
||||
#
|
||||
# mpirun -np $use_cores Estia_baseline.out \
|
||||
# --dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
# omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.000 \
|
||||
# sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
# lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=0 enable_polarizer=1 enable_analyzer=0
|
||||
#
|
||||
# DESTi=$DEST/pol2_10x10_70
|
||||
# if [ -e "$DESTi" ]; then
|
||||
# rm -r "$DESTi"
|
||||
# fi
|
||||
#
|
||||
# mpirun -np $use_cores Estia_baseline.out \
|
||||
# --dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
# omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.000 \
|
||||
# sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
# lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=0 enable_polarizer=2 enable_analyzer=0
|
||||
#
|
||||
|
||||
# ncount=1e8
|
||||
omega=0.0
|
||||
# DESTi=$DEST/pol1_10x50_70
|
||||
# if [ -e "$DESTi" ]; then
|
||||
# rm -r "$DESTi"
|
||||
# fi
|
||||
#
|
||||
# mpirun -np $use_cores Estia_baseline.out \
|
||||
# --dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
# omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.0025 \
|
||||
# sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
# lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=0 enable_polarizer=1 enable_analyzer=0
|
||||
|
||||
DESTi=$DEST/pol2_10x50_70
|
||||
if [ -e "$DESTi" ]; then
|
||||
rm -r "$DESTi"
|
||||
fi
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.0025 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=0 enable_polarizer=2 enable_analyzer=0
|
||||
#
|
||||
# ncount=5e9
|
||||
# sample_length=0.01
|
||||
# omega=1.0
|
||||
# sample_length=0.003
|
||||
# DESTi=$DEST/pol_ref_3x10_10
|
||||
# if [ -e "$DESTi" ]; then
|
||||
# rm -r "$DESTi"
|
||||
# fi
|
||||
#
|
||||
# mpirun -np $use_cores Estia_baseline.out \
|
||||
# --dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
# omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.000 \
|
||||
# sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
# lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=0 enable_polarizer=0 enable_analyzer=0
|
||||
#
|
||||
# DESTi=$DEST/pol1_3x10_10
|
||||
# if [ -e "$DESTi" ]; then
|
||||
# rm -r "$DESTi"
|
||||
# fi
|
||||
#
|
||||
# mpirun -np $use_cores Estia_baseline.out \
|
||||
# --dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
# omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.000 \
|
||||
# sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
# lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=0 enable_polarizer=1 enable_analyzer=0
|
||||
#
|
||||
# DESTi=$DEST/pol2_3x10_10
|
||||
# if [ -e "$DESTi" ]; then
|
||||
# rm -r "$DESTi"
|
||||
# fi
|
||||
#
|
||||
# mpirun -np $use_cores Estia_baseline.out \
|
||||
# --dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
# omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.000 \
|
||||
# sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
# lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=0 enable_polarizer=2 enable_analyzer=0
|
||||
#
|
||||
#
|
||||
# sample_height=0.01
|
||||
# sample_length=0.01
|
||||
# DESTi=$DEST/ana1_10x10_10
|
||||
# if [ -e "$DESTi" ]; then
|
||||
# rm -r "$DESTi"
|
||||
# fi
|
||||
#
|
||||
# mpirun -np $use_cores Estia_baseline.out \
|
||||
# --dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
# omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.000 \
|
||||
# sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
# lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=0 enable_polarizer=0 enable_analyzer=1
|
||||
|
||||
# sample_height=0.001
|
||||
# DESTi=$DEST/ana1_10x1_10
|
||||
# if [ -e "$DESTi" ]; then
|
||||
# rm -r "$DESTi"
|
||||
# fi
|
||||
#
|
||||
# mpirun -np $use_cores Estia_baseline_ana1.out \
|
||||
# --dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
# omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.000 \
|
||||
# sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
# lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=0 enable_polarizer=0 enable_analyzer=1
|
||||
#
|
||||
# sample_height=0.01
|
||||
# DESTi=$DEST/ana2_10x10_10
|
||||
# if [ -e "$DESTi" ]; then
|
||||
# rm -r "$DESTi"
|
||||
# fi
|
||||
#
|
||||
# mpirun -np $use_cores Estia_baseline.out \
|
||||
# --dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
# omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.000 \
|
||||
# sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
# lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=0 enable_polarizer=0 enable_analyzer=2
|
||||
|
||||
# sample_height=0.001
|
||||
# DESTi=$DEST/ana2_10x1_10
|
||||
# if [ -e "$DESTi" ]; then
|
||||
# rm -r "$DESTi"
|
||||
# fi
|
||||
#
|
||||
# mpirun -np $use_cores Estia_baseline_ana2.out \
|
||||
# --dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
# omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.000 \
|
||||
# sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
# lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=0 enable_polarizer=0 enable_analyzer=2
|
||||
|
||||
|
||||
###################### Reference and Ni-layer measurement 1x1mm² sample ####################
|
||||
|
||||
|
||||
Executable
+52
@@ -0,0 +1,52 @@
|
||||
#!/usr/bin/env python
|
||||
#-*- coding: utf8 -*-
|
||||
|
||||
import sys, os
|
||||
from subprocess import call
|
||||
from numpy import *
|
||||
from scipy.optimize import leastsq
|
||||
|
||||
|
||||
CALL='/afs/psi.ch/project/sinq/sl6-64/bin/mpirun -np $SLURM_NPROCS Estia_baseline.out '
|
||||
CALL=+'--dir=../results/selene_geo --ncount=1e9 '
|
||||
CALL=+'omegaa=1.0 sample=4 sample_length=0.001 sample_height=0.01 '
|
||||
CALL=+'lambda_start=2.5 lambda_end=15.0 enable_gravity=1 enable_chopper=1 '
|
||||
CALL=+'lambda_min=3.75 selene1_foot1y=%.4f selene1_foot2y=%.4f '
|
||||
|
||||
def B(x,w):
|
||||
# box function with full width w
|
||||
return float32(abs(x)<=(w/2.))
|
||||
|
||||
def G(x,I0,x0,sigma):
|
||||
# Gaussian with intensity I0, center x0 and standard deviation sigma
|
||||
return I0*exp(-0.5*(x-x0)**2/sigma**2)
|
||||
|
||||
def Intensity(x, p):
|
||||
I0, x0, sigma, w=p
|
||||
return convolve(B(xc,w), G(xc, I0, x0, sigma), mode='same')
|
||||
|
||||
def Beam(x,I0,x0,sigma,w):
|
||||
I=I0*where((x-x0)<(-w/2.), exp(-0.5*(x-x0+w/2.)**2/sigma**2),
|
||||
where((x-x0)>(w/2.), exp(-0.5*(x-x0-w/2.)**2/sigma**2), 1.))
|
||||
return I
|
||||
|
||||
def residuals(p, x, y):
|
||||
I0, x0, sigma, w=p
|
||||
return y-Beam(x, I0, x0, sigma, w)
|
||||
|
||||
def FWHM(pi):
|
||||
rng=xc[where(Beam(xc, *pi)>=(pi[0]/2.))[0]]
|
||||
return rng[-1]-rng[0]
|
||||
|
||||
x=linspace(-10.0005, 10.0005, 20001)
|
||||
xc=(x[1:]+x[:-1])/2.
|
||||
|
||||
def analyze(fi):
|
||||
sim=mr.McSim(fi)
|
||||
det=sim['tof_sample']
|
||||
ignore,y=det.project1d('x', bins=x/1000.)
|
||||
p,res=leastsq(residuals, (y.max(), (xc*y).sum()/y.sum(), 0.01, 0.1), (xc, y))
|
||||
return p,det
|
||||
|
||||
if __name__=='__main__':
|
||||
pass
|
||||
@@ -0,0 +1,86 @@
|
||||
#!/bin/bash
|
||||
|
||||
DEST=../results
|
||||
use_cores=6
|
||||
|
||||
###################### Compile model if necessary ####################
|
||||
bash compile_if_neede.sh
|
||||
|
||||
######## Reference and Ni-layer conventional measurement 10x10mm² sample #############
|
||||
ncount=1e10
|
||||
sample_length=0.01
|
||||
sample_height=0.01
|
||||
lambda_start=3.25
|
||||
lambda_min=3.75
|
||||
lambda_end=12.75
|
||||
sample=1
|
||||
|
||||
omega=1.0
|
||||
DESTi=$DEST/tof_reference_10x10_10
|
||||
if [ -e "$DESTi" ]; then
|
||||
rm -r "$DESTi"
|
||||
fi
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --format=NeXuS --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.01 over_illumination=0.0001 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_min=$lambda_min lambda_end=$lambda_end enable_gravity=1 enable_chopper=1
|
||||
|
||||
# Ni Sample
|
||||
ncount=2e10
|
||||
sample=2
|
||||
omega=0.5
|
||||
DESTi=$DEST/tof_nickle_10x10_05
|
||||
if [ -e "$DESTi" ]; then
|
||||
rm -r "$DESTi"
|
||||
fi
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --format=NeXuS --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.03 over_illumination=0.0001 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_min=$lambda_min lambda_end=$lambda_end enable_gravity=1 enable_chopper=1
|
||||
|
||||
ncount=1e10
|
||||
omega=1.2
|
||||
DESTi=$DEST/tof_nickle_10x10_12
|
||||
if [ -e "$DESTi" ]; then
|
||||
rm -r "$DESTi"
|
||||
fi
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --format=NeXuS --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.03 over_illumination=0.0001 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_min=$lambda_min lambda_end=$lambda_end enable_gravity=1 enable_chopper=1
|
||||
|
||||
ncount=5e9
|
||||
omega=3.0
|
||||
DESTi=$DEST/tof_nickle_10x10_30
|
||||
if [ -e "$DESTi" ]; then
|
||||
rm -r "$DESTi"
|
||||
fi
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --format=NeXuS --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.03 over_illumination=0.0001 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_min=$lambda_min lambda_end=$lambda_end enable_gravity=1 enable_chopper=1
|
||||
|
||||
ncount=2e9
|
||||
omega=7.5
|
||||
DESTi=$DEST/tof_nickle_10x10_75
|
||||
if [ -e "$DESTi" ]; then
|
||||
rm -r "$DESTi"
|
||||
fi
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --format=NeXuS --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.03 over_illumination=0.0001 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_min=$lambda_min lambda_end=$lambda_end enable_gravity=1 enable_chopper=1
|
||||
|
||||
###################### Reference and Ni-layer measurement 1x1mm² sample ####################
|
||||
|
||||
|
||||
@@ -0,0 +1,46 @@
|
||||
#!/bin/bash
|
||||
|
||||
DEST=../results
|
||||
ncount=1e9
|
||||
use_cores=6
|
||||
sample=4
|
||||
omega=0.8
|
||||
sample_length=0.01
|
||||
sample_height=0.01
|
||||
lambda_start=3.0
|
||||
lambda_end=32.0
|
||||
|
||||
bash compile_if_neede.sh
|
||||
|
||||
|
||||
################# Reference and Ni-layer measurement 2 pulse skipping ####################
|
||||
ncount=4e9
|
||||
sample_length=0.01
|
||||
sample_height=0.01
|
||||
lambda_start=3.9
|
||||
lambda_end=24.0
|
||||
sample=1
|
||||
omega=2.0
|
||||
|
||||
DESTi=$DEST/single_skip_reference
|
||||
if [ -e "$DESTi" ]; then
|
||||
rm -r "$DESTi"
|
||||
fi
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --format=NeXuS --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.0002 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=3
|
||||
|
||||
sample=2
|
||||
DESTi=$DEST/single_skip_nickle
|
||||
if [ -e "$DESTi" ]; then
|
||||
rm -r "$DESTi"
|
||||
fi
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --format=NeXuS --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.0002 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=3
|
||||
@@ -0,0 +1,119 @@
|
||||
#!/bin/bash
|
||||
|
||||
DEST=../results
|
||||
ncount=1e9
|
||||
use_cores=6
|
||||
sample=1
|
||||
omega=15.0
|
||||
sample_length=0.02
|
||||
sample_height=0.01
|
||||
|
||||
lambda_start=3.5
|
||||
lambda_min=5.0
|
||||
lambda_end=15.5
|
||||
frame_usage=0.92
|
||||
|
||||
###################### Brilliance Transfer 10x10mm² and 1x1mm² VS ####################
|
||||
bash compile_if_neede.sh
|
||||
|
||||
###################### Reference and Ni-layer measurement 10x10mm² sample ####################
|
||||
|
||||
frame_usage=0.92
|
||||
DESTi=$DEST/chopper_f1_092
|
||||
if [ -e "$DESTi" ]; then
|
||||
rm -r "$DESTi"
|
||||
fi
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --format=NeXuS --ncount=$ncount \
|
||||
omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.0001 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=0 enable_chopper=1 \
|
||||
frame_usage=$frame_usage lambda_min=$lambda_min
|
||||
|
||||
frame_usage=0.94
|
||||
DESTi=$DEST/chopper_f1_094
|
||||
if [ -e "$DESTi" ]; then
|
||||
rm -r "$DESTi"
|
||||
fi
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --format=NeXuS --ncount=$ncount \
|
||||
omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.0001 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=0 enable_chopper=1 \
|
||||
frame_usage=$frame_usage lambda_min=$lambda_min
|
||||
|
||||
frame_usage=0.96
|
||||
DESTi=$DEST/chopper_f1_096
|
||||
if [ -e "$DESTi" ]; then
|
||||
rm -r "$DESTi"
|
||||
fi
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --format=NeXuS --ncount=$ncount \
|
||||
omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.0001 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=0 enable_chopper=1 \
|
||||
frame_usage=$frame_usage lambda_min=$lambda_min
|
||||
|
||||
frame_usage=0.974
|
||||
DESTi=$DEST/chopper_f1_097s
|
||||
if [ -e "$DESTi" ]; then
|
||||
rm -r "$DESTi"
|
||||
fi
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --format=NeXuS --ncount=$ncount \
|
||||
omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.0001 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=0 enable_chopper=1 \
|
||||
frame_usage=$frame_usage lambda_min=$lambda_min
|
||||
|
||||
frame_usage=0.98
|
||||
DESTi=$DEST/chopper_f1_098
|
||||
if [ -e "$DESTi" ]; then
|
||||
rm -r "$DESTi"
|
||||
fi
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --format=NeXuS --ncount=$ncount \
|
||||
omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.0001 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=0 enable_chopper=1 \
|
||||
frame_usage=$frame_usage lambda_min=$lambda_min
|
||||
|
||||
frame_usage=0.99
|
||||
DESTi=$DEST/chopper_f1_099
|
||||
if [ -e "$DESTi" ]; then
|
||||
rm -r "$DESTi"
|
||||
fi
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --format=NeXuS --ncount=$ncount \
|
||||
omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.0001 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=0 enable_chopper=1 \
|
||||
frame_usage=$frame_usage lambda_min=$lambda_min
|
||||
|
||||
|
||||
frame_usage=0.98
|
||||
lambda_start=2.75
|
||||
lambda_min=3.75
|
||||
lambda_end=12.75
|
||||
DESTi=$DEST/chopper_f1_375A
|
||||
if [ -e "$DESTi" ]; then
|
||||
rm -r "$DESTi"
|
||||
fi
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --format=NeXuS --ncount=$ncount \
|
||||
omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.0001 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=0 enable_chopper=1 \
|
||||
frame_usage=$frame_usage lambda_min=$lambda_min
|
||||
|
||||
|
||||
###################### Reference and Ni-layer measurement 1x1mm² sample ####################
|
||||
|
||||
|
||||
Executable
+251
@@ -0,0 +1,251 @@
|
||||
#!/bin/tcsh
|
||||
#SBATCH -J mcEstia
|
||||
#SBATCH -N 6
|
||||
#SBATCH --ntasks-per-node=24
|
||||
#SBATCH --time=2-00:00:00
|
||||
#SBATCH --mail-type=fail
|
||||
#SBATCH --mail-user=artur.glavic@psi.ch
|
||||
|
||||
#SBATCH -o stdout.log
|
||||
#SBATCH -e stderr.log
|
||||
|
||||
#SBATCH --partition=medium
|
||||
|
||||
echo "Starting at `date`"
|
||||
echo "Running on hosts: $SLURM_NODELIST"
|
||||
echo "Running on $SLURM_NNODES nodes."
|
||||
echo "Running on $SLURM_NPROCS processors."
|
||||
echo "Current working directory is `pwd`"
|
||||
|
||||
module load mcstas
|
||||
bash compile_if_needed.sh
|
||||
|
||||
set DEST=../results
|
||||
set ncount=1e10
|
||||
set use_cores=4
|
||||
set sample=1
|
||||
set sample_length=0.03
|
||||
set sample_height=0.01
|
||||
set lambda_start=3.5
|
||||
set lambda_end=14.0
|
||||
|
||||
|
||||
################## Reference ########################
|
||||
set omega=1.0
|
||||
set DESTi=$DEST/solid_liquid_reference
|
||||
if ( -e "$DESTi" ) then
|
||||
rm -r "$DESTi"
|
||||
endif
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.04 over_illumination=0.0025 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=1 enable_polarizer=0 enable_analyzer=0
|
||||
|
||||
################## H2O ########################
|
||||
set ncount=1e11
|
||||
set sample=5
|
||||
set omega=0.2
|
||||
set DESTi=$DEST/solid_liquid_h2o_02
|
||||
if ( -e "$DESTi" ) then
|
||||
rm -r "$DESTi"
|
||||
endif
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.04 over_illumination=0.0025 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=1 enable_polarizer=0 enable_analyzer=0
|
||||
|
||||
set ncount=1e10
|
||||
set omega=0.5
|
||||
set DESTi=$DEST/solid_liquid_h2o_05
|
||||
if ( -e "$DESTi" ) then
|
||||
rm -r "$DESTi"
|
||||
endif
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.04 over_illumination=0.0025 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=1 enable_polarizer=0 enable_analyzer=0
|
||||
|
||||
set ncount=4e9
|
||||
set omega=1.2
|
||||
set DESTi=$DEST/solid_liquid_h2o_12
|
||||
if ( -e "$DESTi" ) then
|
||||
rm -r "$DESTi"
|
||||
endif
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.04 over_illumination=0.0025 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=1 enable_polarizer=0 enable_analyzer=0
|
||||
|
||||
set ncount=2e9
|
||||
set omega=3.1
|
||||
set DESTi=$DEST/solid_liquid_h2o_31
|
||||
if ( -e "$DESTi" ) then
|
||||
rm -r "$DESTi"
|
||||
endif
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.04 over_illumination=0.0025 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=1 enable_polarizer=0 enable_analyzer=0
|
||||
|
||||
set ncount=1e9
|
||||
set omega=7.8
|
||||
set DESTi=$DEST/solid_liquid_h2o_78
|
||||
if ( -e "$DESTi" ) then
|
||||
rm -r "$DESTi"
|
||||
endif
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.04 over_illumination=0.0025 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=1 enable_polarizer=0 enable_analyzer=0
|
||||
|
||||
|
||||
################## D2O ########################
|
||||
set ncount=1e11
|
||||
set sample=6
|
||||
set omega=0.2
|
||||
set DESTi=$DEST/solid_liquid_d2o_02
|
||||
if ( -e "$DESTi" ) then
|
||||
rm -r "$DESTi"
|
||||
endif
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.04 over_illumination=0.0025 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=1 enable_polarizer=0 enable_analyzer=0
|
||||
|
||||
set ncount=1e10
|
||||
set omega=0.5
|
||||
set DESTi=$DEST/solid_liquid_d2o_05
|
||||
if ( -e "$DESTi" ) then
|
||||
rm -r "$DESTi"
|
||||
endif
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.04 over_illumination=0.0025 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=1 enable_polarizer=0 enable_analyzer=0
|
||||
|
||||
set ncount=4e9
|
||||
set omega=1.2
|
||||
set DESTi=$DEST/solid_liquid_d2o_12
|
||||
if ( -e "$DESTi" ) then
|
||||
rm -r "$DESTi"
|
||||
endif
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.04 over_illumination=0.0025 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=1 enable_polarizer=0 enable_analyzer=0
|
||||
|
||||
set ncount=2e9
|
||||
set omega=3.1
|
||||
set DESTi=$DEST/solid_liquid_d2o_31
|
||||
if ( -e "$DESTi" ) then
|
||||
rm -r "$DESTi"
|
||||
endif
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.04 over_illumination=0.0025 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=1 enable_polarizer=0 enable_analyzer=0
|
||||
|
||||
set ncount=1e9
|
||||
set omega=7.8
|
||||
set DESTi=$DEST/solid_liquid_d2o_78
|
||||
if ( -e "$DESTi" ) then
|
||||
rm -r "$DESTi"
|
||||
endif
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.04 over_illumination=0.0025 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=1 enable_polarizer=0 enable_analyzer=0
|
||||
|
||||
################## SMW ########################
|
||||
set ncount=1e11
|
||||
set sample=7
|
||||
set omega=0.2
|
||||
set DESTi=$DEST/solid_liquid_smw_02
|
||||
if ( -e "$DESTi" ) then
|
||||
rm -r "$DESTi"
|
||||
endif
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.04 over_illumination=0.0025 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=1 enable_polarizer=0 enable_analyzer=0
|
||||
|
||||
set ncount=1e10
|
||||
set sample=5
|
||||
set omega=0.5
|
||||
set DESTi=$DEST/solid_liquid_smw_05
|
||||
if ( -e "$DESTi" ) then
|
||||
rm -r "$DESTi"
|
||||
endif
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.04 over_illumination=0.0025 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=1 enable_polarizer=0 enable_analyzer=0
|
||||
|
||||
set ncount=4e9
|
||||
set sample=5
|
||||
set omega=1.2
|
||||
set DESTi=$DEST/solid_liquid_smw_12
|
||||
if ( -e "$DESTi" ) then
|
||||
rm -r "$DESTi"
|
||||
endif
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.04 over_illumination=0.0025 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=1 enable_polarizer=0 enable_analyzer=0
|
||||
|
||||
set ncount=2e9
|
||||
set sample=5
|
||||
set omega=3.1
|
||||
set DESTi=$DEST/solid_liquid_smw_31
|
||||
if ( -e "$DESTi" ) then
|
||||
rm -r "$DESTi"
|
||||
endif
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.04 over_illumination=0.0025 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=1 enable_polarizer=0 enable_analyzer=0
|
||||
|
||||
set ncount=1e9
|
||||
set sample=5
|
||||
set omega=7.8
|
||||
set DESTi=$DEST/solid_liquid_smw_78
|
||||
if ( -e "$DESTi" ) then
|
||||
rm -r "$DESTi"
|
||||
endif
|
||||
|
||||
mpirun -np $use_cores Estia_baseline.out \
|
||||
--dir="$DESTi" --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=1 theta_resolution=0.04 over_illumination=0.0025 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
lambda_start=$lambda_start lambda_end=$lambda_end enable_gravity=1 enable_chopper=1 enable_polarizer=0 enable_analyzer=0
|
||||
|
||||
Executable
+32
@@ -0,0 +1,32 @@
|
||||
#!/bin/bash
|
||||
|
||||
DEST=../results
|
||||
ncount=1e8
|
||||
use_cores=6
|
||||
sample=1
|
||||
omega=0.8
|
||||
sample_length=0.01
|
||||
sample_height=0.015
|
||||
|
||||
if [ Estia_vs.instr -nt Estia_vs.out ]; then
|
||||
rm Estia_vs.c Estia_vs.out
|
||||
mcstas -t -o Estia_vs.c Estia_vs.instr
|
||||
mpicc -O3 -o Estia_vs.out Estia_vs.c -lm -DUSE_MPI -DUSE_NEXUS -lNeXus
|
||||
fi
|
||||
|
||||
|
||||
omega=2.0
|
||||
DESTi=$DEST/vs_data
|
||||
if [ -e "$DESTi" ]; then
|
||||
rm -r "$DESTi"
|
||||
fi
|
||||
|
||||
mpirun -np $use_cores Estia_vs.out \
|
||||
--dir="$DESTi" --format=NeXuS --ncount=$ncount --gravitation \
|
||||
omegaa=$omega operationmode=0 theta_resolution=0.04 over_illumination=0.0025 \
|
||||
sample=$sample sample_length=$sample_length sample_height=$sample_height \
|
||||
enable_gravity=1 enable_chopper=0 source_power=5
|
||||
|
||||
###################### Reference and Ni-layer measurement 1x1mm² sample ####################
|
||||
|
||||
|
||||
@@ -0,0 +1,40 @@
|
||||
#!/bin/bash
|
||||
#SBATCH -J mcEstia_1
|
||||
#SBATCH -N 2
|
||||
#SBATCH --ntasks-per-node=24
|
||||
#SBATCH --time=1-00:00:00
|
||||
#SBATCH --mail-type=fail
|
||||
#SBATCH --mail-user=artur.glavic@psi.ch
|
||||
|
||||
#SBATCH -o stdout.log
|
||||
#SBATCH -e stderr.log
|
||||
|
||||
#SBATCH --partition=ll_long
|
||||
|
||||
echo "Starting at `date`"
|
||||
echo "Running on hosts: $SLURM_NODELIST"
|
||||
echo "Running on $SLURM_NNODES nodes."
|
||||
echo "Running on $SLURM_NPROCS processors."
|
||||
echo "Current working directory is `pwd`"
|
||||
|
||||
/usr/bin/modulecmd tcsh load mcstas
|
||||
bash compile_if_needed.sh
|
||||
|
||||
for y1 in $(seq -0.009 0.001 0.009)
|
||||
do
|
||||
for y2 in $(seq -0.009 0.001 0.009) $(seq -0.09 0.01 -0.01) $(seq 0.01 0.01 0.09) $(seq -0.5 0.1 -0.1) $(seq 0.1 0.1 0.5)
|
||||
do
|
||||
printf -v Y1 %.3f $y1
|
||||
printf -v Y2 %.3f $y2
|
||||
echo $Y1 $Y2
|
||||
/afs/psi.ch/project/sinq/sl6-64/bin/mpirun -np $SLURM_NPROCS Estia_baseline.out \
|
||||
--dir=../results/selene1_geo_$Y1\_$Y2 --ncount=1e9 \
|
||||
omegaa=1.0 sample=4 sample_length=0.00005 sample_height=0.01 \
|
||||
lambda_start=2.5 lambda_end=15.0 enable_gravity=1 enable_chopper=1 \
|
||||
lambda_min=3.75 selene1_foot1y=$Y1 selene1_foot2y=$Y2
|
||||
done
|
||||
done
|
||||
|
||||
|
||||
echo "Program finished with exit code $? at: `date`"
|
||||
|
||||
@@ -0,0 +1,40 @@
|
||||
#!/bin/bash
|
||||
#SBATCH -J mcEstia_2
|
||||
#SBATCH -N 2
|
||||
#SBATCH --ntasks-per-node=24
|
||||
#SBATCH --time=1-00:00:00
|
||||
#SBATCH --mail-type=fail
|
||||
#SBATCH --mail-user=artur.glavic@psi.ch
|
||||
|
||||
#SBATCH -o stdout2.log
|
||||
#SBATCH -e stderr2.log
|
||||
|
||||
#SBATCH --partition=ll_long
|
||||
|
||||
echo "Starting at `date`"
|
||||
echo "Running on hosts: $SLURM_NODELIST"
|
||||
echo "Running on $SLURM_NNODES nodes."
|
||||
echo "Running on $SLURM_NPROCS processors."
|
||||
echo "Current working directory is `pwd`"
|
||||
|
||||
/usr/bin/modulecmd tcsh load mcstas
|
||||
bash compile_if_needed.sh
|
||||
|
||||
for y1 in $(seq -0.09 0.01 -0.01)
|
||||
do
|
||||
for y2 in $(seq -0.009 0.001 0.009) $(seq -0.09 0.01 -0.01) $(seq 0.01 0.01 0.09) $(seq -0.5 0.1 -0.1) $(seq 0.1 0.1 0.5)
|
||||
do
|
||||
printf -v Y1 %.3f $y1
|
||||
printf -v Y2 %.3f $y2
|
||||
echo $Y1 $Y2
|
||||
/afs/psi.ch/project/sinq/sl6-64/bin/mpirun -np $SLURM_NPROCS Estia_baseline.out \
|
||||
--dir=../results/selene1_geo_$Y1\_$Y2 --ncount=1e9 \
|
||||
omegaa=1.0 sample=4 sample_length=0.00005 sample_height=0.01 \
|
||||
lambda_start=2.5 lambda_end=15.0 enable_gravity=1 enable_chopper=1 \
|
||||
lambda_min=3.75 selene1_foot1y=$Y1 selene1_foot2y=$Y2
|
||||
done
|
||||
done
|
||||
|
||||
|
||||
echo "Program finished with exit code $? at: `date`"
|
||||
|
||||
@@ -0,0 +1,40 @@
|
||||
#!/bin/bash
|
||||
#SBATCH -J mcEstia_3
|
||||
#SBATCH -N 2
|
||||
#SBATCH --ntasks-per-node=24
|
||||
#SBATCH --time=1-00:00:00
|
||||
#SBATCH --mail-type=fail
|
||||
#SBATCH --mail-user=artur.glavic@psi.ch
|
||||
|
||||
#SBATCH -o stdout3.log
|
||||
#SBATCH -e stderr3.log
|
||||
|
||||
#SBATCH --partition=ll_long
|
||||
|
||||
echo "Starting at `date`"
|
||||
echo "Running on hosts: $SLURM_NODELIST"
|
||||
echo "Running on $SLURM_NNODES nodes."
|
||||
echo "Running on $SLURM_NPROCS processors."
|
||||
echo "Current working directory is `pwd`"
|
||||
|
||||
/usr/bin/modulecmd tcsh load mcstas
|
||||
bash compile_if_needed.sh
|
||||
|
||||
for y1 in $(seq 0.01 0.01 0.09)
|
||||
do
|
||||
for y2 in $(seq -0.009 0.001 0.009) $(seq -0.09 0.01 -0.01) $(seq 0.01 0.01 0.09) $(seq -0.5 0.1 -0.1) $(seq 0.1 0.1 0.5)
|
||||
do
|
||||
printf -v Y1 %.3f $y1
|
||||
printf -v Y2 %.3f $y2
|
||||
echo $Y1 $Y2
|
||||
/afs/psi.ch/project/sinq/sl6-64/bin/mpirun -np $SLURM_NPROCS Estia_baseline.out \
|
||||
--dir=../results/selene1_geo_$Y1\_$Y2 --ncount=1e9 \
|
||||
omegaa=1.0 sample=4 sample_length=0.00005 sample_height=0.01 \
|
||||
lambda_start=2.5 lambda_end=15.0 enable_gravity=1 enable_chopper=1 \
|
||||
lambda_min=3.75 selene1_foot1y=$Y1 selene1_foot2y=$Y2
|
||||
done
|
||||
done
|
||||
|
||||
|
||||
echo "Program finished with exit code $? at: `date`"
|
||||
|
||||
@@ -0,0 +1,40 @@
|
||||
#!/bin/bash
|
||||
#SBATCH -J mcEstia_4
|
||||
#SBATCH -N 2
|
||||
#SBATCH --ntasks-per-node=24
|
||||
#SBATCH --time=1-00:00:00
|
||||
#SBATCH --mail-type=fail
|
||||
#SBATCH --mail-user=artur.glavic@psi.ch
|
||||
|
||||
#SBATCH -o stdout4.log
|
||||
#SBATCH -e stderr4.log
|
||||
|
||||
#SBATCH --partition=ll_long
|
||||
|
||||
echo "Starting at `date`"
|
||||
echo "Running on hosts: $SLURM_NODELIST"
|
||||
echo "Running on $SLURM_NNODES nodes."
|
||||
echo "Running on $SLURM_NPROCS processors."
|
||||
echo "Current working directory is `pwd`"
|
||||
|
||||
/usr/bin/modulecmd tcsh load mcstas
|
||||
bash compile_if_needed.sh
|
||||
|
||||
for y1 in $(seq -0.5 0.1 -0.1)
|
||||
do
|
||||
for y2 in $(seq -0.009 0.001 0.009) $(seq -0.09 0.01 -0.01) $(seq 0.01 0.01 0.09) $(seq -0.5 0.1 -0.1) $(seq 0.1 0.1 0.5)
|
||||
do
|
||||
printf -v Y1 %.3f $y1
|
||||
printf -v Y2 %.3f $y2
|
||||
echo $Y1 $Y2
|
||||
/afs/psi.ch/project/sinq/sl6-64/bin/mpirun -np $SLURM_NPROCS Estia_baseline.out \
|
||||
--dir=../results/selene1_geo_$Y1\_$Y2 --ncount=1e9 \
|
||||
omegaa=1.0 sample=4 sample_length=0.00005 sample_height=0.01 \
|
||||
lambda_start=2.5 lambda_end=15.0 enable_gravity=1 enable_chopper=1 \
|
||||
lambda_min=3.75 selene1_foot1y=$Y1 selene1_foot2y=$Y2
|
||||
done
|
||||
done
|
||||
|
||||
|
||||
echo "Program finished with exit code $? at: `date`"
|
||||
|
||||
@@ -0,0 +1,40 @@
|
||||
#!/bin/bash
|
||||
#SBATCH -J mcEstia_5
|
||||
#SBATCH -N 2
|
||||
#SBATCH --ntasks-per-node=24
|
||||
#SBATCH --time=1-00:00:00
|
||||
#SBATCH --mail-type=fail
|
||||
#SBATCH --mail-user=artur.glavic@psi.ch
|
||||
|
||||
#SBATCH -o stdout5.log
|
||||
#SBATCH -e stderr5.log
|
||||
|
||||
#SBATCH --partition=ll_long
|
||||
|
||||
echo "Starting at `date`"
|
||||
echo "Running on hosts: $SLURM_NODELIST"
|
||||
echo "Running on $SLURM_NNODES nodes."
|
||||
echo "Running on $SLURM_NPROCS processors."
|
||||
echo "Current working directory is `pwd`"
|
||||
|
||||
/usr/bin/modulecmd tcsh load mcstas
|
||||
bash compile_if_needed.sh
|
||||
|
||||
for y1 in $(seq 0.1 0.1 0.5)
|
||||
do
|
||||
for y2 in $(seq -0.009 0.001 0.009) $(seq -0.09 0.01 -0.01) $(seq 0.01 0.01 0.09) $(seq -0.5 0.1 -0.1) $(seq 0.1 0.1 0.5)
|
||||
do
|
||||
printf -v Y1 %.3f $y1
|
||||
printf -v Y2 %.3f $y2
|
||||
echo $Y1 $Y2
|
||||
/afs/psi.ch/project/sinq/sl6-64/bin/mpirun -np $SLURM_NPROCS Estia_baseline.out \
|
||||
--dir=../results/selene1_geo_$Y1\_$Y2 --ncount=1e9 \
|
||||
omegaa=1.0 sample=4 sample_length=0.00005 sample_height=0.01 \
|
||||
lambda_start=2.5 lambda_end=15.0 enable_gravity=1 enable_chopper=1 \
|
||||
lambda_min=3.75 selene1_foot1y=$Y1 selene1_foot2y=$Y2
|
||||
done
|
||||
done
|
||||
|
||||
|
||||
echo "Program finished with exit code $? at: `date`"
|
||||
|
||||
-551
File diff suppressed because one or more lines are too long
Reference in New Issue
Block a user