Message: Re: LivermoreModel Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: LivermoreModel 

Forum: Physics List
Re: Question LivermoreModel
Date: 27 Jul, 2011
From: Paul Nicholas Colin Gloster <Paul Nicholas Colin Gloster>

Someone <dis@HLL.MPG.De> sent on July 26th, 2011:
|---------------------------------------------------------------------------------------------------------|
|"[..]                                                                                                    |
|                                                                                                         |
|Contained in my PhysicsList:                                                                             |
|                                                                                                         |
|ExPhysicsList::ExPhysicsList() :  G4VUserPhysicsList()                                                   |
|{                                                                                                        |
| .                                                                                                       |
| .                                                                                                       |
| defaultCutValue = 0.0001*mm;                                                                            |
| .                                                                                                       |
|  SetVerboseLevel(2);                                                                                    |
|}                                                                                                        |
| .                                                                                                       |
| .                                                                                                       |
| .                                                                                                       |
| .                                                                                                       |
|[..]                                                                                                     |
|      G4LivermorePhotoElectricModel* livermorephotoelectricmodel = new G4LivermorePhotoElectricModel();  |
|[..]                                                                                                     |
|      G4LivermoreIonisationModel* LivermoreIonisationModel = new G4LivermoreIonisationModel();           |
|[..]                                                                                                     |
|The information I need is contained in the G4LivermoreIonisationModel.cc                                 |
|.                                                                                                        |
|.                                                                                                        |
|.                                                                                                        |
| if (verboseLevel > 0)                                                                                   |
|    {                                                                                                    |
|      G4cout << "Livermore Ionisation model is initialized " << G4endl                                   |
|             << "Energy range: "                                                                         |
|             << LowEnergyLimit() / keV << " keV - "                                                      |
|             << HighEnergyLimit() / GeV << " GeV"                                                        |
|             << G4endl;                                                                                  |
|    }                                                                                                    |
|.                                                                                                        |
|.                                                                                                        |
|.                                                                                                        |
|.                                                                                                        |
|if (verboseLevel > 1)                                                                                    |
|    {                                                                                                    |
|      G4cout << "-----------------------------------------------------------" << G4endl;                 |
|      G4cout << "Energy balance from G4LivermoreIonisation" << G4endl;                                   |
|      G4cout << "Incoming primary energy: " << kineticEnergy/keV << " keV" << G4endl;                    |
|      G4cout << "-----------------------------------------------------------" << G4endl;                 |
|      G4cout << "Outgoing primary energy: " << finalKinEnergy/keV << " keV" << G4endl;                   |
|      G4cout << "Delta ray " << energyDelta/keV << " keV" << G4endl;                                     |
|      G4cout << "Fluorescence: " << (bindingEnergy-theEnergyDeposit)/keV << " keV" << G4endl;            |
|      G4cout << "Local energy deposit " << theEnergyDeposit/keV << " keV" << G4endl;                     |
|      G4cout << "Total final state: " << (finalKinEnergy+energyDelta+bindingEnergy+                      |
|                                          theEnergyDeposit)/keV << " keV" << G4endl;                     |
|      G4cout << "-----------------------------------------------------------" << G4endl;                 |
|    }                                                                                                    |
|.                                                                                                        |
|.                                                                                                        |
|.                                                                                                        |
|                                                                                                         |
|and in G4LivermoreIonisationModel.cc                                                                     |
|.                                                                                                        |
|.                                                                                                        |
|.                                                                                                        |
|if(verboseLevel>0) {                                                                                     |
|    G4cout << "Livermore PhotoElectric is constructed " << G4endl                                        |
|           << "Energy range: "                                                                           |
|           << lowEnergyLimit / eV << " eV - "                                                            |
|           << highEnergyLimit / GeV << " GeV"                                                            |
|           << G4endl;                                                                                    |
|  }                                                                                                      |
|.                                                                                                        |
|.                                                                                                        |
|."                                                                                                       |
|---------------------------------------------------------------------------------------------------------|

Those last few lines are actually only in  G4LivermorePhotoElectricModel.cc

|---------------------------------------------------------------------------------------------------------|
|"By setting the verboseLevel on 2 in the PhysicsList this information                                    |
|should be diplayed after starting the executable. But this is not the                                    |
|case. I get only some information on the display if I set                                                |
|/tracking/verbose 1 ore bigger,                                                                          |
|                                                                                                         |
|*********************************************************************************************************|
|* G4Track Information:   Particle = e-,   Track ID = 35,   Parent ID = 2                                 |
|*********************************************************************************************************|
|                                                                                                         |
|Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process       |
|    0   1.34 mm    947 um   4.54 mm   2.91 keV     0 eV      0 fm      0 fm   Detector     initStep      |
|    1   1.34 mm    947 um   4.54 mm      0 eV   2.91 keV   216 nm    216 nm   Detector     eIoni         |
|                                                                                                         |
|*********************************************************************************************************|
|                                                                                                         |
|but this is not the information listed in G4LivermoreIonisationModel.cc                                  |
|and G4LivermoreIonisationModel.cc.                                                                       |
|                                                                                                         |
|[..]"                                                                                                    |
|---------------------------------------------------------------------------------------------------------|

Good news: I can help you. Bad news: you are using Geant.

Try
LivermoreIonisationModel->SetVerboseLevel(2);
Unfortunately, this does not liberate the logging in
G4LivermorePhotoElectricModel  , which is an illogically programmed
class.  G4LivermorePhotoElectricModel  has no reasonable way to
liberate the logging excerpt which you quoted. One way to liberate it
is to replace
if(verboseLevel>0) {
    G4cout << "Livermore PhotoElectric is constructed " << G4endl
           << "Energy range: "
           << lowEnergyLimit / eV << " eV - "
           << highEnergyLimit / GeV << " GeV"
           << G4endl;
  }
by
  {
    G4cout << "Livermore PhotoElectric is constructed " << G4endl
           << "Energy range: "
           << lowEnergyLimit / eV << " eV - "
           << highEnergyLimit / GeV << " GeV"
           << G4endl;
  }
in  G4LivermorePhotoElectricModel.cc  and then to recompile Geant and
then to recompile your program. (There is one way to liberate the
logging code without modifying  G4LivermorePhotoElectricModel.cc  but
it is unreasonable.)

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

1 Ok: Re: LivermoreModel   (Mike - 28 Jul, 2011)
 Add Message Add Message
to: "Re: LivermoreModel"

 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 ]