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/06 16:26]
asc [Looking at track resolution]
fcs:eic:intro [2016/12/09 16:34] (current)
asc [Hit timing information]
Line 110: Line 110:
  
 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).
 +
 +
 +====== 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 ====== ====== 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. 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":+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 bash> <code bash>
Line 126: Line 140:
  
 Now, let's calculate track resolution. We will use this script {{:fcs:eic:track_resolution.py|track_resolution.py}} which 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): +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> <code bash>
 hs-get pgun_elec_eta4_pt30gev%rfull052 data  2 20 # get 20 files in 2 threads in the directory "data" hs-get pgun_elec_eta4_pt30gev%rfull052 data  2 20 # get 20 files in 2 threads in the directory "data"
Line 133: Line 147:
 </code> </code>
  
 +The final plot is shown below:
  
 {{:fcs:eic:pt_resolution.png?400|}} {{:fcs:eic:pt_resolution.png?400|}}
Line 156: 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.1478449618.txt.gz · Last modified: 2016/11/06 16:26 by asc