Skip to content

Commit 3936a73

Browse files
authored
Merge pull request #1 from OpenSMFS/save_positions
Save particle position for each timestamp
2 parents dfd43d3 + ab82907 commit 3936a73

9 files changed

Lines changed: 1015 additions & 485 deletions

notebooks/PyBroMo - 1. Simulate 3D trajectories - single core.ipynb

Lines changed: 132 additions & 91 deletions
Large diffs are not rendered by default.

notebooks/PyBroMo - 2. Generate smFRET data, including mixtures.ipynb

Lines changed: 323 additions & 18 deletions
Large diffs are not rendered by default.

notebooks/PyBroMo - B.1 Disk-single-core - Generate photon timestamps.ipynb

Lines changed: 29 additions & 51 deletions
Original file line numberDiff line numberDiff line change
@@ -31,9 +31,7 @@
3131
{
3232
"cell_type": "code",
3333
"execution_count": null,
34-
"metadata": {
35-
"collapsed": false
36-
},
34+
"metadata": {},
3735
"outputs": [],
3836
"source": [
3937
"%matplotlib inline\n",
@@ -66,9 +64,7 @@
6664
{
6765
"cell_type": "code",
6866
"execution_count": null,
69-
"metadata": {
70-
"collapsed": false
71-
},
67+
"metadata": {},
7268
"outputs": [],
7369
"source": [
7470
"S = pbm.ParticlesSimulation.from_datafile('016', mode='w')"
@@ -77,9 +73,7 @@
7773
{
7874
"cell_type": "code",
7975
"execution_count": null,
80-
"metadata": {
81-
"collapsed": false
82-
},
76+
"metadata": {},
8377
"outputs": [],
8478
"source": [
8579
"def em_rates_from_E(em_rate_tot, E_values):\n",
@@ -97,9 +91,7 @@
9791
{
9892
"cell_type": "code",
9993
"execution_count": null,
100-
"metadata": {
101-
"collapsed": false
102-
},
94+
"metadata": {},
10395
"outputs": [],
10496
"source": [
10597
"em_rate_tot = 200e3\n",
@@ -112,9 +104,7 @@
112104
{
113105
"cell_type": "code",
114106
"execution_count": null,
115-
"metadata": {
116-
"collapsed": false
117-
},
107+
"metadata": {},
118108
"outputs": [],
119109
"source": [
120110
"# Get the random state at the end of the diffusion simulation\n",
@@ -136,9 +126,7 @@
136126
{
137127
"cell_type": "code",
138128
"execution_count": null,
139-
"metadata": {
140-
"collapsed": false
141-
},
129+
"metadata": {},
142130
"outputs": [],
143131
"source": [
144132
"em_rate_list"
@@ -166,9 +154,7 @@
166154
{
167155
"cell_type": "code",
168156
"execution_count": null,
169-
"metadata": {
170-
"collapsed": false
171-
},
157+
"metadata": {},
172158
"outputs": [],
173159
"source": [
174160
"%%timeit -n1 -r1\n",
@@ -188,9 +174,7 @@
188174
{
189175
"cell_type": "code",
190176
"execution_count": null,
191-
"metadata": {
192-
"collapsed": false
193-
},
177+
"metadata": {},
194178
"outputs": [],
195179
"source": [
196180
"%%timeit -n1 -r1\n",
@@ -203,9 +187,7 @@
203187
{
204188
"cell_type": "code",
205189
"execution_count": null,
206-
"metadata": {
207-
"collapsed": false
208-
},
190+
"metadata": {},
209191
"outputs": [],
210192
"source": [
211193
"for k in S.ts_store.h5file.root.timestamps._v_children.keys():\n",
@@ -216,20 +198,16 @@
216198
{
217199
"cell_type": "code",
218200
"execution_count": null,
219-
"metadata": {
220-
"collapsed": false
221-
},
201+
"metadata": {},
222202
"outputs": [],
223203
"source": [
224-
"ts, ts_par = S.get_timestamps_part('Pop1_P20_Pstart0_max_rate198000cps_BG4000cps_t_1s_rs_8798a6')"
204+
"ts, ts_par, _ = S.get_timestamp_data('Pop1_P20_Pstart0_max_rate198000cps_BG4000cps_t_1s_rs_8798a6')"
225205
]
226206
},
227207
{
228208
"cell_type": "code",
229209
"execution_count": null,
230-
"metadata": {
231-
"collapsed": false
232-
},
210+
"metadata": {},
233211
"outputs": [],
234212
"source": [
235213
"ts[:]"
@@ -245,9 +223,7 @@
245223
{
246224
"cell_type": "code",
247225
"execution_count": null,
248-
"metadata": {
249-
"collapsed": false
250-
},
226+
"metadata": {},
251227
"outputs": [],
252228
"source": [
253229
"bins = np.arange(0, 1, 1e-3)\n",
@@ -266,9 +242,7 @@
266242
{
267243
"cell_type": "code",
268244
"execution_count": null,
269-
"metadata": {
270-
"collapsed": false
271-
},
245+
"metadata": {},
272246
"outputs": [],
273247
"source": [
274248
"group = '/timestamps'\n",
@@ -284,9 +258,7 @@
284258
{
285259
"cell_type": "code",
286260
"execution_count": null,
287-
"metadata": {
288-
"collapsed": false
289-
},
261+
"metadata": {},
290262
"outputs": [],
291263
"source": [
292264
"[t for t in S.timestamp_names if 'BG4000cps' in t]"
@@ -295,9 +267,7 @@
295267
{
296268
"cell_type": "code",
297269
"execution_count": null,
298-
"metadata": {
299-
"collapsed": false
300-
},
270+
"metadata": {},
301271
"outputs": [],
302272
"source": [
303273
"S.ts_store.close()"
@@ -306,9 +276,7 @@
306276
{
307277
"cell_type": "code",
308278
"execution_count": null,
309-
"metadata": {
310-
"collapsed": false
311-
},
279+
"metadata": {},
312280
"outputs": [],
313281
"source": []
314282
},
@@ -338,9 +306,19 @@
338306
"name": "python",
339307
"nbconvert_exporter": "python",
340308
"pygments_lexer": "ipython3",
341-
"version": "3.5.0"
309+
"version": "3.6.7"
310+
},
311+
"toc": {
312+
"nav_menu": {},
313+
"number_sections": true,
314+
"sideBar": true,
315+
"skip_h1_title": false,
316+
"toc_cell": false,
317+
"toc_position": {},
318+
"toc_section_display": "block",
319+
"toc_window_display": false
342320
}
343321
},
344322
"nbformat": 4,
345-
"nbformat_minor": 0
323+
"nbformat_minor": 1
346324
}

0 commit comments

Comments
 (0)