geant4系列講座二(編輯修改稿)
2024-11-17 04:42 本頁面
【文章內(nèi)容簡介】 pe( name=“U235“, iz=92, n=235, a=235.01*g/mole)。 G4Isotope* U8 = new G4Isotope( name=“U238“, iz=92, n=238, a=238.03*g/mole)。 G4Element* elU = new G4Element( name=“enriched Uranium”, symbol=“U”, ncomponents )。 elUAddIsotope(U5, abundance= 90.*perCent)。 elUAddIsotope(U8, abundance= 10.*perCent)。,第十一頁，共二十九頁。,三、Geant4中使用物理量的單位,基本單位： Geant4中用戶可以為選定的物理量選擇各種單位，但是 Geant4內(nèi)核內(nèi)定義了如下基本單位： millimeter (mm) nanosecond (ns) Mega electron Volt (MeV) positron charge (eplus) degree Kelvin (kelvin) the amount of substance (mole) luminous intensity (candela) radian (radian) steradian (steradian) 其它所有單位都以上述單位為基礎(chǔ)得到，如： Millimeter＝ mm ＝1； Meter ＝ m ＝1000*mm；,第十二頁，共二十九頁。,用戶輸入數(shù)據(jù)單位： Geant4中用戶輸入的數(shù)據(jù)必須帶有單位（系統(tǒng)默認(rèn)強(qiáng)烈建議不使用），如： G4double Size = 15*km, KineticEnergy = 90.3*GeV, density = 11*mg/cm3。 同樣，如果數(shù)據(jù)為數(shù)組格式或者從文件讀入，也必須帶有單位，如： for (int j=0, jjmax, j++) CrossSection[j] *= millibarn。 數(shù)據(jù)輸出單位： 帶單位數(shù)據(jù)輸出格式如： G4cout KineticEnergy/keV “ keV“。 G4cout density/(g/cm3) “ g/cm3“。 或者給出Geant4選擇的最優(yōu)化單位： G4cout G4BestUnit(StepSize, “Length“)。,第十三頁，共二十九頁。,四、幾何結(jié)構(gòu)定義,Geant4中以“體”（Volume）的概念定義幾何形狀。模型中最大的“體”稱為“世界體”（World Volume），其它的體都位于世界體內(nèi)部?！笆澜珞w”內(nèi)的“體”之間是包含和被包含關(guān)系，被包含的稱包含它的體為母體（mother Volume）。 Geant4中的體的定義包含三個(gè)層次： 1.定義幾何形狀； 2.定義物理屬性； 3.定義所在母體的位置；,第十四頁，共二十九頁。,定義幾何形狀（solid）： G4double expHall_x = 30.0*m。 G4double expHall_y = 40.0*m。 G4double expHall_z = 60.0*m。 G4Box* experimentalHall_box = new G4Box(“expHall_box“,expHall_x,expHall_y,expHall_z)。 指定幾何形狀的名稱、尺寸等信息。,第十五頁，共二十九頁。,定義物理屬性： G4LogicalVolume* experimentalHall_log = new G4LogicalVolume( experimentalHall_box,Ar,“expHall_log“)。 給定物理屬性指向的幾何體，該幾何體對應(yīng)的物 質(zhì)，該幾何體的名稱等信息。 G4LogicalVolume* tracker_log = new G4LogicalVolume(tracker_tube,Al,“tracker_log“)。 ,第十六頁，共二十九頁。,定義物理位置 G4double trackerPos_x = 1.0*meter。 G4double trackerPos_y = 0.0*meter。 G4double trackerPos_z = 0.0*meter。 G4VPhysicalVolume* tracker_phys = new G4PVPlacement(0, // no rotation G4ThreeVector(trackerPos_x,trackerPos_y,trackerPos_z), // translation position tracker_log, // its logical volume “tracker“, // its
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1