Commit df0afbf6 authored by Thomas Planche's avatar Thomas Planche

plot two masses using REBELOTE

parent 17213688
zgoubi HRS tuto: DIPOLE
'OBJET'
543.701 !543.687 or 543.701 and not 543.714, since we are aiming for a mass resolution of 20000 and not momentum resolution of 20000
543.687 !543.687 or 543.701 and not 543.714, since we are aiming for a mass resolution of 20000 and not momentum resolution of 20000
3.01 !read from particles from
1 3000 1 !IP1; IP2; IPStep (Read particles numbered IT1 to IT2, step ITStep)
1 1 1 !not used here
......@@ -61,7 +60,12 @@ noEspread.beam ! name of the input file
36.32357
'FAISCNL'
mass2.fai
twomasses.fai
'REBELOTE'
1 0.1 0 1
1
1 1 543.701
'END'
......
%!PS-Adobe-2.0 EPSF-2.0
%%Title: twomasses.eps
%%Creator: gnuplot 5.3 patchlevel 0
%%CreationDate: Sat Jul 6 23:47:21 2019
%%CreationDate: Sun Jul 7 12:58:38 2019
%%DocumentFonts: (atend)
%%BoundingBox: 50 50 410 302
%%EndComments
......@@ -481,7 +481,7 @@ SDict begin [
/Creator (gnuplot 5.3 patchlevel 0)
% /Producer (gnuplot)
% /Keywords ()
/CreationDate (Sat Jul 6 23:47:21 2019)
/CreationDate (Sun Jul 7 12:58:38 2019)
/DOCINFO pdfmark
end
} ifelse
......@@ -710,6 +710,8 @@ LCb setrgbcolor
1.000 UL
LTb
0.58 0.00 0.83 C
5277 1434 CircleF
5124 3774 CircleF
5506 904 CircleF
5251 1845 CircleF
5459 1184 CircleF
......@@ -3731,6 +3733,8 @@ LTb
1.000 UL
LTb
0.00 0.62 0.45 C
5672 1434 CircleF
5478 3774 CircleF
5913 904 CircleF
5638 1845 CircleF
5860 1184 CircleF
......
......@@ -10,12 +10,12 @@ cm2mm=10.0
set key at -0.5, 5
set xrange [-1.5:0.5]
set yrange [-60:60]
plot "mass1.fai" u ($10*cm2mm):11 w points pt 7 ps 0.2 lc 1 notitle, \
plot "twomasses.fai" u ($10*cm2mm):($37==543.687 ? $11:1/0) w points pt 7 ps 0.2 lc 1 notitle, \
"" w points pt 7 ps 0.7 lc 1 title "^{238}U mass ", \
"mass2.fai" u ($10*cm2mm):11 w points pt 7 ps 0.2 lc 2 notitle ,\
"" u ($10*cm2mm):($37==543.701 ? $11:1/0) w points pt 7 ps 0.2 lc 2 notitle ,\
"" w points pt 7 ps 0.7 lc 2 title "^{238}U mass*20001/20000"
zgoubi HRS tuto: DIPOLE
'OBJET'
543.687 !543.687 or 543.701 and not 543.714, since we are aiming for a mass resolution of 20000 and not momentum resolution of 20000
3.01 !read from particles from
......@@ -61,7 +60,12 @@ noEspread.beam ! name of the input file
36.32357 ! drift length=80-120*tan(20/deg.)
'FAISCNL'
mass1.fai
twomasses.fai
'REBELOTE'
1 0.1 0 1
1
1 1 543.701
'END'
......
%!PS-Adobe-2.0 EPSF-2.0
%%Title: twomasses.eps
%%Creator: gnuplot 5.3 patchlevel 0
%%CreationDate: Sat Jul 6 23:47:31 2019
%%CreationDate: Sun Jul 7 12:56:36 2019
%%DocumentFonts: (atend)
%%BoundingBox: 50 50 410 302
%%EndComments
......@@ -481,7 +481,7 @@ SDict begin [
/Creator (gnuplot 5.3 patchlevel 0)
% /Producer (gnuplot)
% /Keywords ()
/CreationDate (Sat Jul 6 23:47:31 2019)
/CreationDate (Sun Jul 7 12:56:36 2019)
/DOCINFO pdfmark
end
} ifelse
......
......@@ -15,8 +15,7 @@ set xrange [-1.5:0.5]
set yrange [-60:60]
plot "mass1.fai" u ($10*cm2mm):11 w points pt 7 ps 0.2 lc 1 notitle, \
plot "twomasses.fai" u ($10*cm2mm):($37==543.687 ? $11:1/0) w points pt 7 ps 0.2 lc 1 notitle, \
"" w points pt 7 ps 0.7 lc 1 title "^{238}U mass ", \
"mass2.fai" u ($10*cm2mm):11 w points pt 7 ps 0.2 lc 2 notitle ,\
"" u ($10*cm2mm):($37==543.701 ? $11:1/0) w points pt 7 ps 0.2 lc 2 notitle ,\
"" w points pt 7 ps 0.7 lc 2 title "^{238}U mass*20001/20000"
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment