ProSHADE  0.6.5 (NOV 2018)
Protein Shape Descriptors and Symmetry Detection
All Classes Namespaces Files Functions Variables Pages
useDatabase.py
1 
4 
5 
7 import sys;
8 import os;
9 
10 cwd = os.getcwd ( );
11 proshadePath = cwd + "/../install/python3";
12 
13 sys.path.append ( proshadePath );
14 import proshade3 as proshade;
15 
16 
18 setUp = proshade.ProSHADE_settings ();
19 
20 
22 setUp.energyLevelDist = True;
23 setUp.traceSigmaDist = True;
24 setUp.fullRotFnDist = True;
25 
26 # Settings regarding hierarchical distance computation
27 setUp.enLevelsThreshold = -999.9;
28 setUp.trSigmaThreshold = -999.9;
29 
30 # ... Settings regarding the task
31 setUp.taskToPerform = proshade.Distances;
32 
33 # Settings regarding the database
34 setUp.databaseName = "";
35 setUp.volumeTolerance = 0.4;
36 setUp.clearMapData = False;
37 
38 # ... Settings regarding loudness
39 setUp.verbose = -1;
40 setUp.htmlReport = False;
41 
42 # Get command line info
43 if len(sys.argv) != 3:
44  print ( "Usage: python useDatabase.py [dbName] [filename1] to get distances from [filename1] to all files in the database." );
45  quit ( -1 );
46 else:
47  setUp.databaseName = str( sys.argv[1] );
48  hlpPyStrVec = proshade.StringList ( 1 );
49  hlpPyStrVec[ 0 ] = str( sys.argv[2] );
50 
51 # ... set the file list
52 setUp.structFiles = hlpPyStrVec;
53 
54 
56 runProshade = proshade.ProSHADE ( setUp );
57 
58 
60 crossCorrDists = runProshade.getCrossCorrDists ( );
61 traceSigmaDists = runProshade.getTraceSigmaDists ( );
62 rotationFunctionDists = runProshade.getRotFunctionDists ( );
63 
64 
66 sys.stdout.write ( "ProSHADE module version: " + runProshade.getProSHADEVersion() + "\n" );
67 sys.stdout.write ( "\n" );
68 sys.stdout.flush ( );
69 
70 if len(crossCorrDists) > 0:
71  ccOut = "Cross-correlation distances : " + "%+1.5f" % crossCorrDists[0];
72  for it in range ( 1, len(crossCorrDists) ):
73  ccOut = ccOut + "\t" + "%+1.5f" % crossCorrDists[it];
74  sys.stdout.write ( ccOut + "\n" );
75 sys.stdout.flush ( );
76 
77 if len(traceSigmaDists) > 0:
78  tsOut = "Trace sigma distances : " + "%+1.5f" % traceSigmaDists[0];
79  for it in range ( 1, len(traceSigmaDists) ):
80  tsOut = tsOut + "\t" + "%+1.5f" % traceSigmaDists[it];
81  sys.stdout.write ( tsOut + "\n" );
82 sys.stdout.flush ( );
83 
84 if len(rotationFunctionDists) > 0:
85  rfOut = "Rotation function distances : " + "%+1.5f" % rotationFunctionDists[0];
86  for it in range ( 1, len(rotationFunctionDists) ):
87  rfOut = rfOut + "\t" + "%+1.5f" % rotationFunctionDists[it];
88  sys.stdout.write ( rfOut + "\n" );
89 sys.stdout.flush ( );
90 
91 
93 quit ( );