Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/field/field02/src/F02PrimaryGeneratorAction.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/extended/field/field02/src/F02PrimaryGeneratorAction.cc (Version 11.3.0) and /examples/extended/field/field02/src/F02PrimaryGeneratorAction.cc (Version 9.6.p4)


  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 /// \file field/field02/src/F02PrimaryGenerato     26 /// \file field/field02/src/F02PrimaryGeneratorAction.cc
 27 /// \brief Implementation of the F02PrimaryGen     27 /// \brief Implementation of the F02PrimaryGeneratorAction class
 28 //                                                 28 //
 29 //                                             <<  29 // $Id$
 30 //                                             <<  30 // 
 31 //                                             << 
 32 //....oooOO0OOooo........oooOO0OOooo........oo << 
 33 //....oooOO0OOooo........oooOO0OOooo........oo << 
 34                                                    31 
 35 #include "F02PrimaryGeneratorAction.hh"        <<  32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
                                                   >>  33 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 36                                                    34 
                                                   >>  35 #include "F02PrimaryGeneratorAction.hh"
 37 #include "F02DetectorConstruction.hh"              36 #include "F02DetectorConstruction.hh"
 38 #include "F02PrimaryGeneratorMessenger.hh"         37 #include "F02PrimaryGeneratorMessenger.hh"
 39                                                    38 
 40 #include "G4Event.hh"                              39 #include "G4Event.hh"
 41 #include "G4ParticleDefinition.hh"             << 
 42 #include "G4ParticleGun.hh"                        40 #include "G4ParticleGun.hh"
 43 #include "G4ParticleTable.hh"                      41 #include "G4ParticleTable.hh"
                                                   >>  42 #include "G4ParticleDefinition.hh"
 44 #include "G4PhysicalConstants.hh"                  43 #include "G4PhysicalConstants.hh"
 45 #include "G4SystemOfUnits.hh"                      44 #include "G4SystemOfUnits.hh"
                                                   >>  45 #include "G4ios.hh"
 46 #include "Randomize.hh"                            46 #include "Randomize.hh"
 47                                                    47 
 48 //....oooOO0OOooo........oooOO0OOooo........oo <<  48 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 49                                                <<  49  
 50 G4ParticleDefinition* F02PrimaryGeneratorActio <<  50 G4String F02PrimaryGeneratorAction::fgPrimaryParticleName = "proton" ; 
 51                                                <<  51 
 52 //....oooOO0OOooo........oooOO0OOooo........oo <<  52 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 53                                                <<  53 
 54 F02PrimaryGeneratorAction::F02PrimaryGenerator <<  54 F02PrimaryGeneratorAction::F02PrimaryGeneratorAction(
                                                   >>  55                                             F02DetectorConstruction* det)
                                                   >>  56   : fParticleGun(0),
                                                   >>  57     fDetector(det), 
                                                   >>  58     fMessenger(0),
                                                   >>  59     fRndmFlag("off"),
                                                   >>  60     fXVertex(0.),
                                                   >>  61     fYVertex(0.),
                                                   >>  62     fZVertex(0.),
                                                   >>  63     fVertexDefined(false)
 55 {                                                  64 {
 56   G4int n_particle = 1;                            65   G4int n_particle = 1;
 57   fParticleGun = new G4ParticleGun(n_particle) <<  66   fParticleGun  = new G4ParticleGun(n_particle);
 58                                                <<  67   
 59   // create a messenger for this class             68   // create a messenger for this class
 60   fGunMessenger = new F02PrimaryGeneratorMesse <<  69   //
                                                   >>  70   fMessenger = new F02PrimaryGeneratorMessenger(this);
 61                                                    71 
 62   // default particle kinematic                    72   // default particle kinematic
 63                                                    73 
 64   G4ParticleTable* particleTable = G4ParticleT     74   G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
 65   G4String particleName;                           75   G4String particleName;
 66   G4ParticleDefinition* particle = particleTab <<  76   G4ParticleDefinition* particle
                                                   >>  77     = particleTable->FindParticle(particleName="proton");
 67   fParticleGun->SetParticleDefinition(particle     78   fParticleGun->SetParticleDefinition(particle);
                                                   >>  79   
                                                   >>  80   fgPrimaryParticleName = particle->GetParticleName() ;
 68                                                    81 
 69   fgPrimaryParticle = particle;                <<  82   fParticleGun->SetParticleMomentumDirection(G4ThreeVector(0.,0.,1.));
                                                   >>  83   fParticleGun->SetParticleEnergy(100.*GeV);
 70                                                    84 
 71   fParticleGun->SetParticleMomentumDirection(G <<  85   fZVertex = 0. ; //  -0.5*(fDetector->GetAbsorberThickness());
 72   fParticleGun->SetParticleEnergy(100. * GeV); <<  86   fParticleGun->SetParticlePosition(G4ThreeVector(fXVertex,fYVertex,fZVertex));
 73                                                    87 
 74   fZVertex = 0.;  //  -0.5*(fDetector->GetAbso << 
 75   fParticleGun->SetParticlePosition(G4ThreeVec << 
 76 }                                                  88 }
 77                                                    89 
 78 //....oooOO0OOooo........oooOO0OOooo........oo <<  90 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 79                                                    91 
 80 F02PrimaryGeneratorAction::~F02PrimaryGenerato     92 F02PrimaryGeneratorAction::~F02PrimaryGeneratorAction()
 81 {                                                  93 {
 82   delete fParticleGun;                             94   delete fParticleGun;
 83   delete fGunMessenger;                        <<  95   delete fMessenger;
 84 }                                                  96 }
 85                                                    97 
 86 //....oooOO0OOooo........oooOO0OOooo........oo <<  98 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 87                                                    99 
 88 void F02PrimaryGeneratorAction::GeneratePrimar    100 void F02PrimaryGeneratorAction::GeneratePrimaries(G4Event* anEvent)
 89 {                                                 101 {
 90   // this function is called at the begining o    102   // this function is called at the begining of event
 91   //                                           << 103   // 
 92   fgPrimaryParticle = fParticleGun->GetParticl << 104   fgPrimaryParticleName 
 93                                                << 105     = fParticleGun->GetParticleDefinition()->GetParticleName() ;
 94   G4double x0, y0, z0;                         << 106 
 95   if (fVertexDefined) {                        << 107   G4double x0,y0,z0 ;
 96     x0 = fXVertex;                             << 108   if (fVertexDefined)
 97     y0 = fYVertex;                             << 109   {
 98     z0 = fZVertex;                             << 110     x0 = fXVertex ;
                                                   >> 111     y0 = fYVertex ;
                                                   >> 112     z0 = fZVertex ;
 99   }                                               113   }
100   else {                                       << 114   else
101     x0 = 0.;                                   << 115   {
102     y0 = 0.;                                   << 116     x0 = 0. ;
103     z0 = 0.;  // -0.5*(fDetector->GetWorldSize << 117     y0 = 0. ;
                                                   >> 118     z0 = 0. ; // -0.5*(fDetector->GetWorldSizeZ()) ;
104   }                                               119   }
105                                                   120 
106   G4double r0, phi0;                           << 121   G4double r0,phi0 ;
107   if (fRndmFlag == "on") {                     << 122   if (fRndmFlag == "on")
108     r0 = (fDetector->GetAbsorberRadius()) * st << 123   {
109     phi0 = twopi * G4UniformRand();            << 124       r0 = (fDetector->GetAbsorberRadius())*std::sqrt(G4UniformRand());
110     x0 = r0 * std::cos(phi0);                  << 125       phi0 = twopi*G4UniformRand();
111     y0 = r0 * std::sin(phi0);                  << 126       x0 = r0*std::cos(phi0);
112   }                                            << 127       y0 = r0*std::sin(phi0);
                                                   >> 128   } 
113                                                   129 
114   fParticleGun->SetParticlePosition(G4ThreeVec << 130   fParticleGun->SetParticlePosition(G4ThreeVector(x0,y0,z0));
115   fParticleGun->GeneratePrimaryVertex(anEvent)    131   fParticleGun->GeneratePrimaryVertex(anEvent);
116 }                                                 132 }
117                                                   133 
118 //....oooOO0OOooo........oooOO0OOooo........oo << 134 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
119                                                   135 
120 G4String F02PrimaryGeneratorAction::GetPrimary    136 G4String F02PrimaryGeneratorAction::GetPrimaryName()
121 {                                                 137 {
122   return fgPrimaryParticle->GetParticleName(); << 138    return fgPrimaryParticleName ;
123 }                                                 139 }
124                                                   140 
125 //....oooOO0OOooo........oooOO0OOooo........oo << 141 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
126                                                   142 
127 void F02PrimaryGeneratorAction::SetXVertex(G4d    143 void F02PrimaryGeneratorAction::SetXVertex(G4double x)
128 {                                                 144 {
129   fVertexDefined = true;                       << 145   fVertexDefined = true ;
130   fXVertex = x;                                << 146   fXVertex = x ;
131   G4cout << " X coordinate of the primary vert << 147   G4cout << " X coordinate of the primary vertex = " << fXVertex/mm <<
                                                   >> 148             " mm." << G4endl;
132 }                                                 149 }
133                                                   150 
134 //....oooOO0OOooo........oooOO0OOooo........oo << 151 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
135                                                   152 
136 void F02PrimaryGeneratorAction::SetYVertex(G4d    153 void F02PrimaryGeneratorAction::SetYVertex(G4double y)
137 {                                                 154 {
138   fVertexDefined = true;                       << 155   fVertexDefined = true ;
139   fYVertex = y;                                << 156   fYVertex = y ;
140   G4cout << " Y coordinate of the primary vert << 157   G4cout << " Y coordinate of the primary vertex = " << fYVertex/mm <<
                                                   >> 158             " mm." << G4endl;
141 }                                                 159 }
142                                                   160 
143 //....oooOO0OOooo........oooOO0OOooo........oo << 161 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
144                                                   162 
145 void F02PrimaryGeneratorAction::SetZVertex(G4d    163 void F02PrimaryGeneratorAction::SetZVertex(G4double z)
146 {                                                 164 {
147   fVertexDefined = true;                       << 165   fVertexDefined = true ;
148   fZVertex = z;                                << 166   fZVertex = z ;
149   G4cout << " Z coordinate of the primary vert << 167   G4cout << " Z coordinate of the primary vertex = " << fZVertex/mm <<
                                                   >> 168             " mm." << G4endl;
150 }                                                 169 }
151                                                   170 
152 //....oooOO0OOooo........oooOO0OOooo........oo << 171 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
153                                                   172