Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/HGCal_testbeam/src/SiPMHit.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  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 "SiPMHit.hh"
 27 
 28 #include <cstdlib>
 29 
 30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 31 
 32 SiPMHit::SiPMHit(G4String aVolumeName, G4int aCopyNumSensor,
 33                  G4int aCopyNumCell)
 34     : fVolumeName(aVolumeName), fCopyNumCell(aCopyNumCell),
 35       fCopyNumSensor(aCopyNumSensor) {}
 36 
 37 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 38 
 39 void SiPMHit::Digitise(const G4double aTimeWindow,
 40                        const G4double aToaThreshold) {
 41 
 42   // process energy deposits
 43   if (fEdep.empty()) {
 44     fIsValidHit = false;
 45     return;
 46   }
 47 
 48   std::sort(fEdep.begin(), fEdep.end(),
 49             [](const std::pair<G4double, G4double> &left,
 50                const std::pair<G4double, G4double> &right) {
 51               return left.second < right.second; // second = time
 52             });
 53 
 54   G4double firstHitTime = fEdep[0].second;
 55   fEdepDigi = 0;
 56   for (size_t i = 0; i < fEdep.size(); ++i) {
 57     if (aTimeWindow < 0 || fEdep[i].second < firstHitTime + aTimeWindow)
 58       fEdepDigi += fEdep[i].first;
 59 #ifdef DEBUG
 60     else
 61       std::cout << "Rejecting hit (" << fEdep[i].first << " keV)  at time "
 62                 << fEdep[i].second << "(start: " << firstHitTime << ")"
 63                 << std::endl;
 64 #endif
 65     if ((fTimeOfArrival == -1) && (fEdepDigi > aToaThreshold))
 66       fTimeOfArrival = fEdep[i].second;
 67   }
 68   fIsValidHit = (fEdepDigi > 0);
 69 
 70   // non ionizing part, does not contribute to TOAs
 71   if (fEdepNonIonizing.empty())
 72     return;
 73 
 74   std::sort(fEdepNonIonizing.begin(), fEdepNonIonizing.end(),
 75             [](const std::pair<G4double, G4double> &left,
 76                const std::pair<G4double, G4double> &right) {
 77               return left.second < right.second; // second = time
 78             });
 79 
 80   fEdepNonIonizingDigi = 0;
 81   for (size_t i = 0; i < fEdepNonIonizing.size(); ++i) {
 82     if (aTimeWindow == -1 ||
 83         fEdepNonIonizing[i].second < firstHitTime + aTimeWindow)
 84       fEdepNonIonizingDigi += fEdepNonIonizing[i].first;
 85   }
 86 }
 87