Purpose of the test

Motivation : See Camelia's talk - http://indico.cern.ch/getFile.py/access?contribId=5&resId=0&materialId=slides&confId=110336

Which subset of these are useful enough (and sufficiently smaller than full RECO content of triggered + prescaled MB sample) to produce in conjunction with PromptReco and re-reco?

● We do not have the numbers to make a decision

● Can think at other scenarios (for trigger/cuts combination)

● IsSTA, isGLB; should be extended to isTRK (from pp experience) – check how much increases data volume by doing this

Skimming option

  • Legends
    • cs = central skim
    • sd = secondary dataset
    • pd = primary dataset

  • SD (Secondary dataset)
process.filterSdMuHI = cms.Path(process.hltMuHI)
# jet trigger
hltJetHI = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
hltJetHI.HLTPaths = ["HLT_HIJet35U"]
hltJetHI.throw = False
hltJetHI.andOr = True
# dimuon trigger
hltMuHI = hltJetHI.clone()
#hltMuHI.HLTPaths = ["HLT_HIL1DoubleMuOpen"]
hltMuHI.HLTPaths = ["HLT_L1DoubleMuOpen"]

  • CS (Central Skim)
    • Single muon
muonSelector = cms.EDFilter("MuonSelector",
    src = cms.InputTag("muons"),
    cut = cms.string("(isStandAloneMuon || isGlobalMuon) && pt > 1."),
    filter = cms.bool(True)
    )
muonFilter = cms.EDFilter("MuonCountFilter",
    src = cms.InputTag("muonSelector"),
    minNumber = cms.uint32(1)
    )
muonTrkGlbSelector = cms.EDFilter("MuonSelector",
                                 src = cms.InputTag("muons"),
                                 cut = cms.string("(isTrackerMuon || isStandAloneMuon || isGlobalMuon) && pt > 1."),
                                 filter = cms.bool(True)
                                 )

muonTrkGlbFilter = cms.EDFilter("MuonCountFilter",
                          src = cms.InputTag("muonTrkGlbSelector"),
                          minNumber = cms.uint32(1)
                          )


process.filterCsMuHInoTRK = cms.Path(process.muonSelector * process.muonFilter * process.makeHeavyIonMuons)
process.filterCsMuHIyesTRK = cms.Path(process.muonTrkGlbSelector * process.muonTrkGlbFilter * process.makeHeavyIonMuons)
    • Dimuon
      • Default cut
dimuonsMassCut = cms.EDProducer("CandViewShallowCloneCombiner",
    checkCharge = cms.bool(True),
    cut = cms.string('mass > 70 & mass < 120 & charge=0'),
    decay = cms.string("muonSelector@+ muonSelector@-")
    )

dimuonsMassCutFilter = cms.EDFilter("CandViewCountFilter",
    src = cms.InputTag("dimuonsMassCut"),
    minNumber = cms.uint32(1)
    )
dimuonsCut = cms.EDProducer("CandViewShallowCloneCombiner",
                            checkCharge = cms.bool(True),
                            cut = cms.string('charge=0'),
                            decay = cms.string("muonSelector@+ muonSelector@-")
                            )

dimuonsTrkGlbCut = cms.EDProducer("CandViewShallowCloneCombiner",
                            checkCharge = cms.bool(True),
                            cut = cms.string('charge=0'),
                            decay = cms.string("muonTrkGlbSelector@+ muonTrkGlbSelector@-")
                            )

dimuonsJPsiMassCut = cms.EDProducer("CandViewShallowCloneCombiner",
                                    checkCharge = cms.bool(True),
                                    cut = cms.string('mass > 2.1  & mass < 4.1 & charge=0'),
                                    decay = cms.string("muonTrkGlbSelector@+ muonTrkGlbSelector@-")
                                    )

dimuonsnoTRKJPsiMassCut = cms.EDProducer("CandViewShallowCloneCombiner",
                                    checkCharge = cms.bool(True),
                                    cut = cms.string('mass > 2.1  & mass < 4.1 & charge=0'),
                                    decay = cms.string("muonSelector@+ muonSelector@-"
)
                                    )
dimuonsUpsilonMassCut = cms.EDProducer("CandViewShallowCloneCombiner",
                                       checkCharge = cms.bool(True),
                                       cut = cms.string('mass > 8  & mass < 12  & charge=0'),
                                       decay = cms.string("muonTrkGlbSelector@+ muonTrkGlbSelector@-")
                                       )
dimuonsZ0MassCut = cms.EDProducer("CandViewShallowCloneCombiner",
                                  checkCharge = cms.bool(True),
                                  cut = cms.string('mass > 70  & mass < 110  & charge=0'),
                                  decay = cms.string("muonSelector@+ muonSelector@-")
                                  )

dimuonsTrkGlbZ0MassCut = cms.EDProducer("CandViewShallowCloneCombiner",
                                  checkCharge = cms.bool(True),
                                  cut = cms.string('mass > 70  & mass < 110  & charge=0'),
                                  decay = cms.string("muonTrkGlbSelector@+ muonTrkGlbSelector@-")
                                  )


dimuonsnoTRKUpsilonMassCut = cms.EDProducer("CandViewShallowCloneCombiner",
                                       checkCharge = cms.bool(True),
                                       cut = cms.string('mass > 8  & mass < 12  & charge=0'),
                                       decay = cms.string("muonSelector@+ muonSelector@-")
                                       )



process.filterCsZmumuHInoTRK = cms.Path(process.muonSelector * process.muonFilter
 * process.dimuonsMassCut * process.dimuonsMassCutFilter * process.makeHeavyIonMuons)
process.filterCsZmumuHIyesTRK = cms.Path(process.muonTrkGlbSelector * process.muonTrkGlbFilter
 * process.dimuonsMassCut * process.dimuonsMassCutFilter * process.makeHeavyIonMuons)
process.filterCsZmumuHInoTRKpassHLT = cms.Path(process.hltMuHI * process.muonSelector * process.muonFilter
 * process.dimuonsMassCut * process.dimuonsMassCutFilter * process.makeHeavyIonMuons)
process.filterCsZmumuHIyesTRKpassHLT = cms.Path(process.hltMuHI * process.muonTrkGlbSelector * process.muonTrkGlbFilter
 * process.dimuonsMassCut * process.dimuonsMassCutFilter * process.makeHeavyIonMuons)
      • No cut
process.filterCsdimuonHInoTRK = cms.Path(process.muonSelector * process.muonFilter
 * process.dimuonsCut * process.dimuonsCutFilter * process.makeHeavyIonMuons)
process.filterCsdimuonHIyesTRK = cms.Path(process.muonTrkGlbSelector * process.muonTrkGlbFilter
 * process.dimuonsCut * process.dimuonsCutFilter * process.makeHeavyIonMuons)
process.filterCsdimuonHInoTRKpassHLT = cms.Path(process.hltMuHI * process.muonSelector * process.muonFilter
 * process.dimuonsCut * process.dimuonsCutFilter * process.makeHeavyIonMuons)
process.filterCsdimuonHIyesTRKpassHLT = cms.Path(process.hltMuHI * process.muonTrkGlbSelector * process.muonTrkGlbFilter
 * process.dimuonsCut * process.dimuonsCutFilter * process.makeHeavyIonMuons)
      • new filter for JPsi
process.filterCsdimuonJPsiHInoTRK = cms.Path(process.muonSelector * process.muonFilter
 * process.dimuonsJPsiMassCut * process.dimuonsJPsiMassCutFilter * process.makeHeavyIonMuons)
process.filterCsdimuonJPsiHIyesTRK = cms.Path(process.muonTrkGlbSelector * process.muonTrkGlbFilter
 * process.dimuonsJPsiMassCut * process.dimuonsJPsiMassCutFilter * process.makeHeavyIonMuons)
process.filterCsdimuonJPsiHInoTRKpassHLT = cms.Path(process.hltMuHI * process.muonSelector * process.muonFilter
 * process.dimuonsJPsiMassCut * process.dimuonsJPsiMassCutFilter * process.makeHeavyIonMuons)
process.filterCsdimuonJPsiHIyesTRKpassHLT = cms.Path(process.hltMuHI * process.muonTrkGlbSelector * process.muonTrkGlbFilter
 * process.dimuonsJPsiMassCut * process.dimuonsJPsiMassCutFilter * process.makeHeavyIonMuons)
      • new filter for Upsilon
process.filterCsdimuonUpsilonHInoTRK = cms.Path(process.muonSelector * process.muonFilter
 * process.dimuonsUpsilonMassCut * process.dimuonsUpsilonMassCutFilter * process.makeHeavyIonMuons)
process.filterCsdimuonUpsilonHIyesTRK = cms.Path(process.muonTrkGlbSelector * process.muonTrkGlbFilter
 * process.dimuonsUpsilonMassCut * process.dimuonsUpsilonMassCutFilter * process.makeHeavyIonMuons)
process.filterCsdimuonUpsilonHInoTRKpassHLT = cms.Path(process.hltMuHI * process.muonSelector * process.muonFilter
 * process.dimuonsUpsilonMassCut * process.dimuonsUpsilonMassCutFilter * process.makeHeavyIonMuons)
process.filterCsdimuonUpsilonHIyesTRKpassHLT = cms.Path(process.hltMuHI * process.muonTrkGlbSelector * process.muonTrkGlbFilter
 * process.dimuonsUpsilonMassCut * process.dimuonsUpsilonMassCutFilter * process.makeHeavyIonMuons)
      • new filter for Z0
process.filterCsdimuonZ0HInoTRK = cms.Path(process.muonSelector * process.muonFilter
 * process.dimuonsZ0MassCut * process.dimuonsZ0MassCutFilter * process.makeHeavyIonMuons)
process.filterCsdimuonZ0HIyesTRK = cms.Path(process.muonTrkGlbSelector * process.muonTrkGlbFilter
 * process.dimuonsZ0MassCut * process.dimuonsZ0MassCutFilter * process.makeHeavyIonMuons)
process.filterCsdimuonZ0HInoTRKpassHLT = cms.Path(process.hltMuHI * process.muonSelector * process.muonFilter
 * process.dimuonsZ0MassCut * process.dimuonsZ0MassCutFilter * process.makeHeavyIonMuons)
process.filterCsdimuonZ0HIyesTRKpassHLT = cms.Path(process.hltMuHI * process.muonTrkGlbSelector * process.muonTrkGlbFilter
 * process.dimuonsZ0MassCut * process.dimuonsZ0MassCutFilter * process.makeHeavyIonMuons)

Procedure

Code & Scripts

  • cfg file for skimming
  • modified analysisFilter_cff.py (original location :
  • modified HLT (original location :
  • crab.cfg

Sample

  • Wei's sample (dbs link)
    • MinBias : /Hydjet_Quenched_MinBias_2760GeV_START39V0_GEN_SIM_RAW_390pre5/davidlw-Hydjet_Quenched_MinBias_2760GeV_START39V0_GEN_SIM_RECO_390pre5-6dc75065b07f0b7687f52fc678e8a8c4/USER
    • SingleJPsimumu : /Hydjet_Quenched_MinBias_2760GeV_START39V0_GEN_SIM_RAW_390pre5/davidlw-SingleJPsimumu_FlatPt0to20_START39V1_GEN_SIM_RECO_390pre5-6dc75065b07f0b7687f52fc678e8a8c4/USER
    • SingleUpsimumu : /Hydjet_Quenched_MinBias_2760GeV_START39V0_GEN_SIM_RAW_390pre5/davidlw-SingleUpsilonmumu_FlatPt0to20_START39V1_GEN_SIM_RECO_390pre5-6dc75065b07f0b7687f52fc678e8a8c4/USER
    • SingleZ0mumu : /Hydjet_Quenched_MinBias_2760GeV_START39V0_GEN_SIM_RAW_390pre5/davidlw-SingleZmumu_FlatPt0to20_START39V1_GEN_SIM_RECO_390pre5-6dc75065b07f0b7687f52fc678e8a8c4/USER
  • JEX (dbs link)
    • MinBias : /MinimumBiasHI/kimy-Spring10-JulyAnalysisExercise_MC_37Y_V4-NZS-RECO-prelim-v1-599a3a4efd1a163b86e329cbf44e9f6a/USER

Sotrage of the skimmed file

  • Wei's sample
    • MinBias : /castor/cern.ch/user/h/hckim/TestJEX/crab/v3/MinBias
    • SingleJPsimumu : /castor/cern.ch/user/h/hckim/TestJEX/crab/v3/SingleJPsimumu
    • SingleUpsimumu : /castor/cern.ch/user/h/hckim/TestJEX/crab/v3/SingleUpsimumu
    • SingleZ0mumu : /castor/cern.ch/user/h/hckim/TestJEX/crab/v3/SingleZ0mumu
  • JEX
    • MinBias : /castor/cern.ch/user/h/hckim/TestJEX/crab/v3/JEXMinBias

Result

Number of Event

Skimming option MinBias SingleJPsimumu SingleUpsimumu SinlgleZ0mumu
Total Event 5950 5320 5280 5450
Basic Filter pass HLT Incl. Trk No. of Event skimming Efficiency No. of Event skimming Efficiency No. of Event skimming Efficiency No. of Event skimming Efficiency
SdMuHI Yes ? 0 0 2031 0.382 2662 0.504 3836 0.704
CsMuHI No No 778 0.131 4546 0.855 5021 0.951 5427 0.996
  Yes 2258 0.379 4879 0.917 5113 0.968 5436 0.997
CsZmumuHI No No 0 0 2 0.000 2 0.000 3217 0.590
  Yes 0 0 2 0.000 2 0.000 3217 0.590
Yes No 0 0 1 0.000 2 0.000 3130 0.574
  Yes 0 0 1 0.000 2 0.000 3130 0.574
CsdimuonHI No No 109 0.018 1837 0.345 2578 0.488 3662 0.672
  Yes 109 0.018 1837 0.345 2578 0.488 3662 0.672
Yes No 0 0 1669 0.314 2359 0.447 3540 0.650
  Yes 0 0 1669 0.314 2359 0.447 3540 0.650
CsdimuonJPsi No No 589 0.099 2736 0.514 1646 0.312 1554 0.285
  Yes 1468 0.247 2894 0.544 1685 0.319 1557 0.286
Yes No 0 0 1851 0.348 1035 0.196 1213 0.223
  Yes 0 0 1851 0.348 1035 0.196 1213 0.223
CsdimuonUpsi No No 466 0.078 1223 0.230 3185 0.603 1439 0.264
  Yes 1049 0.176 1318 0.248 3213 0.609 1442 0.265
Yes No 0 0 727 0.137 2446 0.463 1123 0.206
  Yes 0 0 727 0.137 2446 0.463 1123 0.206
CsdimuonZ0 No No 0 0 2 0.000 2 0.000 3186 0.585
  Yes 0 0 2 0.000 2 0.000 3186 0.585
Yes No 0 0 1 0.000 2 0.000 3101 0.569
  Yes 0 0 1 0.000 2 0.000 3101 0.569

Data size

Skimming option MinBias SingleJPsimumu SingleUpsimumu SinlgleZ0mumu
Basic Filter pass HLT Incl. Trk No. of Event skimming Efficiency No. of Event skimming Efficiency No. of Event skimming Efficiency No. of Event skimming Efficiency
SdMuHI Yes ?
CsMuHI No No
  Yes
CsZmumuHI No No
  Yes
Yes No
  Yes
CsdimuonHI No No
  Yes
Yes No
  Yes
CsdimuonJPsi No No
  Yes
Yes No
  Yes
CsdimuonUpsi No No
  Yes
Yes No
  Yes
CsdimuonZ0 No No
  Yes
Yes No
  Yes

-- HyunChulKim - 22 Oct 2010

Edit | Attach | Watch | Print version | History: r4 < r3 < r2 < r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r4 - 2010-10-26 - HyunChulKim
 
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