source: trunk/examples/extended/electromagnetic/TestEm1/src/PrimaryGeneratorAction.cc @ 812

Last change on this file since 812 was 807, checked in by garnier, 16 years ago

update
File size: 4.1 KB
Line 
1//
2// ********************************************************************
3// * License and Disclaimer                                           *
4// *                                                                  *
5// * The  Geant4 software  is  copyright of the Copyright Holders  of *
6// * the Geant4 Collaboration.  It is provided  under  the terms  and *
7// * conditions of the Geant4 Software License,  included in the file *
8// * LICENSE and available at  http://cern.ch/geant4/license .  These *
9// * include a list of copyright holders.                             *
10// *                                                                  *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work  make  any representation or  warranty, express or implied, *
14// * regarding  this  software system or assume any liability for its *
15// * use.  Please see the license in the file  LICENSE  and URL above *
16// * for the full disclaimer and the limitation of liability.         *
17// *                                                                  *
18// * This  code  implementation is the result of  the  scientific and *
19// * technical work of the GEANT4 collaboration.                      *
20// * By using,  copying,  modifying or  distributing the software (or *
21// * any work based  on the software)  you  agree  to acknowledge its *
22// * use  in  resulting  scientific  publications,  and indicate your *
23// * acceptance of all terms of the Geant4 Software license.          *
24// ********************************************************************
25//
26//
27// $Id: PrimaryGeneratorAction.cc,v 1.4 2006/06/29 16:37:21 gunter Exp $
28// GEANT4 tag $Name:  $
29//
30//
31
32//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
33//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
34
35#include "PrimaryGeneratorAction.hh"
36
37#include "DetectorConstruction.hh"
38#include "PrimaryGeneratorMessenger.hh"
39
40#include "G4Event.hh"
41#include "G4ParticleTable.hh"
42#include "G4ParticleDefinition.hh"
43#include "Randomize.hh"
44
45//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
46
47PrimaryGeneratorAction::PrimaryGeneratorAction(
48                                               DetectorConstruction* DC)
49:Detector(DC)                                         
50{
51  particleGun  = new G4ParticleGun(1);
52  SetDefaultKinematic(1);
53  rndmBeam = 0.;
54   
55  //create a messenger for this class
56  gunMessenger = new PrimaryGeneratorMessenger(this); 
57}
58
59//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
60
61PrimaryGeneratorAction::~PrimaryGeneratorAction()
62{
63  delete particleGun;
64  delete gunMessenger; 
65}
66
67//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
68
69void PrimaryGeneratorAction::SetDefaultKinematic(G4int front)
70{
71  G4ParticleDefinition* particle
72           = G4ParticleTable::GetParticleTable()->FindParticle("e-");
73  particleGun->SetParticleDefinition(particle);
74  particleGun->SetParticleMomentumDirection(G4ThreeVector(1.,0.,0.));
75  particleGun->SetParticleEnergy(100*MeV);
76  G4double position = 0.*cm;
77  if (front) position = -0.5*(Detector->GetSize());
78  particleGun->SetParticlePosition(G4ThreeVector(position,0.*cm,0.*cm));
79}
80
81//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
82
83void PrimaryGeneratorAction::GeneratePrimaries(G4Event* anEvent)
84{
85  //this function is called at the begining of event
86  //
87  //randomize the beam, if requested.
88  if (rndmBeam > 0.) 
89    {
90      G4ThreeVector oldPosition = particleGun->GetParticlePosition();   
91      G4double rbeam = 0.5*(Detector->GetSize())*rndmBeam;
92      G4double x0 = oldPosition.x();
93      G4double y0 = oldPosition.y() + (2*G4UniformRand()-1.)*rbeam;
94      G4double z0 = oldPosition.z() + (2*G4UniformRand()-1.)*rbeam;
95      particleGun->SetParticlePosition(G4ThreeVector(x0,y0,z0));
96      particleGun->GeneratePrimaryVertex(anEvent);
97      particleGun->SetParticlePosition(oldPosition);     
98    }
99  else  particleGun->GeneratePrimaryVertex(anEvent); 
100}
101
102//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
103
Note: See TracBrowser for help on using the repository browser.