import random
a = 50
length = 2
dl = 0.05
globalwidth = 1.5
degree = 70
def Start():
   #random.seed(150)
   pass
   
Axiom: SetGuide(treepath,90)_(1)@GsG(0.5)A(1,0.30)
derivation length: 100
production:
A(x,k) :
  if k < 3.5:
    if x <= 7:
      produce /(96)F(0.7,1-k/5)A(x+1,k+0.01)
    elif x <= 14:
      produce /(degree)F(0.7,1-k/5)[SetGuide(treepath2,8)_(0.2)&(90)-(60)B(1,k)]A(x+1,k+0.1)
    else:
      x = 0
      produce A(x+1,k+0.01)
  elif k < 4:
    produce /(degree)F(0.7,1-k/5)[SetGuide(finish,8)_(0.2)&(90)-(60)B(1,k)]A(x+1,k+0.1)
B(x,k):
  r = random.random()
  if x <= 20: 
    produce F(0.8,0.05+0.25/k)[/(r*500)+(20)_(0.02)F(1)C(1,1.0)]B(x+1,k+1)
C(x,k):
  r = random.random()
  if x <= 3:
    produce F(0.6)C(x+1,k+1);(2)&(r*200)Leaf()

Leaf():
  produce  Sweep(path,section,length,dl,globalwidth,leafwidth)
interpretation:
endlsystem
###### INITIALISATION ######

__lpy_code_version__ = 1.1

def __initialiseContext__(context):
	import openalea.plantgl.all as pgl
	Color_1 = pgl.Material("Color_1" , ambient = (66,33,0) , diffuse = 0 , specular = (204,204,204) , )
	Color_1.name = "Color_1"
	context.turtle.setMaterial(1,Color_1)
	Color_2 = pgl.Material("Color_2" , ambient = (6,13,2) , diffuse = 5.76923 , )
	Color_2.name = "Color_2"
	context.turtle.setMaterial(2,Color_2)
	import openalea.plantgl.all as pgl
	path = pgl.BezierCurve2D(	
	    pgl.Point3Array([(-0.5, 0, 1),(-0.326083, -0.0978345, 1),(-0.0956592, 0.0176126, 1),(0.0620721, -0.00621237, 1)]) , 
	    )
	path.name = "path"
	section = pgl.NurbsCurve2D(	
	    ctrlPointList = pgl.Point3Array([(-0.672727, 0, 1),(-0.324715, 0.227106, 1),(0.023332, -0.20146, 1),(0.337225, -0.0436618, 1),(0.695455, 0.00454545, 1)]) , 
	    )
	section.name = "section"
	import openalea.plantgl.all as pgl
	leafwidth = pgl.NurbsCurve2D(	
	    ctrlPointList = pgl.Point3Array([(0, 0, 1),(0.0305872, 0.539066, 1),(0.125225, 0.681196, 1),(0.344686, 0.720503, 1),(0.544193, 0.504133, 1),(0.695718, 0.293404, 1),(0.731597, 0.166532, 1),(0.771218, 0.0592187, 1),(0.897377, 0.030847, 1),(1, -0.000465795, 1)]) , 
	    )
	leafwidth.name = "leafwidth"
	treepath = pgl.NurbsCurve2D(	
	    ctrlPointList = pgl.Point3Array([(-0.5, 0, 1),(-0.238627, 0.00619407, 1),(-0.07298, -0.00805499, 1),(0.0491646, 0.0141531, 1),(0.234358, -9.59484e-005, 1),(0.398942, -0.00805499, 1),(0.565017, -0.00943561, 1)]) , 
	    )
	treepath.name = "treepath"
	finish = pgl.BezierCurve2D(	
	    pgl.Point3Array([(-0.465217, 0.00434783, 1),(-0.310145, 0.0826087, 1),(-0.324638, 0.456522, 1),(-0.408696, 0.669565, 1)]) , 
	    )
	finish.name = "finish"
	treepath2 = pgl.NurbsCurve2D(	
	    ctrlPointList = pgl.Point3Array([(-0.00917746, 0.0133432, 1),(0.72348, 1.57179, 1),(3.10364, 2.97222, 1),(3.50129, 4.69656, 1),(3.82637, 6.13593, 1)]) , 
	    )
	treepath2.name = "treepath2"
	panel_0 = ({'active': True, 'visible': True, 'name': 'Panel 1'},[('Curve2D',path),('Curve2D',section),('Function',leafwidth),('Curve2D',treepath),('Curve2D',finish),('Curve2D',treepath2)])
	parameterset = [panel_0,]
	context["__functions__"] = [('leafwidth',leafwidth),]
	context["__curves__"] = [('path',path),('section',section),('treepath',treepath),('finish',finish),('treepath2',treepath2),]
	context["__parameterset__"] = parameterset
	context["path"] = path
	context["section"] = section
	context["leafwidth"] = pgl.QuantisedFunction(leafwidth)
	context["treepath"] = treepath
	context["finish"] = finish
	context["treepath2"] = treepath2