Geant4 Cross Reference |
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 // 28 // G4PSNofSecondary 29 #include "G4PSNofSecondary.hh" 30 #include "G4VScoreHistFiller.hh" 31 32 // (Description) 33 // This is a primitive scorer class for scoring number of secondaries 34 // in the Cell. 35 // 36 // Created: 2005-11-14 Tsukasa ASO, Akinori Kimura. 37 // Modify: 2011-09-09 T.Aso modify comment in PrintAll(). 38 // 2020-10-06 Use G4VPrimitivePlotter and fill 1-D histo of kinetic 39 // energy of the secondary in MeV (x) vs. track weight (y) 40 // (Makoto Asai) 41 // 42 43 G4PSNofSecondary::G4PSNofSecondary(const G4String& name, G4int depth) 44 : G4VPrimitivePlotter(name, depth) 45 {} 46 47 G4bool G4PSNofSecondary::ProcessHits(G4Step* aStep, G4TouchableHistory*) 48 { 49 //- check for newly produced particle. e.g. first step. 50 if(aStep->GetTrack()->GetCurrentStepNumber() != 1) 51 return false; 52 //- check for this is not a primary particle. e.g. ParentID > 0 . 53 if(aStep->GetTrack()->GetParentID() == 0) 54 return false; 55 //- check the particle if the partifle definition is given. 56 if((particleDef != nullptr) && particleDef != aStep->GetTrack()->GetDefinition()) 57 return false; 58 // 59 //- This is a newly produced secondary particle. 60 G4int index = GetIndex(aStep); 61 G4double weight = 1.0; 62 if(weighted) 63 weight *= aStep->GetPreStepPoint()->GetWeight(); 64 EvtMap->add(index, weight); 65 66 if(!hitIDMap.empty() && hitIDMap.find(index) != hitIDMap.cend()) 67 { 68 auto filler = G4VScoreHistFiller::Instance(); 69 if(filler == nullptr) 70 { 71 G4Exception( 72 "G4PSVolumeFlux::ProcessHits", "SCORER0123", JustWarning, 73 "G4TScoreHistFiller is not instantiated!! Histogram is not filled."); 74 } 75 else 76 { 77 filler->FillH1(hitIDMap[index], 78 aStep->GetPreStepPoint()->GetKineticEnergy(), weight); 79 } 80 } 81 82 return true; 83 } 84 85 void G4PSNofSecondary::SetParticle(const G4String& particleName) 86 { 87 G4ParticleDefinition* pd = 88 G4ParticleTable::GetParticleTable()->FindParticle(particleName); 89 if(pd == nullptr) 90 { 91 G4String msg = "Particle <"; 92 msg += particleName; 93 msg += "> not found."; 94 G4Exception("G4PSNofSecondary::SetParticle", "DetPS0101", FatalException, 95 msg); 96 } 97 particleDef = pd; 98 } 99 100 void G4PSNofSecondary::Initialize(G4HCofThisEvent* HCE) 101 { 102 EvtMap = new G4THitsMap<G4double>(detector->GetName(), GetName()); 103 if(HCID < 0) 104 { 105 HCID = GetCollectionID(0); 106 } 107 HCE->AddHitsCollection(HCID, (G4VHitsCollection*) EvtMap); 108 } 109 110 void G4PSNofSecondary::clear() { EvtMap->clear(); } 111 112 void G4PSNofSecondary::PrintAll() 113 { 114 G4cout << " PrimitiveScorer " << GetName() << G4endl; 115 G4cout << " Number of entries " << EvtMap->entries() << G4endl; 116 for(const auto& [copy, secondaries] : *(EvtMap->GetMap())) 117 { 118 G4cout << " copy no.: " << copy 119 << " num of secondaries: " << *(secondaries) / GetUnitValue() 120 << G4endl; 121 } 122 } 123 124 void G4PSNofSecondary::SetUnit(const G4String& unit) 125 { 126 if(unit.empty()) 127 { 128 unitName = unit; 129 unitValue = 1.0; 130 } 131 else 132 { 133 G4String msg = "Invalid unit [" + unit + "] (Current unit is [" + 134 GetUnit() + "] ) for " + GetName(); 135 G4Exception("G4PSNofSecondary::SetUnit", "DetPS0010", JustWarning, msg); 136 } 137 } 138