macro triangle 
vec/dele * 
vec/cre babar(1)
vec/inp babar 0. 
vec/cre radplus(1) 
vec/cre radminus(1) 
vec/cre cenplus(1) 
vec/cre cenminus(1) 
vec/cre rho(100) 
vec/cre eta(100) 
vec/cre eta_up
vec/inp eta_up 1.

*----  Initialisation 
call contrainte.for(0.)
zone
set ymgu 5
set ymgl 7
null -1. 1. 0 1.



*----  Dmd
call contrainte.for(1.)
igset fais 3
igset faci 2
igset plci 1
igset fasi 345
farea 100 rho eta
pline 100 rho eta
igset plci 2 
dbox -1. -.8 1.5 1.6
text -.75 1.5 '[D]m?d! constraint' .3
*----  Epsilon K 
call contrainte.for(2.)
igset fais 3
igset faci 3
igset plci 1
igset fasi 395
farea 100 rho eta
pline 100 rho eta
dbox -1. -.8 1.3 1.4
text -.75 1.3 '[e]?K! constraint' .3

*----  Vbu/Vcb
call contrainte.for(3.)
igset fais 3
igset faci 7
igset plci 1 
igset fasi 353
farea 100 rho eta
pline 100 rho eta
dbox -1. -.8 1.1 1.2
text -.75 1.1 'V?ub!/V?cb! constraint' .3

if $sigma(babar(1)) = 0 goto nobabar



*----  Sin(2b)
call contrainte.for(4.)
igset fais 3
igset faci 6
igset plci 1 
igset fasi 305
farea 100 rho eta
pline 100 rho eta
dbox -.1 .1 1.5 1.6
text .15 1.5 'sin(2[b]) constraint' .3

vec/cre rho(200) 
vec/cre eta(200) 
*----  Sin(2a)
call contrainte.for(5.)
igset fais 3
igset faci 4
igset plci 1 
igset fasi 305
farea 200 rho eta
graph/primitives/arc .5 cenminus(1) radminus(1) 
graph/primitives/arc .5 cenplus(1) radplus(1) 
dbox -.1 .1 1.3 1.4
text .15 1.3 'sin(2[a]) constraint' .3


nobabar:

atit '[r]' '[c]'
return