Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/optical/LXe/src/LXePMTSD.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 ]

Diff markup

Differences between /examples/extended/optical/LXe/src/LXePMTSD.cc (Version 11.3.0) and /examples/extended/optical/LXe/src/LXePMTSD.cc (Version 10.6.p3)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26 //                                                 26 //
 27 /// \file optical/LXe/src/LXePMTSD.cc              27 /// \file optical/LXe/src/LXePMTSD.cc
 28 /// \brief Implementation of the LXePMTSD clas     28 /// \brief Implementation of the LXePMTSD class
 29 //                                                 29 //
 30 //                                                 30 //
 31 #include "LXePMTSD.hh"                             31 #include "LXePMTSD.hh"
 32                                                << 
 33 #include "LXeDetectorConstruction.hh"          << 
 34 #include "LXePMTHit.hh"                            32 #include "LXePMTHit.hh"
                                                   >>  33 #include "LXeDetectorConstruction.hh"
 35 #include "LXeUserTrackInformation.hh"              34 #include "LXeUserTrackInformation.hh"
 36                                                    35 
                                                   >>  36 #include "G4VPhysicalVolume.hh"
 37 #include "G4LogicalVolume.hh"                      37 #include "G4LogicalVolume.hh"
 38 #include "G4ParticleDefinition.hh"             << 
 39 #include "G4ParticleTypes.hh"                  << 
 40 #include "G4SDManager.hh"                      << 
 41 #include "G4Step.hh"                           << 
 42 #include "G4TouchableHistory.hh"               << 
 43 #include "G4Track.hh"                              38 #include "G4Track.hh"
 44 #include "G4VPhysicalVolume.hh"                <<  39 #include "G4Step.hh"
 45 #include "G4VTouchable.hh"                         40 #include "G4VTouchable.hh"
                                                   >>  41 #include "G4TouchableHistory.hh"
 46 #include "G4ios.hh"                                42 #include "G4ios.hh"
                                                   >>  43 #include "G4ParticleTypes.hh"
                                                   >>  44 #include "G4ParticleDefinition.hh"
 47                                                    45 
 48 //....oooOO0OOooo........oooOO0OOooo........oo     46 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 49                                                    47 
 50 LXePMTSD::LXePMTSD(G4String name) : G4VSensiti <<  48 LXePMTSD::LXePMTSD(G4String name)
                                                   >>  49   : G4VSensitiveDetector(name),fPMTHitCollection(nullptr),
                                                   >>  50     fPMTPositionsX(nullptr),fPMTPositionsY(nullptr),fPMTPositionsZ(nullptr)
 51 {                                                  51 {
 52   collectionName.insert("pmtHitCollection");       52   collectionName.insert("pmtHitCollection");
 53 }                                                  53 }
 54                                                    54 
 55 //....oooOO0OOooo........oooOO0OOooo........oo     55 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 56                                                    56 
 57 LXePMTSD::~LXePMTSD()                          <<  57 LXePMTSD::~LXePMTSD() {}
 58 {                                              << 
 59   delete fPMTPositionsX;                       << 
 60   delete fPMTPositionsY;                       << 
 61   delete fPMTPositionsZ;                       << 
 62 }                                              << 
 63                                                    58 
 64 //....oooOO0OOooo........oooOO0OOooo........oo     59 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 65                                                    60 
 66 void LXePMTSD::SetPmtPositions(const std::vect     61 void LXePMTSD::SetPmtPositions(const std::vector<G4ThreeVector>& positions)
 67 {                                                  62 {
 68   for (size_t i = 0; i < positions.size(); ++i <<  63   for (G4int i=0; i<G4int(positions.size()); ++i) {
 69     if (fPMTPositionsX) fPMTPositionsX->push_b <<  64     if(fPMTPositionsX)fPMTPositionsX->push_back(positions[i].x());
 70     if (fPMTPositionsY) fPMTPositionsY->push_b <<  65     if(fPMTPositionsY)fPMTPositionsY->push_back(positions[i].y());
 71     if (fPMTPositionsZ) fPMTPositionsZ->push_b <<  66     if(fPMTPositionsZ)fPMTPositionsZ->push_back(positions[i].z());
 72   }                                                67   }
 73 }                                                  68 }
 74                                                    69 
 75 //....oooOO0OOooo........oooOO0OOooo........oo     70 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 76                                                    71 
 77 void LXePMTSD::Initialize(G4HCofThisEvent* hit <<  72 void LXePMTSD::Initialize(G4HCofThisEvent* hitsCE){
 78 {                                              <<  73   fPMTHitCollection = new LXePMTHitsCollection
 79   fPMTHitCollection = new LXePMTHitsCollection <<  74                       (SensitiveDetectorName,collectionName[0]);
 80                                                <<  75   //Store collection with event and keep ID
 81   if (fHitCID < 0) {                           <<  76   static G4int hitCID = -1;
 82     fHitCID = G4SDManager::GetSDMpointer()->Ge <<  77   if(hitCID<0){
                                                   >>  78     hitCID = GetCollectionID(0);
 83   }                                                79   }
 84   hitsCE->AddHitsCollection(fHitCID, fPMTHitCo <<  80   hitsCE->AddHitsCollection( hitCID, fPMTHitCollection );
 85 }                                                  81 }
 86                                                    82 
 87 //....oooOO0OOooo........oooOO0OOooo........oo     83 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 88                                                    84 
 89 G4bool LXePMTSD::ProcessHits(G4Step*, G4Toucha <<  85 G4bool LXePMTSD::ProcessHits(G4Step* ,G4TouchableHistory* ){
 90 {                                              << 
 91   return false;                                    86   return false;
 92 }                                                  87 }
 93                                                    88 
 94 //....oooOO0OOooo........oooOO0OOooo........oo     89 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 95                                                    90 
 96 // Generates a hit and uses the postStepPoint' <<  91 //Generates a hit and uses the postStepPoint's mother volume replica number
 97 // PostStepPoint because the hit is generated  <<  92 //PostStepPoint because the hit is generated manually when the photon is
 98 // absorbed by the photocathode                <<  93 //absorbed by the photocathode
 99                                                <<  94 
100 G4bool LXePMTSD::ProcessHits_boundary(const G4 <<  95 G4bool LXePMTSD::ProcessHits_constStep(const G4Step* aStep,
101 {                                              <<  96                                        G4TouchableHistory* ){
102   // need to know if this is an optical photon <<  97 
103   if (aStep->GetTrack()->GetDefinition() != G4 <<  98   //need to know if this is an optical photon
104     return false;                              <<  99   if(aStep->GetTrack()->GetDefinition()
105                                                << 100      != G4OpticalPhoton::OpticalPhotonDefinition()) return false;
106   // User replica number 1 since photocathode  << 101  
107   // to the pmt which was replicated           << 102   //User replica number 1 since photocathode is a daughter volume
108   G4int pmtNumber = aStep->GetPostStepPoint()- << 103   //to the pmt which was replicated
109   G4VPhysicalVolume* physVol = aStep->GetPostS << 104   G4int pmtNumber=
                                                   >> 105     aStep->GetPostStepPoint()->GetTouchable()->GetReplicaNumber(1);
                                                   >> 106   G4VPhysicalVolume* physVol=
                                                   >> 107     aStep->GetPostStepPoint()->GetTouchable()->GetVolume(1);
110                                                   108 
111   // Find the correct hit collection           << 109   //Find the correct hit collection
112   size_t n = fPMTHitCollection->entries();     << 110   G4int n=fPMTHitCollection->entries();
113   LXePMTHit* hit = nullptr;                       111   LXePMTHit* hit = nullptr;
114   for (size_t i = 0; i < n; ++i) {             << 112   for(G4int i=0;i<n;i++){
115     if ((*fPMTHitCollection)[i]->GetPMTNumber( << 113     if((*fPMTHitCollection)[i]->GetPMTNumber()==pmtNumber){
116       hit = (*fPMTHitCollection)[i];           << 114       hit=(*fPMTHitCollection)[i];
117       break;                                      115       break;
118     }                                             116     }
119   }                                               117   }
120                                                << 118  
121   if (hit == nullptr) {  // this pmt wasn't pr << 119   if (hit == nullptr) {//this pmt wasnt previously hit in this event
122     hit = new LXePMTHit();  // so create new h << 120     hit = new LXePMTHit(); //so create new hit
123     hit->SetPMTNumber(pmtNumber);                 121     hit->SetPMTNumber(pmtNumber);
124     hit->SetPMTPhysVol(physVol);                  122     hit->SetPMTPhysVol(physVol);
125     fPMTHitCollection->insert(hit);               123     fPMTHitCollection->insert(hit);
126     hit->SetPMTPos((*fPMTPositionsX)[pmtNumber << 124     hit->SetPMTPos((*fPMTPositionsX)[pmtNumber],(*fPMTPositionsY)[pmtNumber],
127                    (*fPMTPositionsZ)[pmtNumber    125                    (*fPMTPositionsZ)[pmtNumber]);
128   }                                               126   }
129                                                   127 
130   hit->IncPhotonCount();  // increment hit for << 128   hit->IncPhotonCount(); //increment hit for the selected pmt
131                                                << 129  
132   if (!LXeDetectorConstruction::GetSphereOn()) << 130   if(!LXeDetectorConstruction::GetSphereOn()){
133     hit->SetDrawit(true);                         131     hit->SetDrawit(true);
134     // If the sphere is disabled then this hit << 132     //If the sphere is disabled then this hit is automaticaly drawn
135   }                                               133   }
136   else {  // sphere enabled                    << 134   else{//sphere enabled
137     auto trackInfo = (LXeUserTrackInformation* << 135     LXeUserTrackInformation* trackInfo=
138     if (trackInfo->GetTrackStatus() & hitSpher << 136       (LXeUserTrackInformation*)aStep->GetTrack()->GetUserInformation();
139       // only draw this hit if the photon has  << 137     if(trackInfo->GetTrackStatus()&hitSphere)
                                                   >> 138       //only draw this hit if the photon has hit the sphere first
140       hit->SetDrawit(true);                       139       hit->SetDrawit(true);
141   }                                               140   }
142                                                   141 
143   return true;                                    142   return true;
144 }                                                 143 }
                                                   >> 144 
                                                   >> 145 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >> 146 
                                                   >> 147 void LXePMTSD::EndOfEvent(G4HCofThisEvent* ) {}
                                                   >> 148 
                                                   >> 149 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >> 150 
                                                   >> 151 void LXePMTSD::clear() {}
                                                   >> 152 
                                                   >> 153 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >> 154 
                                                   >> 155 void LXePMTSD::DrawAll() {}
                                                   >> 156 
                                                   >> 157 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >> 158 
                                                   >> 159 void LXePMTSD::PrintAll() {}
145                                                   160