Message: proton QMD model -> bad energy non-conservation Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question proton QMD model -> bad energy non-conservation 

Forum: Hadronic Processes
Date: 05 Jan, 2016
From: Cristian Bungau <Cristian Bungau>

Dear all,

I am trying to use the proton QMD model:

	  // Inelastic scattering
	  G4ProtonInelasticProcess* theInelasticProcess = 
	    new G4ProtonInelasticProcess("inelastic");
         theInelasticProcess->RegisterMe(bcModel);
         theInelasticProcess->RegisterMe(qmd);
         theInelasticProcess->AddDataSet(shenXS);
         theInelasticProcess->AddDataSet(tripXS);
         theInelasticProcess->AddDataSet(tripLightXS);

where

  G4BinaryCascade* bcModel = new G4BinaryCascade();
  bcModel->SetMinEnergy(0.*MeV);
  bcModel->SetMaxEnergy(101.*MeV);

  // QMD model
  G4QMDReaction* qmd = new G4QMDReaction;
  qmd->SetMinEnergy(100.*MeV);
  qmd->SetMaxEnergy(10.*GeV);

  G4IonsShenCrossSection* shenXS = new G4IonsShenCrossSection;
  G4TripathiCrossSection* tripXS = new G4TripathiCrossSection;
  G4TripathiLightCrossSection* tripLightXS = new G4TripathiLightCrossSection;

I am guessing the coss sections data sets are not appropriate for protons (not really light ions)? But what cross-section data set should I be using in order to use the QMD model for protons?

The warning messages I get are:

---> Begin of event: 0

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : had012
      issued by : G4HadronicProcess:CheckResult()
Warning: Bad energy non-conservation detected, will re-sample the interaction
 Process / Model: inelastic / QMDModel
 Primary: proton (2212),  E= 1048.44, target nucleus (29, 65)
 E(initial - final) = -1027.03 MeV.

*** This is just a warning message. *** -------- WWWW -------- G4Exception-END --------- WWWW -------

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : had012
      issued by : G4HadronicProcess:CheckResult()
Warning: Bad energy non-conservation detected, will re-sample the interaction
 Process / Model: inelastic / QMDModel
 Primary: proton (2212),  E= 1126.81, target nucleus (26, 56)
 E(initial - final) = -1116.22 MeV.

*** This is just a warning message. *** -------- WWWW -------- G4Exception-END --------- WWWW -------

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : had012
      issued by : G4HadronicProcess:CheckResult()
Warning: Bad energy non-conservation detected, will re-sample the interaction
 Process / Model: inelastic / QMDModel
 Primary: proton (2212),  E= 1052.28, target nucleus (26, 56)
 E(initial - final) = -1040.64 MeV.

*** This is just a warning message. *** -------- WWWW -------- G4Exception-END --------- WWWW -------

---> Begin of event: 10

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : had012
      issued by : G4HadronicProcess:CheckResult()
Warning: Bad energy non-conservation detected, will re-sample the interaction
 Process / Model: inelastic / QMDModel
 Primary: proton (2212),  E= 1137.47, target nucleus (26, 56)
 E(initial - final) = -1134.35 MeV.

*** This is just a warning message. *** -------- WWWW -------- G4Exception-END --------- WWWW -------

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : had012
      issued by : G4HadronicProcess:CheckResult()
Warning: Bad energy non-conservation detected, will re-sample the interaction
 Process / Model: inelastic / QMDModel
 Primary: proton (2212),  E= 1079.62, target nucleus (26, 56)
 E(initial - final) = -1088.18 MeV.

*** This is just a warning message. *** -------- WWWW -------- G4Exception-END --------- WWWW -------

---> Begin of event: 20

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : had012
      issued by : G4HadronicProcess:CheckResult()
Warning: Bad energy non-conservation detected, will re-sample the interaction
 Process / Model: inelastic / QMDModel
 Primary: proton (2212),  E= 1128.33, target nucleus (26, 56)
 E(initial - final) = -1113.01 MeV.

*** This is just a warning message. *** -------- WWWW -------- G4Exception-END --------- WWWW -------

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : had012
      issued by : G4HadronicProcess:CheckResult()
Warning: Bad energy non-conservation detected, will re-sample the interaction
 Process / Model: inelastic / QMDModel
 Primary: proton (2212),  E= 1040.14, target nucleus (26, 56)
 E(initial - final) = -1042.64 MeV.

*** This is just a warning message. *** -------- WWWW -------- G4Exception-END --------- WWWW -------

---> Begin of event: 30

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : had012
      issued by : G4HadronicProcess:CheckResult()
Warning: Bad energy non-conservation detected, will re-sample the interaction
 Process / Model: inelastic / QMDModel
 Primary: proton (2212),  E= 1079.7, target nucleus (29, 63)
 E(initial - final) = -1060.47 MeV.

*** This is just a warning message. *** -------- WWWW -------- G4Exception-END --------- WWWW -------

Thank you,

Cristian

 Add Message Add Message
to: "proton QMD model -> bad energy non-conservation"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]