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 #include "SiliconPixelHit.hh" 27 28 #include "G4RotationMatrix.hh" 29 #include "G4SubtractionSolid.hh" 30 #include "G4Transform3D.hh" 31 #include "G4VVisManager.hh" 32 #include "G4VisAttributes.hh" 33 #include "HGCalTBMaterials.hh" 34 35 #include <cstdlib> 36 37 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 38 39 SiliconPixelHit::SiliconPixelHit(G4String aVolumeName, G4int aCopyNumSensor, 40 G4int aCopyNumCell) 41 : fVolumeName(aVolumeName), fCopyNumCell(aCopyNumCell), 42 fCopyNumSensor(aCopyNumSensor) {} 43 44 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 45 46 void SiliconPixelHit::Digitise(const G4double aTimeWindow, 47 const G4double aToaThreshold) { 48 49 // process energy deposits 50 if (fEdep.empty()) { 51 fIsValidHit = false; 52 return; 53 } 54 55 std::sort(fEdep.begin(), fEdep.end(), 56 [](const std::pair<G4double, G4double> &left, 57 const std::pair<G4double, G4double> &right) { 58 return left.second < right.second; // second = time 59 }); 60 61 G4double firstHitTime = fEdep[0].second; 62 fEdepDigi = 0; 63 for (size_t i = 0; i < fEdep.size(); ++i) { 64 if (aTimeWindow < 0 || fEdep[i].second < firstHitTime + aTimeWindow) 65 fEdepDigi += fEdep[i].first; 66 if (fEdepDigi > aToaThreshold) { 67 if (fTimeOfArrival == -1) 68 fTimeOfArrival = fEdep[i].second; 69 fTimeOfArrivalLast = fEdep[i].second; 70 } 71 } 72 fIsValidHit = (fEdepDigi > 0); 73 74 // non ionizing part, does not contribute to TOAs 75 if (fEdepNonIonizing.empty()) 76 return; 77 78 std::sort(fEdepNonIonizing.begin(), fEdepNonIonizing.end(), 79 [](const std::pair<G4double, G4double> &left, 80 const std::pair<G4double, G4double> &right) { 81 return left.second < right.second; // second = time 82 }); 83 84 fEdepNonIonizingDigi = 0; 85 for (size_t i = 0; i < fEdepNonIonizing.size(); ++i) { 86 if (aTimeWindow == -1 || 87 fEdepNonIonizing[i].second < firstHitTime + aTimeWindow) 88 fEdepNonIonizingDigi += fEdepNonIonizing[i].first; 89 } 90 } 91 92 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 93 94 void SiliconPixelHit::Draw() { 95 G4VVisManager *pVVisManager = G4VVisManager::GetConcreteInstance(); 96 if (!(fEdepDigi > 0)) 97 return; 98 if (pVVisManager) { 99 if (!pVVisManager->FilterHit(*this)) 100 return; 101 G4Transform3D trans(G4RotationMatrix(), 102 G4ThreeVector(fPosX * CLHEP::cm, fPosY * CLHEP::cm, fPosZ * CLHEP::cm)); 103 G4VisAttributes attribs; 104 auto solid = HexagonSolid("dummy", 0.3 * CLHEP::mm, 0.6496345 * CLHEP::cm); 105 G4Colour colour(1, 0, 0); 106 attribs.SetColour(colour); 107 attribs.SetForceSolid(true); 108 pVVisManager->Draw(*solid, attribs, trans); 109 } 110 } 111 112 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 113