geant4系列講座四(參考版)
2024-11-17 04:43本頁面
【正文】 因此散射角劃分可以方便的表示為：,第三十五頁，共三十五頁。G4VprimitiveScorer為每次事件產(chǎn)生G4ThitsMap對象，map的名稱和G4VprimitiveScorer相同。trackerCollectioninsert( newHit )。,內(nèi)容總結(jié),提綱。,因此散射角劃分可以方便的表示為：,其中r是散射光子的方向，r0是入射光子的方向 同樣，把方位角也用余玄表示，但是方位角取值在[0,2Pi]之間，其余玄不單調(diào)，此時就要分情況討論，根據(jù)所處不同象限確定。,對于任意散射角 θs有 因此,均勻劃分散射角的判定條件:,第三十二頁，共三十五頁。,只考慮散射角上的劃分時，對Ф在[0,Pi]上積分，立體角可表示為：,假設(shè)在散射角上把立體劃分為n等分，對應(yīng)每一等分的端點為：,第三十頁，共三十五頁。,第二十八頁，共三十五頁。因此，如果對散射角均勻劃分，那么得不到均勻的立體角（在靠近90度時候最大）。,第二十七頁，共三十五頁。,第二十六頁，共三十五頁。比如以Gamma射線的散射為例，模擬Gamma光子不同立體角上的分布。如： 2008Sep27.16.txt 2008Sep28.16.txt,第二十五頁，共三十五頁。 //****file out~~,這段代碼是將數(shù)據(jù)輸出到一指定路徑的文件下。i++) { alloutGetTheta()GetPhi()GetCounter()G4endl。,std::ofstream allout(name,std::ios::app)。,按自己需要的格式將需要的數(shù)據(jù)結(jié)果輸出：,//****file out~~ //***get local time of the system and set file name time_t date。,復(fù)雜的數(shù)據(jù)輸出，或者輸出有大量數(shù)據(jù)時候可以考慮把數(shù)據(jù)結(jié)果輸出到文件保存。,五、文件輸出技巧,Geant4沒有默認的數(shù)據(jù)輸出，用戶必須自行設(shè)計所有的數(shù)據(jù)輸出。 用預(yù)處理命令查錯。 G4cout“注釋”變量G4endl。 如果在編譯階段就不通過，那肯定是語法什么地方有問題，通過make編譯，不通過時可以看到錯誤出在那個文件中，定義頭文件注意有；號。 },SteppingAction.cc,在SteppingAction.cc中對RunAction中定義的CountProcesses()調(diào)用,第二十一頁，共三十五頁。 G4String procName = endPointGetProcessDefinedStep()GetProcessName()。,void SteppingAction::UserSteppingAction(const G4Step* aStep) { G4StepPoint* prePoint = aStepGetPreStepPoint()。 G4cout “\t“ procName “ = “ count“\n“。i++) { G4String procName = (*ProcCounter)[i]GetName()。 for (size_t i=0。 (*ProcCounter)[i]Count()。 while ((iGetName()!=procName)) i++。,void RunAction::CountProcesses(G4String procName) { // is the process already counted ? size_t nbProc = ProcCountersize()。 ProcCounter = new ProcessesCount。,include “RunAction.hh“ include “G4Run.hh“ include “G4RunManager.hh“ include “Randomize.hh“ include RunAction::RunAction() : ProcCounter(0) { } RunAction::~RunAction() { } void RunAction::BeginOfRunAction(const G4Run* aRun) { G4cout G
合同協(xié)議相關(guān)推薦
鄂ICP備17016276號-1