Montecarlo method for particle simulation
8 views (last 30 days)
Show older comments
Hi all,
I have to simulate a number of particle with a initial random velocity along z axis emitted by a point source and distributed inside a solid angle passing through a pinhole. I would like to "kill" particles that cannot pass through the pinhole becuase of the angular spread and mantain those that can pass.
Is there any example I can use to write this simulation?
0 Comments
Accepted Answer
Image Analyst
on 8 Feb 2012
It's really not that hard. Have a for loop over a bunch of particles. For each particle, set its velocity (speed and angle) vector using the rand() function. Then use an "if" statement to decide whether the particle passed through or not and count those that do. It's like 7-10 lines or so.
0 Comments
More Answers (3)
Francesco
on 9 Feb 2012
1 Comment
Image Analyst
on 9 Feb 2012
You need to figure out what the distance is from the start, and where that particular particle would hit laterally the blocking screen/aperture plane, and decide if that location is within the pinhole locations or not.
Francesco
on 10 Feb 2012
2 Comments
Image Analyst
on 10 Feb 2012
I don't know. It's hard for me to follow since I don't understand Italian.
See Also
Categories
Find more on Particle & Nuclear Physics in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!