Message: antiproton annihilation at rest with CHIPS Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question antiproton annihilation at rest with CHIPS 

Keywords: antiproton, CHIPS
Forum: Hadronic Processes
Date: 13 Mar, 2008
From: Tsuguo <tsuguo@phys.columbia.edu>

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();

   pManager->AddProcess(theapIonisation);
   pManager->AddProcess(theapMultipleScattering);

   // 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.

Best regards,

Tsuguo

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

1 None: Re: antiproton annihilation at rest with CHIPS   (Vladimir IVANTCHENKO - 14 Mar, 2008)
1 None: Re: antiproton annihilation at rest with CHIPS   (Tsuguo - 14 Mar, 2008)
1 None: Re: antiproton annihilation at rest with CHIPS   (Tsuguo - 14 Mar, 2008)
...
3 None: Re: antiproton annihilation at rest with CHIPS   (Tsuguo - 17 Mar, 2008)
1 None: Re: antiproton annihilation at rest with CHIPS   (Vladimir IVANTCHENKO - 18 Mar, 2008)
... 7 Message(s)
2 Question: Re: antiproton annihilation at rest with CHIPS   (Dr_Cards - 08 Feb, 2010)
 Add Message Add Message
to: "antiproton annihilation at rest with CHIPS"

 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 ]