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 5.2.p2)


  1 //                                                  1 
  2 // *******************************************    
  3 // * License and Disclaimer                       
  4 // *                                              
  5 // * The  Geant4 software  is  copyright of th    
  6 // * the Geant4 Collaboration.  It is provided    
  7 // * conditions of the Geant4 Software License    
  8 // * LICENSE and available at  http://cern.ch/    
  9 // * include a list of copyright holders.         
 10 // *                                              
 11 // * Neither the authors of this software syst    
 12 // * institutes,nor the agencies providing fin    
 13 // * work  make  any representation or  warran    
 14 // * regarding  this  software system or assum    
 15 // * use.  Please see the license in the file     
 16 // * for the full disclaimer and the limitatio    
 17 // *                                              
 18 // * This  code  implementation is the result     
 19 // * technical work of the GEANT4 collaboratio    
 20 // * By using,  copying,  modifying or  distri    
 21 // * any work based  on the software)  you  ag    
 22 // * use  in  resulting  scientific  publicati    
 23 // * acceptance of all terms of the Geant4 Sof    
 24 // *******************************************    
 25 //                                                
 26 //                                                
 27 //                                                
 28 // G4PSMinKinEAtGeneration                        
 29 #include "G4PSMinKinEAtGeneration.hh"             
 30 #include "G4UnitsTable.hh"                        
 31 #include "G4VScoreHistFiller.hh"                  
 32                                                   
 33 // (Description)                                  
 34 //   This is a primitive scorer class for scor    
 35 //  newly produced particles in the geometry.     
 36 //                                                
 37 // Created: 2005-11-17  Tsukasa ASO, Akinori K    
 38 // 2010-07-22   Introduce Unit specification.     
 39 // 2011-09-09   Modify comment in PrintAll().     
 40 // 2020-10-06   Use G4VPrimitivePlotter and fi    
 41 //              vs. track weight (y)              
 42 //                                                
 43                                                   
 44 G4PSMinKinEAtGeneration::G4PSMinKinEAtGenerati    
 45   : G4PSMinKinEAtGeneration(name, "MeV", depth    
 46 {}                                                
 47                                                   
 48 G4PSMinKinEAtGeneration::G4PSMinKinEAtGenerati    
 49                                                   
 50                                                   
 51   : G4VPrimitivePlotter(name, depth)              
 52   , HCID(-1)                                      
 53   , EvtMap(nullptr)                               
 54 {                                                 
 55   SetUnit(unit);                                  
 56 }                                                 
 57                                                   
 58 G4bool G4PSMinKinEAtGeneration::ProcessHits(G4    
 59 {                                                 
 60   // ===============                              
 61   //  First Step,                                 
 62   //  Confirm this is a newly produced seconda    
 63   //                                              
 64   //- check for newly produced particle. e.g.     
 65   if(aStep->GetTrack()->GetCurrentStepNumber()    
 66     return false;                                 
 67   //- check for this is not a primary particle    
 68   if(aStep->GetTrack()->GetParentID() == 0)       
 69     return false;                                 
 70                                                   
 71   //==========================================    
 72   //- This is a newly produced secondary parti    
 73   //==========================================    
 74                                                   
 75   // ===============                              
 76   //  Second Step,                                
 77   //  Confirm this track has lower energy than    
 78   //                                              
 79                                                   
 80   // -Kinetic energy of this particle at the s    
 81   G4int index      = GetIndex(aStep);             
 82   G4double kinetic = aStep->GetPreStepPoint()-    
 83                                                   
 84   if(!hitIDMap.empty() && hitIDMap.find(index)    
 85   {                                               
 86     auto filler = G4VScoreHistFiller::Instance    
 87     if(filler == nullptr)                         
 88     {                                             
 89       G4Exception(                                
 90         "G4PSMinKinEAtGeneration::ProcessHits"    
 91         "G4TScoreHistFiller is not instantiate    
 92     }                                             
 93     else                                          
 94     {                                             
 95       filler->FillH1(hitIDMap[index], kinetic,    
 96                      aStep->GetPreStepPoint()-    
 97     }                                             
 98   }                                               
 99                                                   
100   // -Stored value in the current HitsMap.        
101   G4double* mapValue = ((*EvtMap)[index]);        
102                                                   
103   // -                                            
104   // If mapValue exits (e.g not NULL ), compar    
105   // current track's kinetic energy.              
106   if((mapValue != nullptr) && (kinetic > *mapV    
107     return false;                                 
108                                                   
109   // -                                            
110   // Current Track is a newly produced seconda    
111   // kinetic energy than previous one in this     
112   //                                              
113   EvtMap->set(index, kinetic);                    
114   return true;                                    
115 }                                                 
116                                                   
117 void G4PSMinKinEAtGeneration::Initialize(G4HCo    
118 {                                                 
119   EvtMap = new G4THitsMap<G4double>(detector->    
120   if(HCID < 0)                                    
121   {                                               
122     HCID = GetCollectionID(0);                    
123   }                                               
124   HCE->AddHitsCollection(HCID, (G4VHitsCollect    
125 }                                                 
126                                                   
127 void G4PSMinKinEAtGeneration::clear() { EvtMap    
128                                                   
129 void G4PSMinKinEAtGeneration::PrintAll()          
130 {                                                 
131   G4cout << " PrimitiveScorer " << GetName() <    
132   G4cout << " Number of entries " << EvtMap->e    
133   for(const auto& [copy, energy] : *(EvtMap->G    
134   {                                               
135     G4cout << "  copy no.: " << copy              
136            << "  energy: " << *(energy) / GetU    
137            << GetUnit() << "]" << G4endl;         
138   }                                               
139 }                                                 
140                                                   
141 void G4PSMinKinEAtGeneration::SetUnit(const G4    
142 {                                                 
143   CheckAndSetUnit(unit, "Energy");                
144 }                                                 
145