Difference: DPMatterStudies (1 vs. 7)

Revision 72018-09-14 - AlessioSarti

META TOPICPARENT	name="DoseProfiler"

Dose Profiler matter studies page

Dose Profiler matter studies page
Rando Full Simulation studies
The weights definition + machinery
How to run simulations and produce files to apply weights:
The mental ball configurations
- Comments

This page documents the matter studies.

Rando Full Simulation studies

The file chain is the following:

Simulation input: /NAS_arpg/mfischetti/SimRandoLight/notr_profDATA_rando_280_15R.inp

Simulation size: 10^7 12C, 280 MeV, sigma beam as from CNAO.

Detector position: 51 cm from isocentre. Flag for running AnaProfi, check link in CNAO TB page. -ang 58 -dang -6 -phi 1.4

Rando position.

Ana Profi

Added:

>
>

To compute weights, on the selected tracks, run FRED.

/home_arpg/mfischetti/MeMs/MLEM/plugin

To run it

modify runMode = 4
put the correct detector positioning.
make -C plugin&&freddevel
- Run on the file produced by AnaProfi. Should compare with. /NAS_arpg/MeMs/RandoData/rando280_tracks.txt
- Produces as output something: ProfilerWeights2D.txt (check plugin for content).

In order to convert the Z with the rhoZ a map can be used.

/home_arpg/mfischetti/MeMs/MLEM/localout/BP.txt

The weights definition + machinery

Here we document how the weights are defined, where is the code

How to run simulations and produce files to apply weights:

To Run V2:

-reference cilinder, necessary if the target does not start at 0cm (as reference cilinder3 does) -
$FLUKA/flutil/rfluka -e fluka_FILTER_mbcil.exe -N0 -M5 filter_mbcil_1R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbcil.exe -N0 -M5 filter_mbcil_2R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbcil.exe -N0 -M5 filter_mbcil_3R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbcil.exe -N0 -M5 filter_mbcil_4R &

-mb brain+bone radius_ext = 10 cm centered in (0,0,10)
$FLUKA/flutil/rfluka -e fluka_FILTER_mbfull.exe -N0 -M5 filter_mbfull_1R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbfull.exe -N0 -M5 filter_mbfull_2R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbfull.exe -N0 -M5 filter_mbfull_3R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbfull.exe -N0 -M5 filter_mbfull_4R &

cat filter_mbcil*00*dat > filter_mbcil_all_TXT.dat &
cat filter_mbfull*00*dat > filter_mbfull_all_TXT.dat &

./Txt2Root -in filter_mbcil_all_TXT.dat -out filter_mbcil_all.root &

Then in ../ApplyWeights:

see ListApplyWeights _mbfull.sh (the first part) ->to have reference hZ distr

Then:

cd ../V5

ln -s ../V2/filter_mbfull_all_TXT.dat secondary_TXT.dat

and run FLUKA on new ray*inp files (same full geo, different beam (10TeV muon+) and physics card)
NB: check the source for reso!
$FLUKA/flutil/rfluka -e fluka_RAY.exe -N0 -M1 ray_filter_mbfull_reso0 &
$FLUKA/flutil/rfluka -e fluka_RAY.exe -N0 -M1 ray_filter_mbfull_reso1 &
$FLUKA/flutil/rfluka -e fluka_RAY.exe -N0 -M1 ray_filter_mbfull_reso2 &
$FLUKA/flutil/rfluka -e fluka_RAY.exe -N0 -M1 ray_filter_mbfull_reso3 &

./Txt2Root -in ray_filter_mbfull_reso0001_TXT.dat -out ray_filter_mbfull_reso0.root &
./Txt2Root -in ray_filter_mbfull_reso1001_TXT.dat -out ray_filter_mbfull_reso1.root &
./Txt2Root -in ray_filter_mbfull_reso2001_TXT.dat -out ray_filter_mbfull_reso2.root &
./Txt2Root -in ray_filter_mbfull_reso3001_TXT.dat -out ray_filter_mbfull_reso3.root &

Then in ApplyWeights:

see ListApplyWeights _mbfull.sh (second part)

to produce the weighted Z distr

- Compare the distributions with TestWeighs _mb.cc o quello che hai <3

The mental ball configurations

Beam: 12C ions, energy: 300 MeV /u, spatial distribution: gaussian with FWHM_x,y = 0.5 cm

direction: positive z axis, starting point:(0,0,4) cm

1)*mbfull, REFERENCE distr: mbcil*

Mental Ball: sphere with rmin = 9.6 cm, rmax = 10 cm, inner: BRAIN (ICRU), outer: BONE (ICRU)

center:(0,0,15) cm

- fluka region number of detectors (to be implemented in ApplyWeightsLib.h): {3,4,5,6,7,8}

REFERENCE distr: mbcil = MB with same geo, material: water, cut with cilinder of r = 1.5 cm

2)*mbimp, REFERENCE distr: mbimp_cil*

Mental Ball: sphere with rmin = 9.6 cm, rmax = 10 cm, inner: BRAIN (ICRU), outer: BONE (ICRU)

center:(0,0,15) cm

Implant: sphere, r = 1cm, material: AIR, center: (0.,5.,16.)

- fluka region number of detectors (to be implemented in ApplyWeightsLib.h): {0,0,0,3,4,5}

NB: the 60deg detectors are cut in order to have just a portion of ~ 20 x 20 cm2 ; no 90deg detectors in this configuration

REFERENCE distr: mbimp_cil = MB with same geo, material: water, cut with cilinder of r = 1.5 cm (as in sim 1) ); detectors in same config of mbimp sim

3)*mbimp2, REFERENCE distr: mbimp_cil*

Mental Ball: sphere with rmin = 9.6 cm, rmax = 10 cm, inner: BRAIN (ICRU), outer: BONE (ICRU)

center:(0,0,15) cm

Implant: sphere, r = 2cm, material: AIR, center: (0.,5.,16.)

- fluka region number of detectors (to be implemented in ApplyWeightsLib.h): {0,0,0,3,4,5}

NB: the 60deg detectors are cut in order to have just a portion of ~ 20 x 20 cm2 ; no 90deg detectors in this configuration

REFERENCE distr: mbimp_cil = MB with same geo, material: water, cut with cilinder of r = 1.5 cm (as in sim 1) ); detectors in same config of mbimp sim

4)*mbimp3, REFERENCE distr: mbimp_cil*

Mental Ball: sphere with rmin = 9.6 cm, rmax = 10 cm, inner: BRAIN (ICRU), outer: BONE (ICRU)

center:(0,0,15) cm

Implant: sphere, r = 1cm, material: ADIPOSE TISSUE, center: (0.,5.,16.)

- fluka region number of detectors (to be implemented in ApplyWeightsLib.h): {0,0,0,3,4,5}

NB: the 60deg detectors are cut in order to have just a portion of ~ 20 x 20 cm2 ; no 90deg detectors in this configuration

REFERENCE distr: mbimp_cil = MB with same geo, material: water, cut with cilinder of r = 1.5 cm (as in sim 1) ); detectors in same config of mbimp sim

Mental ball with inset on beam

Comments

<--/commentPlugin-->

META FILEATTACHMENT	attachment="MBimpB_5cm_YZ.jpeg" attr="" comment="" date="1519920220" name="MBimpB_5cm_YZ.jpeg" path="MBimpB_5cm_YZ.jpeg" size="24644" user="AlessioSarti" version="2"
META FILEATTACHMENT	attachment="MBcil_0cm_zoom.png" attr="" comment="" date="1519920358" name="MBcil_0cm_zoom.png" path="MBcil_0cm_zoom.png" size="51648" user="AlessioSarti" version="1"
META FILEATTACHMENT	attachment="MBcil_0cm.png" attr="" comment="" date="1519920358" name="MBcil_0cm.png" path="MBcil_0cm.png" size="113631" user="AlessioSarti" version="1"
META FILEATTACHMENT	attachment="MBfull_0cm_zoom.png" attr="" comment="" date="1519920358" name="MBfull_0cm_zoom.png" path="MBfull_0cm_zoom.png" size="66281" user="AlessioSarti" version="1"
META FILEATTACHMENT	attachment="MBfull_0cm.png" attr="" comment="" date="1519920358" name="MBfull_0cm.png" path="MBfull_0cm.png" size="113672" user="AlessioSarti" version="1"

Revision 62018-09-14 - AlessioSarti

META TOPICPARENT	name="DoseProfiler"

Dose Profiler matter studies page

Dose Profiler matter studies page
Rando Full Simulation studies
The weights definition + machinery
How to run simulations and produce files to apply weights:
The mental ball configurations
- Comments

This page documents the matter studies.

Added:

>
>

Rando Full Simulation studies

The file chain is the following:

Simulation input: /NAS_arpg/mfischetti/SimRandoLight/notr_profDATA_rando_280_15R.inp

Simulation size: 10^7 12C, 280 MeV, sigma beam as from CNAO.

Detector position: 51 cm from isocentre. Flag for running AnaProfi, check link in CNAO TB page. -ang 58 -dang -6 -phi 1.4

Rando position.

Ana Profi

The weights definition + machinery

Here we document how the weights are defined, where is the code

How to run simulations and produce files to apply weights:

To Run V2:

cat filter_mbcil*00*dat > filter_mbcil_all_TXT.dat &
cat filter_mbfull*00*dat > filter_mbfull_all_TXT.dat &

./Txt2Root -in filter_mbcil_all_TXT.dat -out filter_mbcil_all.root &

Then in ../ApplyWeights:

see ListApplyWeights _mbfull.sh (the first part) ->to have reference hZ distr

Then:

cd ../V5

ln -s ../V2/filter_mbfull_all_TXT.dat secondary_TXT.dat

Then in ApplyWeights:

see ListApplyWeights _mbfull.sh (second part)

to produce the weighted Z distr

- Compare the distributions with TestWeighs _mb.cc o quello che hai <3

The mental ball configurations

Beam: 12C ions, energy: 300 MeV /u, spatial distribution: gaussian with FWHM_x,y = 0.5 cm

direction: positive z axis, starting point:(0,0,4) cm

1)*mbfull, REFERENCE distr: mbcil*

Mental Ball: sphere with rmin = 9.6 cm, rmax = 10 cm, inner: BRAIN (ICRU), outer: BONE (ICRU)

center:(0,0,15) cm

- fluka region number of detectors (to be implemented in ApplyWeightsLib.h): {3,4,5,6,7,8}

REFERENCE distr: mbcil = MB with same geo, material: water, cut with cilinder of r = 1.5 cm

2)*mbimp, REFERENCE distr: mbimp_cil*

Mental Ball: sphere with rmin = 9.6 cm, rmax = 10 cm, inner: BRAIN (ICRU), outer: BONE (ICRU)

center:(0,0,15) cm

Changed:

<
<

Implant: sphere, r = 1cm, material: AIR, center: (0.,5.,16.)

>
>

Implant: sphere, r = 1cm, material: AIR, center: (0.,5.,16.)

- fluka region number of detectors (to be implemented in ApplyWeightsLib.h): {0,0,0,3,4,5}

NB: the 60deg detectors are cut in order to have just a portion of ~ 20 x 20 cm2 ; no 90deg detectors in this configuration

REFERENCE distr: mbimp_cil = MB with same geo, material: water, cut with cilinder of r = 1.5 cm (as in sim 1) ); detectors in same config of mbimp sim

3)*mbimp2, REFERENCE distr: mbimp_cil*

Mental Ball: sphere with rmin = 9.6 cm, rmax = 10 cm, inner: BRAIN (ICRU), outer: BONE (ICRU)

center:(0,0,15) cm

Changed:

<
<

Implant: sphere, r = 2cm, material: AIR, center: (0.,5.,16.)

>
>

Implant: sphere, r = 2cm, material: AIR, center: (0.,5.,16.)

- fluka region number of detectors (to be implemented in ApplyWeightsLib.h): {0,0,0,3,4,5}

NB: the 60deg detectors are cut in order to have just a portion of ~ 20 x 20 cm2 ; no 90deg detectors in this configuration

REFERENCE distr: mbimp_cil = MB with same geo, material: water, cut with cilinder of r = 1.5 cm (as in sim 1) ); detectors in same config of mbimp sim

4)*mbimp3, REFERENCE distr: mbimp_cil*

Mental Ball: sphere with rmin = 9.6 cm, rmax = 10 cm, inner: BRAIN (ICRU), outer: BONE (ICRU)

center:(0,0,15) cm

Changed:

<
<

Implant: sphere, r = 1cm, material: ADIPOSE TISSUE, center: (0.,5.,16.)

>
>

Implant: sphere, r = 1cm, material: ADIPOSE TISSUE, center: (0.,5.,16.)

- fluka region number of detectors (to be implemented in ApplyWeightsLib.h): {0,0,0,3,4,5}

NB: the 60deg detectors are cut in order to have just a portion of ~ 20 x 20 cm2 ; no 90deg detectors in this configuration

REFERENCE distr: mbimp_cil = MB with same geo, material: water, cut with cilinder of r = 1.5 cm (as in sim 1) ); detectors in same config of mbimp sim

Mental ball with inset on beam

Changed:

<
<

>
>

Comments

<--/commentPlugin-->

META FILEATTACHMENT	attachment="MBimpB_5cm_YZ.jpeg" attr="" comment="" date="1519920220" name="MBimpB_5cm_YZ.jpeg" path="MBimpB_5cm_YZ.jpeg" size="24644" user="AlessioSarti" version="2"
META FILEATTACHMENT	attachment="MBcil_0cm_zoom.png" attr="" comment="" date="1519920358" name="MBcil_0cm_zoom.png" path="MBcil_0cm_zoom.png" size="51648" user="AlessioSarti" version="1"
META FILEATTACHMENT	attachment="MBcil_0cm.png" attr="" comment="" date="1519920358" name="MBcil_0cm.png" path="MBcil_0cm.png" size="113631" user="AlessioSarti" version="1"
META FILEATTACHMENT	attachment="MBfull_0cm_zoom.png" attr="" comment="" date="1519920358" name="MBfull_0cm_zoom.png" path="MBfull_0cm_zoom.png" size="66281" user="AlessioSarti" version="1"
META FILEATTACHMENT	attachment="MBfull_0cm.png" attr="" comment="" date="1519920358" name="MBfull_0cm.png" path="MBfull_0cm.png" size="113672" user="AlessioSarti" version="1"

Revision 52018-03-01 - AlessioSarti

 			 META TOPICPARENT 
			 name="DoseProfiler" 
		
 Dose Profiler matter studies page 


 
  Dose Profiler matter studies page
   Rando Full Simulation studies
   The weights definition + machinery
   How to run simulations and produce files to apply weights:
   The mental ball configurations 
  Comments
 
 


This page documents the matter studies.

 The weights definition + machinery 

Here we document how the weights are defined, where is the code

 How to run simulations and produce files to apply weights: 

To Run V2:

-reference cilinder, necessary if the target does not start at 0cm (as reference cilinder3 does) -
$FLUKA/flutil/rfluka -e fluka_FILTER_mbcil.exe -N0 -M5 filter_mbcil_1R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbcil.exe -N0 -M5 filter_mbcil_2R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbcil.exe -N0 -M5 filter_mbcil_3R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbcil.exe -N0 -M5 filter_mbcil_4R &

-mb brain+bone radius_ext = 10 cm centered in (0,0,10)
$FLUKA/flutil/rfluka -e fluka_FILTER_mbfull.exe -N0 -M5 filter_mbfull_1R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbfull.exe -N0 -M5 filter_mbfull_2R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbfull.exe -N0 -M5 filter_mbfull_3R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbfull.exe -N0 -M5 filter_mbfull_4R &

cat filter_mbcil*00*dat > filter_mbcil_all_TXT.dat &
cat filter_mbfull*00*dat > filter_mbfull_all_TXT.dat &

./Txt2Root -in filter_mbcil_all_TXT.dat -out filter_mbcil_all.root &

Then in ../ApplyWeights:
- META TOPICPARENT
+ name="DoseProfiler"
-<
<
+see ListApplyWeights _mbfull.sh (the first part) ->to have reference hZ distr
->
>
+see ListApplyWeights _mbfull.sh (the first part) ->to have reference hZ distr
 Then:

cd ../V5

ln -s ../V2/filter_mbfull_all_TXT.dat secondary_TXT.dat

and run FLUKA on new ray*inp files (same full geo, different beam (10TeV muon+) and physics card) 
NB: check the source for reso!
$FLUKA/flutil/rfluka -e fluka_RAY.exe -N0 -M1 ray_filter_mbfull_reso0 &
$FLUKA/flutil/rfluka -e fluka_RAY.exe -N0 -M1 ray_filter_mbfull_reso1 &
$FLUKA/flutil/rfluka -e fluka_RAY.exe -N0 -M1 ray_filter_mbfull_reso2 &
$FLUKA/flutil/rfluka -e fluka_RAY.exe -N0 -M1 ray_filter_mbfull_reso3 &

./Txt2Root -in ray_filter_mbfull_reso0001_TXT.dat -out ray_filter_mbfull_reso0.root &
./Txt2Root -in ray_filter_mbfull_reso1001_TXT.dat -out ray_filter_mbfull_reso1.root &
./Txt2Root -in ray_filter_mbfull_reso2001_TXT.dat -out ray_filter_mbfull_reso2.root &
./Txt2Root -in ray_filter_mbfull_reso3001_TXT.dat -out ray_filter_mbfull_reso3.root &

Then in ApplyWeights:

see ListApplyWeights _mbfull.sh (second part)

to produce the weighted Z distr

- Compare the distributions with TestWeighs _mb.cc o quello che hai <3

 The mental ball configurations 

Beam: 12C ions, energy: 300 MeV /u, spatial distribution: gaussian with FWHM_x,y = 0.5 cm

direction: positive z axis, starting point:(0,0,4) cm

1)*mbfull, REFERENCE distr: mbcil*

Mental Ball: sphere with rmin = 9.6 cm, rmax = 10 cm, inner: BRAIN (ICRU), outer: BONE (ICRU)

center:(0,0,15) cm
-<
<
+- fluka region number of detectors (to be implemented in ApplyWeightsLib.h): {3,4,5,6,7,8}
->
>
+- fluka region number of detectors (to be implemented in ApplyWeightsLib.h): {3,4,5,6,7,8}
 REFERENCE distr: mbcil = MB with same geo, material: water, cut with cilinder of r = 1.5 cm

2)*mbimp, REFERENCE distr: mbimp_cil*

Mental Ball: sphere with rmin = 9.6 cm, rmax = 10 cm, inner: BRAIN (ICRU), outer: BONE (ICRU)

center:(0,0,15) cm
-<
<
+Implant: sphere, r = 1cm, material: AIR, center: (0.,5.,16.)
->
>
+Implant: sphere, r = 1cm, material: AIR, center: (0.,5.,16.)
-<
<
+- fluka region number of detectors (to be implemented in ApplyWeightsLib.h): {0,0,0,3,4,5}
->
>
+- fluka region number of detectors (to be implemented in ApplyWeightsLib.h): {0,0,0,3,4,5}
-<
<
+ NB: the 60deg detectors are cut in order to have just a portion of ~ 20 x 20 cm2 ; no 90deg detectors in this configuration
->
>
+NB: the 60deg detectors are cut in order to have just a portion of ~ 20 x 20 cm2 ; no 90deg detectors in this configuration
 REFERENCE distr: mbimp_cil = MB with same geo, material: water, cut with cilinder of r = 1.5 cm (as in sim 1) ); detectors in same config of mbimp sim

3)*mbimp2, REFERENCE distr: mbimp_cil*

Mental Ball: sphere with rmin = 9.6 cm, rmax = 10 cm, inner: BRAIN (ICRU), outer: BONE (ICRU)

center:(0,0,15) cm
-<
<
+Implant: sphere, r = 2cm, material: AIR, center: (0.,5.,16.)
->
>
+Implant: sphere, r = 2cm, material: AIR, center: (0.,5.,16.)
-<
<
+- fluka region number of detectors (to be implemented in ApplyWeightsLib.h): {0,0,0,3,4,5}
->
>
+- fluka region number of detectors (to be implemented in ApplyWeightsLib.h): {0,0,0,3,4,5}
-<
<
+ NB: the 60deg detectors are cut in order to have just a portion of ~ 20 x 20 cm2 ; no 90deg detectors in this configuration
->
>
+NB: the 60deg detectors are cut in order to have just a portion of ~ 20 x 20 cm2 ; no 90deg detectors in this configuration
 REFERENCE distr: mbimp_cil = MB with same geo, material: water, cut with cilinder of r = 1.5 cm (as in sim 1) ); detectors in same config of mbimp sim

4)*mbimp3, REFERENCE distr: mbimp_cil*

Mental Ball: sphere with rmin = 9.6 cm, rmax = 10 cm, inner: BRAIN (ICRU), outer: BONE (ICRU)

center:(0,0,15) cm
-<
<
+Implant: sphere, r = 1cm, material: ADIPOSE TISSUE, center: (0.,5.,16.)
->
>
+Implant: sphere, r = 1cm, material: ADIPOSE TISSUE, center: (0.,5.,16.)
-<
<
+- fluka region number of detectors (to be implemented in ApplyWeightsLib.h): {0,0,0,3,4,5}
->
>
+- fluka region number of detectors (to be implemented in ApplyWeightsLib.h): {0,0,0,3,4,5}
-<
<
+ NB: the 60deg detectors are cut in order to have just a portion of ~ 20 x 20 cm2 ; no 90deg detectors in this configuration
->
>
+NB: the 60deg detectors are cut in order to have just a portion of ~ 20 x 20 cm2 ; no 90deg detectors in this configuration
 REFERENCE distr: mbimp_cil = MB with same geo, material: water, cut with cilinder of r = 1.5 cm (as in sim 1) ); detectors in same config of mbimp sim
->
>
+Mental ball with inset on beam
  Comments 


 
  (Edit - Preview)
 
<--/commentPlugin-->
->
>
+			 META FILEATTACHMENT 
			 attachment="MBimpB_5cm_YZ.jpeg" attr="" comment="" date="1519920220" name="MBimpB_5cm_YZ.jpeg" path="MBimpB_5cm_YZ.jpeg" size="24644" user="AlessioSarti" version="2" 
		
		
			 META FILEATTACHMENT 
			 attachment="MBcil_0cm_zoom.png" attr="" comment="" date="1519920358" name="MBcil_0cm_zoom.png" path="MBcil_0cm_zoom.png" size="51648" user="AlessioSarti" version="1" 
		
		
			 META FILEATTACHMENT 
			 attachment="MBcil_0cm.png" attr="" comment="" date="1519920358" name="MBcil_0cm.png" path="MBcil_0cm.png" size="113631" user="AlessioSarti" version="1" 
		
		
			 META FILEATTACHMENT 
			 attachment="MBfull_0cm_zoom.png" attr="" comment="" date="1519920358" name="MBfull_0cm_zoom.png" path="MBfull_0cm_zoom.png" size="66281" user="AlessioSarti" version="1" 
		
		
			 META FILEATTACHMENT 
			 attachment="MBfull_0cm.png" attr="" comment="" date="1519920358" name="MBfull_0cm.png" path="MBfull_0cm.png" size="113672" user="AlessioSarti" version="1"
- META FILEATTACHMENT
+ attachment="MBimpB_5cm_YZ.jpeg" attr="" comment="" date="1519920220" name="MBimpB_5cm_YZ.jpeg" path="MBimpB_5cm_YZ.jpeg" size="24644" user="AlessioSarti" version="2"
- META FILEATTACHMENT
+ attachment="MBcil_0cm_zoom.png" attr="" comment="" date="1519920358" name="MBcil_0cm_zoom.png" path="MBcil_0cm_zoom.png" size="51648" user="AlessioSarti" version="1"
- META FILEATTACHMENT
+ attachment="MBcil_0cm.png" attr="" comment="" date="1519920358" name="MBcil_0cm.png" path="MBcil_0cm.png" size="113631" user="AlessioSarti" version="1"
- META FILEATTACHMENT
+ attachment="MBfull_0cm_zoom.png" attr="" comment="" date="1519920358" name="MBfull_0cm_zoom.png" path="MBfull_0cm_zoom.png" size="66281" user="AlessioSarti" version="1"
- META FILEATTACHMENT
+ attachment="MBfull_0cm.png" attr="" comment="" date="1519920358" name="MBfull_0cm.png" path="MBfull_0cm.png" size="113672" user="AlessioSarti" version="1"

Revision 42018-02-28 - IlariaMattei

 			 META TOPICPARENT 
			 name="DoseProfiler" 
		
 Dose Profiler matter studies page
- META TOPICPARENT
+ name="DoseProfiler"
-<
<
+  Dose Profiler matter studies page
   Rando Full Simulation studies
   The weights definition + machinery
   How to run simulations and produce files to apply weights:
   The mental ball configurations 
  Comments
->
>
+  Dose Profiler matter studies page
   Rando Full Simulation studies
   The weights definition + machinery
   How to run simulations and produce files to apply weights:
   The mental ball configurations 
  Comments
 This page documents the matter studies.

 The weights definition + machinery 

Here we document how the weights are defined, where is the code

 How to run simulations and produce files to apply weights:
-<
<
+To Run V2:
->
>
+To Run V2:
-<
<
+-reference cilinder, but not necessary-
$FLUKA/flutil/rfluka -e fluka_FILTER_mbcil.exe -N0 -M5 filter_mbcil_1R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbcil.exe -N0 -M5 filter_mbcil_2R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbcil.exe -N0 -M5 filter_mbcil_3R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbcil.exe -N0 -M5 filter_mbcil_4R &
->
>
+-reference cilinder, necessary if the target does not start at 0cm (as reference cilinder3 does) -
$FLUKA/flutil/rfluka -e fluka_FILTER_mbcil.exe -N0 -M5 filter_mbcil_1R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbcil.exe -N0 -M5 filter_mbcil_2R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbcil.exe -N0 -M5 filter_mbcil_3R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbcil.exe -N0 -M5 filter_mbcil_4R &
 -mb brain+bone radius_ext = 10 cm centered in (0,0,10)
$FLUKA/flutil/rfluka -e fluka_FILTER_mbfull.exe -N0 -M5 filter_mbfull_1R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbfull.exe -N0 -M5 filter_mbfull_2R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbfull.exe -N0 -M5 filter_mbfull_3R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbfull.exe -N0 -M5 filter_mbfull_4R &

cat filter_mbcil*00*dat > filter_mbcil_all_TXT.dat &
cat filter_mbfull*00*dat > filter_mbfull_all_TXT.dat &

./Txt2Root -in filter_mbcil_all_TXT.dat -out filter_mbcil_all.root &

Then in ../ApplyWeights:
-<
<
+see ListApplyWeights _mbfull.sh (first part)->to have reference distr (not necessary)
->
>
+see ListApplyWeights _mbfull.sh (the first part) ->to have reference hZ distr
 Then:

cd ../V5

ln -s ../V2/filter_mbfull_all_TXT.dat secondary_TXT.dat

and run FLUKA on new ray*inp files (same full geo, different beam (10TeV muon+) and physics card) 
NB: check the source for reso!
$FLUKA/flutil/rfluka -e fluka_RAY.exe -N0 -M1 ray_filter_mbfull_reso0 &
$FLUKA/flutil/rfluka -e fluka_RAY.exe -N0 -M1 ray_filter_mbfull_reso1 &
$FLUKA/flutil/rfluka -e fluka_RAY.exe -N0 -M1 ray_filter_mbfull_reso2 &
$FLUKA/flutil/rfluka -e fluka_RAY.exe -N0 -M1 ray_filter_mbfull_reso3 &

./Txt2Root -in ray_filter_mbfull_reso0001_TXT.dat -out ray_filter_mbfull_reso0.root &
./Txt2Root -in ray_filter_mbfull_reso1001_TXT.dat -out ray_filter_mbfull_reso1.root &
./Txt2Root -in ray_filter_mbfull_reso2001_TXT.dat -out ray_filter_mbfull_reso2.root &
./Txt2Root -in ray_filter_mbfull_reso3001_TXT.dat -out ray_filter_mbfull_reso3.root &

Then in ApplyWeights:

see ListApplyWeights _mbfull.sh (second part)
->
>
+to produce the weighted Z distr
 - Compare the distributions with TestWeighs _mb.cc o quello che hai <3
-<
<
+ The mental ball configuration
->
>
+ The mental ball configurations
-<
<
+) Mental Ball: rmin = 9.6 cm, rmax = 10 cm, inner: BRAIN (ICRU), outer: BONE (ICRU)
->
>
+Beam: 12C ions, energy: 300 MeV /u, spatial distribution: gaussian with FWHM_x,y = 0.5 cm
-<
<
+center:(0,0,10) cm
->
>
+direction: positive z axis, starting point:(0,0,4) cm
-<
<
+Beam: 12C ions, energy: 300 MeV /u, spatial distribution: gaussian with FWHM_x,y = 0.5 cm
->
>
+)*mbfull, REFERENCE distr: mbcil*
-<
<
+direction: positive z axis, starting point:(0,0,-1) cm
->
>
+Mental Ball: sphere with rmin = 9.6 cm, rmax = 10 cm, inner: BRAIN (ICRU), outer: BONE (ICRU)
-<
<
+--   Alessio Sarti  - 2018-02-22
->
>
+center:(0,0,15) cm
-<
<
+ Comments
->
>
+- fluka region number of detectors (to be implemented in ApplyWeightsLib.h): {3,4,5,6,7,8}
-<
<
+  (Edit - Preview)
 
<--/commentPlugin-->
->
>
+REFERENCE distr: mbcil = MB with same geo, material: water, cut with cilinder of r = 1.5 cm
->
>
+)*mbimp, REFERENCE distr: mbimp_cil*

Mental Ball: sphere with rmin = 9.6 cm, rmax = 10 cm, inner: BRAIN (ICRU), outer: BONE (ICRU)

center:(0,0,15) cm

Implant: sphere, r = 1cm, material: AIR, center: (0.,5.,16.)  

- fluka region number of detectors (to be implemented in ApplyWeightsLib.h): {0,0,0,3,4,5} 

 NB: the 60deg detectors are cut in order to have just a portion of ~ 20 x 20 cm2 ; no 90deg detectors in this configuration

REFERENCE distr: mbimp_cil = MB with same geo, material: water, cut with cilinder of r = 1.5 cm (as in sim 1) ); detectors in same config of mbimp sim

3)*mbimp2, REFERENCE distr: mbimp_cil*

Mental Ball: sphere with rmin = 9.6 cm, rmax = 10 cm, inner: BRAIN (ICRU), outer: BONE (ICRU)

center:(0,0,15) cm

Implant: sphere, r = 2cm, material: AIR, center: (0.,5.,16.)  

- fluka region number of detectors (to be implemented in ApplyWeightsLib.h): {0,0,0,3,4,5} 

 NB: the 60deg detectors are cut in order to have just a portion of ~ 20 x 20 cm2 ; no 90deg detectors in this configuration

REFERENCE distr: mbimp_cil = MB with same geo, material: water, cut with cilinder of r = 1.5 cm (as in sim 1) ); detectors in same config of mbimp sim

4)*mbimp3, REFERENCE distr: mbimp_cil*

Mental Ball: sphere with rmin = 9.6 cm, rmax = 10 cm, inner: BRAIN (ICRU), outer: BONE (ICRU)

center:(0,0,15) cm

Implant: sphere, r = 1cm, material: ADIPOSE TISSUE, center: (0.,5.,16.)  

- fluka region number of detectors (to be implemented in ApplyWeightsLib.h): {0,0,0,3,4,5} 

 NB: the 60deg detectors are cut in order to have just a portion of ~ 20 x 20 cm2 ; no 90deg detectors in this configuration

REFERENCE distr: mbimp_cil = MB with same geo, material: water, cut with cilinder of r = 1.5 cm (as in sim 1) ); detectors in same config of mbimp sim

 Comments 


 
  (Edit - Preview)
 
<--/commentPlugin-->

Revision 32018-02-23 - IlariaMattei

 			 META TOPICPARENT 
			 name="DoseProfiler" 
		
 Dose Profiler matter studies page
- META TOPICPARENT
+ name="DoseProfiler"
-<
<
+  Dose Profiler matter studies page
   Rando Full Simulation studies
   The weights definition + machinery
   How to run simulations and produce files to apply weights:
   The mental ball configurations 
  Comments
->
>
+  Dose Profiler matter studies page
   Rando Full Simulation studies
   The weights definition + machinery
   How to run simulations and produce files to apply weights:
   The mental ball configurations 
  Comments
 This page documents the matter studies.

 The weights definition + machinery 

Here we document how the weights are defined, where is the code
->
>
+ How to run simulations and produce files to apply weights: 

To Run V2:

-reference cilinder, but not necessary-
$FLUKA/flutil/rfluka -e fluka_FILTER_mbcil.exe -N0 -M5 filter_mbcil_1R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbcil.exe -N0 -M5 filter_mbcil_2R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbcil.exe -N0 -M5 filter_mbcil_3R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbcil.exe -N0 -M5 filter_mbcil_4R &

-mb brain+bone radius_ext = 10 cm centered in (0,0,10)
$FLUKA/flutil/rfluka -e fluka_FILTER_mbfull.exe -N0 -M5 filter_mbfull_1R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbfull.exe -N0 -M5 filter_mbfull_2R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbfull.exe -N0 -M5 filter_mbfull_3R &
$FLUKA/flutil/rfluka -e fluka_FILTER_mbfull.exe -N0 -M5 filter_mbfull_4R &

cat filter_mbcil*00*dat > filter_mbcil_all_TXT.dat &
cat filter_mbfull*00*dat > filter_mbfull_all_TXT.dat &

./Txt2Root -in filter_mbcil_all_TXT.dat -out filter_mbcil_all.root &

Then in ../ApplyWeights:

see ListApplyWeights _mbfull.sh (first part)->to have reference distr (not necessary)

Then:

cd ../V5

ln -s ../V2/filter_mbfull_all_TXT.dat secondary_TXT.dat

and run FLUKA on new ray*inp files (same full geo, different beam (10TeV muon+) and physics card) 
NB: check the source for reso!
$FLUKA/flutil/rfluka -e fluka_RAY.exe -N0 -M1 ray_filter_mbfull_reso0 &
$FLUKA/flutil/rfluka -e fluka_RAY.exe -N0 -M1 ray_filter_mbfull_reso1 &
$FLUKA/flutil/rfluka -e fluka_RAY.exe -N0 -M1 ray_filter_mbfull_reso2 &
$FLUKA/flutil/rfluka -e fluka_RAY.exe -N0 -M1 ray_filter_mbfull_reso3 &

./Txt2Root -in ray_filter_mbfull_reso0001_TXT.dat -out ray_filter_mbfull_reso0.root &
./Txt2Root -in ray_filter_mbfull_reso1001_TXT.dat -out ray_filter_mbfull_reso1.root &
./Txt2Root -in ray_filter_mbfull_reso2001_TXT.dat -out ray_filter_mbfull_reso2.root &
./Txt2Root -in ray_filter_mbfull_reso3001_TXT.dat -out ray_filter_mbfull_reso3.root &

Then in ApplyWeights:

see ListApplyWeights _mbfull.sh (second part)

- Compare the distributions with TestWeighs _mb.cc o quello che hai <3
  The mental ball configuration 

1) Mental Ball: rmin = 9.6 cm, rmax = 10 cm, inner: BRAIN (ICRU), outer: BONE (ICRU)
-<
<
+  center:(0,0,10) cm
->
>
+center:(0,0,10) cm
-<
<
+  Beam: 12C ions, energy: 300 MeV/u, spatial distribution: gaussian with FWHM_x,y = 0.5 cm
->
>
+Beam: 12C ions, energy: 300 MeV /u, spatial distribution: gaussian with FWHM_x,y = 0.5 cm
-<
<
+        direction: positive z axis, starting point:(0,0,-1) cm
->
>
+direction: positive z axis, starting point:(0,0,-1) cm
-<
<
+--   Alessio Sarti  - 2018-02-22
->
>
+--   Alessio Sarti  - 2018-02-22
  Comments
-<
<
+  (Edit - Preview)
 
<--/commentPlugin-->
->
>
+  (Edit - Preview)
 
<--/commentPlugin-->
->
>

Revision 22018-02-23 - IlariaMattei

 			 META TOPICPARENT 
			 name="DoseProfiler" 
		
 Dose Profiler matter studies page 

 
  Dose Profiler matter studies page
   Rando Full Simulation studies
   The weights definition + machinery
   How to run simulations and produce files to apply weights:
   The mental ball configurations 
  Comments
 
 


This page documents the matter studies.

 The weights definition + machinery 

Here we document how the weights are defined, where is the code

 The mental ball configuration
- META TOPICPARENT
+ name="DoseProfiler"
-<
<
+What samples are available, what statistics
->
>
+) Mental Ball: rmin = 9.6 cm, rmax = 10 cm, inner: BRAIN (ICRU), outer: BONE (ICRU)
-<
<
+--   Alessio Sarti  - 2018-02-22
->
>
+  center:(0,0,10) cm
->
>
+  Beam: 12C ions, energy: 300 MeV/u, spatial distribution: gaussian with FWHM_x,y = 0.5 cm

        direction: positive z axis, starting point:(0,0,-1) cm

--   Alessio Sarti  - 2018-02-22
  Comments 


 
  (Edit - Preview)
 
<--/commentPlugin-->
-<
<

Revision 12018-02-22 - AlessioSarti

META TOPICPARENT	name="DoseProfiler"

Dose Profiler matter studies page

Dose Profiler matter studies page
Rando Full Simulation studies
The weights definition + machinery
How to run simulations and produce files to apply weights:
The mental ball configurations
- Comments

This page documents the matter studies.

The weights definition + machinery

Here we document how the weights are defined, where is the code

The mental ball configuration

What samples are available, what statistics

-- Alessio Sarti - 2018-02-22

Comments

<--/commentPlugin-->

View topic | History: r7 < r6 < r5 < r4 | More topic actions...