In [1]:

    

import numpy as np
import sys
import matplotlib.pyplot as plt
import os
import scipy.constants as ct
import lxml.etree as lxml
import sqlite3

from mpl_toolkits.mplot3d import Axes3D
%matplotlib inline
Ryd2eV=13.605692
epsilon=3
bohr2m=5.2917721092E-11


    

In [4]:

    

def readSqlall(sqlname):
        sqlstatement = "SELECT Jeff2s FROM pairs"
        con = sqlite3.connect(sqlname)
        with con:
                cur = con.cursor()
                cur.execute(sqlstatement)
                rows = cur.fetchall()

        sqlall=(np.array(rows)[:,0])
        
        return sqlall


    

In [6]:

    

pairs=readSqlall("system_individualmps.sql")


    

In [17]:

    

Jeff=np.sqrt(pairs)
Jeffbig=Jeff[Jeff>0.326099/2.0]
print np.sum(Jeffbig>0)
print np.sum(Jeff>0)


    

    




15
56151

In [18]:

    

n, bins, patches = plt.hist(Jeffbig, 50,  facecolor='green')
plt.xlabel('Coupling strength [eV]')
plt.ylabel('Counts')
plt.show()


    

In [6]:

    




    

In [7]:

    




    

In [13]:

    




    

In [7]:

    




    

In [ ]: