I have built a detector for the detection of the cosmic-ray muon. Till now, I was using point particle gun for shooting high energy muons into the detector. Now I want to create a rectangular plane above the detector to throw thousands on particles in a single event by covering the whole detector plane. It is easy to do in a macro file. But I want to implement it in src file so that I can visualize the entire process. Kindly answer.
Hello, I'm sorry to don't answer your question, but I'm gonna give you a lead towards the current trend in the field. Have you ever heard of Backward Monte-Carlo applied to muon transport* ? Here you'll find a paper explaining the principle, which is basically sending the muons "from" the detector, instead of "to" the detector. It helps improving a lot the computer resources and statistics.
*Article Backward Monte-Carlo applied to muon transport
Hi, if I understood your question currently - you want to create a muon beam that doesn't have a point source but multiple point sources such that the generation area resemble a rectangle? - if that is the case, I suggest two things - look into using CRY (https://nuclear.llnl.gov/simulation/doc_cry_v1.7/cry.pdf) to generate the initial muon flux and then in the PrimaryGeneratorAction.cc file create a nested for loop (2 loops) to create an array of such point sources that would effectively create the rectangular source that you need.
Hints:
1) Define a particle Gun:
% particleGun = new G4ParticleGun();
2) Define a generation point and angle:
% particleGun->SetParticlePosition(G4ThreeVector(x, y, z));
% particleGun->SetParticleMomentumDirection(G4ThreeVector(( xu,yu,zu ))
3) Loop over x and y points to create the rectangle
- Badges
- Science topic
- Similar topics
- Filing