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

From New IAC Wiki
Jump to navigation Jump to search
Line 9: Line 9:
  
  
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");
+
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");
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");
 
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.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");
  
 
[[Simulation_TF_Ncor]]
 
[[Simulation_TF_Ncor]]

Revision as of 23:48, 30 June 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");

Simulation_TF_Ncor