import sys
from JMTucker.Tools.BasicAnalyzer_cfg import cms, process
process.TFileService.fileName = cms.string('evids.root')
process.evids = cms.EDAnalyzer('EventIdRecorder')
process.p = cms.Path(process.evids)
process.source.duplicateCheckMode = cms.untracked.string('noDuplicateCheck')
if __name__ == '__main__' and hasattr(sys, 'argv') and 'submit' in sys.argv:
import JMTucker.Tools.Samples as Samples
bkg_samples = Samples.ttbar_samples + Samples.qcd_samples
from JMTucker.Tools.CRABSubmitter import CRABSubmitter
from JMTucker.Tools.SampleFiles import SampleFiles
def crab_cfg_modifier(sample):
return [('CMSSW', 'lumi_mask', '/uscms/home/tucker/mfvrecipe/HalfMCLists/jsons/%s.json' % sample.name)]
cs = CRABSubmitter('HalfMCLists_TestJsons',
use_ana_dataset = True,
total_number_of_lumis = -1,
lumis_per_job = 3500,
crab_cfg_modifier = crab_cfg_modifier,
)
cs.submit_all(bkg_samples)
import sys
from JMTucker.Tools.BasicAnalyzer_cfg import cms, process
process.TFileService.fileName = cms.string('evids.root')
process.source.fileNames = ['/store/user/tucker/TTJets_HadronicMGDecays_8TeV-madgraph/mfvntuple_v17/f1615d49c4ae9d19e350601d059c4237/ntuple_600_1_BPu.root']
process.evids = cms.EDAnalyzer('EventIdRecorder')
process.veto = cms.EDFilter('EventIdVeto', list_fn = cms.string('ttbarhadronic.txt.gz'), use_run = cms.bool(False))
process.p = cms.Path(~process.veto * process.evids)
process.options.wantSummary = True
process.source.duplicateCheckMode = cms.untracked.string('noDuplicateCheck')
if __name__ == '__main__' and hasattr(sys, 'argv') and 'submit' in sys.argv:
import JMTucker.Tools.Samples as Samples
bkg_samples = Samples.ttbar_samples + Samples.qcd_samples
from JMTucker.Tools.CRABSubmitter import CRABSubmitter
from JMTucker.Tools.SampleFiles import SampleFiles
def pset_modifier(sample):
to_add = []
to_modify = []
to_add.append("process.veto.list_fn = '%s.txt.gz'" % sample.name)
return to_add, to_modify
def crab_cfg_modifier(sample):
return [('USER', 'additional_input_files', '/uscms/home/tucker/mfvrecipe/HalfMCLists/%s.txt.gz' % sample.name)]
cs = CRABSubmitter('HalfMCLists_TestLists',
#!/usr/bin/env python
import sys
from JMTucker.Tools.general import typed_from_argv
from JMTucker.Tools.CMSSWTools import *
from JMTucker.Tools.BasicAnalyzer_cfg import cms, process
process.TFileService.fileName = 'mctruth.root'
process.load('JMTucker.MFVNeutralino.WeightProducer_cfi')
process.p = cms.Path(process.mfvWeight)
process.mfvMovedTree = cms.EDAnalyzer('MFVMovedTracksTreer',
event_src = cms.InputTag('mfvEvent'),
weight_src = cms.InputTag('mfvWeight'),
mover_src = cms.string(''),
vertices_src = cms.InputTag('mfvVerticesAux'),
max_dist2move = cms.double(-1),
apply_presel = cms.bool(False),
njets_req = cms.uint32(0),
nbjets_req = cms.uint32(0),
for_mctruth = cms.bool(True),
)
process.p *= process.mfvMovedTree
file_event_from_argv(process)
if __name__ == '__main__' and hasattr(sys, 'argv') and 'submit' in sys.argv:
from JMTucker.Tools.CRAB3Submitter import CRABSubmitter
import JMTucker.Tools.Samples as Samples
from JMTucker.Tools.BasicAnalyzer_cfg import cms, process, add_analyzer
process.maxEvents.input = 1000
process.source.fileNames = ['file:/uscms/home/tucker/nobackup/fromt3/mfv_neutralino_tau1000um_M0400_jtuple_v6_547d3313903142038335071634b26604_pat_1_1_Dpa.root']
process.TFileService.fileName = 'track_play.root'
#process.load('Configuration.Geometry.GeometryIdeal_cff')
#process.load('Configuration.StandardSequences.MagneticField_AutoFromDBCurrent_cff')
#process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
#process.GlobalTag.globaltag = 'START53_V21::All'
#process.load('TrackingTools.TransientTrack.TransientTrackBuilder_cfi')
ana = process.mfvTrackingPlay = cms.EDAnalyzer('MFVTrackPlay',
track_src = cms.InputTag('generalTracks'),
dxy_cut = cms.double(0),
quality_cut = cms.string(''),
jet_src = cms.InputTag('selectedPatJetsPF'),
jet_pt_cut = cms.double(20),
delta_r_jet_cut = cms.double(1e10),
)
process.p = cms.Path(process.mfvTrackingPlay)
cuts = [
('Qpt0p5eta2p4DRJet0p4', 'pt > 0.5 && abs(eta) < 2.4'),
('Qpt0p5eta2p4DRJet0p4JetPt80', 'pt > 0.5 && abs(eta) < 2.4'),
('Qpt5eta2p4', 'pt > 5 && abs(eta) < 2.4'),
('Qpt5eta2p48hits', 'pt > 5 && abs(eta) < 2.4 && hitPattern.numberOfValidHits >= 8'),
('Qpt5eta2p48hitsdxy0p1', 'pt > 5 && abs(eta) < 2.4 && hitPattern.numberOfValidHits >= 8'),
('Qpt20eta2p48hitsdxy0p1', 'pt > 20 && abs(eta) < 2.4 && hitPattern.numberOfValidHits >= 8'),
]
process.load('CommonTools.ParticleFlow.goodOfflinePrimaryVertices_cfi')
process.goodOfflinePrimaryVertices.filter = cms.bool(True)
process.load('JMTucker.MFVNeutralino.Vertexer_cff')
process.p = cms.EDAnalyzer('LeptonicNBAnalysis',
trigger_results_src = cms.InputTag('TriggerResults', '', 'HLT'),
primary_vertex_src = cms.InputTag('goodOfflinePrimaryVertices'),
secondary_vertex_src = cms.InputTag('mfvSelectedVerticesTight'),
muon_src = cms.InputTag('selectedPatMuonsPF'),
min_nmuons = cms.int32(1),
jet_src = cms.InputTag('selectedPatJetsPF'),
min_jet_pt = cms.double(30.0),
max_jet_eta = cms.double(2.4),
min_njets = cms.int32(6),
b_discriminator_name = cms.string('combinedSecondaryVertexBJetTags'),
bdisc_min = cms.double(0.679),
min_nbtags = cms.int32(1),
muon_stations_min = cms.bool(True),
muon_trackerhits_min = cms.bool(True),
muon_pixelhits_min = cms.bool(True),
muon_dxy_max = cms.bool(True),
muon_dz_max = cms.bool(True),
muon_eta_max = cms.bool(True),
muon_tkpterror_max = cms.bool(True),
muon_trkkink_max = cms.bool(True),
muon_pt_min = cms.bool(True),
muon_iso_max = cms.bool(True),
)
process.no_stations_min = process.p.clone(muon_stations_min = False)
import sys, os
from JMTucker.Tools.BasicAnalyzer_cfg import cms, process
process.TFileService.fileName = 'simple_trigger_efficiency.root'
process.genMus = cms.EDFilter('CandViewSelector', src = cms.InputTag('genParticles'), cut = cms.string('abs(pdgId) == 13 && abs(mother.pdgId) == 24'))
process.genMuCount = cms.EDFilter('CandViewCountFilter', src = cms.InputTag('genMus'), minNumber = cms.uint32(1))
process.genEls = cms.EDFilter('CandViewSelector', src = cms.InputTag('genParticles'), cut = cms.string('abs(pdgId) == 11 && abs(mother.pdgId) == 24'))
process.genElCount = cms.EDFilter('CandViewCountFilter', src = cms.InputTag('genEls'), minNumber = cms.uint32(1))
process.genMusInAcc = cms.EDFilter('CandViewSelector', src = cms.InputTag('genParticles'), cut = cms.string('abs(pdgId) == 13 && abs(mother.pdgId) == 24 && pt > 26 && abs(eta) < 2.1'))
process.genElsInAcc = cms.EDFilter('CandViewSelector', src = cms.InputTag('genParticles'), cut = cms.string('abs(pdgId) == 11 && abs(mother.pdgId) == 24 && pt > 30 && abs(eta) < 2.5'))
process.genMuInAccCount = cms.EDFilter('CandViewCountFilter', src = cms.InputTag('genMusInAcc'), minNumber = cms.uint32(1))
process.genElInAccCount = cms.EDFilter('CandViewCountFilter', src = cms.InputTag('genElsInAcc'), minNumber = cms.uint32(1))
process.RandomNumberGeneratorService = cms.Service('RandomNumberGeneratorService')
process.RandomNumberGeneratorService.SimpleTriggerEfficiency = cms.PSet(initialSeed = cms.untracked.uint32(1219))
process.RandomNumberGeneratorService.SimpleTriggerEfficiencyMu = cms.PSet(initialSeed = cms.untracked.uint32(1220))
process.RandomNumberGeneratorService.SimpleTriggerEfficiencyMuInAcc = cms.PSet(initialSeed = cms.untracked.uint32(1221))
#import prescales
process.SimpleTriggerEfficiency = cms.EDAnalyzer('SimpleTriggerEfficiency',
trigger_results_src = cms.InputTag('TriggerResults', '', 'HLT'),
prescale_paths = cms.vstring(), #*prescales.prescales.keys()),
prescale_values = cms.vuint32(), #*[o for l,h,o in prescales.prescales.itervalues()]),
)
process.SimpleTriggerEfficiencyMu = process.SimpleTriggerEfficiency.clone()
process.SimpleTriggerEfficiencyMuInAcc = process.SimpleTriggerEfficiency.clone()
process.SimpleTriggerEfficiencyEl = process.SimpleTriggerEfficiency.clone()
