'*PLOT*MICE Muon Beamline Model for TURTLE'

(Reference version 19 Feb 2007. Includes TOF0/Ckov scattering patches, but no 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:- 134862 from ref. File2.)
134862

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)

1. -32.0 0.0 0.0 0.0 0.0 0.0 0.258098 ; (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 m/scattering)

16.0 11.0 0.258098 ; (central momentum 0.258098)


-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)
(Dipole Aperture + Fringe Field)
16.0 4.0 33.0 ; (Dipole hz half ap)
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. 1.04995 4.293258339 0. ;
2.0 15.0 ;
20.0 -180.0 ;
-51.0 1 -40 30 1.0 'afB2_X';
-52.0 3 -30 30 1.0 'afB2_Y';


3.0 0.25;
3.0 0.45;
3.0 0.3;

(---- 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';
-51.0 1 -30 30 1.0 'afQ4_3c_X';
-52.0 3 -30 30 1.0 'afQ4_3c_Y';

3. 0.5;

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';
-51.0 1 -30 30 1.0 'afQ5_2c_X';
-52.0 3 -30 30 1.0 'afQ5_2c_Y';

3. 0.5;

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'; 
-51.0 1 -25 25 1.0 'afQ6_1c_X';
-52.0 3 -25 25 1.0 'afQ6_1c_Y';

3.0 0.25;
3.0 0.25;

-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 ; (MAR05 Setting. Not quite acc but ok)

(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)

(Diagnostic after TOF0 device)
-51.0 1 -25 25 1.0 'XYT0';
-52.0 3 -25 25 1.0 'XYT0';
-51.0 1 -30 30 1.0 'afT0_X';
-52.0 3 -25 25 1.0 'afT0_Y';

3.0 0.2;

(---- Ckov1 Module here: consists of 4parts. ----)
(NB Model is for MAR05 setting.)
-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 ; (MAR05 Setting. Not quite acc but ok)

(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 ; (MAR05 Setting. Not quite acc but ok)

(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 ; (MAR05 Setting. Not quite acc but ok)

(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 ; (MAR06 Setting. 7.1pi     Not quite acc but ok)

(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 ; (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)

(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 0.5;

(This is the extra 2.0m since AUG05:- )
3.0 0.5 ;
-51.0 1 -85 85 2.0 'XY1_X';
-52.0 3 -85 85 2.0 'XY1_Y';
3.0 1.0 ;
-51.0 1 -85 85 2.0 'XY2_X';
-52.0 3 -85 85 2.0 'XY2_Y';
3.0 0.5 ;

3. 0.9003;


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

16.0  100.0 23.6 23.6;
16.0  101.0 17.82 0.0;
-51.0 1 -85 85 2.0 'b4Q7_3_X';
-52.0 3 -85 85 2.0 'b4Q7_3_Y';
5.0 0.66 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.5;

16.0  100.0 23.6 23.6;
16.0  101.0 17.82 0.0;
-51.0 1 -50 50 1.0 'b4Q8_4_X';
-52.0 3 -50 50 1.0 'b4Q8_4_Y';
5.0 0.66 -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.5;

16.0  100.0 23.6 23.6;
16.0  101.0 17.82 0.0;
-51.0 1 -30 30 1.0 'b4Q9_5_X';
-52.0 3 -30 30 1.0 'b4Q9_5_Y';
5.0 0.66 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. 1.0;

-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 ; (MAR06 Setting.   Not quite acc but ok)


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


3.0 0.467;


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

(---- Diffuser here: ----)
(Scatterer. 9.5mm )
(Designed for 0.220 Gev/c)
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 25.0 3.0 25.0 ; (A 50cm x 50cm area)
(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)


(-----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