Commit: e716cbc21b8d2c7c5e49f401ddb166b65043fe43 - myprojects-hg-reborn (hg)

Commit MetaInfo

• Revision: e716cbc21b8d2c7c5e49f401ddb166b65043fe43 (tree)
• Zeit: 2010-10-20 21:49:55
• Autor: lorenzo
• Commiter: lorenzo

Log Message

I added a new filter on the tag IDs (useful when dealing with compartimentized
interactions).

Ändern Zusammenfassung

• modified: Python-codes/t_ab_ac_from_scratch.py

Diff

diff -r cdd2c30a8152 -r e716cbc21b8d Python-codes/t_ab_ac_from_scratch.py
--- a/Python-codes/t_ab_ac_from_scratch.py	Wed Oct 20 12:47:00 2010 +0000
+++ b/Python-codes/t_ab_ac_from_scratch.py	Wed Oct 20 12:49:55 2010 +0000
@@ -105,6 +105,27 @@
 
     return (all_times,tags_abc)
     
+
+def p_tab_tac_sel(info,selection):
+    tag_list=s.unique1d(s.ravel(info[:,1:3]))
+    all_times=s.zeros(0).astype("int64")
+    tags_abc=s.zeros((0,4)).astype("int64")
+
+    for i in xrange(len(tag_list)):
+        tag_id=tag_list[i]
+
+        if (tag_id in selection):
+
+            tags_info=single_tag_contact_times(info, tag_id)
+
+            delta_t=tags_info[2]
+            abc=tags_info[3]
+            all_times=s.hstack((all_times, delta_t))
+            tags_abc=s.vstack((tags_abc, abc))
+
+    return (all_times,tags_abc)
+    
+
     
     
 def p_tab_tac_filtered(info):
@@ -168,8 +189,8 @@
 
 #Now I am adding an extra feature (i.e. I am cutting a slice of the data)
 
-t_ini=1275283490
-t_fin=1275456290
+t_ini= 1257417888 #    1275283490
+t_fin= 1258104328                #   1275456290
 
 sel=s.where((info[:,0]>=t_ini) & (info[:,0]<=t_fin) )
 
@@ -184,11 +205,18 @@
 #############################################################################################
 # Now iterate on all tags
 
-results=p_tab_tac(info)
+# results=p_tab_tac(info)
 
 #results=p_tab_tac_filtered(info)
 
 
+sel_p = s.arange(1400,1500)
+
+sel_else = s.arange(1200,1300)
+
+results=p_tab_tac_sel(info, sel_else)
+
+
 n.savetxt("delta_tab_tac.dat",results[0]  ,  fmt='%d')
 
 n.savetxt("tags_ab_ac.dat",results[1]  ,  fmt='%d')