# flair project file Version: 300 Mode: fluka Title: Mouse Irradiator Manual Half Dome md5: 0082c17264e8ffc74104a545574f6684 Input: TITLE Mouse Irradiator Manual Half Dome GLOBAL OLDFLAIR *Seed Wall Thickness cm #define seedWTh 0.006 *Seed Radius cm #define seedR 0.04 *Seed Length cm #define seedL 0.45 *Seed Cavity L cm #define seedCL 0.35 *Seed Tube Wall Thickness #define seedTWth 0.04 *Seed Silver Rod L cm #define AgRodL 0.3 *Seed Silver Rod Radius cm #define AgRodR 0.025 *Radioactive Layer Thickness #define srcTh 1e-07 *Target Width #define TargW 2.5 *Target Heigth #define TargH 2.0 *Target Length #define TargL 9.0 *Table Width #define TableW TargW*1.25 *Table Length #define TableL TargL+0.1 *Table Thickness #define TableTh 5*mm *Table Movement Width #define TargMovW 2*TableW *Table Movement Length #define TargMovL 2*TableL *Table Movement Heigth #define TargMovH TargH *Collimator Thickness #define collTh 3.0*cm *Collimator Hole Radius #define collR 1.*mm *Collimator Hole Longitudinal Spacing Gap #define collHgap seedL/2.0 *Collimator Inner Surface Radius #define collinR TargMovW/2.0+2*cm *Collimator Outer Surface Radius #define colloutR collinR+collTh *Collimator Shield Thickness #define collShTh 0.5*cm *Collimator Opening Clearance #define collOCl 1*mm *Collimator Longitudinal Body Half Angle #define collLBhA acos((collOCl+TargMovW/2.0)/collinR) *Collimator Opening Half Angle #define collOAng 'pi/2 - collLBhA' *Collimator Hole Long Spacing Angle In #define colHLAin '2*asin(2*collR/(2*collinR)) + 2*asin(collHgap/(2*collinR))' *Collimator Longitudinal Body Half Angle End Offset #define collLBEO 'collLBhA - 4.0*deg' *Collimator Longitudinal Hole Spacing Angle Out #define colLHSan '2*asin(seedL/(2*colloutR)) + 2*asin(collHgap/(2*colloutR))' *Collimator Longitudinal Hole Number #define LngHnum 'int( 2*collLBhA / ( 2*colLHSan ) )' *Collimator Transverse Hole Spacing Angle #define colTHSan colLHSan *Collimator Transverse Hole Full Angle #define colTHFan pi *Collimator Transverse Hole Number #define TrnsHnum 'int( colTHFan / ( 2*colTHSan ) )' *Seed Y Position Shift #define seedYshf colloutR+seedTWth+seedR *Iodine seed Strength in U = 1 cGy cm^2/h #define SeedU 0.5 *Iodine 125 Half life #define t12I125 59.49*day *Iodine 125 molecular mass in g/mol #define mI125 125.9126 *Iodine 125 Specific Activity in Bq/g (1 Bq = decay or disintegration/second ) #define SaI125 642e12 *Iodine 125 Source Strength U to Activity Bq Conversion #define I125UtoA 2.912E+7 *Iodine seed Activity (dps) = (Bq) #define I125A SeedU*I125UtoA *Phosphorus 32 Half life #define t12P32 14.3*day *Phosphorus 32 molecular mass in g/mol #define mP32 31.973907 *Phosphorus 325 Specific Activity in Bq/g (1 Bq = decay or disintegration/second ) #define SaP32 2.85518e5/2.7e-11 *Phosphorus 32 Source Strength U to Activity Bq Conversion #define P32UtoA 1 *Phosphorus seed Activity (dps) = (Bq) #define P32A SeedU*P32UtoA *Seed Longitudinal Position (-nL to nL) #define nL 0 *Seed Transverse Position (-nT to nT) #define nT 0 *Use Phosphorus 32 #define P32 *Use Iodine 125 #define I125 *Particle Tracks #define tracks *Do BDX Crossing Scoring #define ScoreBDX *Set the defaults for precision simulations DEFAULTS PRECISIO *Define the beam characteristics BEAM ISOTOPE @select=8 #if P32 *Phosphorus 32 HI-PROPE , 15. 32. *Do Iodine 125 by Default #else HI-PROPE , 53. 125. #endif *Define the beam position BEAMPOS , =-seedYshf*cos(-nL*colLHSan)*sin(-nT*colTHSan) =seedYshf*cos(-nL*colLHSan)*cos(-nT*colTHSan) \ =seedYshf*sin(-nL*colLHSan) *Define the beam position BEAMPOS CYLI-VOL =AgRodR-srcTh =AgRodR+0 =AgRodL-2*srcTh =AgRodL+0 BEAMAXES , =cos(-nT*colTHSan) =sin(-nT*colTHSan) 0 =sin(-nL*colLHSan)*sin(-nT*colTHSan) \ =-sin(-nL*colLHSan)*cos(-nT*colTHSan) =cos(-nL*colLHSan) ROT-DEFI seedLPos 100,,,, =seedYshf+0 ROT-DEFI seedLPos 100,, =nL*colLHSan/deg,, 0. ROT-DEFI seedPos 100,,,, =seedYshf+0 ROT-DEFI seedPos 100,, =nL*colLHSan/deg,, 0. ROT-DEFI seedPos 300,, =nT*colTHSan/deg,, 0. ROT-DEFI seedTPln 300,, =nT*colTHSan/deg,, 0. ROT-DEFI Tm9Rot,,, =-9*colTHSan/deg ROT-DEFI Tm8Rot,,, =-8*colTHSan/deg ROT-DEFI Tm7Rot,,, =-7*colTHSan/deg ROT-DEFI Tm6Rot,,, =-6*colTHSan/deg ROT-DEFI Tm5Rot,,, =-5*colTHSan/deg ROT-DEFI Tm4Rot,,, =-4*colTHSan/deg ROT-DEFI Tm3Rot,,, =-3*colTHSan/deg ROT-DEFI Tm2Rot,,, =-2*colTHSan/deg ROT-DEFI Tm1Rot,,, =-1*colTHSan/deg ROT-DEFI Tp1Rot,,, =1*colTHSan/deg ROT-DEFI Tp2Rot,,, =2*colTHSan/deg ROT-DEFI Tp3Rot,,, =3*colTHSan/deg ROT-DEFI Tp4Rot,,, =4*colTHSan/deg ROT-DEFI Tp5Rot,,, =5*colTHSan/deg ROT-DEFI Tp6Rot,,, =6*colTHSan/deg ROT-DEFI Tp7Rot,,, =7*colTHSan/deg ROT-DEFI Tp8Rot,,, =8*colTHSan/deg ROT-DEFI Tp9Rot,,, =9*colTHSan/deg GEOBEGIN COMBNAME RPP airbod -50. 50. -50. 50 -50. 50 *Void sphere RPP void '=body(airbod,1)-1' '=body(airbod,2)+1' '=body(airbod,3)-1' '=body(airbod,4)+1' \ '=body(airbod,5)-1' '=body(airbod,6)+1' *Black body RPP blkbody '=body(void,1)-1' '=body(void,2)+1' '=body(void,3)-1' '=body(void,4)+1' \ '=body(void,5)-1' '=body(void,6)+1' RPP clipbody -50. 0. 0. 50 -50. 50 RPP clipbox -20. 0. -10. 10. 0.0 20. YZP clippln 0.0 RPP targvol =-TargMovW/2.0 =TargMovW/2.0 =-(TargMovH)/2.0 =(TargMovH)/2.0 =-(TargMovL)/2.0 \ =(TargMovL)/2.0 RPP targbody =-TargW/2.0 =TargW/2.0 =-TargH/2.0 =TargH/2.0 =-TargL/2.0 =TargL/2.0 $start_transform seedLPos *Seed SPH seedendL 0.0 0.0 =-(seedL-2*seedR)/2 =seedR+0 *Seed XYP seedLpla =-(seedL-2*seedR)/2 *Seed SPH seedendR 0.0 0.0 =(seedL-2*seedR)/2 =seedR+0 *Seed XYP seedRpla =(seedL-2*seedR)/2 *Seed RCC seedbody 0.0 0.0 =-(seedL-2*seedR)/2 0.0 0.0 =(seedL-2*seedR) =seedR+0 *Seed RCC seedcav 0.0 0.0 =-seedCL/2 0.0 0.0 =seedCL+0 =seedR-seedWTh *Seed RCC AgRdBdy 0.0 0.0 =-AgRodL/2 0.0 0.0 =AgRodL =AgRodR+0 $end_transform PLA LLpln =cos(1*colTHSan/2.) =sin(1*colTHSan/2.) 0 0 0 0 PLA LRpln =cos(-1*colTHSan/2.) =sin(-1*colTHSan/2.) 0 0 0 0 *Collimator SPH collbody 0 0 0 =collinR+collTh+collShTh *Collimator SPH collinSf 0 0 0 =collinR+0 *Collimator RCC colTbIn 0 0 =-TargMovL/2.0 0 0 =TargMovL+0 =TargMovW/2.0 *Collimator RCC colTbOut 0 0 =-(TargMovL+2.0*collShTh)/2.0 0 0 =TargMovL+2.0*collShTh =TargMovW/2.0+collShTh *Collimator XZP colCtPln =-(TargMovH)/2.0 *Collimator Plate RPP collPlte =-(TargMovL+2.0*collShTh)/2.0-1 =(TargMovL+2.0*collShTh)/2.0+1 \ =-(TargMovH)/2.0-collShTh =-(TargMovH)/2.0 =-(TargMovL+2.0*collShTh)/2.0-1 \ =(TargMovL+2.0*collShTh)/2.0+1 *Collimator Channel SPH ChnlTop 0 0 0 =collinR+collTh+2*seedTWth+2*seedR *Collimator Channel SPH ChnlBot 0 0 0 =collinR+collTh *Collimator Channel TRC ChEndL 0 0 0 0 0 -10.0 0.0 =10.0/tan(collLBEO) *Collimator Channel TRC ChEndR 0 0 0 0 0 10.0 0 =10.0/tan(collLBEO) *Collimator Channel RPP ChCt0 =-(2*seedTWth+2*seedR)/2.0 =(2*seedTWth+2*seedR)/2.0 =-(collinR+collTh+collShTh+5) \ =collinR+collTh+collShTh+5 =-(collinR+collTh+collShTh+5) \ =collinR+collTh+collShTh+5 RCC Lm9T0 0 =(collinR-seedTWth/2.0)*cos(-9*colLHSan) =(collinR-seedTWth/2.0)*sin(-9*colLHSan) \ 0 =(collTh+seedTWth)*cos(-9*colLHSan) =(collTh+seedTWth)*sin(-9*colLHSan) \ =collR+0 RCC Lm8T0 0 =(collinR-seedTWth/2.0)*cos(-8*colLHSan) =(collinR-seedTWth/2.0)*sin(-8*colLHSan) \ 0 =(collTh+seedTWth)*cos(-8*colLHSan) =(collTh+seedTWth)*sin(-8*colLHSan) \ =collR+0 RCC Lm7T0 0 =(collinR-seedTWth/2.0)*cos(-7*colLHSan) =(collinR-seedTWth/2.0)*sin(-7*colLHSan) \ 0 =(collTh+seedTWth)*cos(-7*colLHSan) =(collTh+seedTWth)*sin(-7*colLHSan) \ =collR+0 RCC Lm6T0 0 =(collinR-seedTWth/2.0)*cos(-6*colLHSan) =(collinR-seedTWth/2.0)*sin(-6*colLHSan) \ 0 =(collTh+seedTWth)*cos(-6*colLHSan) =(collTh+seedTWth)*sin(-6*colLHSan) \ =collR+0 RCC Lm5T0 0 =(collinR-seedTWth/2.0)*cos(-5*colLHSan) =(collinR-seedTWth/2.0)*sin(-5*colLHSan) \ 0 =(collTh+seedTWth)*cos(-5*colLHSan) =(collTh+seedTWth)*sin(-5*colLHSan) \ =collR+0 RCC Lm4T0 0 =(collinR-seedTWth/2.0)*cos(-4*colLHSan) =(collinR-seedTWth/2.0)*sin(-4*colLHSan) \ 0 =(collTh+seedTWth)*cos(-4*colLHSan) =(collTh+seedTWth)*sin(-4*colLHSan) \ =collR+0 RCC Lm3T0 0 =(collinR-seedTWth/2.0)*cos(-3*colLHSan) =(collinR-seedTWth/2.0)*sin(-3*colLHSan) \ 0 =(collTh+seedTWth)*cos(-3*colLHSan) =(collTh+seedTWth)*sin(-3*colLHSan) \ =collR+0 RCC Lm2T0 0 =(collinR-seedTWth/2.0)*cos(-2*colLHSan) =(collinR-seedTWth/2.0)*sin(-2*colLHSan) \ 0 =(collTh+seedTWth)*cos(-2*colLHSan) =(collTh+seedTWth)*sin(-2*colLHSan) \ =collR+0 RCC Lm1T0 0 =(collinR-seedTWth/2.0)*cos(-1*colLHSan) =(collinR-seedTWth/2.0)*sin(-1*colLHSan) \ 0 =(collTh+seedTWth)*cos(-1*colLHSan) =(collTh+seedTWth)*sin(-1*colLHSan) \ =collR+0 RCC L0T0 0 =(collinR-seedTWth/2.0) 0 0 =collTh+seedTWth 0 =collR+0 RCC Lp1T0 0 =(collinR-seedTWth/2.0)*cos(1*colLHSan) =(collinR-seedTWth/2.0)*sin(1*colLHSan) \ 0 =(collTh+seedTWth)*cos(1*colLHSan) =(collTh+seedTWth)*sin(1*colLHSan) \ =collR+0 RCC Lp2T0 0 =(collinR-seedTWth/2.0)*cos(2*colLHSan) =(collinR-seedTWth/2.0)*sin(2*colLHSan) \ 0 =(collTh+seedTWth)*cos(2*colLHSan) =(collTh+seedTWth)*sin(2*colLHSan) \ =collR+0 RCC Lp3T0 0 =(collinR-seedTWth/2.0)*cos(3*colLHSan) =(collinR-seedTWth/2.0)*sin(3*colLHSan) \ 0 =(collTh+seedTWth)*cos(3*colLHSan) =(collTh+seedTWth)*sin(3*colLHSan) \ =collR+0 RCC Lp4T0 0 =(collinR-seedTWth/2.0)*cos(4*colLHSan) =(collinR-seedTWth/2.0)*sin(4*colLHSan) \ 0 =(collTh+seedTWth)*cos(4*colLHSan) =(collTh+seedTWth)*sin(4*colLHSan) \ =collR+0 RCC Lp5T0 0 =(collinR-seedTWth/2.0)*cos(5*colLHSan) =(collinR-seedTWth/2.0)*sin(5*colLHSan) \ 0 =(collTh+seedTWth)*cos(5*colLHSan) =(collTh+seedTWth)*sin(5*colLHSan) \ =collR+0 RCC Lp6T0 0 =(collinR-seedTWth/2.0)*cos(6*colLHSan) =(collinR-seedTWth/2.0)*sin(6*colLHSan) \ 0 =(collTh+seedTWth)*cos(6*colLHSan) =(collTh+seedTWth)*sin(6*colLHSan) \ =collR+0 RCC Lp7T0 0 =(collinR-seedTWth/2.0)*cos(7*colLHSan) =(collinR-seedTWth/2.0)*sin(7*colLHSan) \ 0 =(collTh+seedTWth)*cos(7*colLHSan) =(collTh+seedTWth)*sin(7*colLHSan) \ =collR+0 RCC Lp8T0 0 =(collinR-seedTWth/2.0)*cos(8*colLHSan) =(collinR-seedTWth/2.0)*sin(8*colLHSan) \ 0 =(collTh+seedTWth)*cos(8*colLHSan) =(collTh+seedTWth)*sin(8*colLHSan) \ =collR+0 RCC Lp9T0 0 =(collinR-seedTWth/2.0)*cos(9*colLHSan) =(collinR-seedTWth/2.0)*sin(9*colLHSan) \ 0 =(collTh+seedTWth)*cos(9*colLHSan) =(collTh+seedTWth)*sin(9*colLHSan) \ =collR+0 $start_transform Tm9Rot PLA LLplnTm9 =cos(1*colTHSan/2.) =sin(1*colTHSan/2.) 0 0 0 0 $end_transform $start_transform Tm8Rot PLA LLplnTm8 =cos(1*colTHSan/2.) =sin(1*colTHSan/2.) 0 0 0 0 $end_transform $start_transform Tm7Rot PLA LLplnTm7 =cos(1*colTHSan/2.) =sin(1*colTHSan/2.) 0 0 0 0 $end_transform $start_transform Tm6Rot PLA LLplnTm6 =cos(1*colTHSan/2.) =sin(1*colTHSan/2.) 0 0 0 0 $end_transform $start_transform Tm5Rot PLA LLplnTm5 =cos(1*colTHSan/2.) =sin(1*colTHSan/2.) 0 0 0 0 $end_transform $start_transform Tm4Rot PLA LLplnTm4 =cos(1*colTHSan/2.) =sin(1*colTHSan/2.) 0 0 0 0 $end_transform $start_transform Tm3Rot PLA LLplnTm3 =cos(1*colTHSan/2.) =sin(1*colTHSan/2.) 0 0 0 0 $end_transform $start_transform Tm2Rot PLA LLplnTm2 =cos(1*colTHSan/2.) =sin(1*colTHSan/2.) 0 0 0 0 $end_transform $start_transform Tm1Rot PLA LLplnTm1 =cos(1*colTHSan/2.) =sin(1*colTHSan/2.) 0 0 0 0 $end_transform *Postive Channels $start_transform Tp1Rot PLA LRplnTp1 =cos(-1*colTHSan/2.) =sin(-1*colTHSan/2.) 0 0 0 0 $end_transform $start_transform Tp2Rot PLA LRplnTp2 =cos(-1*colTHSan/2.) =sin(-1*colTHSan/2.) 0 0 0 0 $end_transform $start_transform Tp3Rot PLA LRplnTp3 =cos(-1*colTHSan/2.) =sin(-1*colTHSan/2.) 0 0 0 0 $end_transform $start_transform Tp4Rot PLA LRplnTp4 =cos(-1*colTHSan/2.) =sin(-1*colTHSan/2.) 0 0 0 0 $end_transform $start_transform Tp5Rot PLA LRplnTp5 =cos(-1*colTHSan/2.) =sin(-1*colTHSan/2.) 0 0 0 0 $end_transform $start_transform Tp6Rot PLA LRplnTp6 =cos(-1*colTHSan/2.) =sin(-1*colTHSan/2.) 0 0 0 0 $end_transform $start_transform Tp7Rot PLA LRplnTp7 =cos(-1*colTHSan/2.) =sin(-1*colTHSan/2.) 0 0 0 0 $end_transform $start_transform Tp8Rot PLA LRplnTp8 =cos(-1*colTHSan/2.) =sin(-1*colTHSan/2.) 0 0 0 0 $end_transform $start_transform Tp9Rot PLA LRplnTp9 =cos(-1*colTHSan/2.) =sin(-1*colTHSan/2.) 0 0 0 0 $end_transform END *Black hole REGION BLKBODY 5 +blkbody -void *Void around REGION VOID 5 +void -airbod *Target REGION TARGBDY 5 +targbody *Void around REGION EXTAIR 5 +airbod -collbody -colTbOut -colCtPln | +airbod +colCtPln -collPlte alpha=100 *Void around REGION CAVAIR 5 +colTbIn -targbody -colCtPln | +collinSf -targbody -colCtPln alpha=100 *Void around REGION COLLAIR 5 (+ChnlBot +Lm9T0 | +ChnlBot +Lm8T0 | +ChnlBot +Lm7T0 | +ChnlBot +Lm6T0 | +ChnlBot +Lm5T0 | +ChnlBot +Lm4T0 | +ChnlBot +Lm3T0 | +ChnlBot +Lm2T0 | +ChnlBot +Lm1T0 | +ChnlBot +L0T0 | +ChnlBot +Lp1T0 | +ChnlBot +Lp2T0 | +ChnlBot +Lp3T0 | +ChnlBot +Lp4T0 | +ChnlBot +Lp5T0 | +ChnlBot +Lp6T0 | +ChnlBot +Lp7T0 | +ChnlBot +Lp8T0 | +ChnlBot +Lp9T0) -collinSf alpha=100 *Collimator REGION COLLTm9 +LLplnTm8 +collbody -LLplnTm9 -collinSf -ChEndL -ChEndR *Collimator REGION COLLTm8 +LLplnTm7 +collbody -LLplnTm8 -collinSf -ChEndL -ChEndR *Collimator REGION COLLTm7 +LLplnTm6 +collbody -LLplnTm7 -collinSf -ChEndL -ChEndR *Collimator REGION COLLTm6 +LLplnTm5 +collbody -LLplnTm6 -collinSf -ChEndL -ChEndR *Collimator REGION COLLTm5 +LLplnTm4 +collbody -LLplnTm5 -collinSf -ChEndL -ChEndR *Collimator REGION COLLTm4 +LLplnTm3 +collbody -LLplnTm4 -collinSf -ChEndL -ChEndR *Collimator REGION COLLTm3 +LLplnTm2 +collbody -LLplnTm3 -collinSf -ChEndL -ChEndR *Collimator REGION COLLTm2 +LLplnTm1 +collbody -LLplnTm2 -collinSf -ChEndL -ChEndR *Collimator REGION COLLTm1 +LLpln +collbody -LLplnTm1 -collinSf -ChEndL -ChEndR *Collimator REGION COLLPRTO +LRpln +collbody -LLpln -ChnlBot -ChEndL -ChEndR -( +ChnlTop +ChCt0 -colCtPln -ChnlBot -ChEndL -ChEndR) | +LRpln +ChnlBot -LLpln -collinSf -ChEndL -ChEndR -Lm9T0 -Lm8T0 -Lm7T0 -Lm6T0 -Lm5T0 -Lm4T0 -Lm3T0 -Lm2T0 -Lm1T0 -L0T0 -Lp1T0 -Lp2T0 -Lp3T0 -Lp4T0 -Lp5T0 -Lp6T0 -Lp7T0 -Lp8T0 -Lp9T0 *Collimator REGION COLLTp1 +LRplnTp1 +collbody -LRpln -collinSf -ChEndL -ChEndR *Collimator REGION COLLTp2 +LRplnTp2 +collbody -LRplnTp1 -collinSf -ChEndL -ChEndR *Collimator REGION COLLTp3 +LRplnTp3 +collbody -LRplnTp2 -collinSf -ChEndL -ChEndR *Collimator REGION COLLTp4 +LRplnTp4 +collbody -LRplnTp3 -collinSf -ChEndL -ChEndR *Collimator REGION COLLTp5 +LRplnTp5 +collbody -LRplnTp4 -collinSf -ChEndL -ChEndR *Collimator REGION COLLTp6 +LRplnTp6 +collbody -LRplnTp5 -collinSf -ChEndL -ChEndR *Collimator REGION COLLTp7 +LRplnTp7 +collbody -LRplnTp6 -collinSf -ChEndL -ChEndR *Collimator REGION COLLTp8 +LRplnTp8 +collbody -LRplnTp7 -collinSf -ChEndL -ChEndR *Collimator REGION COLLTp9 +LRplnTp9 +collbody -LRplnTp8 -collinSf -ChEndL -ChEndR *Collimator REGION COLLBODY +collbody +LLplnTm9 -collinSf -colTbIn -colCtPln -ChEndL -ChEndR | +collbody -collinSf -colTbIn -colCtPln -ChEndL -ChEndR -LRplnTp9 | +collPlte | +collbody +ChEndL -collinSf -colTbIn -colCtPln -ChEndR | +collbody +ChEndR -collinSf -colTbIn -colCtPln -ChEndL | +colTbOut +ChEndL -collbody -colTbIn -colCtPln -ChEndR | +colTbOut +ChEndR -collbody -colTbIn -colCtPln -ChEndL *Collimator Channel REGION COLLCHNL 5 ( +ChnlTop +ChCt0 -colCtPln -ChnlBot -ChEndL -ChEndR ) -(+seedendL| +seedbody -seedLpla +seedRpla | +seedendR ) *Seed Shell REGION SEEDSHEL 5 +seedendL +seedLpla | +seedbody -seedLpla +seedRpla -seedcav | +seedendR -seedRpla *Seed Cavity REGION SEEDCVTY 5 +seedcav -AgRdBdy *Silver Rod Body REGION AgROD 5 +AgRdBdy LATTICE -Tm9Rot COLLTm9,,, COLLLTm9 LATTICE -Tm8Rot COLLTm8,,, COLLLTm8 LATTICE -Tm7Rot COLLTm7,,, COLLLTm7 LATTICE -Tm6Rot COLLTm6,,, COLLLTm6 LATTICE -Tm5Rot COLLTm5,,, COLLLTm5 LATTICE -Tm4Rot COLLTm4,,, COLLLTm4 LATTICE -Tm3Rot COLLTm3,,, COLLLTm3 LATTICE -Tm2Rot COLLTm2,,, COLLLTm2 LATTICE -Tm1Rot COLLTm1,,, COLLLTm1 LATTICE -Tp1Rot COLLTp1,,, COLLLTp1 LATTICE -Tp2Rot COLLTp2,,, COLLLTp2 LATTICE -Tp3Rot COLLTp3,,, COLLLTp3 LATTICE -Tp4Rot COLLTp4,,, COLLLTp4 LATTICE -Tp5Rot COLLTp5,,, COLLLTp5 LATTICE -Tp6Rot COLLTp6,,, COLLLTp6 LATTICE -Tp7Rot COLLTp7,,, COLLLTp7 LATTICE -Tp8Rot COLLTp8,,, COLLLTp8 LATTICE -Tp9Rot COLLTp9,,, COLLLTp9 END GEOEND MATERIAL ZINC 30,, 7.133 *Brass (typical) *Brass can have a wide variation of Cu/Zn content. Pb is not *always present. MATERIAL Brass,,, 8.52 COMPOUND Brass 62.0 COPPER 35.0 ZINC 3.0 LEAD *..+....1....+....2....+....3....+....4....+....5....+....6....+....7.. ASSIGNMA , BLCKHOLE BLKBODY ASSIGNMA , VACUUM VOID ASSIGNMA , AIR EXTAIR COLLAIR ASSIGNMA , WATER TARGBDY ASSIGNMA , Brass COLLTm9 COLLBODY ASSIGNMA , TITANIUM SEEDSHEL ASSIGNMA , AIR SEEDCVTY ASSIGNMA , SILVER AgROD ASSIGNMA , POLYETHY COLLCHNL *For high threshold (> ~100 keV), secondary electrons are not generated and their contribution to the MCS should be taken in to account with the FUDGEM=1. *For low threshold, (<~5 keV), secondary electrons are explicitly produced and therefore FUDGEM=0 in order not to count them twice in the MCS. *In between, a linear interpolation from 0 to 1 can be used. However, the latter is intended only as rough approximation. EMFCUT PROD-CUT =-10.0*keV =10.0*keV =(10.0-5.0)/(100.0-5.0) VACUUM @LASTMAT *For high threshold (> ~100 keV), secondary electrons are not generated and their contribution to the MCS should be taken in to account with the FUDGEM=1. *For low threshold, (<~5 keV), secondary electrons are explicitly produced and therefore FUDGEM=0 in order not to count them twice in the MCS. *In between, a linear interpolation from 0 to 1 can be used. However, the latter is intended only as rough approximation. !EMFCUT PROD-CUT =-100*keV =100*keV =(100.0-5.0)/(100.0-5.0) Brass RADDECAY , 2 DCYSCORE USRBIN -1,,, SRCCHK doseqbot #if ScoreBDX DCYSCORE USRBDX -1,,, HOLEPRTS USRBDX HOLEPRTS 101 ALL-PART -26 COLLAIR CAVAIR,,,, 1,,, 1 #endif USRBIN SRCCHK 11 ALL-PART -21 =seedR-seedWTh 0.0 =seedCL/2 0.0 0.0 =-seedCL/2 100.0 \ 100. 100. ROTPRBIN ,, -seedPos,, SRCCHK USRBIN FlncXY 10 ALL-PART -22 10. 10. 0.1 -10. -10. -0.1 =20.0/(0.5*mm) =20.0/(0.5*mm) \ 1. USRBIN FlncXZ 10 ALL-PART -22 10. 0.1 10. -10. -0.1 -10. =20.0/(0.5*mm) 1. =20.0/(0.5*mm) USRBIN FlncYZ 10 ALL-PART -22 0.1 10. 10. -0.1 -10. -10. 1. =20.0/(0.5*mm) =20.0/(0.5*mm) USRBIN SdFlncYZ 10 ALL-PART -22 0.1 10. 10. -0.1 -10. -10. 1. =20.0/(0.5*mm) =20.0/(0.5*mm) ROTPRBIN ,, -seedTPln,, SdFlncYZ USRBIN dose 10 DOSE -23 =TargW/2.0 =TargH/2.0 =TargL/2.0/2.0 =-TargW/2.0 =-TargH/2.0 \ =-TargL/2.0/2.0 =TargW/(0.5*mm) =TargH/(0.5*mm) =TargL/(0.5*mm) USRBIN doseqtop 10 DOSE-EQ -24 =.5 =collinR+collTh+collShTh+1.0*cm .5 -.5 =collinR+collTh+collShTh \ -.5 1 1 1 USRBIN doseqbot 10 DOSE-EQ -24 =.5 =-(TargMovH)/2.0-collShTh .5 -.5 =-(TargMovH)/2.0-collShTh-1 \ -.5 1 1 1 *Set the random number seed RANDOMIZ , 1.0 *Set the number of primary histories to be simulated in the run START , 10000. STOP EndInput Notes: |If decays are simulated in semi-analogue mode, detector results are | expressed per unit primary weight (possibly scoring together prompt | and decay particles, if requested by DCYSCORE with WHAT(1) = -1.0). | If decays are instead calculated for requested cooling times, the | results are expressed per unit time (pSv/s, cm-2/s, Bq, etc.) | |When the source is a radioactive isotope (defined by command BEAM with | SDUM = ISOTOPE and by command HI-PROPErt), RADDECAY must be used in | semi-analogue mode (WHAT(1) > 1). The detector results are then | expressed per isotope decay. Note that command DCYSCORE must be issued | with WHAT(1) = -1, and must be applying to all relevant estimators | and detectors. Without DCYSCORE, no scoring will occur (see Note 8 to | command BEAM). | |when the DCYSCORE option is applied to a detector with | WHAT(1) = -1., {b:all quantities are normalized per unit primary | weight, or per decay if the source has been defined as a radioactive | isotope.} | |Energy deposition will be expressed in GeV per cm3 per unit primary | weight. Doses will be expressed in GeV/g per unit primary weight. To | obtain dose in Gy, multiply GeV/g by 1.602176462E-7. | |Bq = decay/s | |Gray/Primary = Gray/disintegration -> Gy/decay*decay/s =Gy/s | |< 2 mrem/hr = <20 uSv/hr for radiation safety | |Dose equivalent will be expressed in pSv per unit primary weight. | |pSv/decay = 10^-6 uSv/decay*decay/s*3600 s/hr = uSv/hr | |1 decay/s = 1 Bq = 2.703×10−11 Ci = 2.703x10-8 mCi => 1 MBq = 0.027 mCi | |Source Strength to Activity Conversion |1 U = 1 mGy m^2/h = 1 cGy cm^2/h = 0.348 mCi = 1.287x10+7 Bq for 137Cs | = 0.243 mCi = 8.99x10+6 Bq for 192Ir | = 0.486 mCi = 1.798x10+7 Bq for 198Au | = 0.787 mCi = 2.912x10+7 Bq for 125I | = 0.773 mCi = 2.860x10+7 Bq for 103Pd Page: Input Calculator: |1.2e-8*642e12*.5*amugr |> ans4 = 6.396400850399999e-18 | |642e12*0.5*amugr |> ans1 = 5.330334042e-10 |amugr |> ans2 = 1.6605402e-24 # Run information Run: Last: 1 End Run: Irr_I125_L0T0_tracks Title: Seed Longitudinal 0 position and Transverse 0 Position without BDX Scoring Define: seedWTh Define: seedR Define: seedL Define: seedCL Define: seedTWth Define: AgRodL Define: AgRodR Define: srcTh Define: TargW Define: TargH Define: TargL Define: TableW Define: TableL Define: TableTh Define: TargMovW Define: TargMovL Define: TargMovH Define: collTh Define: collR Define: collHgap Define: collinR Define: colloutR Define: collShTh Define: collOCl Define: collLBhA Define: collOAng Define: colHLAin Define: collLBEO Define: colLHSan Define: LngHnum Define: colTHSan Define: colTHFan Define: TrnsHnum Define: seedYshf Define: SeedU Define: t12I125 Define: mI125 Define: SaI125 Define: I125UtoA Define: I125A Define: nL=0 Define: nT=0 Define: I125 Rnd: 1015485966 Start: 100000 Last: 1 mode: fluka StartRun: 1591125372 End Run: Irr_I125_Lp5Tm5_tracks Title: Seed Longitudinal +5 position and Transverse -5 Position without BDX Scoring Define: seedWTh Define: seedR Define: seedL Define: seedCL Define: seedTWth Define: AgRodL Define: AgRodR Define: srcTh Define: TargW Define: TargH Define: TargL Define: TableW Define: TableL Define: TableTh Define: TargMovW Define: TargMovL Define: TargMovH Define: collTh Define: collR Define: collHgap Define: collinR Define: colloutR Define: collShTh Define: collOCl Define: collLBhA Define: collOAng Define: colHLAin Define: collLBEO Define: colLHSan Define: LngHnum Define: colTHSan Define: colTHFan Define: TrnsHnum Define: seedYshf Define: SeedU Define: t12I125 Define: mI125 Define: SaI125 Define: I125UtoA Define: I125A Define: nL=+5 Define: nT=-5 Define: I125 Rnd: 362019020 Start: 100000 Last: 1 mode: fluka StartRun: 1591125372 End Run: Irr_I125_L0T0_tracks-BDX Title: Seed Longitudinal 0 position and Transverse 0 Position with BDX Scoring Define: seedWTh Define: seedR Define: seedL Define: seedCL Define: seedTWth Define: AgRodL Define: AgRodR Define: srcTh Define: TargW Define: TargH Define: TargL Define: TableW Define: TableL Define: TableTh Define: TargMovW Define: TargMovL Define: TargMovH Define: collTh Define: collR Define: collHgap Define: collinR Define: colloutR Define: collShTh Define: collOCl Define: collLBhA Define: collOAng Define: colHLAin Define: collLBEO Define: colLHSan Define: LngHnum Define: colTHSan Define: colTHFan Define: TrnsHnum Define: seedYshf Define: SeedU Define: t12I125 Define: mI125 Define: SaI125 Define: I125UtoA Define: I125A Define: nL=0 Define: nT=0 Define: I125 Define: ScoreBDX Rnd: 340118451 Start: 100000 Last: 1 mode: fluka StartRun: 1591125372 End Run: Irr_I125_Lp5Tm5_tracks-BDX Title: Seed Longitudinal +5 position and Transverse -5 Position with BDX Scoring Define: seedWTh Define: seedR Define: seedL Define: seedCL Define: seedTWth Define: AgRodL Define: AgRodR Define: srcTh Define: TargW Define: TargH Define: TargL Define: TableW Define: TableL Define: TableTh Define: TargMovW Define: TargMovL Define: TargMovH Define: collTh Define: collR Define: collHgap Define: collinR Define: colloutR Define: collShTh Define: collOCl Define: collLBhA Define: collOAng Define: colHLAin Define: collLBEO Define: colLHSan Define: LngHnum Define: colTHSan Define: colTHFan Define: TrnsHnum Define: seedYshf Define: SeedU Define: t12I125 Define: mI125 Define: SaI125 Define: I125UtoA Define: I125A Define: nL=+5 Define: nT=-5 Define: I125 Define: ScoreBDX Rnd: 703475129 Start: 100000 Last: 1 mode: fluka StartRun: 1591125372 End # Custom materials information Geometry: Geometry Lock: 1 Layer: Media Options.crosshair: 17 Options.latticelevel: 32 Show.lattice: 1 End Layer: 3D 3D: 1 3D._visible: 1 3D.clip.0: clipbox 3D.edgedetect: 1 3D.xray: 172 Options.latticelevel: 0 Options.viewport: 0 Show: 0 Show.lattice: 1 Show.voxel: 1 End Layer: 3D Tracks 3D: 1 3D._visible: 1 3D.clip.0: clipbox 3D.edgedetect: 1 3D.skip1stblack: 1 3D.xray: 172 Options.viewport: 0 Palette._visible: 1 Palette.alphamin: 1 Show: 0 Show.lattice: 1 Show.voxel: 1 Userdump: 1 Userdump.ELECTRON: 1 Userdump.PHOTON: 1 Userdump.POSITRON: 1 Userdump._visible: 1 Userdump.color.ELECTRON: #ff0000 Userdump.color.PHOTON: #00ff00 Userdump.color.POSITRON: #00ffff Userdump.color.PROTON: #00ff00 Userdump.file: Irr_I125_L0T0_tracks001_draw Userdump.n: 1000 Usrbin: 1 Usrbin._visible: 1 Usrbin.det: 4 Usrbin.file: Irr_I125_L0T0_tracks_22.bnn Usrbin.rotdefi: -seedTPln End Layer: 3D Tracks Off 3D: 1 3D._visible: 1 3D.ambient: 68 3D.clip.0: clipbox 3D.edgedetect: 1 3D.shadows: 0 3D.skip1stblack: 1 3D.xray: 65 Options.latticelevel: 0 Options.viewport: 0 Palette: 1 Palette._visible: 1 Palette.alphamin: 1 Show: 0 Show.lattice: 1 Show.voxel: 1 Userdump: 1 Userdump.ELECTRON: 1 Userdump.PHOTON: 1 Userdump.POSITRON: 1 Userdump._visible: 1 Userdump.color.ELECTRON: #ff0000 Userdump.color.PHOTON: #00ff00 Userdump.color.POSITRON: #00ffff Userdump.color.PROTON: #00ff00 Userdump.file: Irr_I125_Lp5Tm5_tracks001_draw Userdump.n: 1000 Usrbin: 1 Usrbin._visible: 1 Usrbin.det: 4 Usrbin.file: Irr_I125_Lp5Tm5_tracks_22.bnn Usrbin.rotdefi: -seedTPln End End # Plots information Plot: Red Type: Geometry basisu: 0 0 1 coord: Z-Y extends: 13.9161 6.92104 grid: 0 kind: Material layer: Media origin: 0 3.6563 0.154359 ratio: 0.459457975577811 run: Irr_I125_Lp5Tm5_tracks End Plot: Green Type: Geometry coord: X-Y extends: 13.4588 6.70252 grid: 0 kind: Material layer: Media origin: -0.110338 3.86048 0 End Plot: Blue Type: Geometry basisu: 0.841116 0.14322 0.521547 basisv: 0.185069 0.829878 -0.526357 extends: 16.6993 8.28303 grid: 0 kind: Material layer: 3D Tracks origin: -9.48508 10.9583 12.2001 ratio: 1 End Plot: Magenta Type: Geometry basisu: 0 0 1 basisv: 1 0 0 coord: Z-X extends: 20.0125 9.93963 grid: 0 kind: Material layer: Media origin: 1.24369 0.2 -3.59935 ratio: 0.4582299811870859 End Plot: IrrHalf_I125_Seed_All_FluenceXY Type: USRBIN cbtics: 1 datafile: Irr_I125_L0T0_tracks_21.bnn geo: -No- grid: 0 hist: 2D Projection int: 2.9469667919642312E-002 max: 84.4406586 min: 9.88895670E-02 ratio: 1.0 title: Seed Fluence (All) xlabel: X (cm) ylabel: Y (cm) End Plot: IrrHalf_I125_L0T0_All_FluenceXY Type: USRBIN cblog: 1 cbtics: 1 datafile: Irr_I125_L0T0_tracks_22.bnn grid: 0 hist: 2D Projection index: 1 int: 6.4539305103910483E-002 max: 16.8145924 min: 8.18006592E-05 ratio: 3.25 title: All Particle Fluence xlabel: Z (cm) xmax: 2 xmin: -2 ylabel: Y (cm) ymax: 10 ymin: -3 End Plot: IrrHalf_I125_L0T0_All_FluenceYZ Type: USRBIN cblog: 1 cbtics: 1 datafile: Irr_I125_L0T0_tracks_22.bnn det: 3 grid: 0 hist: 2D Projection index: 2 int: 0.10367169859923403 max: 10.1736927 min: 5.41566260E-06 proj: X ratio: 0.5416666666666666 title: All Particle Fluence xlabel: Z (cm) xmax: 12 xmin: -12 ylabel: Y (cm) ymax: 10 ymin: -3 End Plot: IrrHalf_I125_L0T0_All_SeedPlnFluenceYZ Type: USRBIN cblog: 1 cbtics: 1 datafile: Irr_I125_L0T0_tracks_22.bnn det: 4 grid: 0 hist: 2D Projection index: 2 int: 0.10321598791833052 max: 10.9317017 min: 1.97873916E-04 proj: X ratio: 0.5416666666666666 title: All Particle Fluence xlabel: Z (cm) xmax: 12 xmin: -12 ylabel: Y (cm) ymax: 10 ymin: -3 End Plot: IrrHalf_I125_L0T0_All_SeedPlnFluenceYZ01 Type: USRBIN cblog: 1 cbtics: 1 datafile: Irr_I125_Lp5Tm5_tracks_22.bnn det: 4 grid: 0 hist: 2D Projection index: 2 int: 0.10341511789064944 max: 10.9363537 min: 1.07198612E-04 proj: X ratio: 0.5416666666666666 title: All Particle Fluence xlabel: Z (cm) xmax: 12 xmin: -12 ylabel: Y (cm) ymax: 10 ymin: -3 End Plot: IrrHalf_I125_L0T0_DoseYZ Type: USRBIN bxmax: 0.05 bxmin: -0.05 cblabel: Dose Rate (cGy/min) cbtics: 1 commands: set format cb "%.2e" datafile: Irr_I125_L0T0_tracks_23.bnn grid: 0 hist: 2D Projection index: 3 int: 0.0000000000000000 max: 0.00000000 min: 1.00000002E+30 norm: 1.602176462E-7*0.5*2.912E7*60*100 proj: X ratio: 1 title: Longitudinal Dose L0T0 xlabel: Y (cm) ylabel: Z (cm) End Plot: IrrHalf_I125_L0T0_DoseXZ Type: USRBIN bymax: 1 bymin: 0.9 cblabel: Dose Rate (cGy/min) cbtics: 1 commands: set format cb "%.2e" datafile: Irr_I125_L0T0_tracks_23.bnn grid: 0 hist: 2D Projection index: 4 int: 8.6503850751595431E-011 max: 5.32893552E-09 min: 7.94107754E-14 norm: 1.602176462E-7*0.5*2.912E7*60*100 proj: Y ratio: 1 title: Top Dose L0T0 xlabel: X (cm) ylabel: Z (cm) End Plot: IrrHalf_I125_Lp2T0_DoseYZ Type: USRBIN bxmax: 0.1 bxmin: -0.05 cblabel: Dose Rate (cGy/min) cblog: 1 cbtics: 1 commands: set format cb "%.2e" datafile: Irr_I125_Lp2T0_22.bnn det: 3 grid: 0 hist: 2D Projection index: 3 int: 9.1665068862364571E-002 max: 16.5854530 min: 1.05816478E-04 norm: 1.602176462E-7*0.5*2.912E7*60*100 proj: X ratio: 1 title: Longitudinal Dose Lp2T0 xlabel: Y (cm) ylabel: Z (cm) End Plot: IrrHalf_I125_Lm2T0_DoseYZ Type: USRBIN bxmax: 0.1 bxmin: -0.05 cblabel: Dose Rate (cGy/min) cblog: 1 cbtics: 1 commands: set format cb "%.2e" datafile: Irr_I125_Lm2T0_22.bnn det: 3 grid: 0 hist: 2D Projection index: 3 int: 9.4402698980653788E-002 max: 16.5846119 min: 2.79698688E-05 norm: 1.602176462E-7*0.5*2.912E7*60*100 proj: X ratio: 1 title: Longitudinal Dose Lm2T0 xlabel: Y (cm) ylabel: Z (cm) End Plot: IrrHalf_I125_Combo_All_FluenceYZ Type: USRBIN bxmax: 0.1 bxmin: -0.05 cblog: 1 cbtics: 1 commands: set format cb "%.2e" datafile: Irr_I125_Combo_22.bnn det: 3 grid: 0 hist: 2D Projection index: 2 int: 6.8961578629030892E-002 max: 1.71767342 min: 1.27174766E-22 proj: X ratio: 0.5416666666666666 title: All Particle Fluence xlabel: Z (cm) xmax: 12 xmin: -12 ylabel: Y (cm) ymax: 10 ymin: -3 End Plot: IrrHalf_I125_Combo_DoseYZ Type: USRBIN bxmax: 0.05 bxmin: -0.05 cblabel: Dose Rate (cGy/min) cbtics: 1 commands: set format cb "%.2e" datafile: Irr_I125_Combo_23.bnn grid: 0 hist: 2D Projection index: 3 int: 2.8509035331618879E-010 max: 3.02359537E-09 min: 5.03379568E-15 norm: 1.602176462E-7*0.5*2.912E7*60*100 proj: X ratio: 0.4 title: Longitudinal Dose xlabel: Y (cm) ylabel: Z (cm) End