source: trunk/examples/extended/field/field03/src/F03CalorimeterSD.cc@ 1036

Last change on this file since 1036 was 807, checked in by garnier, 17 years ago

update

File size: 4.6 KB
RevLine 
[807]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// $Id: F03CalorimeterSD.cc,v 1.6 2006/06/29 17:19:20 gunter Exp $
28// GEANT4 tag $Name: $
29//
30//
31
32//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
33//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
34
35#include "F03CalorimeterSD.hh"
36
37#include "F03CalorHit.hh"
38#include "F03DetectorConstruction.hh"
39
40#include "G4VPhysicalVolume.hh"
41#include "G4Step.hh"
42#include "G4VTouchable.hh"
43#include "G4TouchableHistory.hh"
44#include "G4SDManager.hh"
45
46#include "G4ios.hh"
47
48//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
49
50F03CalorimeterSD::F03CalorimeterSD(G4String name,
51 F03DetectorConstruction* det)
52:G4VSensitiveDetector(name),Detector(det)
53{
54 collectionName.insert("CalCollection");
55 HitID = new G4int[500];
56}
57
58//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
59
60F03CalorimeterSD::~F03CalorimeterSD()
61{
62 delete [] HitID;
63}
64
65//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
66
67void F03CalorimeterSD::Initialize(G4HCofThisEvent*)
68{
69 CalCollection = new F03CalorHitsCollection(SensitiveDetectorName,
70 collectionName[0]);
71 for (G4int j=0;j<1; j++)
72 {
73 HitID[j] = -1;
74 }
75}
76
77//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
78
79G4bool F03CalorimeterSD::ProcessHits(G4Step* aStep,G4TouchableHistory*)
80{
81 G4double edep = aStep->GetTotalEnergyDeposit();
82
83 G4double stepl = 0.;
84
85 stepl = aStep->GetStepLength();
86
87 if ((edep == 0.) && (stepl == 0.) ) return false;
88
89 G4TouchableHistory* theTouchable
90 = (G4TouchableHistory*)(aStep->GetPreStepPoint()->GetTouchable());
91
92 G4VPhysicalVolume* physVol = theTouchable->GetVolume();
93
94 G4int F03Number = 0 ;
95 if (HitID[F03Number]==-1)
96 {
97 F03CalorHit* calHit = new F03CalorHit();
98 if (physVol == Detector->GetAbsorber()) calHit->AddAbs(edep,stepl);
99 HitID[F03Number] = CalCollection->insert(calHit) - 1;
100 if (verboseLevel>0)
101 G4cout << " New Calorimeter Hit on F03: " << F03Number << G4endl;
102 }
103 else
104 {
105 if (physVol == Detector->GetAbsorber())
106 (*CalCollection)[HitID[F03Number]]->AddAbs(edep,stepl);
107 if (verboseLevel>0)
108 G4cout << " Energy added to F03: " << F03Number << G4endl;
109 }
110
111 return true;
112}
113
114//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
115
116void F03CalorimeterSD::EndOfEvent(G4HCofThisEvent* HCE)
117{
118 static G4int HCID = -1;
119 if(HCID<0)
120 { HCID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]); }
121 HCE->AddHitsCollection(HCID,CalCollection);
122}
123
124//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
125
126void F03CalorimeterSD::clear()
127{}
128
129//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
130
131
132void F03CalorimeterSD::PrintAll()
133{}
134
135//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
136
Note: See TracBrowser for help on using the repository browser.