Message: Re: Antiproton annihilation cross section Not Logged In (login)

Question Re: Antiproton annihilation cross section 

Forum: Hadronic Processes
Re: Question Antiproton annihilation cross section (Kirill Skovpen)
Re: Idea Re: Antiproton annihilation cross section (Vladimir Ivanchenko)
Date: 08 Apr, 2009
From: Kirill Skovpen <Kirill Skovpen>

Dear, Vladimir,

Thank you very much for answering. Well, I was not using default PhysicsLists, I've made up my own. As for antiprotons, defined processes are:

        else if (particleName == "anti_proton")
          {
             pmanager->AddDiscreteProcess(theElasticProcess);

             G4AntiProtonInelasticProcess* theInelasticProcess =
               new G4AntiProtonInelasticProcess();

             G4LEAntiProtonInelastic* theLEInelasticModel =
                    new G4LEAntiProtonInelastic;
             theLEInelasticModel->SetMaxEnergy(20.0*GeV);
             theInelasticProcess->RegisterMe(theLEInelasticModel);

             pmanager->AddDiscreteProcess(theInelasticProcess);

             pmanager->AddRestProcess(new G4AntiProtonAnnihilationAtRest);

             pmanager->AddRestProcess(new G4QCaptureAtRest());
          }

and, of course, EM part:

        pmanager->AddProcess(new G4MultipleScattering(), -1, 1, 1);
        pmanager->AddProcess(new G4hIonisation(), -1, 2, 2);

Or maybe you are suggesting to use another models/processes for antiprotons?

Thank you,

 Kirill Skovpen

Inline Depth:
 0 0
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

Idea Re: Antiproton annihilation cross section  by Vladimir Ivanchenko <Vladimir Ivanchenko>,   12 Apr, 2009
Re: Question Re: Antiproton annihilation cross section (Kirill Skovpen)
Hello Kirill,

Your Physics List is similar to the default which is using parameterized model (GHEISHA, derived from G3). There is a bug - you instantiate 2 AtRest processes - one obsolere G4AntiProtonAnnihilationAtRest and recommanded one G4QCaptureAtRest. I would remove the first one.

To have better anti_proton physics you likely need to use other G4 models. I would let you communicate with original authors of FTF base models, QGS base models, CHIPS models and Bertini. They should clarify if any of above can be used now and at what energy.

VI

1 None: Re: Antiproton annihilation cross section   (Kirill Skovpen - 16 Apr, 2009)
 Add Message Add Message
to: "Re: Antiproton annihilation cross section"

 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 ]