Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/runAndEvent/RE06/src/RE06ParallelWorld.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/runAndEvent/RE06/src/RE06ParallelWorld.cc (Version 11.3.0) and /examples/extended/runAndEvent/RE06/src/RE06ParallelWorld.cc (Version 10.2.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 /// \file RE06/src/RE06ParallelWorld.cc            26 /// \file RE06/src/RE06ParallelWorld.cc
 27 /// \brief Implementation of the RE06ParallelW     27 /// \brief Implementation of the RE06ParallelWorld class
 28 //                                                 28 //
 29 //                                             <<  29 // $Id: RE06ParallelWorld.cc 83381 2014-08-21 14:18:13Z gcosmo $
                                                   >>  30 // 
 30                                                    31 
 31 #include "RE06ParallelWorld.hh"                    32 #include "RE06ParallelWorld.hh"
 32                                                    33 
                                                   >>  34 #include "G4Tubs.hh"
 33 #include "G4LogicalVolume.hh"                      35 #include "G4LogicalVolume.hh"
 34 #include "G4MultiFunctionalDetector.hh"        << 
 35 #include "G4PSEnergyDeposit.hh"                << 
 36 #include "G4PSMinKinEAtGeneration.hh"          << 
 37 #include "G4PSNofSecondary.hh"                 << 
 38 #include "G4PSNofStep.hh"                      << 
 39 #include "G4PSTrackLength.hh"                  << 
 40 #include "G4PVPlacement.hh"                        36 #include "G4PVPlacement.hh"
 41 #include "G4PVReplica.hh"                          37 #include "G4PVReplica.hh"
                                                   >>  38 
 42 #include "G4SDManager.hh"                          39 #include "G4SDManager.hh"
 43 #include "G4SDParticleFilter.hh"               <<  40 #include "G4MultiFunctionalDetector.hh"
 44 #include "G4SystemOfUnits.hh"                  << 
 45 #include "G4Tubs.hh"                           << 
 46 #include "G4VPrimitiveScorer.hh"                   41 #include "G4VPrimitiveScorer.hh"
                                                   >>  42 #include "G4PSEnergyDeposit.hh"
                                                   >>  43 #include "G4PSNofSecondary.hh"
                                                   >>  44 #include "G4PSTrackLength.hh"
                                                   >>  45 #include "G4PSNofStep.hh"
                                                   >>  46 #include "G4PSMinKinEAtGeneration.hh"
 47 #include "G4VSDFilter.hh"                          47 #include "G4VSDFilter.hh"
                                                   >>  48 #include "G4SDParticleFilter.hh"
 48 #include "G4ios.hh"                                49 #include "G4ios.hh"
                                                   >>  50 #include "G4SystemOfUnits.hh"
 49                                                    51 
 50 //....oooOO0OOooo........oooOO0OOooo........oo     52 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 51                                                    53 
 52 G4ThreadLocal G4bool RE06ParallelWorld::fSDCon     54 G4ThreadLocal G4bool RE06ParallelWorld::fSDConstructed = false;
 53                                                    55 
 54 RE06ParallelWorld::RE06ParallelWorld(G4String      56 RE06ParallelWorld::RE06ParallelWorld(G4String worldName)
 55   : G4VUserParallelWorld(worldName),           <<  57 :G4VUserParallelWorld(worldName),
 56     fConstructed(false),                       <<  58  fConstructed(false),
 57     fSerial(false),                            <<  59  fSerial(false),
 58     fTotalThickness(2.0 * m),                  <<  60  fTotalThickness(2.0*m),
 59     fNumberOfLayers(20)                        <<  61  fNumberOfLayers(20)
 60 {                                                  62 {
 61   for (size_t i = 0; i < 3; i++) {             <<  63   for(size_t i=0;i<3;i++)
                                                   >>  64   {
 62     fCalorLogical[i] = 0;                          65     fCalorLogical[i] = 0;
 63     fLayerLogical[i] = 0;                          66     fLayerLogical[i] = 0;
 64     fCalorPhysical[i] = 0;                         67     fCalorPhysical[i] = 0;
 65     fLayerPhysical[i] = 0;                         68     fLayerPhysical[i] = 0;
 66   }                                                69   }
 67   fCalName[0] = "Calor-AP";                        70   fCalName[0] = "Calor-AP";
 68   fCalName[1] = "Calor-BP";                        71   fCalName[1] = "Calor-BP";
 69   fCalName[2] = "Calor-CP";                        72   fCalName[2] = "Calor-CP";
 70 }                                                  73 }
 71                                                    74 
 72 //....oooOO0OOooo........oooOO0OOooo........oo     75 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 73                                                    76 
 74 RE06ParallelWorld::~RE06ParallelWorld()            77 RE06ParallelWorld::~RE06ParallelWorld()
 75 {                                              <<  78 {;}
 76   ;                                            << 
 77 }                                              << 
 78                                                    79 
 79 //....oooOO0OOooo........oooOO0OOooo........oo     80 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 80                                                    81 
 81 void RE06ParallelWorld::Construct()                82 void RE06ParallelWorld::Construct()
 82 {                                                  83 {
 83   if (!fConstructed) {                         <<  84   if(!fConstructed)
                                                   >>  85   { 
 84     fConstructed = true;                           86     fConstructed = true;
 85     SetupGeometry();                               87     SetupGeometry();
 86   }                                                88   }
 87 }                                                  89 }
 88                                                    90 
 89 //....oooOO0OOooo........oooOO0OOooo........oo     91 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 90                                                    92 
 91 void RE06ParallelWorld::ConstructSD()              93 void RE06ParallelWorld::ConstructSD()
 92 {                                                  94 {
 93   if (!fSDConstructed) {                       <<  95   if(!fSDConstructed)
                                                   >>  96   {
 94     fSDConstructed = true;                         97     fSDConstructed = true;
 95     SetupDetectors();                              98     SetupDetectors();
 96   }                                                99   }
 97 }                                                 100 }
 98                                                   101 
 99 //....oooOO0OOooo........oooOO0OOooo........oo    102 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
100 void RE06ParallelWorld::SetupGeometry()           103 void RE06ParallelWorld::SetupGeometry()
101 {                                                 104 {
102   //                                           << 105   //     
103   // World                                        106   // World
104   //                                              107   //
105   G4VPhysicalVolume* ghostWorld = GetWorld();     108   G4VPhysicalVolume* ghostWorld = GetWorld();
106   G4LogicalVolume* worldLogical = ghostWorld->    109   G4LogicalVolume* worldLogical = ghostWorld->GetLogicalVolume();
107                                                << 110   
108   //                                           << 111   //                               
109   // Calorimeter                                  112   // Calorimeter
110   //                                           << 113   //  
111   G4VSolid* calorSolid = new G4Tubs("Calor", 0 << 114   G4VSolid* calorSolid 
                                                   >> 115     = new G4Tubs("Calor",0.0,0.5*m,fTotalThickness/2.,0.0,360.*deg);
112   G4int i;                                        116   G4int i;
113   for (i = 0; i < 3; i++) {                    << 117   for(i=0;i<3;i++)
114     fCalorLogical[i] = new G4LogicalVolume(cal << 118   {
115     if (fSerial) {                             << 119     fCalorLogical[i] = new G4LogicalVolume(calorSolid,0,fCalName[i]);
116       fCalorPhysical[i] =                      << 120     if(fSerial)
117         new G4PVPlacement(0, G4ThreeVector(0., << 121     {
118                           fCalorLogical[i], fC << 122       fCalorPhysical[i] = new G4PVPlacement(0,
                                                   >> 123                  G4ThreeVector(0.,0.,G4double(i-1)*fTotalThickness),
                                                   >> 124                  fCalorLogical[i],fCalName[i],worldLogical,false,i);
119     }                                             125     }
120     else {                                     << 126     else
121       fCalorPhysical[i] = new G4PVPlacement(0, << 127     {
122                                             fC << 128       fCalorPhysical[i] = new G4PVPlacement(0,
                                                   >> 129                  G4ThreeVector(0.,G4double(i-1)*m,0.),
                                                   >> 130                  fCalorLogical[i],fCalName[i],worldLogical,false,i);
123     }                                             131     }
124   }                                               132   }
125                                                << 133  
126   //                                           << 134   //                                 
127   // Layers --- as absorbers                      135   // Layers --- as absorbers
128   //                                              136   //
129   G4VSolid* layerSolid = new G4Tubs("Layer", 0 << 137   G4VSolid* layerSolid 
130   for (i = 0; i < 3; i++) {                    << 138     = new G4Tubs("Layer",0.0,0.5*m,fTotalThickness/2.,0.0,360.*deg);
131     fLayerLogical[i] = new G4LogicalVolume(lay << 139   for(i=0;i<3;i++)
132     fLayerPhysical[i] = new G4PVReplica(fCalNa << 140   {
133                                         kRho,  << 141     fLayerLogical[i] 
                                                   >> 142       = new G4LogicalVolume(layerSolid,0,fCalName[i]+"_LayerLog");
                                                   >> 143     fLayerPhysical[i] 
                                                   >> 144       = new G4PVReplica(fCalName[i]+"_Layer",fLayerLogical[i],fCalorLogical[i],
                                                   >> 145                         kRho,fNumberOfLayers,0.5*m/fNumberOfLayers);
134   }                                               146   }
                                                   >> 147    
135 }                                                 148 }
136                                                   149 
137 //....oooOO0OOooo........oooOO0OOooo........oo    150 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
138                                                   151 
139 void RE06ParallelWorld::SetupDetectors()          152 void RE06ParallelWorld::SetupDetectors()
140 {                                                 153 {
141   G4SDManager::GetSDMpointer()->SetVerboseLeve    154   G4SDManager::GetSDMpointer()->SetVerboseLevel(1);
142   G4String filterName, particleName;              155   G4String filterName, particleName;
143                                                   156 
144   G4SDParticleFilter* gammaFilter =            << 157   G4SDParticleFilter* gammaFilter 
145     new G4SDParticleFilter(filterName = "gamma << 158     = new G4SDParticleFilter(filterName="gammaFilter",particleName="gamma");
146   G4SDParticleFilter* electronFilter =         << 159   G4SDParticleFilter* electronFilter 
147     new G4SDParticleFilter(filterName = "elect << 160     = new G4SDParticleFilter(filterName="electronFilter",particleName="e-");
148   G4SDParticleFilter* positronFilter =         << 161   G4SDParticleFilter* positronFilter 
149     new G4SDParticleFilter(filterName = "posit << 162     = new G4SDParticleFilter(filterName="positronFilter",particleName="e+");
150   G4SDParticleFilter* epFilter = new G4SDParti << 163   G4SDParticleFilter* epFilter 
151   epFilter->add(particleName = "e-");          << 164     = new G4SDParticleFilter(filterName="epFilter");
152   epFilter->add(particleName = "e+");          << 165   epFilter->add(particleName="e-");
153                                                << 166   epFilter->add(particleName="e+");
154   for (G4int i = 0; i < 3; i++) {              << 167 
155     G4String detName = fCalName[i] + "_para";  << 168   for(G4int i=0;i<3;i++)
                                                   >> 169   {
                                                   >> 170     G4String detName = fCalName[i]+"_para";
156     G4MultiFunctionalDetector* det = new G4Mul    171     G4MultiFunctionalDetector* det = new G4MultiFunctionalDetector(detName);
157                                                   172 
158     G4VPrimitiveScorer* primitive;                173     G4VPrimitiveScorer* primitive;
159     primitive = new G4PSEnergyDeposit("eDep");    174     primitive = new G4PSEnergyDeposit("eDep");
160     det->RegisterPrimitive(primitive);            175     det->RegisterPrimitive(primitive);
161     primitive = new G4PSNofSecondary("nGamma")    176     primitive = new G4PSNofSecondary("nGamma");
162     primitive->SetFilter(gammaFilter);            177     primitive->SetFilter(gammaFilter);
163     det->RegisterPrimitive(primitive);            178     det->RegisterPrimitive(primitive);
164     primitive = new G4PSNofSecondary("nElectro    179     primitive = new G4PSNofSecondary("nElectron");
165     primitive->SetFilter(electronFilter);         180     primitive->SetFilter(electronFilter);
166     det->RegisterPrimitive(primitive);            181     det->RegisterPrimitive(primitive);
167     primitive = new G4PSNofSecondary("nPositro    182     primitive = new G4PSNofSecondary("nPositron");
168     primitive->SetFilter(positronFilter);         183     primitive->SetFilter(positronFilter);
169     det->RegisterPrimitive(primitive);            184     det->RegisterPrimitive(primitive);
170     primitive = new G4PSTrackLength("trackLeng    185     primitive = new G4PSTrackLength("trackLength");
171     primitive->SetFilter(epFilter);               186     primitive->SetFilter(epFilter);
172     det->RegisterPrimitive(primitive);            187     det->RegisterPrimitive(primitive);
173     primitive = new G4PSNofStep("nStep");         188     primitive = new G4PSNofStep("nStep");
174     primitive->SetFilter(epFilter);               189     primitive->SetFilter(epFilter);
175     det->RegisterPrimitive(primitive);            190     det->RegisterPrimitive(primitive);
176                                                   191 
177     G4SDManager::GetSDMpointer()->AddNewDetect    192     G4SDManager::GetSDMpointer()->AddNewDetector(det);
178     SetSensitiveDetector(fLayerLogical[i], det    193     SetSensitiveDetector(fLayerLogical[i], det);
179   }                                               194   }
180   G4SDManager::GetSDMpointer()->SetVerboseLeve    195   G4SDManager::GetSDMpointer()->SetVerboseLevel(0);
                                                   >> 196   
181 }                                                 197 }
182                                                   198 
183 //....oooOO0OOooo........oooOO0OOooo........oo    199 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
184                                                   200 
185 void RE06ParallelWorld::SetSerialGeometry(G4bo    201 void RE06ParallelWorld::SetSerialGeometry(G4bool serial)
186 {                                                 202 {
187   if (fSerial == serial) return;               << 203   if(fSerial==serial) return;
188   fSerial = serial;                            << 204   fSerial=serial;
189   if (!fConstructed) return;                   << 205   if(!fConstructed) return;
190   for (G4int i = 0; i < 3; i++) {              << 206   for(G4int i=0;i<3;i++)
191     if (fSerial) {                             << 207   {
192       fCalorPhysical[i]->SetTranslation(G4Thre << 208     if(fSerial)
                                                   >> 209     { 
                                                   >> 210       fCalorPhysical[i]
                                                   >> 211         ->SetTranslation(G4ThreeVector(0.,0.,G4double(i-1)*2.*m)); 
193     }                                             212     }
194     else {                                     << 213     else
195       fCalorPhysical[i]->SetTranslation(G4Thre << 214     { 
                                                   >> 215       fCalorPhysical[i]
                                                   >> 216         ->SetTranslation(G4ThreeVector(0.,G4double(i-1)*m,0.)); 
196     }                                             217     }
197   }                                               218   }
198 }                                                 219 }
199                                                   220 
200 //....oooOO0OOooo........oooOO0OOooo........oo    221 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
201                                                   222