def getdMu(h5pFile):
TINY = 1.0e-8
pIds, Muf = lfmpp.getH5pFin(h5pFile, "Mu")
pIds, Mu0 = lfmpp.getH5pInit(h5pFile, "Mu")
Ind = (Muf > TINY) & (Mu0 > TINY)
Muf = Muf[Ind]
Mu0 = Mu0[Ind]
#dMu = np.abs(Muf-Mu0)/Mu0
dMu = (Muf - Mu0) / Mu0
print("\nCutting %d particles for undefined Mu" % (len(Ind) - Ind.sum()))
print("Min dMu = %f" % (dMu.min()))
print("Min Muf = %f, Min Mu0 = %f" % (Muf.min(), Mu0.min()))
print("\n\n")
return dMu, Ind
def getA0(h5pFile):
pIds, A0 = lfmpp.getH5pInit(h5pFile, "alpha")
return A0
def getdK(h5pFile):
pIds, Kf = lfmpp.getH5pFin(h5pFile, "kev")
pIds, K0 = lfmpp.getH5pInit(h5pFile, "kev")
dK = Kf / K0
return dK