Geant4 Cross Reference

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


  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 //                                                 27 //
 28 //  CaTS (Calorimetry and Tracking Simulation)     28 //  CaTS (Calorimetry and Tracking Simulation)
 29 //                                                 29 //
 30 //  Authors : Hans Wenzel                          30 //  Authors : Hans Wenzel
 31 //            Soon Yung Jun                        31 //            Soon Yung Jun
 32 //            (Fermi National Accelerator Labo     32 //            (Fermi National Accelerator Laboratory)
 33 //                                                 33 //
 34 // History                                         34 // History
 35 //   October 18th, 2021 : first implementation     35 //   October 18th, 2021 : first implementation
 36 //                                                 36 //
 37 // *******************************************     37 // ********************************************************************
 38 //                                                 38 //
 39 /// \file InteractionSD.cc                         39 /// \file InteractionSD.cc
 40 /// \brief Implementation of the CaTS::Interac     40 /// \brief Implementation of the CaTS::InteractionSD class
 41                                                    41 
 42 // Geant4 headers:                                 42 // Geant4 headers:
 43 #include "G4HCofThisEvent.hh"                      43 #include "G4HCofThisEvent.hh"
 44 #include "G4Step.hh"                               44 #include "G4Step.hh"
 45 #include "G4ThreeVector.hh"                        45 #include "G4ThreeVector.hh"
 46 #include "G4SDManager.hh"                          46 #include "G4SDManager.hh"
 47 #include "G4ios.hh"                                47 #include "G4ios.hh"
 48 #include "G4SystemOfUnits.hh"                      48 #include "G4SystemOfUnits.hh"
 49 // project headers                                 49 // project headers
 50 #include "InteractionSD.hh"                        50 #include "InteractionSD.hh"
 51 #include "ParticleChange.hh"                       51 #include "ParticleChange.hh"
 52                                                    52 
 53 InteractionSD::InteractionSD(G4String name)        53 InteractionSD::InteractionSD(G4String name)
 54   : G4VSensitiveDetector(name)                     54   : G4VSensitiveDetector(name)
 55 {                                                  55 {
 56   G4String HCname = name + "_HC";                  56   G4String HCname = name + "_HC";
 57   collectionName.insert(HCname);                   57   collectionName.insert(HCname);
 58   G4cout << collectionName.size() << "   Inter     58   G4cout << collectionName.size() << "   InteractionSD name:  " << name
 59          << " collection Name: " << HCname <<      59          << " collection Name: " << HCname << G4endl;
 60   fHCID = -1;                                      60   fHCID = -1;
 61 }                                                  61 }
 62                                                    62 
 63 InteractionSD::~InteractionSD()                    63 InteractionSD::~InteractionSD()
 64 {                                                  64 {
 65   delete fFirstInter;                              65   delete fFirstInter;
 66   delete fOtherInter;                              66   delete fOtherInter;
 67 }                                                  67 }
 68                                                    68 
 69 void InteractionSD::Initialize(G4HCofThisEvent     69 void InteractionSD::Initialize(G4HCofThisEvent* HCE)
 70 {                                                  70 {
 71   G4cout << "Hits Collection capacity:  " << H     71   G4cout << "Hits Collection capacity:  " << HCE->GetCapacity() << G4endl;
 72   fInteractionHitsCollection =                     72   fInteractionHitsCollection =
 73     new InteractionHitsCollection(SensitiveDet     73     new InteractionHitsCollection(SensitiveDetectorName, collectionName[0]);
 74   if(fHCID < 0)                                    74   if(fHCID < 0)
 75   {                                                75   {
 76     G4cout << "InteractionSD::Initialize:  " <     76     G4cout << "InteractionSD::Initialize:  " << SensitiveDetectorName << "   "
 77            << collectionName[0] << G4endl;         77            << collectionName[0] << G4endl;
 78     fHCID = G4SDManager::GetSDMpointer()->GetC     78     fHCID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
 79   }                                                79   }
 80   HCE->AddHitsCollection(fHCID, fInteractionHi     80   HCE->AddHitsCollection(fHCID, fInteractionHitsCollection);
 81   fFirstInter = new ParticleChange(true);          81   fFirstInter = new ParticleChange(true);
 82   fOtherInter = new ParticleChange();              82   fOtherInter = new ParticleChange();
 83 }                                                  83 }
 84                                                    84 
 85 G4bool InteractionSD::ProcessHits(G4Step* aSte     85 G4bool InteractionSD::ProcessHits(G4Step* aStep, G4TouchableHistory*)
 86 {                                                  86 {
 87   G4int nsc = fFirstInter->GetNumberOfSecondar     87   G4int nsc = fFirstInter->GetNumberOfSecondaries();
 88   if(nsc > 0)                                      88   if(nsc > 0)
 89   {                                                89   {
 90     for(G4int i = 0; i < nsc; i++)                 90     for(G4int i = 0; i < nsc; i++)
 91     {                                              91     {
 92       delete fFirstInter->GetSecondary(i);         92       delete fFirstInter->GetSecondary(i);
 93     }                                              93     }
 94     fFirstInter->Clear();                          94     fFirstInter->Clear();
 95   }                                                95   }
 96   nsc = fOtherInter->GetNumberOfSecondaries();     96   nsc = fOtherInter->GetNumberOfSecondaries();
 97   if(nsc > 0)                                      97   if(nsc > 0)
 98   {                                                98   {
 99     for(G4int i = 0; i < nsc; i++)                 99     for(G4int i = 0; i < nsc; i++)
100     {                                             100     {
101       delete fOtherInter->GetSecondary(i);        101       delete fOtherInter->GetSecondary(i);
102     }                                             102     }
103     fOtherInter->Clear();                         103     fOtherInter->Clear();
104   }                                               104   }
105   const std::vector<const G4Track*>* secs = aS    105   const std::vector<const G4Track*>* secs = aStep->GetSecondaryInCurrentStep();
106   G4int nsec                              = se    106   G4int nsec                              = secs->size();
107   for(G4int i = 0; i < nsec; i++)                 107   for(G4int i = 0; i < nsec; i++)
108   {                                               108   {
109     G4Track* tr = new G4Track(*((*secs)[i]));     109     G4Track* tr = new G4Track(*((*secs)[i]));
110     if(aStep->GetTrack()->GetTrackStatus() !=     110     if(aStep->GetTrack()->GetTrackStatus() != fAlive)  // track looses identity
111     {                                             111     {
112       if(aStep->GetTrack()->GetParentID() == 0    112       if(aStep->GetTrack()->GetParentID() == 0)  // primary track
113       {                                           113       {
114         fFirstInter->AddSecondary(tr);            114         fFirstInter->AddSecondary(tr);
115       }                                           115       }
116       else  // secondary track, and it's also     116       else  // secondary track, and it's also looses identity (re-interaction)
117       {                                           117       {
118         fOtherInter->AddSecondary(tr);            118         fOtherInter->AddSecondary(tr);
119       }                                           119       }
120     }                                             120     }
121   }  // end loop over secondaries                 121   }  // end loop over secondaries
122   G4int NSec = fFirstInter->GetNumberOfSeconda    122   G4int NSec = fFirstInter->GetNumberOfSecondaries();
123   if(NSec > 0)                                    123   if(NSec > 0)
124   {                                               124   {
125     const G4DynamicParticle* sec = 0;             125     const G4DynamicParticle* sec = 0;
126     for(G4int i = 0; i < NSec; i++)               126     for(G4int i = 0; i < NSec; i++)
127     {                                             127     {
128       sec = fFirstInter->GetSecondary(i)->GetD    128       sec = fFirstInter->GetSecondary(i)->GetDynamicParticle();
129       const G4String& pname  = sec->GetDefinit    129       const G4String& pname  = sec->GetDefinition()->GetParticleName();
130       G4double pmom          = (sec->GetTotalM    130       G4double pmom          = (sec->GetTotalMomentum()) / GeV;
131       G4double Ekin          = (sec->GetKineti    131       G4double Ekin          = (sec->GetKineticEnergy()) / GeV;
132       G4double theta         = (sec->GetMoment    132       G4double theta         = (sec->GetMomentum()).theta();
133       InteractionHit* newHit = new Interaction    133       InteractionHit* newHit = new InteractionHit(pname, pmom, Ekin, theta);
134       fInteractionHitsCollection->insert(newHi    134       fInteractionHitsCollection->insert(newHit);
135     }                                             135     }
136   }                                               136   }
137   return true;                                    137   return true;
138 }                                                 138 }
139                                                   139