Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/CaTS/src/CalorimeterSD.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/advanced/CaTS/src/CalorimeterSD.cc (Version 11.3.0) and /examples/advanced/CaTS/src/CalorimeterSD.cc (Version 10.3)


  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 // *******************************************    
 27 //                                                
 28 //  CaTS (Calorimetry and Tracking Simulation)    
 29 //                                                
 30 //  Authors : Hans Wenzel                         
 31 //            Soon Yung Jun                       
 32 //            (Fermi National Accelerator Labo    
 33 //                                                
 34 // History                                        
 35 //   October 18th, 2021 : first implementation    
 36 //                                                
 37 // *******************************************    
 38 //                                                
 39 /// \file CalorimeterSD.cc                        
 40 /// \brief Implementation of the CaTS::Calorim    
 41                                                   
 42 // Geant4 headers                                 
 43 #include "G4HCofThisEvent.hh"                     
 44 #include "G4Step.hh"                              
 45 #include "G4SDManager.hh"                         
 46 #include "G4ios.hh"                               
 47 // project headers                                
 48 #include "CalorimeterSD.hh"                       
 49 #include "ConfigurationManager.hh"                
 50                                                   
 51 CalorimeterSD::CalorimeterSD(G4String name)       
 52   : G4VSensitiveDetector(name)                    
 53 {                                                 
 54   G4String HCname = name + "_HC";                 
 55   collectionName.insert(HCname);                  
 56   G4cout << collectionName.size() << "   Calor    
 57          << " collection Name: " << HCname <<     
 58   fHCID   = -1;                                   
 59   verbose = ConfigurationManager::getInstance(    
 60 }                                                 
 61                                                   
 62 void CalorimeterSD::Initialize(G4HCofThisEvent    
 63 {                                                 
 64   fCalorimeterHitsCollection =                    
 65     new CalorimeterHitsCollection(SensitiveDet    
 66   if(fHCID < 0)                                   
 67   {                                               
 68     if(verbose)                                   
 69       G4cout << "CalorimeterSD::Initialize:  "    
 70              << collectionName[0] << G4endl;      
 71     fHCID = G4SDManager::GetSDMpointer()->GetC    
 72   }                                               
 73   hce->AddHitsCollection(fHCID, fCalorimeterHi    
 74 }                                                 
 75                                                   
 76 G4bool CalorimeterSD::ProcessHits(G4Step* aSte    
 77 {                                                 
 78   // auto can be used for all results from ret    
 79   auto edep = aStep->GetTotalEnergyDeposit() /    
 80   if(edep == 0.)                                  
 81     return false;                                 
 82   auto time       = aStep->GetPreStepPoint()->    
 83   auto touch      = aStep->GetPreStepPoint()->    
 84   auto cellpos    = touch->GetTranslation();      
 85   unsigned int ID = aStep->GetPreStepPoint()->    
 86   auto theTrack   = aStep->GetTrack();            
 87   auto particleType = theTrack->GetDefinition(    
 88   //                                              
 89   //  check if this cell has been hit before      
 90   // fCalorimeterHitsCollection                   
 91   for(unsigned int j = 0; j < fCalorimeterHits    
 92   {                                               
 93     CalorimeterHit* aPreviousHit = (*fCalorime    
 94     if(ID == aPreviousHit->GetId())               
 95     {                                             
 96       aPreviousHit->SetEdep(aStep->GetTotalEne    
 97                             aPreviousHit->GetE    
 98       if((particleType == "e+") || (particleTy    
 99          (particleType == "e-"))                  
100       {                                           
101         aPreviousHit->Setem_Edep(edep + aPrevi    
102       }                                           
103       return true;                                
104     }                                             
105   }                                               
106   //                                              
107   // otherwise create a new hit:                  
108   //                                              
109   CalorimeterHit* newHit;                         
110   if((particleType == "e+") || (particleType =    
111      (particleType == "e-"))                      
112   {                                               
113     newHit = new CalorimeterHit(ID, edep, edep    
114   }                                               
115   else                                            
116   {                                               
117     newHit = new CalorimeterHit(ID, edep, 0.0,    
118   }                                               
119   fCalorimeterHitsCollection->insert(newHit);     
120   return true;                                    
121 }                                                 
122                                                   
123 void CalorimeterSD::EndOfEvent(G4HCofThisEvent    
124 {                                                 
125   G4int NbHits = fCalorimeterHitsCollection->e    
126   if(verbose)                                     
127     G4cout << " Number of CalorimeterHits:  "     
128 }                                                 
129