|Message: Introducing randomization of resonance mass causes problem with decay||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)|
Greetings to All,
I would like to share with you a problem, which I hope can be resolved with your help.
In short, I'm investigating with Geant4 (version 9.6.p02) the detector response to products of decay of the baryonic resonances e.g. Delta(1232), N(1440). I have noticed, that Geant4 doesn't "smear" the mass of the resonances (though it contains all the information such as the width). So, I decided to randomize the mass of primary resonance following Breit-Wigner distribution by setting it through SetMass() method in G4PrimaryParticle class. I did it because I believed (and still believe) that Geant4 uses this information to properly set the kinematics of the decay products, what is crucial in my problem. However, if the mass of resonance is too low the exception is thrown:
G4PhaseSpaceDecayChannel::TwoBodyDecayIt sum of daughter mass is larger than parent mass parent :N(1440)+ 1.02138 daughter 1 :neutron 0.939565 daughter 2:pi+ 0.13957
-------- WWWW ------- G4Exception-START -------- WWWW ------- *** G4Exception : PART112 issued by : G4PhaseSpaceDecayChannel::TwoBodyDecayIt Can not create decay products: sum of daughter mass is larger than parent mass *** This is just a warning message. *** -------- WWWW -------- G4Exception-END --------- WWWW -------
Why Geant4 doesn't check which decay channel is allowed from energy and momentum conservation? How should I proceed to have properly set the kinematics of the products of decay?
I have found, that in Geant4 10.1 some solution to the similiar problem was found (from release notes): "Added new constructor with dynamical mass in G4DynamicParticle. G4PhaseSpaceDecayChannel takes in to account mass width of daughter particles and gives the dynamical mass according to Breit-Wigner distribution. Added daughters_width in G4VDecayChannel. Addressing problem report #1629", but I'm not sure if switching to newer Geant4 version is really the final solution.
Best regards, Rafal