Message: Understanding ProcName Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Understanding ProcName 

Forum: Processes Involving Optical Photons
Date: 21 Aug, 2013
From: <bg92>

Dear all,

I have some hard time to figure out the meaning of ProcName is some of my events. In the G4Track information, I get at a particular step of a pi-:

...
119    -8.52      351 9.01e+03         0        0        0  3.41e+03 physical_box_TEST Scintillation
    :----- List of 2ndaries - #SpawnInStep=  2(Rest= 2,Along= 0,Post= 0), #SpawnTotal= 89 ---------------
    :     -8.52       351  9.01e+03      29.8         anti_nu_mu
    :     -8.52       351  9.01e+03      4.12                mu-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------

Why this process is called scintillation?

The muon decay gives :

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

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0    -8.52      351 9.01e+03      4.12        0        0         0 physical_box_TEST initStep
    1    -8.58      350 9.01e+03         0     4.12    0.794     0.794 physical_box_TEST muIoni
    2    -8.58      350 9.01e+03         0        0        0     0.794 physical_box_TEST Scintillation
    :----- List of 2ndaries - #SpawnInStep=  3(Rest= 3,Along= 0,Post= 0), #SpawnTotal=  3 ---------------
    :     -8.58       350  9.01e+03      45.4              nu_mu
    :     -8.58       350  9.01e+03      11.4          anti_nu_e
    :     -8.58       350  9.01e+03      47.8                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------

Still called scintillation.

In the same idea, looking at Cerenkov photons, I wanted to find them using ->GetProcessName() == "Cerenkov", but it seems that Cerenkov photon are not all produced under the ProcName Cerenkov. For instance, the following step of a Kaon (at 5GeV, so a lot of Cerenkov photons):

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
43    0.499      360 5.34e+03  4.53e+03   0.0141     93.4  5.34e+03 cad_physical_bar Transportation
44    0.501      360 5.35e+03  4.53e+03     2.39     7.22  5.35e+03 cad_physical_bar Cerenkov
    :----- List of 2ndaries - #SpawnInStep=296(Rest= 0,Along= 0,Post=296), #SpawnTotal=338 ---------------
    :       0.5       360  5.34e+03     3e-06      opticalphoton
    :       0.5       360  5.34e+03   2.9e-06      opticalphoton
    :       0.5       360  5.34e+03  2.53e-06      opticalphoton
    :     0.501       360  5.35e+03   2.4e-06      opticalphoton
    :       0.5       360  5.35e+03  3.34e-06      opticalphoton
    .... 
My Optical photons are under the ProcName Cerenkov, so this is ok.

46    0.498      360 5.36e+03  4.52e+03    0.581     2.07  5.36e+03 cad_physical_bar hIoni
    :----- List of 2ndaries - #SpawnInStep= 94(Rest= 0,Along= 0,Post=94), #SpawnTotal= 94 ---------------
    :     0.498       360  5.36e+03      5.41                 e-
    :     0.499       360  5.36e+03  2.54e-06      opticalphoton
    :     0.498       360  5.36e+03  2.83e-06      opticalphoton
    :     0.498       360  5.36e+03  2.76e-06      opticalphoton
    :     0.498       360  5.36e+03  3.13e-06      opticalphoton
    :     0.499       360  5.36e+03  2.52e-06      opticalphoton
    :     0.499       360  5.36e+03   2.5e-06      opticalphoton
    :     0.498       360  5.36e+03  2.08e-06      opticalphoton
    ... 
Here, the Kaon step is under the name hIoni, but I have also optical photon that are Cerenkov photons.

47    0.498      360 5.36e+03  4.52e+03    0.251    0.743  5.36e+03 world_physical Transportation
    :----- List of 2ndaries - #SpawnInStep= 40(Rest= 0,Along= 0,Post=40), #SpawnTotal= 40 ---------------
    :     0.498       360  5.36e+03  2.06e-06      opticalphoton
    :     0.498       360  5.36e+03  4.03e-06      opticalphoton
    :     0.498       360  5.36e+03  2.12e-06      opticalphoton
    :     0.498       360  5.36e+03  2.72e-06      opticalphoton
    :     0.498       360  5.36e+03  3.62e-06      opticalphoton
    ...
And here, the Kaon step is under the name Transportation, where I also have Cerenkov photon.

So it seems to me that the Cerenkov photon can hide under a different process name than Cerenkov, so not a good idea to use ->GetProcessName() == "Cerenkov".

Do you have a clarification to propose, because I am a little bit confused by the ProcName variable.

Thank you very much for any help you could give me.

Best Baptiste

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

1 Feedback: Re: Understanding ProcName   (Gumplinger Peter - 23 Aug, 2013)
 Add Message Add Message
to: "Understanding ProcName"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews