test/expert/runMuonOverlap_NN_Regression_FP.py

0001 # -*- coding: utf-8 -*-
0002 import FWCore.ParameterSet.Config as cms
0003 process = cms.Process("L1TMuonEmulation")
0004 import os
0005 import sys
0006
0007 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0008
0009 verbose = True
0010
0011 version = "_p2DT_v4"
0012
0013 if verbose:
0014     process.MessageLogger = cms.Service("MessageLogger",
0015        #suppressInfo       = cms.untracked.vstring('AfterSource', 'PostModule'),
0016        destinations   = cms.untracked.vstring(
0017                                                #'detailedInfo',
0018                                                #'critical',
0019                                                #'cout',
0020                                                #'cerr',
0021                                                'omtfEventPrint'
0022                     ),
0023        categories        = cms.untracked.vstring('l1tOmtfEventPrint', 'OMTFReconstruction'),
0024        omtfEventPrint = cms.untracked.PSet(
0025                          filename  = cms.untracked.string('log_MuonOverlap_nn' + version),
0026                          extension = cms.untracked.string('.txt'),
0027                          threshold = cms.untracked.string('DEBUG'),
0028                          default = cms.untracked.PSet( limit = cms.untracked.int32(0) ),
0029                          #INFO   =  cms.untracked.int32(0),
0030                          #DEBUG   = cms.untracked.int32(0),
0031                          l1tOmtfEventPrint = cms.untracked.PSet( limit = cms.untracked.int32(1000000000) ),
0032                          OMTFReconstruction = cms.untracked.PSet( limit = cms.untracked.int32(1000000000) )
0033                        ),
0034        debugModules = cms.untracked.vstring('L1MuonAnalyzerOmtf', 'simOmtfPhase2Digis')
0035        #debugModules = cms.untracked.vstring('*')
0036     )
0037
0038     #process.MessageLogger.cerr.FwkReport.reportEvery = cms.untracked.int32(100)
0039 if not verbose:
0040     process.MessageLogger.cerr.FwkReport.reportEvery = cms.untracked.int32(100)
0041     process.options = cms.untracked.PSet(wantSummary = cms.untracked.bool(False),
0042                                          #SkipEvent = cms.untracked.vstring('ProductNotFound')
0043                                      )
0044
0045 process.source = cms.Source('PoolSource',
0046  #fileNames = cms.untracked.vstring('file:/afs/cern.ch/work/g/gflouris/public/SingleMuPt6180_noanti_10k_eta1.root')
0047  #fileNames = cms.untracked.vstring('file:///afs/cern.ch/work/k/kbunkow/private/omtf_data/SingleMu_15_p_1_1_qtl.root')
0048  #fileNames = cms.untracked.vstring('file:///eos/user/k/kbunkow/cms_data/mc/PhaseIIFall17D/SingleMu_PU200_32DF01CC-A342-E811-9FE7-48D539F3863E_dump500Events.root')
0049  fileNames = cms.untracked.vstring("file:///eos/user/k/kbunkow/cms_data/mc/PhaseIITDRSpring19DR/PhaseIITDRSpring19DR_Mu_FlatPt2to100_noPU_v31_dump4000Ev.root")
0050  )
0051
0052 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(500))
0053
0054 # import of standard configurations
0055 process.load('Configuration.StandardSequences.Services_cff')
0056 #process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0057 #process.load('FWCore.MessageService.MessageLogger_cfi')
0058 #process.load('Configuration.EventContent.EventContent_cff')
0059 #process.load('SimGeneral.MixingModule.mixNoPU_cfi')
0060 process.load('Configuration.Geometry.GeometryExtended2026D95Reco_cff')
0061 process.load('Configuration.StandardSequences.MagneticField_cff')
0062 #process.load('Configuration.StandardSequences.RawToDigi_cff')
0063 #process.load('Configuration.StandardSequences.SimL1Emulator_cff')
0064 #process.load('Configuration.StandardSequences.SimPhase2L1GlobalTriggerEmulator_cff')
0065 #process.load('L1Trigger.Configuration.Phase2GTMenus.SeedDefinitions.prototypeSeeds')
0066 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0067 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0068
0069 from Configuration.AlCa.GlobalTag import GlobalTag
0070 process.GlobalTag = GlobalTag(process.GlobalTag, '131X_mcRun4_realistic_v9', '')
0071
0072
0073
0074 process.TFileService = cms.Service("TFileService", fileName = cms.string('omtfAnalysis1.root'), closeFileFast = cms.untracked.bool(True) )
0075
0076
0077 #Calibrate Digis
0078 process.load("L1Trigger.DTTriggerPhase2.CalibratedDigis_cfi")
0079 #process.CalibratedDigis.flat_calib = 325 #turn to 0 to use the DB  , 325 for JM and Jorge benchmark
0080 process.CalibratedDigis.dtDigiTag = "simMuonDTDigis" #turn to 0 to use the DB  , 325 for JM and Jorge benchmark
0081 process.CalibratedDigis.scenario = 0 # 0 for mc, 1 for data, 2 for slice test
0082
0083 #DTTriggerPhase2
0084 process.load("L1Trigger.DTTriggerPhase2.dtTriggerPhase2PrimitiveDigis_cfi")
0085 #process.dtTriggerPhase2PrimitiveDigis.trigger_with_sl = 3  #4 means SL 1 and 3
0086 #for the moment the part working in phase2 format is the slice test
0087 #process.dtTriggerPhase2PrimitiveDigis.p2_df = True
0088 #process.dtTriggerPhase2PrimitiveDigis.filter_primos = True
0089 #for debugging
0090 #process.dtTriggerPhase2PrimitiveDigis.pinta = True
0091 #process.dtTriggerPhase2PrimitiveDigis.min_phinhits_match_segment = 4
0092 #process.dtTriggerPhase2PrimitiveDigis.debug = True
0093 process.dtTriggerPhase2PrimitiveDigis.scenario = 0
0094 process.dtTriggerPhase2PrimitiveDigis.dump = True
0095
0096 ####OMTF Emulator
0097
0098 process.load('L1Trigger.L1TMuonOverlapPhase2.simOmtfPhase2Digis_cfi')
0099
0100 process.simOmtfPhase2Digis.dumpResultToXML = cms.bool(True)
0101 process.simOmtfPhase2Digis.XMLDumpFileName = cms.string("TestEvents_" + version + ".xml")
0102
0103 process.simOmtfPhase2Digis.eventCaptureDebug = cms.bool(True)
0104
0105 process.simOmtfPhase2Digis.rpcMaxClusterSize = cms.int32(3)
0106 process.simOmtfPhase2Digis.rpcMaxClusterCnt = cms.int32(2)
0107 process.simOmtfPhase2Digis.rpcDropAllClustersIfMoreThanMax = cms.bool(True)
0108
0109 process.simOmtfPhase2Digis.lctCentralBx = cms.int32(8);#<<<<<<<<<<<<<<<<!!!!!!!!!!!!!!!!!!!!TODO this was changed in CMSSW 10(?) to 8. if the data were generated with the previous CMSSW then you have to use 6
0110
0111 process.simOmtfPhase2Digis.dropDTPrimitives = cms.bool(True)
0112 process.simOmtfPhase2Digis.usePhase2DTPrimitives = cms.bool(True)
0113
0114 process.simOmtfPhase2Digis.dropRPCPrimitives = cms.bool(True)     #TODO!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
0115 process.simOmtfPhase2Digis.dropCSCPrimitives = cms.bool(True)   #TODO!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
0116
0117 #process.simOmtfPhase2Digis.patternsXMLFile = cms.FileInPath("L1Trigger/L1TMuon/data/omtf_config/Patterns_0x0003.xml")
0118
0119 process.simOmtfPhase2Digis.neuralNetworkFile = cms.FileInPath("L1Trigger/L1TMuon/data/omtf_config/lutNN_omtfRegression_FP_6.xml")
0120
0121 process.dumpED = cms.EDAnalyzer("EventContentAnalyzer")
0122 process.dumpES = cms.EDAnalyzer("PrintEventSetupContent")
0123
0124 process.L1TMuonSeq = cms.Sequence( process.simOmtfPhase2Digis
0125                                    #+ process.dumpED
0126                                    #+ process.dumpES
0127 )
0128
0129 process.L1TMuonPath = cms.Path(process.CalibratedDigis *
0130                             process.dtTriggerPhase2PrimitiveDigis *
0131                             process.L1TMuonSeq)
0132
0133 process.out = cms.OutputModule("PoolOutputModule",
0134    fileName = cms.untracked.string("l1tomtf_superprimitives1.root")
0135 )
0136
0137 #process.output_step = cms.EndPath(process.out)
0138 #process.schedule = cms.Schedule(process.L1TMuonPath)
0139 #process.schedule.extend([process.output_step])