Monte Carlo Particle Validation with HepMCAnalysis

How to install and run HepMCAnalysis _i

1) Install in release 17

Log into an lxplus 5 machine and setup a version:

  • export AtlasSetup =/afs/cern.ch/atlas/software/dist/AtlasSetup
  • alias asetup='source $AtlasSetup/scripts/asetup.sh'
  • asetup 17.2.X.Y-VAL,r5,slc5,here
Check out and compile the code
  • mkdir test_hepmc
  • cd test_hepmc
  • cmt co -r GeneratorUtils -00-01-08 Generators/GeneratorUtils
  • use separate branch for release 17, latest tag: cmt co -r HepMCAnalysis _i-00-00-91-02 Generators/HepMCAnalysis_i
  • cd Generators/GeneratorUtils/cmt
  • cmt config & make
  • cd ../../HepMCAnalysis_i/cmt
  • cmt config & cmt make
Install in release 19

setup a version:

  • export AtlasSetup =/afs/cern.ch/atlas/software/dist/AtlasSetup
  • alias asetup='source $AtlasSetup/scripts/asetup.sh'
  • asetup AtlasSimulation,19.2.4,here
Check out and compile the code
  • mkdir test_hepmc
  • cd test_hepmc
  • check out the latest tag: cmt co -r HepMCAnalysis_i-00-01-05 Generators/HepMCAnalysis_i
  • or check out trunk version: cmt co Generators/HepMCAnalysis_i
  • cd HepMCAnalysis_i/cmt
  • cmt config & cmt make
2) Download EVNT files using dq2 tools

Install by executing:

  • source /afs/cern.ch/project/gd/LCG-share/current/etc/profile.d/grid_env.sh
  • source /afs/cern.ch/atlas/offline/external/GRID/ddm/DQ2Clients/setup.sh
If not yet done, create a voms proxy:
  • voms-proxy-init -voms atlas
Get a file list of MC files. As an example, for a W->mu nu process generated by Sherpa, ID=147775, execute:
  • cd ../ (directory: test_hepmc/Generators/HepMCAnalysis_i)
  • mkdir inputFiles
  • cd inputFiles
  • dq2-ls "mc12_8TeV.147775.Sherpa_CT10_Wmunu.evgen.EVNT.e1434*"
This will give a list of data sets where the names start with "mc12_8TeV.147775.Sherpa_CT10_Wmunu.evgen.EVNT.e1434"

mc12_8TeV.147775.Sherpa_CT10_Wmunu.evgen.EVNT.e1434_tid00927769_00
mc12_8TeV.147775.Sherpa_CT10_Wmunu.evgen.EVNT.e1434_tid00927771_00

Download 1 (or more) files by executing:

  • dq2-get -n 1 mc12_8TeV.147775.Sherpa_CT10_Wmunu.evgen.EVNT.e1434/
For the example, if you download 3 files (by selecting "-n 3" in the command above), you get:

/mc12_8TeV.147775.Sherpa_CT10_Wmunu.evgen.EVNT.e1434_tid00927769_00/EVNT.00927769._000441.pool.root.1

/mc12_8TeV.147775.Sherpa_CT10_Wmunu.evgen.EVNT.e1434_tid00927769_00/EVNT.00927769._000878.pool.root.1

/mc12_8TeV.147775.Sherpa_CT10_Wmunu.evgen.EVNT.e1434_tid00927771_00/EVNT.00927771._002604.pool.root.1

3) Run

Define input files by

  • cd ../share (directory: Generators/HepMCAnalysis_i/share)
  • open the file jobOptions.evgen.py and change the line:
svcMgr.EventSelector.InputCollections = [ '*YOUR EVNT FILE 1*', '*YOUR EVNT FILE 2*', '*YOUR EVNT FILE 3*' ]

  • define the output file in the jobOptions.evgen.py file:
myHepMCAnalysis = HepMCAnalysis _i("HepMCAnalysis_i", file = 'EvGen_test.root')

  • use the following setting in the jobOptions.evgen.py file:
myHepMCAnalysis.JetAnalysis=False
myHepMCAnalysis.WplusJetAnalysis=False
myHepMCAnalysis.ZAnalysis=False
myHepMCAnalysis.ZtautauAnalysis=False
myHepMCAnalysis.WtaunuAnalysis=False
myHepMCAnalysis.ttbarAnalysis=False
myHepMCAnalysis.bbbarAnalysis=False
myHepMCAnalysis.UEAnalysis=False
myHepMCAnalysis.EtmissAnalysis=False
myHepMCAnalysis.ElasScatAnalysis=False
#myHepMCAnalysis.UserAnalysis=False

myHepMCAnalysis.ParticleContentAnalysis=True
myHepMCAnalysis.LeptonJetAnalysis=True
myHepMCAnalysis.PdfAnalysis=True

  • then execute: athena jobOptions.evgen.py
This produces the root file to be used for validation.

Scripts to make comparison plots are available here (list will be updated):

LeptonJetAnalysis: https://twiki.cern.ch/twiki/pub/Sandbox/HepMCAnalysis/LeptonJetPlots.cxx

ParticleContentAnalysis: https://twiki.cern.ch/twiki/pub/Sandbox/HepMCAnalysis/ParticleContentPlots.cxx

Validation analyses

1. LeptonJetAnalysis

  • generale multipurpose analysis
  • generator independent
  • arbitrary number of leptons and jets
  • possibility to produce different sets of plots with different kinematics cuts
  • final state of particles
  • stable, charged leptons
  • antikT jet algorithm
2. ParticleContentAnalysis
  • generale multipurpose analysis
  • decay chain tests of the final state hadrons, tau
  • tau polarization variables
  • particle mass/status/pid checks
3. PdfAnalysis
  • event pdf information checks

Links

HepMCAnalysis Tool

DCube integration within RTT framework

RTT main page

GENSER main page

Topic attachments
I Attachment History Action Size Date Who Comment
Unknown file formatcxx LeptonJetPlots.cxx r1 manage 55.0 K 2013-11-27 - 14:11 NataliiaKondrashova  
Unknown file formatcxx ParticleContentPlots.cxx r1 manage 48.1 K 2013-11-27 - 14:11 NataliiaKondrashova  
Edit | Attach | Watch | Print version | History: r46 < r45 < r44 < r43 < r42 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r46 - 2015-11-09 - MariaMorenoLlacer
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2020 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback