source: trunk/examples/extended/analysis/AnaEx01/src/AnaEx01CalorimeterSD.cc

Last change on this file was 1337, checked in by garnier, 14 years ago

tag geant4.9.4 beta 1 + modifs locales

File size: 5.1 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: AnaEx01CalorimeterSD.cc,v 1.5 2006/06/29 16:33:46 gunter Exp $
28// GEANT4 tag $Name: geant4-09-04-beta-01 $
29//
30//
31
32//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
33//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
34
35#include "AnaEx01CalorimeterSD.hh"
36
37#include "AnaEx01CalorHit.hh"
38#include "AnaEx01DetectorConstruction.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
50AnaEx01CalorimeterSD::AnaEx01CalorimeterSD(G4String name,
51                                   AnaEx01DetectorConstruction* det)
52:G4VSensitiveDetector(name),Detector(det)
53{
54  collectionName.insert("CalCollection");
55  HitID = new G4int[500];
56}
57
58//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
59
60AnaEx01CalorimeterSD::~AnaEx01CalorimeterSD()
61{
62  delete [] HitID;
63}
64
65//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
66
67void AnaEx01CalorimeterSD::Initialize(G4HCofThisEvent*)
68{
69  CalCollection = new AnaEx01CalorHitsCollection
70                      (SensitiveDetectorName,collectionName[0]); 
71  for (G4int j=0;j<Detector->GetNbOfLayers();j++) {HitID[j] = -1;};
72}
73
74//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
75
76G4bool AnaEx01CalorimeterSD::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 LayerNumber = 0;
92  if (Detector->GetNbOfLayers()>1) LayerNumber=theTouchable->GetReplicaNumber(1);
93
94  if (HitID[LayerNumber]==-1)
95    { 
96      AnaEx01CalorHit* calHit = new AnaEx01CalorHit();
97      if (physVol == Detector->GetAbsorber()) calHit->AddAbs(edep,stepl);
98      if (physVol == Detector->GetGap     ()) calHit->AddGap(edep,stepl);
99      HitID[LayerNumber] = CalCollection->insert(calHit) - 1;
100      if (verboseLevel>0)
101        G4cout << " New Calorimeter Hit on layer: " << LayerNumber << G4endl;
102    }
103  else
104    { 
105      if (physVol == Detector->GetAbsorber())
106         (*CalCollection)[HitID[LayerNumber]]->AddAbs(edep,stepl);
107      if (physVol == Detector->GetGap())
108         (*CalCollection)[HitID[LayerNumber]]->AddGap(edep,stepl);
109      if (verboseLevel>0)
110        G4cout << " Energy added to Layer: " << LayerNumber << G4endl; 
111    }
112   
113  return true;
114}
115
116//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
117
118void AnaEx01CalorimeterSD::EndOfEvent(G4HCofThisEvent* HCE)
119{
120  static G4int HCID = -1;
121  if(HCID<0)
122  { HCID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]); }
123  HCE->AddHitsCollection(HCID,CalCollection);
124}
125
126//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
127
128void AnaEx01CalorimeterSD::clear()
129{} 
130
131//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
132
133void AnaEx01CalorimeterSD::DrawAll()
134{} 
135
136//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
137
138void AnaEx01CalorimeterSD::PrintAll()
139{} 
140
141//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
142
Note: See TracBrowser for help on using the repository browser.