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 // Class G4PSVolumeFlux 27 // 28 // Class description: 29 // Scorer that scores number of tracks coming 30 // Optionally number of tracks can be divided 31 // score the volume current and divided by co 32 // where theta is the incident angle 33 // 34 // - Created M. Asai, Sept. 2020 35 // 36 // 37 #include "G4PSVolumeFlux.hh" 38 39 #include "G4SystemOfUnits.hh" 40 #include "G4StepStatus.hh" 41 #include "G4Track.hh" 42 #include "G4VSolid.hh" 43 #include "G4VPhysicalVolume.hh" 44 #include "G4VPVParameterisation.hh" 45 #include "G4UnitsTable.hh" 46 #include "G4GeometryTolerance.hh" 47 #include "G4VScoreHistFiller.hh" 48 49 G4PSVolumeFlux::G4PSVolumeFlux(const G4String& 50 : G4VPrimitivePlotter(name, depth) 51 , fDirection(direction) 52 {} 53 54 G4bool G4PSVolumeFlux::ProcessHits(G4Step* aSt 55 { 56 G4StepPoint* preStepPoint = aStep->GetPreSt 57 G4StepPoint* postStepPoint = aStep->GetPostS 58 G4StepPoint* thisStepPoint = nullptr; 59 if(fDirection == 1) // Score only the inwar 60 { 61 if(preStepPoint->GetStepStatus() == fGeomB 62 { 63 thisStepPoint = preStepPoint; 64 } 65 else 66 { 67 return false; 68 } 69 } 70 else if(fDirection == 2) // Score only the 71 { 72 if(postStepPoint->GetStepStatus() == fGeom 73 { 74 thisStepPoint = postStepPoint; 75 } 76 else 77 { 78 return false; 79 } 80 } 81 82 G4double flux = preStepPoint->GetWeight(); 83 84 if(divare || divcos) 85 { 86 G4VPhysicalVolume* physVol = preStep 87 G4VPVParameterisation* physParam = physVol 88 G4VSolid* solid = nullptr 89 if(physParam != nullptr) 90 { // for parameterized volume 91 auto idx = ((G4TouchableHistory*) (preSt 92 ->GetReplicaNumber(indexDep 93 solid = physParam->ComputeSolid(idx, phy 94 solid->ComputeDimensions(physParam, idx, 95 } 96 else 97 { // for ordinary volume 98 solid = physVol->GetLogicalVolume()->Get 99 } 100 101 if(divare) 102 { 103 flux /= solid->GetSurfaceArea(); 104 } 105 106 if(divcos) 107 { 108 G4TouchableHandle theTouchable = thisSte 109 G4ThreeVector pdirection = thisSte 110 G4ThreeVector localdir = 111 theTouchable->GetHistory()->GetTopTran 112 G4ThreeVector globalPos = thisStepPoint- 113 G4ThreeVector localPos = 114 theTouchable->GetHistory()->GetTopTran 115 G4ThreeVector surfNormal = solid->Surfac 116 G4double cosT = surfNormal.co 117 if(cosT != 0.) 118 flux /= std::abs(cosT); 119 } 120 } 121 122 G4int index = GetIndex(aStep); 123 EvtMap->add(index, flux); 124 125 if(!hitIDMap.empty() && hitIDMap.find(index) 126 { 127 auto filler = G4VScoreHistFiller::Instance 128 if(filler == nullptr) 129 { 130 G4Exception( 131 "G4PSVolumeFlux::ProcessHits", "SCORER 132 "G4TScoreHistFiller is not instantiate 133 } 134 else 135 { 136 filler->FillH1(hitIDMap[index], thisStep 137 } 138 } 139 140 return true; 141 } 142 143 void G4PSVolumeFlux::Initialize(G4HCofThisEven 144 { 145 if(HCID < 0) 146 HCID = GetCollectionID(0); 147 EvtMap = new G4THitsMap<G4double>(GetMultiFu 148 GetName()) 149 HCE->AddHitsCollection(HCID, (G4VHitsCollect 150 } 151 152 void G4PSVolumeFlux::clear() { EvtMap->clear() 153 154 void G4PSVolumeFlux::PrintAll() 155 { 156 G4cout << " MultiFunctionalDet " << detecto 157 G4cout << " PrimitiveScorer" << GetName() << 158 G4cout << " Number of entries " << EvtMap->e 159 for(const auto& [copy, flux] : *(EvtMap->Get 160 { 161 G4cout << " copy no.: " << copy << " flu 162 << G4endl; 163 } 164 } 165