Message: Problem in Muon Capture(No GenericIon) Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Problem in Muon Capture(No GenericIon) 

Forum: Particles
Date: 12 Nov, 2010
From: Yoshi <Yoshi>

Hi,I would like to simulate muon capture mu- beam and Pb Target.
So,I added G4MuMinusCaptureAtRest to ExN02PhysicsList like following.
     .
     .
     .
else if( particleName == "mu+" || 
               particleName == "mu-"    ) {
      //muon  
      pmanager->AddProcess(new G4MuMultipleScattering, -1, 1, 1);
      pmanager->AddProcess(new G4MuIonisation,         -1, 2, 2);
      pmanager->AddProcess(new G4MuBremsstrahlung,     -1, 3, 3);
      pmanager->AddProcess(new G4MuPairProduction,     -1, 4, 4);       
      if ( particleName == "mu-"){
        pmanager->AddProcess(new G4MuonMinusCaptureAtRest(),0,-1,-1);
        }
     .
     .
     .

Then compile error is nothing,but I run ExN02,problem is occurred. (By the way, I changed ExN02 little bit) ====================================================================

### Run 0 start.

*********************************************************************************************************
* G4Track Information:   Particle = mu-,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#     KineE    dEStep   StepLeng  TrakLeng    Volume     ParentID       Particle       Process
    0     40 MeV     0 eV      0 fm      0 fm       World        0              mu-       initStep
    1   39.9 MeV  66.4 keV    26 cm     26 cm       World        0              mu-   Transportation
    2   34.1 MeV  5.81 MeV  1.16 cm   27.2 cm      Sinti1        0              mu-           muIoni
    3   29.2 MeV  4.96 MeV  8.37 mm     28 cm      Sinti1        0              mu-   Transportation
    4   29.1 MeV   110 keV  32.8 cm   60.8 cm       World        0              mu-   Transportation
    5   24.2 MeV  4.86 MeV  7.49 mm   61.6 cm     Sinti2a        0              mu-           muIoni
    6     20 MeV   4.2 MeV  5.89 mm   62.1 cm     Sinti2a        0              mu-           muIoni
    7   19.5 MeV   500 keV   639 um   62.2 cm     Sinti2a        0              mu-   Transportation
    8   19.5 MeV  3.29 keV  8.03 mm     63 cm       World        0              mu-   Transportation
    9   12.8 MeV  6.66 MeV  2.27 mm   63.2 cm     Target1        0              mu-           muIoni
   10   5.32 MeV  7.51 MeV  1.65 mm   63.4 cm     Target1        0              mu-           muIoni
   11      0 eV   5.32 MeV   479 um   63.5 cm     Target1        0              mu-           muIoni

*** G4Exception : No GenericIon
      issued by : G4VUserPhysicsList::AddProcessManager
GenericIon does not exist
*** Run Must Be Aborted 
   12      0 eV      0 eV      0 fm   63.5 cm     Target1        0              mu-            Decay	
    :----- List of 2ndaries - #SpawnInStep= 17(Rest=17,Along= 0,Post= 0), #SpawnTotal= 17 ---------------
    :   3.73 mm  -6.44 mm   13.4 cm    413 keV   neutron
    :   3.73 mm  -6.44 mm   13.4 cm    855 keV     gamma
    :   3.73 mm  -6.44 mm   13.4 cm    911 keV     gamma
    :   3.73 mm  -6.44 mm   13.4 cm    860 keV     gamma
    :   3.73 mm  -6.44 mm   13.4 cm    704 keV     gamma
    :   3.73 mm  -6.44 mm   13.4 cm    216 keV     gamma
    :   3.73 mm  -6.44 mm   13.4 cm    453 keV     gamma
    :   3.73 mm  -6.44 mm   13.4 cm    686 keV     gamma
    :   3.73 mm  -6.44 mm   13.4 cm    266 keV     gamma
    :   3.73 mm  -6.44 mm   13.4 cm   13.4 keV Tl206[0.0]
    :   3.73 mm  -6.44 mm   13.4 cm   81.1 MeV     nu_mu
    :   3.73 mm  -6.44 mm   13.4 cm   96.4 keV        e-
    :   3.73 mm  -6.44 mm   13.4 cm    659 keV     gamma
    :   3.73 mm  -6.44 mm   13.4 cm   9.64 MeV     gamma
    :   3.73 mm  -6.44 mm   13.4 cm   48.3 MeV     nu_mu
    :   3.73 mm  -6.44 mm   13.4 cm   41.8 MeV anti_nu_e
    :   3.73 mm  -6.44 mm   13.4 cm   14.5 MeV        e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
>>> Event 0
    0 trajectories stored in this event.
Idle> exit
Graphics systems deleted.
Visualization Manager deleting...

*** G4Exception : No GenericIon
      issued by : G4VUserPhysicsList::AddProcessManager
GenericIon does not exist
*** Run Must Be Aborted 

Is it not enough only to add G4MuMinusCaptureAtRest? Or, is it not good for simulation of muon capture to use ExN02?? Thanks,

  yoccha

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

1 Idea: Re: Problem in Muon Capture(No GenericIon)   (Vladimir Ivanchenko - 13 Nov, 2010)
 Add Message Add Message
to: "Problem in Muon Capture(No GenericIon)"

 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 ]