Difference between revisions of "TF Simul PhotField IAC 2n"

From New IAC Wiki
Jump to navigation Jump to search
 
(3 intermediate revisions by the same user not shown)
Line 2: Line 2:
  
  
 +
Electrons impinge a 3 x 10 x 3 cm aluminum radiator.  The 10.5 MeV electrons radiate in the air and the aluminum before stopping in the aluminum.
  
root [10] Phot2n->Draw("evt.p1x*evt.p1x+evt.p1y*evt.p1y+evt.p1z*evt.p1z","evt.p1x*evt.p1x+evt.p1y*evt.p1y+evt.p1z*evt.p1z>0");
+
Photons are emitted in all directions. A 6 foot high wall made of 36" thick concrete is located downstream of the Aluminum radiator.
root [11] Phot2n->Draw("sqrt(evt.p1x*evt.p1x+evt.p1y*evt.p1y+evt.p1z*evt.p1z)","evt.p1x*evt.p1x+evt.p1y*evt.p1y+evt.p1z*evt.p1z>0");
+
 
root [12] Phot2n->Draw("sqrt(evt.p2x*evt.p2x+evt.p2y*evt.p2y+evt.p2z*evt.p2z)","evt.p2x*evt.p2x+evt.p2y*evt.p2y+evt.p2z*evt.p2z>0");
+
A 30 cm radius lead collimator is position within the wall.
root [13] Phot2n->Draw("sqrt(evt.p3x*evt.p3x+evt.p3y*evt.p3y+evt.p3z*evt.p3z)","evt.p3x*evt.p3x+evt.p3y*evt.p3y+evt.p3z*evt.p3z>0");
+
 
 +
 
 +
Phot2n->Draw("sqrt(evt.p1x*evt.p1x+evt.p1y*evt.p1y+evt.p1z*evt.p1z)>>PhotonSource(100,0,10)","evt.p1x*evt.p1x+evt.p1y*evt.p1y+evt.p1z*evt.p1z>0");
 +
 
 +
Phot2n->Draw("sqrt(evt.p2x*evt.p2x+evt.p2y*evt.p2y+evt.p2z*evt.p2z)>>PhotonsOnWall(100,0,10)","evt.p2x*evt.p2x+evt.p2y*evt.p2y+evt.p2z*evt.p2z>0");
 +
 
 +
 
 +
Phot2n->Draw("sqrt(evt.p1x*evt.p1x+evt.p1y*evt.p1y+evt.p1z*evt.p1z)>>PhotonSource(100,0,10)","evt.p1x*evt.p1x+evt.p1y*evt.p1y+evt.p1z*evt.p1z>0");
 +
 
 +
 
 +
Phot2n->Draw("sqrt(evt.p2x*evt.p2x+evt.p2y*evt.p2y+evt.p2z*evt.p2z)","evt.p2x*evt.p2x+evt.p2y*evt.p2y+evt.p2z*evt.p2z>0");
 +
 
 +
 
 +
Phot2n->Draw("sqrt(evt.p3x*evt.p3x+evt.p3y*evt.p3y+evt.p3z*evt.p3z)","evt.p3x*evt.p3x+evt.p3y*evt.p3y+evt.p3z*evt.p3z>0");
 +
 
 +
 
 +
Phot2n->Draw("sqrt(evt.p3x*evt.p3x+evt.p3y*evt.p3y+evt.p3z*evt.p3z)>>PhotonsInCell(100,0,10)","evt.p3x*evt.p3x+evtroot [26] 3y+evt.p3z*evt.p3z>0 && abs(evt.pos3x)<1000 && abs(evt.pos3y)<1000");
 +
 
 +
 
 +
root [26] Phot2n->Draw("evt.pos3x:evt.pos3y","evt.p3x*evt.p3x+evt.p3y*evt.p3y+evt.p3z*evt.p3z>0 && abs(evt.pos3x)<1500 && abs(evt.pos3y)<1500");
 +
 
 +
 
 +
==Sim1==
 +
 
 +
Lead Collimator is 15 cm in radius with a 1/4" diameter hole.
 +
Aluminum radiator 2cm thick
  
  
 
[[Simulation_TF_Ncor]]
 
[[Simulation_TF_Ncor]]

Latest revision as of 22:44, 1 July 2014

The goal of this simulation is to determine the photon field in the shielded cell that result from 10.5 MeV electrons hitting 2 cm thick Aluminum radiator.


Electrons impinge a 3 x 10 x 3 cm aluminum radiator. The 10.5 MeV electrons radiate in the air and the aluminum before stopping in the aluminum.

Photons are emitted in all directions. A 6 foot high wall made of 36" thick concrete is located downstream of the Aluminum radiator.

A 30 cm radius lead collimator is position within the wall.


Phot2n->Draw("sqrt(evt.p1x*evt.p1x+evt.p1y*evt.p1y+evt.p1z*evt.p1z)>>PhotonSource(100,0,10)","evt.p1x*evt.p1x+evt.p1y*evt.p1y+evt.p1z*evt.p1z>0");

Phot2n->Draw("sqrt(evt.p2x*evt.p2x+evt.p2y*evt.p2y+evt.p2z*evt.p2z)>>PhotonsOnWall(100,0,10)","evt.p2x*evt.p2x+evt.p2y*evt.p2y+evt.p2z*evt.p2z>0");


Phot2n->Draw("sqrt(evt.p1x*evt.p1x+evt.p1y*evt.p1y+evt.p1z*evt.p1z)>>PhotonSource(100,0,10)","evt.p1x*evt.p1x+evt.p1y*evt.p1y+evt.p1z*evt.p1z>0");


Phot2n->Draw("sqrt(evt.p2x*evt.p2x+evt.p2y*evt.p2y+evt.p2z*evt.p2z)","evt.p2x*evt.p2x+evt.p2y*evt.p2y+evt.p2z*evt.p2z>0");


Phot2n->Draw("sqrt(evt.p3x*evt.p3x+evt.p3y*evt.p3y+evt.p3z*evt.p3z)","evt.p3x*evt.p3x+evt.p3y*evt.p3y+evt.p3z*evt.p3z>0");


Phot2n->Draw("sqrt(evt.p3x*evt.p3x+evt.p3y*evt.p3y+evt.p3z*evt.p3z)>>PhotonsInCell(100,0,10)","evt.p3x*evt.p3x+evtroot [26] 3y+evt.p3z*evt.p3z>0 && abs(evt.pos3x)<1000 && abs(evt.pos3y)<1000");


root [26] Phot2n->Draw("evt.pos3x:evt.pos3y","evt.p3x*evt.p3x+evt.p3y*evt.p3y+evt.p3z*evt.p3z>0 && abs(evt.pos3x)<1500 && abs(evt.pos3y)<1500");


Sim1

Lead Collimator is 15 cm in radius with a 1/4" diameter hole. Aluminum radiator 2cm thick


Simulation_TF_Ncor