ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
UltraPMTSD.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file UltraPMTSD.cc
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 //
27 // --------------------------------------------------------------
28 // GEANT 4 - ULTRA experiment example
29 // --------------------------------------------------------------
30 //
31 // Code developed by:
32 // B. Tome, M.C. Espirito-Santo, A. Trindade, P. Rodrigues
33 //
34 // ****************************************************
35 // * UltraPMTSD.cc
36 // ****************************************************
37 //
38 // Class used to define the Ultra photomultiplier as a sensitive detector.
39 // Hits in this sensitive detector are defined in the UltraOpticalHit class
40 //
41 #include "UltraPMTSD.hh"
42 
43 #include "G4Material.hh"
44 #include "G4Step.hh"
45 #include "G4VTouchable.hh"
46 #include "G4TouchableHistory.hh"
47 #include "G4SDManager.hh"
48 #include "G4OpticalPhoton.hh"
49 #include "G4ParticleDefinition.hh"
50 
51 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
52 
53 UltraPMTSD::UltraPMTSD(G4String name):G4VSensitiveDetector(name)
54 {
55 
56  collectionName.insert("OpticalHitsCollection");
57 
58 }
59 
60 
61 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
62 
63 UltraPMTSD::~UltraPMTSD()
64 {;}
65 
66 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
67 
68 void UltraPMTSD::Initialize(G4HCofThisEvent* HCE)
69 {
70 
71 
72  static int HCID1 = -1;
73 
74 
75  // SensitiveDetectorName and collectionName are data members of G4VSensitiveDetector
76 
77  OpticalHitsCollection =
78  new UltraOpticalHitsCollection(SensitiveDetectorName,collectionName[0]);
79 
80  if(HCID1<0)
81  { HCID1 = GetCollectionID(0); }
82  HCE->AddHitsCollection(HCID1,OpticalHitsCollection);
83 
84 }
85 
86 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
87 
88 G4bool UltraPMTSD::ProcessHits(G4Step* aStep,G4TouchableHistory*)
89 {
90 
91  // Get Material
92 
93  G4String thisVolume = aStep->GetTrack()->GetVolume()->GetName() ;
94  const G4ParticleDefinition* particle = aStep->GetTrack()->GetDefinition();
95 
96 
97  if (thisVolume != "PMT1" && thisVolume != "PMT2")
98  return false;
99  if (particle != G4OpticalPhoton::Definition() )
100  return false;
101 
102  if(particle == G4OpticalPhoton::Definition())
103  aStep->GetTrack()->SetTrackStatus(fStopAndKill);
104 
105  G4double kineticEnergy = aStep->GetTrack()->GetKineticEnergy();
106  G4ThreeVector HitPosition = aStep->GetPreStepPoint()->GetPosition() ;
107 
108  UltraOpticalHit* OpticalHit = new UltraOpticalHit ;
109  OpticalHit->SetEnergy(kineticEnergy);
110  OpticalHit->SetPosition(HitPosition);
111 
112 
113  OpticalHitsCollection->insert(OpticalHit);
114 
115 
116 #ifdef ULTRA_VERBOSE
117  G4cout << "*******************************" << G4endl;
118  G4cout << " PMT HIT " << G4endl;
119  G4cout << " Volume: " << thisVolume << G4endl;
120  G4cout << " Photon energy (eV) : " << kineticEnergy/CLHEP::eV << G4endl;
121  G4cout << " POSITION (mm) : "
122  << HitPosition.x()/CLHEP::mm << " " << HitPosition.y()/CLHEP::mm << " " << HitPosition.z()/CLHEP::mm << G4endl;
123  G4cout << "*******************************" << G4endl;
124 #endif
125 
126 
127  return true;
128 }
129 
130 
131 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
132 
133 void UltraPMTSD::EndOfEvent(G4HCofThisEvent* HCE)
134 {
135  static G4int HCID = -1;
136  if(HCID<0)
137  {
138  HCID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
139  }
140  HCE->AddHitsCollection(HCID,OpticalHitsCollection);
141 
142 }
143 
144 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
145 
146 void UltraPMTSD::clear()
147 {;}
148 
149 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
150 
151 void UltraPMTSD::DrawAll()
152 {;}
153 
154 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
155 
156 void UltraPMTSD::PrintAll()
157 {;}
158 
159 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....