SingleTopFDR2 < Main

TWiki>

Main Web>TWikiUsers>JulienDonini>SingleTopFDR2 (2008-10-28, JulienDonini)

EditAttachPDF

My analysis
Links

My analysis

Results

Summary of single top studies (Single Top meeting 21/10/08): slides

Analysis with ARA 14.2.0

setup

source cmthome_14.2.0/setup.sh -tag=14.2.0

then follow instruction here: AthenaROOTAccess. Do not forget to add all additional 14.2.0 tags in AthenaROOTAccess/cmt/requirement, then do gmake all.

Example of requirement file: requirements.

...

Get AOD with dq2

Get list of datasets from here:

AOD status

DPD status

Then do

   source ~/setup_grid.sh
   cd tmp/jdonini
   dq2_get fdr08_run2.0052304.physics_Muon.merge.AOD.o3_f8_m10 (for example)

an example of grid setup can be found here

   source /afs/cern.ch/project/gd/LCG-share/sl4/etc/profile.d/grid_env.sh
   source /afs/cern.ch/atlas/offline/external/GRID/ddm/DQ2Clients/setup.sh
   source /afs/usatlas.bnl.gov/Grid/Don-Quijote/dq2_user_client/setup.zsh.CERN
   voms-proxy-init -voms atlas

do not forget to put your certificates keys in .globus.

AOD analysis with CBNT_AAN

We give in this section simple tools to produce a root file from AOD's (or DxPD, x=1,2), for single top analysis. This example was largely inspired from here: WorkBookAODAnalysisRel14.

First install code:

Install PhysicsAnalysis/AnalysisCommon/AnalysisExamples package following these instructions.

Get tarball for single top analysis here, then untar it in PhysicsAnalysis/AnalysisCommon/AnalysisExamples directory.

The code used for producing a simple single top root file from AOD is listed below:

   cmt/requirements
   python/SingleTopFlags.py
   AnalysisExamples/DoElectronSingleTop.h
   AnalysisExamples/DoMuonSingleTop.h
   AnalysisExamples/DoParticleJetSingleTop.h
   AnalysisExamples/DoBJetSingleTop.h
   AnalysisExamples/singletop.h   # --> main single top code
   src/DoElectronSingleTop.cxx
   src/DoMuonSingleTop.cxx
   src/DoBJetSingleTop.cxx
   src/DoParticleJetSingleTop.cxx
   src/singletop.cxx   # --> main single top code
   share/DoMuonSingleTop_jobOptions.py
   share/DoElectronSingleTop_jobOptions.py
   share/DoBJetSingleTop_jobOptions.py
   share/DoParticleJetSingleTop_jobOptions.py
   run/SingleTopExample_jobOptions.py   # main jobOption
   run/singletop_fdr08_run2.0052290.physics_Muon.merge.AOD.o3_f8_m10._0007.1.root    # --> an example of output root file

go to cmt/ and compile (gmake)

go to run/ and produce a root file out of the AOD:

   athena.py SingleTopExample_jobOptions.py >& log_singletopexample &

an example of output is given here

   run/singletop_fdr08_run2.0052290.physics_Muon.merge.AOD.o3_f8_m10._0007.1.root

Luminosity calculation with LumiCalc

See more info on LumiCalc here.

Examples (tested for 14.2.0):

   get_files LumiCalc.py
   LumiCalc.py --readoracle --mc --r=52280 --trigger=EF_mu40 --debug
   LumiCalc.py --readoracle --debug --mc --trigger=EF_mu40 /tmp/jdonini/fdr08_run2.0052304.physics_Muon.merge.AOD.o3_f8_m10._0001.1

For FDR2 some lumiblocks are missing. See FDR forum here.

Lum for FDR2 runs (nb-1) in Marcello's D3PD:

Trigger	EF_mu40		EF_e20i
Run	Egamma	Muon	Egamma	Muon
52280	345	165	345
52283	349	179	349
52290	354	108	354
52293	317	107	317
52301		161	301
52304		165	345

Event mixing bug

See useful info here. In brief for single top t-channel process, the production cross section is raised from 80 pb to about 1747 pb for 30-seconds luminosity blocks (runs 4, 5) and to about 913 pb for 1-minutes blocks (runs 1, 2).

Links

Topic attachments
I	Attachment	History	Action	Size	Date	Who	Comment
tgz	AnalysisExamples.tgz	r1	manage	6842.3 K	2008-07-10 - 16:34	JulienDonini
ext	requirements	r1	manage	2.2 K	2008-07-04 - 12:23	JulienDonini

Topic revision: r8 - 2008-10-28 - JulienDonini

Main

Webs

Public webs

Welcome Guest

- Cern Search
- TWiki Search
- Google Search
Main All webs

Copyright &© 2008-2024 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