1 / 8

Stuff on MPC simulations for Ken/Matt

Stuff on MPC simulations for Ken/Matt. R. Seto march 29. eres. med ium. bad. good. define eres=(e_prim-clus_e)/e_prim. ebad. this one must reflect mickeys thrown x and y distribution. emed. egood. this shows the edges of the MPC. define fiducial as 12<r<18. eorig. eres.

ofira
Download Presentation

Stuff on MPC simulations for Ken/Matt

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Stuff on MPC simulationsfor Ken/Matt R. Seto march 29

  2. eres med ium bad good define eres=(e_prim-clus_e)/e_prim ebad this one must reflect mickeys thrown x and y distribution emed egood this shows the edges of the MPC define fiducial as 12<r<18

  3. eorig eres eres vs e for R<r_fiducial • slice up plot into bins of e_prim and fit width • fit the widths to • delta E/E_prim = A+B/sqrt(E_prim) • fit for A and B • One slice fit for E_prim=25 GeV • delta E=0.03 • sqrt(25)=5 • B~15% if A is 0 • Ken and matt should do a full fit slice E_prim=25

  4. x,y prim for e_prim 0-5 GeVand eres in “medium” areashows leakage • to get this plot I looked at data with e=0-5 GeV. I think this better because of the smaller radius of the showers?

  5. delta x =x_prim-x_clus • For egood • centered at 0.5, 0.5 • peak is not in center • Shows structure of MPC modules • RMS ~ 1cm

  6. same plot for emed • peak now in center • RMS ~ 1.5-2 (wider)

  7. snippet of code • float r=sqrt(prim_x*prim_x+prim_y*prim_y); • h_cluse->Fill(clus_e); • h_prime->Fill(prim_e); • h_eres->Fill((prim_e-clus_e)/prim_e); • if(r>12 && r<18)h_eres_vs_e->Fill((prim_e-clus_e)/prim_e,prim_e); • h_xy->Fill(prim_x,prim_y); • h_clus_xy->Fill(clus_x,clus_y); • if( (prim_e-clus_e)/prim_e < 0.1){ • h_xy_egood->Fill(prim_x,prim_y); • h_clus_xy_egood->Fill(clus_x,clus_y); • h_delta_xy_egood->Fill(prim_x-clus_x,prim_y-clus_y); • h_delta_x_egood_vs_e->Fill(prim_x-clus_x,prim_e); • } • if( (prim_e-clus_e)/prim_e > 0.1 && (prim_e-clus_e)/prim_e < 0.9 ){ • h_xy_emed->Fill(prim_x,prim_y); • h_clus_xy_emed->Fill(clus_x,clus_y); • h_delta_xy_emed->Fill(prim_x-clus_x,prim_y-clus_y); • } • if( (prim_e-clus_e)/prim_e > 0.9){ • h_xy_ebad->Fill(prim_x,prim_y); • h_clus_xy_ebad->Fill(clus_x,clus_y); • } • }

  8. notes • I made the pictures just using the output file and root • you should read the data files from mickey as he has set up. The filelist will let you read many files – this is the fun4all mode (which is all set up for you) – the ability to read many files is something we will have to figure out (its not hard) • In my case I used a TChain – I cheated, but it wont allow me to a access stiff if its not in the ttree

More Related