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 // G4PSNofCollision 29 #include "G4PSNofCollision.hh" 30 31 ////////////////////////////////////////////// 32 // (Description) 33 // This is a primitive scorer class for scor 34 // Cell. 35 // 36 // Created: 2007-02-02 Tsukasa ASO, Akinori K 37 // 2010-07-22 Introduce Unit specification. 38 // 39 ////////////////////////////////////////////// 40 41 G4PSNofCollision::G4PSNofCollision(const G4Str 42 : G4VPrimitiveScorer(name, depth) 43 { 44 SetUnit(""); 45 } 46 47 G4bool G4PSNofCollision::ProcessHits(G4Step* a 48 { 49 if(aStep->GetPostStepPoint()->GetStepStatus( 50 return true; 51 52 G4int index = GetIndex(aStep); 53 G4double val = 1.0; 54 if(weighted) 55 val *= aStep->GetPreStepPoint()->GetWeight 56 EvtMap->add(index, val); 57 return true; 58 } 59 60 void G4PSNofCollision::Initialize(G4HCofThisEv 61 { 62 EvtMap = new G4THitsMap<G4double>(detector-> 63 if(HCID < 0) 64 { 65 HCID = GetCollectionID(0); 66 } 67 HCE->AddHitsCollection(HCID, (G4VHitsCollect 68 } 69 70 void G4PSNofCollision::clear() { EvtMap->clear 71 72 void G4PSNofCollision::PrintAll() 73 { 74 G4cout << " MultiFunctionalDet " << detecto 75 G4cout << " PrimitiveScorer " << GetName() < 76 G4cout << " Number of entries " << EvtMap->e 77 for(const auto& [copy, collisions] : *(EvtMa 78 { 79 G4cout << " copy no.: " << copy 80 << " collisions: " << *(collisions 81 << " [collision] " << G4endl; 82 } 83 } 84 85 void G4PSNofCollision::SetUnit(const G4String& 86 { 87 if(unit.empty()) 88 { 89 unitName = unit; 90 unitValue = 1.0; 91 } 92 else 93 { 94 G4String msg = "Invalid unit [" + unit + " 95 GetUnit() + "] ) for " + Ge 96 G4Exception("G4PSNofCollision::SetUnit", " 97 } 98 } 99