User Tools

Site Tools


vlhc:hcal

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
vlhc:hcal [2014/11/21 17:15]
asc [To do next:]
vlhc:hcal [2015/10/17 15:33] (current)
asc
Line 1: Line 1:
  
-====== HCAL studies using a fast simulation ======+====== Calorimeter studies using a fast simulation ======
  
-Use [[http://atlaswww.hep.anl.gov/hepsim/HepSim | HepSim]] repository for truth-level MC and run a fast simulation on the fly to study HCAL segmentation as explained in  [[https://atlaswww.hep.anl.gov/asc/wikidoc/doku.php?id=community:refhepsim |this manual]].  +[[http://atlaswww.hep.anl.gov/hepsim/ | HepSim]] repository provides a number of truth-level MC samples for signal and background that can be used as inputs for a fast simulation. The analysis can be done as explained in  [[https://atlaswww.hep.anl.gov/asc/wikidoc/doku.php?id=community:refhepsim |this manual]].  
-Follow these steps:+To simplify the setup, we have prepared a package that includes [[https://cp3.irmp.ucl.ac.be/projects/delphes | Delphes]], [[http://atlaswww.hep.anl.gov/asc/promc/ | ProMC]] and HS-TOOLs. Follow these steps:
  
 <code bash> <code bash>
-wget http://atlaswww.hep.anl.gov/asc/hepsim/FastHepSim.tgz+wget http://atlaswww.hep.anl.gov/asc/hepsim/soft/FastHepSim.tgz
 tar -zvxf FastHepSim.tgz tar -zvxf FastHepSim.tgz
 cd FastHepSim/libraries/ cd FastHepSim/libraries/
Line 25: Line 25:
 </code> </code>
  
-This compiles the analysis program (analysis.cc) that fills pT and jet mass histograms. Before running this program, we will describe what detector geometry can be studied.+This compiles the analysis program (analysis.cc) that fills pT and jet mass histograms. Before running this program, we will describe what detector geometry can be studied. They are explained in the sections below.
  
 ===== ATLAS-like geometry ===== ===== ATLAS-like geometry =====
Line 34: Line 34:
 </code> </code>
  
-Creates jets from "combined towers" from HCAL and ECAL towers. Then we apply a simple scale correction for jets  1.06+400/(pt*pt*) by looking at resolution plots.+Creates jets from "combined towers", i.e. from HCAL and ECAL towers.  This card also apples a simple scale correction for jets  1.06+400/(pt*pt*) (this was done after a study of jet resolution plots)
  
 For the initial test, we assume the following segmentation: For the initial test, we assume the following segmentation:
Line 40: Line 40:
   * HCAL has 64 modules in Phi and Eta=0.1 in the central region   * HCAL has 64 modules in Phi and Eta=0.1 in the central region
   * ECAL has  x4 better segmentation   * ECAL has  x4 better segmentation
 +  * 60% of energy of charged pions and all baryons are in HCAL 
  
 This  HCAL geometry has a well balanced segmentation in 2D (same angle of 5.625 deg in Phi and Theta). This  HCAL geometry has a well balanced segmentation in 2D (same angle of 5.625 deg in Phi and Theta).
Line 56: Line 57:
 #   -6.0<Eta<4.6      and 4.6<Eta<   Segment: 0.2 eta x 20 deg phi #   -6.0<Eta<4.6      and 4.6<Eta<   Segment: 0.2 eta x 20 deg phi
 </code> </code>
-Electrons and photons (pi0!) do not leak to HCAL. 40% of charged pions and all baryons deposit their energy in HCAL.  +Electrons and photons (pi0!) do not leak to HCAL.  
  
  
Line 150: Line 151:
 ===== How to run ===== ===== How to run =====
  
-We need to bring truth-level data from [[http://atlaswww.hep.anl.gov/hepsim/| HepSim]] and put them somewhere. For this example, +We need to bring truth-level data from [[http://atlaswww.hep.anl.gov/hepsim/| HepSim]] and put them somewhere.  
-we will use [[https://atlaswww.hep.anl.gov/hepsim/info.php?item=57 | TTBar sample]] generated with pT>2.5 TeV and look at the jet mass.+We will use [[https://atlaswww.hep.anl.gov/hepsim/info.php?item=57 | TTBar sample]] generated with pT>2.5 TeV and look at the jet mass.
 Go to "analysis" directory and copy a few truth-level files: Go to "analysis" directory and copy a few truth-level files:
 <code bash> <code bash>
 hs-get http://mc.hep.anl.gov/asc/hepsim/events/pp/100tev/ttbar_pt2500_mg5  data 2 10 hs-get http://mc.hep.anl.gov/asc/hepsim/events/pp/100tev/ttbar_pt2500_mg5  data 2 10
 </code> </code>
-This copies 10 files in 2 threads and put them to the directory "data". Now we are ready to run over these files.+This copies 10 files in 2 threads and put them to the directory "data". Now we are ready to run over these 10 files.
  
 Let us make 2 runs: One run  will use the first geometry layout, and the second run will use the 2nd input card. Let us make 2 runs: One run  will use the first geometry layout, and the second run will use the 2nd input card.
Line 173: Line 174:
 To make plot, run {{:vlhc:mass.py|}} script. This will show 2 jets masses for boosted top  for different HCAL segmentation. To make plot, run {{:vlhc:mass.py|}} script. This will show 2 jets masses for boosted top  for different HCAL segmentation.
  
-{{:vlhc:mass_towers.png}}+ 
 + 
 +{{:vlhc:mass_towers.eps.png?400|}} 
  
 Note: the eflow option has much better resolution for jet masses. Note: the eflow option has much better resolution for jet masses.
  
-===== Using particle flows =====+===== Using particle flows (eFlow) =====
  
  
Line 193: Line 197:
  
  
 +Now let's try to run using "eFlow" cards:
 +
 +{{:vlhc:mass_eflow.eps.png?400|}}
 +
 +{{:vlhc:tau32_eflow.eps.png?400|}}
  
 ===== Next steps ===== ===== Next steps =====
vlhc/hcal.1416590147.txt.gz · Last modified: 2014/11/21 17:15 by asc