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

None Re: Secondary electrons 

Forum: Physics List
Re: None Secondary electrons (Davinder Siwal)
Re: None Re: Secondary electrons (John Allison)
Date: 10 May, 2010
From: Davinder Siwal <Davinder Siwal>

First i would like to thanks to you for reply

This is waht exactly i mean that i am shooting gamma to the germanium detector and interested in the pulse shape induce on the electrodes so i define a simple germaniun detector in cylindrical shape and and shooting 300 kev gammas and i am getting the energy deposition by gamma

so here is the output on my screen

========= Table of registered couples ==============================

Index : 0     used in the geometry : Yes     recalculation needed : No 
 Material : Air
 Range cuts        :  gamma 1 mm     e- 1 mm     e+ 1 mm 
 Energy thresholds :  gamma 990 eV     e- 990 eV     e+ 990 eV 
 Region(s) which use this couple : 
    DefaultRegionForTheWorld

Index : 1     used in the geometry : Yes     recalculation needed : No 
 Material : Germanium
 Range cuts        :  gamma 1 mm     e- 1 mm     e+ 1 mm 
 Energy thresholds :  gamma 20.6936 keV    e- 863.242 keV    e+ 821.764 keV
 Region(s) which use this couple : 
    DefaultRegionForTheWorld

====================================================================

Start closing geometry. G4GeometryManager::ReportVoxelStats -- Voxel Statistics

    Total memory consumed for geometry optimisation:   0 kByte
    Total CPU time elapsed for geometry optimisation: 0 seconds

######### Run Begins 0 START

Start Run processing.

------------> Begin of Event:0 0

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

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0      200        0        0       0.3        0        0         0     expHall initStep
    1      225        0        0       0.3        0       25        25 hpge_detector Transportation
    2      233        0        0         0   0.0111     7.75      32.8 hpge_detector phot
  trackID: 1  energy deposit: 11.103 keV  position: 23.27520990773711 0 0 cm time: 83.39102379953802 ps 

*********************************************************************************************************
* G4Track Information:   Particle = e-,   Track ID = 2,   Parent ID = 1
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng extVolume ProcName
    0      233        0        0     0.289        0        0         0 hpge_detector initStep
    1      233   0.0156 -0.00657         0    0.289    0.235     0.235 hpge_detector eIoni
  trackID: 2  energy deposit: 288.897 keV  position: 23.27879642934676 0.001564232073245215 -0.0006568775037383886 cm time: 109.2492429458352 ps 
Total no. of hits by this event = 2
Total energy deposited by this event = 300 keV
Total events impinges = 1

----------------->End of event: 0
G4VisManager: Using G4TrajectoryDrawByCharge as default trajectory model.
See commands in /vis/modeling/trajectories/ for other options.
HepRepFile writing to G4Data0.heprep
Run terminated.
Run Summary
  Number of events processed : 1
  User=0.01s Real=0s Sys=0s

######### Run Completed ##############

but i don't know how much electron hole pair it is forming for each gamma energy deposition that is i would like to know the no. of electron hole pair it is forming for each interaction

Here is my physics list

#include "PlanarPhysicsList.hh"
#include "G4ParticleTypes.hh"

PlanarPhysicsList::PlanarPhysicsList() {;}

PlanarPhysicsList::~PlanarPhysicsList() {;}

void PlanarPhysicsList::ConstructParticle()
{
  // In this method, static member functions should be called
  // for all particles which you want to use.
  // This ensures that objects of these particle types will be
  // created in the program. 

  ConstructBosons();
  ConstructLeptons();
}

void PlanarPhysicsList::ConstructProcess()
{
  // Define transportation process

  AddTransportation();
  ConstructEM();
}

void PlanarPhysicsList::ConstructBosons()
{
  G4Gamma::GammaDefinition();
}

void PlanarPhysicsList::ConstructLeptons()
{
  G4Electron::ElectronDefinition();
  G4Positron::PositronDefinition();
}

#include "G4ComptonScattering.hh"
#include "G4GammaConversion.hh"
#include "G4PhotoElectricEffect.hh"
#include "G4RegionStore.hh"

#include "G4ProcessManager.hh"
#include "G4StepLimiter.hh"

#include "G4MultipleScattering.hh"

#include "G4eIonisation.hh"
#include "G4eBremsstrahlung.hh"
#include "G4eplusAnnihilation.hh"

#include "G4Gamma.hh"
#include "G4Electron.hh"
#include "G4Positron.hh"

void PlanarPhysicsList::ConstructEM()
{
    G4ProcessManager* pManager =0;// particle->GetProcessManager();
    G4cout<<"REGISTER PARTICLES ARE :-"<<G4endl;     
      // gamma
	pManager = G4Gamma::Gamma()->GetProcessManager();
	G4cout<<"Gamma"<<G4endl;
	pManager->AddDiscreteProcess(new G4PhotoElectricEffect);
	pManager->AddDiscreteProcess(new G4ComptonScattering);
	pManager->AddDiscreteProcess(new G4GammaConversion);

      //electron
	 pManager = G4Electron::Electron()->GetProcessManager();
	 G4cout<<"Electron"<<G4endl;
	 G4VProcess* theeminusMultipleScattering = new G4MultipleScattering();
	 G4VProcess* theeminusIonisation         = new G4eIonisation();
	 G4VProcess* theeminusBremsstrahlung     = new G4eBremsstrahlung();//"eBrem", 10*GeV);

      //  add process
	 pManager->AddProcess(theeminusMultipleScattering);//, -1,1,1);
	 pManager->AddProcess(theeminusIonisation);//, -1,2,2);
	 pManager->AddProcess(theeminusBremsstrahlung);//, -1,3,3);

	 // set ordering for AlongStepDoIt

   pManager->SetProcessOrdering(theeminusMultipleScattering, idxAlongStep,1);
   pManager->SetProcessOrdering(theeminusIonisation,         idxAlongStep,2);
   //pManager->SetProcessOrdering(theeminusBremsstrahlung,     idxAlongStep,3);
   //
   // set ordering for PostStepDoIt
   pManager->SetProcessOrdering(theeminusMultipleScattering, idxPostStep,1);
   pManager->SetProcessOrdering(theeminusIonisation,         idxPostStep,2);
   //pManager->SetProcessOrdering(theeminusBremsstrahlung,     idxPostStep,3);

     //positron
	 pManager = G4Positron::Positron()->GetProcessManager();

	 G4cout<<"Positron"<<G4endl;
	 G4cout<<G4endl;
        G4VProcess* theeplusMultipleScattering = new G4MultipleScattering();
	G4VProcess* theeplusIonisation         = new G4eIonisation();
	G4VProcess* theeplusBremsstrahlung     = new G4eBremsstrahlung("eBrem", 10*GeV);
	G4VProcess* theeplusAnnihilation       = new G4eplusAnnihilation();

	pManager->AddProcess(theeplusMultipleScattering,-1,1,1);
	pManager->AddProcess(theeplusIonisation, -1,2,2);
	pManager->AddProcess(theeplusBremsstrahlung ,-1,3,3);
	pManager->AddProcess(theeplusAnnihilation,0,-1,4);

}

void PlanarPhysicsList::SetCuts() {

    SetCutsWithDefault();   

}

Another question is that i am not able to see these electron hole pair inside the detector so here is my detector visualisation attachment

I will be highly thankful if you guide me in this direction

Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2010/05/10/10.00-95295-Screenshot.png

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

1 None: Re: Secondary electrons   (John Allison - 11 May, 2010)
(_ Ok: Re: Secondary electrons   (Davinder Siwal - 13 May, 2010)
(_ None: Re: Secondary electrons   (John Allison - 14 May, 2010)
(_ Ok: Re: Secondary electrons   (Davinder Siwal - 15 May, 2010)
 Add Message Add Message
to: "Re: Secondary electrons"

 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 ]