notebook.community



In [34]:

    
import IMP
import IMP.pmi.output as output



In [35]:

    
po=IMP.pmi.output.ProcessOutput("stat_replica.0.out")
po.show_keys()









    



ReplicaExchange_CurrentTemp                                      ReplicaExchange_SwapSuccessRatio                                 
ReplicaExchange_MaxTempFrequency                                 score                                                            
ReplicaExchange_MinTempFrequency



In [36]:

    
fs=po.get_fields(["ReplicaExchange_SwapSuccessRatio","score","ReplicaExchange_CurrentTemp"])



In [37]:

    
IMP.pmi.output.plot_fields(fs)

Plot score histograms for all replicas



In [ ]:

    
values=[]
for i in range(0,64):
    po=IMP.pmi.output.ProcessOutput("stat_replica."+str(i)+".out")
    values.append(po.get_fields(["score"])["score"])



In [ ]:

    
output.plot_field_histogram("scores",values,yplotrange=[0,0.02])



In [ ]: