Message: MuonMinusCaptureAtRest and decay products Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None MuonMinusCaptureAtRest and decay products 

Forum: Event and Track Management
Date: 11 Apr, 2013
From: Matthew T <Matthew T>

Dear experts,

I am simulating mu- in doped water. As expected, sometimes the muon decays and sometimes it is captured, however, when I look at the daughter particles at the time of tracking, I see too many neutrinos and e-. For example, for a single mu- decay, I get the following verbose output (where I have suppressed e- and gamma tracking for simplicity). It appears that the decay produces the following particles: numu created by muMinusCaptureAtRest, nuebar created by muMinusCaptureAtRest, e- created by muMinusCaptureAtRest, numu created by Decay, nuebar created by Decay, and e- created by Decay. In this case, the "muMinusCaptureAtRest" particles are created at global time 460 ns and the "Decay" particles are created at global time 3092 ns. In addition, for events where the mu- is captured on a nucleus, extra muon decay products are produced, in addition to the normal beta decay products. Is this caused by a bug in the invocation of certain processes in my physics list (pManager->AddProcess(new G4MuonMinusCaptureAtRest(), 0,-1,-1);)? Does G4MuonMinusCaptureAtRest actually kill the muon after capture like it should, or does the muon capture and then stick around until it decays? Or do the extra particles eventually get killed at some time other than tracking time?

Thank you, Matt

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

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0   -1e+03        0       111        0        0         0  World_phys initStep
    1  -0.0205     -966  -0.0197       111  0.00618     34.2      34.2  World_phys muIoni
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  1 ---------------
    :   -0.0205      -966   -0.0197   0.00183                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
    2   -0.229     -853    -0.15       111   0.0163      113       147  World_phys muIoni
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  2 ---------------
    :    -0.229      -853     -0.15   0.00837                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
    3   -0.262     -836   -0.164       111  0.00168     16.9       164  World_phys muIoni
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  3 ---------------
    :    -0.262      -836    -0.164   0.00261                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
    4   -0.336     -798   -0.196       111  0.00612     37.8       202  World_phys muIoni
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  4 ---------------
    :    -0.336      -798    -0.196    0.0014                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
    5    -0.59     -546    -1.04       111   0.0364      252       454  World_phys muIoni
    6   -0.682     -513    -1.01       111  0.00531       33       487  World_phys muIoni
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  5 ---------------
    :    -0.682      -513     -1.01   0.00423                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
    7   -0.815     -430    -1.01       111   0.0119     82.7       570  World_phys muIoni
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  6 ---------------
    :    -0.815      -430     -1.01   0.00634                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
    8    -1.09     -254   -0.802       111   0.0239      176       746 graphite_phys Transportation
    9   -0.786     -226   -0.962      98.3       12       28       774 graphite_phys muIoni
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  7 ---------------
    :    -0.786      -226    -0.962     0.569                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
   10    0.185     -198    0.324      86.4     11.3     27.8       802 graphite_phys muIoni
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  8 ---------------
    :     0.185      -198     0.324     0.621                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
   11    0.932     -187     1.04      80.6     5.05     11.7       814 graphite_phys muIoni
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  9 ---------------
    :     0.932      -187      1.04     0.665                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
   12     2.38     -173     1.99      73.6     5.85     13.6       827 graphite_phys muIoni
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 10 ---------------
    :      2.38      -173      1.99       1.2                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
   13     2.42     -173     2.01      72.4    0.337    0.454       828 graphite_phys muIoni
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 11 ---------------
    :      2.42      -173      2.01     0.797                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
   14     4.45     -151     2.72        62     10.5     21.8       849 graphite_phys muIoni
   15     6.62     -134     2.84      53.6     8.36     17.4       867 graphite_phys muIoni
   16     6.63     -134     2.84      52.8   0.0509    0.111       867 graphite_phys muIoni
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 12 ---------------
    :      6.63      -134      2.84     0.786                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
   17     7.84     -122     3.18      45.8     6.01     11.4       878 graphite_phys muIoni
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 13 ---------------
    :      7.84      -122      3.18     0.995                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
   18     8.55     -118     3.29      42.4     2.67     4.27       883 graphite_phys muIoni
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 14 ---------------
    :      8.55      -118      3.29     0.702                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
   19       10     -108     3.19      35.9     6.52     10.2       893 graphite_phys muIoni
   20     11.3     -101     2.92        31     4.87     7.09       900  Tyvek_phys Transportation
   21     12.2    -95.2     2.73      28.9     2.12     5.82       906  water_phys Transportation
   22     13.4    -86.4     2.36      24.7     4.19     8.93       915  water_phys muIoni
   23     14.4    -79.3     2.22      20.4     4.28     7.19       922  water_phys muIoni
   24     15.7    -73.9     1.62      16.5     3.87     5.61       927  water_phys muIoni
   25     16.6    -69.6     1.83        13     3.53     4.36       932  water_phys muIoni
   26     17.5    -66.4     2.22      9.85     3.16     3.36       935  water_phys muIoni
   27     18.2    -63.9     2.13      6.89     2.96     2.59       938  water_phys muIoni
   28     18.6    -62.1     1.99      3.96     2.93     1.93       940  water_phys muIoni
   29     18.8      -61     1.82    0.0374     3.92      1.1       941  water_phys muIoni
   30     18.8      -61     1.82         0   0.0374 0.000843       941  water_phys muIoni
   31     18.8      -61     1.82         0        0        0       941  water_phys Scintillation
    :----- List of 2ndaries - #SpawnInStep= 17(Rest=17,Along= 0,Post= 0), #SpawnTotal= 31 ---------------
    :      18.8       -61      1.82  0.000912                 e-
    :      18.8       -61      1.82  0.000329              gamma
    :      18.8       -61      1.82  0.000236                 e-
    :      18.8       -61      1.82   0.00031                 e-
    :      18.8       -61      1.82  0.000419              gamma
    :      18.8       -61      1.82  0.000586                 e-
    :      18.8       -61      1.82  0.000855                 e-
    :      18.8       -61      1.82   0.00132                 e-
    :      18.8       -61      1.82   0.00218              gamma
    :      18.8       -61      1.82   0.00402                 e-
    :      18.8       -61      1.82     0.162              gamma
    :      18.8       -61      1.82      33.7                 e-
    :      18.8       -61      1.82      17.7          anti_nu_e
    :      18.8       -61      1.82      53.6              nu_mu
    :      18.8       -61      1.82      39.6              nu_mu
    :      18.8       -61      1.82      40.9          anti_nu_e
    :      18.8       -61      1.82      24.1                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
mu-  Scintillation  in  water_phys  fStopAndKill
(18.84014325147097,-60.98596336319885,1.821316318296512)

*********************************************************************************************************
* G4Track Information:   Particle = anti_nu_e,   Track ID = 31,   Parent ID = 1
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0     18.8      -61     1.82      40.9        0        0         0  water_phys initStep
anti_nu_e	31	Decay	3092.364507341443	1	40.94581372139056
    1     60.7    -67.5      -96      40.9        0      107       107  Tyvek_phys Transportation
    2     63.1    -67.9     -102      40.9        0      6.1       113 graphite_phys Transportation
    3      128    -78.1     -254      40.9        0      166       279  World_phys Transportation
    4  4.3e+03     -731   -1e+04      40.9        0 1.06e+04  1.09e+04  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = nu_mu,   Track ID = 30,   Parent ID = 1
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0     18.8      -61     1.82      39.6        0        0         0  water_phys initStep
nu_mu	30	Decay	3092.364507341443	1	39.58334677746757
    1    -5.34    -95.9       54      39.6        0     67.3      67.3  Tyvek_phys Transportation
    2    -9.05     -101       62      39.6        0     10.3      77.6 graphite_phys Transportation
    3      -98     -229      254      39.6        0      247       325  World_phys Transportation
    4 -4.61e+03 -6.74e+03    1e+04      39.6        0 1.26e+04  1.29e+04  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = nu_mu,   Track ID = 29,   Parent ID = 1
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0     18.8      -61     1.82      53.6        0        0         0  water_phys initStep
nu_mu	29	muMinusCaptureAtRest	460.6700267250922	1	53.56386169452674
    1    -92.1      -27     84.1      53.6        0      142       142  Tyvek_phys Transportation
    2    -98.4    -25.1     88.8      53.6        0     8.12       150 graphite_phys Transportation
    3     -253     22.2      203      53.6        0      198       348  World_phys Transportation
    4   -1e+04    3e+03 7.43e+03      53.6        0 1.25e+04  1.28e+04  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = anti_nu_e,   Track ID = 28,   Parent ID = 1
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0     18.8      -61     1.82      17.7        0        0         0  water_phys initStep
anti_nu_e	28	muMinusCaptureAtRest	460.6700267250922	1	17.72468457961821
    1     73.5    -61.7    -28.9      17.7        0     62.7      62.7  Tyvek_phys Transportation
    2     80.6    -61.8    -32.9      17.7        0     8.15      70.9 graphite_phys Transportation
    3      246      -64     -126      17.7        0      189       260  World_phys Transportation
    4    1e+04     -192 -5.6e+03      17.7        0 1.12e+04  1.14e+04  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = gamma,   Track ID = 26,   Parent ID = 1
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0     18.8      -61     1.82     0.162        0        0         0  water_phys initStep
    1     18.7    -61.1      1.6     0.127        0    0.293     0.293  water_phys compt
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  1 ---------------
    :      18.7     -61.1       1.6    0.0346                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
    2    -40.2    -54.9     35.1     0.124        0       68      68.3  water_phys compt
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  2 ---------------
    :     -40.2     -54.9      35.1   0.00332                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
    3    -69.5    -66.2     67.5     0.124        0     45.2       113  Tyvek_phys Transportation
    4    -75.1    -68.4     73.8     0.124        0     8.67       122 graphite_phys Transportation
    5    -85.1    -72.3     84.8    0.0899        0     15.4       138 graphite_phys compt
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  3 ---------------
    :     -85.1     -72.3      84.8     0.034                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
    6      -65    -90.8     73.5    0.0682        0     29.6       167 graphite_phys compt
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  4 ---------------
    :       -65     -90.8      73.5    0.0218                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
    7     -106    -78.8      128    0.0555        0     69.2       236 graphite_phys compt
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  5 ---------------
    :      -106     -78.8       128    0.0127                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
    8     -106    -80.9      125    0.0543        0     3.39       240 graphite_phys compt
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  6 ---------------
    :      -106     -80.9       125   0.00116                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
    9    -91.1    -92.3      106     0.047        0     26.8       267 graphite_phys compt
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  7 ---------------
    :     -91.1     -92.3       106   0.00731                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
   10    -90.8    -91.5      107    0.0414        0    0.923       267 graphite_phys compt
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  8 ---------------
    :     -90.8     -91.5       107   0.00564                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
   11     -119    -85.1     54.3      0.04        0     59.9       327 graphite_phys compt
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  9 ---------------
    :      -119     -85.1      54.3   0.00143                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
   12     -131    -91.6       51    0.0347        0       14       341 graphite_phys compt
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 10 ---------------
    :      -131     -91.6        51   0.00529                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
   13     -128    -88.5     52.9    0.0308        0     4.76       346 graphite_phys compt
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 11 ---------------
    :      -128     -88.5      52.9   0.00384                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
   14     -159      -94     44.8    0.0298        0     31.8       378 graphite_phys compt
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 12 ---------------
    :      -159       -94      44.8   0.00107                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
   15     -159    -95.1       43    0.0297        0     2.14       380 graphite_phys compt
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 13 ---------------
    :      -159     -95.1        43  3.18e-05                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
   16     -159    -95.6     41.8    0.0289        0     1.34       381 graphite_phys compt
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 14 ---------------
    :      -159     -95.6      41.8  0.000877                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
   17     -153     -117     38.4    0.0286        0     22.4       404 graphite_phys compt
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 15 ---------------
    :      -153      -117      38.4  0.000217                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
   18     -152     -120     39.6    0.0286        0     3.59       407 graphite_phys compt
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 16 ---------------
    :      -152      -120      39.6  2.78e-05                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------
   19     -152     -121     39.9         0 0.000288     1.05       408 graphite_phys phot
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 17 ---------------
    :      -152      -121      39.9    0.0283                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------

*********************************************************************************************************
* G4Track Information:   Particle = gamma,   Track ID = 24,   Parent ID = 1
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0     18.8      -61     1.82   0.00218        0        0         0  water_phys initStep
    1     18.8      -61     1.81         0 0.000277   0.0124    0.0124  water_phys phot
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  1 ---------------
    :      18.8       -61      1.81   0.00191                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------

*********************************************************************************************************
* G4Track Information:   Particle = gamma,   Track ID = 20,   Parent ID = 1
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0     18.8      -61     1.82  0.000419        0        0         0  water_phys initStep
    1     18.8      -61     1.82         0 0.000353 0.000935  0.000935  water_phys phot
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  1 ---------------
    :      18.8       -61      1.82  6.62e-05                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------

*********************************************************************************************************
* G4Track Information:   Particle = gamma,   Track ID = 17,   Parent ID = 1
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0     18.8      -61     1.82  0.000329        0        0         0  water_phys initStep
    1     18.8      -61     1.82         0 0.000277 0.000385  0.000385  water_phys phot
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  1 ---------------
    :      18.8       -61      1.82  5.22e-05                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------

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

1 Feedback: Re: MuonMinusCaptureAtRest and decay products   (Gumplinger Peter - 26 Apr, 2013)
 Add Message Add Message
to: "MuonMinusCaptureAtRest and decay products"

 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 ]