HGCalValidation/python/HGCalValidator_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 from Validation.HGCalValidation.CaloParticleSelectionForEfficiency_cfi import *
0004 from Validation.HGCalValidation.HGVHistoProducerAlgoBlock_cfi import *
0005
0006 from SimCalorimetry.HGCalAssociatorProducers.LCToCPAssociation_cfi import layerClusterCaloParticleAssociation
0007 from SimCalorimetry.HGCalAssociatorProducers.LCToSCAssociation_cfi import layerClusterSimClusterAssociation
0008
0009 from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
0010
0011 from RecoHGCal.TICL.iterativeTICL_cff import ticlIterLabels, ticlIterLabelsMerge
0012
0013 labelTst = [cms.InputTag("ticlTracksters"+iteration) for iteration in ticlIterLabelsMerge]
0014 labelTst.extend([cms.InputTag("ticlSimTracksters", "fromCPs"), cms.InputTag("ticlSimTracksters")])
0015 lcInputMask = [cms.InputTag("ticlTracksters"+iteration) for iteration in ticlIterLabels]
0016 lcInputMask.extend([cms.InputTag("ticlSimTracksters", "fromCPs"), cms.InputTag("ticlSimTracksters")])
0017 hgcalValidator = DQMEDAnalyzer(
0018     "HGCalValidator",
0019
0020     ### general settings ###
0021     # selection of CP for evaluation of efficiency #
0022     CaloParticleSelectionForEfficiency,
0023
0024     ### reco input configuration ###
0025     #2DLayerClusters, PFClusters, Tracksters
0026     label_lcl = layerClusterCaloParticleAssociation.label_lc,
0027     label_tst = cms.VInputTag(labelTst),
0028     label_simTS = cms.InputTag("ticlSimTracksters"),
0029     label_simTSFromCP = cms.InputTag("ticlSimTracksters", "fromCPs"),
0030
0031     associator = cms.untracked.InputTag("layerClusterCaloParticleAssociationProducer"),
0032
0033     associatorSim = cms.untracked.InputTag("layerClusterSimClusterAssociationProducer"),
0034
0035     #General info on layers etc.
0036     SaveGeneralInfo = cms.untracked.bool(True),
0037     #CaloParticle related plots
0038     doCaloParticlePlots = cms.untracked.bool(True),
0039     #Select caloParticles for efficiency or pass through
0040     doCaloParticleSelection = cms.untracked.bool(True),
0041     #SimCluster related plots
0042     doSimClustersPlots = cms.untracked.bool(True),
0043     label_SimClusters = cms.InputTag("SimClusters"),
0044     label_SimClustersLevel = cms.InputTag("ClusterLevel"),
0045     #Layer Cluster related plots
0046     doLayerClustersPlots = cms.untracked.bool(True),
0047     label_layerClusterPlots = cms.InputTag("hgcalLayerClusters"),
0048     label_LCToCPLinking = cms.InputTag("LCToCP_association"),
0049     #Trackster related plots
0050     doTrackstersPlots = cms.untracked.bool(True),
0051     label_TS = cms.string("Morphology"),
0052     label_TSToCPLinking = cms.string("TSToCP_linking"),
0053     label_TSToSTSPR = cms.string("TSToSTS_patternRecognition"),
0054
0055     #The cumulative material budget in front of each layer. To be more specific, it
0056     #is the material budget just in front of the active material (not including it).
0057     #This file is created using the official material budget code.
0058     cummatbudinxo = cms.FileInPath('Validation/HGCalValidation/data/D41.cumulative.xo'),
0059
0060     ### sim input configuration ###
0061     label_cp_effic = layerClusterCaloParticleAssociation.label_cp,
0062     label_cp_fake = cms.InputTag("mix","MergedCaloTruth"),
0063     #simClusters
0064     label_scl = layerClusterSimClusterAssociation.label_scl,
0065
0066     simVertices = cms.InputTag("g4SimHits"),
0067
0068     LayerClustersInputMask = cms.VInputTag(lcInputMask),
0069
0070     #Total number of layers of HGCal that we want to monitor
0071     #Could get this also from HGCalImagingAlgo::maxlayer but better to get it from here
0072     totallayers_to_monitor = cms.int32(52),
0073     #Thicknesses we want to monitor. -1 is for scintillator
0074     thicknesses_to_monitor = cms.vint32(120,200,300,-1),
0075
0076     # HistoProducerAlgo. Defines the set of plots to be booked and filled
0077     histoProducerAlgoBlock = HGVHistoProducerAlgoBlock,
0078
0079     ### output configuration
0080     dirName = cms.string('HGCAL/HGCalValidator/')
0081
0082 )
0083
0084 from Configuration.ProcessModifiers.premix_stage2_cff import premix_stage2
0085 premix_stage2.toModify(hgcalValidator,
0086     label_cp_fake = "mixData:MergedCaloTruth"
0087 )
0088
0089 from Configuration.Eras.Modifier_phase2_hgcalV10_cff import phase2_hgcalV10
0090 phase2_hgcalV10.toModify(hgcalValidator, totallayers_to_monitor = cms.int32(50))
0091
0092 from Configuration.Eras.Modifier_phase2_hgcalV16_cff import phase2_hgcalV16
0093 phase2_hgcalV16.toModify(hgcalValidator, totallayers_to_monitor = cms.int32(47))