User Tools

Site Tools


fcs:eic:intro

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
fcs:eic:intro [2016/11/04 19:50]
asc [Reconstructed events]
fcs:eic:intro [2016/12/09 16:34] (current)
asc [Hit timing information]
Line 111: Line 111:
 Note that the plot is done for all electrons starting from 1 GeV (a region where the efficiency is quite small). Note that the plot is done for all electrons starting from 1 GeV (a region where the efficiency is quite small).
  
-====== Looking at other samples ======+ 
 +====== Hit timing information ====== 
 + 
 + 
 + 
 +You can also look at timing of silicon tracker using this script: {{:fcs:eic:mc_hits_time.py|eic:mc_hits_time.py}} 
 +<code bash> 
 +hs-get pgun_elec_eta4_pt30gev%rfull052 data # 2 100  get 100 files in 2 threads in the directory "data" 
 +fpad eic:mc_hits_time.py 
 +</code> 
 +The plot is shown here: 
 + 
 +{{:fcs:eic:electrons_hit_timing.png?400|}} 
 + 
 +====== Looking at track resolution ======
  
 You can also look at single particles with flat pT, Eta and Phi. Here is an example this using a simple validation script {{:fcs:eic:mc_truth_flat.py|}} You can also look at single particles with flat pT, Eta and Phi. Here is an example this using a simple validation script {{:fcs:eic:mc_truth_flat.py|}}
-for the MCparticle container inside LCIO file:+for the MCparticle container inside LCIO file. We will use the "flat pT" sample called [[http://atlaswww.hep.anl.gov/hepsim/info.php?item=242|pgun_elec_eta4_pt30gev]] reconstructed with the tag "rfull052" (SiEIC with 1T solenoid):
  
-<code>+<code bash>
 hs-get pgun_elec_eta4_pt30gev%rfull052 data  2 100 # get 100 files in 2 threads in the directory "data" hs-get pgun_elec_eta4_pt30gev%rfull052 data  2 100 # get 100 files in 2 threads in the directory "data"
 fpad mc_truth_flat.py fpad mc_truth_flat.py
Line 123: Line 137:
 The resulting image is: The resulting image is:
  
-{{:fcs:eic:mc_truth_pt_ele.png?400|}}+{{:fcs:eic:mc_truth_pt_ele.png?500|}} 
 + 
 +Now, let's calculate track resolution. We will use this script {{:fcs:eic:track_resolution.py|track_resolution.py}} which 
 +makes fit of pT(reco)/pT(truth) and plots the resolution as a function of pT(truth). 
 +This command shown below makes resolution plots, fits with a Gaussian function, and plot resolution as a function of pT: 
 +<code bash> 
 +hs-get pgun_elec_eta4_pt30gev%rfull052 data  2 20 # get 20 files in 2 threads in the directory "data" 
 +fpad track_resolution.py 
 +</code> 
 + 
 +The final plot is shown below: 
 + 
 +{{:fcs:eic:pt_resolution.png?400|}} 
  
 ====== How to continue ====== ====== How to continue ======
Line 135: Line 162:
   * fake rate (cannot be done with this sample!)   * fake rate (cannot be done with this sample!)
   * calculate Q2 and x_bj, check efficiency and resolutions for these variable   * calculate Q2 and x_bj, check efficiency and resolutions for these variable
-===== Displaying reconstructed events =====+====== Displaying reconstructed events ======
  
  
Line 144: Line 171:
  
 <code bash> <code bash>
-wget http://mc.hep.anl.gov/asc/hepsim/events/misc/pgun/pgun_eta45_p30_elec//rfull051/pgun_elec30gev_001_pandora.slcio+wget http://mc.hep.anl.gov/asc/hepsim/events/misc/pgun/pgun_eta45_p30_elec/rfull051/pgun_elec30gev_001_pandora.slcio
 </code> </code>
  
fcs/eic/intro.1478289015.txt.gz · Last modified: 2016/11/04 19:50 by asc