Geant4 Cross Reference |
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 /// \file electromagnetic/TestEm16/src/Steppin 27 /// \brief Implementation of the SteppingActio 28 // 29 // 30 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oo 32 33 #include "SteppingAction.hh" 34 35 #include "HistoManager.hh" 36 #include "Run.hh" 37 38 #include "G4EmProcessSubType.hh" 39 #include "G4ParticleTypes.hh" 40 #include "G4RunManager.hh" 41 #include "G4SteppingManager.hh" 42 #include "G4SystemOfUnits.hh" 43 #include "G4UnitsTable.hh" 44 #include "G4VProcess.hh" 45 46 //....oooOO0OOooo........oooOO0OOooo........oo 47 48 SteppingAction::SteppingAction() : G4UserStepp 49 50 //....oooOO0OOooo........oooOO0OOooo........oo 51 52 void SteppingAction::UserSteppingAction(const 53 { 54 Run* run = static_cast<Run*>(G4RunManager::G 55 56 const G4VProcess* process = aStep->GetPostSt 57 if (process == nullptr) return; 58 59 G4int eventID = G4RunManager::GetRunManager( 60 G4Track* trk = aStep->GetTrack(); 61 62 thread_local G4int iCalled = 0; 63 const G4int nprint = 10; // set to 10 to ge 64 G4AnalysisManager* analysisManager = G4Analy 65 66 if (fSynchrotronRadiation == process->GetPro 67 const G4StepPoint* PrePoint = aStep->GetPr 68 const G4TrackVector* secondary = fpSteppin 69 size_t lp = (*secondary).size(); 70 if (lp) { 71 G4double Egamma = (*secondary)[lp - 1]-> 72 run->f_n_gam_sync++; 73 run->f_e_gam_sync += Egamma; 74 run->f_e_gam_sync2 += Egamma * Egamma; 75 if (Egamma > run->f_e_gam_sync_max) run- 76 run->f_lam_gam_sync += aStep->GetStepLen 77 if (iCalled < nprint) { 78 G4double Eelec = PrePoint->GetTotalEne 79 G4ThreeVector Pelec = PrePoint->GetMom 80 G4ThreeVector PGamma = (*secondary)[lp 81 G4bool IsGamma = ((*secondary)[lp - 1] 82 const G4int wid = 8; 83 if (analysisManager->GetVerboseLevel() 84 G4cout << __FILE__ << " line " << __ 85 << " processName=" << process 86 << " Step Length=" << std::se 87 << G4BestUnit(aStep->GetStepL 88 << G4BestUnit(Eelec, "Energy" 89 << " IsGamma=" << IsGamma << 90 << G4BestUnit(Egamma, "Energy 91 << G4BestUnit(PGamma, "Energy 92 } 93 analysisManager->FillH1(1, Egamma); 94 analysisManager->FillH1(2, Egamma, 95 Egamma / keV); 96 analysisManager->FillH1(3, aStep->GetSte 97 } 98 } 99 100 thread_local G4ThreeVector IncomingPhotonDir 101 if (fGammaReflection == process->GetProcessS 102 const G4TrackVector* secondary = fpSteppin 103 size_t lp = (*secondary).size(); 104 if (lp) { 105 G4double Egamma = (*secondary)[lp - 1]-> 106 ++run->f_n_Xray_Refl; 107 analysisManager->FillH1(1, Egamma, 1. / 108 if (analysisManager->GetVerboseLevel() > 109 G4cout << __FILE__ << " line " << __LI 110 << " iCalled=" << std::setw(3) 111 << " f_n_Xray_Refl=" << run->f_ 112 << " ProcessName=" << process-> 113 << " direction=" << trk->GetMom 114 << " keV" << G4endl; 115 IncomingPhotonDirection = trk->GetMoment 116 } 117 } 118 119 const G4VProcess* creator = trk->GetCreatorP 120 G4int CreatorProcessSubType = 0; 121 if (creator != nullptr) CreatorProcessSubTyp 122 if (CreatorProcessSubType == fGammaReflectio 123 if (IncomingPhotonDirection != G4ThreeVect 124 G4double cos_angle = 125 IncomingPhotonDirection * trk->GetMome 126 G4double IncidentAngle = std::acos(cos_a 127 analysisManager->FillH1(4, IncidentAngle 128 if (analysisManager->GetVerboseLevel() > 129 G4cout << __FILE__ << " line " << __LI 130 << " iCalled=" << std::setw(3) 131 << " CreatorProcname=" << creat 132 << " ProcessName=" << process-> 133 << " IncomingPhotonDirection=" 134 << " outgoing PhotonDirection= 135 << " IncidentAngle=" << Inciden 136 << " NumberOfEventToBeProcessed 137 << " RunID=" << run->GetRunID() 138 } 139 } 140 141 ++iCalled; 142 } 143 144 //....oooOO0OOooo........oooOO0OOooo........oo 145