This project has moved and is read-only. For the latest updates, please go here.

mc.exe timeout/crash for long-running process

Jun 1, 2011 at 8:20 PM
Edited Jun 1, 2011 at 8:23 PM

A summary of email conversation with Rob:

Rob started 10M photons it at ~5:15pm yesterday, expecting the simulation to take 25 hours. Here is a snapshot:

Jun 1, 2011 at 9:33 PM
Edited Jun 1, 2011 at 9:43 PM

An update from Carole:

After scanning the input file, I see one potential problem. The input file specifies the creation of two databases and those file are probably getting *very* large for 10M.

I have edited your 10K input file and made that option nil and attached it (Note to self, change default infile to not specify the creation of databases).

Please update N and give this a try for 1M and 10M.


Aside: not sure how to deal with attachments and images in this discussion tab - maybe we should adopt a different online service.

In the meantime, I recommend using and for code and screenshots, respectvely.

Here's Carole's XML code:

Jun 1, 2011 at 11:16 PM

I started a new launch of 1M photons at 2:15pm (with Carole's new .xml file) and it seems on track to finish around 4:45pm (it just announced 40% completion).  I'll post back soon when/if it finishes.

Jun 1, 2011 at 11:19 PM

Also, in response to Carole... the two databases left behind after the 10M launch are 95MB and 78MB.

I'm not sure if that's what you guys call "*very* large"... I think that's a few albums worth of .mp3s though.

Jun 2, 2011 at 2:54 AM

The simulation completed successfully and I started a 10M photon launch.  In the meantime, I attempted to utilize David's vts_parser (v2) and you can see the results here.  I'm not sure what I'm looking at... does this look correct?


Also, I recommend for uploading pictures ( wouldn't upload my .jpgs for some reason) :D

Jun 2, 2011 at 4:37 AM
Thanks for the tip!

You're looking at 2D maps of R(rho,omega). Vertical axis is the first variable, rho; horizontal axis is the second, omega. The amplitude map is plotted
in log(A), phase is in radians. You can look at all the results generated by typing 'whos' in Matlab. Note the 'results' variable.
Typing 'results' will show you all outputs (detectors) generated by your simulation run, including the full 2D arrays. The midpoints of the
bins are pre-calculated for you, so that you can access values at a constant rho or omega, for example,

'figure; plot(results.ROfRhoAndOmega.Omega_Midpoints, results.ROfRhoAndOmega.Amplitude(10, :));'

will plot the amplitude at the 10th rho bin, which would be at


I typed that from my iPhone, so caveat emptor...let me know how it goes.
Jun 2, 2011 at 5:23 PM
This discussion has been copied to a work item. Click here to go to the work item and continue the discussion.