v/del *
set *

set xval .2                             ! y-label to axis
g/SET LWID 8             | The line width is 6
g/set pass 2

null -10 10 10 35                       ! frame
text -4 8.3 'z vertex position (cm)' .4
text -11.8 12.7 'muon arm acceptance (angle in degrees)' .4 90
v/create left r -9.95

vector/create vertex(5) r left(1) -5. 0. 5. 10.

vector/create end_4_lo(5) r  7.4  8.7 10.6 13.4 18.3
vector/create end_4_hi(5) r 20.8 23.7 26.4 29.6 33.6
vector/create end_3_lo(5) r  6.2  7.0  8.2  9.8 12.3
vector/create end_3_hi(5) r 21.6 24.0 29.7 37.2 43.2
vector/create end_1_hi(5) r 25.5 28.6 32.3 37.1 48.9

v/crea pix_2_lo(5) r 15.1 19.2 26.2 39.7 69.4
v/crea pix_1_lo(5) r  8.0 10.2 14.3 23.2 54.1

v/crea strip_2_lo(5) r 28.3 32.7 38.5 46.2 56.6
v/crea strip_1_lo(5) r 23.1 27.5 33.7 42.9 56.8


*--------------- (2) Barrel pixel layers: ----------------------------*

   SET PLCI 4             | The line colour 2=red 4=blue
   SET LTYP 1             | The line type 2=dashed 1=solid
 pline 5 vertex pix_2_lo
 pline 5 vertex pix_1_lo

vector/create pix_2_y(5) r 
vector/create pix_2_x(5) r 
v/in pix_2_x(1)  vertex(1)
v/in pix_2_x(2)  vertex(2)
v/in pix_2_x(3)  vertex(3)
v/in pix_2_x(4)  vertex(4)
v/in pix_2_x(5)  left(1)

*v/in pix_2_y(1)  pix_2_lo(1)
*v/in pix_2_y(2)  pix_2_lo(2)
*v/in pix_2_y(3)  pix_2_lo(3)
*v/in pix_2_y(4)  pix_2_lo(4)
*v/in pix_2_y(5) 35.0

v/in pix_2_y(1)  end_3_hi(1)
v/in pix_2_y(2)  end_3_hi(2)
v/in pix_2_y(3)  end_3_hi(3)
v/in pix_2_y(4)  end_3_hi(4)
v/in pix_2_y(5) 35.0

graphics/set fais 3     ! fill area interior style 3=hatched
graphics/set faci 4     !           color index
graphics/set fasi 577                 style index  (?)
farea 5 pix_2_x pix_2_y

*--------------- (3) Barrel strip layers: -------------------------------*

   SET PLCI 7             | The line colour 2=red 4=blue
 pline 5 vertex strip_2_lo
 pline 5 vertex strip_1_lo

vector/create strip_2_y(5) r 
vector/create strip_2_x(5) r 
v/in strip_2_x(1)  vertex(1)
v/in strip_2_x(2)  vertex(2)
v/in strip_2_x(3)  vertex(3)
v/in strip_2_x(4)  vertex(1)

v/in strip_2_y(1)  strip_2_lo(1)
v/in strip_2_y(2)  strip_2_lo(2)
v/in strip_2_y(3)  strip_2_lo(3)
v/in strip_2_y(4)  35

graphics/set fais 3     ! fill area interior style 3=hatched
graphics/set faci 7     !           color index
graphics/set fasi 343                 style index  (?)
farea 4 strip_2_x strip_2_y

*--------------- (1) Muon endcaps: -----------------------------------*

   SET PLCI 6             | The line colour 2=red 4=blue
   SET LTYP 1             | The line type 2=dashed 1=solid
 pline 5 vertex end_4_lo
 pline 5 vertex end_4_hi


set plci 2
pline 5 vertex end_3_lo
pline 5 vertex end_3_hi
set ltyp 2
pline 5 vertex end_1_hi

vector/create end_4_y(10) r 
vector/create end_4_x(10) r 
v/in end_4_y(1)  end_4_hi(1)
v/in end_4_y(2)  end_4_hi(2)
v/in end_4_y(3)  end_4_hi(3)
v/in end_4_y(4)  end_4_hi(4)
v/in end_4_y(5)  end_4_hi(5)
v/in end_4_y(6)  end_4_lo(5)
v/in end_4_y(7)  end_4_lo(4)
v/in end_4_y(8)  end_4_lo(3)
v/in end_4_y(9)  end_4_lo(2)
v/in end_4_y(10)  10.0

v/in end_4_x(1)  vertex(1) 
v/in end_4_x(2)  vertex(2)
v/in end_4_x(3)  vertex(3)
v/in end_4_x(4)  vertex(4)
v/in end_4_x(5)  vertex(5)
v/in end_4_x(6)  vertex(5)
v/in end_4_x(7)  vertex(4)
v/in end_4_x(8)  vertex(3)
v/in end_4_x(9)  vertex(2)
v/in end_4_x(10)  left(1)

graphics/set fais 3     ! fill area interior style 3=hatched
graphics/set faci 6     !           color index
graphics/set fasi 344                 style index  (?)
farea 10 end_4_x end_4_y

v/del rx*
v/del ry*
v/crea rx1(1) r 4.0
v/crea ry1(1) r 18
v/crea rx2(1) r
v/crea ry2(1) r
vbias rx1(1) 4.2 rx2(1)
vbias ry1(1) 0.8 ry2(1)
v/crea rx(4) r rx1(1) rx2(1) rx2(1) rx1(1)
v/crea ry(4) r ry1(1) ry1(1) ry2(1) ry2(1)
g/set fais 1
g/set faci 8
g/set fasi 0
farea 4 rx ry
g/set txci 6
g/set pass 3
text rx(1)+0.2 ry(1)+0.2 '4 endcaps hit' .3

v/del rx*
v/del ry*
v/crea rx1(1) r 4.0
v/crea ry1(1) r 13.5
v/crea rx2(1) r
v/crea ry2(1) r
vbias rx1(1) 0.8 rx2(1)
vbias ry1(1) 0.8 ry2(1)
v/crea rx(4) r rx1(1) rx2(1) rx2(1) rx1(1)
v/crea ry(4) r ry1(1) ry1(1) ry2(1) ry2(1)
g/set fais 1
g/set faci 8
g/set fasi 0
farea 4 rx ry
g/set txci 6
g/set pass 3
text rx(1)+0.2 ry(1)+0.2 '4' .3

v/del rx*
v/del ry*
v/crea rx1(1) r 1.5
v/crea ry1(1) r 26.4
v/crea rx2(1) r
v/crea ry2(1) r
vbias rx1(1) 0.8 rx2(1)
vbias ry1(1) 0.8 ry2(1)
v/crea rx(4) r rx1(1) rx2(1) rx2(1) rx1(1)
v/crea ry(4) r ry1(1) ry1(1) ry2(1) ry2(1)
g/set fais 1
g/set faci 8
g/set fasi 0
farea 4 rx ry
g/set txci 6
g/set pass 3
text rx(1)+0.2 ry(1)+0.2 '4' .3

*----- endcaps between 3 hit ans 4 hit, lower bounds: ----------------*
v/crea end_3_x(6) r
v/crea end_3_y(6) r
v/in end_3_x(1) vertex(2) 
v/in end_3_x(2) vertex(3) 
v/in end_3_x(3) vertex(4) 
v/in end_3_x(4) vertex(5)
v/in end_3_x(5) vertex(5) 
v/in end_3_x(6) vertex(4) 

v/in end_3_y(1) end_4_lo(2)
v/in end_3_y(2) end_4_lo(3)
v/in end_3_y(3) end_4_lo(4)
v/in end_3_y(4) end_4_lo(5)
v/in end_3_y(5) end_3_lo(5)
v/in end_3_y(6) end_3_lo(4)

graphics/set fais 3     ! fill area interior style 3=hatched
graphics/set faci 2     !           color index
graphics/set fasi 344                 style index  (?)
farea 10 end_3_x end_3_y

*------ endcaps between 3 and 4 hits, upper bound: ----------------------*

*v/crea end_4_x(10) r
*v/crea end_4_y(10) r
v/in end_4_x(1) vertex(1) 
v/in end_4_x(2) vertex(2) 
v/in end_4_x(3) vertex(3) 
v/in end_4_x(4) vertex(4)
v/in end_4_x(5) vertex(5)

v/in end_4_x(6) vertex(5)
v/in end_4_x(7) vertex(4) 
v/in end_4_x(8) vertex(3) 
v/in end_4_x(9) vertex(2) 
v/in end_4_x(10) vertex(1) 

v/in end_4_y(1) end_4_hi(1)
v/in end_4_y(2) end_4_hi(2)
v/in end_4_y(3) end_4_hi(3)
v/in end_4_y(4) end_4_hi(4)
v/in end_4_y(5) end_4_hi(5)

v/in end_4_y(6) 35.0
v/in end_4_y(7) end_3_hi(4)
v/in end_4_y(8) end_3_hi(3)
v/in end_4_y(9) end_3_hi(2)
v/in end_4_y(10) end_3_hi(1)

graphics/set fais 3     ! fill area interior style 3=hatched
graphics/set faci 2     !           color index
graphics/set fasi 344                 style index  (?)
farea 10 end_4_x end_4_y

v/del rx*
v/del ry*
v/crea rx1(1) r -9.0
v/crea ry1(1) r 18.1
v/crea rx2(1) r
v/crea ry2(1) r
vbias rx1(1) 2.6 rx2(1)
vbias ry1(1) 0.8 ry2(1)
v/crea rx(4) r rx1(1) rx2(1) rx2(1) rx1(1)
v/crea ry(4) r ry1(1) ry1(1) ry2(1) ry2(1)
g/set fais 1
g/set faci 8
g/set fasi 0
farea 4 rx ry
g/set txci 4
g/set pass 3
text rx(1)+0.2 ry(1)+0.2 '2 pix hit' .3

v/del rx*
v/del ry*
v/crea rx1(1) r -5.5
v/crea ry1(1) r 20.0
v/crea rx2(1) r
v/crea ry2(1) r
vbias rx1(1) 0.8 rx2(1)
vbias ry1(1) 0.8 ry2(1)
v/crea rx(4) r rx1(1) rx2(1) rx2(1) rx1(1)
v/crea ry(4) r ry1(1) ry1(1) ry2(1) ry2(1)
g/set fais 1
g/set faci 8
g/set fasi 0
farea 4 rx ry
g/set txci 4
g/set pass 3
text rx(1)+0.2 ry(1)+0.2 '2' .3

v/del rx*
v/del ry*
v/crea rx1(1) r -4.0
v/crea ry1(1) r 12.0
v/crea rx2(1) r
v/crea ry2(1) r
vbias rx1(1) 0.8 rx2(1)
vbias ry1(1) 0.8 ry2(1)
v/crea rx(4) r rx1(1) rx2(1) rx2(1) rx1(1)
v/crea ry(4) r ry1(1) ry1(1) ry2(1) ry2(1)
g/set fais 1
g/set faci 8
g/set fasi 0
farea 4 rx ry
g/set txci 4
g/set pass 3
text rx(1)+0.2 ry(1)+0.2 '1' .3

v/del rx*
v/del ry*
v/crea rx1(1) r -9.0
v/crea ry1(1) r 33.0
v/crea rx2(1) r
v/crea ry2(1) r
vbias rx1(1) 3.0 rx2(1)
vbias ry1(1) 0.8 ry2(1)
v/crea rx(4) r rx1(1) rx2(1) rx2(1) rx1(1)
v/crea ry(4) r ry1(1) ry1(1) ry2(1) ry2(1)
g/set fais 1
g/set faci 8
g/set fasi 0
farea 4 rx ry
g/set txci 7
g/set pass 3
text rx(1)+0.2 ry(1)+0.2 '2 strip hit' .3

v/del rx*
v/del ry*
v/crea rx1(1) r -8.5
v/crea ry1(1) r 30.5
v/crea rx2(1) r
v/crea ry2(1) r
vbias rx1(1) 0.8 rx2(1)
vbias ry1(1) 0.8 ry2(1)
v/crea rx(4) r rx1(1) rx2(1) rx2(1) rx1(1)
v/crea ry(4) r ry1(1) ry1(1) ry2(1) ry2(1)
g/set fais 1
g/set faci 8
g/set fasi 0
farea 4 rx ry
g/set txci 7
g/set pass 3
text rx(1)+0.2 ry(1)+0.2 '2' .3

v/del rx*
v/del ry*
v/crea rx1(1) r -8.5
v/crea ry1(1) r 25.0
v/crea rx2(1) r
v/crea ry2(1) r
vbias rx1(1) 0.8 rx2(1)
vbias ry1(1) 0.8 ry2(1)
v/crea rx(4) r rx1(1) rx2(1) rx2(1) rx1(1)
v/crea ry(4) r ry1(1) ry1(1) ry2(1) ry2(1)
g/set fais 1
g/set faci 8
g/set fasi 0
farea 4 rx ry
g/set txci 7
g/set pass 3
text rx(1)+0.2 ry(1)+0.2 '1' .3

v/del rx*
v/del ry*
v/crea rx1(1) r -1.0
v/crea ry1(1) r 30.7
v/crea rx2(1) r
v/crea ry2(1) r
vbias rx1(1) 0.8 rx2(1)
vbias ry1(1) 0.8 ry2(1)
v/crea rx(4) r rx1(1) rx2(1) rx2(1) rx1(1)
v/crea ry(4) r ry1(1) ry1(1) ry2(1) ry2(1)
g/set fais 1
g/set faci 8
g/set fasi 0
farea 4 rx ry
g/set txci 6
g/set pass 3
text rx(1)+0.2 ry(1)+0.2 '1' .3


v/del rx*
v/del ry*
v/crea rx1(1) r 8.0
v/crea ry1(1) r 12.0
v/crea rx2(1) r
v/crea ry2(1) r
vbias rx1(1) 0.8 rx2(1)
vbias ry1(1) 0.8 ry2(1)
v/crea rx(4) r rx1(1) rx2(1) rx2(1) rx1(1)
v/crea ry(4) r ry1(1) ry1(1) ry2(1) ry2(1)
g/set fais 1
g/set faci 8
g/set fasi 0
farea 4 rx ry
g/set txci 6
g/set pass 3
text rx(1)+0.2 ry(1)+0.2 '3' .3

v/del rx*
v/del ry*
v/crea rx1(1) r -1
v/crea ry1(1) r 27.5
v/crea rx2(1) r
v/crea ry2(1) r
vbias rx1(1) 0.8 rx2(1)
vbias ry1(1) 0.8 ry2(1)
v/crea rx(4) r rx1(1) rx2(1) rx2(1) rx1(1)
v/crea ry(4) r ry1(1) ry1(1) ry2(1) ry2(1)
g/set fais 1
g/set faci 8
g/set fasi 0
farea 4 rx ry
g/set txci 6
g/set pass 3
text rx(1)+0.2 ry(1)+0.2 '3' .3

v/crea all_2_x(6) r
v/crea all_2_y(6) r
v/crea pointx(1) r -2.6
v/crea pointy(1) r 30.5

v/in all_2_x(1) pointx(1)
v/in all_2_x(2) vertex(3)
v/in all_2_x(3) vertex(4)

v/in all_2_x(4) vertex(4)
v/in all_2_x(5) vertex(3)
v/in all_2_x(6) vertex(2)

v/in all_2_y(1) pointy(1)
v/in all_2_y(2) end_1_hi(3)
v/in all_2_y(3) end_1_hi(4)

v/in all_2_y(4) strip_1_lo(4)
v/in all_2_y(5) strip_1_lo(3)
v/in all_2_y(6) strip_1_lo(2)

graphics/set fais 3     ! fill area interior style 3=hatched
graphics/set faci 1     !           color index
graphics/set fasi 177                 style index  (?)
farea 6 all_2_x all_2_y