Esempio n. 1
0
process.TFileService.fileName = 'cosmic_muons.root'

from HLTrigger.HLTfilters.hltHighLevel_cfi import hltHighLevel
process.triggerFilter = hltHighLevel.clone()
process.triggerFilter.HLTPaths = ['HLT_QuadJet50_v*']
process.triggerFilter.andOr = True # = OR

process.CosmicMuons = cms.EDAnalyzer('CosmicMuons',
                                     track_src = cms.InputTag('globalCosmicMuons'),
                                     primary_vertex_src = cms.InputTag('goodOfflinePrimaryVertices'),
                                     beamspot_src = cms.InputTag('offlineBeamSpot'),
                                     general_track_src = cms.InputTag('generalTracks'),
                                     gen_particle_src = cms.InputTag('genParticles'),
                                     min_pt = cms.double(1),
                                     max_eta = cms.double(2.4),
                                     min_npxhits = cms.int32(1),
                                     min_nsthits = cms.int32(6),
                                     min_nmuhits = cms.int32(2),
                                     min_dxy = cms.double(0.01),
                                     max_relpterr = cms.double(0.5),
                                     )

process.CosmicMuonsPt1Nsthits6Dxy0 = process.CosmicMuons.clone(min_dxy = 0)
process.CosmicMuonsPt1Nsthits6Dxy100 = process.CosmicMuons.clone()
process.CosmicMuonsPt1Nsthits8Dxy0 = process.CosmicMuons.clone(min_nsthits = 8, min_dxy = 0)
process.CosmicMuonsPt1Nsthits8Dxy100 = process.CosmicMuons.clone(min_nsthits = 8)
process.CosmicMuonsPt5Nsthits6Dxy0 = process.CosmicMuons.clone(min_pt = 5, min_dxy = 0)
process.CosmicMuonsPt5Nsthits6Dxy100 = process.CosmicMuons.clone(min_pt = 5)
process.CosmicMuonsPt5Nsthits8Dxy0 = process.CosmicMuons.clone(min_pt = 5, min_nsthits = 8, min_dxy = 0)
process.CosmicMuonsPt5Nsthits8Dxy100 = process.CosmicMuons.clone(min_pt = 5, min_nsthits = 8)
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')

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),
Esempio n. 3
0
from JMTucker.Tools.BasicAnalyzer_cfg import cms, process

process.source.fileNames = ['file:gensimhlt.root']
process.TFileService.fileName = 'gencheck.root'

process.genCheck = cms.EDAnalyzer('MFVGenCheck',
                                  gen_src = cms.InputTag('genParticles'),
                                  lsp_id = cms.int32(1000021),
                                  allowed_dau_ids = cms.PSet(
                                      dummy = cms.vint32(1,-1),
                                      ),
                                  )

process.printList = cms.EDAnalyzer('JMTParticleListDrawer',
                                   src = cms.InputTag('genParticles'),
                                   printVertex = cms.untracked.bool(True),
                                   )

process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')

process.p = cms.Path(process.genCheck * process.printList)
process.p = cms.Path(process.genCheck)
Esempio n. 4
0
import sys
from JMTucker.Tools.BasicAnalyzer_cfg import cms, process
from JMTucker.Tools import SampleFiles

#SampleFiles.setup(process, 'MFVNtupleV18', 'mfv_neutralino_tau1000um_M0400', 1000)
process.source.fileNames = ['/store/user/tucker/mfv_neutralino_tau1000um_M0400/mfvmergentuple_v18/2f015b87b0f183c443d811e67eceec2e/merge_1_1_S9a.root','/store/user/tucker/mfv_neutralino_tau1000um_M0400/mfvmergentuple_v18/2f015b87b0f183c443d811e67eceec2e/merge_2_1_Qi0.root']
process.TFileService.fileName = 'resolutions.root'

process.load('JMTucker.MFVNeutralino.VertexSelector_cfi')
process.load('JMTucker.MFVNeutralino.AnalysisCuts_cfi')

mfvResolutions = cms.EDAnalyzer('MFVResolutions',
                                vertex_src = cms.InputTag('mfvSelectedVerticesTight'),
                                mevent_src = cms.InputTag('mfvEvent'),
                                which_mom = cms.int32(0),
                                max_dr = cms.double(-1),
                                max_dist = cms.double(0.005),
                                )

process.p = cms.Path(process.mfvSelectedVerticesSeq)

process.mfvResolutionsByDistCutTrks = mfvResolutions.clone()
process.p *= process.mfvResolutionsByDistCutTrks

process.mfvResolutionsByDistCutJets = mfvResolutions.clone(which_mom = 1)
process.p *= process.mfvResolutionsByDistCutJets

process.mfvResolutionsByDistCutTrksJets = mfvResolutions.clone(which_mom = 2)
process.p *= process.mfvResolutionsByDistCutTrksJets

process.p *= process.mfvAnalysisCuts
import os, sys
from JMTucker.Tools.BasicAnalyzer_cfg import cms, process

process.source.fileNames = ['/store/user/jchu/mfv_neutralino_tau1000um_M0400/jtuple_v7/5d4c2a74c85834550d3f9609274e8548/pat_1_1_hdB.root']
process.source.secondaryFileNames = cms.untracked.vstring('/store/user/tucker/mfv_neutralino_tau1000um_M0400/mfv_neutralino_tau1000um_M0400/a6ab3419cb64660d6c68351b3cff9fb0/aodpat_891_1_sZ9.root','/store/user/tucker/mfv_neutralino_tau1000um_M0400/mfv_neutralino_tau1000um_M0400/a6ab3419cb64660d6c68351b3cff9fb0/aodpat_948_2_lgB.root')
process.TFileService.fileName = 'hadronic_nb_analysis.root'

process.p = cms.EDAnalyzer('HadronicNBAnalysis',
                           trigger_results_src = cms.InputTag('TriggerResults', '', 'HLT'),
                           jet_src = cms.InputTag('selectedPatJetsPF'),
                           min_njets = cms.int32(6),
                           min_jet_pt = cms.double(50.0),
                           min_1st_jet_pt = cms.double(150.0),
                           b_discriminator_name = cms.string('combinedSecondaryVertexBJetTags'),
                           min_nbtags = cms.int32(2),
                           bdisc_min = cms.double(0.898),
                           muon_src = cms.InputTag('selectedPatMuonsPF'),
                           electron_src = cms.InputTag('selectedPatElectronsPF'),
                           )
process.p0 = cms.Path(process.p)

def de_mfv():
    if hasattr(process, 'mfvGenParticleFilter'):
        process.mfvGenParticleFilter.cut_invalid = False
    if hasattr(process, 'mfvGenVertices'):
        process.mfvGenVertices.is_mfv = False

if 'qcdht0100' in sys.argv:
    de_mfv()
    process.source.fileNames = '''/store/user/jchu/QCD_HT-100To250_TuneZ2star_8TeV-madgraph-pythia/jtuple_no_skimming_cuts_v8/a95b67af6a62c008520600808d84ac2d/pat_1_1_Shq.root
/store/user/jchu/QCD_HT-100To250_TuneZ2star_8TeV-madgraph-pythia/jtuple_no_skimming_cuts_v8/a95b67af6a62c008520600808d84ac2d/pat_2_1_RWX.root