valid_roket_files Namespace Reference

Functions
def	validfile (filename)

Variables
string	datapath = "/home/fferreira/Data/correlation/"
	filenames = glob.glob(datapath + "roket_8m*.h5")
	l = len(filenames)
int	ind = 0
	end

Function Documentation

validfile()

def valid_roket_files.validfile

(

filename

)

Definition at line 6 of file valid_roket_files.py.

def validfile(filename):
    f = h5py.File(filename)
    if (list(f.attrs.keys()).count("target.Lambda")):
        Lambda = f.attrs["target.Lambda"][0]
    else:
        Lambda = 1.65
    nactus = f["noise"][:].shape[0]
    niter = f["noise"][:].shape[1]
    P = f["P"][:]
    nmodes = P.shape[0]
    data = np.zeros((nmodes, niter))
    error_list = [
            "noise", "aliasing", "tomography", "filtered modes", "bandwidth",
            "non linearity"
    ]
 
    for i in error_list:
        data += np.dot(P, f[i][:])
 
    data = np.var(data, axis=1)
    data = np.sum(data)
    data = np.exp(-data * (2 * np.pi / Lambda)**2)
    data *= np.exp(-f["fitting"].value)
 
    SR2 = f["SR2"].value
    SR = f["SR"].value
 
    if (np.abs(data - SR) < 0.05 or np.abs(data - SR2) < 0.05):
        f.attrs["validity"] = True
 
 

Variable Documentation

datapath

string valid_roket_files.datapath = "/home/fferreira/Data/correlation/"

Definition at line 37 of file valid_roket_files.py.

end

valid_roket_files.end

Definition at line 44 of file valid_roket_files.py.

filenames

valid_roket_files.filenames = glob.glob(datapath + "roket_8m*.h5")

Definition at line 38 of file valid_roket_files.py.

ind

int valid_roket_files.ind = 0

Definition at line 40 of file valid_roket_files.py.

l

valid_roket_files.l = len(filenames)

Definition at line 39 of file valid_roket_files.py.