def makeDataComponent(self, name, dataset, user, pattern, json): component = cfg.DataComponent(name=name, files=self.getFiles( dataset, user, pattern), intLumi=1, triggers=[], json=json) return component
def makeMCComponent(self, name, dataset, user, pattern): component = cfg.MCComponent( dataset=dataset, name=name, files=self.getFiles(dataset, user, pattern), xSection=1, nGenEvents=1, triggers=[], effCorrFactor=1, # skimEfficiency = self.getSkimEfficiency(dataset,user) ) # print 'Skim Efficiency for ',name,'=', component.skimEfficiency return component
def makePrivateMCComponent(self, name, dataset, files): if len(files) == 0: raise RuntimeError, "Trying to make a component %s with no files" % name # prefix filenames with dataset unless they start with "/" dprefix = dataset + "/" if files[0][0] != "/" else "" component = cfg.MCComponent( dataset=dataset, name=name, files=['root://eoscms//eos/cms%s%s' % (dprefix, f) for f in files], xSection=1, nGenEvents=1, triggers=[], effCorrFactor=1, # skimEfficiency = self.getSkimEfficiency(dataset,user) ) # print 'Skim Efficiency for ',name,'=', component.skimEfficiency return component
import copy import os import CMGTools.RootTools.fwlite.Config as cfg from CMGTools.RootTools.fwlite.Config import printComps from CMGTools.RootTools.RootTools import * # this analyzer finds the initial events before the skim skimAnalyzer = cfg.Analyzer('skimAnalyzerCount') jsonAna = cfg.Analyzer('JSONAnalyzer', ) triggerAna = cfg.Analyzer( #'TriggerAnalyzer', 'triggerBitFilter', ) # this analyzer is just there to select a list of good primary vertices. ttHVertexAna = cfg.Analyzer('VertexAnalyzer', goodVertices='offlinePrimaryVertices', vertexWeight=None, fixedWeight=1, verbose=False) # Gen Info Analyzer ttHGenAna = cfg.Analyzer( 'ttHGenLevelAnalyzer', filterHiggsDecays=[0, 15, 23, 24], verbose=False, PDFWeights=[], ) # Lepton Analyzer
#ttHLepAna.loose_muon_absIso5= 10 ttHLepAna.loose_electron_pt = 7 ttHLepAna.loose_electron_relIso = 0.4 #ttHLepAna.loose_electron_absIso = 10 ttHLepAna.ele_isoCorr = "rhoArea" # --- LEPTON SKIMMING --- ttHLepSkim.minLeptons = 1 ttHLepSkim.maxLeptons = 999 # --- JET-LEPTON CLEANING --- ttHJetAna.minLepPt = 10 # Event Analyzer for susy multi-lepton (at the moment, it's the TTH one) // Do we need had W and Top? ttHEventAna = cfg.Analyzer( 'ttHLepEventAnalyzer', minJets25=0, ) ttHIsoTrackAna.setOff = False susyCoreSequence.insert(susyCoreSequence.index(ttHCoreEventAna), ttHSVAnalyzer) susyCoreSequence.insert(susyCoreSequence.index(ttHCoreEventAna), ttHHeavyFlavourHadronAnalyzer) from CMGTools.TTHAnalysis.samples.samples_8TeV_v517 import triggers_1mu, triggers_1muHT, triggers_1eleHT # need to update the trigger MET pr HTMET? # Tree Producer treeProducer = cfg.Analyzer( 'treeProducerSusySingleSoftLepton', vectorTree=True, saveTLorentzVectors=
mc_tauEffWeight = 'effTau2011A' mc_muEffWeight = 'effMu2011A' elif period == 'Period_2011B': mc_vertexWeight = 'vertexWeight2011B' mc_tauEffWeight = 'effTau2011B' mc_muEffWeight = 'effMu2011B' elif period == 'Period_2011AB': mc_vertexWeight = 'vertexWeight2011AB' mc_tauEffWeight = 'effTau2011AB' mc_muEffWeight = 'effMu2011AB' ZMuMuAna = cfg.Analyzer('ZMuMuAnalyzer', pt1=20, pt2=20, iso1=0.1, iso2=0.1, eta1=2, eta2=2, m_min=0, m_max=200) triggerAna = cfg.Analyzer('TriggerAnalyzer') jetAna = cfg.Analyzer('SimpleJetAnalyzer', ptCut=0) effMuAna = cfg.Analyzer( 'EfficiencyAnalyzer', # recsel = 'cuts_vbtfmuon' genPdgId=13) vertexAna = cfg.Analyzer('VertexAnalyzer',
# vertexFileDir = os.environ['CMSSW_BASE'] + '/src/CMGTools/RootTools/data/Reweight/2012/Vertices' # vertexFileData = '/'.join([vertexFileDir, 'vertices_data_2012A_2012B_start_195947.root']) mc_vertexWeight = None # mc_tauEffWeight_mc = 'effTau2012MC53X' # mc_muEffWeight_mc = 'eff_2012_Rebecca_TauMu_IsoMu1753XMC' # mc_tauEffWeight = 'effTau2012ABC' # mc_muEffWeight = 'effMu2012_Rebecca_TauMu_ABC' mc_tauEffWeight_mc = 'effTau_muTau_MC_2012ABCDSummer13' mc_muEffWeight_mc = 'effMu_muTau_MC_2012ABCD' mc_tauEffWeight = 'effTau_muTau_Data_2012ABCDSummer13' mc_muEffWeight = 'effMu_muTau_Data_2012ABCDSummer13' eventSelector = cfg.Analyzer('EventSelector', toSelect=[105104]) jsonAna = cfg.Analyzer('JSONAnalyzer', ) triggerAna = cfg.Analyzer('TriggerAnalyzer') vertexAna = cfg.Analyzer( 'VertexAnalyzer', goodVertices='goodPVFilter', vertexWeight=mc_vertexWeight, fixedWeight=1, verbose=False, ) embedWeighter = cfg.Analyzer('EmbedWeighter', isRecHit=False, verbose=False)
tauScaleShift = 1.0 syncntuple = True mc_vertexWeight = 'vertexWeightFall112011AB' mc_tauEffWeight = None mc_muEffWeight = None mc_tauEffWeight_mc = 'effMediumIsoTau20MC' mc_eleEffWeight_mc = 'effEle18MC' mc_tauEffWeight = 'effTau2011AB_TauEle' mc_eleEffWeight = 'effEle2011AB' jsonAna = cfg.Analyzer( 'JSONAnalyzer', ) triggerAna = cfg.Analyzer( 'TriggerAnalyzer', verbose = False ) vertexAna = cfg.Analyzer( 'VertexAnalyzer', goodVertices = 'goodPVFilter', vertexWeight = mc_vertexWeight, # fixedWeight = 1, verbose = False )
from CMGTools.H2TauTau.proto.samples.data_diTau import * # from CMGTools.H2TauTau.proto.samples.embed import * from CMGTools.H2TauTau.proto.samples.ewk import * import copy import re import CMGTools.RootTools.fwlite.Config as cfg #from CMGTools.RootTools.yellowreport.YRParser import yrparser Higgsgg110 = cfg.MCComponent(name='Higgsgg110', files=[], xSection=19.84 * 6.37e-2, nGenEvents=199988, triggers=[], effCorrFactor=1) Higgsgg115 = cfg.MCComponent(name='Higgsgg115', files=[], xSection=18.13 * 6.37e-2, nGenEvents=196002, triggers=[], effCorrFactor=1) Higgsgg120 = cfg.MCComponent(name='Higgsgg120', files=[], xSection=16.63 * 6.37e-2, nGenEvents=199981, triggers=[], effCorrFactor=1) Higgsgg125 = cfg.MCComponent(name='Higgsgg125',
puFileMC = '/afs/cern.ch/user/a/agilbert/public/HTT_Pileup/28-09-12/MC_Summer12_PU_S10-600bins.root' puFileData = '/afs/cern.ch/user/a/agilbert/public/HTT_Pileup/28-09-12/Data_Pileup_2012_HCP-600bins.root' vertexFileDir = os.environ[ 'CMSSW_BASE'] + '/src/CMGTools/RootTools/data/Reweight/2012/Vertices' vertexFileData = '/'.join( [vertexFileDir, 'vertices_data_2012A_2012B_start_195947.root']) mc_vertexWeight = None mc_tauEffWeight_mc = 'effTau2012MC53X' mc_muEffWeight_mc = 'eff_2012_Rebecca_TauMu_IsoMu1753XMC' mc_tauEffWeight = 'effTau2012ABC' mc_muEffWeight = 'effMu2012_Rebecca_TauMu_ABC' eventSelector = cfg.Analyzer('EventSelector', toSelect=[105104]) jsonAna = cfg.Analyzer('JSONAnalyzer', ) triggerAna = cfg.Analyzer('TriggerAnalyzer') vertexAna = cfg.Analyzer( 'VertexAnalyzer', goodVertices='goodPVFilter', vertexWeight=mc_vertexWeight, fixedWeight=1, verbose=False, ) embedWeighter = cfg.Analyzer('EmbedWeighter', verbose=False)
mc_tauEffWeight = None mc_tauEffWeight_mc = 'effLooseTau15MC' if period == 'Period_2011A': mc_vertexWeight = 'vertexWeightFall112invfb' mc_tauEffWeight = 'effTau2011A' elif period == 'Period_2011B': mc_vertexWeight = 'vertexWeightFall112011B' mc_tauEffWeight = 'effTau2011B' elif period == 'Period_2011AB': mc_vertexWeight = 'vertexWeightFall112011AB' mc_tauEffWeight = 'effTau2011AB' if runOnMC: mc_tauEffWeight = 'effTau1fb' triggerAna = cfg.Analyzer('TriggerAnalyzer') TauTauAna = cfg.Analyzer( 'TauTauAnalyzer', pt1=20, pt2=20, iso1=999, iso2=999, eta1=999, eta2=999, m_min=10, m_max=99999, diLeptonCutString='cuts_baseline', # No trigger matching for the moment since MC only contains too high threshold trigger #triggerMap = pathsAndFilters )
return lepton.relIso(0.5)<0.1 def idIsoMuon(muon): return idMuon(muon) and isoLepton(muon) def passLepton(lepton): return True trigMap = { 'HLT_IsoMu15_v5':'hltSingleMuIsoL3IsoFiltered15', 'HLT_IsoMu15_v14':'hltSingleMuIsoL3IsoFiltered15' } effMuAnaStd = cfg.Analyzer( 'EfficiencyAnalyzer_std', # recsel = 'cuts_vbtfmuon', recselFun = isoLepton, # refselFun = idMuon, triggerMap = trigMap, instance = 'cmgMuonSelStdLep', type = 'std::vector<cmg::Muon>', genPdgId = 13 ) effMuAnaPF = cfg.Analyzer( 'EfficiencyAnalyzer_pf', # recsel = 'cuts_vbtfmuon', recselFun = isoLepton, # refselFun = idMuon, triggerMap = trigMap, instance = 'cmgMuonSel', type = 'std::vector<cmg::Muon>', genPdgId = 13 )
def createTreeProducer(ana): tp = cfg.Analyzer('_'.join(['TauTauLLTreeProducer', 'TauTauLLAnalyzer']), anaName='TauTauLLAnalyzer') return tp
import copy import os import CMGTools.RootTools.fwlite.Config as cfg tautaullAna = cfg.Analyzer('TauTauLLAnalyzer', ) def createTreeProducer(ana): tp = cfg.Analyzer('_'.join(['TauTauLLTreeProducer', 'TauTauLLAnalyzer']), anaName='TauTauLLAnalyzer') return tp # for debugging stopper = cfg.Analyzer('Stopper') sequence = cfg.Sequence([ tautaullAna, createTreeProducer(tautaullAna), ])
puFileMC = '/afs/cern.ch/user/a/agilbert/public/HTT_Pileup/12-09-12/MC_Summer12_PU_S7.root' puFileData = '/afs/cern.ch/user/a/agilbert/public/HTT_Pileup/12-09-12/Data_Pileup_2012.root' vertexFileDir = os.environ[ 'CMSSW_BASE'] + '/src/CMGTools/RootTools/data/Reweight/2012/Vertices' vertexFileData = '/'.join( [vertexFileDir, 'vertices_data_2012A_2012B_start_195947.root']) mc_vertexWeight = None mc_tauEffWeight_mc = 'eff2012Tau20MC_TauEle' mc_eleEffWeight_mc = 'eff2012Ele20MC' mc_tauEffWeight = 'effTau2012AB_TauEle' mc_eleEffWeight = 'effEle2012AB' jsonAna = cfg.Analyzer('JSONAnalyzer', ) triggerAna = cfg.Analyzer('TriggerAnalyzer', # verbose = True ) vertexAna = cfg.Analyzer('VertexAnalyzer', goodVertices='goodPVFilter', vertexWeight=mc_vertexWeight, fixedWeight=1, verbose=False) embedWeighter = cfg.Analyzer('EmbedWeighter', verbose=False) pileUpAna = cfg.Analyzer('PileUpAnalyzer', true=True)
########################################################## ## SUSY CONFIGURATION FOR LEPTON ID STUDIES ## ## makes trees with one entry per lepton, not per event ## ########################################################## import CMGTools.RootTools.fwlite.Config as cfg from CMGTools.RootTools.fwlite.Config import printComps from CMGTools.RootTools.RootTools import * PDFWeights = [] #PDFWeights = [ ("CT10",53), ("MSTW2008lo68cl",41), ("NNPDF21_100",101) ] #PDFWeights = [ ("cteq61",41) ] # this analyzer finds the initial events before the skim skimAnalyzer = cfg.Analyzer('skimAnalyzerCount') eventSelector = cfg.Analyzer( 'EventSelector', toSelect=[ # here put the event numbers (actual event numbers from CMSSW) ]) jsonAna = cfg.Analyzer('JSONAnalyzer', ) # this analyzer is just there to select a list of good primary vertices. ttHVertexAna = cfg.Analyzer('VertexAnalyzer', vertexWeight=None, fixedWeight=1, verbose=False) # this analyzer actually does the pile-up reweighting.
import CMGTools.RootTools.fwlite.Config as cfg from CMGTools.RootTools.fwlite.Config import printComps from CMGTools.RootTools.RootTools import * reader = cfg.Analyzer( 'ObjectReader', # muonCol = ('cmgMuonSel', 'vector<cmg::Muon>'), # electronCol = ('cmgElectronSel', 'vector<cmg::Electron>'), muonCol = None, electronCol = None, # only if you have PFJets jetCol = ('ak5PFJets', 'std::vector<reco::PFJet>'), # comment above and uncomment below when doing gensim # jetCol = None, simJetCol = ('jets', 'std::vector<reco::LeafCandidate>'), genJetCol = ('ak5GenJets', 'std::vector<reco::GenJet>'), # comment above and uncomment below when doing gensim # genJetCol = ('genJets', 'std::vector<reco::LeafCandidate>'), genPartCol = ('genParticles', 'std::vector<reco::GenParticle>'), simPartCol = ('pfsim', 'std::vector<reco::LeafCandidate>'), recPartCol = None ) jetAna = cfg.Analyzer( 'PFSimJetAnalyzer', jetPt = 10., jetEta = 5.0, btagSFseed = 123456, relaxJetId = True,
import copy import os import CMGTools.RootTools.fwlite.Config as cfg from CMGTools.H2TauTau.proto.samples.getFiles import getFiles DYJets = cfg.MCComponent(name='DYJets', files=[], xSection=3048., nGenEvents=1, triggers=[], effCorrFactor=1) WJets = cfg.MCComponent(name='WJets', files=[], xSection=31314., nGenEvents=1, triggers=[], effCorrFactor=1) TTJets = cfg.MCComponent(name='TTJets', files=[], xSection=165.8, nGenEvents=1, triggers=[], effCorrFactor=1) mc_basicEWK = [DYJets, WJets, TTJets]
useParked = False MC_trigger_matching = True puFileMC = '/afs/cern.ch/user/a/agilbert/public/HTT_Pileup/07-01-13/MC_Summer12_PU_S10-600bins.root' puFileData = '/afs/cern.ch/user/a/agilbert/public/HTT_Pileup/07-01-13/Data_Pileup_2012_Moriond-600bins.root' mc_vertexWeight = None mc_tauEffWeight_mc = 'effLooseTau15MC' mc_tauEffWeight = 'eff2012IsoTau19fb' mc_jetEffWeight = 'eff2012Jet19fb' if use1prong or useParked: mc_tauEffWeight='eff2012IsoTau1prong19fb' triggerAna = cfg.Analyzer( 'TriggerAnalyzer', verbose = False ) TauTauAna = cfg.Analyzer( 'TauTauAnalyzer', pt1 = 0, pt2 = 0, iso1 = 999, iso2 = 999, eta1 = 999, eta2 = 999, m_min = 0, m_max = 99999, diLeptonCutString = '', # No trigger matching for the moment since MC only contains too high threshold trigger triggerMap = pathsAndFilters,
import copy import os import CMGTools.RootTools.fwlite.Config as cfg from CMGTools.RootTools.fwlite.Config import printComps from CMGTools.RootTools.RootTools import * fourJetAna = cfg.Analyzer( 'FourJetAnalyzer', npfj=(3, 5), # at least 8 PF particle in 3 of the 4 jets ntkj=(4, 1), # at least 1 PF charged hadron in 4 of the 4 jets minM=0.5, # No jet with a jet mass smaller than 0.5 GeV mVis=180., # total visible mass in excess of 180 GeV, chi2= 1000., # total chi**2 of the beta4 "fit" smaller than 1000. (i.e., all fit energies > 0.) pair_mass=0., # No jet pair should have a mass smaller than 30 GeV pair_cijkl= 1000., # Tue smallest sum of cos(ij)+cos(kl) should be larger than -1.1 pair_sumtet= 0., # The sum of the four smallest jet-jet angles must be larger than 350 degrees h_mass=105., # The Higgs candidate mass should be larger than 90 GeV z_mass=( 80., 110. ), # The Z candidate jet pair should have a mass closer to the Z mass than 100 GeV ) def createTreeProducer(ana): tp = cfg.Analyzer('_'.join(['FourJetTreeProducer', 'FourJetAnalyzer']), anaName='FourJetAnalyzer') return tp
import glob import CMGTools.RootTools.fwlite.Config as cfg from CMGTools.RootTools.fwlite.Config import printComps from CMGTools.RootTools.RootTools import * reader = cfg.Analyzer( 'ObjectReader', # muonCol = ('cmgMuonSel', 'vector<cmg::Muon>'), # electronCol = ('cmgElectronSel', 'vector<cmg::Electron>'), muonCol = None, electronCol = None, # only if you have PFJets jetCol = None, # jetCol = ('ak5PFJets', 'std::vector<reco::PFJet>'), simJetCol = ('jets', 'std::vector<reco::LeafCandidate>'), # genJetCol = ('ak5GenJets', 'std::vector<reco::GenJet>'), genJetCol = ('genJets', 'std::vector<reco::LeafCandidate>'), genPartCol = ('genParticles', 'std::vector<reco::GenParticle>'), simPartCol = ('pfsim', 'std::vector<reco::LeafCandidate>'), ) jetAna = cfg.Analyzer( 'PFSimJetAnalyzer', jetPt = 10., jetEta = 5.0, btagSFseed = 123456, relaxJetId = True,
'HLT_Ele15_{blah}_LooseIsoPFTau15_v4'.format(blah=blah), 'HLT_Ele15_{blah}_LooseIsoPFTau20_v6'.format(blah=blah), 'HLT_Ele15_{blah}_LooseIsoPFTau20_v8'.format(blah=blah), 'HLT_Ele15_{blah}_LooseIsoPFTau20_v9'.format(blah=blah), 'HLT_Ele15_{blah}_TightIsoPFTau20_v2'.format(blah=blah), 'HLT_Ele18_{blah}_MediumIsoPFTau20_v1'.format(blah=blah) ] data_triggers_11B = [ 'HLT_Ele18_{blah}_MediumIsoPFTau20_v1'.format(blah=blah), 'HLT_Ele20_{blah}_MediumIsoPFTau20_v5'.format(blah=blah), 'HLT_Ele20_{blah}_MediumIsoPFTau20_v6'.format(blah=blah) ] triggerAna = cfg.Analyzer( 'TriggerAnalyzer' ) tauEleAna = cfg.Analyzer( 'TauEleAnalyzer', pt1 = 20, pt2 = 20, iso1 = 999, iso2 = 0.1, eta1 = 999, eta2 = 2.1, # cutString2 = 'cuts_vbtf80ID', # eta2 = 1.4, m_min = 10, m_max = 99999, diLeptonCutString = 'cuts_baseline'
import copy import os import CMGTools.RootTools.fwlite.Config as cfg from CMGTools.H2TauTau.proto.samples.getFiles import getFiles jsondir = '/afs/cern.ch/cms/CAF/CMSCOMM/COMM_DQM/certification/Collisions11/7TeV' promptjson = '/'.join([ jsondir, 'Prompt/Cert_160404-180252_7TeV_PromptReco_Collisions11_JSON.txt' ]) embed_Run2011A_May10ReReco_v1 = cfg.EmbedComponent( name='embed_Run2011A_May10ReReco_v1', files=[], triggers=[], json='/'.join([ jsondir, 'Reprocessing/Cert_160404-163869_7TeV_May10ReReco_Collisions11_JSON_v3.txt' ])) embed_Run2011A_PromptReco_v4 = cfg.EmbedComponent( name='embed_Run2011A_PromptReco_v4', files=[], triggers=[], json=promptjson) embed_Run2011A_05Aug2011_v1 = cfg.EmbedComponent( name='embed_Run2011A_05Aug2011_v1', files=[], triggers=[], json='/'.join([ jsondir,
import copy import os import CMGTools.RootTools.fwlite.Config as cfg from CMGTools.RootTools.fwlite.Config import printComps from CMGTools.WMass.triggerMap import pathsAndFilters #from CMGTools.WMass.analyzers.WAnalyzer import WAnalyzer #from CMGTools.WMass.analyzers.ZAnalyzer_notrig import ZAnalyzer_notrig jsonAna = cfg.Analyzer('JSONAnalyzer', ) triggerAna = cfg.Analyzer('TriggerAnalyzer', keepFailingEvents=True) vertexAna = cfg.Analyzer('VertexAnalyzer', fixedWeight=1., keepFailingEvents=True) WAna = cfg.Analyzer('WAnalyzer', recoilcut=1000, pfmetcut=0, jetptcut=1000, pt=30, eta=2.1, iso=0.5, savegenp=True, verbose=True, triggerMap=pathsAndFilters) WtreeProducer = cfg.Analyzer('WTreeProducer') ZAna = cfg.Analyzer('ZAnalyzer', recoilcut=1000,
import CMGTools.RootTools.fwlite.Config as cfg import copy # the luminosities for these components correspond to the # HLT triggers used in the H->tau tau analysis # they might differ from the luminosities you want for a # Z-> mu mu or Z-> ee analysis. data_Run2012A = cfg.DataComponent( name='data_Run2012A', files=[], intLumi=890., # 806.194 from Colin triggers=[], json=None) data_Run2012B = cfg.DataComponent(name='data_Run2012B', files=[], intLumi=4429., triggers=[], json=None) data_Run2012C = cfg.DataComponent( name='data_Run2012C', files=[], # | HLT_IsoMu17_eta2p1_LooseIsoPFTau20_v3 | 6719 | 495.003(/pb) | 495.003(/pb) | intLumi=7153., triggers=[], json=None) data_Run2012D = cfg.DataComponent( name='data_Run2012D',
import copy import os import CMGTools.RootTools.fwlite.Config as cfg fourJetAna = cfg.Analyzer( 'FourJetAnalyzer', npfj=(3, 5), # at least 8 PF particle in 3 of the 4 jets ntkj=(4, 1), # at least 1 PF charged hadron in 4 of the 4 jets minM=0.5, # No jet with a jet mass smaller than 2.5 GeV mVis=180., # total visible mass in excess of 180 GeV, chi2= 1000., # total chi**2 of the beta4 "fit" smaller than 1000. (i.e., all fit energies > 0.) pair_mass=0., # No jet pair should have a mass smaller than 30 GeV pair_cijkl= 1000., # Tue smallest sum of cos(ij)+cos(kl) should be larger than -1.1 pair_sumtet= 0., # The sum of the four smallest jet-jet angles must be larger than 350 degrees h_mass=105., # The Higgs candidate mass should be larger than 90 GeV z_mass=( 80., 110. ), # The Z candidate jet pair should have a mass closer to the Z mass than 100 GeV ) def createTreeProducer(ana): tp = cfg.Analyzer('_'.join(['FourJetTreeProducer', 'FourJetAnalyzer']), anaName='FourJetAnalyzer') return tp # for debugging
import CMGTools.RootTools.fwlite.Config as cfg # from CMGTools.H2TauTau.proto.samples.sampleShift import sampleShift DYJets = cfg.MCComponent( name = 'DYJets', files = [], xSection = 2008.4*3., nGenEvents = 1, triggers = [], effCorrFactor = 1 ) DY1Jets = cfg.MCComponent( name = 'DY1Jets', files = [], xSection = 561.0, nGenEvents = 1, triggers = [], effCorrFactor = 1 ) DY2Jets = cfg.MCComponent( name = 'DY2Jets', files = [], xSection = 181., nGenEvents = 1, triggers = [], effCorrFactor = 1 ) DY3Jets = cfg.MCComponent( name = 'DY3Jets', files = [], xSection = 51.1,
def createTreeProducer(ana): tp = cfg.Analyzer('_'.join(['FourJetTreeProducer', 'FourJetAnalyzer']), anaName='FourJetAnalyzer') return tp
import copy import os import CMGTools.RootTools.fwlite.Config as cfg from CMGTools.RootTools.fwlite.Config import printComps from CMGTools.RootTools.RootTools import * PDFWeights = [] #PDFWeights = [ ("CT10",53), ("MSTW2008lo68cl",41), ("NNPDF21_100",101) ] #PDFWeights = [ ("cteq61",41) ] # this analyzer finds the initial events before the skim skimAnalyzer = cfg.Analyzer( 'skimAnalyzerCount', useLumiBlocks = False, # since we run on a non-CMG dataset, there's no info on lumi blocks, and this way we use less memory ) eventSelector = cfg.Analyzer( 'EventSelector', toSelect = [ # here put the event numbers (actual event numbers from CMSSW) ] ) jsonAna = cfg.Analyzer( 'JSONAnalyzer', ) triggerAna = cfg.Analyzer( #'TriggerAnalyzer', 'triggerBitFilter', )
import copy import os import CMGTools.RootTools.fwlite.Config as cfg from CMGTools.RootTools.fwlite.Config import printComps from CMGTools.WMass.triggerMap import triggers_mu jsonAna = cfg.Analyzer('JSONAnalyzer', ) triggerAna = cfg.Analyzer('triggerBitFilter', # keepFailingEvents = False ) vertexAna = cfg.Analyzer( 'VertexAnalyzer', allVertices='slimmedPrimaryVertices', goodVertices='slimmedPrimaryVertices', vertexWeight=None, fixedWeight=1, verbose=False, ) WAna = cfg.Analyzer('WAnalyzer', recoilcut=1000, pfmetcut=0, jetptcut=1000, pt=30, eta=2.1, iso=0.5, savegenp=False, verbose=True, triggerBits={'SingleMu': triggers_mu},