Skip to content

Re-miniaod #27

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
sidnarayanan wants to merge 3 commits into
base: master
Choose a base branch
from
Open

Re-miniaod #27

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
6f0a086
testing smear
Feb 9, 2017
069e6c2
setting stuff up for EG and muon corrections
Feb 10, 2017
ecfa118
updating runNtupler.py
Feb 10, 2017
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
3 changes: 3 additions & 0 deletions Ntupler/interface/ElectronFiller.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -24,6 +24,9 @@ class ElectronFiller : virtual public BaseFiller
edm::EDGetTokenT<pat::ElectronCollection> el_token;
edm::Handle<pat::ElectronCollection> el_handle;

edm::EDGetTokenT<pat::ElectronCollection> el_unsmeared_token;
edm::Handle<pat::ElectronCollection> el_unsmeared_handle;

edm::EDGetTokenT<edm::ValueMap<bool> > el_vetoid_token;
edm::Handle<edm::ValueMap<bool> > el_vetoid_handle;

Expand Down
16 changes: 8 additions & 8 deletions Ntupler/interface/METFiller.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -23,16 +23,17 @@ class METFiller : virtual public BaseFiller
virtual inline string name(){return "METFiller";};
void init(TTree *t);

bool rerun=false;

edm::EDGetTokenT<pat::METCollection> met_token;
edm::Handle<pat::METCollection> met_handle;

// rerun
edm::EDGetTokenT<reco::PFMETCollection> remet_token;
edm::Handle<reco::PFMETCollection> remet_handle;
edm::EDGetTokenT<reco::PFMETCollection> remetuncorr_token;
edm::Handle<reco::PFMETCollection> remetuncorr_handle;
edm::EDGetTokenT<pat::METCollection> cleanmu_met_token;
edm::Handle<pat::METCollection> cleanmu_met_handle;

edm::EDGetTokenT<pat::METCollection> cleaneg_met_token;
edm::Handle<pat::METCollection> cleaneg_met_handle;

edm::EDGetTokenT<pat::METCollection> unclean_met_token;
edm::Handle<pat::METCollection> unclean_met_handle;

// pf cands - used to recalc some raw METs
edm::EDGetTokenT<pat::PackedCandidateCollection> pat_token;
Expand All @@ -43,7 +44,6 @@ class METFiller : virtual public BaseFiller

pfcandtype which_cand=kPat;

void fillMETs(std::vector<const reco::Candidate*> pfs);
bool minimal=true;

private:
Expand Down
3 changes: 3 additions & 0 deletions Ntupler/interface/PhotonFiller.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -20,6 +20,9 @@ class PhotonFiller : virtual public BaseFiller
edm::EDGetTokenT<pat::PhotonCollection> pho_token;
edm::Handle<pat::PhotonCollection> pho_handle;

edm::EDGetTokenT<pat::PhotonCollection> pho_unsmeared_token;
edm::Handle<pat::PhotonCollection> pho_unsmeared_handle;

edm::EDGetTokenT<edm::ValueMap<bool> > pho_looseid_token;
edm::Handle<edm::ValueMap<bool> > pho_looseid_handle;

Expand Down
55 changes: 29 additions & 26 deletions Ntupler/plugins/Ntupler.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -80,37 +80,39 @@ Ntupler::Ntupler(const edm::ParameterSet& iConfig)
// MET FILLERS -----------------------------------------------
METFiller *pfmet = new METFiller("pfmet");
pfmet->skipEvent = skipEvent;
pfmet->rerun = false;
pfmet->met_token = consumes<pat::METCollection>(iConfig.getParameter<edm::InputTag>("pfmet"));
pfmet->cleanmu_met_token = consumes<pat::METCollection>(iConfig.getParameter<edm::InputTag>("pfmet_cleanmu"));
pfmet->cleaneg_met_token = consumes<pat::METCollection>(iConfig.getParameter<edm::InputTag>("pfmet_cleaneg"));
pfmet->unclean_met_token = consumes<pat::METCollection>(iConfig.getParameter<edm::InputTag>("pfmet_unclean"));
pfmet->pat_token = consumes<pat::PackedCandidateCollection>(iConfig.getParameter<edm::InputTag>("chsPFCands")); // these are not actually CHS, so don't worry
pfmet->which_cand = METFiller::kPat;
pfmet->which_cand = METFiller::kPat;
pfmet->minimal = false;
obj.push_back(pfmet);

METFiller *puppimet = new METFiller("puppimet");
puppimet->skipEvent = skipEvent;
puppimet->rerun = false;
puppimet->met_token = consumes<pat::METCollection>(iConfig.getParameter<edm::InputTag>("puppimet"));
puppimet->minimal = true;
obj.push_back(puppimet);


// LEPTON FILLERS --------------------------------------------
MuonFiller *muon = new MuonFiller("muon");
muon->skipEvent = skipEvent;
muon->evt = event;
muon->mu_token = consumes<pat::MuonCollection>(iConfig.getParameter<edm::InputTag>("muons"));
MuonFiller *muon = new MuonFiller("muon");
muon->skipEvent = skipEvent;
muon->evt = event;
muon->mu_token = consumes<pat::MuonCollection>(iConfig.getParameter<edm::InputTag>("muons"));
obj.push_back(muon);

ElectronFiller *electron = new ElectronFiller("electron");
electron->evt = event;
electron->skipEvent = skipEvent;
electron->el_token = consumes<pat::ElectronCollection>(iConfig.getParameter<edm::InputTag>("electrons"));
electron->el_vetoid_token = consumes<edm::ValueMap<bool> >(iConfig.getParameter<edm::InputTag>("eleVetoIdMap"));
electron->el_looseid_token = consumes<edm::ValueMap<bool> >(iConfig.getParameter<edm::InputTag>("eleLooseIdMap"));
electron->el_mediumid_token = consumes<edm::ValueMap<bool> >(iConfig.getParameter<edm::InputTag>("eleMediumIdMap"));
electron->el_tightid_token = consumes<edm::ValueMap<bool> >(iConfig.getParameter<edm::InputTag>("eleTightIdMap"));
electron->el_hltid_token = consumes<edm::ValueMap<bool> >(iConfig.getParameter<edm::InputTag>("eleHLTIdMap"));
ElectronFiller *electron = new ElectronFiller("electron");
electron->evt = event;
electron->skipEvent = skipEvent;
electron->el_token = consumes<pat::ElectronCollection>(iConfig.getParameter<edm::InputTag>("electrons"));
electron->el_unsmeared_token = consumes<pat::ElectronCollection>(edm::InputTag("slimmedElectrons"));
electron->el_vetoid_token = consumes<edm::ValueMap<bool> >(iConfig.getParameter<edm::InputTag>("eleVetoIdMap"));
electron->el_looseid_token = consumes<edm::ValueMap<bool> >(iConfig.getParameter<edm::InputTag>("eleLooseIdMap"));
electron->el_mediumid_token = consumes<edm::ValueMap<bool> >(iConfig.getParameter<edm::InputTag>("eleMediumIdMap"));
electron->el_tightid_token = consumes<edm::ValueMap<bool> >(iConfig.getParameter<edm::InputTag>("eleTightIdMap"));
electron->el_hltid_token = consumes<edm::ValueMap<bool> >(iConfig.getParameter<edm::InputTag>("eleHLTIdMap"));
electron->effArea.reset(
new EffectiveAreas(
edm::FileInPath(iConfig.getParameter<std::string>("eleEA")).fullPath()
Expand All @@ -125,16 +127,17 @@ Ntupler::Ntupler(const edm::ParameterSet& iConfig)


// PHOTON FILLER --------------------------------------------
PhotonFiller *photon = new PhotonFiller("photon");
photon->skipEvent = skipEvent;
photon->pho_token = consumes<pat::PhotonCollection>(iConfig.getParameter<edm::InputTag>("photons"));
photon->pho_looseid_token = consumes<edm::ValueMap<bool>>(iConfig.getParameter<edm::InputTag>("phoLooseIdMap"));
photon->pho_mediumid_token = consumes<edm::ValueMap<bool>>(iConfig.getParameter<edm::InputTag>("phoMediumIdMap"));
photon->pho_tightid_token = consumes<edm::ValueMap<bool>>(iConfig.getParameter<edm::InputTag>("phoTightIdMap"));
photon->iso_ch_token = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("phoChargedIsolation"));
photon->iso_nh_token = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("phoNeutralHadronIsolation"));
photon->iso_pho_token = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("phoPhotonIsolation"));
photon->iso_wch_token = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("phoWorstChargedIsolation"));
PhotonFiller *photon = new PhotonFiller("photon");
photon->skipEvent = skipEvent;
photon->pho_token = consumes<pat::PhotonCollection>(iConfig.getParameter<edm::InputTag>("photons"));
photon->pho_unsmeared_token = consumes<pat::PhotonCollection>(edm::InputTag("slimmedPhotons"));
photon->pho_looseid_token = consumes<edm::ValueMap<bool>>(iConfig.getParameter<edm::InputTag>("phoLooseIdMap"));
photon->pho_mediumid_token = consumes<edm::ValueMap<bool>>(iConfig.getParameter<edm::InputTag>("phoMediumIdMap"));
photon->pho_tightid_token = consumes<edm::ValueMap<bool>>(iConfig.getParameter<edm::InputTag>("phoTightIdMap"));
photon->iso_ch_token = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("phoChargedIsolation"));
photon->iso_nh_token = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("phoNeutralHadronIsolation"));
photon->iso_pho_token = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("phoPhotonIsolation"));
photon->iso_wch_token = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("phoWorstChargedIsolation"));
obj.push_back(photon);


Expand Down
12 changes: 7 additions & 5 deletions Ntupler/python/PandaProd_cfi.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -21,27 +21,29 @@

# offline skimming
doJetSkim = cms.bool(False),
doRecoilFilter = cms.bool(True),
doRecoilFilter = cms.bool(False),

# jet toggles
savePuppiCands = cms.bool(False),
saveCHSCands = cms.bool(False),
doCHSAK4 = cms.bool(True),
doPuppiAK4 = cms.bool(True),
doPuppiCA15 = cms.bool(True),
doCHSCA15 = cms.bool(False),
doCHSCA15 = cms.bool(True),
doPuppiAK8 = cms.bool(True),
doCHSAK8 = cms.bool(False),
doCHSAK8 = cms.bool(True),

chsAK4 = cms.InputTag("slimmedJets"),
puppiAK4 = cms.InputTag("slimmedJetsPuppi"),
#puppiAK4 = cms.InputTag("patJetsPFAK4Puppi"),
chsAK8 = cms.InputTag("packedPatJetsPFchsAK8"),
puppiAK8 = cms.InputTag("packedPatJetsPFpuppiAK8"),
chsCA15 = cms.InputTag("packedPatJetsPFchsCA15"),
puppiCA15 = cms.InputTag("packedPatJetsPFpuppiCA15"),

pfmet = cms.InputTag("slimmedMETs"),
pfmet = cms.InputTag("slimmedMETsMuEGClean"),
pfmet_cleanmu = cms.InputTag('slimmedMETs'),
pfmet_cleaneg = cms.InputTag('slimmedMETsEGClean'),
pfmet_unclean = cms.InputTag('slimmedMETsUncorrected'),
puppimet = cms.InputTag("slimmedMETsPuppi"),

puppiPFCands = cms.InputTag("puppi"),
Expand Down
97 changes: 52 additions & 45 deletions Ntupler/python/makeFatJets_cff.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -18,7 +18,7 @@ def initFatJets(process,isData):
## INITIAL SETUP ##
########################################

## Load standard PAT objects
## load chs and genparticles ##

process.pfCHS = cms.EDFilter("CandPtrSelector", src = cms.InputTag("packedPFCandidates"), cut = cms.string("fromPV"))
addGenNoNu=False
Expand All @@ -34,6 +34,13 @@ def initFatJets(process,isData):
process.fatjetInitSequence = cms.Sequence()
process.fatjetInitSequence += process.pfCHS

## load basic b-tag stuff ##
process.load("RecoBTag.ImpactParameter.impactParameter_EventSetup_cff")
process.load("RecoBTag.SecondaryVertex.secondaryVertex_EventSetup_cff")
process.load("RecoBTag.SoftLepton.softLepton_EventSetup_cff")
process.load("RecoBTag.Combined.combinedMVA_EventSetup_cff")
process.load("RecoBTag.CTagging.cTagging_EventSetup_cff")

return process.fatjetInitSequence

def makeFatJets(process,isData,pfCandidates,algoLabel,jetRadius):
Expand Down Expand Up @@ -103,6 +110,9 @@ def makeFatJets(process,isData,pfCandidates,algoLabel,jetRadius):

bTagInfosSubjets=['None']
bTagDiscriminatorsSubjets=['None']

bTagInfos = ['None']
bTagDiscriminators = ['None']

### jet clustering ###
########################################
Expand All @@ -114,40 +124,40 @@ def makeFatJets(process,isData,pfCandidates,algoLabel,jetRadius):
if not(hasattr(process,"genJetsNoNu"+rLabel)) and isMC:
addingGenJets = True
setattr(process,"genJetsNoNu"+rLabel, ak4GenJets.clone(
jetAlgorithm = cms.string(jetAlgo),
rParam = cms.double(jetRadius),
src = cms.InputTag("packedGenParticlesForJetsNoNu")
)
jetAlgorithm = cms.string(jetAlgo),
rParam = cms.double(jetRadius),
src = cms.InputTag("packedGenParticlesForJetsNoNu")
)
)
setattr(process,"PFJets"+customLabel, ak4PFJets.clone(
jetAlgorithm = cms.string(jetAlgo),
rParam = cms.double(jetRadius),
src = cms.InputTag(pfCandidates),
jetPtMin = cms.double(100)
)
jetAlgorithm = cms.string(jetAlgo),
rParam = cms.double(jetRadius),
src = cms.InputTag(pfCandidates),
jetPtMin = cms.double(100)
)
)
if not(hasattr(process,"genJetsNoNuSoftDrop"+rLabel)) and isMC:
addingGenJets = True
setattr(process,"genJetsNoNuSoftDrop"+rLabel, getattr(process,'genJetsNoNu'+rLabel).clone(
R0 = cms.double(jetRadius),
useSoftDrop = cms.bool(True),
zcut = cms.double(sdZcut),
beta = cms.double(sdBeta),
writeCompound = cms.bool(True),
useExplicitGhosts = cms.bool(True),
jetCollInstanceName=cms.string("SubJets")
)
R0 = cms.double(jetRadius),
useSoftDrop = cms.bool(True),
zcut = cms.double(sdZcut),
beta = cms.double(sdBeta),
writeCompound = cms.bool(True),
useExplicitGhosts = cms.bool(True),
jetCollInstanceName=cms.string("SubJets")
)
)
setattr(process,"PFJets"+"SoftDrop"+customLabel, getattr(process,'PFJets'+customLabel).clone(
useSoftDrop = cms.bool(True),
R0 = cms.double(jetRadius),
zcut = cms.double(sdZcut),
beta = cms.double(sdBeta),
writeCompound = cms.bool(True),
useExplicitGhosts = cms.bool(True),
jetCollInstanceName=cms.string("SubJets"),
jetPtMin = cms.double(100)
)
useSoftDrop = cms.bool(True),
R0 = cms.double(jetRadius),
zcut = cms.double(sdZcut),
beta = cms.double(sdBeta),
writeCompound = cms.bool(True),
useExplicitGhosts = cms.bool(True),
jetCollInstanceName=cms.string("SubJets"),
jetPtMin = cms.double(100)
)
)
process.dump = cms.EDAnalyzer("EventContentAnalyzer")
if addingGenJets and not(isData):
Expand All @@ -162,20 +172,20 @@ def makeFatJets(process,isData,pfCandidates,algoLabel,jetRadius):
#######################################

setattr(process,customLabel+'Njettiness',
Njettiness.clone(
src = cms.InputTag('PFJets'+customLabel),
R0 = cms.double(jetRadius),
Njets = cms.vuint32(1,2,3,4)
)
Njettiness.clone(
src = cms.InputTag('PFJets'+customLabel),
R0 = cms.double(jetRadius),
Njets = cms.vuint32(1,2,3,4)
)
)

setattr(process,customLabel+'SDKinematics',
cms.EDProducer('RecoJetDeltaRValueMapProducer',
src = cms.InputTag('PFJets'+customLabel),
matched = cms.InputTag('PFJets'+"SoftDrop"+customLabel),
distMax = cms.double(1.5),
values = cms.vstring('mass'),
valueLabels = cms.vstring('Mass'),
src = cms.InputTag('PFJets'+customLabel),
matched = cms.InputTag('PFJets'+"SoftDrop"+customLabel),
distMax = cms.double(1.5),
values = cms.vstring('mass'),
valueLabels = cms.vstring('Mass'),
)
)

Expand All @@ -197,10 +207,10 @@ def makeFatJets(process,isData,pfCandidates,algoLabel,jetRadius):

setattr(process,customLabel+'PFImpactParameterTagInfos',
pfImpactParameterTagInfos.clone(
jets = cms.InputTag('PFJets'+"SoftDrop"+customLabel,'SubJets'),
maxDeltaR = cms.double(0.4),
primaryVertex = cms.InputTag('offlineSlimmedPrimaryVertices'),
candidates = cms.InputTag('packedPFCandidates')
jets = cms.InputTag('PFJets'+"SoftDrop"+customLabel,'SubJets'),
maxDeltaR = cms.double(0.4),
primaryVertex = cms.InputTag('offlineSlimmedPrimaryVertices'),
candidates = cms.InputTag('packedPFCandidates')
)
)

Expand All @@ -222,7 +232,7 @@ def makeFatJets(process,isData,pfCandidates,algoLabel,jetRadius):
newSeq += getattr(process,customLabel+'PFInclusiveSecondaryVertexFinderTagInfos')
newSeq += getattr(process,customLabel+'PFCombinedInclusiveSecondaryVertexV2BJetTags')

### doubleb###
### doubleb ###
setattr(process, customLabel+'PFImpactParameterDoubleSVTagInfos',
pfImpactParameterTagInfos.clone(
jets = cms.InputTag("PFJets"+customLabel),
Expand Down Expand Up @@ -271,9 +281,6 @@ def makeFatJets(process,isData,pfCandidates,algoLabel,jetRadius):
newSeq += getattr(process,customLabel+'PFBoostedDoubleSVTagInfos')
newSeq += getattr(process,customLabel+'PFBoostedDoubleSecondaryVertexBJetTags')


bTagInfos = ['None']
bTagDiscriminators = ['None']

########################################
## MAKE PAT JETS ##
Expand Down
11 changes: 8 additions & 3 deletions Ntupler/src/ElectronFiller.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -32,19 +32,20 @@ int ElectronFiller::analyze(const edm::Event& iEvent){
iEvent.getByToken(el_mediumid_token,el_mediumid_handle);
iEvent.getByToken(el_tightid_token,el_tightid_handle);
iEvent.getByToken(el_hltid_token,el_hltid_handle);
iEvent.getByToken(el_unsmeared_token, el_unsmeared_handle);

unsigned int iE=-1;
for (const pat::Electron& el : *el_handle) {

++iE;

//PDebug("ElectronFiller",TString::Format("Considering electron with pT=%.3f, eta=%.3f",el.pt(),el.eta()));
// PDebug("ElectronFiller",TString::Format("Considering electron with pT=%.3f, eta=%.3f",el.pt(),el.eta()));

if (el.pt()<minPt || fabs(el.eta())>maxEta || !(el.passConversionVeto()) ) {
continue;
}

//PDebug("ElectronFiller",TString::Format("Passed kinematics and conv veto"));
// PDebug("ElectronFiller",TString::Format("Passed kinematics and conv veto"));

edm::RefToBase<pat::Electron> ref ( edm::Ref< pat::ElectronCollection >(el_handle, iE) ) ;

Expand All @@ -58,7 +59,8 @@ int ElectronFiller::analyze(const edm::Event& iEvent){
continue;
}

//PDebug("ElectronFiller",TString::Format("Passed veto"));
// PDebug("ElectronFiller",TString::Format("Passed veto"));


// compute isolation
float chiso = el.pfIsolationVariables().sumChargedHadronPt;
Expand Down Expand Up @@ -92,6 +94,9 @@ int ElectronFiller::analyze(const edm::Event& iEvent){
electron->id |= (unsigned(tight)*PElectron::kTight);
electron->id |= (unsigned(hltsafe)*PElectron::kHLTPresel);

// get the unsmeared momentum
electron->pt_unsmeared = (*el_unsmeared_handle)[iE].pt();

data->push_back(electron);
}

Expand Down
Loading