source: trunk/examples/extended/electromagnetic/TestEm10/src/Em10CalorimeterSD.cc @ 1281

Last change on this file since 1281 was 1230, checked in by garnier, 14 years ago

update to geant4.9.3

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