Difference between revisions of "TF Simul PhotField IAC 2n"
| Line 9: | Line 9: | ||
| − | + | 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]] | [[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");