Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2023-10-25 09:55:18

0001 from __future__ import absolute_import
0002 import FWCore.ParameterSet.Config as cms
0003 
0004 import SimCalorimetry.HGCalSimProducers.hgcalDigitizer_cfi as digiparam
0005 import RecoLocalCalo.HGCalRecProducers.HGCalUncalibRecHit_cfi as recoparam
0006 import RecoLocalCalo.HGCalRecProducers.HGCalRecHit_cfi as recocalibparam
0007 import L1Trigger.L1THGCal.hgcalLayersCalibrationCoefficients_cfi as layercalibparam
0008 
0009 
0010 fcPerMip = recoparam.HGCalUncalibRecHit.HGCEEConfig.fCPerMIP
0011 keV2fC = digiparam.hgceeDigitizer.digiCfg.keV2fC
0012 layerWeights = layercalibparam.TrgLayer_dEdX_weights
0013 thicknessCorrections = recocalibparam.HGCalRecHit.thicknessCorrection
0014 
0015 ntuple_event = cms.PSet(
0016     NtupleName = cms.string('HGCalTriggerNtupleEvent')
0017 )
0018 
0019 from FastSimulation.Event.ParticleFilter_cfi import ParticleFilterBlock
0020 PartFilterConfig = ParticleFilterBlock.ParticleFilter.copy()
0021 PartFilterConfig.protonEMin = cms.double(100000)
0022 PartFilterConfig.etaMax = cms.double(3.1)
0023 
0024 ntuple_gen = cms.PSet(
0025     NtupleName = cms.string('HGCalTriggerNtupleGen'),
0026     GenParticles = cms.InputTag('genParticles'),
0027     GenPU = cms.InputTag('addPileupInfo'),
0028     MCEvent = cms.InputTag('generatorSmeared'),
0029     SimTracks = cms.InputTag('g4SimHits'),
0030     SimVertices = cms.InputTag('g4SimHits'),
0031     particleFilter = PartFilterConfig
0032 )
0033 
0034 ntuple_gentau = cms.PSet(
0035     NtupleName = cms.string('HGCalTriggerNtupleGenTau'),
0036     GenParticles = cms.InputTag('genParticles'),
0037     isPythia8 = cms.bool(False)
0038 )
0039 
0040 ntuple_genjet = cms.PSet(
0041     NtupleName = cms.string('HGCalTriggerNtupleGenJet'),
0042     GenJets = cms.InputTag('ak4GenJetsNoNu')
0043 )
0044 
0045 ntuple_digis = cms.PSet(
0046     NtupleName = cms.string('HGCalTriggerNtupleHGCDigis'),
0047     HGCDigisEE = cms.InputTag('simHGCalUnsuppressedDigis:EE'),
0048     HGCDigisFH = cms.InputTag('simHGCalUnsuppressedDigis:HEfront'),
0049     HGCDigisBH = cms.InputTag('simHGCalUnsuppressedDigis:HEback'),
0050     eeSimHits = cms.InputTag('g4SimHits:HGCHitsEE'),
0051     fhSimHits = cms.InputTag('g4SimHits:HGCHitsHEfront'),
0052     bhSimHits = cms.InputTag('g4SimHits:HGCHitsHEback'),
0053     isSimhitComp = cms.bool(False),
0054     digiBXselect = cms.vuint32(2)
0055 )
0056 
0057 ntuple_triggercells = cms.PSet(
0058     NtupleName = cms.string('HGCalTriggerNtupleHGCTriggerCells'),
0059     TriggerCells = cms.InputTag('l1tHGCalConcentratorProducer:HGCalConcentratorProcessorSelection'),
0060     Multiclusters = cms.InputTag('l1tHGCalBackEndLayer2Producer:HGCalBackendLayer2Processor3DClustering'),
0061     eeSimHits = cms.InputTag('g4SimHits:HGCHitsEE'),
0062     fhSimHits = cms.InputTag('g4SimHits:HGCHitsHEfront'),
0063     bhSimHits = cms.InputTag('g4SimHits:HGCHitsHEback'), 
0064     FillSimEnergy = cms.bool(False),
0065     FillTruthMap = cms.bool(False),
0066     fcPerMip = fcPerMip,
0067     keV2fC = keV2fC,
0068     layerWeights = layerWeights,
0069     thicknessCorrections = thicknessCorrections,
0070     FilterCellsInMulticlusters = cms.bool(False)
0071 )
0072 
0073 ntuple_triggersums = cms.PSet(
0074     NtupleName = cms.string('HGCalTriggerNtupleHGCTriggerSums'),
0075     TriggerSums = cms.InputTag('l1tHGCalConcentratorProducer:HGCalConcentratorProcessorSelection'),
0076 )
0077 
0078 ntuple_econdata = cms.PSet(
0079     NtupleName = cms.string('HGCalTriggerNtupleHGCConcentratorData'),
0080     ConcentratorData = cms.InputTag('l1tHGCalConcentratorProducer:HGCalConcentratorProcessorSelection'),
0081 )
0082 
0083 ntuple_clusters = cms.PSet(
0084     NtupleName = cms.string('HGCalTriggerNtupleHGCClusters'),
0085     Clusters = cms.InputTag('l1tHGCalBackEndLayer1Producer:HGCalBackendLayer1Processor2DClustering'),
0086     Multiclusters = cms.InputTag('l1tHGCalBackEndLayer2Producer:HGCalBackendLayer2Processor3DClustering'),
0087     FilterClustersInMulticlusters = cms.bool(False)
0088 )
0089 
0090 from L1Trigger.L1THGCal.egammaIdentification import egamma_identification_histomax
0091 ntuple_multiclusters = cms.PSet(
0092     NtupleName = cms.string('HGCalTriggerNtupleHGCMulticlusters'),
0093     Multiclusters = cms.InputTag('l1tHGCalBackEndLayer2Producer:HGCalBackendLayer2Processor3DClustering'),
0094     EGIdentification = egamma_identification_histomax.clone(),
0095     FillLayerInfo = cms.bool(False),
0096     FillInterpretationInfo = cms.bool(True)
0097 )
0098 
0099 ntuple_towers = cms.PSet(
0100     NtupleName = cms.string('HGCalTriggerNtupleHGCTowers'),
0101     Towers = cms.InputTag('l1tHGCalTowerProducer:HGCalTowerProcessor')
0102 )
0103 
0104 l1tHGCalTriggerNtuplizer = cms.EDAnalyzer(
0105     "HGCalTriggerNtupleManager",
0106     Ntuples = cms.VPSet(
0107         ntuple_event,
0108         ntuple_gen,
0109         ntuple_genjet,
0110         ntuple_gentau,
0111         ntuple_digis,
0112         ntuple_triggercells,
0113         ntuple_triggersums,
0114         ntuple_multiclusters,
0115         ntuple_towers
0116     )
0117 )