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 // ------------------------------------------------------------ 29 // GEANT 4 class implementation 30 // 31 // ---------------- GFlashHitMaker ---------------- 32 // 33 // Authors: E.Barberio & Joanna Weng 34 // ------------------------------------------------------------ 35 36 #include "GFlashHitMaker.hh" 37 38 #include "G4GFlashSpot.hh" 39 #include "G4Step.hh" 40 #include "G4StepPoint.hh" 41 #include "G4TouchableHandle.hh" 42 #include "G4TransportationManager.hh" 43 #include "G4VGFlashSensitiveDetector.hh" 44 #include "G4VSensitiveDetector.hh" 45 #include "G4ios.hh" 46 47 GFlashHitMaker::GFlashHitMaker() 48 { 49 fTouchableHandle = new G4TouchableHistory(); // talk to ?@@@ 50 fpNavigator = new G4Navigator(); 51 fNaviSetup = false; 52 fWorldWithSdName = ""; 53 fpSpotS = new G4Step(); 54 fpSpotP = new G4StepPoint(); 55 // N.B. Pre and Post step points are common. 56 fpSpotS->SetPreStepPoint(fpSpotP); 57 fpSpotS->SetPostStepPoint(fpSpotP); 58 } 59 60 GFlashHitMaker::~GFlashHitMaker() 61 { 62 delete fpNavigator; 63 delete fpSpotP; 64 fpSpotS->ResetPreStepPoint(); 65 fpSpotS->ResetPostStepPoint(); 66 delete fpSpotS; 67 } 68 69 void GFlashHitMaker::make(GFlashEnergySpot* aSpot, const G4FastTrack* aT) 70 { 71 // Locate the spot 72 if (!fNaviSetup) { 73 // Choose the world volume that contains the sensitive detector based on its name (empty name 74 // for mass geometry) 75 G4VPhysicalVolume* worldWithSD = nullptr; 76 if (fWorldWithSdName.empty()) { 77 worldWithSD = G4TransportationManager::GetTransportationManager() 78 ->GetNavigatorForTracking() 79 ->GetWorldVolume(); 80 } 81 else { 82 worldWithSD = 83 G4TransportationManager::GetTransportationManager()->GetParallelWorld(fWorldWithSdName); 84 } 85 fpNavigator->SetWorldVolume(worldWithSD); 86 fpNavigator->LocateGlobalPointAndUpdateTouchable(aSpot->GetPosition(), fTouchableHandle(), 87 false); 88 fNaviSetup = true; 89 } 90 else { 91 fpNavigator->LocateGlobalPointAndUpdateTouchable(aSpot->GetPosition(), fTouchableHandle()); 92 } 93 94 //-------------------------------------- 95 // Produce Hits 96 // call sensitive part: taken/adapted from the stepping: 97 // Send G4Step information to Hit/Dig if the volume is sensitive 98 //--------------G4TouchableHistory---------------------------------------- 99 100 G4VPhysicalVolume* pCurrentVolume = fTouchableHandle()->GetVolume(); 101 G4VSensitiveDetector* pSensitive; 102 if (pCurrentVolume != 0) { 103 pSensitive = pCurrentVolume->GetLogicalVolume()->GetSensitiveDetector(); 104 G4VGFlashSensitiveDetector* gflashSensitive = 105 dynamic_cast<G4VGFlashSensitiveDetector*>(pSensitive); 106 if (gflashSensitive) { 107 // set spot information: 108 G4GFlashSpot theSpot(aSpot, aT, fTouchableHandle); 109 gflashSensitive->Hit(&theSpot); 110 } 111 else if (pSensitive) { 112 fpSpotS->SetTotalEnergyDeposit(aSpot->GetEnergy()); 113 fpSpotS->SetTrack(const_cast<G4Track*>(aT->GetPrimaryTrack())); 114 fpSpotP->SetWeight(aT->GetPrimaryTrack()->GetWeight()); 115 fpSpotP->SetPosition(aSpot->GetPosition()); 116 fpSpotP->SetGlobalTime(aT->GetPrimaryTrack()->GetGlobalTime()); 117 fpSpotP->SetLocalTime(aT->GetPrimaryTrack()->GetLocalTime()); 118 fpSpotP->SetProperTime(aT->GetPrimaryTrack()->GetProperTime()); 119 fpSpotP->SetTouchableHandle(fTouchableHandle); 120 fpSpotP->SetProcessDefinedStep(fpProcess); 121 fpSpotP->SetStepStatus(fUserDefinedLimit); 122 pSensitive->Hit(fpSpotS); 123 } 124 } 125 else { 126 #ifdef GFLASH_DEBUG 127 G4cout << "GFlashHitMaker::Out of volume " << G4endl; 128 #endif 129 } 130 } 131