ATC: ATLAS Center at ANL

User Tools

Site Tools

Trace:
community:refhepsim

Differences

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

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision		Previous revision
community:refhepsim [2016/03/21 22:19] – asccommunity:refhepsim [2016/04/28 18:23] (current) – [On the fly reconstruction] asc
Line 1: Line 1:
 +[[http://atlaswww.hep.anl.gov/hepsim/|<< back to HepSim]]
  
 ======  =HepSim= ====== ======  =HepSim= ======
Line 11: Line 12:
  
 The Monte Carlo files  can be accessed from the [[http://atlaswww.hep.anl.gov/hepsim/ | HepSim repository]] link.  The Monte Carlo files  can be accessed from the [[http://atlaswww.hep.anl.gov/hepsim/ | HepSim repository]] link. 
 +
 +
 +======  Physics and detector studies  ======
  
  
Line 16: Line 20:
    
  
-  * [[fcs:fcchh:intro|FCC-hh detector studies]] - explains how to analyses data for FCC-hh detector studies+  * [[fcs:fcchh:|FCC-hh detector studies]] - explains how to analyses data for FCC-hh detector studies
   * [[fcs:sid:intro|SID detector studies]] - explains how to analyses data for the SiD detector (ILC)   * [[fcs:sid:intro|SID detector studies]] - explains how to analyses data for the SiD detector (ILC)
   * [[fcs:cepc:intro|CEPC detector studies]] - shows some results with full simulations for CEPC   * [[fcs:cepc:intro|CEPC detector studies]] - shows some results with full simulations for CEPC
Line 22: Line 26:
   * [[vlhc:hcal|HCAL studies]] explains how to analyse ROOT data after fast detector simulations used for FCC studies   * [[vlhc:hcal|HCAL studies]] explains how to analyse ROOT data after fast detector simulations used for FCC studies
  
-Currently these studies are based on  [[https://atlaswww.hep.anl.gov/asc/jas4pp|Jas4pp data-analysis environment]].+Currently these studies are based on  [[http://atlaswww.hep.anl.gov/asc/jas4pp|Jas4pp data-analysis environment]]. 
 +Geometry files for full detector simulations can be found in the  [[http://atlaswww.hep.anl.gov/hepsim/soft/detectors/|detector repository]].
  
 This manual describes how to work with truth-level samples posted on HepSim. This manual describes how to work with truth-level samples posted on HepSim.
-[[https://atlaswww.hep.anl.gov/asc/wikidoc/doku.php?id=community:refhepsim_analysis|HepSim Python/Java analysis page]] explains how to write code to read truth-level files using Python on the Java platform.+[[community:refhepsim_analysis|HepSim Python/Java analysis page]] explains how to write code to read truth-level files using Python on the Java platform. 
 + 
 +<fc #008000>We are working on a  more extensive manual [[community:hepsim|here]]</fc>.
  
 ======  Data formats ====== ======  Data formats ======
Line 63: Line 70:
 </code> </code>
  
-<note tip>The programs included to **hs-toolkit** can be used for searching and downloading ProMC files, and other files after fast and full simulations. +<note tip>The programs included in **hs-toolkit** can be used for searching and downloading ProMC files, and other files after fast and full simulations. 
-If you need to analyze and display reconstructed events, use the [[https://atlaswww.hep.anl.gov/asc/wikidoc/doku.php?id=fcs:fpad|FPadSoft]] program which includes **hs-toolkit**, but also adds packages for processing SLCIO files.+If you need to analyse and display reconstructed events in the LCIO format, use the [[https://atlaswww.hep.anl.gov/asc/jas4pp|Jas4PP program]] which includes **hs-toolkit**, but also adds packages for processing SLCIO files. This program can also process ProMC files.
 </note> </note>
  
Line 87: Line 94:
 You will see a pop-up GUI browser which displays the MC record. You can search for a given particle name, view data layouts and log files using the [Menu]: You will see a pop-up GUI browser which displays the MC record. You can search for a given particle name, view data layouts and log files using the [Menu]:
  
 +<hidden>
 {{:community:promc_browser.png| ProMC browser}} {{:community:promc_browser.png| ProMC browser}}
 +</hidden>
  
 This works for full parton-shower simulations with detailed information on particles. This works for full parton-shower simulations with detailed information on particles.
Line 228: Line 237:
 If you are interested in a specific reconstruction tag, use "%" to separate the search string and the tag name. If you are interested in a specific reconstruction tag, use "%" to separate the search string and the tag name.
 Example: Example:
-<code>+<code bash>
 hs-find pythia%rfast001 hs-find pythia%rfast001
 </code> </code>
Line 279: Line 288:
 the names are usually have the substring "pt", followed by the pT cuts. In this case, one can download such files as: the names are usually have the substring "pt", followed by the pT cuts. In this case, one can download such files as:
  
-<code>+<code bash>
 hs-get tev13_higgs_pythia8_ptbins data 2 5 pt100_ hs-get tev13_higgs_pythia8_ptbins data 2 5 pt100_
 </code> </code>
Line 304: Line 313:
 fast simulation, version 001). To download the reconstructed events for the reconstruction tag "rfast001", use this syntax: fast simulation, version 001). To download the reconstructed events for the reconstruction tag "rfast001", use this syntax:
  
-<code>+<code bash>
 hs-ls  tev100_ttbar_mg5%rfast001      # list reco files with the tag "rfast001" hs-ls  tev100_ttbar_mg5%rfast001      # list reco files with the tag "rfast001"
 hs-get tev100_ttbar_mg5%rfast001 data # download to the "data" directory hs-get tev100_ttbar_mg5%rfast001 data # download to the "data" directory
Line 314: Line 323:
 </code> </code>
 As before, one can also download the files using the URL: As before, one can also download the files using the URL:
-<code>+<code bash>
 hs-ls http://mc.hep.anl.gov/asc/hepsim/events/pp/100tev/ttbar_mg5/rfast001/ # list all files hs-ls http://mc.hep.anl.gov/asc/hepsim/events/pp/100tev/ttbar_mg5/rfast001/ # list all files
 hs-get http://mc.hep.anl.gov/asc/hepsim/events/pp/100tev/ttbar_mg5/rfast001/ data hs-get http://mc.hep.anl.gov/asc/hepsim/events/pp/100tev/ttbar_mg5/rfast001/ data
Line 333: Line 342:
  You can create plots using a number of programming languages, Java, Python, C++, Ruby, Groovy etc. Plots can be done on any platform, without modifying your system. C++ analysis programs require ROOT and Linux.  You can create plots using a number of programming languages, Java, Python, C++, Ruby, Groovy etc. Plots can be done on any platform, without modifying your system. C++ analysis programs require ROOT and Linux.
  
-Belows we will discuss how to analyse HepSim data using Java, since this approach works on any+Below we will discuss how to analyse HepSim data using Java, since this approach works on any
 platform (Linux, Mac, Windows) and does not require installation of any platform-specific program.  platform (Linux, Mac, Windows) and does not require installation of any platform-specific program. 
 As before, make sure that [[http://java.com/en/download/index.jsp| Java 7]] and above is installed (check it as "java -version"). As before, make sure that [[http://java.com/en/download/index.jsp| Java 7]] and above is installed (check it as "java -version").
Line 345: Line 354:
 hs-run ttbar_mg5.py hs-run ttbar_mg5.py
 </code> </code>
-Another approach is to use [[http://jwork.org/dmelt/ | DataMelt]] that gives more flexibility and more libraries for analysis. 
-In this example, we will run a Python script that downloads data from URL into the computer memory  and runs it in a batch mode. First, you will need an analysis code from HepSim. Look at ttbar sample from Madgraph: [[http://atlaswww.hep.anl.gov/hepsim/info.php?item=15 | ttbar_mg5]].  Find the URL location of the analysis script ("ttbar_mg5.py") located at the bottom of this page. Then copy the URL link of the file *.py using the right mouse button ("Copy URL Location"). Let us make a calculation of differential ttbar cross section: 
  
 +Another approach is to use [[http://atlaswww.hep.anl.gov/asc/jas4pp | Jas4pp]] or  [[http://jwork.org/dmelt/ | DataMelt]]. Such programs 
 +give more flexibility and more libraries for analysis.
 +In this example, we will run a Python script that downloads data from URL into the computer memory  and runs it in a batch mode. First, you will need an analysis code from HepSim. Look at ttbar sample from Madgraph: [[http://atlaswww.hep.anl.gov/hepsim/info.php?item=15 | ttbar_mg5]].  Find the URL location of the analysis script ("ttbar_mg5.py") located at the bottom of this page. Then copy the URL link of the file *.py using the right mouse button ("Copy URL Location"). Let us make a calculation of differential ttbar cross section. 
 +
 +Here is how to process the analysis using  [[http://atlaswww.hep.anl.gov/asc/jas4pp | Jas4pp]]:
 +
 +<code bash>
 +wget http://atlaswww.hep.anl.gov/asc/jas4pp/download/current.php -O jas4pp.tgz
 +tar -zvxf jas4pp.tgz
 +cd jas4pp
 +source ./setup.sh # takes 5 sec for first-time optimization
 +wget http://mc.hep.anl.gov/asc/hepsim/events/pp/100tev/ttbar_mg5/macros/ttbar_mg5.py # get the HepSim script
 +fpad ttbar_mg5.py # process it in a batch mode.
 +</code>
 +
 +
 +
 +Similarly, you can use a more complex [[http://jwork.org/dmelt/ | DataMelt]]:
 <code bash> <code bash>
 wget -O dmelt.zip http://jwork.org/dmelt/download/current.php                        # get DataMelt wget -O dmelt.zip http://jwork.org/dmelt/download/current.php                        # get DataMelt
Line 384: Line 409:
 </code> </code>
  
 +Similarly,  you can use  [[http://atlaswww.hep.anl.gov/asc/jas4pp | Jas4pp]].
 ===== Method III. Running in a GUI mode  ==== ===== Method III. Running in a GUI mode  ====
  
Line 393: Line 419:
 Run this code  by pressing "run". This approach uses a light-weight editor built-in inside the hs-tool package.   Run this code  by pressing "run". This approach uses a light-weight editor built-in inside the hs-tool package.  
  
-You can also bring up a full-featured GUI editor as this:+For Jas4pp, you can  start an editor, correct the script, and run it: 
 + 
 +<code bash> 
 +./jaspp ttbar_mg5.py # Open the script in the editor 
 +</code> 
 +Then, use the right mouse button and select "Run Python". You will see the output.  
 + 
 + 
 +For the DMelt IDE, you can also bring up a full-featured GUI editor as this:
  
 <code bash> <code bash>
Line 402: Line 436:
  
 ===== Method IV. Running in a GUI mode using URL dialog  ==== ===== Method IV. Running in a GUI mode using URL dialog  ====
 +
 +
  
 If you use DMelt, you can run this code using If you use DMelt, you can run this code using
Line 475: Line 511:
 The execution of this script will bring-up a window with the lego plot. The execution of this script will bring-up a window with the lego plot.
  
 +<hidden>
 {{:community:etaphi.png|Eta-Phi}} {{:community:etaphi.png|Eta-Phi}}
 +</hidden>
  
 You can also look at a simple "tracking"  for a given event: You can also look at a simple "tracking"  for a given event:
  
 +<hidden>
 {{:community:3dview.png| 3D view}} {{:community:3dview.png| 3D view}}
 +</hidden>
  
 The code written in Python is attached: The code written in Python is attached:
  
 +<hidden>
 {{:community:view3d.py| Python code to show event in 3D}} {{:community:view3d.py| Python code to show event in 3D}}
 +</hidden>
  
 ====== Reading NLO predictions ====== ====== Reading NLO predictions ======
Line 508: Line 548:
 On the left panel, click on the event and then look at "Event info".  It shows the integer values (idata) that encode the PDF uncertainties, while the float array ("fdata") shows other information. The first element in the float array is the weight of the event. The particle information is shown as usual (but without mother ID etc.).  On the left panel, click on the event and then look at "Event info".  It shows the integer values (idata) that encode the PDF uncertainties, while the float array ("fdata") shows other information. The first element in the float array is the weight of the event. The particle information is shown as usual (but without mother ID etc.). 
  
 +<hidden>
 {{:community:browser.png?|NLO Browser}} {{:community:browser.png?|NLO Browser}}
 +</hidden>
  
 The scripts that reconstruct cross sections are attached to the HepSim event repository. The scripts that reconstruct cross sections are attached to the HepSim event repository.
Line 525: Line 566:
 NLO event record includes 4-momenta of particles and event weights (double values). In addition, deviations form central weights are included as an array of integer values as: NLO event record includes 4-momenta of particles and event weights (double values). In addition, deviations form central weights are included as an array of integer values as:
  
- +{{:community:screenshot-26.png?300}}
-{{:community:screenshot-26.png}}+
  
 You can calculate differential cross sections using online files using this example: You can calculate differential cross sections using online files using this example:
Line 666: Line 706:
  
 ===== On the fly reconstruction ===== ===== On the fly reconstruction =====
-HepSim ca be used o create ROOT files after fast detector simulation, or +HepSim can be used to create ROOT files after fast detector simulation, or 
 one can analyse events after the Delphes fast simulation program on the fly. one can analyse events after the Delphes fast simulation program on the fly.
 The latter approach allows to make changes to  the detector geometry by the end-users and, at the the same time, perform an analysis. The latter approach allows to make changes to  the detector geometry by the end-users and, at the the same time, perform an analysis.
Line 780: Line 820:
 Currently, the full detector simulation can be done using the SLIC software. You will need to convert ProMC files to LCIO files and Currently, the full detector simulation can be done using the SLIC software. You will need to convert ProMC files to LCIO files and
 use these files for the "sclic" Geant4-based program. use these files for the "sclic" Geant4-based program.
 +
 +The physics performance studies are  [[fcs:start|listed here]]. They are based on  [[http://atlaswww.hep.anl.gov/asc/jas4pp|Jas4pp program]].
  
 How to convert ProMC to LCIO files is described in How to convert ProMC to LCIO files is described in
 [[community:refhepsim&#converting_to_lcio]] [[community:refhepsim&#converting_to_lcio]]
  
-You can analyse SLCIO files as described in [[https://atlaswww.hep.anl.gov/asc/wikidoc/doku.php?id=fcs:fpad|FPadSoft environment]] explains how to analyse LCIO data after full detector simulations 
 ====== Pileups mixing  ====== ====== Pileups mixing  ======
  
community/refhepsim.txt · Last modified: 2016/04/28 18:23 by asc