an example with plots added

dream/simulation/Examples/TwoServersPlots.py

+from SimPy.Simulation import simulate, activate, initialize
+from simulation.Machine import Machine
+from simulation.Source import Source
+from simulation.Exit import Exit
+from simulation.Part import Part
+from simulation.Repairman import Repairman
+from simulation.Queue import Queue
+from simulation.Globals import G
+
+#import Graphs
+from outputanalysis.Plots import Graphs
+
+#define the objects of the model
+R=Repairman('R1', 'Bob') 
+S=Source('S1','Source', mean=0.5, item=Part)
+M1=Machine('M1','Machine1', mean=0.25, failureDistribution='Fixed', MTTF=60, MTTR=5, repairman=R)
+Q=Queue('Q1','Queue')
+M2=Machine('M2','Machine2', mean=1.5, failureDistribution='Fixed', MTTF=40, MTTR=10,repairman=R)
+E=Exit('E1','Exit')  
+
+G.ObjList=[S,M1,M2,E,Q]   #add all the objects in G.ObjList so that they can be easier accessed later
+
+#define predecessors and successors for the objects    
+S.defineRouting([M1])
+M1.defineRouting([S],[Q])
+Q.defineRouting([M1],[M2])
+M2.defineRouting([Q],[E])
+E.defineRouting([M2])
+
+initialize()                        #initialize the simulation (SimPy method)
+    
+#initialize all the objects    
+R.initialize()
+
+for object in G.ObjList:
+    object.initialize()
+
+#activate all the objects 
+for object in G.ObjList:
+    activate(object, object.run())
+
+G.maxSimTime=1440.0     #set G.maxSimTime 1440.0 minutes (1 day)
+    
+simulate(until=G.maxSimTime)    #run the simulation
+
+#carry on the post processing operations for every object in the topology       
+for object in G.ObjList:
+    object.postProcessing()
+R.postProcessing()
+
+#print the results
+print "the system produced", E.numOfExits, "parts"
+print "the blockage ratio of", M1.objName,  "is", (M1.totalBlockageTime/G.maxSimTime)*100, "%"
+print "the working ratio of", R.objName,"is", (R.totalWorkingTime/G.maxSimTime)*100, "%"
+
+#calculate the percentages for the pie
+repairmanWorkingRatio=R.totalWorkingTime/G.maxSimTime*100
+repairmanWaitingRatio=R.totalWaitingTime/G.maxSimTime*100
+
+#create a graph object
+graph=Graphs()
+#create the pie
+graph.Pie([repairmanWorkingRatio,repairmanWaitingRatio], "repairman.jpg")
+
+
+