Luận văn thạc sĩ khảo sát đáp ứng của detector hpge cho phóng xạ môi trường bằng phần mềm geant4

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TP HỒ CHÍ MINH Nguyễn Thị Bích Hậu KHẢO SÁT ĐÁP ỨNG CỦA DETECTOR HPGe CHO PHĨNG XẠ MƠI TRƯỜNG BẰNG PHẦN MỀM GEANT4 LUẬN VĂN THẠC SĨ VẬT LÍ Thành phố Hồ Chí Minh - 2013 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TP HỒ CHÍ MINH Nguyễn Thị Bích Hậu KHẢO SÁT ĐÁP ỨNG CỦA DETECTOR HPGe CHO PHĨNG XẠ MƠI TRƯỜNG BẰNG PHẦN MỀM GEANT4 Chuyên ngành: Vật lí nguyên tử Mã số: 60.44.01.06 LUẬN VĂN THẠC SĨ: VẬT LÍ NGƯỜI HƯỚNG DẪN KHOA HỌC: TS Võ Hồng Hải Thành phố Hồ Chí Minh – 2013 LỜI CẢM ƠN Trong suốt trình thực luận văn, tơi nhận nhiều giúp đỡ tận tình, chu đáo với tinh thần khoa học trách nhiệm cao Thầy/Cô khoa Vật lý trường Đại Học Sư Phạm Thành Phố Hồ Chí Minh mơn Vật Lý Hạt Nhân trường Đại Học Khoa Học Tự Nhiên Thành Phố Hồ Chí Minh Tơi xin gửi lời biết ơn chân thành kính trọng đến: TS Võ Hồng Hải khơng Thầy hướng dẫn khoa học tận tình hướng dẫn, động viên truyền đạt kinh nghiệm quý báu nghiên cứu khoa học Bên cạnh đó, Thầy tạo điều kiện tốt để sớm hoàn thành luận văn PGS TS Châu Văn Tạo Thầy/Cô môn Vật Lý Hạt Nhân tận tình giúp đỡ tạo điều kiện thuận lợi để thực luận văn Ths Nguyễn Quốc Hùng hướng dẫn nhiều bước làm quen với chương trình mơ Geant4 Tơi cũng bày tỏ lòng biết ơn sâu sắc đến quý Thầy/Cô truyền đạt kiến thức trình học, Phịng Sau Đại Học, Khoa Vật Lý trường Đại Học Sư Phạm TP Hồ Chí Minh tạo điều kiện thuận lợi giúp đỡ nhiệt tình suốt thời gian thực luận văn Cảm ơn bạn bè động viên Cuối cùng, xin gửi lời biết ơn chân thành tới ba mẹ, gia đình tạo điều kiện thuận lợi để thực việc học làm luận văn tốt MỤC LỤC LỜI CẢM ƠN MỤC LỤC DANH MỤC CÁC KÝ HIỆU, CÁC CHỮ VIẾT TẮT MỞ ĐẦU CHƯƠNG 1: TỔNG QUAN 1.1 Tương tác xạ gamma với vật chất 1.1.1 Hiệu ứng quang điện 1.1.2 Tán xạ Compton 1.1.3 Hiệu ứng tạo cặp 10 1.1.4 Hiệu ứng Rayleigh 12 1.2 Tương tác xạ muon với vật chất 14 1.3 Giới thiệu phóng xạ mơi trường 16 1.3.1 Các chuỗi phóng xạ nguyên thủy tự nhiên 17 1.3.2 Phóng xạ xạ vũ trụ muon 21 1.4 Giới thiệu hệ phổ kế HPGe GC2018 22 1.4.1 Giới thiệu 22 1.4.2 Mô tả đầu dò HPGe GC2018 23 1.5 Giới thiệu chương trình mơ Geant4 25 1.5.1 Giới thiệu chung 25 1.5.2 Cấu trúc chương trình Geant4 25 CHƯƠNG 2: BỐ TRÍ MƠ PHỎNG HỆ ĐO HPGE 27 2.1 Bố trí mơ 27 2.2 Chương trình Geant4 28 2.2.1 Kết cấu hình học tính chất đối tượng 28 2.2.2 Tương tác vật lý 37 2.2.3 Các thiết lập sở ban đầu 38 2.3 Chương trình tác động độ phân giải 40 2.4 Chương trình lưu liệu vẽ phổ 41 2.4.1 Chương trình lưu liệu 41 2.4.2 Chương trình vẽ phổ 42 CHƯƠNG 3: KẾT QUẢ VÀ NHẬN XÉT 43 3.1 Dạng phổ đồng vị K-40 nhận xét 43 3.2 Dạng phổ đồng vị U-238 nhận xét 46 3.3 Dạng phổ đồng vị Th-232 nhận xét 49 3.4 Dạng phổ muon nhận xét 52 KẾT LUẬN VÀ KIẾN NGHỊ 53 TÀI LIỆU THAM KHẢO 55 PHỤ LỤC 57 DANH MỤC CÁC KÝ HIỆU, CÁC CHỮ VIẾT TẮT Ký hiệu Diễn giải h Hằng số Planck 6,626.10-34J.s c Vận tốc ánh sáng chân không 3.108m/s Na Số Avogadro 6,02.1023 hạt/mol m0 Khối lượng nghỉ electron 9,1.10-31kg re Bán kính quỹ đạo cổ điển electron 2,8179.10-15m σ0 Độ lệch chuẩn Ee Động quang electron E Năng lượng gamma tới EB Năng lượng liên kết electron lớp vỏ Ek Năng lượng liên kết electron lớp K Z Số hiệu nguyên tử Chữ viết tắt Diễn giải Geant4 Geometry and Tracking CERN Conseil Européen pour la Recherche Nucléaire (European Organization for Nuclear Research) FWHM Full Width Hafl Maximum HPGe Hight Pure Germanium Tp.HCM Thành phố Hồ Chí Minh LNHB Laboratoire National Henri Becquerel MCNP Monte Carlo N Particle MỞ ĐẦU Trong môi trường có nhân phóng xạ nguyên thủy phổ biến U-238, Th-232, U-235, K-40 Các nhân phóng xạ có chu kỳ bán rã lớn so sánh với tuổi trái đất Trong nồng độ U-238 có giá trị trung bình tồn giới 25 Bq/kg Nồng độ Th-232 đất dải từ đến 300 Bq/kg tương tự U-238 có giá trị trung bình tồn giới 40 Bq/kg Nồng độ K-40 đất dải từ 37 đến 1100 Bq/kg có giá trị trung bình 400 Bq/kg Cịn lượng nhân phóng xạ U-235 chiếm 0,72% tổng lượng uranium có tự nhiên nên có mơi trường đất Việc nghiên cứu phóng xạ môi trường khả ghi nhận đầu dị, phân tích phổ cho phép xác định hoạt độ đồng vị phóng xạ Để đánh giá hoạt độ mẫu môi trường thường dựa vào nhiều tham số hiệu suất ghi đầu dò, độ phân giải lượng đặt biệt phơng phóng xạ Đối với phơng nền, nguyên nhân phông Compton, phông xạ vũ trụ Sự tự hấp thụ mẫu nguyên nhân Để nắm rõ ảnh hưởng phông khả ghi nhận, tính chất hoạt động đầu dị, bên cạnh thực nghiệm, việc mơ tính chất đầu dị quan trọng Hiện nay, phần mềm có độ tin cậy cao dùng việc mô gồm GEANT4, MCNP, phần mềm chuyên dụng khác Đã có số cơng trình ngồi nước thực mơ số tính chất cho đầu dị HPGe Chẳng hạn số cơng trình tiêu biểu: “Nghiên cứu hiệu suất ghi đầu dò bán dẫn siêu tinh khiết (HPGe) hệ phổ kế gamma phương pháp Monte Carlo thuật toán di truyền” luận án tiến sĩ Võ Xuân Ân năm 2008; “Mô Monte Carlo đường cong hiệu suất đỉnh đầu dị HPGe hệ phổ kế gamma mơi trường sử dụng chương trình MCNP4C2” báo đăng tạp chí khoa học phát triển cơng nghệ nhóm tác giả Trương Thị Hồng Loan, Đặng Nguyên Phương, Trần Ái Khanh Mai Văn Nhơn; “Dead-layer thickness effect for gamma spectra measured in an HPGe p-type detector” tác giả Huy, N.Q (2011); “A detailed investigation of HPGe detector response for improved Monte Carlo efficiency calculations” nhóm tác giả Boson J., Agren G., Johansson L (2008); … Trong luận văn này, thực mô khả đáp ứng đầu dò HPGe cho cho số đồng vị phóng xạ có mơi trường xạ vũ trụ muon Cụ thể khảo sát cho đồng vị U-238, Th-232, K-40, với dạng nguồn điểm nguồn có kích thước Sự tự hấp thụ quan tâm đề tài Đầu dị chúng tơi quan tâm loại HPGe đặt buồng chì Các thơng số kích thước, vật liệu cấu tạo đầu dò buồng chì lấy từ nhà sản xuất Chúng tơi mơ tả cụ thể cho đầu dò HPGe mà sử dụng Bộ môn Vật lý Hạt nhân, trường đại học Khoa học Tự nhiên – TpHCM Phần mềm sử dụng mô Geant4 Trong luận văn tiến hành công việc cụ thể sau: 1) Mơ hệ đầu dị HPGe theo kích thước thực 2) Mơ nguồn theo dạng nguồn điểm, nguồn có kích thước phân bố Nguồn phát xạ đẳng hướng pi Các giá trị lượng K-40, U-238, Th-232 tham khảo từ sở liệu phịng thí nghiệm quốc gia Pháp (LNHB) Riêng xạ vũ trụ muon có lượng 450 MeV, dạng tia, hướng bắn vng góc với bề mặt dầu dị 3) Phổ lượng để lại đầu dò ghi lại, phổ lượng có ảnh hưởng độ phân giải theo phương trình FWHM = a + b * sqrt(E), với hệ số a, b lấy từ số liệu thực nghiệm Phòng thí nghiệm Bộ mơn Vật lý hạt nhân, Trường Đại Học Khoa Học Tự nhiên Thành phố Hồ Chí Minh: a = 0,6982239678*10-3; b = 0,902784899*10-3 4) Phân tích, đánh giá kết thu Phần bố trí luận văn phân bổ sau: Chương - Tổng quan: Trong chương giới thiệu khái quát tương tác xạ gamma với vật chất, tương tác xạ muon với vật chất, phóng xạ mơi trường, hệ phổ kế HPGe GC2018 chương trình mơ Geant4 Chương - Bố trí mơ hệ đo HPGe: Trong phần trình bày bố trí hệ đo mơ phỏng, chương trình Geant4, chương trình tác động độ phân giải chương trình lưu liệu vẽ phổ Chương - Kết nhận xét: Trong chương chúng tơi trình bày kết dạng phổ nhận xét kết thu đồng vị: K-40, U-238, Th-232 muon CHƯƠNG 1: TỔNG QUAN 1.1 Tương tác xạ gamma với vật chất 1.1.1 Hiệu ứng quang điện Khi gamma va chạm với electron nguyên tử, lượng gamma truyền tồn cho electron dạng động Electron bay khỏi nguyên tử gọi quang electron Đây gọi hiệu ứng quang điện Quang electron nhận động Ee hiệu số lượng gamma tới E lượng liên kết EB electron lớp vỏ trước bị bứt [1],[4],[8] Germanium e- Năng lượng photon (MeV) (a) (b) Hình 1.1 a) Hiệu ứng quang điện b) Tiết diện khối hiệu ứng quang điện phụ thuộc lượng photon Ee = E - EB (1.1) Theo công thức (1.1) lượng gamma tới phải lớn lượng liên kết electron hiệu ứng quang điện xảy Tương tác xảy với xác suất lớn lượng gamma vừa vượt qua lượng liên kết, đặc biệt lớp (Hình 1.1b) Do lượng liên kết thay đổi theo số nguyên tử Z nên tiết diện tương tác quang điện phụ thuộc vào Z Xác suất tổng cộng hiệu ứng quang điện tất electron quỹ đạo có lượng E ≥ Ek , Ek lượng liên kết electron lớp K, tuân theo quy luật E theo phương trình[16]: σ photo = 4α Với σ =  me c   Z σ   hυ  7/2 8πre2 = 6,651.10− 25 cm ; (1.2) α= tiết diện tương tác Thomson số 137 cấu trúc tinh tế Các công thức cho thấy hiệu ứng quang điện xảy với tiết diện lớn nguyên tử nặng (chẳng hạn chì) vùng lượng cao, nguyên tử nhẹ hiệu ứng quang điện xuất vùng lượng thấp Hiệu ứng quang điện làm ion hóa nguyên tử đưa nguyên tử đến trạng thái kích thích Sau electron tầng cao nhảy xuống lấp lỗ trống khiếm khuyết tầng điện li, trình phát tia X đặc trưng electron Auger Do lượng thấp nên chúng bị hấp thụ hoàn toàn vùng nhạy đầu dò Phần lớn lượng tia gamma chuyển thành động electron ghi nhận Xung ghi hiệu ứng quang điện hiệu ứng bề mặt tia X đóng góp vào số đếm đỉnh lượng tồn phần 1.1.2 Tán xạ Compton e- (a) (b) Hình 1.2 a) Hiệu ứng Compton; b) Sơ đồ tán xạ gamma lên electron tự G4VPhysicalVolume* AltwoPhys; G4VPhysicalVolume* Cu1Phys; G4VPhysicalVolume* Cu2Phys; G4VPhysicalVolume* Cu3Phys; G4VPhysicalVolume* Cu4Phys; G4VPhysicalVolume* Pb1Phys; G4VPhysicalVolume* Pb2Phys; G4VPhysicalVolume* Pb3Phys; G4VPhysicalVolume* Sn1Phys; G4VPhysicalVolume* Sn2Phys; G4VPhysicalVolume* Sn3Phys; G4VPhysicalVolume* Sn4Phys; G4VPhysicalVolume* VacuumPhys; G4VPhysicalVolume* Vacuum3Phys; G4VPhysicalVolume* Vacuum4Phys; G4VPhysicalVolume* airPhys; G4bool constructed; private: void DefineMaterials(); void SetupGeometry(); private: G4Material* air; G4Material* dat; G4Material* Ge; G4Material* Pb; G4Material* Cu; G4Material* Sn; G4Material* Al; G4Material* vacuum; G4VPhysicalVolume* worldPhys; G4VPhysicalVolume* datPhys; G4VPhysicalVolume* GePhys; G4VPhysicalVolume* AlPhys; }; #endif 64 PHỤ LỤC - CODE TƯƠNG TÁC VẬT LÝ PHYSICSLIST.CC #include "PhysicsList.hh" #include "G4ParticleDefinition.hh" #include "G4ParticleTypes.hh" #include "G4ProcessManager.hh" #include "G4LeptonConstructor.hh" #include "G4BosonConstructor.hh" //#include "G4MesonConstructor.hh" #include "G4IonConstructor.hh" #include "G4BaryonConstructor.hh" #include "G4StepLimiter.hh" #include "G4Region.hh" #include "G4EmProcessOptions.hh" PhysicsList::PhysicsList(): G4VUserPhysicsList() { SetVerboseLevel(1); } PhysicsList::~PhysicsList() { void PhysicsList::ConstructParticle() { G4BosonConstructor pBosonConstructor; pBosonConstructor.ConstructParticle(); G4LeptonConstructor pLeptonConstructor; pLeptonConstructor.ConstructParticle(); G4BaryonConstructor pBaryonConstructor; pBaryonConstructor.ConstructParticle(); G4IonConstructor pIonConstructor; pIonConstructor.ConstructParticle(); ConstructBosons(); ConstructLeptons(); ConstructBarions(); ConstructIons(); } void PhysicsList::ConstructBosons() { void PhysicsList::ConstructIons() { // Ions G4IonConstructor ions; ions.ConstructParticle(); } void PhysicsList::ConstructProcess() { AddTransportation(); ConstructEM(); } // gamma from standard #include "G4PhotoElectricEffect.hh" #include "G4ComptonScattering.hh" #include "G4GammaConversion.hh" #include "G4RayleighScattering.hh" // gamma from Lowenergy #include "G4LivermorePhotoElectricModel.hh" #include "G4LivermoreComptonModel.hh" #include "G4LivermorePolarizedComptonModel.hh" // alternative for polarized photons #include "G4LivermoreGammaConversionModel.hh" #include "G4LivermoreRayleighModel.hh" #include "G4LivermorePolarizedRayleighModel.hh" // alternative for polarized photons // e+ - e- from standard #include "G4eMultipleScattering.hh" #include "G4eIonisation.hh" #include "G4eBremsstrahlung.hh" #include "G4UniversalFluctuation.hh" #include "G4eplusAnnihilation.hh" // e+ - e- from Lowenergy #include "G4LivermoreIonisationModel.hh" #include "G4LivermoreBremsstrahlungModel.hh" //muon #include "G4MuMultipleScattering.hh" #include "G4MuIonisation.hh" #include "G4MuBremsstrahlung.hh" #include "G4MuPairProduction.hh" // gamma G4Gamma::GammaDefinition(); } // oooOO0OOooo oooOO0OOooo oooOO0OOooo oooO O0OOooo void PhysicsList::ConstructLeptons() { // leptons G4Electron::ElectronDefinition(); G4Positron::PositronDefinition(); //G4mu::muDefinition(); } void PhysicsList::ConstructBarions() { G4BaryonConstructor baryon; baryon.ConstructParticle(); G4Proton::ProtonDefinition(); } // protons from standard #include "G4hLowEnergyIonisation.hh" #include "G4hMultipleScattering.hh" // options #include "G4LossTableManager.hh" #include "G4EmProcessOptions.hh" #include "G4hIonisation.hh" // oooOO0OOooo oooOO0OOooo oooOO0OOooo oooO O0OOooo void PhysicsList::ConstructEM() { theParticleIterator->reset(); while( (*theParticleIterator)() ) G4ParticleDefinition* particle = theParticleIterator->value(); 65 G4ProcessManager* pmanager = particle->GetProcessManager(); G4String particleName = particle->GetParticleName(); if (particleName == "gamma") { // gamma G4ComptonScattering* theComptonScattering = new G4ComptonScattering(); G4LivermoreComptonModel* theLivermoreComptonModel = new G4LivermoreComptonModel(); theComptonScattering->SetModel(theLivermoreComptonModel); pmanager->AddDiscreteProcess(theComptonScattering); G4PhotoElectricEffect* thePhotoElectricEffect = new G4PhotoElectricEffect(); G4LivermorePhotoElectricModel* theLivermorePhotoElectricModel = new G4LivermorePhotoElectricModel(); theLivermorePhotoElectricModel->ActivateAuger(true); thePhotoElectricEffect->SetModel(theLivermorePhotoElectricModel); pmanager->AddDiscreteProcess(thePhotoElectricEffect); G4GammaConversion* theGammaConversion = new G4GammaConversion(); G4LivermoreGammaConversionModel* theLivermoreGammaConversionModel = new G4LivermoreGammaConversionModel(); theGammaConversion>SetModel(theLivermoreGammaConversionModel); pmanager->AddDiscreteProcess(theGammaConversion); G4RayleighScattering* theRayleigh = new G4RayleighScattering(); G4LivermoreRayleighModel* theRayleighModel = new G4LivermoreRayleighModel(); theRayleigh->SetModel(theRayleighModel); pmanager->AddDiscreteProcess(theRayleigh); pmanager->AddProcess(new G4StepLimiter(), -1, -1, 5); } else if (particleName == "e-") { } else if( particleName == "mu+" || particleName == "mu-" ) { //muon pmanager->AddProcess(new G4MuMultipleScattering,-1, 1, 1); pmanager->AddProcess(new G4MuIonisation, -1, 2, 2); pmanager->AddProcess(new G4MuBremsstrahlung, -1, 3, 3); pmanager->AddProcess(new G4MuPairProduction, -1, 4, 4); } else if (particleName == "proton") { //proton G4hLowEnergyIonisation* hIoni = new G4hLowEnergyIonisation(); hIoni->SetFluorescence(true); hIoni->SelectShellIonisationCS("analytical"); pmanager->AddProcess(new G4hMultipleScattering,-1,1,1); pmanager->AddProcess(hIoni,-1, 2,2); } else if ( particleName == "alpha" ) { pmanager->AddProcess(new G4hMultipleScattering,-1,1,1); G4hLowEnergyIonisation* hIoni = new G4hLowEnergyIonisation(); hIoni->SetFluorescence(true); pmanager->AddProcess(hIoni,-1,2,2); } } //############################### GEANT4 4.9.4 Atomic Deexcitation Process ###############################// G4EmProcessOptions emOptions; emOptions.SetDeexcitationActive(true); } void PhysicsList::SetCuts() { defaultCutValue = 1.*mm; const G4double cutForGamma = 10e-6*mm; const G4double cutForElectron = defaultCutValue; const G4double cutForPositron = defaultCutValue; SetCutValue(cutForGamma, "gamma"); SetCutValue(cutForElectron, "e-"); SetCutValue(cutForPositron, "e+"); G4double lowLimit = 250 * eV; G4double highLimit = 100 * GeV; G4ProductionCutsTable::GetProductionCutsTable() >SetEnergyRange(lowLimit, highLimit); if (verboseLevel>0) DumpCutValuesTable(); } //electron G4eMultipleScattering* msc = new G4eMultipleScattering(); msc->SetStepLimitType(fUseDistanceToBoundary); pmanager->AddProcess(msc, -1, 1, 1); G4eIonisation* eIoni = new G4eIonisation(); G4LivermoreIonisationModel* theLivermoreIonisationModel = new G4LivermoreIonisationModel(); theLivermoreIonisationModel->ActivateAuger(true); eIoni->SetEmModel(theLivermoreIonisationModel); // eIoni->SetStepFunction(0.2, 100*um); // pmanager->AddProcess(eIoni, -1, 2, 2); G4eBremsstrahlung* eBrem = new G4eBremsstrahlung(); eBrem->SetEmModel(new G4LivermoreBremsstrahlungModel()); pmanager->AddProcess(eBrem, -1,-3, 3); pmanager->AddProcess(new G4StepLimiter(), -1, -1, 4); PHYSICSLIST.HH #ifndef PhysicsList_h #define PhysicsList_h #include "G4VUserPhysicsList.hh" #include "globals.hh" class PhysicsList: public G4VUserPhysicsList { public: PhysicsList(); // Ionisation } else if (particleName == "e+") { //positron G4eMultipleScattering* msc = new G4eMultipleScattering(); msc->SetStepLimitType(fUseDistanceToBoundary); pmanager->AddProcess(msc, -1, 1, 1); pmanager->AddProcess(new G4eIonisation(), -1, 2, 2); pmanager->AddProcess(new G4eBremsstrahlung, -1,-1,3); pmanager->AddProcess(new G4eplusAnnihilation, 0,-1,4); 66 ~PhysicsList(); protected: // Construct particle and physics void ConstructParticle(); void ConstructProcess(); // these methods Construct particles virtual void ConstructBosons(); virtual void ConstructLeptons(); virtual void ConstructBarions(); virtual void ConstructIons(); void SetCuts(); private: // These methods construct physics processes and register them void ConstructEM(); }; #endif 67 PHỤ LỤC - CODE NGUỒN – HƯỚNG BẮN – NĂNG LƯỢNG particleGun->GeneratePrimaryVertex(anEvent); GPSPRIMARYGENERATORACTION.CC } #include "GPSPrimaryGeneratorAction.hh" PRIMARYGENERATORACTION.HH #include "G4Event.hh" #ifndef GPSPrimaryGeneratorAction_h #include "G4GeneralParticleSource.hh" #define GPSPrimaryGeneratorAction_h #include "G4ParticleGun.hh" #include "G4VUserPrimaryGeneratorAction.hh" #include "Randomize.hh" #include "globals.hh" #include "G4ParticleTable.hh" class G4GeneralParticleSource; #include "G4ParticleDefinition.hh" class G4Event; //vi du nguon diem classGPSPrimaryGeneratorAction:public GPSPrimaryGeneratorAction::GPSPrimaryGeneratorAction() G4VUserPrimaryGeneratorAction { { particleGun = new G4GeneralParticleSource(); public: } GPSPrimaryGeneratorAction(); GPSPrimaryGeneratorAction::~GPSPrimaryGeneratorAction() virtual ~GPSPrimaryGeneratorAction(); { public: virtual void GeneratePrimaries(G4Event* anEvent); delete particleGun; private: } void G4GeneralParticleSource* particleGun; GPSPrimaryGeneratorAction::GeneratePrimaries(G4Event* anEvent) }; { #endif // Random generation of the impulse direction of primary particles Năng lượng Nguồn điểm # K40 1uCi /gps/hist/point 0.02001 /gps/hist/point 0.08331 0.061 gps/source/intensity 37000 /gps/hist/point 0.02001 0.0051 /gps/hist/point 0.08702 /gps/particle gamma /gps/hist/point 0.02950 /gps/hist/point 0.08702 0.018 /gps/pos/type Point /gps/hist/point 0.02950 0.00123 /gps/hist/point 0.09238 /gps/pos/centre 0 mm /gps/hist/point 0.043498 /gps/hist/point 0.09238 2.18 /gps/ang/type iso /gps/hist/point 0.043498 0.00186 /gps/hist/point 0.09280 /gps/hist/point 0.1135 /gps/hist/point 0.09280 2.15 gps/ene/type User /gps/hist/point 0.1135 0.0174 /gps/hist/point 0.09986 /gps/hist/type energy /gps/hist/point 0.04955 /gps/hist/point 0.09986 0.0054 /gps/hist/point 1.4602 /gps/hist/point 0.04955 0.0697 /gps/hist/point 0.10335 /gps/hist/point 1.4602 10.55 /gps/hist/point 0.06270 /gps/hist/point 0.10335 0.0032 /gps/hist/point 0.06270 0.0039 /gps/hist/point 0.105 # U238 1uCi /gps/hist/point 0.06288 /gps/hist/point 0.105 0.0114 gps/source/intensity 37000 /gps/hist/point 0.06288 0.0164 /gps/hist/point 0.11281 /gps/particle gamma /gps/hist/point 0.0633 /gps/hist/point 0.11282 0.215 /gps/pos/type Point /gps/hist/point 0.0633 3.75 /gps/hist/point 0.12546 /gps/pos/centre 0 mm /gps/hist/point 0.07392 /gps/hist/point 0.12546 0.00133 /gps/ang/type iso /gps/hist/point 0.07392 0.170 /gps/hist/point 0.13130 /gps/hist/point 0.07392 /gps/hist/point 0.13130 0.0306 gps/ene/type User /gps/hist/point 0.07392 0.017 /gps/hist/point 0.14015 /gps/hist/type energy /gps/hist/point 0.08331 /gps/hist/point 0.14015 0.00128 68 /gps/hist/point 0.15271 /gps/hist/point 0.46810 0.00228 /gps/hist/point 0.73339 /gps/hist/point 0.15271 0.0102 /gps/hist/point 0.47550 /gps/hist/point 0.73339 0.0117 /gps/hist/point 0.15948 /gps/hist/point 0.47550 0.00285 /gps/hist/point 0.73800 /gps/hist/point 0.15948 0.00111 /gps/hist/point 0.50675 /gps/hist/point 0.73800 0.00196 /gps/hist/point 0.1847 /gps/hist/point 0.50675 0.00219 /gps/hist/point 0.73995 /gps/hist/point 0.1847 0.00171 /gps/hist/point 0.50750 /gps/hist/point 0.73995 0.0110 /gps/hist/point 0.1848 /gps/hist/point 0.50750 0.0016 /gps/hist/point 0.74280 /gps/hist/point 0.1848 0.012 /gps/hist/point 0.50920 /gps/hist/point 0.74280 0.096 /gps/hist/point 0.18615 /gps/hist/point 0.50920 0.0021 /gps/hist/point 0.74281 /gps/hist/point 0.18615 0.00299 /gps/hist/point 0.51340 /gps/hist/point 0.74281 00350 /gps/hist/point 0.20097 /gps/hist/point 0.51340 0.00192 /gps/hist/point 0.75500 /gps/hist/point 0.20097 0.00153 /gps/hist/point 0.52140 /gps/hist/point 0.75500 0.00207 /gps/hist/point 0.20312 /gps/hist/point 0.52140 0.00126 /gps/hist/point 0.76030 /gps/hist/point 0.20312 0.00309 /gps/hist/point 0.54410 /gps/hist/point 0.76030 0.0016 /gps/hist/point 0.2099 /gps/hist/point 0.54410 0.0037 /gps/hist/point 0.76637 /gps/hist/point 0.2099 0.00128 /gps/hist/point 0.56520 /gps/hist/point 0.76637 0.316 /gps/hist/point 0.22650 /gps/hist/point 0.56520 0.00175 /gps/hist/point 0.78040 /gps/hist/point 0.22650 0.0071 /gps/hist/point 0.56890 /gps/hist/point 0.78040 0.00153 /gps/hist/point 0.22725 /gps/hist/point 0.56890 0.0061 /gps/hist/point 0.78137 /gps/hist/point 0.22725 0.0099 /gps/hist/point 0.56950 /gps/hist/point 0.78137 0.0075 /gps/hist/point 0.24537 /gps/hist/point 0.56950 0.0139 /gps/hist/point 0.78627 /gps/hist/point 0.24537 0.00128 /gps/hist/point 0.62566 /gps/hist/point 0.78627 0.0570 /gps/hist/point 0.24779 /gps/hist/point 0.62566 0.0014 /gps/hist/point 0.79490 /gps/hist/point 0.24779 0.0010 /gps/hist/point 0.64770 /gps/hist/point 0.79490 0.00114 /gps/hist/point 0.24922 /gps/hist/point 0.64770 0.0016 /gps/hist/point 0.79610 /gps/hist/point 0.24922 0.0043 /gps/hist/point 0.64912 /gps/hist/point 0.79610 0.0044 /gps/hist/point 0.25819 /gps/hist/point 0.64912 0.00110 /gps/hist/point 0.80410 /gps/hist/point 0.25819 0.0754 /gps/hist/point 0.65530 /gps/hist/point 0.80410 0.00105 /gps/hist/point 0.27228 /gps/hist/point 0.65530 0.0014 /gps/hist/point 0.80574 /gps/hist/point 0.27228 0.00184 /gps/hist/point 0.66650 /gps/hist/point 0.80574 0.0043 /gps/hist/point 0.29379 /gps/hist/point 0.66650 0.00197 /gps/hist/point 0.80580 /gps/hist/point 0.29379 0.0051 /gps/hist/point 0.66970 /gps/hist/point 0.80580 0.0043 /gps/hist/point 0.33040 /gps/hist/point 0.66970 0.00168 /gps/hist/point 0.80820 /gps/hist/point 0.33040 0.00131 /gps/hist/point 0.69110 /gps/hist/point 0.80820 0.0030 /gps/hist/point 0.33810 /gps/hist/point 0.69110 0.0078 /gps/hist/point 0.81820 /gps/hist/point 0.33810 0.00114 /gps/hist/point 0.69260 /gps/hist/point 0.81820 0.0010 /gps/hist/point 0.36950 /gps/hist/point 0.69260 0.00211 /gps/hist/point 0.81920 /gps/hist/point 0.36950 0.0043 /gps/hist/point 0.69550 /gps/hist/point 0.81920 0.00320 /gps/hist/point 0.37200 /gps/hist/point 0.69550 0.0016 /gps/hist/point 0.82420 /gps/hist/point 0.37200 0.00207 /gps/hist/point 0.69903 /gps/hist/point 0.82420 0.00211 /gps/hist/point 0.38794 /gps/hist/point 0.69903 0.0069 /gps/hist/point 0.82510 /gps/hist/point 0.38794 0.0014 /gps/hist/point 0.70205 /gps/hist/point 0.82510 0.00320 /gps/hist/point 0.45093 /gps/hist/point 0.70205 0.0077 /gps/hist/point 0.82560 /gps/hist/point 0.45093 0.00301 /gps/hist/point 0.70590 /gps/hist/point 0.82560 0.0014 /gps/hist/point 0.45360 /gps/hist/point 0.70590 0.0079 /gps/hist/point 0.83150 /gps/hist/point 0.45360 0.0024 /gps/hist/point 0.73090 /gps/hist/point 0.83150 0.0070 /gps/hist/point 0.45868 /gps/hist/point 0.73090 0.00107 /gps/hist/point 0.84410 /gps/hist/point 0.45868 0.00192 /gps/hist/point 0.73250 /gps/hist/point 0.84410 0.00181 /gps/hist/point 0.46810 /gps/hist/point 0.73250 0.00128 /gps/hist/point 0.85157 69 /gps/hist/point 0.85157 0.0063 /gps/hist/point 1.06190 /gps/hist/point 1.75981 0.0023 /gps/hist/point 0.86680 /gps/hist/point 1.06190 0.0020 /gps/hist/point 1.76540 /gps/hist/point 0.86680 0.00110 /gps/hist/point 1.12060 /gps/hist/point 1.76540 0.0087 /gps/hist/point 0.87600 /gps/hist/point 1.12060 0.0017 /gps/hist/point 1.80900 /gps/hist/point 0.87600 0.00428 /gps/hist/point 1.12570 /gps/hist/point 1.80900 0.0043 /gps/hist/point 0.88050 /gps/hist/point 1.12570 0.0040 /gps/hist/point 1.81969 /gps/hist/point 0.88050 0.0177 /gps/hist/point 1.17420 /gps/hist/point 1.81969 0.00114 /gps/hist/point 0.88090 /gps/hist/point 1.17420 0.00190 /gps/hist/point 1.8315 /gps/hist/point 0.88090 0.00384 /gps/hist/point 1.19377 /gps/hist/point 1.8315 0.0167 /gps/hist/point 0.88324 /gps/hist/point 1.19377 0.0131 /gps/hist/point 1.86309 /gps/hist/point 0.88324 0.0199 /gps/hist/point 1.22000 /gps/hist/point 1.86309 0.00121 /gps/hist/point 0.88728 /gps/hist/point 1.22000 0.0010 /gps/hist/point 1.86768 /gps/hist/point 0.88728 0.0074 /gps/hist/point 1.23740 /gps/hist/point 1.86768 0.0075 /gps/hist/point 0.89867 /gps/hist/point 1.23740 0.00512 /gps/hist/point 1.87485 /gps/hist/point 0.89867 0.0054 /gps/hist/point 1.35290 /gps/hist/point 1.87485 0.0078 /gps/hist/point 0.92170 /gps/hist/point 1.35290 0.00196 /gps/hist/point 1.89350 /gps/hist/point 0.92170 0.0123 /gps/hist/point 1.39270 /gps/hist/point 1.89350 0.0021 /gps/hist/point 0.92500 /gps/hist/point 1.39270 0.00156 /gps/hist/point 1.91117 /gps/hist/point 0.92500 0.0133 /gps/hist/point 1.39390 /gps/hist/point 1.91117 0.0053 /gps/hist/point 0.92600 /gps/hist/point 1.39390 0.00350 /gps/hist/point 1.93701 /gps/hist/point 0.92600 0.0031 /gps/hist/point 1.41360 /gps/hist/point 1.93701 0.0030 /gps/hist/point 0.92672 /gps/hist/point 1.41369 0.0021 /gps/hist/point 0.92672 0.0122 /gps/hist/point 1.43413 /gps/hist/point 0.92672 /gps/hist/point 1.43413 0.0083 # Th232 1uCi /gps/hist/point 0.92672 0.00124 /gps/hist/point 1.45270 gps/source/intensity 37000 /gps/hist/point 0.93630 /gps/hist/point 1.45270 0.00136 /gps/particle gamma /gps/hist/point 0.93630 0.00185 /gps/hist/point 1.45850 /gps/pos/type Point /gps/hist/point 0.94190 /gps/hist/point 1.45850 0.00185 /gps/pos/centre 0 mm /gps/hist/point 0.94190 0.0030 /gps/hist/point 1.5105 /gps/ang/type iso /gps/hist/point 0.94591 /gps/hist/point 1.5105 0.0130 /gps/hist/point 0.94591 0.0111 /gps/hist/point 1.52721 gps/ene/type User /gps/hist/point 0.94600 /gps/hist/point 1.52721 0.0021 /gps/hist/type energy /gps/hist/point 0.94600 0.0228 /gps/hist/point 1.55000 /gps/hist/point 0.00665 /gps/hist/point 0.94770 /gps/hist/point 1.55000 0.00185 /gps/hist/point 0.00665 0.0165 /gps/hist/point 0.94770 0.00275 /gps/hist/point 1.55374 /gps/hist/point 0.029373 /gps/hist/point 0.98030 /gps/hist/point 1.55384 0.0090 /gps/hist/point 0.029373 2.17 /gps/hist/point 0.98030 0.0076 /gps/hist/point 1.57068 /gps/hist/point 0.057100 /gps/hist/point 0.98160 /gps/hist/point 1.57068 0.00122 /gps/hist/point 0.057100 0.0498 /gps/hist/point 0.98160 0.00122 /gps/hist/point 1.59370 /gps/hist/point 0.063920 /gps/hist/point 0.98420 /gps/hist/point 1.59370 0.00384 /gps/hist/point 0.063920 0.0007 /gps/hist/point 0.98420 0.00275 /gps/hist/point 1.66840 /gps/hist/point 0.070490 /gps/hist/point 0.99610 /gps/hist/point 1.66840 0.00129 /gps/hist/point 0.070490 0.0007 /gps/hist/point 0.99610 0.0041 /gps/hist/point 1.69380 /gps/hist/point 0.074510 /gps/hist/point 1.00103 /gps/hist/point 1.69380 0.00117 /gps/hist/point 0.074510 0.0402 /gps/hist/point 1.00103 0.839 /gps/hist/point 1.7325 /gps/hist/point 0.086477 /gps/hist/point 1.04170 /gps/hist/point 1.7325 0.0212 /gps/hist/point 0.086477 1.843 /gps/hist/point 1.04170 0.0142 /gps/hist/point 1.73773 /gps/hist/point 0.087990 /gps/hist/point 1.05940 /gps/hist/point 1.73773 0.0202 /gps/hist/point 0.087990 0.1698 /gps/hist/point 1.05940 0.00110 /gps/hist/point 1.75981 /gps/hist/point 0.094650 70 /gps/hist/point 0.094650 0.775 /gps/hist/point 0.216540 /gps/hist/point 0.440940 0.1912 /gps/hist/point 0.105200 /gps/hist/point 0.216540 0.0130 /gps/hist/point 0.447762 /gps/hist/point 0.105200 0.041 /gps/hist/point 0.226100 /gps/hist/point 0.447762 0.1043 /gps/hist/point 0.108500 /gps/hist/point 0.226100 0.0171 /gps/hist/point 0.454200 /gps/hist/point 0.108500 0.0006 /gps/hist/point 0.237860 /gps/hist/point 0.454200 0.04 /gps/hist/point 0.115140 /gps/hist/point 0.237860 0.0019 /gps/hist/point 0.459222 /gps/hist/point 0.115140 0.003 /gps/hist/point 0.242300 /gps/hist/point 0.459222 0.989 /gps/hist/point 0.117692 /gps/hist/point 0.242300 0.0029 /gps/hist/point 0.464800 /gps/hist/point 0.117692 0.0029 /gps/hist/point 0.246140 /gps/hist/point 0.464800 0.0026 /gps/hist/point 0.131101 /gps/hist/point 0.246140 0.0041 /gps/hist/point 0.467400 /gps/hist/point 0.131101 0.0508 /gps/hist/point 0.250650 /gps/hist/point 0.467400 0.0144 /gps/hist/point 0.134285 /gps/hist/point 0.250650 0.0047 /gps/hist/point 0.473900 /gps/hist/point 0.134285 0.0018 /gps/hist/point 0.252780 /gps/hist/point 0.473900 0.0033 /gps/hist/point 0.143230 /gps/hist/point 0.252780 0.0066 /gps/hist/point 0.490800 /gps/hist/point 0.143230 0.0114 /gps/hist/point 0.257300 /gps/hist/point 0.490800 0.1078 /gps/hist/point 0.147500 /gps/hist/point 0.257300 0.0524 /gps/hist/point 0.497100 /gps/hist/point 0.147500 0.0018 /gps/hist/point 0.278700 /gps/hist/point 0.497100 0.0128 /gps/hist/point 0.151409 /gps/hist/point 0.278700 0.0047 /gps/hist/point 0.499020 /gps/hist/point 0.151409 0.0067 /gps/hist/point 0.285240 /gps/hist/point 0.499020 0.1576 /gps/hist/point 0.153490 /gps/hist/point 0.285240 0.0154 /gps/hist/point 0.505500 /gps/hist/point 0.153490 0.0407 /gps/hist/point 0.309900 /gps/hist/point 0.505500 0.0055 /gps/hist/point 0.155239 /gps/hist/point 0.309900 0.0032 /gps/hist/point 0.513400 /gps/hist/point 0.155239 0.00023 /gps/hist/point 0.316100 /gps/hist/point 0.513400 0.0133 /gps/hist/point 0.162504 /gps/hist/point 0.316100 0.0037 /gps/hist/point 0.517000 /gps/hist/point 0.162504 0.1674 /gps/hist/point 0.347640 /gps/hist/point 0.517000 0.0046 /gps/hist/point 0.162504 /gps/hist/point 0.347640 0.0145 /gps/hist/point 0.526690 /gps/hist/point 0.162504 0.16 /gps/hist/point 0.359740 /gps/hist/point 0.526690 0.0463 /gps/hist/point 0.169162 /gps/hist/point 0.359740 0.0869 /gps/hist/point 0.531800 /gps/hist/point 0.169162 0.251 /gps/hist/point 0.361285 /gps/hist/point 0.531800 0.0070 /gps/hist/point 0.170600 /gps/hist/point 0.361285 0.0218 /gps/hist/point 0.552210 /gps/hist/point 0.170600 0.507 /gps/hist/point 0.367920 /gps/hist/point 0.552210 0.0165 /gps/hist/point 0.179050 /gps/hist/point 0.367920 0.0037 /gps/hist/point 0.553700 /gps/hist/point 0.179050 0.0278 /gps/hist/point 0.377270 /gps/hist/point 0.553700 0.0030 /gps/hist/point 0.180760 /gps/hist/point 0.377270 0.0275 /gps/hist/point 0.554900 /gps/hist/point 0.180760 0.00011 /gps/hist/point 0.383500 /gps/hist/point 0.554900 0.0031 /gps/hist/point 0.186800 /gps/hist/point 0.383500 0.0019 /gps/hist/point 0.562930 /gps/hist/point 0.186800 0.0209 /gps/hist/point 0.398800 /gps/hist/point 0.562930 0.0545 /gps/hist/point 0.190552 /gps/hist/point 0.398800 0.0111 /gps/hist/point 0.573700 /gps/hist/point 0.190552 0.0861 /gps/hist/point 0.408800 /gps/hist/point 0.573700 0.0332 /gps/hist/point 0.194970 /gps/hist/point 0.408800 0.0005 /gps/hist/point 0.578700 /gps/hist/point 0.194970 0.1073 /gps/hist/point 0.412500 /gps/hist/point 0.578700 0.0017 /gps/hist/point 0.201620 /gps/hist/point 0.412500 0.0083 /gps/hist/point 0.583200 /gps/hist/point 0.201620 0.0221 /gps/hist/point 0.418400 /gps/hist/point 0.583200 0.0016 /gps/hist/point 0.210670 /gps/hist/point 0.418400 0.0091 /gps/hist/point 0.595390 /gps/hist/point 0.210670 0.0178 /gps/hist/point 0.430900 /gps/hist/point 0.595390 0.1178 /gps/hist/point 0.211300 /gps/hist/point 0.430900 0.0178 /gps/hist/point 0.599300 /gps/hist/point 0.211300 0.0202 /gps/hist/point 0.433200 /gps/hist/point 0.599300 0.0294 /gps/hist/point 0.212340 /gps/hist/point 0.433200 0.0117 /gps/hist/point 0.610000 /gps/hist/point 0.212340 0.0065 /gps/hist/point 0.440940 /gps/hist/point 0.610000 0.0567 71 /gps/hist/point 0.642400 /gps/hist/point 0.774000 /gps/hist/point 0.942800 /gps/hist/point 0.642400 0.0202 /gps/hist/point 0.774000 0.0108 /gps/hist/point 0.942800 0.0019 /gps/hist/point 0.663300 /gps/hist/point 0.783200 /gps/hist/point 0.948300 /gps/hist/point 0.663300 0.0037 /gps/hist/point 0.783200 0.0056 /gps/hist/point 0.948300 0.006 /gps/hist/point 0.669900 /gps/hist/point 0.784200 /gps/hist/point 0.955000 /gps/hist/point 0.669900 0.0018 /gps/hist/point 0.784200 0.0022 /gps/hist/point 0.955000 0.0002 /gps/hist/point 0.669901 /gps/hist/point 0.805000 /gps/hist/point 0.960800 /gps/hist/point 0.669901 0.504 /gps/hist/point 0.805000 0.0214 /gps/hist/point 0.960800 0.0041 /gps/hist/point 0.678040 /gps/hist/point 0.806400 /gps/hist/point 0.962800 /gps/hist/point 0.678040 0.0647 /gps/hist/point 0.806400 0.0123 /gps/hist/point 0.962800 0.0015 /gps/hist/point 0.681200 /gps/hist/point 0.811600 /gps/hist/point 0.968200 /gps/hist/point 0.681200 0.0143 /gps/hist/point 0.811600 0.0060 /gps/hist/point 0.968200 0.0083 /gps/hist/point 0.690000 /gps/hist/point 0.815900 /gps/hist/point 0.978200 /gps/hist/point 0.690000 0.0021 /gps/hist/point 0.815900 0.0195 /gps/hist/point 0.978200 0.0058 /gps/hist/point 0.698500 /gps/hist/point 0.817000 /gps/hist/point 0.984800 /gps/hist/point 0.698500 0.0106 /gps/hist/point 0.817000 0.0095 /gps/hist/point 0.984800 0.0102 /gps/hist/point 0.703700 /gps/hist/point 0.832000 /gps/hist/point 0.994000 /gps/hist/point 0.703700 0.0091 /gps/hist/point 0.832000 0.0075 /gps/hist/point 0.994000 0.0006 /gps/hist/point 0.707800 /gps/hist/point 0.846800 /gps/hist/point 1.001000 /gps/hist/point 0.707800 0.0091 /gps/hist/point 0.846800 0.0013 /gps/hist/point 1.001000 0.0008 /gps/hist/point 0.717000 /gps/hist/point 0.849500 /gps/hist/point 1.007000 /gps/hist/point 0.717000 0.0421 /gps/hist/point 0.849500 0.0039 /gps/hist/point 1.007000 0.0014 /gps/hist/point 0.725100 /gps/hist/point 0.870700 /gps/hist/point 1.011000 /gps/hist/point 0.725100 0.0633 /gps/hist/point 0.870700 0.0031 /gps/hist/point 1.011000 0.0019 /gps/hist/point 0.727800 /gps/hist/point 0.874000 /gps/hist/point 1.026500 /gps/hist/point 0.727800 0.0029 /gps/hist/point 0.874000 0.00120 /gps/hist/point 1.026500 0.0075 /gps/hist/point 0.741100 /gps/hist/point 0.880900 /gps/hist/point 1.092500 /gps/hist/point 0.741100 0.0236 /gps/hist/point 0.880900 0.0097 /gps/hist/point 1.092500 0.006 /gps/hist/point 0.744900 /gps/hist/point 0.890100 /gps/hist/point 1.132100 /gps/hist/point 0.744900 0.0053 /gps/hist/point 0.890100 0.1052 /gps/hist/point 1.132100 0.0006 /gps/hist/point 0.751600 /gps/hist/point 0.898300 /gps/hist/point 1.139100 /gps/hist/point 0.751600 0.0023 /gps/hist/point 0.898300 0.0022 /gps/hist/point 1.139100 0.0004 /gps/hist/point 0.757900 /gps/hist/point 0.918900 /gps/hist/point 1.144000 /gps/hist/point 0.757900 0.0324 /gps/hist/point 0.918900 0.006 /gps/hist/point 1.144000 0.0027 /gps/hist/point 0.764550 /gps/hist/point 0.935200 /gps/hist/point 1.201000 /gps/hist/point 0.764550 0.0891 /gps/hist/point 0.935200 0.0369 /gps/hist/point 1.201000 0.006 /gps/hist/point 0.767500 /gps/hist/point 0.941900 /gps/hist/point 0.767500 0.0032 /gps/hist/point 0.941900 0.0048 Nguồn có kích thước /gps/pos/centre 0 mm gps/source/intensity 37000 /gps/ang/type iso /gps/particle gamma gps/ene/type User /gps/pos/type Plane /gps/hist/type energy /gps/pos/shape Circle /gps/hist/point … /gps/pos/halfx 50 mm /gps/pos/halfy 50 mm Nguồn muon /gps/pos/halfz 10 mm gps/source/intensity 37000 /gps/pos/radius 50 mm /gps/particle mu- /gps/pos/inner_radius mm /gps/pos/type Beam 72 /gps/pos/shape Circle /gps/ang/type beam2d /gps/pos/radius mm /gps/ene/type Mono /gps/pos/sigma_r mm /gps/ene/mono 450 MeV /gps/pos/centre 0 -500 mm /gps/direction 0 73 PHỤ LỤC - CODE TÁC ĐỘNG ĐỘ PHÂN GIẢI RESOLUTION #include #include #include #include #include using namespace std; const double PI = 3.14159265358979323846; const int evtNb = 50000000; float edep_stdgaus(void); float edep_resolution(float &edep,float &edepstdgaus); int main () { float edep; float edepstdgaus, edepresolution; ofstream filesave; char input[100]; char output[100]; cout input; ifstream openfile(input); cout output; filesave.open (output); /* Initialize random seed: */ srand ((unsigned int)time(NULL)); // Open energy deposit save file int numberofentries = 0; for (int i=0;i> edep; } else coutx; if (!in.good()) break; h1->Fill(x); } // Ve histogram h1->Draw(); in.close(); } 78 ... TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TP HỒ CHÍ MINH Nguyễn Thị Bích Hậu KHẢO SÁT ĐÁP ỨNG CỦA DETECTOR HPGe CHO PHĨNG XẠ MƠI TRƯỜNG BẰNG PHẦN MỀM GEANT4 Chuyên ngành: Vật lí nguyên tử Mã số: 60.44.01.06 LUẬN... Agren G., Johansson L (2008); … Trong luận văn này, thực mơ khả đáp ứng đầu dị HPGe cho cho số đồng vị phóng xạ có môi trường xạ vũ trụ muon Cụ thể khảo sát cho đồng vị U-238, Th-232, K-40, với... 4,5.109 99,274 Chuỗi phóng xạ U-238 chia thành chuỗi phóng xạ hoạt tính phóng xạ nhân phóng xạ đầu chuỗi chi phối hoạt tính phóng xạ nhân phóng xạ khác có chuỗi Đó chuỗi phóng xạ: U-238 > U-234;