Geant4  v4-10.4-release
 모두 클래스 네임스페이스들 파일들 함수 변수 타입정의 열거형 타입 열거형 멤버 Friends 매크로 그룹들 페이지들
extended/electromagnetic/TestEm6/src/SteppingAction.cc
이 파일의 문서화 페이지로 가기
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 //
28 //
29 // $Id: SteppingAction.cc 83428 2014-08-21 15:46:01Z gcosmo $
30 //
31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
33 
34 #include "SteppingAction.hh"
35 #include "RunAction.hh"
36 #include "G4SteppingManager.hh"
37 #include "G4VProcess.hh"
38 #include "G4ParticleTypes.hh"
39 
40 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
41 
43 :G4UserSteppingAction(),fRunAction(RuAct)
44 {
46 }
47 
48 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
49 
51 { }
52 
53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
54 
56 {
57  const G4VProcess* process = aStep->GetPostStepPoint()->GetProcessDefinedStep();
58  if (process == 0) return;
59  G4String processName = process->GetProcessName();
60  fRunAction->CountProcesses(processName); //count processes
61 
62  if (processName != "GammaToMuPair") return;
63 
64  G4StepPoint* PrePoint = aStep->GetPreStepPoint();
65  G4double EGamma = PrePoint->GetTotalEnergy();
66  G4ThreeVector PGamma = PrePoint->GetMomentum();
67 
68  G4double Eplus(0), Eminus(0);
69  G4ThreeVector Pplus , Pminus;
70  const G4TrackVector* secondary = fpSteppingManager->GetSecondary();
71  for (size_t lp=0; lp<(*secondary).size(); lp++) {
72  if ((*secondary)[lp]->GetDefinition()==G4MuonPlus::MuonPlusDefinition()) {
73  Eplus = (*secondary)[lp]->GetTotalEnergy();
74  Pplus = (*secondary)[lp]->GetMomentum();
75  } else {
76  Eminus = (*secondary)[lp]->GetTotalEnergy();
77  Pminus = (*secondary)[lp]->GetMomentum();
78  }
79  }
80 
81  G4double xPlus = Eplus/EGamma, xMinus = Eminus/EGamma;
82  G4double thetaPlus = PGamma.angle(Pplus), thetaMinus = PGamma.angle(Pminus);
83  G4double GammaPlus = Eplus/fMuonMass;
84  G4double GammaMinus= Eminus/fMuonMass;
85 
86  G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
87 
88  if(0.0 == thetaPlus || 0.0 == thetaMinus) {
89  G4cout << "SteppingAction: "
90  << "thetaPlus= " << thetaPlus << " thetaMinus= " << thetaMinus
91  << " gamPlus= " << GammaPlus << " gamMinus= " << GammaMinus
92  << " " << thetaPlus *GammaPlus - thetaMinus*GammaMinus << G4endl;
93  return;
94  }
95  analysisManager->FillH1(1,1./(1.+std::pow(thetaPlus*GammaPlus,2)));
96  analysisManager->FillH1(2,std::log10(thetaPlus*GammaPlus));
97 
98  analysisManager->FillH1(3,std::log10(thetaMinus*GammaMinus));
99  analysisManager->FillH1(4,std::log10(std::fabs(thetaPlus *GammaPlus
100  -thetaMinus*GammaMinus)));
101 
102  analysisManager->FillH1(5,xPlus);
103  analysisManager->FillH1(6,xMinus);
104 }
105 
106 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
107 
108 
G4double GetTotalEnergy() const
G4StepPoint * GetPreStepPoint() const
#define G4endl
Definition: G4ios.hh:61
G4double GetPDGMass() const
double angle(const Hep3Vector &) const
double G4double
Definition: G4Types.hh:76
const G4TrackVector * GetSecondary() const
const G4String & GetProcessName() const
Definition: G4VProcess.hh:411
G4CsvAnalysisManager G4AnalysisManager
Definition: g4csv_defs.hh:77
Definition: G4Step.hh:76
G4StepPoint * GetPostStepPoint() const
void UserSteppingAction(const G4Step *)
G4SteppingManager * fpSteppingManager
std::vector< G4Track * > G4TrackVector
static G4MuonPlus * MuonPlus()
Definition: G4MuonPlus.cc:99
G4GLOB_DLL std::ostream G4cout
G4ThreeVector GetMomentum() const
static G4MuonPlus * MuonPlusDefinition()
Definition: G4MuonPlus.cc:94
const G4VProcess * GetProcessDefinedStep() const