Pattuple production for A->Zh an other

CMSSW version

CMSSW version 5_3_11

showtags


Test Release based on: CMSSW_5_3_11
Base Release in: /cvmfs/cms.cern.ch/slc5_amd64_gcc462/cms/cmssw/CMSSW_5_3_11
Your Test release in: /lustre/cmswork/lacaprar/Ana/CMSSW_5_3_11
--- Tag ---    --- RelTag --- -------- Package --------                        
NoCVS                         AZh/Analysis                                     
NoCVS                         AZh/Analysis_old                                 
NoCVS                         AZh/Limits                                       
NoCVS                         AZh/Ntuple                                       
NoCVS                         AZh/PDFUncert                                    
NoCVS                         AZh/Pattupler                                    
NoCVS                         AZh/Test                                         
NoCVS                         AZh/lhapdf-5.9.1                                 
V00-03-04                     CMGTools/External                                
V00-01-04      V00-01-04      CommonTools/RecoUtils                            
V02-05-11      V02-05-12      DataFormats/CaloRecHit                           
V03-03-11-01   V03-03-11      DataFormats/METReco                              
V15-03-04-01   V15-03-04-01   DataFormats/ParticleFlowCandidate                
V06-05-06-12   V06-05-06-11   DataFormats/PatCandidates                        
V00-02-14      V00-02-13-01   DataFormats/StdDictionaries                      
V10-02-02-01   V10-02-02-01   DataFormats/TrackReco                            
V02-00-04-01   V02-00-04-01   DataFormats/VertexReco                           
V00-00-08                     EgammaAnalysis/ElectronTools                     
CMSSW_5_3_11   V00-04-01      ElectroWeakAnalysis/Utilities                    
V00-00-70      V00-00-70      FWCore/GuiBrowsers                               
V05-01-09      V05-01-06      JetMETCorrections/Modules                        
V04-06-09-02   V04-06-09-01   JetMETCorrections/Type1MET                       
muscle_v4_2_0                 MuScleFit/Calibration                            
V08-09-62      V08-09-59      PhysicsTools/PatAlgos                            
V03-09-28      V03-09-28      PhysicsTools/PatUtils                            
V00-07-05      V00-06-02      RecoBTag/Configuration                           
V01-04-09-01   V01-04-08      RecoBTag/ImpactParameter                         
V01-03-21      V01-03-19      RecoBTag/PerformanceDB                           
V01-10-06      V01-08-01      RecoBTag/SecondaryVertex                         
V05-09-11      V05-08-01      RecoBTag/SoftLepton                              
V02-03-02      V02-02-07      RecoBTau/JetTagComputer                          
V09-00-01      V08-11-22      RecoEgamma/EgammaTools                           
V00-05-cmsdas2013-02 V00-07-02-03   RecoJets/JetAnalyzers                            
NoCVS          V04-02-01      RecoLuminosity/LumiDB                            
V04-02-07                     RecoLuminosity/LumiDB_old                        
V00-00-08                     RecoMET/METAnalyzers                             
V00-00-13-01   V00-00-13-01   RecoMET/METFilters                               
V03-03-12-02   V03-03-12-02   RecoMET/METProducers                             
V15-02-06      V15-02-05-01   RecoParticleFlow/PFProducer                      
V01-04-14      V01-03-04      RecoTauTag/Configuration                         
V01-04-26      V01-04-12-00   RecoTauTag/RecoTau                               
---------------------------------------
total packages: 41 (41 displayed)

Hand made changes

MuSCLE Fit requires this cvs co -r muscle_v4_2_0 -d MuScleFit/Calibration UserCode/scasasso/MuScleFit/Calibration

cd EgammaAnalysis/ElectronTools/data/
cat download.url | xargs wget 

EgammaAnalysis/ElectronTools/plugins/CalibratedPatElectronProducer.cc l82

case 1: pathToDataCorr = "./scalesMoriond.csv";
was this incredible piece of shit. In a public release! Can you imagine it?
case 1: pathToDataCorr = "../data/scalesMoriond.csv";
the csv file is in EgammaAnalysis/ElectronTools/data/scalesMoriond.csv No way it would works if you run it anywhere else than EgammaAnalysis/ElectronTools/test Unbelievable!

MuSCLE Fit requires this cvs co -r muscle_v4_2_0 -d MuScleFit/Calibration UserCode/scasasso/MuScleFit/Calibration

Information sources

Configuration file

https://twiki.cern.ch/twiki/pub/Sandbox/StefanoLacapraraPat2012/patTuple_PF2PAT_cfg.py.txt

I'm not 100% sure about the top-projection stuff, not sure what it does, what it is supposed to do (remove muon and electron and taus recognized as such from the jets collection?) or if it works or not. Or if we need it, in the first place. It was there in hbb analysis, will keep it for the time being.

What works so far:

Electrons
    • Ele calibration, regression and smearing (for MC)
    • Electron ID:
ID: eidLoose 15
ID: eidRobustHighEnergy 4
ID: eidRobustLoose 7
ID: eidRobustTight 4
ID: eidTight 15
ID: mvaNonTrigV0 0.926894
ID: mvaTrigV0 0.931236
ID: simpleEleId60cIso 6
ID: simpleEleId60relIso 6
ID: simpleEleId70cIso 6
ID: simpleEleId70relIso 6
ID: simpleEleId80cIso 6
ID: simpleEleId80relIso 6
ID: simpleEleId85cIso 6
ID: simpleEleId85relIso 6
ID: simpleEleId90cIso 7
ID: simpleEleId90relIso 7
ID: simpleEleId95cIso 7
ID: simpleEleId95relIso 7
ID: pf_evspi 0.837197
ID: pf_evsmu -999
ID: pf_evspi 0.647363
ID: pf_evsmu -999 
    • PF Isolation

Muons
    • Now are non-isolated muons: should we keep only the isolated?
    • MuScle calibration (
    • PF Isolation

Trigger Path filter
    • HLT_Mu17_Mu8_*
    • HLT_Ele17_CaloIdT_CaloIsoVL_TrkIdVL_TrkIsoVL_Ele8_CaloIdT_CaloIsoVL_TrkIdVL_TrkIsoVL_*
    • HLT_Mu17_Ele8_CaloIdT_CaloIsoVL_TrkIdVL_TrkIsoVL_*, HLT_Mu8_Ele17_CaloIdT_CaloIsoVL_TrkIdVL_TrkIsoVL_*
    • HLT_IsoMu24_*, HLT_IsoMu24_eta2p1_v*
    • HLT_Ele27_WP80*, HLT_Ele25_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_TriCentralPFNoPUJet30_v*

Jets

GenParticles
* Selection: all pdgId<50 or (pt>3 and eta<5) * Warning: to get a correct match with reco object, a set of input source need to be changed
    process.genParticlesPt3 = cms.EDProducer(
        "GenParticlePruner",
        src = cms.InputTag("genParticles"),
        select = cms.vstring(
        "keep (abs(pdgId) <50 ) | pt > 3 & abs(eta) < 5"
        )
    )

getattr(process,"muonMatch"+postfix).matched = cms.InputTag("genParticlesPt3")

getattr(process,"electronMatch"+postfix).matched = cms.InputTag("genParticlesPt3")

getattr(process,"photonMatch"+postfix).matched = cms.InputTag("genParticlesPt3")

getattr(process,"tauGenJets"+postfix).GenParticles = cms.InputTag("genParticlesPt3")
getattr(process,"tauMatch"+postfix).matched = cms.InputTag("genParticlesPt3")

getattr(process,"patJetPartonMatch"+postfix).matched = cms.InputTag("genParticlesPt3")
getattr(process,"patJetPartons"+postfix).matched = cms.InputTag("genParticlesPt3")

  • Added New B-tagger
D: jetBProbabilityBJetTags = 3.15469
D: jetProbabilityBJetTags = 0.132742
D: trackCountingHighPurBJetTags = -0.647959
D: trackCountingHighEffBJetTags = 0.380164
D: simpleSecondaryVertexHighEffBJetTags = -1
D: simpleSecondaryVertexHighPurBJetTags = -1
D: combinedSecondaryVertexBJetTags = 0.275713
D: combinedSecondaryVertexV1BJetTags = 0.319692
D: combinedSecondaryVertexSoftPFLeptonV1BJetTags = 0.313491

To be tested:

  • Conversion

Jet Energy Scale issue

* Apparently the Z->bb peak is misplaced (as well as the H->bb one), whle in kostya pattuples it looked fine. Why is that?

Changing the GT from START53_V23 to START52_V11 (suggested by Mia: 5.2????) does not change significantly the peak position, so I kept START53_V23

getattr(process,"pfPileUp"+postfix).checkClosestZVertex = False
improves the H->jj peak (from 96 to 103 GeV).

trying

getattr(process,"pfJets"+postfix).doAreaFastjet = True

The results are visible in these plots. The mass peak shift toward the correct direction, although not to the nominal value. At least, the GenJets and RecoJets peak agrees reasonably well.

To be noted that if the pt of the H/Z/W is boosted, then the agreement is more like the one typically found on VH analysis.

H(125)->jj

can_Hjj.png To be noted that in the H->jj related plots for the GenJets, the last row shows the results when using all GenParticle (except neutrinos) instead of GenParticle with pt>3 GeV.

Z->jj

can_Zjj.png

W->jj

can_Wjj.png

---+++++ ToDo

  • Check that all the variable needed for bb mass regression are available
  • Try to include SuperCombinedTagger and/or CSV retrained. Failed so far.
  • ...

Production

I keep track of the production of these PATTUPLES in the following spreadsheet from google drive

https://docs.google.com/spreadsheet/ccc?key=0ArGe-qLD__-KdHh0ZTI2WjNWN3RNV1NWa1ZGVmZQeFE&usp=sharing

in Tab AZh and "MC for AZh" for Data and MC, respectively

The list of pattuples are available on

/lustre/cmswork/lacaprar/Ana/PATuple_2012/Prod6/

DoubleElectron_Run2012A-22Jan2013-v1_AZh_Prod6.list
DoubleElectron_Run2012B-22Jan2013-v1_AZh_Prod6.list
DoubleElectron_Run2012C-22Jan2013-v1_AZh_Prod6.list
DoubleElectron_Run2012D-22Jan2013-v1_AZh_Prod6.list

DoubleMuParked_Run2012B-22Jan2013-v1_AZh_Prod6.list
DoubleMuParked_Run2012C-22Jan2013-v1_AZh_Prod6.list
DoubleMuParked_Run2012D-22Jan2013-v1_AZh_Prod6.list
DoubleMu_Run2012A-22Jan2013-v1_AZh_Prod6.list

MuEG_Run2012A-22Jan2013-v1_AZh_Prod6.list
MuEG_Run2012B-22Jan2013-v1_AZh_Prod6.list
MuEG_Run2012C-22Jan2013-v1_AZh_Prod6.list
MuEG_Run2012D-22Jan2013-v1_AZh_Prod6.list

SingleElectron_Run2012A-22Jan2013-v1_AZh_Prod6.list
SingleElectron_Run2012B-22Jan2013-v1_AZh_Prod6.list
SingleElectron_Run2012C-22Jan2013-v1_AZh_Prod6.list
SingleElectron_Run2012D-22Jan2013-v1_AZh_Prod6.list

SingleMu_Run2012A-22Jan2013-v1_AZh_Prod6.list
SingleMu_Run2012B-22Jan2013-v1_AZh_Prod6.list
SingleMu_Run2012C-22Jan2013-v1_AZh_Prod6.list
SingleMu_Run2012D-22Jan2013-v1_AZh_Prod6.list

DY1JetsToLL_M-50_TuneZ2Star_8TeV-madgraph_AZh_Prod6.list
DY2JetsToLL_M-50_TuneZ2Star_8TeV-madgraph_AZh_Prod6.list
DY3JetsToLL_M-50_TuneZ2Star_8TeV-madgraph_AZh_Prod6.list
DY4JetsToLL_M-50_TuneZ2Star_8TeV-madgraph_AZh_Prod6.list
DYJetsToLL_HT-200To400_TuneZ2Star_8TeV-madgraph_AZh_Prod6.list
DYJetsToLL_HT-200To400_TuneZ2Star_8TeV-madgraph_ext_AZh_Prod6.list
DYJetsToLL_HT-400ToInf_TuneZ2Star_8TeV-madgraph_AZh_Prod6.list
DYJetsToLL_HT-400ToInf_TuneZ2Star_8TeV-madgraph_ext_AZh_Prod6.list
DYJetsToLL_M-10To50_TuneZ2Star_8TeV-madgraph_AZh_Prod6.list
DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball_AZh_Prod6.list
DYJetsToLL_PtZ-100_TuneZ2star_8TeV-madgraph_AZh_Prod6.list
DYJetsToLL_PtZ-100_TuneZ2star_8TeV_ext-madgraph-tarball_AZh_Prod6.list
DYJetsToLL_PtZ-180_TuneZ2star_8TeV-madgraph_AZh_Prod6.list
DYJetsToLL_PtZ-50To70_TuneZ2star_8TeV-madgraph-tarball_AZh_Prod6.list
DYJetsToLL_PtZ-50To70_TuneZ2star_8TeV_ext-madgraph-tarball_AZh_Prod6.list
DYJetsToLL_PtZ-70To100_TuneZ2star_8TeV-madgraph-tarball_AZh_Prod6.list
DYJetsToLL_PtZ-70To100_TuneZ2star_8TeV_ext-madgraph-tarball_AZh_Prod6.list

QCD_Pt_20_MuEnrichedPt_15_TuneZ2star_8TeV_pythia6_AZh_Prod6.list
TTJets_FullLeptMGDecays_8TeV-madgraph-tauola_AZh_Prod6.list
TTJets_SemiLeptMGDecays_8TeV-madgraph-tauola_AZh_Prod6.list
TTWJets_8TeV-madgraph_Summer12_DR53X-PU_S10_START53_V7A-v1_AZh_Prod6.list
TTZJets_8TeV-madgraph_v2_Summer12_DR53X-PU_S10_START53_V7A-v1_AZh_Prod6.list
T_s-channel_TuneZ2star_8TeV-powheg-tauola_AZh_Prod6.list
T_t-channel_TuneZ2star_8TeV-powheg-tauola_AZh_Prod6.list
T_tW-channel-DR_TuneZ2star_8TeV-powheg-tauola_AZh_Prod6.list
Tbar_s-channel_TuneZ2star_8TeV-powheg-tauola_AZh_Prod6.list
Tbar_t-channel_TuneZ2star_8TeV-powheg-tauola_AZh_Prod6.list
Tbar_tW-channel-DR_TuneZ2star_8TeV-powheg-tauola_AZh_Prod6.list
WJetsToLNu_TuneZ2Star_8TeV-madgraph-tarball_AZh_Prod6.list
WWWJets_8TeV-madgraph_AZh_Prod6.list
WW_TuneZ2star_8TeV_pythia6_tauola_AZh_Prod6.list
WZ_TuneZ2star_8TeV_pythia6_tauola_AZh_Prod6.list
ZH_ZToLL_HToBB_M-125_8TeV-powheg-herwigpp_AZh_Prod6.list
ZZ_TuneZ2star_8TeV_pythia6_tauola_AZh_Prod6.list

GluGluToAToZhToLLBB_mA-250_mh-125_8TeV-pythia6-tauola_AZh_Prod6.list
GluGluToAToZhToLLBB_mA-300_mh-125_8TeV-pythia6-tauola_AZh_Prod6.list
GluGluToAToZhToLLBB_mA-350_mh-125_8TeV-pythia6-tauola_AZh_Prod6.list
SeesawTo3Lminus_M-140_FDS_TuneZ2_8TeV-madgraph_AZh_Prod6.list
SeesawTo3Lminus_M-180_FDS_TuneZ2_8TeV-madgraph_AZh_Prod6.list
SeesawTo3Lminus_M-200_FDS_TuneZ2_8TeV-madgraph_AZh_Prod6.list
SeesawTo3Lminus_M-220_FDS_TuneZ2_8TeV-madgraph_AZh_Prod6.list
SeesawTo3Lminus_M-240_FDS_TuneZ2_8TeV-madgraph_AZh_Prod6.list
SeesawTo3Lminus_M-260_FDS_TuneZ2_8TeV-madgraph_AZh_Prod6.list
SeesawTo3Lminus_M-280_FDS_TuneZ2_8TeV-madgraph_AZh_Prod6.list
SeesawTo3Lminus_M-300_FDS_TuneZ2_8TeV-madgraph_AZh_Prod6.list
SeesawTo3Lminus_M-320_FDS_TuneZ2_8TeV-madgraph_AZh_Prod6.list
SeesawTo3Lminus_M-340_FDS_TuneZ2_8TeV-madgraph_AZh_Prod6.list
SeesawTo3Lplus_M-140_FDS_TuneZ2_8TeV-madgraph_AZh_Prod6.list
SeesawTo3Lplus_M-180_FDS_TuneZ2_8TeV-madgraph_AZh_Prod6.list
SeesawTo3Lplus_M-200_FDS_TuneZ2_8TeV-madgraph_AZh_Prod6.list
SeesawTo3Lplus_M-220_FDS_TuneZ2_8TeV-madgraph_AZh_Prod6.list
SeesawTo3Lplus_M-240_FDS_TuneZ2_8TeV-madgraph_AZh_Prod6.list
SeesawTo3Lplus_M-260_FDS_TuneZ2_8TeV-madgraph_AZh_Prod6.list
SeesawTo3Lplus_M-280_FDS_TuneZ2_8TeV-madgraph_AZh_Prod6.list
SeesawTo3Lplus_M-300_FDS_TuneZ2_8TeV-madgraph_AZh_Prod6.list
SeesawTo3Lplus_M-320_FDS_TuneZ2_8TeV-madgraph_AZh_Prod6.list
SeesawTo3Lplus_M-340_FDS_TuneZ2_8TeV-madgraph_AZh_Prod6.list

---++++ TODO production

  • Rerun the production on the same datasetpath (5.2 for data and 5.3 for MC) with the latest patch

---+++ Test on 53X ReReco

Goal: validate the configuration on 5.3 ReReco data

  • the cfg does run, out of the box, on ReReco (tested on DoubleMu RunA)
  • Should I keep the Electron and Muon calibration?

Luminosity

Twiki https://twiki.cern.ch/twiki/bin/view/CMS/LumiCalc

Usage:

pixelLumiCalc.py -i /lustre/cmswork/lacaprar/Ana/CMSSW_5_3_9/src/AZh/Pattupler/test/DoubleMu_Run2012A-13Jul2012-v1/res/lumiSummary.json overview

CRAB and tools

go to CRAB wtiki page

Topic attachments
I Attachment History Action Size Date Who Comment
Unix shell scriptsh EdmCountsEvents.sh r1 manage 1.7 K 2013-05-27 - 16:43 UnknownUser  
PDFpdf can_Hjj.pdf r1 manage 38.6 K 2013-07-01 - 15:00 UnknownUser Peak position for H, Z and W ->jj
PNGpng can_Hjj.png r1 manage 24.3 K 2013-07-01 - 15:00 UnknownUser Peak position for H, Z and W ->jj
PDFpdf can_Wjj.pdf r1 manage 39.1 K 2013-07-01 - 15:00 UnknownUser Peak position for H, Z and W ->jj
PNGpng can_Wjj.png r1 manage 24.0 K 2013-07-01 - 15:00 UnknownUser Peak position for H, Z and W ->jj
PDFpdf can_Zjj.pdf r1 manage 40.4 K 2013-07-01 - 15:00 UnknownUser Peak position for H, Z and W ->jj
PNGpng can_Zjj.png r1 manage 24.7 K 2013-07-01 - 15:00 UnknownUser Peak position for H, Z and W ->jj
Texttxt contaEventiCrab.py.txt r1 manage 2.2 K 2013-05-27 - 16:41 UnknownUser  
Unix shell scriptsh creaListeFromLNL.sh r1 manage 2.0 K 2013-05-27 - 16:36 UnknownUser  
Unix shell scriptsh getLumi.sh r1 manage 0.8 K 2013-05-27 - 17:00 UnknownUser  
Texttxt patTuple_PF2PAT_cfg.py.txt r6 r5 r4 r3 r2 manage 18.8 K 2013-07-01 - 14:15 UnknownUser Better GenJets and Jets
Edit | Attach | Watch | Print version | History: r27 < r26 < r25 < r24 < r23 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r27 - 2013-11-12 - unknown
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2020 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