Bugfix for HA2018 and HA2025 Pion Absorption#503
Conversation
…state energy budget for the multinucleon absorption mechanism
nusense
left a comment
nusense
left a comment
There was a problem hiding this comment.
No review here about the physics, just a couple of comments about the code implementation.
| double probKE = p->P4()->E() -probM; | ||
| double clusKE = probKE * (1./5.); | ||
| double probEAvail; | ||
| if ( pdgc==kPdgPiP || pdgc==kPdgPi0 || pdgc==kPdgPiM) probEAvail = p->P4()->E() + BE_correction_per_nucleon*5; //pion probe -- includes mass, BE correction for 4 nucleons |
There was a problem hiding this comment.
You might consider the use of genie::pdg::IsPion(pgdc) from Framework/ParticleData/PDGUtils.h here and elsewhere.
There was a problem hiding this comment.
Good call -- I'll update that
| TLorentzVector clusP4(pP3,clusKE); //no mass | ||
| double BE_correction_per_nucleon = .005; //GeV | ||
| probM -= BE_correction_per_nucleon*5; // BE correction | ||
| TVector3 pP3 = p->P4()->Vect() * (1./5.); |
There was a problem hiding this comment.
The indent for this line seems borked. Potentially due to use of tabs perhaps.
| if ( pdgc==kPdgPiP || pdgc==kPdgPi0 || pdgc==kPdgPiM) probEAvail = p->P4()->E() + probBE; //pion probe -- includes mass, BE correction for 4 nucleons | ||
| else probEAvail = p->P4()->E() - (probM - probBE); //proton or neutron probe -- does not include mass | ||
| double clusKE = probEAvail; // + np*0.9383 + nn*.9396; | ||
| TLorentzVector clusP4(pP3,clusKE); //no mass is correct |
There was a problem hiding this comment.
Another weird indent. Also you might want to re-word "no mass is correct" -- I assume it means "using no mass here is correct"
There was a problem hiding this comment.
I'll go in and fix the indents. I can update the comment, but I will note that for my bugfix, I did not change the line with the comment as written (I think I only changed the indentation on that line by accident).
2 participants
Bugfix: For pion absorption in HA2018 and HA2025, I add the pion mass to the final state energy budget for the multinucleon absorption mechanism.