Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-04-06 12:15:38

0001 import FWCore.ParameterSet.Config as cms
0002 
0003 from PhysicsTools.PatAlgos.tools.helpers import getPatAlgosToolsTask
0004 
0005 process = cms.Process("bphAnalysis")
0006 
0007 patAlgosToolsTask = getPatAlgosToolsTask(process)
0008 
0009 #process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
0010 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(1000) )
0011 
0012 process.load("Configuration.Geometry.GeometryRecoDB_cff")
0013 process.load("Configuration.StandardSequences.MagneticField_cff")
0014 process.load('Configuration.StandardSequences.Services_cff')
0015 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0016 process.load('FWCore.MessageService.MessageLogger_cfi')
0017 process.load('Configuration.EventContent.EventContent_cff')
0018 process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
0019 process.load('Configuration.StandardSequences.MagneticField_AutoFromDBCurrent_cff')
0020 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0021 patAlgosToolsTask.add(process.MEtoEDMConverter)
0022 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0023 process.load("TrackingTools/TransientTrack/TransientTrackBuilder_cfi")
0024 
0025 process.MessageLogger.cerr.FwkReport.reportEvery = 100
0026 process.options = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) )
0027 
0028 process.source = cms.Source("PoolSource",fileNames = cms.untracked.vstring(
0029 #
0030 ### use this to access the nearest copy of the input file, querying the catalog
0031     '/store/mc/RunIISpring15DR74/BsToJpsiPhi_BMuonFilter_TuneCUEP8M1_13TeV-pythia8-evtgen/AODSIM/Asympt25nsRaw_MCRUN2_74_V9-v1/50000/D0D90725-1D61-E511-B812-0025907277CE.root'
0032 #
0033 ### use this to access the input file if by any reason you want to specify 
0034 ### the data server
0035 #    'root://xrootd-cms.infn.it//store/mc/RunIISpring15DR74/BsToJpsiPhi_BMuonFilter_TuneCUEP8M1_13TeV-pythia8-evtgen/AODSIM/Asympt25nsRaw_MCRUN2_74_V9-v1/50000/D0D90725-1D61-E511-B812-0025907277CE.root'
0036 #
0037 ### use this to access an input file locally available
0038 #    'file:/...complete_file_path.../D0D90725-1D61-E511-B812-0025907277CE.root'
0039 ))
0040 
0041 from Configuration.AlCa.GlobalTag import GlobalTag
0042 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:run2_mc', '')
0043 
0044 process.load('PhysicsTools.PatAlgos.producersLayer1.patCandidates_cff')
0045 patAlgosToolsTask.add(process.patCandidatesTask)
0046 process.load('PhysicsTools.PatAlgos.selectionLayer1.selectedPatCandidates_cff')
0047 patAlgosToolsTask.add(process.selectedPatCandidatesTask)
0048 process.load('PhysicsTools.PatAlgos.cleaningLayer1.cleanPatCandidates_cff')
0049 patAlgosToolsTask.add(process.cleanPatCandidatesTask)
0050 
0051 process.selectedPatMuons.cut = cms.string('muonID(\"TMOneStationTight\")'
0052     ' && abs(innerTrack.dxy) < 0.3'
0053     ' && abs(innerTrack.dz)  < 20.'
0054     ' && innerTrack.hitPattern.trackerLayersWithMeasurement > 5'
0055     ' && innerTrack.hitPattern.pixelLayersWithMeasurement > 0'
0056     ' && innerTrack.quality(\"highPurity\")'
0057 )
0058 
0059 #make patTracks
0060 from PhysicsTools.PatAlgos.tools.trackTools import makeTrackCandidates
0061 makeTrackCandidates(process,
0062     label        = 'TrackCands',                  # output collection
0063     tracks       = cms.InputTag('generalTracks'), # input track collection
0064     particleType = 'pi+',                         # particle type (for assigning a mass)
0065     preselection = 'pt > 0.7',                    # preselection cut on candidates
0066     selection    = 'pt > 0.7',                    # selection on PAT Layer 1 objects
0067     isolation    = {},                            # isolations to use (set to {} for None)
0068     isoDeposits  = [],
0069     mcAs         = None                           # replicate MC match as the one used for Muons
0070 )
0071 process.patTrackCands.embedTrack = True
0072 
0073 process.testBPHSpecificDecay = cms.EDAnalyzer('TestBPHSpecificDecay',
0074     patMuonLabel = cms.string('selectedPatMuons'),
0075     pfCandsLabel = cms.string('particleFlow'),
0076     outDump = cms.string('dump_full.txt'),
0077     outHist = cms.string('hist_full.root')
0078 )
0079 
0080 process.p = cms.Path(
0081     process.testBPHSpecificDecay,
0082     patAlgosToolsTask
0083 )
0084