Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/event/src/G4ParticleGun.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /event/src/G4ParticleGun.cc (Version 11.3.0) and /event/src/G4ParticleGun.cc (Version 4.0.p1)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                    <<   3 // * DISCLAIMER                                                       *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th <<   5 // * The following disclaimer summarizes all the specific disclaimers *
  6 // * the Geant4 Collaboration.  It is provided <<   6 // * of contributors to this software. The specific disclaimers,which *
  7 // * conditions of the Geant4 Software License <<   7 // * govern, are listed with their locations in:                      *
  8 // * LICENSE and available at  http://cern.ch/ <<   8 // *   http://cern.ch/geant4/license                                  *
  9 // * include a list of copyright holders.      << 
 10 // *                                                9 // *                                                                  *
 11 // * Neither the authors of this software syst     10 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     11 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     12 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     13 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file  <<  14 // * use.                                                             *
 16 // * for the full disclaimer and the limitatio << 
 17 // *                                               15 // *                                                                  *
 18 // * This  code  implementation is the result  <<  16 // * This  code  implementation is the  intellectual property  of the *
 19 // * technical work of the GEANT4 collaboratio <<  17 // * GEANT4 collaboration.                                            *
 20 // * By using,  copying,  modifying or  distri <<  18 // * By copying,  distributing  or modifying the Program (or any work *
 21 // * any work based  on the software)  you  ag <<  19 // * based  on  the Program)  you indicate  your  acceptance of  this *
 22 // * use  in  resulting  scientific  publicati <<  20 // * statement, and all its terms.                                    *
 23 // * acceptance of all terms of the Geant4 Sof << 
 24 // *******************************************     21 // ********************************************************************
 25 //                                                 22 //
 26 // G4ParticleGun class implementation          << 
 27 //                                                 23 //
 28 // Author: Makoto Asai, 1997                   <<  24 // $Id: G4ParticleGun.cc,v 1.6 2001/11/20 23:21:41 asaim Exp $
 29 // ------------------------------------------- <<  25 // GEANT4 tag $Name: geant4-04-00 $
                                                   >>  26 //
 30                                                    27 
                                                   >>  28 // G4ParticleGun
 31 #include "G4ParticleGun.hh"                        29 #include "G4ParticleGun.hh"
 32 #include "G4SystemOfUnits.hh"                  << 
 33 #include "G4PrimaryParticle.hh"                    30 #include "G4PrimaryParticle.hh"
 34 #include "G4ParticleGunMessenger.hh"               31 #include "G4ParticleGunMessenger.hh"
 35 #include "G4Event.hh"                              32 #include "G4Event.hh"
 36 #include "G4ios.hh"                                33 #include "G4ios.hh"
 37                                                    34 
 38 G4ParticleGun::G4ParticleGun()                     35 G4ParticleGun::G4ParticleGun()
 39 {                                                  36 {
 40   SetInitialValues();                              37   SetInitialValues();
 41 }                                                  38 }
 42                                                    39 
 43 G4ParticleGun::G4ParticleGun(G4int numberofpar     40 G4ParticleGun::G4ParticleGun(G4int numberofparticles)
 44 {                                                  41 {
 45   SetInitialValues();                              42   SetInitialValues();
 46   NumberOfParticlesToBeGenerated = numberofpar     43   NumberOfParticlesToBeGenerated = numberofparticles;
 47 }                                                  44 }
 48                                                    45 
 49 G4ParticleGun::G4ParticleGun(G4ParticleDefinit <<  46 G4ParticleGun::G4ParticleGun
 50                              G4int numberofpar <<  47     (G4ParticleDefinition * particleDef, G4int numberofparticles)
 51 {                                                  48 {
 52   SetInitialValues();                              49   SetInitialValues();
 53   NumberOfParticlesToBeGenerated = numberofpar     50   NumberOfParticlesToBeGenerated = numberofparticles;
 54   SetParticleDefinition( particleDef );        <<  51   particle_definition = particleDef;
 55 }                                                  52 }
 56                                                    53 
 57 void G4ParticleGun::SetInitialValues()             54 void G4ParticleGun::SetInitialValues()
 58 {                                                  55 {
 59   NumberOfParticlesToBeGenerated = 1;              56   NumberOfParticlesToBeGenerated = 1;
 60   particle_definition = nullptr;               <<  57   particle_definition = 0;
 61   G4ThreeVector zero;                              58   G4ThreeVector zero;
 62   particle_momentum_direction = (G4ParticleMom     59   particle_momentum_direction = (G4ParticleMomentum)zero;
 63   particle_energy = 0.0;                           60   particle_energy = 0.0;
 64   particle_momentum = 0.0;                     << 
 65   particle_position = zero;                        61   particle_position = zero;
 66   particle_time = 0.0;                             62   particle_time = 0.0;
 67   particle_polarization = zero;                    63   particle_polarization = zero;
 68   particle_charge = 0.0;                           64   particle_charge = 0.0;
 69   theMessenger = new G4ParticleGunMessenger(th     65   theMessenger = new G4ParticleGunMessenger(this);
 70 }                                                  66 }
 71                                                    67 
 72 G4ParticleGun::~G4ParticleGun()                    68 G4ParticleGun::~G4ParticleGun()
 73 {                                                  69 {
 74   delete theMessenger;                             70   delete theMessenger;
 75 }                                                  71 }
 76                                                    72 
 77 void G4ParticleGun::                           <<  73 void G4ParticleGun::SetParticleDefinition
 78 SetParticleDefinition(G4ParticleDefinition* aP <<  74                  (G4ParticleDefinition * aParticleDefinition)
 79 {                                                  75 { 
 80   if(aParticleDefinition == nullptr)           << 
 81   {                                            << 
 82     G4Exception("G4ParticleGun::SetParticleDef << 
 83                 FatalException, "Null pointer  << 
 84   }                                            << 
 85   if(aParticleDefinition->IsShortLived())      << 
 86   {                                            << 
 87     if(aParticleDefinition->GetDecayTable() == << 
 88     {                                          << 
 89       G4ExceptionDescription ED;               << 
 90       ED << "G4ParticleGun does not support sh << 
 91          << "particle without a valid decay ta << 
 92       ED << "G4ParticleGun::SetParticleDefinit << 
 93          << aParticleDefinition->GetParticleNa << 
 94       G4Exception("G4ParticleGun::SetParticleD << 
 95                   JustWarning, ED);            << 
 96       return;                                  << 
 97     }                                          << 
 98   }                                            << 
 99   particle_definition = aParticleDefinition;       76   particle_definition = aParticleDefinition; 
100   particle_charge = particle_definition->GetPD     77   particle_charge = particle_definition->GetPDGCharge();
101   if(particle_momentum>0.0)                    << 
102   {                                            << 
103     G4double mass = particle_definition->GetPD << 
104     particle_energy =                          << 
105       std::sqrt(particle_momentum*particle_mom << 
106   }                                            << 
107 }                                              << 
108                                                << 
109 void G4ParticleGun::SetParticleEnergy(G4double << 
110 {                                              << 
111   particle_energy = aKineticEnergy;            << 
112   if(particle_momentum>0.0)                    << 
113   {                                            << 
114     if(particle_definition != nullptr)         << 
115     {                                          << 
116       G4cout << "G4ParticleGun::" << particle_ << 
117              << G4endl;                        << 
118     }                                          << 
119     else                                       << 
120     {                                          << 
121       G4cout << "G4ParticleGun::" << " " << G4 << 
122     }                                          << 
123     G4cout << " was defined in terms of Moment << 
124            << particle_momentum/GeV << "GeV/c" << 
125     G4cout << " is now defined in terms of Kin << 
126            << particle_energy/GeV   << "GeV"   << 
127     particle_momentum = 0.0;                   << 
128   }                                            << 
129 }                                              << 
130                                                << 
131 void G4ParticleGun::SetParticleMomentum(G4doub << 
132 {                                              << 
133   if(particle_energy>0.0)                      << 
134   {                                            << 
135     if(particle_definition != nullptr)         << 
136     {                                          << 
137       G4cout << "G4ParticleGun::" << particle_ << 
138              << G4endl;                        << 
139     }                                          << 
140     else                                       << 
141     {                                          << 
142       G4cout << "G4ParticleGun::" << " " << G4 << 
143     }                                          << 
144     G4cout << " was defined in terms of Kineti << 
145            << particle_energy/GeV << "GeV"   < << 
146     G4cout << " is now defined in terms Moment << 
147            << aMomentum/GeV       << "GeV/c" < << 
148   }                                            << 
149   if(particle_definition == nullptr)           << 
150   {                                            << 
151     G4cout << "Particle Definition not defined << 
152            << G4endl;                          << 
153     G4cout << "Zero Mass is assumed" << G4endl << 
154     particle_momentum = aMomentum;             << 
155     particle_energy = aMomentum;               << 
156   }                                            << 
157   else                                         << 
158   {                                            << 
159     G4double mass = particle_definition->GetPD << 
160     particle_momentum = aMomentum;             << 
161     particle_energy =                          << 
162       std::sqrt(particle_momentum*particle_mom << 
163   }                                            << 
164 }                                                  78 }
165                                                    79  
166 void G4ParticleGun::SetParticleMomentum(G4Part     80 void G4ParticleGun::SetParticleMomentum(G4ParticleMomentum aMomentum)
167 {                                                  81 {
168   if(particle_energy>0.0)                      <<  82   if(particle_definition==0)
169   {                                                83   {
170     if(particle_definition != nullptr)         <<  84     G4cout <<"Particle Definition not defined yet for G4ParticleGun"<< G4endl;
171     {                                          <<  85     G4cout <<"Zero Mass is assumed"<<G4endl;
172       G4cout << "G4ParticleGun::" << particle_ <<  86     particle_momentum_direction =  aMomentum.unit();
173              << G4endl;                        << 
174     }                                          << 
175     else                                       << 
176     {                                          << 
177       G4cout << "G4ParticleGun::" << " " << G4 << 
178     }                                          << 
179     G4cout << " was defined in terms of Kineti << 
180            << particle_energy/GeV << "GeV"   < << 
181     G4cout << " is now defined in terms Moment << 
182            << aMomentum.mag()/GeV << "GeV/c" < << 
183   }                                            << 
184   if(particle_definition == nullptr)           << 
185   {                                            << 
186     G4cout << "Particle Definition not defined << 
187            << G4endl;                          << 
188     G4cout << "Zero Mass is assumed" << G4endl << 
189     particle_momentum_direction = aMomentum.un << 
190     particle_momentum = aMomentum.mag();       << 
191     particle_energy = aMomentum.mag();             87     particle_energy = aMomentum.mag();
192   }                                                88   } 
193   else                                             89   else 
194   {                                                90   {
195     G4double mass =  particle_definition->GetP     91     G4double mass =  particle_definition->GetPDGMass();
196     particle_momentum = aMomentum.mag();       <<  92     G4double p = aMomentum.mag();
197     particle_momentum_direction = aMomentum.un <<  93     particle_momentum_direction =  aMomentum.unit();
198     particle_energy =                          <<  94     if ((particle_energy>0.0)&&(abs(particle_energy+mass-sqrt(p*p+mass*mass))>keV))
199       std::sqrt(particle_momentum*particle_mom <<  95     {
                                                   >>  96       G4cout << "G4ParticleGun::" << particle_definition->GetParticleName() << G4endl;
                                                   >>  97       G4cout << "  KineticEnergy and Momentum could be inconsistent" << G4endl;
                                                   >>  98       G4cout << " (Momentum:" << p/GeV << " GeV/c";
                                                   >>  99       G4cout << "  Mass:" << mass/GeV << " GeV/c/c)" << G4endl;
                                                   >> 100       G4cout << "  KineticEnergy is overwritten!! ";
                                                   >> 101       G4cout << particle_energy/GeV << "->";
                                                   >> 102       G4cout << (sqrt(p*p+mass*mass)-mass)/GeV << "GeV" << G4endl;
                                                   >> 103     }
                                                   >> 104     particle_energy = sqrt(p*p+mass*mass)-mass;
200   }                                               105   }
201 }                                                 106 }
202                                                   107 
203 void G4ParticleGun::GeneratePrimaryVertex(G4Ev    108 void G4ParticleGun::GeneratePrimaryVertex(G4Event* evt)
204 {                                                 109 {
205   if(particle_definition == nullptr)           << 110   if(particle_definition==0) return;
206   {                                            << 
207     G4ExceptionDescription ED;                 << 
208     ED << "Particle definition is not defined. << 
209     ED << "G4ParticleGun::SetParticleDefinitio << 
210        << G4endl;                              << 
211     G4Exception("G4ParticleGun::GeneratePrimar << 
212                 FatalException, ED);           << 
213     return;                                    << 
214   }                                            << 
215                                                   111 
216   // Create a new vertex                       << 112   // create a new vertex
217   //                                           << 113   G4PrimaryVertex* vertex = 
218   auto* vertex =                               << 
219     new G4PrimaryVertex(particle_position,part    114     new G4PrimaryVertex(particle_position,particle_time);
220                                                   115 
221   // Create new primaries and set them to the  << 116   // create new primaries and set them to the vertex
222   //                                           << 
223   G4double mass =  particle_definition->GetPDG    117   G4double mass =  particle_definition->GetPDGMass();
224   for( G4int i=0; i<NumberOfParticlesToBeGener << 118   G4double energy = particle_energy + mass;
                                                   >> 119   G4double pmom = sqrt(energy*energy-mass*mass);
                                                   >> 120   G4double px = pmom*particle_momentum_direction.x();
                                                   >> 121   G4double py = pmom*particle_momentum_direction.y();
                                                   >> 122   G4double pz = pmom*particle_momentum_direction.z();
                                                   >> 123   for( G4int i=0; i<NumberOfParticlesToBeGenerated; i++ )
225   {                                               124   {
226     auto* particle =                           << 125     G4PrimaryParticle* particle =
227       new G4PrimaryParticle(particle_definitio << 126       new G4PrimaryParticle(particle_definition,px,py,pz);
228     particle->SetKineticEnergy( particle_energ << 
229     particle->SetMass( mass );                    127     particle->SetMass( mass );
230     particle->SetMomentumDirection( particle_m << 
231     particle->SetCharge( particle_charge );       128     particle->SetCharge( particle_charge );
232     particle->SetPolarization(particle_polariz    129     particle->SetPolarization(particle_polarization.x(),
233                               particle_polariz << 130                                particle_polarization.y(),
234                               particle_polariz << 131                                particle_polarization.z());
235     particle->SetWeight( particle_weight );    << 
236     vertex->SetPrimary( particle );               132     vertex->SetPrimary( particle );
237   }                                               133   }
                                                   >> 134 
238   evt->AddPrimaryVertex( vertex );                135   evt->AddPrimaryVertex( vertex );
239 }                                                 136 }
                                                   >> 137 
                                                   >> 138 
240                                                   139