def customise(process): # From python/simEmtfDigis_cfi.py if hasattr(process, 'simEmtfDigis'): process.simEmtfDigis.spPCParams16.ZoneBoundaries = [0,36,54,96,127] process.simEmtfDigis.spPCParams16.UseNewZones = True process.simEmtfDigis.RPCEnable = True process.simEmtfDigis.GEMEnable = True process.simEmtfDigis.IRPCEnable = True process.simEmtfDigis.ME0Enable = True process.simEmtfDigis.TTEnable = False process.simEmtfDigis.ME0Input = cms.InputTag('fakeSimMuonME0PadDigis') process.simEmtfDigis.Era = cms.string('Phase2C2') process.simEmtfDigis.spPAParams16.PtLUTVersion = cms.int32(7) # From python/fakeEmtfParams_cff.py if hasattr(process, 'emtfParams'): process.emtfParams.PtAssignVersion = cms.int32(7) if hasattr(process, 'emtfForestsDB'): process.emtfForestsDB = cms.ESSource( "EmptyESSource", recordName = cms.string('L1TMuonEndCapForestRcd'), iovIsRunNotTime = cms.bool(True), firstValid = cms.vuint32(1) ) process.emtfForests = cms.ESProducer( "L1TMuonEndCapForestESProducer", PtAssignVersion = cms.int32(7), bdtXMLDir = cms.string("2017_v7") ) return process
def run36xOnReRecoMC( process ): """ ------------------------------------------------------------------ running GenJets for ak5 and ak7 process : process genJets : which gen jets to run ------------------------------------------------------------------ """ print "*********************************************************************" print "NOTE TO USER: when running on 31X samples re-recoed in 3.5.6 " print " with this CMSSW version of PAT " print " it is required to re-run the GenJet production for " print " anti-kT since that is not part of the re-reco " print "*********************************************************************" process.load("RecoJets.Configuration.GenJetParticles_cff") process.load("RecoJets.JetProducers.ak5GenJets_cfi") process.ak7GenJets = process.ak5GenJets.clone( rParam = 0.7 ) process.makePatJets.replace( process.patJetCharge, process.genParticlesForJets+process.ak5GenJets+process.ak7GenJets+process.patJetCharge) #-- Remove changes for GenJets ------------------------------------------------ process.genParticlesForJets.ignoreParticleIDs = cms.vuint32(1000022, 1000012, 1000014, 1000016, 2000012, 2000014, 2000016, 1000039, 5100039, 4000012, 4000014, 4000016, 9900012, 9900014, 9900016, 39) process.genParticlesForJets.excludeResonances = True
def custom_triggercellselect_supertriggercell(process, stcSize = cms.vuint32(4,4,4) ): parameters = process.hgcalConcentratorProducer.ProcessorParameters parameters.Method = cms.string('superTriggerCellSelect') parameters.stcSize = stcSize return process
def create_supertriggercell(process, inputs, stcSize = cms.vuint32(4,4,4) ): producer = process.hgcalConcentratorProducer.clone() producer.ProcessorParameters.Method = cms.string('superTriggerCellSelect') producer.ProcessorParameters.stcSize = stcSize producer.InputTriggerCells = cms.InputTag('{}:HGCalVFEProcessorSums'.format(inputs)) producer.InputTriggerSums = cms.InputTag('{}:HGCalVFEProcessorSums'.format(inputs)) return producer
def custom_geometry_ZoltanSplit_V7(process): process.hgcalTriggerGeometryESProducer.TriggerGeometry.TriggerGeometryName = cms.string('HGCalTriggerGeometryHexLayerBasedImp1') process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TCellsMapping = cms.FileInPath("L1Trigger/L1THGCal/data/triggercell_mapping_8inch_aligned_192_432_V7_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TModulesMapping = cms.FileInPath("L1Trigger/L1THGCal/data/panel_mapping_tdr_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TCellNeighborsMapping = cms.FileInPath("L1Trigger/L1THGCal/data/triggercell_neighbor_mapping_8inch_aligned_192_432_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TCellsBHMapping = cms.FileInPath("L1Trigger/L1THGCal/data/triggercell_mapping_BH_3x3_30deg_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TCellNeighborsBHMapping = cms.FileInPath("L1Trigger/L1THGCal/data/triggercell_neighbor_mapping_BH_3x3_30deg_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.DisconnectedModules = cms.vuint32(0) process.hgcalConcentratorProducer.ProcessorParameters.MaxCellsInModule = cms.uint32(288) return process
def custom_geometry_ZoltanSplit_V8(process): process.hgcalTriggerGeometryESProducer.TriggerGeometry.TriggerGeometryName = cms.string('HGCalTriggerGeometryHexLayerBasedImp1') process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TCellsMapping = cms.FileInPath("L1Trigger/L1THGCal/data/triggercell_mapping_8inch_aligned_192_432_V8_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TModulesMapping = cms.FileInPath("L1Trigger/L1THGCal/data/panel_mapping_60deg_6mod_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TCellNeighborsMapping = cms.FileInPath("L1Trigger/L1THGCal/data/triggercell_neighbor_mapping_8inch_aligned_192_432_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TCellsBHMapping = cms.FileInPath("L1Trigger/L1THGCal/data/triggercell_mapping_BH_3x3_30deg_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TCellNeighborsBHMapping = cms.FileInPath("L1Trigger/L1THGCal/data/triggercell_neighbor_mapping_BH_3x3_30deg_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.DisconnectedModules = cms.vuint32(0) process.hgcalTriggerPrimitiveDigiProducer.FECodec.MaxCellsInModule = cms.uint32(288) process.hgcalTriggerPrimitiveDigiProducer.BEConfiguration.algorithms[0].FECodec.MaxCellsInModule = cms.uint32(288) return process
def custom_geometry_V9(process): process.hgcalTriggerGeometryESProducer.TriggerGeometry.TriggerGeometryName = cms.string('HGCalTriggerGeometryV9Imp1') process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TCellsMapping = cms.FileInPath("L1Trigger/L1THGCal/data/triggercell_mapping_8inch_aligned_192_432_V9_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TWafersMapping = cms.FileInPath("L1Trigger/L1THGCal/data/wafer_mapping_V9_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TModulesMapping = cms.FileInPath("L1Trigger/L1THGCal/data/panel_mapping_tdr_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TCellNeighborsMapping = cms.FileInPath("L1Trigger/L1THGCal/data/triggercell_neighbor_mapping_8inch_aligned_192_432_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TCellsSciMapping = cms.FileInPath("L1Trigger/L1THGCal/data/triggercell_mapping_sci_2x2_V9_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TCellNeighborsSciMapping = cms.FileInPath("L1Trigger/L1THGCal/data/triggercell_neighbor_mapping_sci_2x2_V9_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.DisconnectedModules = cms.vuint32(0) process.hgcalConcentratorProducer.ProcessorParameters.MaxCellsInModule = cms.uint32(288) return process
def addSemiLepKinFitMuon(process, isData=False) : ## std sequence to produce the kinematic fit for semi-leptonic events process.load("TopQuarkAnalysis.TopKinFitter.TtSemiLepKinFitProducer_Muons_cfi") #apply selections on muon process.cleanPatMuons.preselection = cms.string("pt>20 && abs(eta)<2.4"+ " && isGlobalMuon && isPFMuon && isTrackerMuon" + " && globalTrack.isNonnull "+ " && globalTrack.normalizedChi2<10"+ " && globalTrack.hitPattern.numberOfValidMuonHits>0"+ " && numberOfMatchedStations>1"+ " && innerTrack.hitPattern.numberOfValidPixelHits>0"+ " && track.hitPattern.trackerLayersWithMeasurement > 5"+ " && dB() < 0.2"+ " && (pfIsolationR04.sumChargedHadronPt+ max(0.,pfIsolationR04.sumNeutralHadronEt+pfIsolationR04.sumPhotonEt-0.5*pfIsolationR04.sumPUPt))/pt < 0.15" ) #clean jets from muons process.cleanPatJets.checkOverlaps.muons.requireNoOverlaps = cms.bool(True) process.cleanPatJets.preselection = cms.string("pt>15 && abs(eta)<2.5") #change constraints on kineFit process.kinFitTtSemiLepEvent.constraints = cms.vuint32(3, 4) process.kinFitTtSemiLepEvent.maxNJets = cms.int32(5) process.kinFitTtSemiLepEvent.jets = cms.InputTag("cleanPatJets") process.kinFitTtSemiLepEvent.leps = cms.InputTag("cleanPatMuons") #process.kinFitTtSemiLepEvent.mets = cms.InputTag("pfType1CorrectedMet") process.kinFitTtSemiLepEvent.mets = cms.InputTag("patMETsPF") process.kinFitTtSemiLepEvent.udscResolutions = udscResolutionPF.functions process.kinFitTtSemiLepEvent.bResolutions = bjetResolutionPF.functions process.kinFitTtSemiLepEvent.lepResolutions = muonResolution.functions process.kinFitTtSemiLepEvent.metResolutions = metResolutionPF.functions process.kinFitTtSemiLepEvent.metResolutions[0].eta = "9999" if not isData : process.kinFitTtSemiLepEvent.jetEnergyResolutionScaleFactors = cms.vdouble ( 1.052, 1.057, 1.096, 1.134, 1.288 ) process.kinFitTtSemiLepEvent.jetEnergyResolutionEtaBinning = cms.vdouble( 0.0, 0.5, 1.1, 1.7, 2.3, -1. ) # Add JES Up and Down and Rerun the KineFitter process.scaledJetEnergyUp = scaledJetEnergy.clone() process.scaledJetEnergyUp.inputJets = "cleanPatJets" process.scaledJetEnergyUp.inputMETs = "patMETsPF" process.scaledJetEnergyUp.scaleType = "jes:up" process.kinFitTtSemiLepEventJESUp = process.kinFitTtSemiLepEvent.clone() process.kinFitTtSemiLepEventJESUp.jets = cms.InputTag("scaledJetEnergyUp:cleanPatJets") process.kinFitTtSemiLepEventJESUp.mets = cms.InputTag("scaledJetEnergyUp:patMETsPF") process.scaledJetEnergyDown = scaledJetEnergy.clone() process.scaledJetEnergyDown.inputJets = "cleanPatJets" process.scaledJetEnergyDown.inputMETs = "patMETsPF" process.scaledJetEnergyDown.scaleType = "jes:down" process.kinFitTtSemiLepEventJESDown = process.kinFitTtSemiLepEvent.clone() process.kinFitTtSemiLepEventJESDown.jets = cms.InputTag("scaledJetEnergyDown:cleanPatJets") process.kinFitTtSemiLepEventJESDown.mets = cms.InputTag("scaledJetEnergyDown:patMETsPF") process.kinFitSequence = cms.Sequence(process.kinFitTtSemiLepEvent * process.scaledJetEnergyUp * process.kinFitTtSemiLepEventJESUp * process.scaledJetEnergyDown * process.kinFitTtSemiLepEventJESDown)
def custom_geometry_V9(process): process.hgcalTriggerGeometryESProducer.TriggerGeometry.TriggerGeometryName = cms.string('HGCalTriggerGeometryV9Imp1') process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TCellsMapping = cms.FileInPath("L1Trigger/L1THGCal/data/triggercell_mapping_8inch_aligned_192_432_V9_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TWafersMapping = cms.FileInPath("L1Trigger/L1THGCal/data/wafer_mapping_V9_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TModulesMapping = cms.FileInPath("L1Trigger/L1THGCal/data/panel_mapping_tdr_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TCellNeighborsMapping = cms.FileInPath("L1Trigger/L1THGCal/data/triggercell_neighbor_mapping_8inch_aligned_192_432_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TCellsSciMapping = cms.FileInPath("L1Trigger/L1THGCal/data/triggercell_mapping_sci_2x2_V9_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.L1TCellNeighborsSciMapping = cms.FileInPath("L1Trigger/L1THGCal/data/triggercell_neighbor_mapping_sci_2x2_V9_0.txt") process.hgcalTriggerGeometryESProducer.TriggerGeometry.DisconnectedModules = cms.vuint32(0) process.hgcalTriggerPrimitiveDigiProducer.FECodec.MaxCellsInModule = cms.uint32(288) process.hgcalTriggerPrimitiveDigiProducer.BEConfiguration.algorithms[0].FECodec.MaxCellsInModule = cms.uint32(288) return process
def add_tau_variables(process, sequence, pfjet_tag, patjet_tag, algo, taus): tag = "taus" + patjet_tag.replace("patJets", "") tau_calculator = Njettiness.clone( src=cms.InputTag(pfjet_tag), Njets=cms.vuint32(taus), # variables for measure definition: measureDefinition = cms.uint32(0), # CMS default is normalized measure beta=cms.double(1.0), # CMS default is 1 R0=cms.double(algo.r), # CMS default is jet cone size Rcutoff=cms.double(999.0), # not used by default # variables for axes definition: axesDefinition=cms.uint32(6), # CMS default is 1-pass KT axes nPass=cms.int32(999), # not used by default akAxesR0=cms.double(-999.0), # not used by default ) setattr(process, tag, tau_calculator) for tau in taus: getattr(process, patjet_tag).userData.userFloats.src += ['{}:tau{}'.format(tag, tau)] sequence += getattr(process, tag) return tag
import FWCore.ParameterSet.Config as cms EventNumberArrayFilter = cms.EDFilter( "EventNumberArrayFilter", # select event numbers selectedEventNumbers=cms.vuint32())
#define a parameter set to be passed to all modules that utilize GenTauDecayID for signal taus commonGenTauDecayIDPSet = cms.PSet(momPDGID=cms.vint32(A_PDGID), chargedHadronPTMin=cms.double(0.0), neutralHadronPTMin=cms.double(0.0), chargedLeptonPTMin=cms.double(0.0), totalPTMin=cms.double(0.0)) #define a parameter set for the W-->munu selector WMuNuPSet = commonGenTauDecayIDPSet.clone() WMuNuPSet.momPDGID = cms.vint32(W_PDGID) #only proceed if event is a true W-->munu event process.genWMuNuSelector = cms.EDFilter( 'GenObjectProducer', genParticleTag=cms.InputTag('genParticles'), absMatchPDGIDs=cms.vuint32(MU_PDGID), sisterAbsMatchPDGID=cms.uint32(NUMU_PDGID), genTauDecayIDPSet=WMuNuPSet, primaryTauDecayType=cms.uint32(TAU_ALL), sisterTauDecayType=cms.uint32(TAU_ALL), primaryTauPTRank=cms.int32(ANY_PT_RANK), primaryTauHadronicDecayType=cms.int32(TAU_ALL_HAD), sisterHadronicDecayType=cms.int32(TAU_ALL_HAD), primaryTauAbsEtaMax=cms.double(-1.0), primaryTauPTMin=cms.double(-1.0), countSister=cms.bool(False), applyPTCuts=cms.bool(False), countKShort=cms.bool(True), minNumGenObjectsToPassFilter=cms.uint32(0), makeAllCollections=cms.bool(False))
#define VertexValidationVertices for the vertex DQM validation process.VertexValidationVertices= cms.EDAnalyzer("HLTVertexPerformanceAnalyzer", SimVertexCollection = cms.InputTag("g4SimHits"), TriggerResults = cms.InputTag('TriggerResults','',fileini.processname), HLTPathNames = cms.vstring(fileini.vertex_pathes), Vertex = fileini.vertex_modules, ) #define bTagValidation for the b-tag DQM validation (distribution plot) process.bTagValidation = cms.EDAnalyzer("HLTBTagPerformanceAnalyzer", TriggerResults = cms.InputTag('TriggerResults','',fileini.processname), HLTPathNames = cms.vstring(fileini.btag_pathes), JetTag = fileini.btag_modules, MinJetPT = cms.double(20), mcFlavours = cms.PSet( light = cms.vuint32(1, 2, 3, 21), # udsg c = cms.vuint32(4), b = cms.vuint32(5), g = cms.vuint32(21), uds = cms.vuint32(1, 2, 3) ), mcPartons = cms.InputTag("hltBtagJetsbyValAlgo") ) #define bTagPostValidation for the b-tag DQM validation (efficiency and mistagrate plot) process.bTagPostValidation = cms.EDAnalyzer("HLTBTagHarvestingAnalyzer", HLTPathNames = fileini.btag_pathes, histoName = fileini.btag_modules_string, minTag = cms.double(0.6), # MC stuff mcFlavours = cms.PSet(
# import of standard configurations process.load('Configuration/StandardSequences/Services_cff') process.load('FWCore/MessageService/MessageLogger_cfi') process.load('Configuration/EventContent/EventContent_cff') # module to select events based on HLT Trigger bits process.load('HLTrigger/HLTfilters/hltHighLevelDev_cfi') # Loading "hltHighLevelDev_cfi" defines an EDFilter called hltHighLevelDev # now we can configure it process.hltHighLevelDev.TriggerResultsTag = cms.InputTag("TriggerResults","","HLT8E29") process.hltHighLevelDev.HLTPaths = ("HLT_Mu9",) process.hltHighLevelDev.andOr = True # True = OR, False = AND process.hltHighLevelDev.HLTPathsPrescales = cms.vuint32(1,) process.hltHighLevelDev.HLTOverallPrescale = cms.uint32 (20) # Once the hltHighLevelDev filter has been defined I can clone it and configure as needed process.hltHighLevelDev2 = process.hltHighLevelDev.clone(andOr = True) process.hltHighLevelDev2.HLTPaths = ("HLT_L1Mu20","HLT_L2Mu9","HLT_L2Mu11","HLT_L1Mu14_L1SingleEG10","HLT_L1Mu14_L1SingleJet6U","HLT_L1Mu14_L1ETM30") process.hltHighLevelDev2.andOr = True # True = OR, False = AND process.hltHighLevelDev2.HLTPathsPrescales = cms.vuint32(1,1,1,1,1,1) process.hltHighLevelDev2.HLTOverallPrescale = cms.uint32 (20) ############# pt_hat Filter ##### #process.ptHat_filter = cms.EDFilter("MCProcessFilter", # ProcessID = cms.untracked.vint32(0), # MinPthat = cms.untracked.vdouble(50), # MaxPthat = cms.untracked.vdouble(150.)
dataTier = cms.untracked.string('RAW-RECO'), filterName = cms.untracked.string('MuonDPG_skim')), SelectEvents = cms.untracked.PSet( SelectEvents = cms.vstring("l1MuBitsSkim","dtHLTSkim","dtSkim","cscHLTSkim","cscSkimAlone","rpcTecSkim","muonTracksSkim") ) ) #################################################################################### #################################logerrorharvester############################################ process.load("FWCore.Modules.logErrorFilter_cfi") from Configuration.StandardSequences.RawToDigi_Data_cff import gtEvmDigis process.gtEvmDigis = gtEvmDigis.clone() process.stableBeam = cms.EDFilter("HLTBeamModeFilter", L1GtEvmReadoutRecordTag = cms.InputTag("gtEvmDigis"), AllowedBeamMode = cms.vuint32(11), saveTags = cms.bool(False) ) process.logerrorpath=cms.Path(process.gtEvmDigis+process.stableBeam+process.logErrorFilter) process.outlogerr = cms.OutputModule("PoolOutputModule", outputCommands = process.FEVTEventContent.outputCommands, fileName = cms.untracked.string('/tmp/azzi/logerror_filter.root'), dataset = cms.untracked.PSet( dataTier = cms.untracked.string('RAW-RECO'), filterName = cms.untracked.string('Skim_logerror')), SelectEvents = cms.untracked.PSet( SelectEvents = cms.vstring("logerrorpath") ))
process.load('FWCore/MessageService/MessageLogger_cfi') process.load('Configuration/EventContent/EventContent_cff') process.load('JetMETCorrections.Configuration.jecHLTFilters_cfi') # module to select events based on HLT Trigger bits process.load('HLTrigger/HLTfilters/hltHighLevelDev_cfi') process.load('Configuration/StandardSequences/MagneticField_38T_cff') # Loading "hltHighLevelDev_cfi" defines an EDFilter called hltHighLevelDev # now we can configure it process.hltHighLevelDev.TriggerResultsTag = cms.InputTag("TriggerResults","","HLT8E29") process.hltHighLevelDev.HLTPaths = ("HLT_Jet30U","HLT_QuadJet15U","HLT_FwdJet20U") process.hltHighLevelDev.andOr = True # True = OR, False = AND process.hltHighLevelDev.HLTPathsPrescales = cms.vuint32(10,1,10) process.hltHighLevelDev.HLTOverallPrescale =cms.uint32(1) # Once the hltHighLevelDev filter has been defined I can clone it and configure as needed process.hltHighLevelDev2 = process.hltHighLevelDev.clone(andOr = True) process.hltHighLevelDev2.HLTPaths = ("HLT_Jet50U",) process.hltHighLevelDev2.andOr = True # True = OR, False = AND process.hltHighLevelDev2.HLTPathsPrescales = cms.vuint32(1,) process.hltHighLevelDev2.HLTOverallPrescale = cms.uint32(1) # Once the hltHighLevelDev filter has been defined I can clone it and configure as needed process.hltHighLevelDev3 = process.hltHighLevelDev.clone(andOr = True) process.hltHighLevelDev3.HLTPaths = ("HLT_DiJetAve30U_8E29",) process.hltHighLevelDev3.andOr = True # True = OR, False = AND process.hltHighLevelDev3.HLTPathsPrescales = cms.vuint32(1,) process.hltHighLevelDev3.HLTOverallPrescale = cms.uint32(1)
def customise_gpu_ecal(process): # FIXME replace the Sequences with empty ones to avoid exanding them during the (re)definition of Modules and EDAliases process.HLTDoFullUnpackingEgammaEcalMFSequence = cms.Sequence() process.HLTDoFullUnpackingEgammaEcalWithoutPreshowerSequence = cms.Sequence( ) process.HLTDoFullUnpackingEgammaEcalSequence = cms.Sequence() # Event Setup process.load( "EventFilter.EcalRawToDigi.ecalElectronicsMappingGPUESProducer_cfi") process.load( "RecoLocalCalo.EcalRecProducers.ecalGainRatiosGPUESProducer_cfi") process.load( "RecoLocalCalo.EcalRecProducers.ecalPedestalsGPUESProducer_cfi") process.load( "RecoLocalCalo.EcalRecProducers.ecalPulseCovariancesGPUESProducer_cfi") process.load( "RecoLocalCalo.EcalRecProducers.ecalPulseShapesGPUESProducer_cfi") process.load( "RecoLocalCalo.EcalRecProducers.ecalSamplesCorrelationGPUESProducer_cfi" ) process.load( "RecoLocalCalo.EcalRecProducers.ecalTimeBiasCorrectionsGPUESProducer_cfi" ) process.load( "RecoLocalCalo.EcalRecProducers.ecalTimeCalibConstantsGPUESProducer_cfi" ) # Modules and EDAliases process.hltEcalDigisGPU = cms.EDProducer( "EcalRawToDigiGPU", InputLabel=cms.InputTag("rawDataCollector"), FEDs=cms.vint32(601, 602, 603, 604, 605, 606, 607, 608, 609, 610, 611, 612, 613, 614, 615, 616, 617, 618, 619, 620, 621, 622, 623, 624, 625, 626, 627, 628, 629, 630, 631, 632, 633, 634, 635, 636, 637, 638, 639, 640, 641, 642, 643, 644, 645, 646, 647, 648, 649, 650, 651, 652, 653, 654), digisLabelEB=cms.string('ebDigis'), digisLabelEE=cms.string('eeDigis'), maxChannels=cms.uint32(4 * 20000)) process.hltEcalDigis = cms.EDProducer( "EcalCPUDigisProducer", digisInLabelEB=cms.InputTag("hltEcalDigisGPU", "ebDigis"), digisInLabelEE=cms.InputTag("hltEcalDigisGPU", "eeDigis"), digisOutLabelEB=cms.string('ebDigis'), digisOutLabelEE=cms.string('eeDigis'), produceDummyIntegrityCollections=cms.bool(True)) process.hltEcalUncalibRecHitGPU = cms.EDProducer( "EcalUncalibRecHitProducerGPU", digisLabelEB=cms.InputTag("hltEcalDigisGPU", "ebDigis"), digisLabelEE=cms.InputTag("hltEcalDigisGPU", "eeDigis"), recHitsLabelEB=cms.string("EcalUncalibRecHitsEB"), recHitsLabelEE=cms.string("EcalUncalibRecHitsEE"), EBamplitudeFitParameters=cms.vdouble(1.138, 1.652), EBtimeConstantTerm=cms.double(0.6), EBtimeFitLimits_Lower=cms.double(0.2), EBtimeFitLimits_Upper=cms.double(1.4), EBtimeFitParameters=cms.vdouble(-2.015452, 3.130702, -12.3473, 41.88921, -82.83944, 91.01147, -50.35761, 11.05621), EBtimeNconst=cms.double(28.5), EEamplitudeFitParameters=cms.vdouble(1.89, 1.4), EEtimeConstantTerm=cms.double(1.0), EEtimeFitLimits_Lower=cms.double(0.2), EEtimeFitLimits_Upper=cms.double(1.4), EEtimeFitParameters=cms.vdouble(-2.390548, 3.553628, -17.62341, 67.67538, -133.213, 140.7432, -75.41106, 16.20277), EEtimeNconst=cms.double(31.8), amplitudeThresholdEB=cms.double(10.0), amplitudeThresholdEE=cms.double(10.0), outOfTimeThresholdGain12mEB=cms.double(5.0), outOfTimeThresholdGain12mEE=cms.double(1000.0), outOfTimeThresholdGain12pEB=cms.double(5.0), outOfTimeThresholdGain12pEE=cms.double(1000.0), outOfTimeThresholdGain61mEB=cms.double(5.0), outOfTimeThresholdGain61mEE=cms.double(1000.0), outOfTimeThresholdGain61pEB=cms.double(5.0), outOfTimeThresholdGain61pEE=cms.double(1000.0), kernelMinimizeThreads=cms.vuint32(32, 1, 1), maxNumberHits=cms.uint32(4 * 20000), shouldRunTimingComputation=cms.bool(False)) process.hltEcalUncalibRecHitSoA = cms.EDProducer( "EcalCPUUncalibRecHitProducer", containsTimingInformation=cms.bool(False), recHitsInLabelEB=cms.InputTag("hltEcalUncalibRecHitGPU", "EcalUncalibRecHitsEB"), recHitsInLabelEE=cms.InputTag("hltEcalUncalibRecHitGPU", "EcalUncalibRecHitsEE"), recHitsOutLabelEB=cms.string('EcalUncalibRecHitsEB'), recHitsOutLabelEE=cms.string('EcalUncalibRecHitsEE')) process.hltEcalUncalibRecHit = cms.EDProducer( "EcalUncalibRecHitConvertGPU2CPUFormat", recHitsLabelGPUEB=cms.InputTag("hltEcalUncalibRecHitSoA", "EcalUncalibRecHitsEB"), recHitsLabelGPUEE=cms.InputTag("hltEcalUncalibRecHitSoA", "EcalUncalibRecHitsEE"), recHitsLabelCPUEB=cms.string("EcalUncalibRecHitsEB"), recHitsLabelCPUEE=cms.string("EcalUncalibRecHitsEE")) # Sequences process.HLTDoFullUnpackingEgammaEcalMFSequence = cms.Sequence( process.hltEcalDigisGPU + process.hltEcalDigis + process.hltEcalPreshowerDigis + process.hltEcalUncalibRecHitGPU + process.hltEcalUncalibRecHitSoA + process.hltEcalUncalibRecHit + process.hltEcalDetIdToBeRecovered + process.hltEcalRecHit + process.hltEcalPreshowerRecHit) process.HLTDoFullUnpackingEgammaEcalWithoutPreshowerSequence = cms.Sequence( process.hltEcalDigisGPU + process.hltEcalDigis + process.hltEcalUncalibRecHitGPU + process.hltEcalUncalibRecHitSoA + process.hltEcalUncalibRecHit + process.hltEcalDetIdToBeRecovered + process.hltEcalRecHit) process.HLTDoFullUnpackingEgammaEcalSequence = cms.Sequence( process.hltEcalDigisGPU + process.hltEcalDigis + process.hltEcalPreshowerDigis + process.hltEcalUncalibRecHitGPU + process.hltEcalUncalibRecHitSoA + process.hltEcalUncalibRecHit + process.hltEcalDetIdToBeRecovered + process.hltEcalRecHit + process.hltEcalPreshowerRecHit) # done return process
import FWCore.ParameterSet.Config as cms hltHcalSimpleRecHitFilter = cms.EDFilter( 'HLTHcalSimpleRecHitFilter', saveTags=cms.bool(True), HFRecHitCollection=cms.InputTag('hltHfreco'), threshold=cms.double(3), minNHitsNeg=cms.int32(1), minNHitsPos=cms.int32(1), doCoincidence=cms.bool(True), maskedChannels=cms.vuint32())
process.load('SimCalorimetry.HcalTrigPrimProducers.hcaltpdigi_cff' ) # needed to read HCal TPs process.load('CalibCalorimetry.CaloTPG.CaloTPGTranscoder_cfi') process.load('Configuration.StandardSequences.SimL1Emulator_cff') process.load('L1Trigger.TrackTrigger.TrackTrigger_cff') process.load("L1Trigger.TrackFindingTracklet.L1HybridEmulationTracks_cff") process.load("L1Trigger.TrackerDTC.ProducerES_cff") process.load("L1Trigger.TrackerDTC.ProducerED_cff") process.load("RecoVertex.BeamSpotProducer.BeamSpot_cfi") process.l1ctLayer1Barrel9 = process.l1ctLayer1Barrel.clone() process.l1ctLayer1Barrel9.puAlgo.nFinalSort = 32 process.l1ctLayer1Barrel9.regions[0].etaBoundaries = [-1.5, -0.5, 0.5, 1.5] process.l1ctLayer1Barrel9.boards = cms.VPSet( cms.PSet(regions=cms.vuint32( *[0 + 9 * ie + i for ie in range(3) for i in range(3)])), cms.PSet(regions=cms.vuint32( *[3 + 9 * ie + i for ie in range(3) for i in range(3)])), cms.PSet(regions=cms.vuint32( *[6 + 9 * ie + i for ie in range(3) for i in range(3)])), ) process.pfInputsTask = cms.Task(process.TTClustersFromPhase2TrackerDigis, process.TTStubsFromPhase2TrackerDigis, process.TrackerDTCProducer, process.offlineBeamSpot, process.TTTracksFromTrackletEmulation, process.SimL1EmulatorTask) process.runPF = cms.Path(process.l1ctLayer1Barrel + process.l1ctLayer1Barrel9 + process.l1ctLayer1HGCal + process.l1ctLayer1HGCalNoTK + process.l1ctLayer1HF)
#no consideration of pT rank ANY_PT_RANK = -1 #define a parameter set to be passed to all modules that utilize GenTauDecayID commonGenTauDecayIDPSet = cms.PSet(momPDGID=cms.vint32(A_PDGID), chargedHadronPTMin=cms.double(0.0), neutralHadronPTMin=cms.double(0.0), chargedLeptonPTMin=cms.double(0.0), totalPTMin=cms.double(0.0)) #di-mu+X filter process.genDiTauMuTauXSelector = cms.EDFilter( 'GenObjectProducer', genParticleTag=cms.InputTag('genParticles'), absMatchPDGIDs=cms.vuint32(TAU_PDGID), sisterAbsMatchPDGID=cms.uint32(TAU_PDGID), genTauDecayIDPSet=commonGenTauDecayIDPSet, primaryTauDecayType=cms.uint32(TAU_MU), sisterTauDecayType=cms.uint32(TAU_ALL), primaryTauPTRank=cms.int32(ANY_PT_RANK), primaryTauHadronicDecayType=cms.int32(TAU_ALL_HAD), sisterHadronicDecayType=cms.int32(TAU_ALL_HAD), primaryTauAbsEtaMax=cms.double(-1.0), primaryTauPTMin=cms.double(-1.0), countSister=cms.bool(False), applyPTCuts=cms.bool(False), countKShort=cms.bool(True), minNumGenObjectsToPassFilter=cms.uint32(2), makeAllCollections=cms.bool(False))
# this module process.load('L1Trigger.GlobalTriggerAnalyzer.l1GtBeamModeFilter_cfi') # replacing arguments for l1GtBeamModeFilter # input tag for input tag for ConditionInEdm products #process.l1GtBeamModeFilter.CondInEdmInputTag = cms.InputTag("conditionsInEdm"), # input tag for the L1 GT EVM product #process.l1GtBeamModeFilter.L1GtEvmReadoutRecordTag = cms.InputTag("gtEvmDigis") # vector of allowed beam modes (see enumeration in header file for implemented values) # default value: 11 (STABLE) #process.l1GtBeamModeFilter.AllowedBeamMode = cms.vuint32(11) process.l1GtBeamModeFilter.AllowedBeamMode = cms.vuint32(9, 10, 11) # return the inverted result, to be used instead of NOT #process.l1GtBeamModeFilter.InvertResult = False # path to be run # for RAW data, run first the RAWTODIGI if (dataType == 'RAW') and not (useRelValSample): process.load('Configuration/StandardSequences/RawToDigi_Data_cff') process.gtEvmDigis.UnpackBxInEvent = cms.int32(1) # set EVM unpacker to verbose process.gtEvmDigis.Verbosity = cms.untracked.int32(1) process.p = cms.Path(process.RawToDigi * process.l1GtBeamModeFilter) elif (dataType == 'RAW') and (useRelValSample):
process.load("Configuration.TotemCommon.RandomNumbers_cfi") process.load("TotemAlignment.T2TrkBasedInternalAlignment.T2TrkIntAlign_cfi") process.T2TrkBasedIntAl.OutputFile=cms.untracked.string('50000IPMuons_T2DigisH07x7yDispl_H17x7yDispl_PhiFlatHighEta_conf5_0AndShiftAndTilt_testMille.root') process.T2TrkBasedIntAl.simufile=cms.bool(True) process.T2TrkBasedIntAl.DoALign=cms.bool(True) process.T2TrkBasedIntAl.UseMillePede=cms.bool(True) #Two reference detector per qurter is assumed in the Millepede algorithm process.T2TrkBasedIntAl.IdreferencedetMille=cms.vuint32(1,9,11,19,21,29,31,39) process.T2TrkBasedIntAl.XreferencedetMille=cms.vdouble(0.,0.,0.,0.,0.,0.,0.,0.) process.T2TrkBasedIntAl.YreferencedetMille=cms.vdouble(0.,0.,0.,0.,0.,0.,0.,0.) #process.T2TrkBasedIntAl.DetEffRWind =cms.double(20.0) #process.T2TrkBasedIntAl.DetEffPhiWind =10.0 #process.T2TrkBasedIntAl.Numhittrkeff =4 # num cl1 hit in other detectors for efficiency calculation #process.T2TrkBasedIntAl.Testedcamera = cms.vuint32(0,2,4,6,8) process.T2TrkBasedIntAl.HitNumb4Align=10 process.T2TrkBasedIntAl.MeasuredResol_X=0.01 #Stop Condition internal alignment process.T2TrkBasedIntAl.MeasuredResol_Y=0.01 process.T2TrkBasedIntAl.SHIFTprescale=1.0 process.T2TrkBasedIntAl.MaxStepalignstep=50 process.T2TrkBasedIntAl.Idreferencedet=cms.vuint32(1,11) #force fixed ref det for alignment process.T2TrkBasedIntAl.MaxEvents=49900 process.T2TrkBasedIntAl.ratioXYAlignTrkcut=0.3
# Filter setup # ------------------------ # https://github.com/cms-sw/cmssw/blob/CMSSW_8_0_X/PhysicsTools/HepMCCandAlgos/python/genParticles_cfi.py tmpGenParticles = cms.EDProducer( "GenParticleProducer", saveBarCodes=cms.untracked.bool(True), src=cms.InputTag("generator", "unsmeared"), abortOnUnknownPDGCode=cms.untracked.bool(False)) # https://github.com/cms-sw/cmssw/blob/CMSSW_8_0_X/RecoJets/Configuration/python/GenJetParticles_cff.py # https://github.com/cms-sw/cmssw/blob/CMSSW_8_0_X/RecoMET/Configuration/python/GenMETParticles_cff.py tmpGenParticlesForJetsNoNu = cms.EDProducer( "InputGenJetsParticleSelector", src=cms.InputTag("tmpGenParticles"), ignoreParticleIDs=cms.vuint32(1000022, 1000012, 1000014, 1000016, 2000012, 2000014, 2000016, 1000039, 5100039, 4000012, 4000014, 4000016, 9900012, 9900014, 9900016, 39, 12, 14, 16), partonicFinalState=cms.bool(False), excludeResonances=cms.bool(False), excludeFromResonancePids=cms.vuint32(12, 13, 14, 16), tausAsJets=cms.bool(False)) # https://github.com/cms-sw/cmssw/blob/CMSSW_8_0_X/RecoJets/JetProducers/python/AnomalousCellParameters_cfi.py AnomalousCellParameters = cms.PSet(maxBadEcalCells=cms.uint32(9999999), maxRecoveredEcalCells=cms.uint32(9999999), maxProblematicEcalCells=cms.uint32(9999999), maxBadHcalCells=cms.uint32(9999999), maxRecoveredHcalCells=cms.uint32(9999999), maxProblematicHcalCells=cms.uint32(9999999)) # https://github.com/cms-sw/cmssw/blob/CMSSW_8_0_X/RecoJets/JetProducers/python/GenJetParameters_cfi.py
import FWCore.ParameterSet.Config as cms Njettiness = cms.EDProducer( "NjettinessAdder", src=cms.InputTag("ak8PFJetsCHS"), Njets=cms.vuint32(1, 2, 3, 4), # compute 1-, 2-, 3-, 4- subjettiness # variables for measure definition : measureDefinition=cms.uint32(0), # CMS default is normalized measure beta=cms.double(1.0), # CMS default is 1 R0=cms.double(0.8), # CMS default is jet cone size Rcutoff=cms.double(-999.0), # not used by default # variables for axes definition : axesDefinition=cms.uint32(6), # CMS default is 1-pass KT axes nPass=cms.int32(-999), # not used by default akAxesR0=cms.double(-999.0) # not used by default )
process.btaggingAK5PF = cms.Sequence( process.impactParameterTagInfosAK5PF +process.secondaryVertexTagInfosAK5PF +process.simpleSecondaryVertexHighEffBJetTagsAK5PF +process.simpleSecondaryVertexHighPurBJetTagsAK5PF +process.combinedSecondaryVertexBJetTagsAK5PF ) ################################################################ ### MC specific config ################################################################ if not isData: process.load("RecoJets.Configuration.GenJetParticles_cff") process.genParticlesForJets.ignoreParticleIDs = cms.vuint32( 1000022, 2000012, 2000014, 2000016, 1000039, 5000039, 4000012, 9900012, 9900014, 9900016, 39, 12, 14, 16 ) process.load("RecoJets.JetProducers.ak5GenJets_cfi") process.ak5GenJetsNoNuBSM = process.ak5GenJets.clone(); ## Configure Parton Jets #process.genParticlesForPartonJets = process.genParticlesForJets.clone() #process.genParticlesForPartonJets.partonicFinalState = True #process.genParticlesForPartonJets.excludeFromResonancePids = cms.vuint32(11, 12, 13, 14, 15, 16) ################################################################ ### General Config ################################################################ process.MessageLogger.cerr.FwkReport.reportEvery = 10000 process.MessageLogger.suppressError = cms.untracked.vstring( 'patTrigger','HLTConfigData' )
cms.PSet( importerName = cms.string("SuperClusterImporter"), source_eb = cms.InputTag("particleFlowSuperClusterECAL:particleFlowSuperClusterECALBarrel"), source_ee = cms.InputTag("particleFlowSuperClusterECAL:particleFlowSuperClusterECALEndcapWithPreshower"), source_towers = cms.InputTag("towerMaker"), maximumHoverE = cms.double(0.5), minSuperClusterPt = cms.double(10.0), minPTforBypass = cms.double(100.0), superClustersArePF = cms.bool(True) ), # all secondary track importers cms.PSet( importerName = cms.string("GeneralTracksImporter"), source = cms.InputTag("pfTrack"), muonSrc = cms.InputTag("hiMuons1stStep"), useIterativeTracking = cms.bool(False), DPtOverPtCuts_byTrackAlgo = cms.vdouble(-1.0,-1.0,-1.0, 1.0,1.0), NHitCuts_byTrackAlgo = cms.vuint32(3,3,3,3,3) ), # to properly set SC based links you need to run ECAL importer # after you've imported all SCs to the block cms.PSet( importerName = cms.string("ECALClusterImporter"), source = cms.InputTag("particleFlowClusterECAL"), BCtoPFCMap = cms.InputTag('particleFlowSuperClusterECAL:PFClusterAssociationEBEE') ), cms.PSet( importerName = cms.string("GenericClusterImporter"), source = cms.InputTag("particleFlowClusterHCAL") ), cms.PSet( importerName = cms.string("GenericClusterImporter"), source = cms.InputTag("particleFlowClusterHO") ), cms.PSet( importerName = cms.string("GenericClusterImporter"), source = cms.InputTag("particleFlowClusterHF") ), cms.PSet( importerName = cms.string("GenericClusterImporter"), source = cms.InputTag("particleFlowClusterPS") ) )
errortype = cms.string('whole'), detid = cms.uint32(302123024) ), #Lay3 # BPix_BpI_SEC1_LYR3_LDR3_MOD1 cms.PSet( errortype = cms.string('whole'), detid = cms.uint32(302188820) ), #FPix_BpO_D1_BLD5_PNL2_PLQ2_ROC6 (352398856) cms.PSet( errortype = cms.string('none'), detid = cms.uint32(352398856), badroclist = cms.vuint32(6) ), #--------------------------------------------------------------------------- ), Record = cms.string('SiPixelQualityFromDbRcd'), SinceAppendMode = cms.bool(True), IOVMode = cms.string('Run'), printDebug = cms.untracked.bool(True), doStoreOnDB = cms.bool(True) ) #process.print = cms.OutputModule("AsciiOutputModule") process.p = cms.Path(process.prod) #process.ep = cms.EndPath(process.print)
process.schedule = cms.Schedule(process.patPath, process.fakPath, process.countOverlaps, process.scrap, process.outpath) elif doFakeRates == 'only': process.out.SelectEvents = cms.untracked.PSet( SelectEvents=cms.vstring('fakPath')) process.schedule = cms.Schedule(process.fakPath, process.scrap, process.outpath) if doNoFilter: ##Extra GenJets process.load('RecoJets.Configuration.RecoGenJets_cff') process.load('RecoJets.Configuration.GenJetParticles_cff') process.genParticlesForJetsNoElNoMuNoNu = process.genParticlesForJetsNoMuNoNu.clone( ) process.genParticlesForJetsNoElNoMuNoNu.ignoreParticleIDs += cms.vuint32( 11) process.ak5GenJetsNoElNoMuNoNu = process.ak5GenJetsNoMuNoNu.clone() process.ak5GenJetsNoElNoMuNoNu.src = cms.InputTag( "genParticlesForJetsNoElNoMuNoNu") process.extraGenJets = cms.Sequence( process.genParticlesForJetsNoNu + process.ak5GenJetsNoNu + process.genParticlesForJetsNoMuNoNu + process.ak5GenJetsNoMuNoNu + process.genParticlesForJetsNoElNoMuNoNu + process.ak5GenJetsNoElNoMuNoNu) process.allPath += process.extraGenJets #No need to run patPath (with preYieldFilter) process.out.SelectEvents = cms.untracked.PSet( SelectEvents=cms.vstring('allPath')) process.schedule = cms.Schedule(process.allPath, process.scrap, process.outpath)
fileNames = cms.untracked.vstring( '/store/relval/CMSSW_3_1_2/RelValQCD_Pt_80_120/GEN-SIM-RECO/STARTUP31X_V2-v1/0007/DCAE40E8-CA78-DE11-8F20-001D09F2305C.root', '/store/relval/CMSSW_3_1_2/RelValQCD_Pt_80_120/GEN-SIM-RECO/STARTUP31X_V2-v1/0007/D099AB85-CA78-DE11-9A5E-001D09F2503C.root', '/store/relval/CMSSW_3_1_2/RelValQCD_Pt_80_120/GEN-SIM-RECO/STARTUP31X_V2-v1/0007/B2E190DF-CA78-DE11-9F35-001D09F2532F.root', '/store/relval/CMSSW_3_1_2/RelValQCD_Pt_80_120/GEN-SIM-RECO/STARTUP31X_V2-v1/0007/90A07B7B-CA78-DE11-A4B7-001D09F26C5C.root', '/store/relval/CMSSW_3_1_2/RelValQCD_Pt_80_120/GEN-SIM-RECO/STARTUP31X_V2-v1/0007/4C5F86D1-CB78-DE11-8650-000423D6B42C.root', '/store/relval/CMSSW_3_1_2/RelValQCD_Pt_80_120/GEN-SIM-RECO/STARTUP31X_V2-v1/0007/2247873A-B378-DE11-8F5B-001D09F24664.root', ), ) process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(10) ) # Identify GenParticles to be used to build GenJets (ie, no neutrinos or BSM) process.load("RecoJets.Configuration.GenJetParticles_cff") process.genParticlesForJets.ignoreParticleIDs = cms.vuint32( 1000022, 2000012, 2000014, 2000016, 1000039, 5000039, 4000012, 9900012, 9900014, 9900016, 39, 12, 14, 16 ) process.genParticlesForJets.excludeFromResonancePids = cms.vuint32(12, 14, 16) # Build reco::GenJets from GenParticles from RecoJets.JetProducers.iterativeCone5GenJets_cff import iterativeCone5GenJets process.iterativeCone5GenJetsNoNuBSM = iterativeCone5GenJets.clone() # Jet-track association from RecoJets.JetAssociationProducers.iterativeCone5JTA_cff import* ZSPiterativeCone5JetTracksAssociatorAtVertex = iterativeCone5JetTracksAssociatorAtVertex.clone() ZSPiterativeCone5JetTracksAssociatorAtVertex.jets = cms.InputTag("ZSPJetCorJetIcone5") ZSPiterativeCone5JetTracksAssociatorAtCaloFace = iterativeCone5JetTracksAssociatorAtCaloFace.clone() ZSPiterativeCone5JetTracksAssociatorAtCaloFace.jets = cms.InputTag("ZSPJetCorJetIcone5") ZSPiterativeCone5JetExtender = iterativeCone5JetExtender.clone()
options.outputFile = 'this_is_not_used' options.maxEvents = 10 # parse command-line options options.parseArguments() # minimum of logs process.MessageLogger = cms.Service( "MessageLogger", statistics=cms.untracked.vstring(), destinations=cms.untracked.vstring('cerr'), cerr=cms.untracked.PSet(threshold=cms.untracked.string('WARNING'))) # raw data source process.source = cms.Source("PoolSource", fileNames=cms.untracked.vstring( options.inputFiles)) process.maxEvents = cms.untracked.PSet( input=cms.untracked.int32(options.maxEvents)) # frame analyzer process.totemVFATFrameAnalyzer = cms.EDAnalyzer( "TotemVFATFrameAnalyzer", rawDataTag=cms.InputTag("rawDataCollector"), fedIds=cms.vuint32(578, 579, 580), RawUnpacking=cms.PSet()) # execution configuration process.p = cms.Path(process.totemVFATFrameAnalyzer)
import FWCore.ParameterSet.Config as cms singleLepCalc = cms.PSet( dataType = cms.string('All'), isMc = cms.bool(True), pvCollection = cms.InputTag("offlineSlimmedPrimaryVertices"), genParticles = cms.InputTag("prunedGenParticles"), genJets_it = cms.InputTag("slimmedGenJets"), triggerCollection = cms.InputTag("TriggerResults::HLT"), triggerSummary = cms.InputTag("selectedPatTrigger"), keepFullMChistory = cms.bool(True), keepPDGID = cms.vuint32(1, 2, 3, 4, 5, 21, 11, 12, 13, 14, 15, 16), keepMomPDGID = cms.vuint32(6, 24), rhoSrc = cms.InputTag("fixedGridRhoAll"), isWJets = cms.bool(False), isTB = cms.bool(False), )
directory=cms.string(options.input_directory), is_mc=cms.bool(options.is_mc), only_sample=cms.int32(options.only_sample), check_for_wrong_sample=cms.bool(options.check_for_wrong_sample), filename=cms.string(options.input_fn), min_muon_hits=cms.int32(options.min_muon_hits), min_pixel_layers=cms.int32(options.min_pixel_layers), min_strip_layers=cms.int32(options.min_strip_layers), max_tpfms_station=cms.int32(-1), no_dt_allowed=cms.bool(options.muon_subdet == 'csconly'), no_csc_allowed=cms.bool(options.muon_subdet == 'dtonly'), check_tksta_dphi=cms.bool(False), check_shared_hits=cms.bool(options.check_shared_hits), use_unpropagated_values=cms.bool(True), pp_reco_mode=cms.bool(False), force_run_list=cms.vuint32(run_list), force_event_list=cms.vuint32(), require_tt25=cms.bool(options.require_tt25), require_not_tt25=cms.bool(options.require_not_tt25), min_bfield=cms.double(3.7), copy_selected_events=cms.bool(options.copy_selected_events), make_scatterplots=cms.bool(options.make_scatterplots), ) if options.bin_by_run: if not run_list: raise ValueError('bin_by_run specified but run_list is empty') bins = make_bins('run', run_list) elif options.bin_by_phi: bins = make_bins('phi', [x / 100. for x in range(-320, 321, 32)]) else:
FudgeLayer = cms.vuint32( 104, # Pixel Barrel services 104, # Pixel Barrel services 104, # Pixel Barrel services 104, # Pixel Barrel services 104, # Pixel Barrel services 104, # Pixel Barrel services 104, # Pixel Barrel services 106, # Pixel wall 106, # Pixel wall 107, # Pixel endcap services 107, # Pixel endcap services 107, # Pixel endcap services 6, # TIB1 services 7, # TIB2 services 8, # TIB3 services 9, # TIB4 services 10, # TID Layer 1 11, # TID Layer 2 12, # TID Layer 3 110, # TID outside services 110, # TID outside services 111, # TOB inside services 111, # TOB inside services 13, # TOB Layer1 13, # TOB Layer1 13, # TOB Layer1 13, # TOB Layer1 14, # TOB Layer2 14, # TOB Layer2 14, # TOB Layer2 14, # TOB Layer2 15, # TOB Layer3 15, # TOB Layer3 15, # TOB Layer3 15, # TOB Layer3 16, # TOB Layer4 16, # TOB Layer4 16, # TOB Layer4 16, # TOB Layer4 17, # TOB Layer5 17, # TOB Layer5 17, # TOB Layer5 17, # TOB Layer5 18, # TOB Layer6 18, # TOB Layer6 18, # TOB Layer6 18, # TOB Layer6 112, # TOB services 112, # TOB services 112, # TOB services 19, # TEC Layer 1 19, # TEC Layer 1 19, # TEC Layer 1 20, # TEC Layer 2 20, # TEC Layer 2 20, # TEC Layer 2 21, # TEC Layer 3 21, # TEC Layer 3 21, # TEC Layer 3 22, # TEC Layer 4 22, # TEC Layer 4 22, # TEC Layer 4 23, # TEC Layer 5 23, # TEC Layer 5 23, # TEC Layer 5 24, # TEC Layer 6 24, # TEC Layer 6 24, # TEC Layer 6 25, # TEC Layer 7 25, # TEC Layer 7 26, # TEC Layer 8 26, # TEC Layer 8 27, # TEC Layer 9 27, # TEC Layer 9 113, # Barrel Wall 114, # Endcap Wall : 4.86<eta<4.91 114, # Endcap Wall : 4.82<eta<4.86 114, # Endcap Wall : 4.40<eta<4.82 114, # Endcap Wall : 4.00<eta<4.40 114, # Endcap Wall : 3.71<eta<4.00 114, # Endcap Wall : 3.61<eta<3.71 114, # Endcap Wall : 3.30<eta<3.61 114, # Endcap Wall : 3.05<eta<3.30 114, # Endcap Wall : 2.95<eta<3.05 114, # Endcap Wall : 2.75<eta<2.95 114 # Endcap Wall ),
# import of standard configurations process.load('Configuration/StandardSequences/Services_cff') process.load('FWCore/MessageService/MessageLogger_cfi') process.load('Configuration/EventContent/EventContent_cff') # module to select events based on HLT Trigger bits process.load('HLTrigger/HLTfilters/hltHighLevelDev_cfi') # Loading "hltHighLevelDev_cfi" defines an EDFilter called hltHighLevelDev # now we can configure it process.hltHighLevelDev.TriggerResultsTag = cms.InputTag("TriggerResults","","HLT8E29") process.hltHighLevelDev.HLTPaths = ("HLT_Mu9",) process.hltHighLevelDev.andOr = True # True = OR, False = AND process.hltHighLevelDev.HLTPathsPrescales = cms.vuint32(1,) process.hltHighLevelDev.HLTOverallPrescale = cms.uint32 (100) # Once the hltHighLevelDev filter has been defined I can clone it and configure as needed process.hltHighLevelDev2 = process.hltHighLevelDev.clone(andOr = True) process.hltHighLevelDev2.HLTPaths = ("HLT_L1Mu20","HLT_L2Mu9","HLT_L2Mu11","HLT_L1Mu14_L1SingleEG10","HLT_L1Mu14_L1SingleJet6U","HLT_L1Mu14_L1ETM30") process.hltHighLevelDev2.andOr = True # True = OR, False = AND process.hltHighLevelDev2.HLTPathsPrescales = cms.vuint32(1,1,1,1,1,1) process.hltHighLevelDev2.HLTOverallPrescale = cms.uint32 (100) # Once the hltHighLevelDev filter has been defined I can clone it and configure as needed process.hltHighLevelDev3 = process.hltHighLevelDev.clone(andOr = True) process.hltHighLevelDev3.HLTPaths = ("HLT_L1DoubleMuOpen","HLT_DoubleMu0","HLT_DoubleMu3") process.hltHighLevelDev3.andOr = True # True = OR, False = AND process.hltHighLevelDev3.HLTPathsPrescales = cms.vuint32(5,1,1) process.hltHighLevelDev3.HLTOverallPrescale = cms.uint32 (100)
TriggerResults=cms.InputTag('TriggerResults', '', 'HLT'), HLTPathNames=cms.vstring('HLT_PFMET120_PFMHT120_IDTight_v', 'HLT_PFHT300PT30_QuadPFJet_75_60_45_40_v', 'HLT_PFHT380_SixPFJet32_DoublePFBTagCSV_', 'HLT_PFHT380_SixPFJet32_DoublePFBTagDeepCSV_', 'HLT_IsoMu24_eta2p1_v'), JetTag=cms.VInputTag( cms.InputTag("hltCombinedSecondaryVertexBJetTagsCalo"), cms.InputTag("hltCombinedSecondaryVertexBJetTagsCalo"), cms.InputTag("hltCombinedSecondaryVertexBJetTagsCalo"), cms.InputTag("hltCombinedSecondaryVertexBJetTagsPF"), cms.InputTag("hltCombinedSecondaryVertexBJetTagsPF"), ), MinJetPT=cms.double(20), mcFlavours=cms.PSet( light=cms.vuint32(1, 2, 3, 21), # udsg c=cms.vuint32(4), b=cms.vuint32(5), g=cms.vuint32(21), uds=cms.vuint32(1, 2, 3)), mcPartons=cms.InputTag("hltBtagJetsbyValAlgo")) #put all in a path hltbtagValidationSequence = cms.Sequence( # remove noisy warnings # hltBtagTriggerSelection + hltBtagJetMCTools + HltVertexValidationVertices + hltbTagValidation) # fastsim customs from Configuration.Eras.Modifier_fastSim_cff import fastSim fastSim.toModify(HltVertexValidationVertices,
"HLT_IsoMu15_L1ETM20_v1", "HLT_IsoMu15_L1ETM20_v2", "HLT_IsoMu15_L1ETM20_v3", "HLT_IsoMu15_L1ETM20_v4", "HLT_IsoPFTau35_Trk20_MET45_v1", "HLT_IsoPFTau35_Trk20_MET45_v2", "HLT_IsoPFTau35_Trk20_MET45_v4", "HLT_IsoPFTau35_Trk20_MET45_v6", "HLT_IsoPFTau35_Trk20_v2", "HLT_IsoPFTau35_Trk20_v3", "HLT_IsoPFTau35_Trk20_v4", "HLT_IsoPFTau35_Trk20_v6", "HLT_IsoPFTau35_Trk20_MET60_v2", "HLT_IsoPFTau35_Trk20_MET60_v3", "HLT_IsoPFTau35_Trk20_MET60_v4", "HLT_IsoPFTau35_Trk20_MET60_v6", "HLT_IsoPFTau45_Trk20_MET60_v2", "HLT_IsoPFTau45_Trk20_MET60_v3", "HLT_IsoPFTau45_Trk20_MET60_v4", "HLT_IsoPFTau35_Trk20_MET70_v2", "HLT_MediumIsoPFTau35_Trk20_v1", "HLT_MediumIsoPFTau35_Trk20_v5", "HLT_MediumIsoPFTau35_Trk20_v6", "HLT_MediumIsoPFTau35_Trk20_MET60_v1", "HLT_MediumIsoPFTau35_Trk20_MET60_v5", "HLT_MediumIsoPFTau35_Trk20_MET60_v6", "HLT_MediumIsoPFTau35_Trk20_MET70_v1", "HLT_MediumIsoPFTau35_Trk20_MET70_v5", "HLT_MediumIsoPFTau35_Trk20_MET70_v6", ), L1TauTriggerSource = cms.InputTag("tteffL1GTSeed"), L1JetMatchingCone = cms.double(0.5), L1JetMatchingMode = cms.string("nearestDR"), # "nearestDR", "highestEt" L1IsolationThresholds = cms.vuint32(1,2,3,4), # count regions with "et() < threshold", these are in GeV L2AssociationCollection = cms.InputTag("openhltL2TauIsolationProducer"), EERecHits = cms.untracked.InputTag("ecalRecHit","EcalRecHitsEE"), EBRecHits = cms.untracked.InputTag("ecalRecHit","EcalRecHitsEB"), CaloTowers = cms.untracked.InputTag("towerMaker"), outerCone = cms.untracked.double(0.5), innerCone = cms.untracked.double(0.15), crystalThresholdEB = cms.untracked.double(0.15), crystalThresholdEE = cms.untracked.double(0.45), L2matchingDeltaR = cms.double(0.2), l25JetSource = cms.InputTag("hltPFTauTagInfo"), l25PtCutSource = cms.InputTag("hltPFTaus"), l3IsoSource = cms.InputTag("hltL3TauIsolationSelector"), #obsolet: L25/L3 merged? l25MatchingCone = cms.double(0.3), l3IsoQualityCuts = process.hltPFTauLooseIsolationDiscriminator.qualityCuts.isolationQualityCuts.clone(), MCMatchingCone = cms.double(0.2),
process.EvFDaqDirector = cms.Service("EvFDaqDirector", buBaseDir = cms.untracked.string(options.buBaseDir), baseDir = cms.untracked.string(options.dataDir), directorIsBU = cms.untracked.bool(False ), testModeNoBuilderUnit = cms.untracked.bool(False), runNumber = cms.untracked.uint32(options.runNumber) ) process.PrescaleService = cms.Service( "PrescaleService", lvl1DefaultLabel = cms.string( "B" ), lvl1Labels = cms.vstring( 'A', 'B' ), prescaleTable = cms.VPSet( cms.PSet( pathName = cms.string( "p1" ), prescales = cms.vuint32( 0, 10) ), cms.PSet( pathName = cms.string( "p2" ), prescales = cms.vuint32( 0, 100) ) )) process.source = cms.Source("FedRawDataInputSource", getLSFromFilename = cms.untracked.bool(True), testModeNoBuilderUnit = cms.untracked.bool(False), eventChunkSize = cms.untracked.uint32(128), numBuffers = cms.untracked.uint32(2), eventChunkBlock = cms.untracked.uint32(128), useL1EventID=cms.untracked.bool(True) )
) process.source = cms.Source("PoolSource", fileNames = cms.untracked.vstring( 'file:reco_9_1_ZAm.root' ) ) # 'dcap://pnfs/cms/WAX/resilient/spadhi/CMS/spadhi/PhysicsProcesses_TopologyT2_38xFall10/PhysicsProcesses_TopologyT2_38xFall10/f959c379445f9d7540d41cf9ffa87a96/ROOTFILE' process.source.fileNames = [ 'dcap:/ROOTFILE' ] process.load("RecoJets.Configuration.GenJetParticles_cff") process.genParticlesForJets.ignoreParticleIDs.extend(cms.vuint32(5000039)) process.load("RecoJets.JetProducers.ak5GenJets_cfi") process.printGenParticles = cms.EDAnalyzer("ParticleListDrawer", src = cms.InputTag("genParticles"), # src = cms.InputTag("ak5GenJets"), maxEventsToPrint = cms.untracked.int32(1) ) process.MessageLogger = cms.Service("MessageLogger", cout2 = cms.untracked.PSet( default = cms.untracked.PSet( limit = cms.untracked.int32(0)
import FWCore.ParameterSet.Config as cms # parameters for HIPAlignmentAlgorithm HIPAlignmentAlgorithm = cms.PSet( algoName=cms.string('HIPAlignmentAlgorithm'), debug=cms.bool(False), verbosity=cms.bool(False), checkDbAlignmentValidity=cms.bool(False), isCollision=cms.bool(True), UsePreSelection=cms.bool(False), multiIOV=cms.bool(False), IOVrange=cms.vuint32(1, 99999999), minRelParameterError=cms.double(0), maxRelParameterError=cms.double(-1), # -1 for no cut minimumNumberOfHits=cms.int32(1), maxAllowedHitPull=cms.double(-1), # -1 for no cut applyCutsPerComponent=cms.bool( False), # Overrides settings above for the specified detectors cutsPerComponent=cms.VPSet( cms.PSet( Selector=cms.PSet(alignParams=cms.vstring( "AllAlignables,000000" # Obligatory second string ) # can use "selected" for the already-specified alignables ), # Parameter cuts minRelParError=cms.double(0), maxRelParError=cms.double(-1), # -1 for no cut # Hit cuts minNHits=cms.int32(0), maxHitPull=cms.double(-1), # -1 for no cut applyPixelProbCut=cms.bool(False),
# module to select events based on HLT Trigger bits process.load('HLTrigger/HLTfilters/hltHighLevelDev_cfi') # Loading "hltHighLevelDev_cfi" defines an EDFilter called hltHighLevelDev # now we can configure it # All events from MinBias PD process.hltHighLevelDev.TriggerResultsTag = cms.InputTag( "TriggerResults", "", "HLT") process.hltHighLevelDev.HLTPaths = ('HLT_PhysicsDeclared', ) process.hltHighLevelDev.andOr = True # True = OR, False = AND process.hltHighLevelDev.throw = False # throw exception on unknown path names process.hltHighLevelDev.HLTPathsPrescales = cms.vuint32( 1, #'HLT_PhysicsDeclared', ) process.hltHighLevelDev.HLTOverallPrescale = cms.uint32(1) # All events from MinBias PD, prescaled by a factor 10 process.hltHighLevelDev2 = process.hltHighLevelDev.clone(andOr=True) process.hltHighLevelDev2.HLTPaths = ('HLT_PhysicsDeclared', ) process.hltHighLevelDev2.andOr = True # True = OR, False = AND process.hltHighLevelDev2.throw = False # throw exception on unknown path names process.hltHighLevelDev2.HLTPathsPrescales = cms.vuint32( 1, #'HLT_PhysicsDeclared', ) process.hltHighLevelDev2.HLTOverallPrescale = cms.uint32(10) # All interesting events.
import FWCore.ParameterSet.Config as cms l1RPCHwConfig = cms.ESProducer("RPCTriggerHwConfig", disableTowers = cms.vint32(), disableCrates = cms.vint32(), disableTowersInCrates = cms.vint32(), enableTowers = cms.vint32(), enableCrates = cms.vint32(), enableTowersInCrates = cms.vint32(), disableAll = cms.bool(False) ) rpchwconfsrc = cms.ESSource("EmptyESSource", recordName = cms.string('L1RPCHwConfigRcd'), iovIsRunNotTime = cms.bool(True), firstValid = cms.vuint32(1) )
# ThePhotons = cms.InputTag("slimmedPhotons"), # TheJets = cms.InputTag("slimmedJets"), # L1Maps = cms.string(relBase+"/src/L1Prefiring/EventWeightProducer/files/L1PrefiringMaps_new.root"), # L1Maps = cms.string("L1PrefiringMaps_new.root"), # update this line with the location of this file #L1Maps = cms.string("CMSSW_8_0_32/src/L1Prefiring/EventWeightProducer/data/L1PrefiringMaps_new.root"), # DataEra = cms.string("2016BtoH"), #Use 2016BtoH for 2016 # UseJetEMPt = cms.bool(False), #can be set to true to use jet prefiring maps parametrized vs pt(em) instead of pt # PrefiringRateSystematicUncty = cms.double(0.2) #Minimum relative prefiring uncty per object # ) process.load('RecoMET.METFilters.ecalBadCalibFilter_cfi') baddetEcallist = cms.vuint32([ 872439604, 872422825, 872420274, 872423218, 872423215, 872416066, 872435036, 872439336, 872420273, 872436907, 872420147, 872439731, 872436657, 872420397, 872439732, 872439339, 872439603, 872422436, 872439861, 872437051, 872437052, 872420649, 872422436, 872421950, 872437185, 872422564, 872421566, 872421695, 872421955, 872421567, 872437184, 872421951, 872421694, 872437056, 872437057, 872437313 ]) process.ecalBadCalibReducedMINIAODFilter = cms.EDFilter( "EcalBadCalibFilter", EcalRecHitSource=cms.InputTag("reducedEgamma:reducedEERecHits"), ecalMinEt=cms.double(50.), baddetEcal=baddetEcallist, taggingMode=cms.bool(True), debug=cms.bool(False)) process.treeDumper = cms.EDAnalyzer( "PKUTreeMaker", originalNEvents=cms.int32(1),
CA8QGTaggerPuppi = QGTagger.clone() CA8QGTaggerPuppi.srcJets = cms.InputTag('CA8PFJetsPuppi') CA8QGTaggerPuppi.jetsLabel = cms.string('QGL_AK4PFchs') CA8QGTaggerPuppi.jec = cms.InputTag("ak8PuppiL1FastL2L3Corrector") # NOTE: use "ca8PFPuppiL1FastL2L3Corrector" for MC / "ca8PFPuppiL1FastL2L3ResidualCorrector" for Data CA8QGTaggerSubJetsPuppi = CA8QGTaggerPuppi.clone() CA8QGTaggerSubJetsPuppi.srcJets = cms.InputTag('CA8caPFJetsSoftDropPuppi','SubJets') CA8QGTaggerSubJetsPuppi.jetsLabel = cms.string('QGL_AK4PFchs') CA8QGTaggerSubJetsPuppi.jec = cms.InputTag("ak8PuppiL1FastL2L3Corrector") # NOTE: use "ca8PFPuppiL1FastL2L3Corrector" for MC / "ca8PFPuppiL1FastL2L3ResidualCorrector" for Data # N-subjettiness from RecoJets.JetProducers.nJettinessAdder_cfi import Njettiness CA8NjettinessPuppi = Njettiness.clone( src = cms.InputTag('CA8PFJetsPuppi'), R0 = cms.double(0.8), Njets = cms.vuint32(1,2,3,4) ) CA8genjetsequencePuppi = cms.Sequence( CA8GenJetsPuppi* CA8FlavorPuppi ) # # Define sequences # CA8jetsequencePuppi = cms.Sequence( CA8PFJetsPuppi* CA8caPFJetsPrunedPuppi* CA8caPFJetsTrimmedPuppi* CA8caPFJetsSoftDropPuppi*
L1MuonCollectionTag=cms.InputTag("hltL1extraParticles"), L1UseL1TriggerObjectMaps=cms.bool(True), L1UseAliasesForSeeding=cms.bool(True), L1GtReadoutRecordTag=cms.InputTag("hltGtDigis"), L1CollectionsTag=cms.InputTag("hltL1extraParticles"), L1NrBxInEvent=cms.int32(3), L1GtObjectMapTag=cms.InputTag("hltL1GtObjectMap"), L1TechTriggerSeeding=cms.bool(False)) fragment.hltPreZeroBias = cms.EDFilter( "HLTPrescaler", L1GtReadoutRecordTag=cms.InputTag("hltGtDigis"), offset=cms.uint32(0)) fragment.hltFEDSelector = cms.EDProducer( "EvFFEDSelector", inputTag=cms.InputTag("rawDataCollector"), fedList=cms.vuint32(1023)) fragment.hltTriggerSummaryAOD = cms.EDProducer("TriggerSummaryProducerAOD", processName=cms.string("@")) fragment.hltTriggerSummaryRAW = cms.EDProducer("TriggerSummaryProducerRAW", processName=cms.string("@")) fragment.hltPreHLTAnalyzerEndpath = cms.EDFilter( "HLTPrescaler", L1GtReadoutRecordTag=cms.InputTag("hltGtDigis"), offset=cms.uint32(0)) fragment.hltL1GtTrigReport = cms.EDAnalyzer( "L1GtTrigReport", PrintVerbosity=cms.untracked.int32(10), UseL1GlobalTriggerRecord=cms.bool(False), PrintOutput=cms.untracked.int32(3), L1GtRecordInputTag=cms.InputTag("hltGtDigis")) fragment.hltTrigReport = cms.EDAnalyzer(
connect = cms.string(dbconnection) ) process.CondDBCommon.connect = cms.string( dbconnection ) process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(maxevts) ) process.source = cms.Source("PoolSource", fileNames = cms.untracked.vstring( inputfile ) ) process.rpcTechnicalTrigger = cms.EDProducer('RPCTechnicalTrigger', TriggerMode = cms.int32(1), RPCDigiLabel = cms.InputTag("muonRPCDigis"), UseDatabase = cms.untracked.int32(0), ConfigFile = cms.untracked.string("hardware-pseudoconfig.txt"), BitNumbers=cms.vuint32(24,25,26,27,28), BitNames=cms.vstring('L1Tech_rpcBit1', 'L1Tech_rpcBit2', 'L1Tech_rpcBit3', 'L1Tech_rpcBit4', 'L1Tech_rpcBit5') ) process.out = cms.OutputModule("PoolOutputModule", fileName = cms.untracked.string('rpcttbits.root'), outputCommands = cms.untracked.vstring('drop *','keep L1GtTechnicalTriggerRecord_*_*_*') ) process.p = cms.Path(process.rpcTechnicalTrigger) process.e = cms.EndPath(process.out)
import FWCore.ParameterSet.Config as cms process = cms.Process("HLTFake2") process.HLTConfigVersion = cms.PSet( tableName=cms.string('/dev/CMSSW_11_3_0/Fake2/V8')) process.streams = cms.PSet(A=cms.vstring('InitialPD')) process.datasets = cms.PSet(InitialPD=cms.vstring( 'HLT_Physics_v1', 'HLT_Random_v1', 'HLT_ZeroBias_v1')) process.GlobalParametersRcdSource = cms.ESSource( "EmptyESSource", iovIsRunNotTime=cms.bool(True), recordName=cms.string("L1TGlobalParametersRcd"), firstValid=cms.vuint32(1)) process.GlobalTag = cms.ESSource( "PoolDBESSource", globaltag=cms.string("80X_dataRun2_HLT_v12"), RefreshEachRun=cms.untracked.bool(False), snapshotTime=cms.string(""), toGet=cms.VPSet(), DBParameters=cms.PSet( authenticationPath=cms.untracked.string("."), connectionRetrialTimeOut=cms.untracked.int32(60), idleConnectionCleanupPeriod=cms.untracked.int32(10), messageLevel=cms.untracked.int32(0), enablePoolAutomaticCleanUp=cms.untracked.bool(False), enableConnectionSharing=cms.untracked.bool(True), enableReadOnlySessionOnUpdateConnection=cms.untracked.bool(False), connectionTimeOut=cms.untracked.int32(0),
process.patElectronsWithTrigger = cms.EDProducer("PATTriggerMatchElectronEmbedder", src = cms.InputTag("patElectrons"), matches = cms.VInputTag(cms.InputTag('eleTriggerMatchHLT')) ) switchOnTriggerMatching( process, ['eleTriggerMatchHLT' ],sequence ='patDefaultSequence', hltProcess = '*' ) ##################### # # # GEN JET PROP # # # ##################### genParticlesForJetsNoNu = process.genParticlesForJets.clone() process.genParticlesForJetsNoNu.ignoreParticleIDs += cms.vuint32( 12,14,16) GenJetParameters = cms.PSet( src = cms.InputTag("genParticlesForJets"), srcPVs = cms.InputTag(''), jetType = cms.string('GenJet'), jetPtMin = cms.double(3.0), inputEtMin = cms.double(0.0), inputEMin = cms.double(0.0), doPVCorrection = cms.bool(False), # pileup with offset correction doPUOffsetCorr = cms.bool(False), # if pileup is false, these are not read: nSigmaPU = cms.double(1.0), radiusPU = cms.double(0.5),
process.RPStraightTrackAligner.resolveRotZ = True process.RPStraightTrackAligner.resolveShZ = False process.RPStraightTrackAligner.resolveRPShZ = False process.RPStraightTrackAligner.algorithms = cms.vstring('Ideal', 'Jan') process.RPStraightTrackAligner.constraintsType = "final" process.RPStraightTrackAligner.useExtendedRotZConstraint = True process.RPStraightTrackAligner.useZeroThetaRotZConstraint = True process.RPStraightTrackAligner.useExtendedShZConstraints = True process.RPStraightTrackAligner.useExtendedRPShZConstraint = True process.RPStraightTrackAligner.homogeneousConstraints.RotZ_values = cms.vdouble( 0, 0, 0, 0) process.RPStraightTrackAligner.fixedDetectorsConstraints.ShR.ids = cms.vuint32( 1200, 1201, 1248, 1249) process.RPStraightTrackAligner.fixedDetectorsConstraints.RotZ.ids = cms.vuint32( 1200, 1201) process.RPStraightTrackAligner.fixedDetectorsConstraints.ShZ.ids = cms.vuint32( 1200, 1201, 1248, 1249) process.RPStraightTrackAligner.fixedDetectorsConstraints.RPShZ.ids = cms.vuint32( 1200, 1240) process.RPStraightTrackAligner.fixedDetectorsConstraints.RPShZ.values = cms.vdouble( 0, 0) process.RPStraightTrackAligner.JanAlignmentAlgorithm.stopOnSingularModes = False result_dir = '.' process.RPStraightTrackAligner.saveIntermediateResults = False process.RPStraightTrackAligner.taskDataFileName = result_dir + "/task_data.root"
from RecoJets.Configuration.RecoPFJets_cff import * from RecoJets.JetProducers.AnomalousCellParameters_cfi import * from RecoJets.JetProducers.PFJetParameters_cfi import * process.chs = cms.EDFilter("CandPtrSelector", src = cms.InputTag('packedPFCandidates'), cut = cms.string('fromPV') ) process.ak4PFJetsCHS = ak4PFJetsCHS.clone( src = 'chs' ) process.ak4PFJetsCHS.doAreaFastjet = True process.ak8CHSJets = ak8PFJetsCHS.clone( src = 'chs', jetPtMin = fatjet_ptmin ) process.NjettinessAK8 = cms.EDProducer("NjettinessAdder", src = cms.InputTag("ak8CHSJets"), Njets = cms.vuint32(1, 2, 3, 4), # variables for measure definition : measureDefinition = cms.uint32( 0 ), # CMS default is normalized measure beta = cms.double(1.0), # CMS default is 1 R0 = cms.double( 0.8 ), # CMS default is jet cone size Rcutoff = cms.double( 999.0), # not used by default # variables for axes definition : axesDefinition = cms.uint32( 6 ), # CMS default is 1-pass KT axes nPass = cms.int32(999), # not used by default akAxesR0 = cms.double(-999.0) # not used by default ) if config["DOAK10TRIMMEDRECLUSTERING"]: process.ECFAK10 = cms.EDProducer("ECFAdder", src = cms.InputTag("ak10CHSJetsTrimmed"), Njets = cms.vuint32(1, 2, 3),
maxMass=cms.string('1.6'), minMass=cms.double(0.9), nCharged=cms.uint32(3), nChargedPFCandsMin=cms.uint32(1), nPiZeros=cms.uint32(1), nTracksMin=cms.uint32(2))), matchingCone=cms.double(0.5), minPixelHits=cms.int32(1), minTauPt=cms.double(0.0), requireTauChargedHadronsToBeChargedPFCands=cms.bool(False)) process.hpsPFTauProducerSansRefs = cms.EDProducer( "RecoBaseTauCleaner", cleaners=cms.VPSet( cms.PSet(name=cms.string('Charge'), nprongs=cms.vuint32(1, 3), passForCharge=cms.int32(1), plugin=cms.string('RecoBaseTauChargeCleanerPlugin'), selectionFailValue=cms.double(0)), cms.PSet(name=cms.string('HPS_Select'), plugin=cms.string('RecoBaseTauDiscriminantCleanerPlugin'), src=cms.InputTag("hpsSelectionDiscriminator")), cms.PSet( minTrackPt=cms.double(5.0), name=cms.string('killSoftTwoProngTaus'), plugin=cms.string('RecoBaseTauSoftTwoProngTausCleanerPlugin')), cms.PSet(name=cms.string('ChargedHadronMultiplicity'), plugin=cms.string( 'RecoBaseTauChargedHadronMultiplicityCleanerPlugin')), cms.PSet(name=cms.string('Pt'), plugin=cms.string('RecoBaseTauStringCleanerPlugin'),
## ---- Services ---- process.load("DQM.SiStripCommon.MessageLogger_cfi") process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1)) ## Global tag - see http://twiki.cern.ch/twiki/bin/view/CMS/SWGuideFrontierConditions process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff") process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:run2_data', '') process.load("Configuration.Geometry.GeometryRecoDB_cff") # --- The unpacking configuration --- process.load('DQM.SiStripMonitorHardware.SiStripSpyUnpacker_cfi') process.load('DQM.SiStripMonitorHardware.SiStripSpyDigiConverter_cfi') ## * Scope digi settings process.SiStripSpyUnpacker.FEDIDs = cms.vuint32() #use a subset of FEDs or leave empty for all. #process.SiStripSpy.FEDIDs = cms.vuint32(50, 187, 260, 356) #one from each partition process.SiStripSpyUnpacker.InputProductLabel = cms.InputTag('source') process.SiStripSpyUnpacker.AllowIncompleteEvents = True process.SiStripSpyUnpacker.StoreCounters = True process.SiStripSpyUnpacker.StoreScopeRawDigis = cms.bool(True) # Note - needs to be True for use in other modules. ## * Module digi settings process.SiStripSpyDigiConverter.InputProductLabel = cms.InputTag('SiStripSpyUnpacker','ScopeRawDigis') process.SiStripSpyDigiConverter.StorePayloadDigis = True process.SiStripSpyDigiConverter.StoreReorderedDigis = True process.SiStripSpyDigiConverter.StoreModuleDigis = True process.SiStripSpyDigiConverter.StoreAPVAddress = True process.SiStripSpyDigiConverter.MinDigiRange = 100 process.SiStripSpyDigiConverter.MaxDigiRange = 1024 process.SiStripSpyDigiConverter.MinZeroLight = 0 process.SiStripSpyDigiConverter.MaxZeroLight = 1024
"isPFMuon && innerTrack.validFraction >= 0.49 && ( isGlobalMuon && globalTrack.normalizedChi2 < 3 && combinedQuality.chi2LocalPosition < 12 && combinedQuality.trkKink < 20 && segmentCompatibility >= 0.303 || segmentCompatibility >= 0.451 )" ), electronSel=cms.string(""), tauSel=cms.string(""), photonSel=cms.string(""), jetName=cms.string("Jet"), muonName=cms.string("Muon"), electronName=cms.string("Electron"), tauName=cms.string("Tau"), photonName=cms.string("Photon")) genWeightsTable = cms.EDProducer( "GenWeightsTableProducer", genEvent=cms.InputTag("generator"), lheInfo=cms.InputTag("externalLHEProducer"), preferredPDFs=cms.vuint32(91400, 260001), namedWeightIDs=cms.vstring(), namedWeightLabels=cms.vstring(), lheWeightPrecision=cms.int32(14), maxPdfWeights=cms.uint32( 50), # for NNPDF, keep only the first 50 replicas (save space) debug=cms.untracked.bool(False), ) lheInfoTable = cms.EDProducer( "LHETablesProducer", lheInfo=cms.InputTag("externalLHEProducer"), ) l1bits = cms.EDProducer("L1TriggerResultsConverter", src=cms.InputTag("gtStage2Digis"), legacyL1=cms.bool(False))
ObjectsFromDigis = [ "2d" , "pixelray", "globalgeometry" , "bpphel" ] # Set Stub production pt threshold thresh = 2.0 # Set Tracklet flags trackletonlyherm = cms.bool(True) # True creates only hermetic-FNAL style Tracklets tracklethelixfit = cms.bool(True) # True fits the Tracklet with helicoidal style in R-Z, creating dedicated collection trackletbeamspot = cms.bool(True) # True corrects the Tracklet fit for the beam position, creating dedicated collection # Set Track flags trackletseedtype = "ShortTrackletsVtxBSHelFit" # Recommended because the only algorithm available is bpphel! # NOTE be self consistent with available seeds # and with available fit strategies trackhelixfit = cms.bool(True) # True applies helicoidal fit to Tracks trackalsoseed = cms.bool(True) # True includes seed vertex within points used for the fit tracknoduplic = cms.bool(False) # True removes duplicate candidate Tracks # NOTE only one collection is created by the L1TkTrackBuilder trackdsvector = cms.vuint32( ) # empty to allow all of them index = 1 puindex = 111 if os.getenv('pgun_index'): index = str(os.getenv('pgun_index')) puindex = int(os.getenv('pgun_index')) theSeedValue = 12345 if os.getenv('pgun_index'): theSeedValue = int(12345+1000*int(os.getenv('pgun_index'))*int(os.getenv('pgun_index'))) print 'Seed: '+str(theSeedValue)
process.puppi.useExistingWeights = True process.puppi.candName = cms.InputTag('packedPFCandidates') process.puppi.vertexName = cms.InputTag('offlineSlimmedPrimaryVertices') process.ak8PFJetsCHS = ak8PFJetsCHS.clone(src='puppi', jetPtMin=100.0) process.ak8PFJetsCHSPruned = ak8PFJetsCHSPruned.clone(src='puppi', jetPtMin=100.0) process.ak8PFJetsCHSPrunedMass = ak8PFJetsCHSPrunedMass.clone() process.ak8PFJetsCHSSoftDrop = ak8PFJetsCHSSoftDrop.clone(src='puppi', jetPtMin=100.0) process.ak8PFJetsCHSSoftDropMass = ak8PFJetsCHSSoftDropMass.clone() process.NjettinessAK8 = cms.EDProducer( "NjettinessAdder", src=cms.InputTag("ak8PFJetsCHS"), Njets=cms.vuint32(1, 2, 3, 4), # variables for measure definition : measureDefinition=cms.uint32(0), # CMS default is normalized measure beta=cms.double(1.0), # CMS default is 1 R0=cms.double(0.8), # CMS default is jet cone size Rcutoff=cms.double(999.0), # not used by default # variables for axes definition : axesDefinition=cms.uint32(6), # CMS default is 1-pass KT axes nPass=cms.int32(0), # not used by default akAxesR0=cms.double(-999.0) # not used by default ) process.substructureSequence = cms.Sequence() process.substructureSequence += process.puppi process.substructureSequence += process.ak8PFJetsCHS process.substructureSequence += process.NjettinessAK8
process.options = cms.untracked.PSet(wantSummary = cms.untracked.bool(True)) process.options.allowUnscheduled = cms.untracked.bool(False) process.puppiNoLeptons = process.puppi.clone() process.puppiNoLeptons.candName = cms.InputTag('pfNoLeptons') process.ElectronIsolation = cms.EDProducer("CITKPFIsolationSumProducer", srcToIsolate = cms.InputTag("slimmedElectrons"), srcForIsolationCone = cms.InputTag('packedPFCandidates'), isolationConeDefinitions = cms.VPSet( cms.PSet( isolationAlgo = cms.string('ElectronPFIsolationWithConeVeto'), coneSize = cms.double(0.3), VetoConeSizeEndcaps = cms.double(0.015), VetoConeSizeBarrel = cms.double(0.0), isolateAgainst = cms.string('h+'), miniAODVertexCodes = cms.vuint32(2,3) ), cms.PSet( isolationAlgo = cms.string('ElectronPFIsolationWithConeVeto'), coneSize = cms.double(0.3), VetoConeSizeEndcaps = cms.double(0.0), VetoConeSizeBarrel = cms.double(0.0), isolateAgainst = cms.string('h0'), miniAODVertexCodes = cms.vuint32(2,3) ), cms.PSet( isolationAlgo = cms.string('ElectronPFIsolationWithConeVeto'), coneSize = cms.double(0.3), VetoConeSizeEndcaps = cms.double(0.08), VetoConeSizeBarrel = cms.double(0.0), isolateAgainst = cms.string('gamma'), miniAODVertexCodes = cms.vuint32(2,3) ) ) )
import FWCore.ParameterSet.Config as cms disconnectedTriggerLayers = [ 2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28 ] geometry = cms.PSet( TriggerGeometryName=cms.string('HGCalTriggerGeometryHexLayerBasedImp1'), L1TCellsMapping=cms.FileInPath( "L1Trigger/L1THGCal/data/triggercell_mapping_8inch_aligned_192_432_V8_0.txt" ), L1TModulesMapping=cms.FileInPath( "L1Trigger/L1THGCal/data/panel_mapping_tdr_0.txt"), L1TCellNeighborsMapping=cms.FileInPath( "L1Trigger/L1THGCal/data/triggercell_neighbor_mapping_8inch_aligned_192_432_0.txt" ), L1TCellsBHMapping=cms.FileInPath( "L1Trigger/L1THGCal/data/triggercell_mapping_BH_3x3_30deg_0.txt"), L1TCellNeighborsBHMapping=cms.FileInPath( "L1Trigger/L1THGCal/data/triggercell_neighbor_mapping_BH_3x3_30deg_0.txt" ), DisconnectedModules=cms.vuint32(0), DisconnectedLayers=cms.vuint32(disconnectedTriggerLayers)) hgcalTriggerGeometryESProducer = cms.ESProducer( 'HGCalTriggerGeometryESProducer', TriggerGeometry=geometry, isV9Geometry=cms.bool(True))
import FWCore.ParameterSet.Config as cms pfMassDecorrelatedDeepBoostedJetPreprocessParams = cms.PSet( input_names = cms.vstring( 'pfcand', 'sv' ), pfcand = cms.PSet( input_shape = cms.vuint32(1, 36, 100, 1), var_infos = cms.PSet( pfcand_VTX_ass = cms.PSet( median = cms.double(7.0), norm_factor = cms.double(1.0) ), pfcand_abseta = cms.PSet( median = cms.double(0.632648706436), norm_factor = cms.double(1.59032225958) ), pfcand_btagEtaRel = cms.PSet( median = cms.double(1.19703966379), norm_factor = cms.double(0.521026991705) ), pfcand_btagJetDistVal = cms.PSet( median = cms.double(-0.000215483247302), norm_factor = cms.double(161.385119349) ), pfcand_btagPParRatio = cms.PSet( median = cms.double(0.839023888111), norm_factor = cms.double(1.19186117841) ), pfcand_btagPtRatio = cms.PSet(
"RecoPixelVertexing/PixelLowPtUtilities/data/pixelShapePhase1_noL1.par" )) #SiStrip Local Reco process.load("CalibTracker.SiStripCommon.TkDetMapESProducer_cfi") #---- for P5 (online) DB access process.load("DQM.Integration.config.FrontierCondition_GT_cfi") # Condition for lxplus: change and possibly customise the GT #from Configuration.AlCa.GlobalTag import GlobalTag as gtCustomise #process.GlobalTag = gtCustomise(process.GlobalTag, 'auto:run2_data', '') process.hltESSHcalSeverityLevel = cms.ESSource( "EmptyESSource", iovIsRunNotTime=cms.bool(True), recordName=cms.string("HcalSeverityLevelComputerRcd"), firstValid=cms.vuint32(1)) process.hcalRecAlgos = cms.ESProducer( "HcalRecAlgoESProducer", RecoveredRecHitBits=cms.vstring('TimingAddedBit', 'TimingSubtractedBit'), SeverityLevels=cms.VPSet( cms.PSet(ChannelStatus=cms.vstring(), RecHitFlags=cms.vstring(), Level=cms.int32(0)), cms.PSet(ChannelStatus=cms.vstring('HcalCellCaloTowerProb'), RecHitFlags=cms.vstring(), Level=cms.int32(1)), cms.PSet( ChannelStatus=cms.vstring('HcalCellExcludeFromHBHENoiseSummary'), RecHitFlags=cms.vstring('HSCP_R1R2', 'HSCP_FracLeader', 'HSCP_OuterEnergy', 'HSCP_ExpFit', 'ADCSaturationBit', 'HBHEIsolatedNoise',
L1MuonCollectionTag = cms.InputTag( "hltL1extraParticles" ), L1UseL1TriggerObjectMaps = cms.bool( True ), L1UseAliasesForSeeding = cms.bool( True ), L1GtReadoutRecordTag = cms.InputTag( "hltGtDigis" ), L1CollectionsTag = cms.InputTag( "hltL1extraParticles" ), L1NrBxInEvent = cms.int32( 3 ), L1GtObjectMapTag = cms.InputTag( "hltL1GtObjectMap" ), L1TechTriggerSeeding = cms.bool( False ) ) fragment.hltPreZeroBias = cms.EDFilter( "HLTPrescaler", L1GtReadoutRecordTag = cms.InputTag( "hltGtDigis" ), offset = cms.uint32( 0 ) ) fragment.hltFEDSelector = cms.EDProducer( "EvFFEDSelector", inputTag = cms.InputTag( "rawDataCollector" ), fedList = cms.vuint32( 1023 ) ) fragment.hltTriggerSummaryAOD = cms.EDProducer( "TriggerSummaryProducerAOD", moduleLabelPatternsToSkip = cms.vstring( ), processName = cms.string( "@" ), moduleLabelPatternsToMatch = cms.vstring( 'hlt*' ), throw = cms.bool( False ) ) fragment.hltTriggerSummaryRAW = cms.EDProducer( "TriggerSummaryProducerRAW", processName = cms.string( "@" ) ) fragment.hltPreHLTAnalyzerEndpath = cms.EDFilter( "HLTPrescaler", L1GtReadoutRecordTag = cms.InputTag( "hltGtDigis" ), offset = cms.uint32( 0 ) ) fragment.hltL1GtTrigReport = cms.EDAnalyzer( "L1GtTrigReport",
fileName=cms.untracked.string('/tmp/azzi/WZMuStream.root'), dataset=cms.untracked.PSet(dataTier=cms.untracked.string('RAW-RECO'), filterName=cms.untracked.string('WZMuFilter')), SelectEvents=cms.untracked.PSet(SelectEvents=cms.vstring( 'ZMuSkimPath', 'WtcMetSkimPath', 'WpfMetSkimPath'))) #################################logerrorharvester############################################ process.load("FWCore.Modules.logErrorFilter_cfi") from Configuration.StandardSequences.RawToDigi_Data_cff import gtEvmDigis process.gtEvmDigis = gtEvmDigis.clone() process.stableBeam = cms.EDFilter( "HLTBeamModeFilter", L1GtEvmReadoutRecordTag=cms.InputTag("gtEvmDigis"), AllowedBeamMode=cms.vuint32(11)) process.logerrorpath = cms.Path(process.gtEvmDigis + process.stableBeam + process.logErrorFilter) process.outlogerr = cms.OutputModule( "PoolOutputModule", outputCommands=process.FEVTEventContent.outputCommands, fileName=cms.untracked.string('/tmp/azzi/logerror_filter.root'), dataset=cms.untracked.PSet( dataTier=cms.untracked.string('RAW-RECO'), filterName=cms.untracked.string('Skim_logerror')), SelectEvents=cms.untracked.PSet(SelectEvents=cms.vstring("logerrorpath"))) #====================== ######################################TPG Performance SKIMS#####################################