Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống
1
/ 80 trang
THÔNG TIN TÀI LIỆU
Thông tin cơ bản
Định dạng
Số trang
80
Dung lượng
1,42 MB
Nội dung
luan van thac si su pham, thac si giao duc ,ths xa hoi1 of 138 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 luan van thac si su pham, thac si giao duc ,ths xa hoi1 of 138 luan van thac si su pham, thac si giao duc ,ths xa hoi2 of 138 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 luan van thac si su pham, thac si giao duc ,ths xa hoi2 of 138 luan van thac si su pham, thac si giao duc ,ths xa hoi3 of 138 LỜI CẢM ƠN Trong suốt trình thực luận văn, 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 luan van thac si su pham, thac si giao duc ,ths xa hoi3 of1 138 luan van thac si su pham, thac si giao duc ,ths xa hoi4 of 138 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 luan van thac si su pham, thac si giao duc ,ths xa hoi4 of2 138 luan van thac si su pham, thac si giao duc ,ths xa hoi5 of 138 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 luan van thac si su pham, thac si giao duc ,ths xa hoi5 of3 138 luan van thac si su pham, thac si giao duc ,ths xa hoi6 of 138 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 luan van thac si su pham, thac si giao duc ,ths xa hoi6 of4 138 luan van thac si su pham, thac si giao duc ,ths xa hoi7 of 138 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 toà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 toà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 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ể luan van thac si su pham, thac si giao duc ,ths xa hoi7 of5 138 luan van thac si su pham, thac si giao duc ,ths xa hoi8 of 138 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ò 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 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 vuông 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 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 luan van thac si su pham, thac si giao duc ,ths xa hoi8 of6 138 luan van thac si su pham, thac si giao duc ,ths xa hoi9 of 138 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 toà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ỹ luan van thac si su pham, thac si giao duc ,ths xa hoi9 of7 138 luan van thac si su pham, thac si giao duc ,ths xa hoi10 of 138 đạ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 thoát tia X đóng góp vào số đếm đỉnh lượng toà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ự luan van thac si su pham, thac si giao duc ,ths xa hoi10 of 138 luan van thac si su pham, thac si giao duc ,ths xa hoi66 of 138 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; 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; }; #endif luan van thac si su pham, thac si giao duc ,ths xa hoi6664 of 138 luan van thac si su pham, thac si giao duc ,ths xa hoi67 of 138 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() { // 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(); } 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" // 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(); luan van thac si su pham, thac si giao duc ,ths xa hoi6765 of 138 luan van thac si su pham, thac si giao duc ,ths xa hoi68 of 138 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-") { //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); // 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); } 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(); } PHYSICSLIST.HH #ifndef PhysicsList_h #define PhysicsList_h #include "G4VUserPhysicsList.hh" #include "globals.hh" class PhysicsList: public G4VUserPhysicsList { public: PhysicsList(); luan van thac si su pham, thac si giao duc ,ths xa hoi6866 of 138 luan van thac si su pham, thac si giao duc ,ths xa hoi69 of 138 ~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 luan van thac si su pham, thac si giao duc ,ths xa hoi6967 of 138 luan van thac si su pham, thac si giao duc ,ths xa hoi70 of 138 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 GPSPrimaryGeneratorAction::GeneratePrimaries(G4Event* G4GeneralParticleSource* particleGun; 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 luan van thac si su pham, thac si giao duc ,ths xa hoi7068 of 138 luan van thac si su pham, thac si giao duc ,ths xa hoi71 of 138 /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 luan van thac si su pham, thac si giao duc ,ths xa hoi7169 of 138 luan van thac si su pham, thac si giao duc ,ths xa hoi72 of 138 /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 luan van thac si su pham, thac si giao duc ,ths xa hoi7270 of 138 luan van thac si su pham, thac si giao duc ,ths xa hoi73 of 138 /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 luan van thac si su pham, thac si giao duc ,ths xa hoi7371 of 138 luan van thac si su pham, thac si giao duc ,ths xa hoi74 of 138 /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 luan van thac si su pham, thac si giao duc ,ths xa hoi7472 of 138 luan van thac si su pham, thac si giao duc ,ths xa hoi75 of 138 /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 luan van thac si su pham, thac si giao duc ,ths xa hoi7573 of 138 luan van thac si su pham, thac si giao duc ,ths xa hoi76 of 138 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(); } luan van thac si su pham, thac si giao duc ,ths xa hoi8078 of 138 ... of 138 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... 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;... thường có cân phóng xạ Điều đồng nghĩa với hoạt độ phóng xạ nhân phóng xạ có chuỗi với hoạt độ phóng xạ nhân bắt đầu chuỗi Bảng 1.1 Độ giàu đồng vị nhân phóng xạ nguyên thuỷ [9] Nhân phóng xạ Thời