18
|
1 from commons.core.sql.DbFactory import DbFactory
|
|
2 from commons.core.utils.FileUtils import FileUtils
|
|
3 from commons.pyRepetUnit.convCoord.ConvSetChr2Chunk import ConvSetChr2Chunk
|
|
4 import unittest
|
|
5 import os
|
|
6
|
|
7 class Test_ConvSetChr2Chunk(unittest.TestCase):
|
|
8
|
|
9 def setUp(self):
|
|
10 self._setfilename="dummy.set"
|
|
11 self._settablename="dummy_set"
|
|
12
|
|
13 self._chunkmapfilename="chunk.map"
|
|
14 self._chunkmaptable="test_chunk_map"
|
|
15
|
|
16 self._resfilename="dummy_res.set"
|
|
17
|
|
18 map4convData="chunk1\tchr1\t1\t100\n"
|
|
19 map4convData+="chunk2\tchr1\t90\t200\n"
|
|
20 map4convData+="chunk3\tchr1\t190\t300\n"
|
|
21
|
|
22 map4conv=open(self._chunkmapfilename,"w")
|
|
23 map4conv.write(map4convData)
|
|
24 map4conv.close()
|
|
25
|
|
26 self._db = DbFactory.createInstance()
|
|
27 self._db.createTable(self._chunkmaptable, "map", self._chunkmapfilename, True)
|
|
28
|
|
29 set2convData="1\treg1\tchr1\t10\t80\n"
|
|
30 set2convData+="1\treg1\tchr1\t90\t110\n"
|
|
31 set2convData+="2\treg2\tchr1\t150\t160\n"
|
|
32
|
|
33 set2convData+="3\treg3\tchr1\t70\t210\n"
|
|
34 set2convData+="4\treg4\tchr1\t85\t100\n"
|
|
35 set2convData+="5\treg5\tchr1\t85\t110\n"
|
|
36
|
|
37 set2convData+="6\treg6\tchr1\t80\t10\n"
|
|
38 set2convData+="6\treg6\tchr1\t110\t90\n"
|
|
39 set2convData+="7\treg7\tchr1\t160\t150\n"
|
|
40 set2convData+="8\treg8\tchr1\t190\t300\n"
|
|
41
|
|
42 set2convData+="9\treg9\tchr1\t210\t70\n"
|
|
43 set2convData+="10\treg10\tchr1\t100\t85\n"
|
|
44 set2convData+="11\treg11\tchr1\t110\t85\n"
|
|
45 set2convData+="12\treg12\tchr1\t300\t190\n"
|
|
46
|
|
47 set2conv=open(self._setfilename,"w")
|
|
48 set2conv.write(set2convData)
|
|
49 set2conv.close()
|
|
50
|
|
51 self._db.createTable(self._settablename, "set", self._setfilename, True)
|
|
52
|
|
53 setResData="1\treg1\tchunk1\t10\t80\n"
|
|
54 setResData+="1\treg1\tchunk1\t90\t100\n"
|
|
55 setResData+="1\treg1\tchunk2\t1\t21\n"
|
|
56 setResData+="2\treg2\tchunk2\t61\t71\n"
|
|
57
|
|
58 setResData+="3\treg3\tchunk1\t70\t100\n"
|
|
59 setResData+="3\treg3\tchunk2\t1\t111\n"
|
|
60 setResData+="3\treg3\tchunk3\t1\t21\n"
|
|
61
|
|
62 setResData+="4\treg4\tchunk1\t85\t100\n"
|
|
63 setResData+="4\treg4\tchunk2\t1\t11\n"
|
|
64 setResData+="5\treg5\tchunk1\t85\t100\n"
|
|
65 setResData+="5\treg5\tchunk2\t1\t21\n"
|
|
66
|
|
67 setResData+="6\treg6\tchunk1\t80\t10\n"
|
|
68 setResData+="6\treg6\tchunk1\t100\t90\n"
|
|
69 setResData+="6\treg6\tchunk2\t21\t1\n"
|
|
70 setResData+="7\treg7\tchunk2\t71\t61\n"
|
|
71 setResData+="8\treg8\tchunk2\t101\t111\n"
|
|
72 setResData+="8\treg8\tchunk3\t1\t111\n"
|
|
73
|
|
74 setResData+="9\treg9\tchunk1\t100\t70\n"
|
|
75 setResData+="9\treg9\tchunk2\t111\t1\n"
|
|
76 setResData+="9\treg9\tchunk3\t21\t1\n"
|
|
77
|
|
78 setResData+="10\treg10\tchunk1\t100\t85\n"
|
|
79 setResData+="10\treg10\tchunk2\t11\t1\n"
|
|
80 setResData+="11\treg11\tchunk1\t100\t85\n"
|
|
81 setResData+="11\treg11\tchunk2\t21\t1\n"
|
|
82 setResData+="12\treg12\tchunk2\t111\t101\n"
|
|
83 setResData+="12\treg12\tchunk3\t111\t1\n"
|
|
84
|
|
85 setRes=open(self._resfilename,"w")
|
|
86 setRes.write(setResData)
|
|
87 setRes.close()
|
|
88
|
|
89 def tearDown(self):
|
|
90 self._db.dropTable(self._chunkmaptable)
|
|
91 self._db.dropTable(self._settablename)
|
|
92 self._db.close()
|
|
93 os.remove(self._resfilename)
|
|
94 os.remove(self._setfilename)
|
|
95 os.remove(self._chunkmapfilename)
|
|
96
|
|
97 def test_convert(self):
|
|
98 conv = ConvSetChr2Chunk(self._db,self._settablename, self._chunkmaptable, self._settablename+"_onchk")
|
|
99 conv.convert()
|
|
100 obsFileName = "dummy_set.onchk"
|
|
101 self._db.exportDataToFile("dummy_set_onchk", obsFileName)
|
|
102 self._db.dropTable("dummy_set_onchk")
|
|
103 self.assertTrue(FileUtils.are2FilesIdentical(self._resfilename, obsFileName))
|
|
104 os.remove(obsFileName)
|
|
105
|
|
106 test_suite = unittest.TestSuite()
|
|
107 test_suite.addTest( unittest.makeSuite( Test_ConvSetChr2Chunk ) )
|
|
108 if __name__ == '__main__':
|
|
109 unittest.TextTestRunner(verbosity=2).run( test_suite ) |