FDsimG4LensSD.cc
Go to the documentation of this file.
1 #include "FDsimG4LensSD.hh"
2 #include "G4RunManager.hh"
3 #include "G4Material.hh"
4 #include "G4Step.hh"
5 #include "G4VTouchable.hh"
6 #include "G4TouchableHistory.hh"
7 #include "G4SDManager.hh"
8 #include "G4ios.hh"
9 #include <sstream>
10 
11 using namespace TelescopeSimulatorLX ;
12 
13 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
14 FDsimG4LensSD::FDsimG4LensSD(G4String name):G4VSensitiveDetector(name)
15 {
16  collectionName.insert("LensHits");
17 
18  G4cerr << "Created " << name << " Sensitive Detector." << G4endl ;
19 }
20 
21 
22 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
23 
24 FDsimG4LensSD::~FDsimG4LensSD()
25 {;}
26 
27 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
28 
29 void FDsimG4LensSD::Initialize(G4HCofThisEvent* HCE)
30 {
31 
32  static int HCID1 = -1;
33 
34  // SensitiveDetectorName e collectionName sao data members de G4VSensitiveDetector
35 
36  LensHitsCollection =
37  new FDsimG4OpticalHitsCollection(SensitiveDetectorName,collectionName[0]);
38 
39  if(HCID1<0)
40  { HCID1 = GetCollectionID(0); }
41  HCE->AddHitsCollection(HCID1,LensHitsCollection);
42 
43 }
44 
45 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
46 
47 G4bool FDsimG4LensSD::ProcessHits(G4Step* aStep,G4TouchableHistory*){
48 
49  G4TouchableHistory* theTouchable
50  = (G4TouchableHistory*)(aStep->GetPreStepPoint()->GetTouchable());
51 
52 
53  // Get Material
54 
55  G4String thisVolume = theTouchable->GetVolume()->GetName() ;
56  G4String particleName = aStep->GetTrack()->GetDefinition()->GetParticleName();
57 
58 
59 #ifdef VERBOSE_USE
60  G4cerr << "Volume name: " << thisVolume << G4endl ;;
61 #endif
62 
63 
64  G4int HitID ;
65 
66  G4VPhysicalVolume *fPhysVolume = theTouchable->GetVolume() ;
67  G4String VolumeName = fPhysVolume->GetName() ;
68  G4int copyNo = fPhysVolume->GetCopyNo();
69 
70 
71 
72  if(particleName == "opticalphoton"){HitID=1;}else{HitID=10;}
73 
74  G4double fGlobalTOF = aStep->GetTrack()->GetGlobalTime() ;
75  G4double kineticEnergy = aStep->GetTrack()->GetKineticEnergy() ;
76  G4double Wavelength = h_Planck*c_light/kineticEnergy ;
77 
78 
79  // G4StepPoint* preStepPoint ; // unused (except in #def VERBOSE_USE)
80  G4StepPoint* posStepPoint ;
81  // preStepPoint = aStep->GetPreStepPoint(); // unused (except in #def VERBOSE_USE)
82  posStepPoint = aStep->GetPostStepPoint();
83 
84  G4ThreeVector worldPos = posStepPoint->GetPosition();
85  G4ThreeVector localPos = theTouchable->GetHistory()->GetTopTransform().TransformPoint(worldPos);
86 
87  G4ThreeVector HitPosition = localPos;
88 
89 
90  G4ThreeVector HitDirection = posStepPoint->GetMomentumDirection() ;
91  G4ThreeVector SDposition = fPhysVolume->GetObjectTranslation() ;
92 
93  if (HitDirection.z() < 0.){
94 
95  FDsimG4OpticalHit* Hit = new FDsimG4OpticalHit ;
96 
97  Hit->SetEnergy(kineticEnergy);
98  Hit->SetWavelength(Wavelength);
99  Hit->SetTime(fGlobalTOF);
100  Hit->SetPosition(HitPosition);
101  Hit->SetWorldPosition(worldPos);
102  Hit->SetDirection(HitDirection);
103  Hit->SetPixelPosition(SDposition);
104  Hit->SetHitID(HitID);
105 
106  Hit->SetPMTid(copyNo) ;
107 
108 
109  LensHitsCollection->insert(Hit);
110 
111  }
112 
113 #ifdef VERBOSE_USE
114  G4cerr << "Volume name: " << VolumeName << " "
115  << "Pre step volume " << aStep->GetPreStemPoint()->GetPhysicalVolume()->GetName() << " "
116  << "Pos step volume " << posStepPoint->GetPhysicalVolume()->GetName()<< G4endl ;
117 
118 
119  G4cerr << fPhysVolume->GetObjectTranslation() << " "
120  << fPhysVolume->GetFrameTranslation() << " "
121  << fPhysVolume->GetTranslation() << G4endl ;
122 
123  G4cerr << fPhysVolume->GetMultiplicity() << " "
124  << fPhysVolume->GetCopyNo() << G4endl ;
125 
126 #endif
127  return true ;
128 
129 }
130 
131 
132 
133 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
134 
135 void FDsimG4LensSD::EndOfEvent(G4HCofThisEvent* HCE)
136 {
137  static G4int HCID = -1;
138  if(HCID<0)
139  {
140  HCID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
141  }
142  HCE->AddHitsCollection(HCID,LensHitsCollection);
143 
144 }
145 
146 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
147 
148 void FDsimG4LensSD::clear()
149 {;}
150 
151 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
152 
153 void FDsimG4LensSD::DrawAll()
154 {;}
155 
156 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
157 
158 void FDsimG4LensSD::PrintAll()
159 {;}
160 
161 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
TelescopeSimulatorLX::FDsimG4OpticalHitsCollection
G4THitsCollection< FDsimG4OpticalHit > FDsimG4OpticalHitsCollection
Definition: FDsimG4OpticalHit.hh:62
TelescopeSimulatorLX::FDsimG4LensSD::Initialize
void Initialize(G4HCofThisEvent *)
Definition: FDsimG4LensSD.cc:29
TelescopeSimulatorLX::FDsimG4LensSD::ProcessHits
G4bool ProcessHits(G4Step *, G4TouchableHistory *)
Definition: FDsimG4LensSD.cc:47
TelescopeSimulatorLX::FDsimG4LensSD::LensHitsCollection
FDsimG4OpticalHitsCollection * LensHitsCollection
Definition: FDsimG4LensSD.hh:32
TelescopeSimulatorLX::FDsimG4LensSD::EndOfEvent
void EndOfEvent(G4HCofThisEvent *)
Definition: FDsimG4LensSD.cc:135
, generated on Tue Sep 26 2023.