#! /usr/bin/python2.7

import os,subprocess,shlex

if __name__ == "__main__":
	inFile=open("PDBsequences.csv",'r')
	outFile=open("ExpNumNbor_300",'w')
	outFile.write("pdb\tUniform_MS1\tUniform_MS1_norm\tUniform_MS2\tUniform_MS2_norm\tTurner_MS1\tTurner_MS1_norm\tTurner_MS2\tTurner_MS2_norm\n")
	for line in inFile.readlines():
		rna = line.strip().split('\t')[1]
		pdbId = line.strip().split('\t')[0]
		print rna
		cmdUniformMS1 = "/home/amir/ExpNumNhborswithShift/Cprogram_ExpNumNborsAddRemShiftUnifProb/nbor0 %s" %rna
		cmdUniformMS2 = "/home/amir/ExpNumNhborswithShift/Cprogram_ExpNumNborsAddRemShiftUnifProb/nbor1 %s" %rna
		cmdTurnerMS1 = "/home/amir/ExpNumNhborswithShift/srcRNAexpNumNbors_Cprograms/a.out %s" %rna
		cmdTurnerMS2 = "/home/amir/ExpNumNhborswithShift/Turner-Shift/nbor %s" %rna
		UMS1=subprocess.Popen(shlex.split(cmdUniformMS1),stdout=subprocess.PIPE)
		outUMS1,errMS1 = UMS1.communicate()
		UMS2 = subprocess.Popen(shlex.split(cmdUniformMS2),stdout=subprocess.PIPE)
		outUMS2,errMSt= UMS2.communicate()
		TMS1 = subprocess.Popen(shlex.split(cmdTurnerMS1),stdout=subprocess.PIPE)
		outTMS1,errTMS1 = TMS1.communicate()
		TMS2 = subprocess.Popen(shlex.split(cmdTurnerMS2),stdout = subprocess.PIPE)
		outTMS2,errTMS2 = TMS2.communicate()
		outFile.write("%s\t%s\t%s\t%s\t%s\n"%(pdbId,outUMS1.strip(),outUMS2.strip(),outTMS1.strip(),outTMS2.strip()))
	inFile.close()
	outFile.close()