SLArSimulationPage

Introduction

This page is aimed to provide information about sLAr upgrade simulation work in Tokyo group. Plan in Tokyo group

Presentation

Date Meeting Slides
2012.07.13 LAr Upgrade Simulation sLAr Simulation Plan in Tokyo
2012.08.28 LAr Simulation L1 Jet Trigger Study for LAr Upgrade
2012.09.24 LAr sLHC Simulation Update of L1 Jet Trigger Study
2012.09.27 --- L1 Jet Trigger Study for High-Luminosity Upgrade
2012.12.04 LAr sLHC Simulation L1 Jet Study for sLAr Upgrade
2012.12.04 --- Summary of J17 Sample Investigation
2013.03.13 LAr sLHC Simulation Summary of Present Status and Plans for Jet Trigger
2013.05.15 LAr sLHC Simulation Gaussian Jets for L1 Jet Trigger
2013.07.04 LAr sLHC Simulation Trigger Rate Study for L1 Jet
2013.09.02 LAr Weekly Meeting Gaussian Filter Jets for L1 Jet Trigger
2013.11.05 LAr sLHC Simulation Jet Trigger Results for Phase I upgrade

Document

Outline of Jet Trigger Performance Chapter in TDR

The TeX files for the TDR document can be found in SVN area.

Introduction

  • Motivation
  • Problem of present L1 jet
  • Goal

Analysis framework

  • Basics about input to L1 jet (e.g, trigger tower, digitization)
  • Description about Gaussian filtering, sliding window and anti-kt algorithms

Inclusive jet trigger performance

  • Performance evaluation for different energy digitization using QCD dijet samples (fixed mu=60 or 80?)
    • Single jet trigger turn-on (e.g, J20) with energy digitization of 1024 MeV and a fixed mu for sliding window (SW), anti-kT R=0.4 (AKt) and Gaussian filter (GF) sigma=0.4 (or 0.4/sqrt(2)) jets
    • Single jet rigger turn-on resolution (fixed mu) as a function of energy digitization for SW, AKt and GF jets
  • Performance for different pileup conditions (mu=20-80, maybe up to 140-200?) with fixed energy digitization of 250 MeV
    • Single jet trigger turn-on (fixed mu) for SW, AKt and GF jets
    • Single jet trigger turn-on resolution as a function of mu for SW, AKt and GF jets
    • Table of single jet trigger rates vs mu for SW and GF jets

Multijet trigger performance

  • Performance for multijet triggers (fixed mu=60 or 80?) with fixed energy digitization of 250 MeV
    • Multijet trigger (e.g, 4J30) turn-on as a function of 4th offline jet pT for SW, AKt and GF sigma=0.1 and 0.4/sqrt(2) jets (pT_jet1,2,3 > 100 GeV)
    • Modification to GF algorithm : GF(sigma=0.1) for seeding and Gaussian weighted pT sum for energy measurement
    • Multijet trigger turn-on as a function of 4th offline jet pT for SW and modified GF jets (pT_jet1,2,3 > 100 GeV)
    • Table of multijet trigger rates vs mu for SW and modified GF jets (not pT_jet1,2,3 requirement)

Forward jet performance

Performance for new physics signature (if possible)

  • Trigger performance for e.g, full-hadronic ttbar, SUSY (squark, gluino), Z' to ttbar events:
    • Single jet trigger turn-on (fixed mu) for SW and GF jets
    • Multijet trigger turn-on (fixed mu) as a function of 4th offline jet pT for SW and GF jets (pT_jet1,2,3 > 100 GeV)

L1 Jet Trigger

D3PD analysis

Simulation framework using D3PD by J. Philipp Grohs:

Analysis code from J. Philipp Grohs located at /data/maxi39/atljphys/terashi/sLArSimulation/D3PD SVN

  • Info about the code : 1, 2
    • cfg: xml configuration scripts, you can steer the simulation with these files
      • Granularity of super cell, trigger towers can be changed if you want by editing xml file.
    • doc: code documentation, you have to execute the doxygen script first, most important classes are documented reasonably well
    • exe_simulation: executes the simulations from lib_analyses
    • GNUmakefile: global make file, build all libraries and executables
    • initProject.sh: requirements for the code compilation
    • job: bsub scripts, also for je dumping, look at this files to learn how to start the executable
    • lib_analyses: actual analysis code, MET calculation, je dumper, ...
    • lib_analysis_utilities: basic analysis utilities
    • lib_calo_sim: actual simulation code
  • Installation tips:(condition: gcc.4.3.5, ROOT 5.30/02)
    • Requires Boost C++, It is obtained locally in /data/data9/zp/yamanaka/opt/boost_1_50_0
      • Set environment variables as follows:
        export BOOST_INC=/data/data9/zp/yamanaka/opt/boost_1_50_0
        export BOOST_LIB=/data/data9/zp/yamanaka/opt/boost_1_50_0/lib

      • Change the following lines (Maybe around line 103) in lib_analysis_utilities/src/Declarations.h
        #if !(defined __GNUG__) || (defined __OPEN64__)
        # define MY_BOOST_SCOPE_EXIT BOOST_SCOPE_EXIT
        #else
        # define MY_BOOST_SCOPE_EXIT BOOST_SCOPE_EXIT_TPL
        #endif

        to the following
        # define MY_BOOST_SCOPE_EXIT BOOST_SCOPE_EXIT
      • Change the line 286 of lib_analyses/src/ETMissSim.cpp to the following:
        CUpgradeSimBase<CEGammaReader>::AfterLoop();
      • Here is the fixed version CaloSim-00-01-01
  • JetElementDumper with reconstructed jets CaloSim-00-01-05 (-j option)
    • Usage: bin/exe_simulation.exe -a JetElementDumper -o output.root -i egamma.root -p -t -m -w -j -s cfg/DefaultThreshold.xml

Procedure

  • Generate D3PD with calorimeter cell using D3PDMaker from existing ESD
    • Branch: cc_*
    • Local samples: /data/maxi39/atljphys/terashi/sLArSimulation/D3PD/Samples/user.jgrohs.mc11_14TeV.106043.PythiaWenu_no_filter.recon.D3PD.e940_s1321_d602_r2887.Jul12.10h/user.jgrohs.002830.StreamNTUP_EGAMMA._00*.root
    • Event size -> too large to store all events in this format?
      • 106043.PythiaWenu mu=69? xsec=18774 pb @14TeV: 36.6MB/event
      • 105802.JF17_pythia_jetjet mu=80 xsec=3.5272E+09 pb @14TeV, filt. eff.=0.0709: 36.5MB/event

  • Calorimeter simulation: create super cells, trigger tower from calorimeter cell information
    • Example command to generate trigger tower D3PD from calorimeter cell D3PD
      bin/exe_simulation.exe -a TriggerTowerDumper -o output.root -i egamma.root -p -t -s cfg/DefaultThreshold.xml
    • Event size:
      • 106043.PythiaWenu mu=69? xsec=18774 pb @14TeV: 73.9 kB/event
      • 105802.JF17_pythia_jetjet mu=80 xsec=3.5272E+09 pb @14TeV, filt. eff.=0.0709: 72.6 kB/event

Sample

  • mc11_14TeV.105802.JF17_pythia_jetjet_filter
    • Samples in local directory: /data/maxi42/atljphys/yamanaka/mc11_14TeV/NTUP_JE
    • Various mu (=0, 12.5, 17.5, 27.5, 41, 55, 80) samples are available and it is generated in 14 TeV.
      μ reco tag
      0 r2962
      12.5 r3088
      27.5 r3089
      41 r3090
      55 r3091
      80 r3092
  • mc11_7TeV.1050xx.Jx_pythia_jetjet wiki
    • Various mu (=0, 10, 20, 30, 40) samples are available although they are lower than 14 TeV samples.
    • ECM=7 TeV samples are only available but it will be enough for jet studying
    • Sample in local directory: /data/maxi42/atljphys/yamanaka/mc11_7TeV/NTUP_JE
  • mc12_8TeV samples?
  • Signal sample
    • mSUGRA with large m0
      • Current limit on mSUGRA with tanβ=10, A0=0, μ>0 at 8 TeV, 5.8 fb-1 by ATLAS CONF.
      • At large m0, sensitivity is almost constant respect to m1/2
    • stop or top pair production, hadronic decay

Athena super-cell

Denis's code to simulate super-cell in athena
  1. Get RDO
  2. Get the right nigtly: e.g. asetup dev,rel_1,here,setup
  3. Get the jobOptions from ~damazio/public/TestNtuple4.py (or ~damazio/public/TestNtuple3.py),
    then edit the line : PoolRDOFile='myRDO_withAddCaloDigi.pool.root'
  4. Run the jobOptions;
    athena TestNtuple3.py

Super cell simulation

D3PD

Jet reconstruction

Gaussian filter jet

Gaussian filtering code (based on SCons) by B. Cole and Y. S. Lai located at /data/maxi39/atljphys/terashi/sLArSimulation/GaussianFilter
  • Installation tips:
    • Get SCons and install to /usr/bin or set PATH to the installed directory.
    • Get ccache and install if not installed as default.
    • In gaussian-filter-jet directory, execute scons, then the executable of example command is made in build/src/example/randomjet.
  • External use:
    • Source files need to be included in compile.
      • -I (gaussian-filter-jet)/src/libjetrec
      • -I (gaussian-filter-jet)/src/libjetbase
      • -I (gaussian-filter-jet)/src/libjetapprox
    • Three libraries are generated. They need to be linked in compile.
      • -L (gaussian-filter-jet)/build/src/libjetrec -ljetrec
      • -L (gaussian-filter-jet)/build/src/libjetbase -ljetbase
      • -L (gaussian-filter-jet)/build/src/libjetapprox -jetapprox
    • In compile, -fopenmp option is required.
    • Other compile options are required? -Wall -fPIC -Wno-unused -Werror -O2 -ftree-vectorize -msse -msse2 -msse3 -mssse3 -mfpmath=sse -DHAVE_SSE -DHAVE_SSE2 -DHAVE_SSE3 -DHAVE_SSSE3 -D_LARGEFILE64_SOURCE
    • Architecture of main source code and gaussian-filter-jet libraries need to be matched.
      • Architecture of scons libraries can be changed by replacing the following in the SConscript
        env.Replace(ARCH = ARGUMENTS.get('arch', os.uname()[4]))
        by
        env.Replace(ARCH = 'something')
        'something' will be 'x86_64', 'i386' and so on.

Detail of the code

  • reconstruction_filter_iir_t() is used, which is defined in libjetrec/jetrec/rec.h
  • Functions are defined in recfilter.c
  • jet reconstruction is performed by reconstruct()
    • In the function, if refine is true (default), refine_jet_maximum() and ghost_cut() are performed.
    • refine_jet_maximum() is defined in libjetrec/recfilteroptimize.cc as a member function of reconstruction_filter_t()
    • ghost_cut() is defined in libjetrec/recfilter.cc

More info about the filtering algorithm found at

Fastjet

To compare with the current jet algorithm used in ATLAS (anti-kt), fastjet is also required.
  • After installed, example command is generated under example directory.
  • Example of anti-kt algorithm is
    ./01-basic < data/single-event.dat
  • External use:
    • Include header files in compile
      • -I(fastjet)/include
      • -I(fastjet)/tools
    • Link library (since there are several version of libraries, they will conflict each other unless one of them is specified directory)
      • (fastjet)/src/.libs/libfastjet.a

L1 jet algorithm

L1 jet algorithm currently used in ATLAS is also needed for comparison, TrigT1CaloUtils/JetAlgorithm.cxx .

L1 jet trigger study

Condition

  • Use central jets (|eta|<3.2) only.
  • As offline jets, anti-kt dR=0.4 TopoCluster jets in EM scale is used.

Study items

  • Jet resolution expected to be similar between different jet algorithms if input is the same.
  • While jet separation is expected to be better for Gaussian filter jets.

Reference

-- KojiTerashi - 24-Jul-2012

Topic attachments
I Attachment History Action Size Date Who Comment
Texttxt Grohs1.txt r1 manage 2.4 K 2012-07-24 - 11:29 KojiTerashi  
Texttxt Grohs2.txt r1 manage 5.7 K 2012-07-24 - 11:30 KojiTerashi  
PDFpdf L1JetStudySummary_120927.pdf r2 r1 manage 547.5 K 2012-10-08 - 22:08 TakashiYamanaka  
PDFpdf SummaryOfJJ17_121204.pdf r1 manage 400.1 K 2012-12-07 - 20:27 TakashiYamanaka  
JPEGjpg http___iktp.tu-dresden.jpg r1 manage 73.0 K 2012-07-24 - 23:03 TakashiYamanaka  
Edit | Attach | Watch | Print version | History: r73 < r72 < r71 < r70 < r69 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r73 - 2014-05-20 - TakashiYamanaka
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

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