source: trunk/source/digits_hits/scorer/src/G4PSCellCharge.cc @ 1340

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

update ti head

File size: 4.2 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// $Id: G4PSCellCharge.cc,v 1.2 2010/07/22 07:23:45 taso Exp $
27// GEANT4 tag $Name:  $
28//
29// G4PSCellCharge
30#include "G4PSCellCharge.hh"
31#include "G4Track.hh"
32
33
34///////////////////////////////////////////////////////////////////////////////
35// (Description)
36//   This is a primitive scorer class for scoring cell charge.
37//   The Cell Charge is defined by  a sum of deposited charge inside the cell.
38//
39// Created: 2006-08-20  Tsukasa ASO
40// 2010-07-22   Introduce Unit specification.
41//
42///////////////////////////////////////////////////////////////////////////////
43
44G4PSCellCharge::G4PSCellCharge(G4String name, G4int depth)
45    :G4VPrimitiveScorer(name,depth),HCID(-1)
46{
47    SetUnit("e+");
48}
49
50G4PSCellCharge::G4PSCellCharge(G4String name, const G4String& unit, 
51                               G4int depth)
52    :G4VPrimitiveScorer(name,depth),HCID(-1)
53{
54    SetUnit(unit);
55}
56
57G4PSCellCharge::~G4PSCellCharge()
58{;}
59
60G4bool G4PSCellCharge::ProcessHits(G4Step* aStep,G4TouchableHistory*)
61{
62
63    // Enter or First step of primary.
64    if( aStep->GetPreStepPoint()->GetStepStatus() == fGeomBoundary
65        || ( aStep->GetTrack()->GetParentID() == 0 &&
66             aStep->GetTrack()->GetCurrentStepNumber() == 1 ) ){
67        G4double CellCharge = aStep->GetPreStepPoint()->GetCharge();
68        CellCharge *= aStep->GetPreStepPoint()->GetWeight();
69        G4int index = GetIndex(aStep);
70        EvtMap->add(index,CellCharge);
71    }
72
73    // Exit
74    if( aStep->GetPostStepPoint()->GetStepStatus() == fGeomBoundary){
75        G4double CellCharge = aStep->GetPreStepPoint()->GetCharge();
76        CellCharge *= aStep->GetPreStepPoint()->GetWeight();
77        G4int index = GetIndex(aStep);
78        CellCharge *= -1.0;
79        EvtMap->add(index,CellCharge);
80    }
81
82  return TRUE;
83}
84
85void G4PSCellCharge::Initialize(G4HCofThisEvent* HCE)
86{
87  EvtMap = new G4THitsMap<G4double>(detector->GetName(),
88                                    GetName());
89  if ( HCID < 0 ) HCID = GetCollectionID(0);
90  HCE->AddHitsCollection(HCID,EvtMap);
91}
92
93void G4PSCellCharge::EndOfEvent(G4HCofThisEvent*)
94{;}
95
96void G4PSCellCharge::clear(){
97  EvtMap->clear();
98}
99
100void G4PSCellCharge::DrawAll()
101{;}
102
103void G4PSCellCharge::PrintAll()
104{
105  G4cout << " MultiFunctionalDet  " << detector->GetName() << G4endl;
106  G4cout << " PrimitiveScorer " << GetName() <<G4endl; 
107  G4cout << " Number of entries " << EvtMap->entries() << G4endl;
108  std::map<G4int,G4double*>::iterator itr = EvtMap->GetMap()->begin();
109  for(; itr != EvtMap->GetMap()->end(); itr++) {
110    G4cout << "  copy no.: " << itr->first
111           << "  cell charge : " << *(itr->second)/GetUnitValue()
112           << " [" << GetUnit() << "]"
113           << G4endl;
114  }
115}
116
117void G4PSCellCharge::SetUnit(const G4String& unit)
118{
119    CheckAndSetUnit(unit,"Electric charge");
120}
121
122
Note: See TracBrowser for help on using the repository browser.