Difference between revisions of "TF Simul PhotField IAC 2n"
(→Sim1) |
|||
(One intermediate revision by the same user not shown) | |||
Line 27: | Line 27: | ||
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"); | 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