| Message: Re: QGSP_BERT_HP + precompound -> Nucleus A, Z = 0 0 Projectile was a Aborted | Not Logged In (login) |
|
Hello Vladimir,
Thank you very much for your reply. It is actually my own QGSP_INCL_ABLA_HP (not yet released physics list). Its code is based on QGSP_BERT_HP, with the appropriate modifications. Because of the energy range of validity of the INCL model (150 MeV - 3 GeV), I though the precompound model has to be included in addition to it. I included it simply like this:
theNeutrons=new G4NeutronBuilder; theNeutrons->RegisterMe(theQGSPNeutron=new G4QGSPNeutronBuilder(QuasiElastic, ProjectileDiffraction)); theNeutrons->RegisterMe(theLEPNeutron=new G4LEPNeutronBuilder); theLEPNeutron->SetMinEnergy(19.9*MeV); theLEPNeutron->SetMinInelasticEnergy(9.5*GeV); theLEPNeutron->SetMaxInelasticEnergy(25*GeV);
theNeutrons->RegisterMe(theBertiniNeutron=new G4BertiniNeutronBuilder); theBertiniNeutron->SetMinEnergy(2.9*GeV); theBertiniNeutron->SetMaxEnergy(9.9*GeV);
theNeutrons->RegisterMe(theInclAblaNeutron=new G4InclAblaNeutronBuilder); theInclAblaNeutron->SetMinEnergy(0.15*GeV); theInclAblaNeutron->SetMaxEnergy(3.0*GeV);
theNeutrons->RegisterMe(thePrecoNeutron=new G4PrecoNeutronBuilder); thePrecoNeutron->SetMinEnergy(19.9*MeV);
theNeutrons->RegisterMe(theHPNeutron=new G4NeutronHPBuilder); theHPNeutron->SetMinEnergy(0.0*MeV); I assumed that the precompound model is not part of the inclabla model. Please let me know if it is and I will remove it from the above code. Thank you, Cristian
|
| Inline Depth: | Outline Depth: | Add message: |
|
to: |