Higgs To Four Leptons Analysis Example With Release 10.0.1

Purpose to give a working example for release 10.0.1 in order to analyse the already
existing Rome Initial Layout samples. Based on the Tutorial http://www.usatlas.bnl.gov/PAT/tutorial904.html

Make sure that CMT is setup correctly

1. Check that the requirements file is set for 10.0.1

########################Atlas requirements file #######################
# Setup Environment
#
#
# Several alternate settings for ATLAS-AREA have to be considered according
# to the CMTSITE
#
# The macros ATLAS_DIST_AREA and ATLAS_AREA are mandatory.
#
set CMTSITE CERN
#set SITEROOT /afs/cern.ch

macro ATLAS_DIST_AREA "/afs/cern.ch/atlas/software/dist"

#macro ATLAS_DIST_AREA "${SITEROOT}/atlas/software/dist"

# Select releases as they are available

macro ATLAS_RELEASE "10.0.1"

use AtlasSettings v* $(ATLAS_DIST_AREA)

# Prepend your test area(s) according to your choices

path_remove CMTPATH /afs/cern.ch/user/l/lagouri/MyTest/$(ATLAS_RELEASE)
path_prepend CMTPATH /afs/cern.ch/user/l/lagouri/MyTest/$(ATLAS_RELEASE)

############################################

2. Source

  • > source /afs/cern.ch/sw/contrib/CMT/v1r18p20050501/mgr/setup.sh
  • > cmt config

3. Create the working directory

  • > mkdir MyTest/10.0.1

4. Setup CMT

  • > source setup.sh -tag=opt

5. Check that CMT path is set correctly

  • > echo $CMTPATH | tr ":" "\n"
  • Something like that is printed to the screen:
/afs/cern.ch/user/l/lagouri/MyTest/10.0.1
/afs/cern.ch/atlas/software/dist/10.0.1
/afs/cern.ch/atlas/offline/external/Gaudi/0.14.6.14-pool201
/afs/cern.ch/atlas/offline/external/LCGCMT/LCGCMT_33

  • The path should at least contain your working directory (MyTest), the ATLAS software distribution (dist, Gaudi), and the external libraries (LCGCMT) .

5. Change Directory to your working area for this tutorial

  • > cd 10.0.1

AOD Analysis

Check out the UserAnalysis Package

  • > cd ~/MyTest/10.0.1

  • > cmt co -r UserAnalysis-00-02-02 PhysicsAnalysis/AnalysisCommon/UserAnalysis

Compile and build the UserAnalysis Package

  • > cd PhysicsAnalysis/AnalysisCommon/UserAnalysis/UserAnalysis-00-02-02/cmt/

  • > cmt config
  • > source setup.sh
  • > cmt broadcast gmake

Create your own analysis algorithm

  • A package, UserAnalysis, is provided to help you get started with the development of your analysis codes. This package consists of a simple Gaudi Algorithm, the AnalysisSkeleton. The user may check out this package and start implementing his/her analysis code.

  • > cd ../run
  • download AnalysisSkeleton_jobOptions.py and rename it to AnalysisSkeleton_jobOptions.py
  • pool_insertFileToCatalog "rfio:/castor/cern.ch/user/k/ketevi/10.0.1/zee/aod/rome.initial.zee.esd1001.aod1001.0001.pool.root"
  • > athena.py -b AnalysisSkeleton_jobOptions.py

  • The source code for this example Z->2 ee is found in /MyTest/10.0.1/PhysicsAnalysis/AnalysisCommon/UserAnalysis/UserAnalysis-00-02-02/src/AnalysisSkeleton.cxx

  • It should be successful, producing a ROOT output file called AnalysisSkeleton.root. Open this file in ROOT and browse the histograms.

  • The data file data you are using is Z -> ee, (/castor/cern.ch/user/k/ketevi/10.0.1/zee/aod/rome.initial.zee.esd1001.aod1001.0001.pool.root), simulated with the Rome initial layout. You will retrieve a container AOD electrons from the persistency storage. You will loop over the electrons in the container, and for each electron, you will plot the electron transverse momentum and pseudo-rapidity and the E/P. You can calculate the ee pair invariant mass and put it in a histogram.

Analysis Examples

How to run the example Higgs to four leptons

Check out the AnalysisExample Package

  • > cd ~/MyTest/10.0.1

  • > cmt co -r AnalysisExamples-00-09-18 PhysicsAnalysis/AnalysisCommon/AnalysisExamples

  • You should have the following .cxx files in your PhysicsAnalysis/AnalysisCommon/AnalysisExamples/AnalysisExamples-00-09-18/src/
    • DoBJet.cxx
    • DoMuon.cxx
    • DoParticleJet.cxx
    • DoTauJet.cxx
    • H4l.cxx
    • MiscellaneousExamples.cxx
    • TrackExample.cxx
    • ZeeOnESD.cxx
    • Ztautau.cxx
    • DoElectron.cxx
    • DoParticle.cxx
    • DoPhoton.cxx
    • FilterExample.cxx
    • Hbb.cxx
    • ReadTruthParticles.cxx
    • ttbar.cxx
    • ZeeZmmOnAOD.cxx

Compile and build the AnalysisExample Package

  • I run this example with 2524 events
  • Input files: Roma Initial Layout H(130)->4l (e,mu)
(/castor/cern.ch/grid/atlas/datafiles/rome/recov10/rome.003045.recov10.H1_130_4l/*.root)

  • Electron selection cuts are controlled by DoElectron_jobOptions.py :

# run the DoElectron algorithm
AOD_ANA.members += [ "DoElectron" ]

# electron analysis parameters
DoElectron = Algorithm( "DoElectron" )
DoElectron.DeltaRMatchCut = 0.2
DoElectron.MaxDeltaR = 0.9999
DoElectron.AODContainerName = "ElectronCollection"
DoElectron.TruthContainerName = "SpclMC"
DoElectron.UserContainerName = "MyPreSelectedElectrons"
# Electron Pt > Pt cut
DoElectron.ElectronEtCut = 6.0*GeV
DoElectron.ElectronEtaCut = 2.5
DoElectron.ElectronCone = 0.9
DoElectron.NtupleLocID = "/FILE1/Electron/t100"
DoElectron.NtuplePrefix = "Electron"
DoElectron.HistDirectoryName = "Electron"

  • Muon selection cuts are controlled by DoMuon_jobOptions.py

# run the DoMuon algorithm
AOD_ANA.members += [ "DoMuon" ]

# muon analysis parameter - muons from MuID
DoMuon = Algorithm( "DoMuon" )
DoMuon.DeltaRMatchCut = 0.2
DoMuon.MaxDeltaR = 0.9999
DoMuon.AODContainerName = "MuonCollection"
DoMuon.TruthContainerName = "SpclMC"
DoMuon.UserContainerName = "MyPreSelectedMuons"
DoMuon.MuonCone = 0.9
# Muon Pt > Pt cut
DoMuon.MuonEtCut = 6.0*GeV
DoMuon.MuonEtaCut = 2.7
# Chi2 > Chi2 cut
DoMuon.MuonChi2Cut = 20.0
DoMuon.NtupleLocID = "/FILE1/Muon/t300"
DoMuon.NtuplePrefix = "Muon"
DoMuon.HistDirectoryName = "Muon"

  • H4l cuts are controlled by H4lExample_jobOptions.py :

# particle pre selection job options
if AnalysisFlags.DoElectron:
include( "./DoElectron_jobOptions.py" )
if AnalysisFlags.DoMuon:
include( "./DoMuon_jobOptions.py" )

# User analysis steering algorithm
AOD_ANA.members += [ "H4l" ]
H4l = Algorithm( "H4l" )
H4l.PreSelectedElectrons = "MyPreSelectedElectrons"
H4l.PreSelectedMuons = "MyPreSelectedMuons"
H4l.doH4e = TRUE
H4l.doH4m = TRUE
H4l.doH2e2m = TRUE
# m12-Mz < deltaMz, 12 leading leptons
H4l.deltaMz = 15.0*GeV
# m34 > minMz, 34 other leptons
H4l.minMz = 20.0*GeV
# 2 leading leptons Pt > Pt cut
H4l.LeadingLeptonsPtCut = 20.0*GeV
# 2 other leptons Pt > Pt cut
H4l.OtherLeptonsPtCut = 7.0*GeV
H4l.isBothZOnShell = FALSE
H4l.ZMassResolution = 2.8*GeV
H4l.NtupleLocID = "/FILE1/H4l/t800"
H4l.NtuplePrefix = "H4l"
H4l.OutputLevel = ERROR

  • When it run successful, it produces the root files: H4l.ntuple.root and H4l.hist.root.
  • Open these files in ROOT and browse the histograms.

Useful Information:

  • Higgs Roma Data Samples in castor area
    • /castor/cern.ch/grid/atlas/datafiles/rome/recov10/
      • 1. H(130)
        • /castor/cern.ch/grid/atlas/datafiles/rome/recov10/rome.003045.recov10.H1_130_4l
      • 2. H(150)
        • /castor/cern.ch/grid/atlas/datafiles/rome/recov10/ rome.004310.recov10.H10_150_4l
      • 3. H(180)
        • /castor/cern.ch/grid/atlas/datafiles/rome/recov10/rome.003046.recov10.H2_180_4l
  • Higgs Background Roma Data Samples in castor area
    • /castor/cern.ch/grid/atlas/datafiles/rome/recov10/
      • 1. Zbbar
        • /castor/cern.ch/grid/atlas/datafiles/rome/recov10/rome.004321.recov10.AcerMC_Zbb_4l
      • 2. ttbar
        • /castor/cern.ch/grid/atlas/datafiles/rome/recov10/rome.004322.recov10.MCatNLO_tt_4l

Useful links:

-- TheodotaLagouri - 28 Oct 2006

Topic attachments
I Attachment History Action Size Date Who Comment
Texttxt AnalysisSkeleton_jobOptions.py.txt r1 manage 1.9 K 2006-10-28 - 18:34 TheodotaLagouri AnalysisSkeleton.py
Texttxt DoElectron_jobOptions.py.txt r1 manage 0.7 K 2006-10-29 - 12:09 TheodotaLagouri DoElectron_jobOptions.py
Texttxt DoMuon_jobOptions.py.txt r1 manage 0.6 K 2006-10-29 - 12:10 TheodotaLagouri DoMuon_jobOptions.py
Postscriptps H2e2m.ps r1 manage 14.4 K 2006-10-29 - 12:17 TheodotaLagouri H2e2m.ps
Postscriptps H4e.ps r1 manage 14.9 K 2006-10-29 - 12:17 TheodotaLagouri H4e.ps
Texttxt H4lExample_jobOptions.py.txt r1 manage 6.2 K 2006-10-29 - 11:13 TheodotaLagouri H4lExample_jobOptions.py
Postscriptps H4m.ps r1 manage 14.5 K 2006-10-29 - 12:16 TheodotaLagouri H4m.ps
XMLxml PoolFileCatalog.xml r1 manage 8.2 K 2006-10-29 - 12:11 TheodotaLagouri PoolFileCatalog.xml
Edit | Attach | Watch | Print version | History: r6 < r5 < r4 < r3 < r2 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r6 - 2006-11-24 - TheodotaLagouri
 
    • 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