Muon L1 emulation high rates in heavy ion

Related Twiki to HI Dilepton group

MC Samples

  • CMSSW 362 (up to DIGI)
    • /castor/cern.ch/user/d/dmoon/cms362/DIGI/
  • CMSSW 370 (up to DIGI)
    • /castor/cern.ch/user/d/dmoon/cms370/DIGI/
  • CMSSW 381 (up to HLTDEBUG)
    • /store/relval/CMSSW_3_8_1/RelValHydjetQ_MinBias_2760GeV/GEN-SIM-DIGI-RAW-HLTDEBUG/MC_38Y_V8-v1/0011/

Summary of the rates

In heavy ion simulations with the uptodate csctf emulator

trigger_on_ME1a = cms.bool(False),
trigger_on_ME1b = cms.bool(False),

initializeFromPSet = cms.bool(True)

    • to make sure one gets the latest tags :
cvs co -r HEAD DataFormats/L1CSCTrackFinder
cvs co -r V00-19-04 L1Trigger/CSCTrackFinder
//which has to go together with
cvs co -r V02-12-00 L1Trigger/DTTrackFinder 
cvs co -r V01-00-06 L1TriggerConfig/CSCTFConfigProducers (just to be sure)

Algorithm Key Efficiency Stat error (%)
L1_DoubleMu3 1.80 +- 0.59
L1_DoubleMuOpen 2.20 +- 0.66
L1_SingleMu0 5.80 +- 1.05
L1_SingleMu10 0.60 +- 0.35
L1_SingleMu14 0.40 +- 0.28
L1_SingleMu20 0.20 +- 0.20
L1_SingleMu3 4.40 +- 0.92
L1_SingleMu5 2.00 +- 0.63
L1_SingleMu7 1.20 +- 0.49
L1_SingleMuOpen 10.40 +- 1.37

In heavy ion simulations : digging into high rates

  • Reported by Christof and Frank, and including the other options (such as Without CSC, DT)

CMSSW 362     370     381   Without CSC 362     370     381   Without DT 362     370     381  
L1_DoubleMu3 2.60 +- 0.71   17.78 +- 1.72   15.80 +- 1.63   0.00 +- 0.00   0.20 +- 0.20   0.00 +- 0.00   2.40 +- 0.68   18.99 +- 1.76      
L1_DoubleMuOpen 3.60 +- 0.83   18.59 +- 1.75   16.40 +- 1.66   0.20 +- 0.20   0.20 +- 0.20   0.00 +- 0.00   3.20 +- 0.79   19.60 +- 1.78      
L1_SingleMu0 7.60 +- 1.19   31.52 +- 2.09   5.60 +- 1.03   2.40 +- 0.68   1.01 +- 0.45   1.00 +- 0.44   6.80 +- 1.13   32.32 +- 2.10      
L1_SingleMu10 0.60 +- 0.35   30.51 +- 2.07   0.40 +- 0.28   0.40 +- 0.28   0.20 +- 0.20   0.00 +- 0.00   0.60 +- 0.35   31.11 +- 2.08      
L1_SingleMu14 0.60 +- 0.35   30.51 +- 2.07   0.20 +- 0.20   0.40 +- 0.28   0.00 +- 0.00   0.00 +- 0.00   0.60 +- 0.35   31.11 +- 2.08      
L1_SingleMu20 0.00 +- 0.00   30.51 +- 2.07   0.20 +- 0.20   0.20 +- 0.20   0.00 +- 0.00   0.00 +- 0.00   0.00 +- 0.00   31.11 +- 2.08      
L1_SingleMu3 6.20 +- 1.08   31.31 +- 2.08   4.20 +- 0.90   2.40 +- 0.68   1.01 +- 0.45   1.00 +- 0.44   5.80 +- 1.05   32.12 +- 2.10      
L1_SingleMu5 2.60 +- 0.71   30.91 +- 2.08   1.80 +- 0.59   0.60 +- 0.35   0.61 +- 0.35   0.60 +- 0.35   2.00 +- 0.63   31.72 +- 2.09      
L1_SingleMu7 1.20 +- 0.49   30.51 +- 2.07   1.00 +- 0.44   0.40 +- 0.28   0.20 +- 0.20   0.00 +- 0.00   1.20 +- 0.49   31.31 +- 2.08      
L1_SingleMuOpen 10.40 +- 1.37   33.94 +- 2.13   31.80 +- 2.08   3.20 +- 0.79   1.01 +- 0.45   1.40 +- 0.53   10.40 +- 1.37   34.55 +- 2.14      

  • None option test : CSC "none" and DT "none"

CSC NONE 362     370     381   DT NONE 362     370     381  
L1_DoubleMu3 0.00 +- 0.00   0.00 +- 0.00   0.00 +- 0.00   2.40 +- 0.68   19.75 +- 1.82   15.60 +- 1.62
L1_DoubleMuOpen 0.20 +- 0.20   0.00 +- 0.00   0.00 +- 0.00   3.20 +- 0.79   20.38 +- 1.85   16.20 +- 1.65
L1_SingleMu0 2.40 +- 0.68   1.05 +- 0.47   1.00 +- 0.44   7.00 +- 1.14   32.14 +- 2.14   5.40 +- 1.01
L1_SingleMu10 0.40 +- 0.28   0.21 +- 0.21   0.00 +- 0.00   0.60 +- 0.35   31.30 +- 2.13   0.40 +- 0.28
L1_SingleMu14 0.40 +- 0.28   0.00 +- 0.00   0.00 +- 0.00   0.60 +- 0.35   31.30 +- 2.13   0.20 +- 0.20
L1_SingleMu20 0.20 +- 0.20   0.00 +- 0.00   0.00 +- 0.00   0.00 +- 0.00   31.30 +- 2.13   0.20 +- 0.20
L1_SingleMu3 2.40 +- 0.68   1.05 +- 0.47   1.00 +- 0.44   6.00 +- 1.06   32.14 +- 2.14   4.00 +- 0.88
L1_SingleMu5 0.60 +- 0.35   0.42 +- 0.30   0.60 +- 0.35   2.00 +- 0.63   31.72 +- 2.13   2.00 +- 0.63
L1_SingleMu7 0.40 +- 0.28   0.21 +- 0.21   0.00 +- 0.00   1.20 +- 0.49   31.51 +- 2.13   1.00 +- 0.44
L1_SingleMuOpen 3.20 +- 0.79   1.05 +- 0.47   1.40 +- 0.53   10.40 +- 1.37   34.24 +- 2.17   31.80 +- 2.08

  • Results of Without DT, when using the sample with Quarkonia and Weak Bosons on at 362
CMSSW 362  
L1_DoubleMu3 2.80 +- 0.74
L1_DoubleMuOpen 3.40 +- 0.81
L1_SingleMu0 4.40 +- 0.92
L1_SingleMu10 0.40 +- 0.28
L1_SingleMu14 0.40 +- 0.28
L1_SingleMu20 0.20 +- 0.20
L1_SingleMu3 3.60 +- 0.83
L1_SingleMu5 1.40 +- 0.53
L1_SingleMu7 0.80 +- 0.40
L1_SingleMuOpen 10.20 +- 1.35

In p+p simulations

  • rates in p+p
    • CMSSW_3_6_2/RelValMinBias/GEN-SIM-DIGI-RAW-HLTDEBUG/MC_36Y_V10-v1
    • CMSSW_3_7_0/RelValMinBias/GEN-SIM-DIGI-RAW-HLTDEBUG/MC_37Y_V4-v1
    • CMSSW_3_8_1/RelValMinBias/GEN-SIM-DIGI-RAW-HLTDEBUG/MC_38Y_V8-v1

p+p 362     370     381
L1_DoubleMu3 0.01 +- 0.01   0.03 +- 0.02   0.02 +- 0.02
L1_DoubleMuOpen 0.01 +- 0.01   0.03 +- 0.02   0.02 +- 0.02
L1_SingleMu0 0.06 +- 0.02   0.69 +- 0.09   0.10 +- 0.03
L1_SingleMu10 0.00 +- 0.00   0.66 +- 0.09   0.00 +- 0.00
L1_SingleMu14 0.00 +- 0.00   0.64 +- 0.08   0.00 +- 0.00
L1_SingleMu20 0.00 +- 0.00   0.63 +- 0.08   0.00 +- 0.00
L1_SingleMu3 0.03 +- 0.02   0.69 +- 0.09   0.07 +- 0.03
L1_SingleMu5 0.01 +- 0.01   0.69 +- 0.09   0.03 +- 0.02
L1_SingleMu7 0.00 +- 0.00   0.66 +- 0.09   0.00 +- 0.00
L1_SingleMuOpen 0.09 +- 0.03   0.73 +- 0.09   0.91 +- 0.10

L1 muon information

  • CMSSW 381
    • among 304 L1 muons, 243 L1muons have 140 GeV/c in the high eta region. This is quite weird. Removing 140 GeV/c L1muons, the number of L1 muos is only 60 muons.
    • We found those 140 GeV/c L1 muons in the high eta region. Is this misassigned pt information on the L1 emulation ?
pt of L1 muons eta of L1 muons
pt vs eta of L1 muons pt vs eta of L1 muons, excluding 140 GeV muons at L1
    • gen and sim muons info
pt vs eta of gen muons pt vs eta of sim muons pt vs eta of sim muons > 1 GeV/c

  • CMSSW 370

  • CMSSW 362

Comparison

  • Detail comparison given by Camelia

Details of the investigation

  • Compare the effect of CSCTFinder with 3 kinds of version : 362, 370, 381
    • Running cfg files : cfg
  • Test options
    • (1) : Remove CSC trigger in the L1 emulation at 362(cfg), 370(cfg), 381(cfg).
      process.simGmtDigis.CSCCandidates = cms.InputTag('none')
      process.simDttfDigis.CSCStub_Source = cms.InputTag('none')
      process.SimL1MuTriggerPrimitives.remove(process.simCscTriggerPrimitiveDigis)
      process.SimL1MuTrackFinders.remove(process.simCsctfTrackDigis)
      process.SimL1MuTrackFinders.remove(process.simCsctfDigis)
    • (2) : Remove DT trigger in the L1 emulation at 362(cfg), 370(cfg), 381(cfg).
      process.simGmtDigis.DTCandidates = cms.InputTag('none')
      process.simCsctfTrackDigis.DTproducer = cms.untracked.InputTag('none')
      process.simCsctfTrackDigis.useDT = cms.bool(False)
      process.simDtTriggerPrimitiveDigis.digiTag =  cms.InputTag('none')
      process.simDttfDigis.DTDigi_Source = cms.InputTag('none')
      process.SimL1MuTriggerPrimitives.remove(process.simDtTriggerPrimitiveDigis)
      process.SimL1MuTrackFinders.remove(process.simDttfDigis)
    • (3) : Test "none" option of CSC and DT at 362, 370, 381.
      process.simGmtDigis.CSCCandidates = cms.InputTag('none')
      process.simGmtDigis.DTCandidates = cms.InputTag('none')
    • (4) : Remove RPC trigger in the L1 emulation at 362, 370, 381.

Study logs

  • Dongho is investigating the fake rate roughly and the dependence of several kinematics (eta, pt)
  • Catherine is testing the effect of number of muons in the gen level, because from the 370 the hydjet sample has Quarkonia and Weak Bosons.
  • According to study so far(2010/09/22), it has tendency for the case of "Without CSC" to reduce more than the case of "Without DT", which means the huge rate of L1 should be affected by CSC trigger.
  • Testing the option (1), (2) in the local
    • (1) : 362, 370, 381 are fine and finished
    • (2) : 362, 370 look fine but 381 has error cfg
      getByLabel: Found zero products matching all criteria
      Looking for type: L1MuDTChambPhContainer
      Looking for module label: simDtTriggerPrimitiveDigis
      Looking for productInstanceName: 

      cms::Exception going through module CSCTFTrackProducer/simCsctfTrackDigis run: 1 lumi: 666672 event: 31
      If you wish to continue processing events after a ProductNotFound exception,
  • After back to normal process (not removing anythings), process.simGmtDigis.CSCCandidates = cms.InputTag('none') test is ongoing in the Batch(362, 370)
  • option input contents to drop L1 Contents
         process.source.inputCommands = cms.untracked.vstring(
                      'keep *', 
                      'drop *_simCscTriggerPrimitiveDigis_*_*', 
                      'drop *_simDtTriggerPrimitiveDigis_*_*', 
                      'drop *_simRpcTriggerDigis_*_*', 
                      'drop *_simRctDigis_*_*', 
                      'drop *_simCsctfDigis_*_*', 
                      'drop *_simCsctfTrackDigis_*_*', 
                      'drop *_simDttfDigis_*_*', 
                      'drop *_simGctDigis_*_*', 
                      'drop *_simGmtDigis_*_*', 
                      'drop *_simGtDigis_*_*', 
                      'drop *_cscTriggerPrimitiveDigis_*_*', 
                      'drop *_dtTriggerPrimitiveDigis_*_*', 
                      'drop *_rpcTriggerDigis_*_*', 
                      'drop *_rctDigis_*_*', 
                      'drop *_csctfDigis_*_*', 
                      'drop *_csctfTrackDigis_*_*', 
                      'drop *_dttfDigis_*_*', 
                      'drop *_gctDigis_*_*', 
                      'drop *_gmtDigis_*_*', 
                      'drop *_gtDigis_*_*',
                      'drop *_gtEvmDigis_*_*',
                      'drop *_l1GtRecord_*_*', 
                      'drop *_l1GtTriggerMenuLite_*_*',
                      'drop *_conditionsInEdm_*_*',
                      'drop *_l1extraParticles_*_*',
                      'drop LumiDetails_lumiProducer_*_*',
                      'drop LumiSummary_lumiProducer_*_*'
          )

  • In the case of 381, when I run L1 and L1Analyzer, it works well, although it doesn't have CSC part. But the result is exactly same as what it has CSC part, which means the CSC doesn't seem to affect on the trigger rate.
  • the used cfg file here
    • the results
                      L1_DoubleMu3      15.80 +-       1.63
                   L1_DoubleMuOpen      16.40 +-       1.66
                      L1_SingleMu0       5.60 +-       1.03
                     L1_SingleMu10       0.40 +-       0.28
                     L1_SingleMu14       0.20 +-       0.20
                     L1_SingleMu20       0.20 +-       0.20
                      L1_SingleMu3       4.20 +-       0.90
                      L1_SingleMu5       1.80 +-       0.59
                      L1_SingleMu7       1.00 +-       0.44
                   L1_SingleMuOpen      31.80 +-       2.08

  • the test of "NONE" option of CSC Candidates applied. cfg
        process.simGmtDigis.CSCCandidates = cms.InputTag('none')
    • the results :
                      L1_DoubleMu3       0.00 +-       0.00
                   L1_DoubleMuOpen       0.00 +-       0.00
                      L1_SingleMu0       1.00 +-       0.44
                     L1_SingleMu10       0.00 +-       0.00
                     L1_SingleMu14       0.00 +-       0.00
                     L1_SingleMu20       0.00 +-       0.00
                      L1_SingleMu3       1.00 +-       0.44
                      L1_SingleMu5       0.60 +-       0.35
                      L1_SingleMu7       0.00 +-       0.00
                   L1_SingleMuOpen       1.40 +-       0.53

  • In the case of 370 (solved)
    • When running from L1 without CSC part, the error of RPCDigis happended like followings cfg
         %MSG-e RPCTechnicalTrigger:  RPCTechnicalTrigger:simRpcTechTrigDigis
         22-Sep-2010 01:15:42 CEST Run: 1 Event: 1
         can't find RPCDigiCollection with label: InputTag:  label =
         simMuonRPCDigis, instance =
    • When running from pdigi to avoid RPCDigis error, the other error came up like following cfg
         22-Sep-2010 01:10:15 CEST  Closed file
         rfio:/castor/cern.ch/cms/store/relval/CMSSW_3_7_0/RelValHydjetQ_MinBias_2760GeV/GEN-SIM-RAW/MC_37Y_V4-v1/0025
         /00C18D3F-6869-DF11-8687-0026189438AB.root
         %MSG-s CMSException:  AfterFile 22-Sep-2010 01:10:15 CEST PostEndRun
         cms::Exception caught in cmsRun
         ---- EventProcessorFailure BEGIN
         EventProcessingStopped
         ---- ScheduleExecutionFailure BEGIN
         ProcessingStopped
         ---- ProductNotFound BEGIN
         getByLabel: Found zero products matching all criteria
         Looking for type: CSCTriggerContainer
         Looking for module label: simCsctfTrackDigis
         Looking for productInstanceName:

         cms::Exception going through module DTTrackFinder/simDttfDigis run: 1
         lumi: 666705 event: 200
         If you wish to continue processing events after a ProductNotFound exception,
         add "SkipEvent = cms.untracked.vstring('ProductNotFound')" to the
         "options" PSet in the configuration.
         ---- ProductNotFound END
         Exception going through path L1simulation_step
         ---- ScheduleExecutionFailure END
         an exception occurred during current event processing
         cms::Exception caught in EventProcessor and rethrown
         ---- EventProcessorFailure END

-- DongHoMoon - 21 Sep 2010

Topic attachments
I Attachment History Action Size Date Who Comment
Unknown file formatcc HLTrgAna.cc r2 r1 manage 27.8 K 2010-09-29 - 23:11 UnknownUser  
Texttxt csctfTrackDigis_cfi.py.txt r1 manage 8.6 K 2010-09-29 - 15:44 UnknownUser  
Texttxt l1_Ana_381_final.py.txt r1 manage 7.9 K 2010-09-29 - 15:45 UnknownUser  
Edit | Attach | Watch | Print version | History: r22 < r21 < r20 < r19 < r18 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r22 - 2010-11-09 - DongHoMoon
 
KOREA-CMS
This site is powered by the TWiki collaboration platform Powered by Perl This site is powered by the TWiki collaboration platformCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding KoreaCmsWiki? Send feedback