Ejemplo n.º 1
0
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
Ejemplo n.º 2
0
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
Ejemplo n.º 3
0
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
Ejemplo n.º 4
0
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
Ejemplo n.º 5
0
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
Ejemplo n.º 6
0
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
Ejemplo n.º 7
0
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
Ejemplo n.º 8
0
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) 
Ejemplo n.º 9
0
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
Ejemplo n.º 11
0
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))
Ejemplo n.º 13
0
#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(
Ejemplo n.º 14
0
# 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.)
Ejemplo n.º 15
0
    	      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")
    ))
Ejemplo n.º 16
0
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())
Ejemplo n.º 19
0
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):
Ejemplo n.º 22
0
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
Ejemplo n.º 23
0
#     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
Ejemplo n.º 24
0
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' )
Ejemplo n.º 26
0
 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)
Ejemplo n.º 28
0
    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)
Ejemplo n.º 29
0
    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() 
Ejemplo n.º 30
0
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)
Ejemplo n.º 31
0
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),
                         )
Ejemplo n.º 32
0
    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:
Ejemplo n.º 33
0
 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
 ),
Ejemplo n.º 34
0
# 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)
Ejemplo n.º 35
0
    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,
Ejemplo n.º 36
0
	    "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),
Ejemplo n.º 37
0
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)
                            )
Ejemplo n.º 38
0
)

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)
Ejemplo n.º 39
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.
Ejemplo n.º 41
0
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)
)

Ejemplo n.º 42
0
#                                 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),
Ejemplo n.º 43
0
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*
Ejemplo n.º 44
0
    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(
Ejemplo n.º 45
0
                                      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)

Ejemplo n.º 46
0
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),  
Ejemplo n.º 48
0
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"
Ejemplo n.º 49
0
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'),
Ejemplo n.º 51
0
## ---- 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
Ejemplo n.º 52
0
        "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))
Ejemplo n.º 53
0
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)

Ejemplo n.º 54
0
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(
Ejemplo n.º 58
0
        "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',
Ejemplo n.º 59
0
    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",
Ejemplo n.º 60
0
    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#####################################