|Message: Re: Problems with Radioactive Decay: no decay lines||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
Thanks for you response.
In fact, that's not exactly what I meant.
I'm trying to simulate the decay of an isomer in 210Ra to radon (not the full chain) and get the response of the detectors to the particles from the decay.
I coded in the simplified version of the detector (only the silicon double-sided strip detector for the moment) and a gun that generates my decaying nuclei in the volume of the detector. I have tested my code with a simpler gun (shooting the alpha-particles of a given energy of 7016. keV) G4int nofParticles = 1; fParticleGun = new G4ParticleGun(nofParticles);
// Define particle properties G4String particleName = "alpha"; G4ThreeVector position(0, 0, -5.*um); G4ThreeVector momentum(0, 0, 0.*keV);
// Default particle kinematics G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable(); G4ParticleDefinition* particle = particleTable->FindParticle(particleName); fParticleGun->SetParticleDefinition(particle); fParticleGun->SetParticleMomentum(momentum); fParticleGun->SetParticlePosition(position); … //randomized position G4double x= 0.*cm, d_x= 20.*cm; G4double y= 0.*cm, d_y= 20.*cm; x += d_x*(G4UniformRand()-0.5)*2; y += d_y*(G4UniformRand()-0.5)*2; G4ThreeVector position(x, y, -2.*um); //implantation depth fParticleGun->SetParticlePosition(position);
// randomized direction G4double dtheta = 360.*deg; G4double dphi = 360*deg; G4double theta = G4UniformRand()*dtheta; G4double phi = G4UniformRand()*dphi; fParticleGun->SetParticleMomentumDirection(G4ThreeVector(sin(theta)*sin(phi), sin(theta)*cos(phi), cos(theta))); fParticleGun->SetParticleEnergy(7016.*keV); fParticleGun->GeneratePrimaryVertex(event);
The response of my "detector" is exactly what I expect it to be - a 7016 line and some not-full-energy-deposing particles: http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2014/06/16/08.03-28433-4-06-16_at_5.03.14_PM.png
But with the Radioactive decay gun (defined as in my first post) I do not get any full-energy alpha peak:
For sure, this has to be wrong, so I guess that I'm not implementing either the gun, or something else related to the radioactive decay correctly.
Could you please help me out with this one?
|Inline Depth:||Outline Depth:||Add message:|