Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/human_phantom/src/G4HumanPhantomPrimaryGeneratorAction.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 /examples/advanced/human_phantom/src/G4HumanPhantomPrimaryGeneratorAction.cc (Version 11.3.0) and /examples/advanced/human_phantom/src/G4HumanPhantomPrimaryGeneratorAction.cc (Version 9.6.p3)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26 // Previous authors: G. Guerrieri, S. Guatelli <<  26 // Authors: S. Guatelli and M. G. Pia, INFN Genova, Italy
 27 // Authors (since 2007): S. Guatelli, Universi <<  27 // 
                                                   >>  28 // Based on code developed by the undergraduate student G. Guerrieri 
                                                   >>  29 // Note: this is a preliminary beta-version of the code; an improved 
                                                   >>  30 // version will be distributed in the next Geant4 public release, compliant
                                                   >>  31 // with the design in a forthcoming publication, and subject to a 
                                                   >>  32 // design and code review.
 28 //                                                 33 //
 29 #include "G4HumanPhantomPrimaryGeneratorAction     34 #include "G4HumanPhantomPrimaryGeneratorAction.hh"
                                                   >>  35 #include "G4HumanPhantomConstruction.hh"
                                                   >>  36 #include "G4HumanPhantomPrimaryGeneratorMessenger.hh"
                                                   >>  37 
 30 #include "globals.hh"                              38 #include "globals.hh"
 31 #include "G4SystemOfUnits.hh"                      39 #include "G4SystemOfUnits.hh"
 32 #include "G4Event.hh"                              40 #include "G4Event.hh"
 33 #include "G4GeneralParticleSource.hh"          <<  41 #include "G4ParticleGun.hh"
 34 #include "G4ParticleTable.hh"                      42 #include "G4ParticleTable.hh"
 35 #include "G4ParticleDefinition.hh"                 43 #include "G4ParticleDefinition.hh"
 36                                                    44 
 37 #include "G4RunManager.hh"                         45 #include "G4RunManager.hh"
 38 #include "G4ios.hh"                                46 #include "G4ios.hh"
 39                                                    47 
 40 G4HumanPhantomPrimaryGeneratorAction::G4HumanP     48 G4HumanPhantomPrimaryGeneratorAction::G4HumanPhantomPrimaryGeneratorAction()
                                                   >>  49   :beamKind("beamAlongZ"),worldLength(200.)
 41 {                                                  50 {
 42   particleGun = new G4GeneralParticleSource(); <<  51   G4int n_particle = 1;
                                                   >>  52 
                                                   >>  53   messenger= new G4HumanPhantomPrimaryGeneratorMessenger(this);
                                                   >>  54 
                                                   >>  55   particleGun = new G4ParticleGun(n_particle);
                                                   >>  56 
                                                   >>  57   probability.push_back(0.1667);
                                                   >>  58   probability.push_back(0.1667);
                                                   >>  59   probability.push_back(0.1667);
                                                   >>  60   probability.push_back(0.1667);
                                                   >>  61   probability.push_back(0.1666);
                                                   >>  62   probability.push_back(0.1666);
                                                   >>  63 
                                                   >>  64   // Default Particle
                                                   >>  65   G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
                                                   >>  66   
                                                   >>  67 particleGun->SetParticleDefinition(particleTable->FindParticle("geantino"));
                                                   >>  68   particleGun->SetParticleEnergy(100.*MeV);
 43 }                                                  69 }
 44                                                    70 
 45 G4HumanPhantomPrimaryGeneratorAction::~G4Human     71 G4HumanPhantomPrimaryGeneratorAction::~G4HumanPhantomPrimaryGeneratorAction()
 46 {                                                  72 {
 47   delete particleGun;                              73   delete particleGun;
                                                   >>  74   delete messenger;
 48 }                                                  75 }
 49                                                    76 
 50 void G4HumanPhantomPrimaryGeneratorAction::Gen     77 void G4HumanPhantomPrimaryGeneratorAction::GeneratePrimaries(G4Event* anEvent)
 51 {                                                  78 {
 52   particleGun -> GeneratePrimaryVertex(anEvent <<  79   if (beamKind == "beamAlongZ") GenerateBeamAlongZ();
                                                   >>  80   if (beamKind == "beamAlongY") GenerateBeamAlongY();
                                                   >>  81   if (beamKind == "beamAlongX") GenerateBeamAlongX();
                                                   >>  82   if (beamKind == "isotropicFlux") GenerateIsotropicFlux();
                                                   >>  83 
                                                   >>  84   particleGun->SetParticlePosition(G4ThreeVector(x0, y0,z0));
                                                   >>  85   particleGun->GeneratePrimaryVertex(anEvent);
                                                   >>  86 }
                                                   >>  87 
                                                   >>  88 void G4HumanPhantomPrimaryGeneratorAction::GenerateBeamAlongZ()
                                                   >>  89 {
                                                   >>  90   z0 = 0.5*(worldLength)*cm;
                                                   >>  91   y0 = (worldLength)*(G4UniformRand()-0.5)*cm;
                                                   >>  92   x0 = (worldLength)*(G4UniformRand()-0.5)*cm;
                                                   >>  93  
                                                   >>  94   //z0 = 0.5*(worldLength)*cm;
                                                   >>  95   //y0 = 50. * cm + 5.*(G4UniformRand()-0.5)*cm;
                                                   >>  96   //x0 = 10.*(G4UniformRand()-0.5)*cm;
                                                   >>  97 
                                                   >>  98    G4ThreeVector direction(0.,0.,-1.);
                                                   >>  99 
                                                   >> 100 
                                                   >> 101   particleGun->SetParticleMomentumDirection(direction);
                                                   >> 102 }
                                                   >> 103 
                                                   >> 104 void G4HumanPhantomPrimaryGeneratorAction::GenerateBeamAlongX()
                                                   >> 105 {
                                                   >> 106   x0 = -0.5*(worldLength)*cm;
                                                   >> 107   y0 = (worldLength)*(G4UniformRand()-0.5)*cm;
                                                   >> 108   z0 = (worldLength)*(G4UniformRand()-0.5)*cm;
                                                   >> 109  
                                                   >> 110   G4ThreeVector direction(1.,0.,0.);
                                                   >> 111   particleGun->SetParticleMomentumDirection(direction);
                                                   >> 112 }
                                                   >> 113 void G4HumanPhantomPrimaryGeneratorAction::GenerateBeamAlongY()
                                                   >> 114 {
                                                   >> 115   y0 = 0.5*(worldLength)*cm;
                                                   >> 116   x0 = (worldLength)*(G4UniformRand()-0.5)*cm;
                                                   >> 117   z0 = (worldLength)*(G4UniformRand()-0.5)*cm;
                                                   >> 118  
                                                   >> 119   G4ThreeVector direction(0.,-1.,0.);
                                                   >> 120   particleGun->SetParticleMomentumDirection(direction);
                                                   >> 121 }
                                                   >> 122 
                                                   >> 123 void G4HumanPhantomPrimaryGeneratorAction::GenerateIsotropicFlux()
                                                   >> 124 {
                                                   >> 125   G4double random = G4UniformRand();
                                                   >> 126   G4double sum = 0.;
                                                   >> 127   G4int i = 0;
                                                   >> 128 
                                                   >> 129   while(sum<random){sum += probability[i]; i++;}
                                                   >> 130   
                                                   >> 131   if(i==1) 
                                                   >> 132     {
                                                   >> 133       z0 = -0.5*(worldLength-2.)*cm;
                                                   >> 134       y0 = (worldLength)*(G4UniformRand()-0.5)*cm;
                                                   >> 135       x0 = (worldLength)*(G4UniformRand()-0.5)*cm;
                                                   >> 136     }
                                                   >> 137 
                                                   >> 138   if(i==2) 
                                                   >> 139     {
                                                   >> 140       y0 = -0.5*(worldLength-2.)*cm;
                                                   >> 141       z0 = (worldLength)*(G4UniformRand()-0.5)*cm;
                                                   >> 142       x0 = (worldLength)*(G4UniformRand()-0.5)*cm;
                                                   >> 143     }
                                                   >> 144 
                                                   >> 145   if(i==3) 
                                                   >> 146     {
                                                   >> 147       x0 = -0.5*(worldLength-2.)*cm;
                                                   >> 148       z0 = (worldLength)*(G4UniformRand()-0.5)*cm;
                                                   >> 149       y0 = (worldLength)*(G4UniformRand()-0.5)*cm;
                                                   >> 150     }
                                                   >> 151 
                                                   >> 152   if (i==4)
                                                   >> 153     {
                                                   >> 154       z0 = 0.5*(worldLength-2.)*cm;
                                                   >> 155       y0 = (worldLength)*(G4UniformRand()-0.5)*cm;
                                                   >> 156       x0 = (worldLength)*(G4UniformRand()-0.5)*cm;
                                                   >> 157     }
                                                   >> 158 
                                                   >> 159  if(i==5) 
                                                   >> 160     {
                                                   >> 161       y0 = 0.5*(worldLength-2.)*cm;
                                                   >> 162       z0 = (worldLength)*(G4UniformRand()-0.5)*cm;
                                                   >> 163       x0 = (worldLength)*(G4UniformRand()-0.5)*cm;
                                                   >> 164     }
                                                   >> 165 
                                                   >> 166   if(i==6) 
                                                   >> 167     {
                                                   >> 168       x0 = 0.5*(worldLength-2.)*cm;
                                                   >> 169       z0 = (worldLength)*(G4UniformRand()-0.5)*cm;
                                                   >> 170       y0 = (worldLength)*(G4UniformRand()-0.5)*cm;
                                                   >> 171     } 
                                                   >> 172 
                                                   >> 173   G4double a,b,c;
                                                   >> 174   G4double n;
                                                   >> 175   do{
                                                   >> 176     a = (G4UniformRand()-0.5)/0.5;
                                                   >> 177     b = (G4UniformRand()-0.5)/0.5; 
                                                   >> 178     c = (G4UniformRand()-0.5)/0.5;
                                                   >> 179     n = a*a+b*b+c*c;
                                                   >> 180   }while(n > 1 || n == 0.0);
                                                   >> 181   n = std::sqrt(n);
                                                   >> 182   a /= n;
                                                   >> 183   b /= n;
                                                   >> 184   c /= n;
                                                   >> 185 
                                                   >> 186   G4ThreeVector direction(a,b,c);
                                                   >> 187   particleGun->SetParticleMomentumDirection(direction);  
                                                   >> 188 }
                                                   >> 189 
                                                   >> 190 void G4HumanPhantomPrimaryGeneratorAction::SetBeam(G4String beam)
                                                   >> 191 {
                                                   >> 192   if((beam == "beamAlongZ")||(beam == "beamAlongX")||
                                                   >> 193   (beam == "beamAlongY")||(beam == "isotropicFlux")) beamKind = beam;
                                                   >> 194   
                                                   >> 195   else G4cout<<"This option is not valid "<<
                                                   >> 196                "---> beamAlongZ/beamAlongY/beamAlongX/isotropicFlux"
                                                   >> 197              <<G4endl;
 53 }                                                 198 }
 54                                                   199 
 55                                                   200