Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/digits_hits/scorer/src/G4PSMinKinEAtGeneration.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/G4PSMinKinEAtGeneration.cc (Version 11.3.0) and /digits_hits/scorer/src/G4PSMinKinEAtGeneration.cc (Version 11.1.1)


  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 // G4PSMinKinEAtGeneration                         28 // G4PSMinKinEAtGeneration
 29 #include "G4PSMinKinEAtGeneration.hh"              29 #include "G4PSMinKinEAtGeneration.hh"
 30 #include "G4UnitsTable.hh"                         30 #include "G4UnitsTable.hh"
 31 #include "G4VScoreHistFiller.hh"                   31 #include "G4VScoreHistFiller.hh"
 32                                                    32 
 33 // (Description)                                   33 // (Description)
 34 //   This is a primitive scorer class for scor     34 //   This is a primitive scorer class for scoring minimum energy of
 35 //  newly produced particles in the geometry.      35 //  newly produced particles in the geometry.
 36 //                                                 36 //
 37 // Created: 2005-11-17  Tsukasa ASO, Akinori K     37 // Created: 2005-11-17  Tsukasa ASO, Akinori Kimura.
 38 // 2010-07-22   Introduce Unit specification.      38 // 2010-07-22   Introduce Unit specification.
 39 // 2011-09-09   Modify comment in PrintAll().      39 // 2011-09-09   Modify comment in PrintAll().  T.Aso
 40 // 2020-10-06   Use G4VPrimitivePlotter and fi     40 // 2020-10-06   Use G4VPrimitivePlotter and fill 1-D histo of kinetic energy (x)
 41 //              vs. track weight (y)               41 //              vs. track weight (y)                 (Makoto Asai)
 42 //                                                 42 //
 43                                                    43 
 44 G4PSMinKinEAtGeneration::G4PSMinKinEAtGenerati <<  44 G4PSMinKinEAtGeneration::G4PSMinKinEAtGeneration(G4String name, G4int depth)
 45   : G4PSMinKinEAtGeneration(name, "MeV", depth     45   : G4PSMinKinEAtGeneration(name, "MeV", depth) 
 46 {}                                                 46 {}
 47                                                    47 
 48 G4PSMinKinEAtGeneration::G4PSMinKinEAtGenerati <<  48 G4PSMinKinEAtGeneration::G4PSMinKinEAtGeneration(G4String name,
 49                                                    49                                                  const G4String& unit,
 50                                                    50                                                  G4int depth)
 51   : G4VPrimitivePlotter(name, depth)               51   : G4VPrimitivePlotter(name, depth)
 52   , HCID(-1)                                       52   , HCID(-1)
 53   , EvtMap(nullptr)                                53   , EvtMap(nullptr)
 54 {                                                  54 {
 55   SetUnit(unit);                                   55   SetUnit(unit);
 56 }                                                  56 }
 57                                                    57 
 58 G4bool G4PSMinKinEAtGeneration::ProcessHits(G4     58 G4bool G4PSMinKinEAtGeneration::ProcessHits(G4Step* aStep, G4TouchableHistory*)
 59 {                                                  59 {
 60   // ===============                               60   // ===============
 61   //  First Step,                                  61   //  First Step,
 62   //  Confirm this is a newly produced seconda     62   //  Confirm this is a newly produced secondary.
 63   //                                               63   //
 64   //- check for newly produced particle. e.g.      64   //- check for newly produced particle. e.g. Step number is 1.
 65   if(aStep->GetTrack()->GetCurrentStepNumber()     65   if(aStep->GetTrack()->GetCurrentStepNumber() != 1)
 66     return false;                                  66     return false;
 67   //- check for this is not a primary particle     67   //- check for this is not a primary particle. e.g. ParentID != 0 .
 68   if(aStep->GetTrack()->GetParentID() == 0)        68   if(aStep->GetTrack()->GetParentID() == 0)
 69     return false;                                  69     return false;
 70                                                    70 
 71   //==========================================     71   //===============================================
 72   //- This is a newly produced secondary parti     72   //- This is a newly produced secondary particle.
 73   //==========================================     73   //===============================================
 74                                                    74 
 75   // ===============                               75   // ===============
 76   //  Second Step,                                 76   //  Second Step,
 77   //  Confirm this track has lower energy than     77   //  Confirm this track has lower energy than previous one.
 78   //                                               78   //
 79                                                    79 
 80   // -Kinetic energy of this particle at the s     80   // -Kinetic energy of this particle at the starting point.
 81   G4int index      = GetIndex(aStep);              81   G4int index      = GetIndex(aStep);
 82   G4double kinetic = aStep->GetPreStepPoint()-     82   G4double kinetic = aStep->GetPreStepPoint()->GetKineticEnergy();
 83                                                    83 
 84   if(!hitIDMap.empty() && hitIDMap.find(index)     84   if(!hitIDMap.empty() && hitIDMap.find(index) != hitIDMap.cend())
 85   {                                                85   {
 86     auto filler = G4VScoreHistFiller::Instance     86     auto filler = G4VScoreHistFiller::Instance();
 87     if(filler == nullptr)                          87     if(filler == nullptr)
 88     {                                              88     {
 89       G4Exception(                                 89       G4Exception(
 90         "G4PSMinKinEAtGeneration::ProcessHits"     90         "G4PSMinKinEAtGeneration::ProcessHits", "SCORER0123", JustWarning,
 91         "G4TScoreHistFiller is not instantiate     91         "G4TScoreHistFiller is not instantiated!! Histogram is not filled.");
 92     }                                              92     }
 93     else                                           93     else
 94     {                                              94     {
 95       filler->FillH1(hitIDMap[index], kinetic,     95       filler->FillH1(hitIDMap[index], kinetic,
 96                      aStep->GetPreStepPoint()-     96                      aStep->GetPreStepPoint()->GetWeight());
 97     }                                              97     }
 98   }                                                98   }
 99                                                    99 
100   // -Stored value in the current HitsMap.        100   // -Stored value in the current HitsMap.
101   G4double* mapValue = ((*EvtMap)[index]);        101   G4double* mapValue = ((*EvtMap)[index]);
102                                                   102 
103   // -                                            103   // -
104   // If mapValue exits (e.g not NULL ), compar    104   // If mapValue exits (e.g not NULL ), compare it with
105   // current track's kinetic energy.              105   // current track's kinetic energy.
106   if((mapValue != nullptr) && (kinetic > *mapV    106   if((mapValue != nullptr) && (kinetic > *mapValue))
107     return false;                                 107     return false;
108                                                   108 
109   // -                                            109   // -
110   // Current Track is a newly produced seconda    110   // Current Track is a newly produced secondary and has lower
111   // kinetic energy than previous one in this     111   // kinetic energy than previous one in this geometry.
112   //                                              112   //
113   EvtMap->set(index, kinetic);                    113   EvtMap->set(index, kinetic);
114   return true;                                    114   return true;
115 }                                                 115 }
116                                                   116 
117 void G4PSMinKinEAtGeneration::Initialize(G4HCo    117 void G4PSMinKinEAtGeneration::Initialize(G4HCofThisEvent* HCE)
118 {                                                 118 {
119   EvtMap = new G4THitsMap<G4double>(detector->    119   EvtMap = new G4THitsMap<G4double>(detector->GetName(), GetName());
120   if(HCID < 0)                                    120   if(HCID < 0)
121   {                                               121   {
122     HCID = GetCollectionID(0);                    122     HCID = GetCollectionID(0);
123   }                                               123   }
124   HCE->AddHitsCollection(HCID, (G4VHitsCollect    124   HCE->AddHitsCollection(HCID, (G4VHitsCollection*) EvtMap);
125 }                                                 125 }
126                                                   126 
127 void G4PSMinKinEAtGeneration::clear() { EvtMap    127 void G4PSMinKinEAtGeneration::clear() { EvtMap->clear(); }
128                                                   128 
129 void G4PSMinKinEAtGeneration::PrintAll()          129 void G4PSMinKinEAtGeneration::PrintAll()
130 {                                                 130 {
131   G4cout << " PrimitiveScorer " << GetName() <    131   G4cout << " PrimitiveScorer " << GetName() << G4endl;
132   G4cout << " Number of entries " << EvtMap->e    132   G4cout << " Number of entries " << EvtMap->entries() << G4endl;
133   for(const auto& [copy, energy] : *(EvtMap->G    133   for(const auto& [copy, energy] : *(EvtMap->GetMap()))
134   {                                               134   {
135     G4cout << "  copy no.: " << copy              135     G4cout << "  copy no.: " << copy
136            << "  energy: " << *(energy) / GetU    136            << "  energy: " << *(energy) / GetUnitValue() << " ["
137            << GetUnit() << "]" << G4endl;         137            << GetUnit() << "]" << G4endl;
138   }                                               138   }
139 }                                                 139 }
140                                                   140 
141 void G4PSMinKinEAtGeneration::SetUnit(const G4    141 void G4PSMinKinEAtGeneration::SetUnit(const G4String& unit)
142 {                                                 142 {
143   CheckAndSetUnit(unit, "Energy");                143   CheckAndSetUnit(unit, "Energy");
144 }                                                 144 }
145                                                   145