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 // 27 /// \file Par02TrackingAction.cc 28 /// \brief Implementation of the Par02TrackingAction class 29 30 #include "Par02TrackingAction.hh" 31 32 #include "Par02EventInformation.hh" 33 #include "Par02Output.hh" 34 #include "Par02PrimaryParticleInformation.hh" 35 36 #include "G4EventManager.hh" 37 #include "G4RunManager.hh" 38 #include "G4SystemOfUnits.hh" 39 #include "G4ThreeVector.hh" 40 #include "G4TrackingManager.hh" 41 #include "Randomize.hh" 42 43 #include <iomanip> 44 #include <vector> 45 46 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 47 48 Par02TrackingAction::Par02TrackingAction() : G4UserTrackingAction() {} 49 50 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 51 52 Par02TrackingAction::~Par02TrackingAction() = default; 53 54 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 55 56 void Par02TrackingAction::PreUserTrackingAction(const G4Track* aTrack) 57 { 58 // Kill the tracks that have a small transverse momentum or that are not 59 // in the central region. 60 if (aTrack->GetMomentum().perp() < 1.0 * MeV 61 || std::abs(aTrack->GetMomentum().pseudoRapidity()) > 5.5) 62 { 63 ((G4Track*)aTrack)->SetTrackStatus(fStopAndKill); 64 } 65 } 66 67 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 68 69 void Par02TrackingAction::PostUserTrackingAction(const G4Track* aTrack) 70 { 71 if (aTrack->GetTrackStatus() == fStopAndKill && aTrack->GetParentID() == 0) { 72 auto info = (Par02PrimaryParticleInformation*)aTrack->GetDynamicParticle() 73 ->GetPrimaryParticle() 74 ->GetUserInformation(); 75 // info->Print(); 76 Par02Output::Instance()->SaveTrack(Par02Output::eSaveMC, info->GetPartID(), info->GetPDG(), 77 info->GetMCMomentum() / MeV); 78 Par02Output::Instance()->SaveTrack(Par02Output::eSaveTracker, info->GetPartID(), info->GetPDG(), 79 info->GetTrackerMomentum() / MeV, 80 info->GetTrackerResolution(), info->GetTrackerEfficiency()); 81 Par02Output::Instance()->SaveTrack(Par02Output::eSaveEMCal, info->GetPartID(), info->GetPDG(), 82 info->GetEMCalPosition() / mm, info->GetEMCalResolution(), 83 info->GetEMCalEfficiency(), info->GetEMCalEnergy() / MeV); 84 Par02Output::Instance()->SaveTrack(Par02Output::eSaveHCal, info->GetPartID(), info->GetPDG(), 85 info->GetHCalPosition() / mm, info->GetHCalResolution(), 86 info->GetHCalEfficiency(), info->GetHCalEnergy() / MeV); 87 } 88 } 89 90 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 91