
System PSLC (vms 4.7)

Username: XRAYOP
Password: 
                    Welcome to VAX/VMS Version 4.7
                  University of Wisconsin - Madison
                     Physical Sciences Laboratory

    Last interactive login on Tuesday, 19-JUL-1988 22:33
    Last non-interactive login on Saturday, 16-JUL-1988 16:11

PSLA will be unavailable Wednesday, 20-July, from 6:30am to 8:00am.
  
       DECnet node tables were updated 18-Jul-1988.

  See BULL "important user messages" for important announcements concerning:
    + IMSL support on PSLB being replaced by CERN library support,
    + Announced charges for ethernet connections are taking effect 1-Sep.
    + New charges for bitnet use are taking effect 1-Sep,
    + HEPnet Decnet nodes have been moved from area 42 to 46

  In order for PSL to provide the maximum amount of uptime possible:
  PSLA   will be down 27-Jul 07:00-08:30 AM for preventive maintenance
  PSLOAS will be down 27-Jul 08:30-10:00 AM for preventive maintenance


    You are logged into _RTA2:





 
 
 
 
 
Defining Logical Names and Environmental Parameters
Graphic device set at TT:. Enter $SETPLOT to modify.





All SHADOW logical names and command defined.
Type $NEWS for more info on SHADOW.


$ SD [.UNDUL]
$ CREATE/DIR [.DEMO]
$ SD [.DEMO]
$ DIR
%DIRECT-W-NOFILES, no files found
$ EPATH
Enter number of periods : 
10
Enter wavelength of undulator (m) : 
.061
Enter magnetic field (tesla) : 
.07587
Enter electron energy (GeV) : 
0.8
Enter number of points along undulator axis: 
101
Enter angle along horizontal plane (mrads) : 
0
Enter angle with vertical (mrads) : 
0
Trajectory file for input to NPHOTON : 
NPHOTON
Parameters file-name : 
EPATH.PAR
Write out data for ERAD ? 
Y
Enter unformatted file name to be used by ERAD:
ERAD.DAT
Do you want a typable file (1 yes, 0 no)?
0
$ DIR

Directory D13:[XRAYOP.UNDUL.DEMO]

EPATH.PAR;1               2  19-JUL-1988 22:42  (RWED,RWED,RE,RE)
ERAD.DAT;1               16  19-JUL-1988 22:42  (RWED,RWED,RE,RE)
NPHOTON.DAT;1           290  19-JUL-1988 22:42  (RWED,RWED,RE,RE)

Total of 3 files, 308 blocks.
$ TY EPATH.PAR
Parameters used for run creating
ERAD.DAT                                                                        
------   INPUT    ---------
Number of periods               =           10
Undulator Wavelength      [ m ] =   6.1000000000000000E-02
Beam Energy             [ GeV ] =   0.8000000000000000    
Peak Magnetic field   [ Tesla ] =   7.5870000000000000E-02
Phi_E ( x )           [ mrads ] =   0.0000000000000000E+00
The_E ( y )           [ mrads ] =   0.0000000000000000E+00
------   OUTPUT    ---------
Gamma                           =    1565.375526427099    
Fundamental              [ Ev ] =    91.11653276369593    
Fundamental       [ Angstroms ] =    136.0891066029992    
K                               =   0.4320960741393696    
Equivalent SR C.W.     [ Angs ] =    478.8206801107157    
              C.E.       [ eV ] =    25.89386907251613    
$ PRE_RAD
Enter number of energy points   (odd, max=51) :
7
Enter number of points in theta (odd, max=31) :
31
Enter number of points in phi   (odd, max=31) :
31
Enter filename for input  (from EPATH   ):
ERAD.DAT
Enter filename for parameter info       ):
SRC.PAR
Trajectory computed by EPATH with following parameters:
Number of periods      :          10
Number of points       :         101
Wavelength of undulator:  6.1000000000000000E-02meters
Fundamental energy     :   91.11653276369593    eV
Fundamental wavelength :   136.0891066029992    Angstrom
Deflection parameter   :  0.4320960741393696    
Horz. traj. angle      :  0.0000000000000000E+00mrad
Vert. traj. angle      :  0.0000000000000000E+00mrad
Gamma		       :   1565.375526427099    
Initial beta along x   :  0.0000000000000000E+00C units
... 	     along y   :  0.0000000000000000E+00C units
...  	     along z   :  0.9999997959517054    C units
Field B0	       :  7.5870000000000000E-02tesla
Electron energy	       :  0.8000000000000000    GeV
How many periods to use ? 
1
Enter starting energy [ eV ] : 
91.0
 ...  final    energy [ eV ] : 
91.25
Electron current      [ A  ] : 
0.1
Enter [ 0 ] for constant bandpass
      [ 1 ]              energy interval
Then ? 
0
Bandpass [ e.g., for 1% enter 0.01 ] ? 
0.001
WARNING::
For CARTESIAN angles, A file for SHADOW will be generated
containing (binary CDFs and polarization).  You will have
the option to create an ASCII file containing the photon
flux at a chosen energy interval.
 
WARNING::
For POLAR angles, A file for SHADOW will not be created.
You will have the option to create an ASCII file containing
the photon flux (phi,the,ener),photon flux(ener,the) or
photon flux(ener) at chosen energy interval or constant
bandpass.
POLAR ANGLES = 1, OR CARTESIAN = 2  ? 
2
NOTE:: If you chose a pinhole or rectangular aperture
       the number of points in theta and phi means the
       number in 1 quadrant.  for pinhole, enter
       phimin = 0.0 phimax = pi/2 in mrads.
Enter 3 for centered rectangle or 4 for rectangle
APERTURE TYPE ?  
3
Enter initial angle phi (mrad):
0
Enter final   angle phi (mrad):
1
Enter initial angle theta (mrad):
0
Enter final   angle theta (mrad):
1
Begin calculations.
BEGIN INTEGRATION OVER  PHI
BEGIN INTEGRATION OVER THETA
Do you want to READ or WRIT RNs? 
WRIT
Enter filename for RNis 
RN.DAT
Do you want to READ or WRIT CDFs? 
WRIT
Enter file name to write CDFs and pol. to SHADOW 
CDF.DAT
DATA WRITTEN TO FILE      
CDF.DAT                                                                         
PARAMETER INFO WRITTEN TO 
SRC.PAR                                                                         
$ DIR

Directory D13:[XRAYOP.UNDUL.DEMO]

CDF.DAT;1               215  19-JUL-1988 22:47  (RWED,RWED,RE,RE)
EPATH.PAR;1               2  19-JUL-1988 22:42  (RWED,RWED,RE,RE)
ERAD.DAT;1               16  19-JUL-1988 22:42  (RWED,RWED,RE,RE)
FOR070.DAT;1            657  19-JUL-1988 22:45  (RWED,RWED,RE,RE)
NPHOTON.DAT;1           290  19-JUL-1988 22:42  (RWED,RWED,RE,RE)
RN.DAT;1                109  19-JUL-1988 22:47  (RWED,RWED,RE,RE)
SRC.PAR;1                 5  19-JUL-1988 22:47  (RWED,RWED,RE,RE)

Total of 7 files, 1294 blocks.
$ GO SOURCE



  :::::::  :::    ::            :::    :::::::        :::::     :::  ::  ::
 ::::::::  :::    ::           ::::    ::::::::      :::::::    :::  ::  ::
:::        :::    ::          ::: :    :::   :::    :::    ::   :::  ::  ::
:::        :::    ::         ::: ::    :::    ::    :::    ::   :::  ::  ::
 :::::::   :::::::::        :::  ::    :::    ::    :::    ::   :::  ::  ::
  :::::::  :::::::::       :::   ::    :::    ::    :::    ::   :::  ::  ::
       ::  :::    ::      :::    ::    :::    ::    :::    ::   :::  ::  ::
      :::  :::    ::     :::     ::    :::   :::    :::   :::   ::: ::: :::
::::::::   :::    ::    :::      ::    ::::::::      :::::::     :::::::::
:::::::    :::    ::   :::       ::    :::::::        :::::       ::: :::




VERSION     88.2                                               1 June  1988




Hit <RETURN> when ready: 


SOURCE selected. Begin procedure.


This procedure generates a SOURCE for SHADOW.
Mode selected [ ? <ret> for HELP ] ? : PROMPT
Call to INPUT




---------------------------------- S H A D O W  ------------------------------- 
Jun 1988                                         F.Cerrina        ECE - UW      

Defining source : 
When prompted for a yes/no answer, you may enter:
for   YES  answer      Y, 1
for   NO   answer      anything else




Do you want a verbose [ 1 ] or terse [ 0 ] output ?
1

------------- SOURCE SPECS ------------------
Options available:
Random in  BOTH REAL and MOMENTUM space		      0
Grid       BOTH					      1
Grid       REAL, random  MOMENTUM	              2
Random     REAL,  grid   MOMENTUM		      3
Ellipses in PHASE space, random around each ellipse   4
Ellipses in PHASE space,   grid around each ellipse   5
Source modelling type [ 0-5 ] ? 
0
How many rays [ 1 - 5 000 ] ? 
3000
Seed [ odd, 1000 - 1 000 000 ] ? 
12357
Do you want to optimize the source ? 
0

Source type : [ 0 ] regular source
              [ 1 ] wiggler
              [ 2 ] undulator
Then ? 
2
File containing the CDF's ? 
CDF.DAT
Sigma along X ? 
1.4
            Z ? 
0.14
Electron beam emittance.  Units are : rads*[ units of length used so far ]
Beam emittances in X [ at waist ] ? 
0
Distance of insertion device's center from X waist [ signed ] ? 
0
Beam emittances in Z [ at waist ] ? 
0
Distance of insertion device's center from Z waist [ signed ] ? 
0

Polarization component of interest. Enter 
parallel polarization    1
perpendicular            2
total                    3
                   then ? 
Polarization Selected [ 1-3 ] ? 
3

Do you want to store the optical paths (OPD) [Y/N] ? 
0
Do you want to generate the A vectors (electric field) [Y/N] ?
0
Exit from INPUT_SOURCE


Generated          250 rays out of         3000
                   500
                   750
                  1000
                  1250
                  1500
                  1750
                  2000
                  2250
                  2500
                  2750
                  3000
Exit from SOURCE
SOURCE => Source has been succesfully generated.

SOURCE procedure completed.

Shadow:: EXIT



Exit to DCL

$ PLOTXY


PLOT> Input file? 
BEGIN.DAT
PLOT> Options --- Enter
PLOT>   0   for excluding the losses
PLOT>   1   for including only the losses
PLOT>   2   for including all the rays.
PLOT> Then ? 

PLOT> Comment for plot [ 80 char ] ?
SRC undulator radiation pattern

PLOT> File read OK. Full specifications:
     D13:[XRAYOP.UNDUL.DEMO]BEGIN.DAT;1                          
     Was created : 19-JUL-1988 22:49
 
PLOT>  Found         3000 good points out of        3000
PLOT>  The following columns are defined for each ray :
       1)  the regular columns [1-12]



Col Par Minimum:       Maximum:       Center:        St. Dev.:

 1  X    -5.4948         5.1509       -0.38320E-01     1.4023    
 2  Y    0.00000E+00    0.00000E+00    0.00000E+00    0.00000E+00
 3  Z   -0.48262        0.50482       -0.58630E-02    0.13997    
 4  X'  -0.99946E-03    0.99816E-03    0.17575E-04    0.51717E-03
 5  Y'    1.0000         1.0000         1.0000        0.24414E-03
 6  Z'  -0.99929E-03    0.99995E-03    0.12254E-04    0.54643E-03
11   Photon Energy  (eV)      91.000         91.250    
20   Numerical Aperture      0.15206E-04    0.14043E-02
 
PLOT> Options. You may plot any two rows from the above
      list versus each other. You may also plot any of
      them versus the ray Numerical Aperture.
      N.A. -- enter 20.
 
PLOT> Rows to use for plot :
PLOT>   for horizontal axis ? 
4
PLOT>   for vertical axis   ? 
6
 
PLOT> Scaling options. Enter 
 
       0   For automatic scaling
       1   For cartesian scaling
       2   For external limits
PLOT> Then ? 
2
PLOT> Enter limits. 
PLOT> Horizontal min.: 
-.001
PLOT> Horizontal max.: 
.001
PLOT> Vertical min.  : 
-.001
PLOT> Vertical max.  : 
.001
PLOT> Plotting options : 
       0   For scattered plot
       1   For connected plot
       2   For contour plot
PLOT> Then ? 
0
PLOT> Hairline at [ 0,0 ] ? 
0
PLOT> Overlay a mirror/slit ? 
0
PLOT> Ready for histograms. Enter:
     -1         to skip 
      0     for same limits as plot
      1                     3*stdev
      2                     external
PLOT> ? 
-1
PLOT> Do you want to save the TD command file ? 
1
PLOT> Terminal type: [ 0 ] VT240, [ 1 ] TEK 40xx, [ 2 ] HIREZ, [ 3 ] TEK 4107 : 
PLOT> Terminal type: [ 0 ] VT240, [ 1 ] TEK 40xx, [ 2 ] HIREZ, [ 3 ] TEK 4107 : 
