Electrons

Code

CMSSW:

in ../crab: mc2pat_386_cfg.py (from dataset page) etc. needs:

~/scratch/top/CMSSW_3_8_6/src% cvs co -r V00-02-01 ElectroWeakAnalysis/WENu
checkdeps -a   (?)
 scram b -j6

Configuration files:

Electron selection: starts from fullLeptonicMuonElectronSelectionMarkus_cff.py (in my public afs area)
~/scratch/top/CMSSW_3_8_6/src/TopAnalysis/TopFilter/python/sequences/fullLeptonElectronElectronSelectionMC_cff.py . NEW! now with PF jets ID corrected
fullLepSynchronisation_cfg.py : for sync files, cmsRun fullLepSynchronisation_cfg.py NEW! , now with vertex cut and trigger match
~/scratch/top/CMSSW_3_8_6/src/MyAnalysis/Configuration/crabnaf/TTEE/fullLepttbar_cfg.py works on MC ttbar samples. NEW! all ready as of 10.2.11
fullLep_data.py under test on Run2010A, needs fullLeptonicElectronElectronSelectionDATA_cff.py NEW! all ready as of 10.2.11

Code added:

For the mass cut:
~/scratch/top/CMSSW_3_8_6/src/TopAnalysis/TopFilter/plugins/DiElectronFilter.cc
~/scratch/top/CMSSW_3_8_6/src/TopAnalysis/TopFilter/plugins/DiElectronFilter.h
~/scratch/top/CMSSW_3_8_6/src/TopAnalysis/TopFilter/python/filters/DiElectronFilter_cfi.py
~/scratch/top/CMSSW_3_8_6/src/TopAnalysis/TopFilter/plugins/SealModule.cc needs to be modified

For the electrons histograms  NEW histograms updated 9.2.11 :
~/scratch/top/CMSSW_3_8_6/src/MyAnalysis/Analyzers/src/ElectronAnalyzer.cc
~/scratch/top/CMSSW_3_8_6/src/MyAnalysis/Analyzers/interface/ElectronAnalyzer.h
~/scratch/top/CMSSW_3_8_6/src/MyAnalysis/Analyzers/python/ElectronAnalyzer_cfi.py
To be fixed: supercluster energy histogram

For the dielectron histograms:
~/scratch/top/CMSSW_3_8_6/src/MyAnalysis/Analyzers/src/DiElectronAnalyzer.cc
~/scratch/top/CMSSW_3_8_6/src/MyAnalysis/Analyzers/interface/DiElectronAnalyzer.h
~/scratch/top/CMSSW_3_8_6/src/MyAnalysis/Analyzers/python/DiElectronAnalyzer_cfi.py

For the electrons ntuple:
~/scratch/top/CMSSW_3_8_6/src/MyAnalysis/Analyzers/src/ElectronNtupleAnalyzer.cc
~/scratch/top/CMSSW_3_8_6/src/MyAnalysis/Analyzers/interface/ElectronNtupleAnalyzer.h
~/scratch/top/CMSSW_3_8_6/src/MyAnalysis/Analyzers/python/ElectronNtupleAnalyzer_cfi.py

For the electron vertex:    NEW added 10.2.11 
~/scratch/top/CMSSW_3_8_6/src/MyAnalysis/Analyzers/src/ElectronVertexAnalyzer.cc
~/scratch/top/CMSSW_3_8_6/src/MyAnalysis/Analyzers/interface/ElectronVertexAnalyzer.h
~/scratch/top/CMSSW_3_8_6/src/MyAnalysis/Analyzers/python/ElectronVertexAnalyzer_cfi.py


For the trigger matching with an object, from Benjamin:
~/scratch/top/CMSSW_3_8_6/src/MyAnalysis/Analyzers/src/TriggerObjectFilter.cc
~/scratch/top/CMSSW_3_8_6/src/MyAnalysis/Analyzers/python/TriggerObjectFilter_cfi.py
Needs in the cfg file:
## produce pat trigger content
process.load("PhysicsTools.PatAlgos.triggerLayer1.triggerProducer_cff")
process.patTriggerEvent.patTriggerMatches = []
## filter trigger (taken trigger for MC on twiki page )
process.load("TopAnalysis.TopFilter.filters.TriggerFilter_cfi") 
process.filterTrigger.TriggerResults = cms.InputTag('TriggerResults','','HLT')
process.filterTrigger.hltPaths  = cms.vstring('HLT_Ele10_SW_L1R')

# match with trigger object with pt>15 GeV
process.load("MyAnalysis.Analyzers.TriggerObjectFilter_cfi")
process.filterTriggerObject.trigger_event=cms.InputTag("patTriggerEvent")
process.filterTriggerObject.hltPath  = cms.string('HLT_Ele10_SW_L1R')
process.filterTriggerObject.ptMin = cms.double(15)

process.p1 = cms.Path(process.generatorTopFilter *
                      # process.generatorZFilter *
                      process.patDefaultSequence *
                      process.filterTrigger *
                      process.patTriggerDefaultSequence *
                      process.filterTriggerObject *
 OK tested 


For the electron vertex cut:  NEW! 
~/scratch/top/CMSSW_3_8_6/src/TopAnalysis/TopFilter/plugins/ElectronVertexFilter.cc
~/scratch/top/CMSSW_3_8_6/src/TopAnalysis/TopFilter/plugins/ElectronVertexFilter.h
~/scratch/top/CMSSW_3_8_6/src/TopAnalysis/TopFilter/python/filters/ElectronVertexFilter_cfi.py


For plotting:
for histograms .L ElectronPlots.C+; .x plotelectrons.C
for tnuples: .L ElectronClass.C+; ElectronClass aa ; aa.Loop()

Sync test Electrons

           step0 filterTrigger  step2 filterOppositeCharge MassVeto TwoHardJets MET OneBtag

ttbar PF     779   760(only EleSw10 used)         291       234         159     129     122        
files:   ttbar_electron_pf_sync_test.root,  Event list       
ttbar PF     779   744(with matching)             291       234         159     129     122        
files:   ttbar_electron_pf_sync_test.root,  Event list       

zjets       17084    14470                 5805              425          11       1       0
zjets       17084    13767(with matching)  5805              425          11       1       0
zjets       17084    13767(with matching)  5768(with vertex) 422          11       1       0
                                           5767              426 (also with beamspot and superCluster.eta cut)
files:    Event list

Still running/todo:

Data top ntuple

DATAA 3.18 pb-1, ready in ../crabnaf/EEDATAA/data_electron_pf_Run2010A_topntuple.root,
DATAB 32.96 pb-1, 69 jobs, ready in ../crabnaf/EEDATAB/data_electron_pf_Run2010B_topntuple.root

             Step1      Step4f     Step5pf
DATAA          1121        2           1
DATAB         10959       20          14       Event list Step4f
DATAB+DATAA   12080       22          15

Data Z ntuple ready in ../crabnaf/eedatainsideZntuple.root
Data top ntuple ready in ../crabnaf/eedatatopntuple.root
TTbar top ntuple 71 jobs ready in ../crabnaf/TT/ttbar_electron_pf_D6T_topntuple.root
TTbar Z ntuple 71 jobs ready in ../crabnaf/TT/ttbar_electron_pf_D6T_insideZntuple.root
DY top ntuple ready in ../crabnaf/DY/dy_electron_pf_D6T_topntuple.root
DY Z ntuple ready in ../crabnaf/DY/dy_electron_pf_D6T_insideZntuple.root

Plots comparison data-MC HERE

Muons
Slides MET: met.pdf , met2.pdf

Muon Ntuples:

DATA Z-selection: ../crabnaf/datainsideZntuple.root
DATAA READY ../crabnaf/DATAA/dataAinsideZntuple.root (Step11 missing)
DATAB1 READY ../crabnaf/DATAB1/dataB1insideZntuple.root
DATAB2 READY ../crabnaf/DATAB2/dataB2insideZntuple.root
DATA top-selection: ../crabnaf/datatopntuple.root
DATAA READY ../crabnaf/DATAA/dataAtopntuple.root -> event list (jets wrongly filled, are not final jets)
DATAB1 READY ../crabnaf/DATAB1/dataB1topntuple.root -> event list (jets wrongly filled)
DATAB2 READY ../crabnaf/DATAB2/dataB2topntuple.root -> event list (jets wrongly filled)
MC Z-selection
DY D6T READY ../crabnaf/DY/dyD6TinsideZntuple.root
TTbar D6T 1306182 events READY ../crabnaf/ttbarD6TinsideZntuple.root
DY Z2
TTbar Z2
MC top-selection
DY D6T READY ../crabnaf/DY/dyD6Ttopntuple.root
TTbar D6T running crab_0_110114_101712 1 job(3%) missing from ../crabnaf/TT/ttbarD6Ttopntuple_temporary.root
DY Z2
TTbar Z2
CODE

======================================================================================================== Sync test Muons

           step0 filterTrigger  step2 filterOppositeCharge MassVeto TwoHardJets MET OneBtag

ttbar PF    727    698                        338               268        197    162   152
files:       fullLepSynchronisation_cfg.py, muon_ttjets_pf_sync.log,  muon_ttjets_pf.root  

ttbarJPT    727    698                        338               268        207    167   152
files:      fullLepSynchronisation_ttbar_maria.py, muon_ttjets_jpt_sync.log,  muon_ttjets_jpt_moreplots.root  

zjets       17291   14010                     6827             554         17      0     0      
files: muon_dy_pf_sync.log, muon_dy_pf.root

==============================================================================================

             Step1         Step4f       Step5f

Run2010B2     11658           22           18                     crabnaf/DATAB2/dataB2new.root
Run2010B1      2161            3            3                     crabnaf/DATAB1/dataB1.root
Run2010A       1278            3            3                     crabnaf/DATAA/dataA.root
Sum           15097           28           24



now with vertex cut
Run2010B2     11625           22           18                     dataB2vertex.root
Run2010B1      2156            3            3                     dataB1vertex.root
Run2010A       1276            3            3                     dataAvertex.root
Sum           15057           28           24                     macros/Nov22datavertex.root

                                                                  plots in plots/topwindow

Z selection on muons and MET


=============================================================================================
Zselection on same samples   
==========

               Step11        Step4MC=Z     Step5MC=2jets   Step6MC=MET     Step7MC=Btag
            (dile.+QCDveto)   


Run2010B2       11658          10538             318            37            18        dataB2insideZ.root
Run2010B1        2161           1940              57             5             1        dataB1insideZ.root
Run2010A         1278           1147              36             1             1        dataAinsideZ.root
-----------------------------------------------------------------------------------
Sum             15097          13625             411            43            20


with vertex and ntuple                      

Run2010B1                                                                               dataB1insideZntuple.root
Run2010A                                                                                dataA_insideZntuple.root

                                                                                      macros/Nov22DatainsideZ.root
                                                                                      plots in plots/Zwindow


TTbarZ2                         9.35             7.13          6.05                   macros/ttbar_muon_pf_Z2_insideZ.root
DYZ2                           13897             347           12.6                   macros/dy_muon_pf_Z2_insideZ.root

Smeared MET


MET smeared tests Z selection (weighted)
============================

             Step11          Step4MC       Step5MC       MET>30(smeared)        BTag

DY 50K                        13882          441                                      crabnaf/DY/dy_muon_pf_test_insideZ.root
TTbar 50K unwe.                     (85)           (64)           (54)          (50)  crabnaf/TT/ttbar_muon_pf_test_insideZ.root


ET smeared tests top selection (weighted)
==============================

DY 50K                          1161        26.4        0 in this small sample        crabnaf/DY/dy_muon_pf_test.root 
TTbar 50K                         33        23.9       19.6                           crabnaf/TT/ttbar_muon_pf_test.root
         

METet>30 sono 172 eventi unweighted at Step 5scp *.gif gallo@zenith.desy.de:~/public/topplots/.
METsmeared2>30 sono 177 eventi, cioe' 5/172 ~ 3% di sistematico allo step 6 (cioe' senza Btag)


                                                                            plots in plots/
Nov22datavertex.root            1463        83         28

Code, NAF documentation

source /afs/desy.de/project/glite/UI/etc/profile.d/grid-env-3.1.sh source /afs/cern.ch/project/gd/LCG-share/current/external/etc/profile.d/grid-env.sh grid-proxy-init -rfc gsissh -X cms.naf.desy.de ini glite32 grid-proxy-init -rfc ini crab ini cmssw cmsenv (nella versione giusta) # dcache: /store/user/egallo srmls "srm://dcache-se-cms.desy.de:8443/srm/managerv2?SFN=/pnfs/desy.de/cms/tier2/store/user/egallo/" # to checkout code fare prima: echo $CVSROOT kinit egallo@CERN.CH # to import code: cmsrel CMSSW_3_8_6 cd CMSSW_3_8_6/src cmsenv addpkg CommonTools/RecoAlgos cvs update -A -r V00-02-15 CommonTools/RecoAlgos/plugins/HBHENoiseFilterResultProducer.cc cvs update -r V00-02-15 CommonTools/RecoAlgos/python/HBHENoiseFilterResultProducer_cfi.py cvs update -r V00-02-15 CommonTools/RecoAlgos/BuildFile cvs co -rV00-03-14-02 RecoEgamma/ElectronIdentification checkdeps -a scram b -j6 cvs co -d TopAnalysis UserCode/Bromo/TopAnalysis cvs update -r 1.15 TopAnalysis/TopAnalyzer/src/JetQuality.cc scram b -j6 # mkdir MyAnalysis commented out not to create a double directory # cd MyAnalysis cvs co -d MyAnalysis/Analyzers UserCode/DirkDammann/Analyzers cvs co -d MyAnalysis/Configuration UserCode/DirkDammann/Configuration scram b -j6 ------------------------------------------------------------ py in: MyAnalysis/Configuration/data/fullLepSynchronisation_cfg.py cmsRun fullLepSynchronisation_cfg.py >& muon_tt_pf.log files in: TopAnalysis/Configuration/python/samples/ZJets_DilepSyncro_cff.py dilepton selection in: TopFilter/python/sequences/fullLeptonicMuonMuonSelection_cff.py dimuon mass filter TopAnalysis/TopFilter/python/filters/DiMuonFilter_cfi.py muon analyzer MyAnalysis/Analyzers/python/MuonAnalyzer_cfi.py cff_py for the different samples in TopAnalysis/Configuration/python/samples/ for instance for TTjets sync TopAnalysis/Configuration/python/samples/TTJets_DilepSyncro_cff.py -------------------------------------------------------------------- Per runnare sync test on 50000 MC events: in MyAnalysis/Configuration/data/ cmsRun fullLepSynchronisation_cfg.py >& muon_tt_pf.log & Risultati: muon_dy_pf_sync.log on DY muon_ttjets_pf_sync.log on ttjets Per runnare sync test sui dati: in MyAnalysis/Configuration/data/ qsub runInbatch.sh l'output e' nella HOME. oppure in MyAnalysis/Configuration/crabnaf con crab sulla NAF (cercare il dataset prima sul DBS di DESY) - FUNZIONA crab cfgfile prima per vedere se funziona hadd output.root directory/*.root funziona -------------------------------------------------------------------- Per runnare crab data2pat sui dati: in: MyAnalysis/Configuration/crab fare prima (?): voms-proxy-init -voms cms sge=glite e io non sono il gruppo dcms crab -create -cfg crabNov22.cfg FUNZIONA con il vecchio data2pat crab -submit ------------------------------------------------------------------------- in macros: .L METPlots.C+ .x plotall.C -> e metterci quello di cui ho bisogno METPlotsSmeared.C confronto con smeared distributions .x count.C per le Ntuples .x metresolutions.C -> per fare i plots delle MET Sigmas .x vertex.C per la dipendenza dal vertice METClass.h, METClass.C per giocare con MET: (../crabnaf/TT/METClassTTbar per il top) .L METClass.C+ METClass aa aa.Loop() .x vertex.C (per riconfrontare ai dati) ../crabnaf/TT/METClassTTbar.C+ per calcolare l'effetto sistematico su TTbar per copiare, dopo klog gallo@DESY.DE: scp ../plots/*.gif gallo@zenith.desy.de:~/public/topplots/. oppure cp ../plots/*.gif /afs/desy.de/user/g/gallo/public/topplots/.