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, 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 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ể 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 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ỹ đạ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ự 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 [...]... các quá trình biến đổi môi trường gây ra sự mất cân bằng phóng xạ thì các chuỗi phóng xạ này thường có cân bằng phóng xạ Điều này cũng đồng nghĩa với hoạt độ phóng xạ của các nhân phóng xạ có trong mỗi chuỗi là bằng nhau và bằng với hoạt độ phóng xạ của nhân bắt đầu mỗi chuỗi Bảng 1.1 Độ giàu đồng vị của các nhân phóng xạ nguyên thuỷ [9] Nhân phóng xạ Thời gian bán huỷ (năm) Độ giàu đồng vị (%) K-40... đối của electron ne là số electron trên một đơn vị thể tích của môi trường ne = ZρNa/A Na là số Avogadro ρ là mật độ của môi trường vật chất 15 (1.22) A là số khối của môi trường Z là số bậc nguyên tử của môi trường m0 là khối lượng của electron v là vận tốc của hạt I là năng lượng ion hóa trung bình Sự phát bức xạ hãm: Khi hạt mang điện đi vào môi trường vật chất, do ảnh hưởng của trường Coulomb của. .. Geant4 được sử dụng một cách linh hoạt tương ứng với mục đích sử dụng của người dùng Trong chương trình này, chúng tôi tập trung phát triển chương trình mô phỏng tương tác của bức xạ gamma phát ra bởi các đồng vị phóng xạ môi trường gồm U-238, Th-232, K-40, tương tác của hạt muon với bán dẫn siêu tinh khiết germanium nhằm nghiên cứu đáp ứng của detector HPGe siêu tinh khiết bằng phần mềm mô phỏng Geant4. .. lập được nguyên tố phóng xạ nhân tạo đầu tiên Po-210 Người ta nhận thấy rằng tính phóng xạ của nguyên tố phóng xạ tự nhiên và nguyên tố phóng xạ nhân tạo là như nhau 16 1.3.1 Các chuỗi phóng xạ nguyên thủy trong tự nhiên Khắp nơi trên trái đất của chúng ta đâu đâu cũng có phóng xạ từ nước, không khí, đất, động vật và thực vật đến cơ thể người Phóng xạ tự nhiên: Là những đồng vị phóng xạ tồn tại trong... 7,1.108 0,72 U-238 4,5.109 99,274 Chuỗi phóng xạ U-238 có thể chia thành các chuỗi phóng xạ con trong đó hoạt tính phóng xạ của nhân phóng xạ đầu chuỗi sẽ chi phối hoạt tính phóng xạ của các nhân phóng xạ khác có trong chuỗi Đó là các chuỗi phóng xạ: U-238 > U-234; Th-230; Ra-226; Rn-222 > Po- 214 và Pb-210 >Po-210 Uranium khá phổ biến trong tự nhiên, về độ giàu nó ứng hàng thứ 38 trong số các nguyên... electron Z số bậc nguyên tử của môi trường 1.3 Giới thiệu phóng xạ môi trường Từ năm 1895 với sự phát hiện ra tia X của nhà bác học người Đức W.Roentgen, và sau khám phá của nhà bác học Henri Becquerel về sự phát xạ tự nhiên đã tạo tiền đề cho sự phát hiện ra nguyên tố phóng xạ Polonium của ông bà Piere Curie và Marie Curie năm 1898, và tiếp theo sau đó là sự tìm ra nguyên tố phóng xạ thứ hai tồn tại trong... trào K-40 Có trong đất = 1-30 pCi/g (0,037-1,1 Bq/g) Từ việc phân tích độ giàu của các đồng vị phóng xạ có trong môi trường, trong luận văn này, chúng tôi đã chọn nghiên cứu trên các đồng vị phóng xạ phổ biến của môi trường là U238, Th-232, K-40 và bức xạ vũ trụ muon 1.3.2 Phóng xạ bức xạ vũ trụ muon Tùy thuộc vào nguồn gốc của chúng mà tia vũ trụ được chia thành hai loại, tia vũ trụ sơ cấp và tia vũ... nghỉ của electron Góc tán xạ θ của electron sau tán xạ liên hệ với góc tán xạ của gamma ϕ như sau: cos θ = (1 + γ ) tan ϕ (1.5) 2 Theo (1.3) và (1.4) góc tán xạ của gamma sau tán xạ càng lớn thì năng lượng của nó E’ càng bé Nghĩa là gamma càng mất nhiều năng lượng Gamma chuyển phần năng lượng lớn nhất cho electron sau tán xạ bay ra một góc 1800, tức là tán xạ giật lùi Khi đó năng lượng giật lùi của. .. trong tự nhiên là Radium của bà Marie Curie Từ đó, việc nghiên cứu về phóng xạ được mở rộng Năm 1899, Rutherford đã nhận thấy một phần tia phóng xạ bị lệch khi xuyên qua từ trường của hạt nhân Heli ( α ) và electron ( β − ) Năm 1903, Rutherford và Soddy kết hợp giữa sự phân rã ( α ) và ( β − ) với sự thay đổi bậc số của nguyên tố phóng xạ đã thiết lập định luật phân rã phóng xạ Năm 1934, Irene và Frederic... phóng xạ, còn U-238, Th-232, U-235 có sản phẩm phân rã là hạt nhân không bền, tiếp tục phân rã đến khi hạt nhân con cháu cuối cùng là hạt nhân bền, và sự phân rã liên tiếp của các đồng vị trong một họ đó tạo thành chuỗi phóng xạ tự nhiên Ngày nay có ba chuỗi phóng xạ tự nhiên phổ biến mà U238, Th-232, U-235 là các nhân bắt đầu của mỗi chuỗi Trong môi trường nếu không có các quá trình biến đổi môi trường ...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;... 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