Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/electromagnetic/TestEm10/src/DetectorHarris73.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/electromagnetic/TestEm10/src/DetectorHarris73.cc (Version 11.3.0) and /examples/extended/electromagnetic/TestEm10/src/DetectorHarris73.cc (Version 10.3.p2)


  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 electromagnetic/TestEm10/src/Detecto     26 /// \file electromagnetic/TestEm10/src/DetectorHarris73.cc
 27 /// \brief Implementation of the DetectorHarri     27 /// \brief Implementation of the DetectorHarris73 class
 28 //                                                 28 //
 29 //                                                 29 //
                                                   >>  30 // $Id: DetectorHarris73.cc 72517 2013-07-24 23:01:29Z gum $
 30 //                                                 31 //
 31 //                                                 32 //
 32                                                    33 
 33 #include "DetectorHarris73.hh"                     34 #include "DetectorHarris73.hh"
 34                                                << 
 35 #include "Materials.hh"                        << 
 36 #include "SensitiveDetector.hh"                    35 #include "SensitiveDetector.hh"
                                                   >>  36 #include "Materials.hh"
 37                                                    37 
                                                   >>  38 #include "G4Material.hh"
 38 #include "G4Box.hh"                                39 #include "G4Box.hh"
 39 #include "G4FieldManager.hh"                   << 
 40 #include "G4LogicalVolume.hh"                      40 #include "G4LogicalVolume.hh"
 41 #include "G4Material.hh"                       << 
 42 #include "G4PVPlacement.hh"                        41 #include "G4PVPlacement.hh"
 43 #include "G4Region.hh"                         << 
 44 #include "G4SDManager.hh"                      << 
 45 #include "G4SystemOfUnits.hh"                  << 
 46 #include "G4TransportationManager.hh"          << 
 47 #include "G4UniformMagField.hh"                    42 #include "G4UniformMagField.hh"
                                                   >>  43 #include "G4FieldManager.hh"
                                                   >>  44 #include "G4TransportationManager.hh"
                                                   >>  45 #include "G4SDManager.hh"
                                                   >>  46 
                                                   >>  47 #include "G4Region.hh"
                                                   >>  48 
 48 #include "G4UnitsTable.hh"                         49 #include "G4UnitsTable.hh"
                                                   >>  50 #include "G4SystemOfUnits.hh"
 49 #include "G4ios.hh"                                51 #include "G4ios.hh"
 50                                                    52 
 51 //....oooOO0OOooo........oooOO0OOooo........oo     53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 52                                                    54 
 53 DetectorHarris73::DetectorHarris73() : fRadiat <<  55 DetectorHarris73::DetectorHarris73()
                                                   >>  56   : fRadiatorDescription(0) 
                                                   >>  57 {}
 54                                                    58 
 55 //....oooOO0OOooo........oooOO0OOooo........oo     59 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 56                                                    60 
 57 DetectorHarris73::~DetectorHarris73()              61 DetectorHarris73::~DetectorHarris73()
 58 {                                                  62 {
 59   // delete fRadiatorDescription;                  63   // delete fRadiatorDescription;
 60   // the description is deleted in detector co <<  64         // the description is deleted in detector construction
 61 }                                                  65 }
 62                                                    66 
 63 //....oooOO0OOooo........oooOO0OOooo........oo     67 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 64                                                    68 
 65 G4VPhysicalVolume* DetectorHarris73::Construct     69 G4VPhysicalVolume* DetectorHarris73::Construct()
 66 {                                                  70 {
 67   // Geometry parameters                           71   // Geometry parameters
 68   //                                               72   //
 69                                                    73 
 70   G4cout << "DetectorHarris73 setup" << G4endl     74   G4cout << "DetectorHarris73 setup" << G4endl;
 71                                                    75 
 72   G4double worldSizeZ = 400. * cm;             <<  76   G4double worldSizeZ = 400.*cm;
 73   G4double worldSizeR = 20. * cm;              <<  77   G4double worldSizeR = 20.*cm;
 74                                                    78 
 75   // Radiator and detector parameters              79   // Radiator and detector parameters
 76                                                    80 
 77   G4double radThickness = 0.0127 * mm;         <<  81   G4double radThickness = 0.0127*mm;
 78   G4double gasGap = 0.762 * mm;                <<  82   G4double gasGap       = 0.762*mm;
 79   G4double foilGasRatio = radThickness / (radT <<  83   G4double foilGasRatio = radThickness/(radThickness+gasGap);
 80   G4double foilNumber = 100;                   <<  84   G4double foilNumber   = 100;
 81                                                <<  85 
 82   G4double absorberThickness = 15.0 * mm;      <<  86   G4double absorberThickness = 15.0*mm;
 83   G4double absorberRadius = 100. * mm;         <<  87   G4double absorberRadius   = 100.*mm;
 84                                                <<  88 
 85   G4double windowThick = 51.0 * micrometer;    <<  89   G4double windowThick    = 51.0*micrometer;
 86   G4double electrodeThick = 10.0 * micrometer; <<  90   G4double electrodeThick = 10.0*micrometer;
 87   G4double gapThick = 10.0 * cm;               <<  91   G4double gapThick       =  10.0*cm;
 88   G4double detGap = 0.01 * mm;                 <<  92   G4double detGap         =   0.01*mm;
 89                                                    93 
 90   G4double startZ = 100.0 * mm;                <<  94   G4double startZ       = 100.0*mm;
 91                                                    95 
 92   // Materials                                     96   // Materials
 93   //                                               97   //
 94                                                    98 
 95   // Change to create materials using NIST         99   // Change to create materials using NIST
 96   G4Material* air = Materials::GetInstance()-> << 100   G4Material* air   = Materials::GetInstance()->GetMaterial("Air");
 97   G4Material* mylar = Materials::GetInstance()    101   G4Material* mylar = Materials::GetInstance()->GetMaterial("Mylar");
 98   G4Material* kr7ch4 = Materials::GetInstance(    102   G4Material* kr7ch4 = Materials::GetInstance()->GetMaterial("Kr7CH4");
 99                                                   103 
100   // Preparation of mixed radiator material    << 104  // Preparation of mixed radiator material
101                                                   105 
102   G4double foilDensity = mylar->GetDensity();     106   G4double foilDensity = mylar->GetDensity();
103   G4double gasDensity = air->GetDensity();     << 107   G4double gasDensity  = air->GetDensity();  
104   G4double totDensity = foilDensity * foilGasR << 108   G4double totDensity  = foilDensity*foilGasRatio 
                                                   >> 109                        + gasDensity*(1.0-foilGasRatio);
105                                                   110 
106   G4double fractionFoil = foilDensity * foilGa << 111   G4double fractionFoil =  foilDensity*foilGasRatio/totDensity;
107   G4double fractionGas = gasDensity * (1.0 - f << 112   G4double fractionGas  =  gasDensity*(1.0-foilGasRatio)/totDensity;
108   G4Material* radiatorMat = new G4Material("ra    113   G4Material* radiatorMat = new G4Material("radiatorMat", totDensity, 2);
109   radiatorMat->AddMaterial(mylar, fractionFoil    114   radiatorMat->AddMaterial(mylar, fractionFoil);
110   radiatorMat->AddMaterial(air, fractionGas);     115   radiatorMat->AddMaterial(air, fractionGas);
111                                                   116 
112   // Radiator description                         117   // Radiator description
113   fRadiatorDescription = new RadiatorDescripti    118   fRadiatorDescription = new RadiatorDescription;
114   fRadiatorDescription->fFoilMaterial = mylar; << 119   fRadiatorDescription->fFoilMaterial  = mylar; // CH2; // Kapton; // Mylar ; // Li ; // CH2 ;
115   fRadiatorDescription->fGasMaterial = air;  / << 120   fRadiatorDescription->fGasMaterial   = air; // CO2; // He; //
116   fRadiatorDescription->fFoilThickness = radTh    121   fRadiatorDescription->fFoilThickness = radThickness;
117   fRadiatorDescription->fGasThickness = gasGap << 122   fRadiatorDescription->fGasThickness  = gasGap;
118   fRadiatorDescription->fFoilNumber = foilNumb    123   fRadiatorDescription->fFoilNumber = foilNumber;
119                                                   124 
120   G4Material* worldMaterial = air;  // CO2;    << 125   G4Material* worldMaterial    = air; // CO2;
121   G4Material* absorberMaterial = kr7ch4;          126   G4Material* absorberMaterial = kr7ch4;
122                                                   127 
123   // Volumes                                      128   // Volumes
124   //                                              129   //
                                                   >> 130  
                                                   >> 131   G4VSolid* solidWorld 
                                                   >> 132     = new G4Box("World", worldSizeR, worldSizeR, worldSizeZ/2.);
                                                   >> 133  
                                                   >> 134   G4LogicalVolume* logicWorld 
                                                   >> 135     = new G4LogicalVolume(solidWorld,  worldMaterial,  "World");
125                                                   136 
126   G4VSolid* solidWorld = new G4Box("World", wo << 137   G4VPhysicalVolume* physicsWorld 
127                                                << 138     = new G4PVPlacement(0, G4ThreeVector(), "World", logicWorld, 0,  false, 0);
128   G4LogicalVolume* logicWorld = new G4LogicalV << 
129                                                << 
130   G4VPhysicalVolume* physicsWorld =            << 
131     new G4PVPlacement(0, G4ThreeVector(), "Wor << 
132                                                   139 
133   // TR radiator envelope                         140   // TR radiator envelope
134                                                   141 
135   G4double radThick = foilNumber * (radThickne << 142   G4double radThick = foilNumber*(radThickness + gasGap) - gasGap + detGap;
136   G4double radZ = startZ + 0.5 * radThick;     << 143   G4double radZ = startZ + 0.5*radThick;
137                                                << 
138   G4VSolid* solidRadiator =                    << 
139     new G4Box("Radiator", 1.1 * absorberRadius << 
140                                                   144 
141   G4LogicalVolume* logicRadiator = new G4Logic << 145   G4VSolid* solidRadiator 
                                                   >> 146     = new G4Box("Radiator", 1.1*absorberRadius, 1.1*absorberRadius, 0.5*radThick);
142                                                   147 
                                                   >> 148   G4LogicalVolume* logicRadiator 
                                                   >> 149     = new G4LogicalVolume(solidRadiator, radiatorMat, "Radiator");
                                                   >> 150  
143   fRadiatorDescription->fLogicalVolume = logic    151   fRadiatorDescription->fLogicalVolume = logicRadiator;
144                                                << 152  
145   new G4PVPlacement(0, G4ThreeVector(0, 0, rad << 153   new G4PVPlacement(0, G4ThreeVector(0, 0, radZ),
146                     0);                        << 154                     "Radiator", logicRadiator, physicsWorld, false, 0 );
147                                                   155 
148   // create region for radiator                   156   // create region for radiator
149                                                   157 
150   G4Region* radRegion = new G4Region("XTRradia    158   G4Region* radRegion = new G4Region("XTRradiator");
151   radRegion->AddRootLogicalVolume(logicRadiato    159   radRegion->AddRootLogicalVolume(logicRadiator);
152                                                   160 
153   G4double windowZ = startZ + radThick + windo << 161   G4double windowZ = startZ + radThick + windowThick/2. + 15.0*mm;
154                                                   162 
155   G4double gapZ = windowZ + windowThick / 2. + << 163   G4double gapZ = windowZ + windowThick/2. + gapThick/2. + 0.01*mm;
156                                                   164 
157   G4double electrodeZ = gapZ + gapThick / 2. + << 165   G4double electrodeZ = gapZ + gapThick/2. + electrodeThick/2. + 0.01*mm;
158                                                   166 
159   // Absorber                                     167   // Absorber
160                                                   168 
161   G4double absorberZ = electrodeZ + electrodeT << 169   G4double absorberZ = electrodeZ + electrodeThick/2. 
                                                   >> 170                      + absorberThickness/2. + 0.01*mm;
162                                                   171 
163   G4VSolid* solidAbsorber =                    << 172   G4VSolid* solidAbsorber 
164     new G4Box("Absorber", absorberRadius, abso << 173     = new G4Box("Absorber", absorberRadius, absorberRadius, absorberThickness/2.);
165                                                   174 
166   G4LogicalVolume* logicAbsorber = new G4Logic << 175   G4LogicalVolume* logicAbsorber 
                                                   >> 176     = new G4LogicalVolume(solidAbsorber, absorberMaterial, "Absorber");
167                                                   177 
168   new G4PVPlacement(0, G4ThreeVector(0., 0., a << 178   new G4PVPlacement(0, G4ThreeVector(0., 0., absorberZ),
169                     false, 0);                 << 179                     "Absorber", logicAbsorber, physicsWorld, false, 0);
170                                                   180 
171   G4Region* regGasDet = new G4Region("XTRdEdxD    181   G4Region* regGasDet = new G4Region("XTRdEdxDetector");
172   regGasDet->AddRootLogicalVolume(logicAbsorbe    182   regGasDet->AddRootLogicalVolume(logicAbsorber);
173                                                   183 
174   // Sensitive Detectors: Absorber                184   // Sensitive Detectors: Absorber
175                                                   185 
176   SensitiveDetector* sd = new SensitiveDetecto    186   SensitiveDetector* sd = new SensitiveDetector("AbsorberSD");
177   G4SDManager::GetSDMpointer()->AddNewDetector << 187   G4SDManager::GetSDMpointer()->AddNewDetector(sd );
178   logicAbsorber->SetSensitiveDetector(sd);        188   logicAbsorber->SetSensitiveDetector(sd);
179                                                   189 
180   // Print geometry parameters                    190   // Print geometry parameters
181                                                   191 
182   G4cout << "\n The  WORLD   is made of " << w << 192   G4cout << "\n The  WORLD   is made of "
183          << worldMaterial->GetName();          << 193          << worldSizeZ/mm << "mm of " << worldMaterial->GetName();
184   G4cout << ", the transverse size (R) of the  << 194   G4cout << ", the transverse size (R) of the world is " 
185   G4cout << " The ABSORBER is made of " << abs << 195          << worldSizeR/mm << " mm. " << G4endl;
186          << absorberMaterial->GetName();       << 196   G4cout << " The ABSORBER is made of "
187   G4cout << ", the transverse size (R) is " << << 197          << absorberThickness/mm << "mm of " << absorberMaterial->GetName();
188   G4cout << " Z position of the (middle of the << 198   G4cout << ", the transverse size (R) is " 
                                                   >> 199          << absorberRadius/mm << " mm. " << G4endl;
                                                   >> 200   G4cout << " Z position of the (middle of the) absorber " 
                                                   >> 201          << absorberZ/mm << "  mm." << G4endl;
189                                                   202 
190   G4cout << "radZ = " << radZ / mm << " mm" << << 203   G4cout << "radZ = " << radZ/mm << " mm" << G4endl;
191   G4cout << "startZ = " << startZ / mm << " mm << 204   G4cout << "startZ = " << startZ/mm<< " mm" << G4endl;
192                                                   205 
193   G4cout << "fRadThick = " << radThick / mm << << 206   G4cout << "fRadThick = " << radThick/mm << " mm"<<G4endl;
194   G4cout << "fFoilNumber = " << foilNumber <<     207   G4cout << "fFoilNumber = " << foilNumber << G4endl;
195   G4cout << "fRadiatorMat = " << radiatorMat->    208   G4cout << "fRadiatorMat = " << radiatorMat->GetName() << G4endl;
196   G4cout << "WorldMaterial = " << worldMateria    209   G4cout << "WorldMaterial = " << worldMaterial->GetName() << G4endl;
197   G4cout << G4endl;                               210   G4cout << G4endl;
198                                                   211 
199   return physicsWorld;                            212   return physicsWorld;
200 }                                                 213 }
201                                                   214 
202 //....oooOO0OOooo........oooOO0OOooo........oo    215 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >> 216 
203                                                   217