| Geant4 Cross Reference |
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 /// \file eventgenerator/HepMC/HepMCEx01/src/ExN04CalorimeterROGeometry.cc
27 /// \brief Implementation of the ExN04CalorimeterROGeometry class
28 //
29 //
30
31 #include "ExN04CalorimeterROGeometry.hh"
32
33 #include "ExN04DummySD.hh"
34
35 #include "G4Box.hh"
36 #include "G4LogicalVolume.hh"
37 #include "G4Material.hh"
38 #include "G4PVPlacement.hh"
39 #include "G4PVReplica.hh"
40 #include "G4SDManager.hh"
41 #include "G4SystemOfUnits.hh"
42 #include "G4ThreeVector.hh"
43 #include "G4Tubs.hh"
44 #include "G4VPhysicalVolume.hh"
45
46 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
47 ExN04CalorimeterROGeometry::ExN04CalorimeterROGeometry() : G4VReadOutGeometry()
48 {
49 #include "ExN04DetectorParameterDef.icc"
50 }
51
52 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
53 ExN04CalorimeterROGeometry::ExN04CalorimeterROGeometry(G4String aString)
54 : G4VReadOutGeometry(aString)
55 {
56 #include "ExN04DetectorParameterDef.icc"
57 }
58
59 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
60 ExN04CalorimeterROGeometry::~ExN04CalorimeterROGeometry() {}
61
62 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
63 G4VPhysicalVolume* ExN04CalorimeterROGeometry::Build()
64 {
65 // A dummy material is used to fill the volumes of the readout geometry.
66 // ( It will be allowed to set a NULL pointer in volumes of such virtual
67 // division in future, since this material is irrelevant for tracking.)
68 G4Material* dummyMat = new G4Material(name = "dummyMat", 1., 1. * g / mole, 1. * g / cm3);
69
70 // Builds the ReadOut World:
71 G4Box* ROWorldBox = new G4Box("ROWorldBox", fexpHall_x, fexpHall_y, fexpHall_z);
72 G4LogicalVolume* ROWorldLog =
73 new G4LogicalVolume(ROWorldBox, dummyMat, "ROWorldLogical", 0, 0, 0);
74 G4PVPlacement* ROWorldPhys =
75 new G4PVPlacement(0, G4ThreeVector(), "ROWorldPhysical", ROWorldLog, 0, false, 0);
76 // Calorimeter volume:
77 G4VSolid* caloROtub = new G4Tubs("caloROtub", fcaloTubs_rmin, fcaloTubs_rmax, fcaloTubs_dz,
78 fcaloTubs_sphi, fcaloTubs_dphi);
79 G4LogicalVolume* caloROlog = new G4LogicalVolume(caloROtub, dummyMat, "caloROlogical", 0, 0, 0);
80 G4VPhysicalVolume* caloROphys =
81 new G4PVPlacement(0, G4ThreeVector(), "calROphysical", caloROlog, ROWorldPhys, false, 0);
82
83 // -------------------------------
84 // Calorimeter readout division:
85 // -------------------------------
86 // Phi division first: 48 sectors
87 G4VSolid* caloROphiDivisionTub = new G4Tubs("caloROphiDivision", fcaloCell_rmin, fcaloCell_rmax,
88 fcaloCell_dz, fcaloCell_sphi, fcaloCell_dphi);
89 G4LogicalVolume* caloROphiDivisionLog =
90 new G4LogicalVolume(caloROphiDivisionTub, dummyMat, "caloROphiDivisionLogical", 0, 0, 0);
91 G4VPhysicalVolume* caloROphiDivisionPhys =
92 new G4PVReplica("caloROphiDivisionPhysical", caloROphiDivisionLog, caloROphys, kPhi,
93 fsegmentsinPhi, fcaloCell_dphi);
94 // then z division: 20 slices:
95 G4VSolid* caloROcellTub = new G4Tubs("caloROcellTub", fcaloRing_rmin, fcaloRing_rmax,
96 fcaloRing_dz, fcaloRing_sphi, fcaloRing_dphi);
97 G4LogicalVolume* caloROcellLog =
98 new G4LogicalVolume(caloROcellTub, dummyMat, "caloROcellLogical", 0, 0, 0);
99 // G4VPhysicalVolume * caloROcellPhys =
100 new G4PVReplica("caloROcellPhysical", caloROcellLog, caloROphiDivisionPhys, kZAxis, fsegmentsinZ,
101 2. * fcaloRing_dz);
102
103 // Flags the cells as sensitive .The pointer here serves
104 // as a flag only to check for sensitivity.
105 // (Could we make it by a simple cast of a non-NULL value ?)
106 ExN04DummySD* dummySensi = new ExN04DummySD;
107 caloROcellLog->SetSensitiveDetector(dummySensi);
108
109 return ROWorldPhys;
110 }
111