L1TMuonOverlapPhase1/test/runMuonOverlap_2023withExtrapol_mc.py

0001 # -*- coding: utf-8 -*-
0002 import FWCore.ParameterSet.Config as cms
0003 process = cms.Process("L1TMuonEmulation")
0004 import os
0005 import sys
0006 import re
0007 from os import listdir
0008 from os.path import isfile, join
0009 import fnmatch
0010
0011 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0012
0013 verbose = True
0014 #version = 't14_extrapolSimpl_displ_allfiles'
0015 #version = 't16_extrapolSimpl_displ_test'
0016 #version = 'ExtraplMB1nadMB2SimplifiedFP_t17_v11_test_valueP1Scale'
0017 #version = 'ExtraplMB1nadMB2SimplifiedFP_t19_v16_test_bits'
0018 #version = 'Patterns_0x00012_t17_v11_extr_off_test_bits'
0019 #version = 'ExtraplMB1nadMB2SimplifiedFP_t19_v16_test_bits_MH-1000_MFF-150_CTau-1000mm'
0020 version = 'ExtraplMB1nadMB2SimplifiedFP_t19_v16_test'
0021
0022 if verbose:
0023     process.MessageLogger = cms.Service("MessageLogger",
0024        #suppressInfo       = cms.untracked.vstring('AfterSource', 'PostModule'),
0025        destinations   = cms.untracked.vstring(
0026                                                #'detailedInfo',
0027                                                #'critical',
0028                                                #'cout',
0029                                                #'cerr',
0030                                                'omtfEventPrint'
0031                     ),
0032        categories        = cms.untracked.vstring('l1tOmtfEventPrint', 'OMTFReconstruction'),
0033        omtfEventPrint = cms.untracked.PSet(
0034                          filename  = cms.untracked.string('log_MuonOverlap_newPats_t' + version),
0035                          extension = cms.untracked.string('.txt'),
0036                          threshold = cms.untracked.string('INFO'),
0037                          default = cms.untracked.PSet( limit = cms.untracked.int32(0) ),
0038                          #INFO   =  cms.untracked.int32(0),
0039                          #DEBUG   = cms.untracked.int32(0),
0040                          l1tOmtfEventPrint = cms.untracked.PSet( limit = cms.untracked.int32(1000000000) ),
0041                          OMTFReconstruction = cms.untracked.PSet( limit = cms.untracked.int32(1000000000) )
0042                        ),
0043        debugModules = cms.untracked.vstring('L1MuonAnalyzerOmtf', 'simOmtfDigis')
0044        #debugModules = cms.untracked.vstring('*')
0045     )
0046
0047     #process.MessageLogger.cerr.FwkReport.reportEvery = cms.untracked.int32(100)
0048 if not verbose:
0049     process.MessageLogger.cerr.FwkReport.reportEvery = cms.untracked.int32(1000)
0050     process.options = cms.untracked.PSet(wantSummary = cms.untracked.bool(False),
0051                                          #SkipEvent = cms.untracked.vstring('ProductNotFound')
0052                                      )
0053
0054
0055 # PostLS1 geometry used
0056 process.load('Configuration.Geometry.GeometryExtended2015Reco_cff')
0057 process.load('Configuration.Geometry.GeometryExtended2015_cff')
0058 ############################
0059 #process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_condDBv2_cff')
0060 #from Configuration.AlCa.GlobalTag_condDBv2 import GlobalTag
0061 #process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:run2_mc', '')
0062
0063
0064 # import of standard configurations
0065 process.load('Configuration.StandardSequences.Services_cff')
0066 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0067 process.load('Configuration.EventContent.EventContent_cff')
0068 process.load('SimGeneral.MixingModule.mixNoPU_cfi')
0069 #process.load('Configuration.Geometry.GeometryExtended2023D41Reco_cff')
0070 #process.load('Configuration.Geometry.GeometryExtended2023D41_cff')
0071 process.load('Configuration.StandardSequences.MagneticField_cff')
0072 #process.load('Configuration.StandardSequences.SimL1Emulator_cff')
0073 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0074 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0075
0076 from Configuration.AlCa.GlobalTag import GlobalTag
0077 #process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:upgradePLS3', '')
0078 process.GlobalTag = GlobalTag(process.GlobalTag, '103X_upgrade2023_realistic_v2', '')
0079
0080
0081 process.source = cms.Source('PoolSource',
0082  #fileNames = cms.untracked.vstring('file:/afs/cern.ch/work/g/gflouris/public/SingleMuPt6180_noanti_10k_eta1.root')
0083  #fileNames = cms.untracked.vstring('file:///afs/cern.ch/work/k/kbunkow/private/omtf_data/SingleMu_15_p_1_1_qtl.root')
0084  #fileNames = cms.untracked.vstring('file:///afs/cern.ch/work/k/kbunkow/private/omtf_data/SingleMu_20_p_100_2_B7Z.root')
0085  #fileNames = cms.untracked.vstring('file:///afs/cern.ch/work/k/kbunkow/private/omtf_data/SingleMu_5_p_1_1_Meh.root')
0086  #fileNames = cms.untracked.vstring('file:///afs/cern.ch/work/k/kbunkow/private/omtf_data/SingleMu_7_p_1_1_DkC.root')
0087  #fileNames = cms.untracked.vstring('file:///afs/cern.ch/work/k/kbunkow/private/omtf_data/SingleMu_18_p_1_1_2KD.root')
0088  #fileNames = cms.untracked.vstring('file:///afs/cern.ch/work/k/kbunkow/public/CMSSW/cmssw_10_x_x_l1tOfflinePhase2/CMSSW_10_6_0_pre4/src/L1Trigger/L1TMuonBayes/test/expert/DisplacedMuonGun_Pt30To100_Dxy_0_1000_E68C6334-7F62-E911-8AA5-0025905B8610_dump2000Ev.root')
0089
0090  #fileNames = cms.untracked.vstring('file:///eos/user/k/kbunkow/cms_data/SingleMuFullEta/721_FullEta_v4/SingleMu_31_p_101_2_DzU.root')
0091  #fileNames = cms.untracked.vstring('file:///eos/user/a/akalinow/Data/SingleMu/9_3_14_FullEta_v2/SingleMu_6_p_1.root'),
0092  #fileNames = cms.untracked.vstring("file:///eos/user/k/kbunkow/cms_data/mc/PhaseIITDRSpring19DR/PhaseIITDRSpring19DR_Mu_FlatPt2to100_noPU_v31_E0D5C6A5-B855-D14F-9124-0B2C9B28D0EA_dump4000Ev.root"),
0093  fileNames = cms.untracked.vstring("file:///eos/user/k/kbunkow/cms_data/mc/PhaseIITDRSpring19DR/JPsiToMuMu_Pt0to100_NoPU_FDA71CB6-4C3B-4540-99EB-803077C6EC2D_dump4000Ev.root",
0094                                    "file:///eos/user/k/kbunkow/cms_data/mc/PhaseIITDRSpring19DR/PhaseIITDRSpring19DR_Mu_FlatPt2to100_noPU_v31_E0D5C6A5-B855-D14F-9124-0B2C9B28D0EA_dump4000Ev.root"
0095                                    ),
0096  #fileNames = cms.untracked.vstring("file:///eos/user/k/kbunkow/cms_data/mc/PhaseIITDRSpring19DR/HSCPppstau_M_871_PU200_v3-v2_1ADE9D9E-8C0C-1948-A405-5DFDA1AF5172_dump100Ev.root"),
0097  #fileNames = cms.untracked.vstring("file:///eos/user/k/kbunkow/cms_data/mc/PhaseIIFall17D/SingleMu_PU200_32DF01CC-A342-E811-9FE7-48D539F3863E_dump500Events.root"),
0098  #fileNames = cms.untracked.vstring("file:///eos/user/k/kbunkow/cms_data/mc/PhaseIITDRSpring19DR/Nu_E10-pythia8-gun_PU250_v3_ext2-v1_FFE07316-3810-6E44-97A1-5753A3070D12_dump100Ev.root"),
0099  #fileNames = cms.untracked.vstring("file:///afs/cern.ch/work/k/kbunkow/public/CMSSW/cmssw_10_x_x_l1tOfflinePhase2/CMSSW_10_6_1_patch2/src/L1Trigger/L1TMuonBayes/test/expert/Nu_E10-pythia8-gun_PU250_v3_ext2-v1_FFE07316-3810-6E44-97A1-5753A3070D12_dump100Ev.root"),
0100
0101
0102
0103
0104  #fileNames = cms.untracked.vstring('file:///eos/home-k/konec/FFCFF986-ED0B-B74F-B253-C511D19B8249.root'),
0105  #fileNames = cms.untracked.vstring('file:///afs/cern.ch/user/k/konec/work/CMSSW_10_6_1_patch2.displaced/src/UserCode/OmtfAnalysis/jobs/FFCFF986-ED0B-B74F-B253-C511D19B8249.root'),
0106
0107
0108         inputCommands=cms.untracked.vstring(
0109         'keep *',
0110         'drop l1tEMTFHit2016Extras_simEmtfDigis_CSC_HLT',
0111         'drop l1tEMTFHit2016Extras_simEmtfDigis_RPC_HLT',
0112         'drop l1tEMTFHit2016s_simEmtfDigis__HLT',
0113         'drop l1tEMTFTrack2016Extras_simEmtfDigis__HLT',
0114         'drop l1tEMTFTrack2016s_simEmtfDigis__HLT')
0115 )
0116
0117 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1))
0118
0119
0120 ####Event Setup Producer
0121 process.load('L1Trigger.L1TMuonOverlapPhase1.fakeOmtfParams_cff')
0122 process.omtfParams.configXMLFile = cms.FileInPath("L1Trigger/L1TMuon/data/omtf_config/hwToLogicLayer_0x0009.xml")
0123
0124 # reding config and patterns from the DB
0125 process.esProd = cms.EDAnalyzer("EventSetupRecordDataGetter",
0126    toGet = cms.VPSet(
0127       cms.PSet(record = cms.string('L1TMuonOverlapParamsRcd'),
0128                data = cms.vstring('L1TMuonOverlapParams'))
0129                    ),
0130    verbose = cms.untracked.bool(False)
0131 )
0132
0133
0134 ####OMTF Emulator
0135 process.load('L1Trigger.L1TMuonOverlapPhase1.simOmtfDigis_extrapolSimple_cfi')
0136
0137 process.simOmtfDigis.patternsXMLFile = cms.FileInPath("L1Trigger/L1TMuonOverlapPhase1/test/expert/omtf/Patterns_ExtraplMB1nadMB2SimplifiedFP_t17_classProb17_recalib2.xml")
0138
0139 process.simOmtfDigis.dumpResultToXML = cms.bool(True)
0140 process.simOmtfDigis.eventCaptureDebug = cms.bool(True)
0141
0142 process.simOmtfDigis.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
0143
0144 #process.simOmtfDigis.rpcSimHitsInputTag = cms.InputTag("g4SimHits", "MuonRPCHits")
0145
0146 #process.dumpED = cms.EDAnalyzer("EventContentAnalyzer")
0147 #process.dumpES = cms.EDAnalyzer("PrintEventSetupContent")
0148
0149
0150 process.L1TMuonSeq = cms.Sequence( #process.esProd  +
0151                                    process.simOmtfDigis
0152                                    #+ process.dumpED
0153                                    #+ process.dumpES
0154 )
0155
0156 process.L1TMuonPath = cms.Path(process.L1TMuonSeq)
0157
0158 #process.load("L1Trigger.TrackFindingTracklet.L1TrackletTracks_cff")
0159 #process.TTTracks = cms.Path(process.L1TrackletTracks)
0160 #process.TTTracksWithTruth = cms.Path(process.L1TrackletTracksWithAssociators)
0161
0162 process.schedule = cms.Schedule(process.L1TMuonPath)
0163
0164 #process.out = cms.OutputModule("PoolOutputModule",
0165 #   fileName = cms.untracked.string("l1tomtf_superprimitives1.root")
0166 #)
0167
0168 #process.output_step = cms.EndPath(process.out)
0169 #process.schedule = cms.Schedule(process.L1TMuonPath)
0170 #process.schedule.extend([process.output_step])