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