Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/digits_hits/scorer/src/G4PSPassageCellCurrent.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /digits_hits/scorer/src/G4PSPassageCellCurrent.cc (Version 11.3.0) and /digits_hits/scorer/src/G4PSPassageCellCurrent.cc (Version 11.1.3)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26 //                                                 26 //
 27 //                                                 27 //
 28 // G4PSPassageCellCurrent                          28 // G4PSPassageCellCurrent
 29 #include "G4PSPassageCellCurrent.hh"               29 #include "G4PSPassageCellCurrent.hh"
 30 #include "G4StepStatus.hh"                         30 #include "G4StepStatus.hh"
 31 #include "G4Track.hh"                              31 #include "G4Track.hh"
 32 #include "G4VSolid.hh"                             32 #include "G4VSolid.hh"
 33 #include "G4UnitsTable.hh"                         33 #include "G4UnitsTable.hh"
 34 #include "G4VScoreHistFiller.hh"                   34 #include "G4VScoreHistFiller.hh"
 35                                                    35 
 36 //////////////////////////////////////////////     36 ////////////////////////////////////////////////////////////////////////////////
 37 // (Description)                                   37 // (Description)
 38 //   This is a primitive scorer class for scor     38 //   This is a primitive scorer class for scoring number of tracks,
 39 //   where only tracks passing through the geo     39 //   where only tracks passing through the geometry are taken
 40 //  into account.                                  40 //  into account.
 41 //                                                 41 //
 42 // Created: 2005-11-14  Tsukasa ASO, Akinori K     42 // Created: 2005-11-14  Tsukasa ASO, Akinori Kimura.
 43 // 2010-07-22   Introduce Unit specification.      43 // 2010-07-22   Introduce Unit specification.
 44 // 2010-07-22   Add weighted option                44 // 2010-07-22   Add weighted option
 45 // 2020-10-06   Use G4VPrimitivePlotter and fi     45 // 2020-10-06   Use G4VPrimitivePlotter and fill 1-D histo of kinetic energy (x)
 46 //              vs. current * track weight (y)     46 //              vs. current * track weight (y)               (Makoto Asai)
 47 //                                                 47 //
 48 //////////////////////////////////////////////     48 ///////////////////////////////////////////////////////////////////////////////
 49                                                    49 
 50 G4PSPassageCellCurrent::G4PSPassageCellCurrent <<  50 G4PSPassageCellCurrent::G4PSPassageCellCurrent(G4String name, G4int depth)
 51   : G4VPrimitivePlotter(name, depth)               51   : G4VPrimitivePlotter(name, depth)
 52 {                                                  52 {
 53   SetUnit("");                                     53   SetUnit("");
 54 }                                                  54 }
 55                                                    55 
 56 G4bool G4PSPassageCellCurrent::ProcessHits(G4S     56 G4bool G4PSPassageCellCurrent::ProcessHits(G4Step* aStep, G4TouchableHistory*)
 57 {                                                  57 {
 58   if(IsPassed(aStep))                              58   if(IsPassed(aStep))
 59   {                                                59   {
 60     fCurrent = 1.;                                 60     fCurrent = 1.;
 61     if(weighted)                                   61     if(weighted)
 62       fCurrent = aStep->GetPreStepPoint()->Get     62       fCurrent = aStep->GetPreStepPoint()->GetWeight();
 63     G4int index = GetIndex(aStep);                 63     G4int index = GetIndex(aStep);
 64     EvtMap->add(index, fCurrent);                  64     EvtMap->add(index, fCurrent);
 65                                                    65 
 66     if(!hitIDMap.empty() && hitIDMap.find(inde     66     if(!hitIDMap.empty() && hitIDMap.find(index) != hitIDMap.cend())
 67     {                                              67     {
 68       auto filler = G4VScoreHistFiller::Instan     68       auto filler = G4VScoreHistFiller::Instance();
 69       if(filler == nullptr)                        69       if(filler == nullptr)
 70       {                                            70       {
 71         G4Exception(                               71         G4Exception(
 72           "G4PSVolumeFlux::ProcessHits", "SCOR     72           "G4PSVolumeFlux::ProcessHits", "SCORER0123", JustWarning,
 73           "G4TScoreHistFiller is not instantia     73           "G4TScoreHistFiller is not instantiated!! Histogram is not filled.");
 74       }                                            74       }
 75       else                                         75       else
 76       {                                            76       {
 77         filler->FillH1(hitIDMap[index],            77         filler->FillH1(hitIDMap[index],
 78                        aStep->GetPreStepPoint(     78                        aStep->GetPreStepPoint()->GetKineticEnergy(), fCurrent);
 79       }                                            79       }
 80     }                                              80     }
 81   }                                                81   }
 82                                                    82 
 83   return true;                                     83   return true;
 84 }                                                  84 }
 85                                                    85 
 86 G4bool G4PSPassageCellCurrent::IsPassed(G4Step     86 G4bool G4PSPassageCellCurrent::IsPassed(G4Step* aStep)
 87 {                                                  87 {
 88   G4bool Passed = false;                           88   G4bool Passed = false;
 89                                                    89 
 90   G4bool IsEnter = aStep->GetPreStepPoint()->G     90   G4bool IsEnter = aStep->GetPreStepPoint()->GetStepStatus() == fGeomBoundary;
 91   G4bool IsExit  = aStep->GetPostStepPoint()->     91   G4bool IsExit  = aStep->GetPostStepPoint()->GetStepStatus() == fGeomBoundary;
 92                                                    92 
 93   G4int trkid = aStep->GetTrack()->GetTrackID(     93   G4int trkid = aStep->GetTrack()->GetTrackID();
 94                                                    94 
 95   if(IsEnter && IsExit)                            95   if(IsEnter && IsExit)
 96   {  // Passed at one step                         96   {  // Passed at one step
 97     Passed = true;                                 97     Passed = true;
 98   }                                                98   }
 99   else if(IsEnter)                                 99   else if(IsEnter)
100   {                         // Enter a new geo    100   {                         // Enter a new geometry
101     fCurrentTrkID = trkid;  // Resetting the c    101     fCurrentTrkID = trkid;  // Resetting the current track.
102   }                                               102   }
103   else if(IsExit)                                 103   else if(IsExit)
104   {  // Exit a current geometry                   104   {  // Exit a current geometry
105     if(fCurrentTrkID == trkid)                    105     if(fCurrentTrkID == trkid)
106     {                                             106     {
107       Passed = true;  // if the track is same     107       Passed = true;  // if the track is same as entered.
108     }                                             108     }
109   }                                               109   }
110   else                                            110   else
111   {  // Inside geometry                           111   {  // Inside geometry
112     if(fCurrentTrkID == trkid)                    112     if(fCurrentTrkID == trkid)
113     {  // Adding the track length to current o    113     {  // Adding the track length to current one ,
114     }                                             114     }
115   }                                               115   }
116   return Passed;                                  116   return Passed;
117 }                                                 117 }
118                                                   118 
119 void G4PSPassageCellCurrent::Initialize(G4HCof    119 void G4PSPassageCellCurrent::Initialize(G4HCofThisEvent* HCE)
120 {                                                 120 {
121   fCurrentTrkID = -1;                             121   fCurrentTrkID = -1;
122                                                   122 
123   EvtMap = new G4THitsMap<G4double>(detector->    123   EvtMap = new G4THitsMap<G4double>(detector->GetName(), GetName());
124   if(HCID < 0)                                    124   if(HCID < 0)
125     HCID = GetCollectionID(0);                    125     HCID = GetCollectionID(0);
126   HCE->AddHitsCollection(HCID, EvtMap);           126   HCE->AddHitsCollection(HCID, EvtMap);
127 }                                                 127 }
128                                                   128 
129 void G4PSPassageCellCurrent::clear() { EvtMap-    129 void G4PSPassageCellCurrent::clear() { EvtMap->clear(); }
130                                                   130 
131 void G4PSPassageCellCurrent::PrintAll()           131 void G4PSPassageCellCurrent::PrintAll()
132 {                                                 132 {
133   G4cout << " MultiFunctionalDet  " << detecto    133   G4cout << " MultiFunctionalDet  " << detector->GetName() << G4endl;
134   G4cout << " PrimitiveScorer " << GetName() <    134   G4cout << " PrimitiveScorer " << GetName() << G4endl;
135   G4cout << " Number of entries " << EvtMap->e    135   G4cout << " Number of entries " << EvtMap->entries() << G4endl;
136   for(const auto& [copy, current] : *(EvtMap->    136   for(const auto& [copy, current] : *(EvtMap->GetMap()))
137   {                                               137   {
138     G4cout << "  copy no.: " << copy              138     G4cout << "  copy no.: " << copy
139            << "  cell current : " << *(current    139            << "  cell current : " << *(current) << " [tracks] " << G4endl;
140   }                                               140   }
141 }                                                 141 }
142                                                   142 
143 void G4PSPassageCellCurrent::SetUnit(const G4S    143 void G4PSPassageCellCurrent::SetUnit(const G4String& unit)
144 {                                                 144 {
145   if(unit.empty())                                145   if(unit.empty())
146   {                                               146   {
147     unitName  = unit;                             147     unitName  = unit;
148     unitValue = 1.0;                              148     unitValue = 1.0;
149   }                                               149   }
150   else                                            150   else
151   {                                               151   {
152     G4String msg = "Invalid unit [" + unit + "    152     G4String msg = "Invalid unit [" + unit + "] (Current  unit is [" +
153                    GetUnit() + "] ) for " + Ge    153                    GetUnit() + "] ) for " + GetName();
154     G4Exception("G4PSPassageCellCurrent::SetUn    154     G4Exception("G4PSPassageCellCurrent::SetUnit", "DetPS0012", JustWarning,
155                 msg);                             155                 msg);
156   }                                               156   }
157 }                                                 157 }
158                                                   158