Message: Re: Problems with Optical Photons [OpAbsorption] Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Problems with Optical Photons [OpAbsorption] 

Forum: Processes Involving Optical Photons
Re: Question Problems with Optical Photons [OpAbsorption] (Ana)
Re: Feedback Re: Problems with Optical Photons [OpAbsorption] (Gumplinger Peter)
Date: 14 Mar, 2011
From: Ana <Ana>

Hi Peter,

I had also tried with:

if (particleName == "opticalphoton") {
      pmanager->AddDiscreteProcess(new G4OpAbsorption());
      pmanager->AddDiscreteProcess(new G4OpRayleigh());
      pmanager->AddDiscreteProcess(theBoundaryProcess);
    }

and I have had also problems.

The ouput is:

*********************************************************************************************************
* G4Track Information:   Particle = opticalphoton,   Track ID = 3,   Parent ID = 1
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0     0.97        0        0  2.95e-06        0        0         0   caloBlock initStep
  G4ParticleChange::CheckIt  : the Momentum Change is not unit vector !!
  Difference:  1
 G4ParticleChange::CheckIt 
      -----------------------------------------------
        G4ParticleChange Information  
      -----------------------------------------------
        # of 2ndaries       :                    0
      -----------------------------------------------
        Energy Deposit (MeV):                    0
        Non-ionizing Energy Deposit (MeV):                    0
        Track Status        :                Alive
        True Path Length (mm) :                 18.1
        Stepping Control      :                    0
    Last Step In the voulme  : 
        Mass (GeV)   :                    0
        Charge (eplus)   :                    0
        MagneticMoment   :                    0
                :  =                    0*[e hbar]/[2 m]
        Position - x (mm)   :                -3.32
        Position - y (mm)   :                -14.5
        Position - z (mm)   :                  -10
        Time (ns)           :                 -nan
        Proper Time (ns)    :                 -nan
        Momentum Direct - x :                    0
        Momentum Direct - y :                    0
        Momentum Direct - z :                    0
        Kinetic Energy (MeV):             2.95e-06
        Polarization - x    :                    0
        Polarization - y    :                    0
        Polarization - z    :                    0

*** G4Exception : 200
      issued by : G4ParticleChange::CheckIt
momentum, energy, and/or time was illegal
*** Event Must Be Aborted 
    1    -3.32    -14.5      -10  2.95e-06        0     18.1      18.1     expHall Transportation
Track (trackID 3, parentID 1) is processed with stopping code 3
NULL returned from G4StackManager.
Terminate current event processing.
Run terminated.
Run Summary
  Number of events processed : 1
  User=0.5600000000000001s Real=0.57s Sys=0.01s
G4 kernel has come to Quit state.

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

And about your question, yes before I have worked with geant4.9.2.p02 and with my simulation without problems, and I have obtained good results. When I say good results I mean results with physical sense.

Now, I want to run my simulation in geant4.9.4.p01 and I had my first problem with the "LowEnergy" packages. I thought it should be easy to change, nowadays I have tried with the Penelope package and with the PhysicList that I have posted in my first message.

I have also made I new simulation much easier, and I have the same problem.

I attach this last simulation, if somebody have time to see it and tell me something about the problem I would be very grateful.

Thanks Ana

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2011/03/14/06.04-70486-01DetectorConstruction.cc
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2011/03/14/06.04-76576-ExN01PhysicsList.cc
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2011/03/14/06.04-19298-PrimaryGeneratorAction.cc
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2011/03/14/06.04-11144-01DetectorConstruction.hh
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2011/03/14/06.04-51650-ExN01PhysicsList.hh
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2011/03/14/06.04-42208-PrimaryGeneratorAction.hh
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2011/03/14/06.04-84947-exampleN01.cc
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2011/03/14/06.05-33819-exampleN01.cc

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

1 Feedback: Re: Problems with Optical Photons [OpAbsorption]   (Gumplinger Peter - 15 Mar, 2011)
2 Warning: Re: Problems with Optical Photons [OpAbsorption]   (Gumplinger Peter - 15 Mar, 2011)
1 Sad: Re: Problems with Optical Photons [OpAbsorption]   (Ana - 22 Mar, 2011)
 Add Message Add Message
to: "Re: Problems with Optical Photons [OpAbsorption]"

 Subscribe Subscribe

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