|Message: antiproton annihilation at rest with CHIPS||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)|
I tried to simulate the antiproton annihilation at rest with GEANT4.9.0 and CHIPS model, but I was not able to get the same result as Mikhail Kossov's paper (CERN), "Simulation of Antiproton-Nuclear Annihilation at Rest" (2005), where GEANT4 and CHIPS were used.
Can anyone give me some advice for this problem? (Did he use the special code for that?)
Here is my physics list for antiproton.
---------- // anti-proton Physics pManager = G4AntiProton::AntiProton()->GetProcessManager();
// G4QCaptureAtRest G4QCaptureAtRest* captureProc = new G4QCaptureAtRest(); pManager->AddRestProcess(captureProc);
// add process G4HadronElasticProcess* theapElasticProcess = new G4HadronElasticProcess(); G4LElastic* theapElasticModel = new G4LElastic(); theapElasticProcess->RegisterMe(theapElasticModel); pManager->AddDiscreteProcess(theapElasticProcess);
G4AntiProtonInelasticProcess* theAntiProtonInelasticProcess = new G4AntiProtonInelasticProcess();
G4LEAntiProtonInelastic* theAntiProtonLEPModel = new G4LEAntiProtonInelastic(); G4HEAntiProtonInelastic* theAntiProtonHEPModel = new G4HEAntiProtonInelastic(); theAntiProtonInelasticProcess->RegisterMe(theAntiProtonLEPModel); theAntiProtonInelasticProcess->RegisterMe(theAntiProtonHEPModel); pManager->AddDiscreteProcess(theAntiProtonInelasticProcess);
G4AntiProtonAnnihilationAtRest* theAntiProtonAnnihilation = new G4AntiProtonAnnihilationAtRest(); pManager->AddRestProcess(theAntiProtonAnnihilation);
G4VProcess* theapMultipleScattering = new G4MultipleScattering(); G4VProcess* theapIonisation = new G4hIonisation();
// set ordering for AlongStepDoIt pManager->SetProcessOrdering(theapMultipleScattering, idxAlongStep,1); pManager->SetProcessOrdering(theapIonisation, idxAlongStep,2); // // set ordering for PostStepDoIt pManager->SetProcessOrdering(theapMultipleScattering, idxPostStep,1); pManager->SetProcessOrdering(theapIonisation, idxPostStep,2); ----------
Thanks in advance.
|Inline Depth:||Outline Depth:||Add message:|