semileptonic ttbar studies with electrons

analysis setup

recipe for dummies

• copied from the 371 official recipe

cmsrel CMSSW_3_7_1
cd CMSSW_3_7_1/src
cmsenv
cvs co -r  Colin_Sept14_histos -d   AnalysisDataFormats/PFAnalyses UserCode/PFAnalyses/AnalysisDataFormats/PFAnalyses      
cvs co -r  mnguyen_Sept29 -d   PFAnalyses/CommonTools   UserCode/PFAnalyses/CommonTools   
cvs co -r  Colin_Sept28_newSamplesInterface -d   PFAnalyses/RootTools  UserCode/PFAnalyses/RootTools                          
cvs co -r  Artur_Sept28_prod  -d   PFAnalyses/TemplateAnalysis  UserCode/PFAnalyses/TemplateAnalysis  
cvs co -r Michal_July6_CMSSW_3_7_0_patch4 PhysicsTools/PatAlgos
cvs co -r V04-06-01 PhysicsTools/PFCandProducer
cvs up -r 1.46.2.2 PhysicsTools/PatAlgos/python/tools/pfTools.py 
cvs up -r 1.57.4.1 PhysicsTools/PatAlgos/python/tools/jetTools.py
cvs co -r V00-02-16 CommonTools/Utils
cvs co -r V00-02-17 PhysicsTools/SelectorUtils           
scram b -j 4            

cvs co -d PFAnalyses/PGTTbarSemiEle UserCode/Govoni/PGTTbarSemiEle
scram b -j 4
cp PFAnalyses/PGTTbarSemiEle/python/otherPackages/electronProducer_cff.py PhysicsTools/PatAlgos/python/producersLayer1/electronProducer_cff.py

files description

• DumpEventID.py list the event and run number for each event in a root file
• pickup_events.py extract a subset of events from a root file
• stdPat_creator_selected.py produce std PAT-tuples from RECO
• patTuple_PF2PAT_cfg.py produce std and PFlow PAT-tuples from RECO
• PATTrivialPlots_both.py plot basic quantities for both std PAT objects and PFlow PAT objects
• PATTrivialTree_both.py put basic quantities into TNtuples for each type of physics objects, for both std PAT objects and PFlow PAT objects
• PATTrivialPlots_std.py plot basic quantities for std PAT objects

pat-tuples production

• edit ../../../PhysicsTools/PatAlgos/python/producersLayer1/electronProducer_cff.py to add full electron ID and ISO

useful things

links

• PFlow analyses twiki
• PFlow ttbar analysis twiki
• useful instructions within CMSSW
• PF2PAT twiki
• FWlite beginners twiki
• link to myt personal repository
• CMG group meetings indico page
• PhysicsTools/PATExamples link to the code

extract the first N events from a data file

• see also here

import FWCore.ParameterSet.Config as cms
process = cms.Process("PICKEVENTS")
process.source = cms.Source("PoolSource",
    fileNames = cms.untracked.vstring( 'file:/tmp/govoni/6E7C4631-9D47-DF11-96CE-003048C69288.root' )
#    eventsToProcess = cms.untracked.VEventRange( '1:909002-1:909767' ) # THIS LINE SHOULD WORK
)
process.Out = cms.OutputModule("PoolOutputModule",
    outputCommands = cms.untracked.vstring('keep *'),
    fileName = cms.untracked.string('requested_events.root')
)
process.e = cms.EndPath(process.Out)
process.maxEvents = cms.untracked.PSet(
    input = cms.untracked.int32(100)
)

-- PietroGovoni - 30-Sep-2010