'*PLOT*MICE Muon Beamline Model for TURTLE'

(Reference version 23 May 2007. Includes TOF0/Ckov/TOF1 scattering patches and air)

(FILE3 - Beam from just before Dipole 2 to end of Pb Scatterer)
(CM14-7.1-200 version)

(want to try to maximise the flux of the 10 dp cut particles)

(Uses Quad Aps 23.6 outer circle)

(TTL only:- 134589 from ref. File2.)
134589

16.0 190. 0 100 ; (dump hists to for100.dat)
13 10 ; (TTL apertures resp by parents.)
-16.0 20.0 0.0; (TTL, +16 20 is for homogeneous distn.)
16.0 30.0 1.0 ; (TTL only)

16. 162.0 ; (Slits made rectangular)

1. -32.0 0.0 0.0 0.0 0.0 0.0 0.258094 ; (pick up data from for032.dat)

17.0 3.0; (1st Order 2-2nd 3-3rd.)
16.0 3.0 206.75 /Mass/ ; (mass term for muon multiple-cattering)

16.0 11.0 0.258094 ; (central momentum 0.258094)


-51.0 1 -40 30 1.0 'b4B2_X';
-52.0 3 -30 30 1.0 'b4B2_Y';
(---- Bend Magnet ----)
(Muon Spectrometer. Type I Dipole)
16.0 4.0 55.88 ;  (Dipole hz half ap - now 111.76cm/2 )
16.0 5.0 10.0 ; (Vert half ap.)
-16.0 7.0 0.0 ; (No fringe fields)
16.0 7.0 0.65348;
16.0 8.0 2.8;
20.0 180.0 ;
2.0 15.0 ;
4.0 1.04995 4.293258339 0.0 ; (4.245726 for 255.240, 4.293258339 for 258.094)
2.0 15.0 ;
20.0 -180.0 ;
-51.0 1 -40 40 1.0 'afB2_X';
-52.0 3 -30 30 1.0 'afB2_Y';

(---- Start of muon extraction straight ----)

(BlockIII)
16.0 165.0 ; (Enable m/s)
16.0 160.0 ; (Slits are inverted ie. targets.)
6.0 1.0 25.0 3.0 25.0 ; (A 50cm x 50cm area)
(For 17.0cm 255.240 Mev/c ~ Rhox4 or so )
(Note the 400.0 ASN is a guessed number, but maynot make much diff here..)
1.0 3.208 9.67 37750.0 1.000 400.0 0.0000456 0.0002135 0.0 ;
3.0 0.17 ; (length)
16.0 166.0 ; (Disable m/s)
16.0 161.0 ; (Slits are taken normal again.)
16.0 11.0 0.257832 ; (For cm14 B2 p0) 

3.0 0.435065 ;
3.0 0.032065 'D71';      (* MAR05: 0.435065. For cm14 layout: 0.435065-0.403 = 0.032065)

(BlockII.1)
16.0 165.0 ; (Enable m/s)
16.0 160.0 ; (Slits are inverted ie. targets.)
6.0 1.0 25.0 3.0 25.0 ; (A 50cm x 50cm area)
(For 17.0cm 254.974 Mev/c ~ Rhox4 or so )
(AIR. model at 254.974: from MuSCAT siga=7.0, sige=0.01 )
(Note the 400.0 ASN is a guessed number, but maynot make much diff here..)
1.0 3.509 9.84 32290.0 1.000 400.0 0.0000534 0.0002496 0.0 ;
3.0 0.17 ; (length)
16.0 166.0 ; (Disable m/s)
16.0 161.0 ; (Slits are taken normal again.)
16.0 11.0 0.257521 ;(For cm14 B2 p0)


3.0 0.193; (MP)

(---- Second Triplet Quad Set ----)

16.0  100.0 23.6 23.6;
16.0  101.0 17.82 0.0;
-51.0 1 -70 60 2.0 'b4Q4_3_X';
-52.0 3 -30 30 1.0 'b4Q4_3_Y';
5.0 0.66 1.76356892 17.82 'q4';
-5.0 0.33 1.76356892 17.82 'q4'; ( 3.326489 for 254.663 )
-51.0 1 -30 30 1.0 'inQ4_3c_X';
-52.0 3 -30 30 1.0 'inQ4_3c_Y';
-5.0 0.33 1.76356892 17.82 'q4'; 
-51.0 1 -30 30 1.0 'afQ4_3c_X';
-52.0 3 -30 30 1.0 'afQ4_3c_Y';

3.0 0.193; (MP)

(BlockII.2)
16.0 165.0 ; (Enable m/s)
16.0 160.0 ; (Slits are inverted ie. targets.)
6.0 1.0 25.0 3.0 25.0 ; (A 50cm x 50cm area)
(Note the 400.0 ASN is a guessed number, but maynot make much diff here..)
1.0 3.509 9.84 32290.0 1.000 400.0 0.0000534 0.0002496 0.0 ;
3.0 0.17 ; (length)
16.0 166.0 ; (Disable m/s)
16.0 161.0 ; (Slits are taken normal again.)
16.0 11.0 0.25721 ; (For cm14 B2 p0)

3.0 0.137; (MP)

16.0  100.0 23.6 23.6;
16.0  101.0 17.82 0.0;
-51.0 1 -30 30 1.0 'b4Q5_2_X';
-52.0 3 -35 35 1.0 'b4Q5_2_Y';
5.0 0.66 -2.36224363 17.82 'q5';
-5.0 0.33 -2.36224363 17.82 'q5'; ( -3.068433 for 254.352 )
-51.0 1 -30 30 1.0 'inQ5_2c_X';
-52.0 3 -30 30 1.0 'inQ5_2c_Y';
-5.0 0.33 -2.36224363 17.82 'q5';
-51.0 1 -30 30 1.0 'afQ5_2c_X';
-52.0 3 -30 30 1.0 'afQ5_2c_Y';

3.0 0.193; (MP)

(BlockII.3)
16.0 165.0 ; (Enable m/s)
16.0 160.0 ; (Slits are inverted ie. targets.)
6.0 1.0 25.0 3.0 25.0 ; (A 50cm x 50cm area)
(Note the 400.0 ASN is a guessed number, but maynot make much diff here..)
1.0 3.509 9.84 32290.0 1.000 400.0 0.0000534 0.0002496 0.0 ;
3.0 0.17 ; (length)
16.0 166.0 ; (Disable m/s)
16.0 161.0 ; (Slits are taken normal again.)
16.0 11.0 0.256899 ; (For cm14 B2 p0) 


3.0 0.137; (MP)

16.0  100.0  23.6 23.6;
16.0  101.0  17.82 0.0;
-51.0 1 -35 30 1.0 'b4Q6_1_X';
-52.0 3 -25 25 1.0 'b4Q6_1_Y';
5.0 0.66 1.567434686 17.82 'q6';
-5.0 0.33 1.567434686 17.82 'q6'; ( 3.157349 for 254.041 )
-51.0 1 -25 25 1.0 'inQ6_1c_X';
-52.0 3 -25 25 1.0 'inQ6_1c_Y';
-5.0 0.33 1.567434686 17.82 'q6'; 
-51.0 1 -25 25 1.0 'afQ6_1c_X';
-52.0 3 -25 25 1.0 'afQ6_1c_Y';

3.0 0.023 ; (Part of MP)

(Next part of MP is AirBlockI:-)
16.0 165.0 ; (Enable m/s)
16.0 160.0 ; (Slits are inverted ie. targets.)
6.0 1.0 25.0 3.0 25.0 ; (A 50cm x 50cm area)
(For 17.0cm 254.041 Mev/c ~ Rhox4 or so )
(AIR. model at 254.041: from MuSCAT siga=7.0, sige=0.01 )
(Note the 400.0 ASN is a guessed number, but maynot make much diff here..)
1.0 2.323 9.02 67180.0 1.000 400.0 0.0000257 0.0001200 0.0 ;
3.0 0.17 ; (length)
16.0 166.0 ; (Disable m/s)
16.0 161.0 ; (Slits are taken normal again.)
16.0 11.0 0.25675 ; (For cm14 B2 p0) 
(End of MP)
(3.0 0.193;) (MP)

3.0 0.05965 ; (0.11045-0.0508)
3.0 0.247 ; 

-51.0 1 -30 30 1.0 'b4T0_X';
-52.0 3 -25 25 1.0 'b4T0_Y';
(---- TOF0 here: ----)
(NB Model is for MAR05 setting.)
16.0 165.0 ; (Enable m/s)
16.0 160.0 ; (Slits are inverted ie. targets.)
9.0 10.0; (Repeat section fllwing x5)
6.0 1.0 24.0 3.0 24.0 ; (A 48cm x 48cm area)
(For 5.08mm 253.892Mev/c. This is for the 0.89817 g cm-3 density.)
(C9H10. model at 253.892: from MuSCAT siga=7.0, sige=0.01. 800 ASN is a guess note)
1.0 5.996 11.12 614.197 0.999 800.0 0.0001897 0.0009170 0.0 ;
3.0 0.00508 ; (length)
9.0 0.0;
16.0 166.0 ; (Disable m/s)
16.0 161.0 ; (Slits are taken normal again.)
16.0 11.0 0.24511 ; (For cm14 B2 p0)

(TOF0 Random Scattering Addition to patch up to g4bl m/s angle)
(ttl ~ 16.96, g4bl ~ 21.79. hence agl 13.68, & for card make 2x13.68 = 27.36)
1. 0.0 27.36 0.0 27.36 0.0 0.0 0.0 0.0 ; (the x' & y' s are max amplitudes / scatterng)

-51.0 1 -30 30 1.0 'afT0_X';
-52.0 3 -25 25 1.0 'afT0_Y';

(AIR_TOF0-Ckov1)
3.0 0.0 ;   (* MAR05: 0.62965. For cm14 layout: 0.62965- #0.62965 part of 1.569# . So 0.0 )
16.0 165.0 ; (Enable m/s)
16.0 160.0 ; (Slits are inverted ie. targets.)
(Note 1x air block only in code means exact energy following is a little out....)
9.0 1.0 ;    (* MAR05: 9.0 3.0. For cm14 layout: need another 1.569-2x0.6295 removed. So rmv & add bit )
6.0 1.0 25.0 3.0 25.0 ; (A 50cm x 50cm area)
(For 17.0cm 242.252 Mev/c ~ Rhox7 or so )
(AIR. model at 242.252: from MuSCAT siga=7.0, sige=0.01 )
(Note the 400.0 ASN is a guessed number, but maynot make much diff here..)
1.0 2.406 9.00 69670.0 1.000 400.0 0.0000251 0.0001167 0.0 ;
3.0 0.17 ; (length)
9.0 0.0 ;
16.0 166.0 ; (Disable m/s)
16.0 161.0 ; (Slits are taken normal again.)
16.0 11.0 0.24467 ;  (For cm14 B2 p0) 
3.0 0.0 ;  (* MAR05: 0.62965. For cm14 layout: 0.62965- #0.62965 part of 1.569# . So 0.0 )

3.0 0.0303 ; (Since removed 0.62965 + 0.62965 + 0.17x2 = too much. Add 0.0303 to have rmvd 1.569)

(---- Ckov1 Module here: consists of 4parts. ----)
(This Ck module has been brought back about 2metres or so from mar05)
-51.0 1 -50 50 1.0 'b4CKV1_X';
-52.0 3 -50 50 1.0 'b4CKV1_Y';

(CkovAlFront)
16.0 165.0 ; (Enable m/s)
16.0 160.0 ; (Slits are inverted ie. targets.)
6.0 1.0 25.0 3.0 25.0 ; (A 50cm x 50cm area)
(For 1.0mm 241.812Mev/c)
(Al. model at 207.9: from MuSCAT siga=0.01, sige=0.01 )
1.0 6.587 10.57 191.001 1.000 331.6 0.0000878 0.0003875 0.0 ;
3.0 0.001 ; (length)
16.0 166.0 ; (Disable m/s)
16.0 161.0 ; (Slits are taken normal again.)
16.0 11.0 0.24418 ; (For cm14 B2 p0)

(Ckov_Radiator)
16.0 165.0 ; (Enable m/s)
16.0 160.0 ; (Slits are inverted ie. targets.)
9.0 4.0; (Repeat section fllwing x4)
6.0 1.0 25.0 3.0 25.0 ; (A 50cm x 50cm area)
(For 5.0mm 241.322Mev/c)
(CKOVI: Quartz window part: from MuSCAT siga=0.01, sige=0.01. 400 ASN is a guess note)
1.0 10.038 11.75 555.459 1.000 400.0 0.0002814 0.0012814 0.0 ;
3.0 0.005 ; (length)
9.0 0.0;
16.0 166.0 ; (Disable m/s)
16.0 161.0 ; (Slits are taken normal again.)
16.0 11.0 0.23765 ; (For cm14 B2 p0)

(Ckov_QuartzWindow here)
16.0 165.0 ; (Enable m/s)
16.0 160.0 ; (Slits are inverted ie. targets.)
6.0 1.0 25.0 3.0 25.0 ; (A 50cm x 50cm area)
(For 2.0mm 234.792Mev/c)
(CKOVI: Quartz window part: from MuSCAT siga=0.01, sige=0.01. 400 ASN is a guess note)
1.0 8.327 11.13 316.183 1.000 400.0 0.0001528 0.0006592 0.0 ;
3.0 0.002 ; (length)
16.0 166.0 ; (Disable m/s)
16.0 161.0 ; (Slits are taken normal again.)
16.0 11.0 0.23681 ; (For cm14 B2 p0)

(CkovLight Air Bit)
3.0 0.154 ;

(AIR CkovLight here)
16.0 165.0 ; (Enable m/s)
16.0 160.0 ; (Slits are inverted ie. targets.)
6.0 1.0 25.0 3.0 25.0 ; (A 50cm x 50cm area)
(For 17.0cm 233.952 Mev/c ~ Rhox4 or so )
(AIR. model at 233.952: from MuSCAT siga=7.0, sige=0.01 )
(Note the 400.0 ASN is a guessed number, but maynot make much diff here..)
1.0 2.217 8.77 86310.0 1.000 400.0 0.0000205 0.0000948 0.0 ;
3.0 0.17 ; (length)
16.0 166.0 ; (Disable m/s)
16.0 161.0 ; (Slits are taken normal again.)
16.0 11.0 0.23669 ; (For cm14 B2 p0)

(Ckv1 Random Scattering Addition to patch up to g4bl m/s angle)
(ttl ~ 20.02, g4bl ~ 23.6. hence agl 12.49, & for card make 2x12.49 = 24.99)
1. 0.0 24.99 0.0 24.99 0.0 0.0 0.0 0.0 ; (the x' & y' s are max amplitudes / scatterng)

3.0 0.154 'D91'; (0.34285-0.07444=0.26841)

-51.0 1 -80 80 2.0 'afCKV1_X';
-52.0 3 -85 85 2.0 'afCKV1_Y';
(End of Ckov1 Module)

(---- Nominal Mid Point---)

(Detailed Diagnostics on Muon beam here)
-51.3 1 -25 25 1.0 'XXD'; (First column for rij mtrx)
-52.3 2 -100 100 5 'XXD';
-51.3 1 -25 25 1.0 'XY';
-52.3 3 -25 25 1.0 'XY';
-51.3 1 -25 25 1.0 'XYD';
-52.3 4 -100 100 5 'XYD';
-51.3 1 -25 25 1.0 'XDL';
-52.3 6 -50 200 5 'XDL';

-51.3 2 -100 100 5 'XDY';(Second column for rij mtrx)
-52.3 3 -25 25 1.0 'XDY';
-51.3 2 -100 100 5 'XDYD';
-52.3 4 -100 100 5 'XDYD';
-51.3 2 -100 100 5 'XDDL';
-52.3 6 -50 200 5 'XDDL';

-51.3 3 -25 25 1.0 'YDY';(Third column for rij mtrx)
-52.3 4 -100 100 5 'YDY';
-51.3 3 -25 25 1.0 'YDEL';
-52.3 6 -50 200 5 'YDEL'; 

-51.3 4 -100 100 5 'YDDL';(Fourth column for rij mtrx)
-52.3 6 -50 200 5 'YDDL';

-50.3 6 -50 200 5.0 ;
-50.3 11 0.1 0.42 0.005 'DSPM';


(3.0 1.30856 'D114';) (0.80533 from itn1)
3.0 0.26 'D114';
(AIR Ckov_TOF1.Block2 here)
16.0 165.0 ; (Enable m/s)
16.0 160.0 ; (Slits are inverted ie. targets.)
6.0 1.0 25.0 3.0 25.0 ; (A 50cm x 50cm area)
(For 17.0cm 233.832 Mev/c ~ Rhox4 or so )
(AIR. model at 233.832: from MuSCAT siga=7.0, sige=0.01 )
(Note the 400.0 ASN is a guessed number, but maynot make much diff here..)
1.0 3.065 9.40 49090.0 1.000 400.0 0.0000361 0.0001667 0.0 ;
3.0 0.17 ; (length)
16.0 166.0 ; (Disable m/s)
16.0 161.0 ; (Slits are taken normal again.)
16.0 11.0 0.236479 ; (For cm14 B2 p0)
3.0 0.175 ;


3.0 0.085 ;
3.0 0.523 ;
-51.0 1 -85 85 2.0 'XY1_X';
-52.0 3 -85 85 2.0 'XY1_Y';

(AIR Ckov_TOF1.Block3.1 here - ctr is 0.608 fr ctr/Q7.)
16.0 165.0 ; (Enable m/s)
16.0 160.0 ; (Slits are inverted ie. targets.)
6.0 1.0 25.0 3.0 25.0 ; (A 50cm x 50cm area)
(AIR. model for 233.621incoming.. from MuSCAT siga=7.0, sige=0.01 )
(Note the 400.0 ASN is a guessed number, but maynot make much diff here..)
1.0 4.231 10.04 27840.0 1.000 400.0 0.0000636 0.0002939 0.0 ;
3.0 0.17 ; (length)
16.0 166.0 ; (Disable m/s)
16.0 161.0 ; (Slits are taken normal again.)
16.0 11.0 0.236107 ; (For cm14 B2 p0)

(In reality could have done with the extra 2x0.17 air in here etc)
3.0 0.830 ;
-51.0 1 -85 85 2.0 'XY2_X';
-52.0 3 -85 85 2.0 'XY2_Y';
3.0 0.994 ;
3.0 0.193; (MP)

(---- Third Triplet Quad Set ----)

16.0  100.0 23.6 23.6;
16.0  101.0 17.82 0.0;
-51.0 1 -90 90 2.0 'b4Q7_3_X';
-52.0 3 -90 90 2.0 'b4Q7_3_Y';
5.0 0.66 1.581672331 17.82 'q7';
-5.0 0.33 1.581672331 17.82 'q7'; ( 2.477666 for 233.249 )
-51.0 1 -30 30 1.0 'inQ7_3c_X';
-52.0 3 -30 30 1.0 'inQ7_3c_Y';
-5.0 0.33 1.581672331 17.82 'q7';
-51.0 1 -30 30 1.0 'afQ7_3c_X';
-52.0 3 -30 30 1.0 'afQ7_3c_Y';

3.0 0.193; (MP) 

(AIR Ckov_TOF1.Block3.2 here:-)
16.0 165.0 ; (Enable m/s)
16.0 160.0 ; (Slits are inverted ie. targets.)
6.0 1.0 25.0 3.0 25.0 ; (A 50cm x 50cm area)
(AIR. model.. from MuSCAT siga=7.0, sige=0.01 )
(Note the 400.0 ASN is a guessed number, but maynot make much diff here..)
1.0 4.231 10.04 27840.0 1.000 400.0 0.0000636 0.0002939 0.0 ;
3.0 0.17 ; (length)
16.0 166.0 ; (Disable m/s)
16.0 161.0 ; (Slits are taken normal again.)
16.0 11.0 0.235735 ; (For cm14 B2 p0)

3.0 0.137; (MP)   (* MAR05: 0.193. For cm14 layout: 0.193-0.056 = 0.137. Eats into MP rep!)

16.0  100.0 23.6 23.6;
16.0  101.0 17.82 0.0;
-51.0 1 -40 40 1.0 'b4Q8_4_X';
-52.0 3 -30 30 1.0 'b4Q8_4_Y';
5.0 0.66 -2.39278459 17.82 'q8';
-5.0 0.33 -2.39278459 17.82 'q8'; ( -2.63130 for 232.877 )
-51.0 1 -30 30 1.0 'inQ8_4c_X';
-52.0 3 -30 30 1.0 'inQ8_4c_Y';
-5.0 0.33 -2.39278459 17.82 'q8';
-51.0 1 -30 30 1.0 'afQ8_4c_X';
-52.0 3 -30 30 1.0 'afQ8_4c_Y';

3.0   0.193; (MP)

(AIR Ckov_TOF1.Block3.3 here:-)
16.0 165.0 ; (Enable m/s)
16.0 160.0 ; (Slits are inverted ie. targets.)
6.0 1.0 25.0 3.0 25.0 ; (A 50cm x 50cm area)
(AIR. model.. from MuSCAT siga=7.0, sige=0.01 )
(Note the 400.0 ASN is a guessed number, but maynot make much diff here..)
1.0 4.231 10.04 27840.0 1.000 400.0 0.0000636 0.0002939 0.0 ;
3.0 0.17 ; (length)
16.0 166.0 ; (Disable m/s)
16.0 161.0 ; (Slits are taken normal again.)
16.0 11.0 0.235363 ; (For cm14 B2 p0)

3.0 0.137; (MP)   (* MAR05: 0.193. For cm14 layout: 0.193-0.056 = 0.137. Eats into MP rep!)

16.0  100.0 23.6 23.6;
16.0  101.0 17.82 0.0;
-51.0 1 -40 40 1.0 'b4Q9_5_X';
-52.0 3 -30 30 1.0 'b4Q9_5_Y';
5.0 0.66 2.044825918 17.82 'q9';
-5.0 0.33 2.044825918 17.82 'q9'; ( 3.02917 for 232.505 )
-51.0 1 -30 30 1.0 'inQ9_5c_X';
-52.0 3 -30 30 1.0 'inQ9_5c_Y';
-5.0 0.33 2.044825918 17.82 'q9';
-51.0 1 -30 30 1.0 'afQ9_5c_X';
-52.0 3 -30 30 1.0 'afQ9_5c_Y';

3.0  0.193; (MP)

(AIR Ckov_TOF1.Block4 here:-)
16.0 165.0 ; (Enable m/s)
16.0 160.0 ; (Slits are inverted ie. targets.)
6.0 1.0 25.0 3.0 25.0 ; (A 50cm x 50cm area)
(For 17.0cm 232.505Mev/c ~ Rhox6.5 or so )
(AIR. model at 232.505: from MuSCAT siga=7.0, sige=0.01 )
(Note the 400.0 ASN is a guessed number, but maynot make much diff here..)
1.0 4.064 9.95 30200.0 1.000 400.0 0.0000587 0.0002713 0.0 ;
3.0 0.17 ; (length)
16.0 166.0 ; (Disable m/s)
16.0 161.0 ; (Slits are taken normal again.)

16.0 11.0 0.23502 ; (For cm14 B2 p0)


3.0 0.4291 ; (Should be at 0.54835 - but interfers with MP thus mvd d/str)

3.0 0.108 ;

-51.0 1 -25 25 1.0 'b4T1_X';
-52.0 3 -25 25 1.0 'b4T1_Y';

(---- TOF1 here: ----)
(NB Model is for MAR05 settings)
16.0 165.0 ; (Enable m/s)
16.0 160.0 ; (Slits are inverted ie. targets.)
9.0 10.0; (Repeat section fllwing x10)
6.0 1.0 24.0 3.0 24.0 ; (A 48cm x 48cm area)
(For 5.08mm 232.162Mev/c. This is for the 0.89817 g cm-3 density.)
(C9H10. model at 232.162: from MuSCAT siga=7.0, sige=0.01. 800 ASN is a guess note)
1.0 6.691 11.16 609.952 0.999 800.0 0.0001965 0.0009387 0.0 ;
3.0 0.00508 ; (length)
9.0 0.0;
16.0 166.0 ; (Disable m/s)
16.0 161.0 ; (Slits are taken normal again.)

16.0 11.0 0.22287 ; (For cm14 B2 p0)

(TOF1 Random Scattering Addition to patch up to g4bl m/s angle)
(Used LBNL CM11 presentation scatters from Feb05)
(ttl ~ 20.66, g4bl ~ 24.76. hence agl 13.64, & for card make 2x13.64 = 27.29)
(TOF0 at present)
1. 0.0 27.29 0.0 27.29 0.0 0.0 0.0 0.0 ; (the x' & y' s are max amplitudes / scatterng)

-51.0 1 -25 25 1.0 'afT1_X';
-52.0 3 -25 25 1.0 'afT1_Y';

3.0 0.177 ; (0.2024-0.0254)
3.0 0.0215 ; (0.1065-0.17/2)

(AIRTOF1_diffuser here:-)
16.0 165.0 ; (Enable m/s)
16.0 160.0 ; (Slits are inverted ie. targets.)
6.0 1.0 25.0 3.0 25.0 ; (A 50cm x 50cm area)
(For 17.0cm 220.012Mev/c Rhox3. )
(AIR. model at 220.012: from MuSCAT siga=7.0, sige=0.01 )
(Note the 400.0 ASN is a guessed number, but maynot make much diff here..)
1.0 2.678 9.02 70700.0 1.000 400.0 0.0000256 0.0001173 0.0 ;
3.0 0.17 ; (length)
16.0 166.0 ; (Disable m/s)
16.0 161.0 ; (Slits are taken normal again.)

16.0 11.0 0.22272 ; (For cm14 B2 p0)

3.0 0.1985 ; 
3.0 0.00377; (Movement of Diffuser 3.77mm downstream)

(53.3 ;) (Marker flag for muons reaching here.)

-51.0 1 -25 25 1.0 'b4Diff_X';
-52.0 3 -25 25 1.0 'b4Diff_Y';

(---- Diffuser here: ----)
(Scatterer. 9.5mm )
(Designed for 0.220 Gev/c)
16. 163.0 ; (Slits made round again)
16.0 165.0 ; (Enable m/s)
16.0 160.0 ; (Slits are inverted ie. targets.)
9.0 5.0; (Repeat section fllwing x5 )
6.0 1.0 15.0 3.0 15.0 ; (A 150mm radius)
(For 1.9mm 220Mev/c)
(Pb. model at 220: from MuSCAT siga=0.01, sige=0.01 )
1.0 55.198 12.21 28.561 1.000 184.1 0.0006893 0.0025247 0.0 ;
3.0 0.0019 ; (length)
9.0 0.0;
16.0 166.0 ; (Disable m/s)
16.0 161.0 ; (Slits are taken normal again.)

16. 11. 0.208582; (Final total momentum into Tracker)     

-6.0 6 10.0 ; (+/-10% Cut on momentum)

16.0 181.0 334 0; (Output muon particle vectors to for334.dat)

(-----Diagnostics for beam after Scatterer -----------------------------------)
(REAL SPACE)
51.0 1 -30 30 1.0 'XYmi';
52.0 3 -30 30 1.0 'XYmi';

(H PHASE SPACE)
-51.0 1 -30 30 1.0 'XDX'; (First column for rij matrix)
-52.0 2 -300 300 10.0 'XDX';
(V PHASE SPACE)
-51.0 3 -30 30 1.0 'YDY';
-52.0 4 -300 300 10 'YDY';
-50.0 6 -50 100 5.0; (plot dp/p)

SENTINEL
SENTINEL