1 |
scasasso |
1.1 |
#################
|
2 |
|
|
#PATMuons analysis
|
3 |
|
|
##################
|
4 |
|
|
|
5 |
|
|
import FWCore.ParameterSet.Config as cms
|
6 |
|
|
|
7 |
|
|
process = cms.Process("TEST")
|
8 |
|
|
# Messages
|
9 |
|
|
process.load("FWCore.MessageService.MessageLogger_cfi")
|
10 |
|
|
process.MessageLogger.destinations = ['cout', 'cerr']
|
11 |
|
|
process.MessageLogger.cerr.FwkReport.reportEvery = 1000
|
12 |
|
|
|
13 |
|
|
process.load("Configuration.StandardSequences.Geometry_cff")
|
14 |
|
|
process.load("Configuration.StandardSequences.MagneticField_38T_cff")
|
15 |
|
|
process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
|
16 |
|
|
|
17 |
|
|
process.GlobalTag.globaltag = 'GR_R_39X_V5::All'
|
18 |
|
|
|
19 |
|
|
|
20 |
|
|
process.source = cms.Source("PoolSource",
|
21 |
|
|
fileNames = cms.untracked.vstring(
|
22 |
|
|
#'rfio:/castor/cern.ch/cms/store/user/emiglior/test/zzSkim/test/PAT_03Mar11/Mu_ZbbSkimDec22ReReco_PAT397_03Mar11.root'
|
23 |
|
|
'rfio:/castor/cern.ch/cms/store/user/emiglior/zbb/ZbbSkimDec22ReReco_PAT397_17Mar11/Mu/MergedOutputFile_1_2_Ala.root'
|
24 |
|
|
)
|
25 |
|
|
)
|
26 |
|
|
|
27 |
|
|
|
28 |
|
|
|
29 |
|
|
##################
|
30 |
|
|
|
31 |
|
|
# TRIGGER FILTER #
|
32 |
|
|
##################
|
33 |
|
|
|
34 |
|
|
|
35 |
|
|
### Trigger path FIRED in the event
|
36 |
|
|
|
37 |
|
|
import HLTrigger.HLTfilters.hltHighLevel_cfi
|
38 |
|
|
process.hltFilter = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
|
39 |
|
|
process.hltFilter.HLTPaths = [ "HLT_DoubleMu3_v*" ]
|
40 |
|
|
process.hltFilter.TriggerResultsTag = cms.InputTag("TriggerResults","","HLT")
|
41 |
|
|
process.hltFilter.throw = cms.bool(False)
|
42 |
|
|
|
43 |
|
|
|
44 |
|
|
process.maxEvents = cms.untracked.PSet(
|
45 |
|
|
input = cms.untracked.int32(-1)
|
46 |
|
|
)
|
47 |
|
|
|
48 |
|
|
process.TFileService=cms.Service('TFileService',
|
49 |
|
|
fileName=cms.string('Mu_2010_Plots_FromPAT_03Mar11.root')
|
50 |
|
|
)
|
51 |
|
|
###############################################################
|
52 |
|
|
#First Analyzer: general spectra for all PATMuons + filter studies
|
53 |
|
|
################################################################
|
54 |
|
|
process.MuAllSpectra = cms.EDAnalyzer("MuAnalyzer",
|
55 |
|
|
MuonCollection = cms.untracked.string('patMuonsWithTrigger'),
|
56 |
|
|
JetCollectionPF = cms.untracked.string('cleanPatJetsPF'),
|
57 |
|
|
ZllCollection = cms.untracked.string('zMMCand'),
|
58 |
|
|
numEventsNames = cms.untracked.vstring('TotalEventCounter','AfterPVFilterCounter', 'AfterNSFilterCounter', 'AfterPATCounter', 'AfterCandidatesCounter', 'AfterJetsCounter')
|
59 |
|
|
)
|
60 |
|
|
|
61 |
|
|
|
62 |
|
|
####################################################
|
63 |
|
|
#Second Analyzer: spectra about preselection analysis
|
64 |
|
|
####################################################
|
65 |
|
|
process.MuPreselSpectra = cms.EDAnalyzer("PreselAnalyzer",
|
66 |
|
|
MuonCollection = cms.untracked.string('patMuonsWithTrigger'),
|
67 |
|
|
ZllCollection = cms.untracked.string('zMMCand'),
|
68 |
|
|
)
|
69 |
|
|
|
70 |
|
|
###############################################################
|
71 |
|
|
#Third Analyzer: Jet collections studies
|
72 |
|
|
################################################################
|
73 |
|
|
process.JetLepton = cms.EDAnalyzer("JetLeptonCleanAnalyzer",
|
74 |
|
|
MuonCollection = cms.untracked.string('patMuonsWithTrigger'),
|
75 |
|
|
ZllCollection = cms.untracked.string('zMMCand'),
|
76 |
|
|
JetCollection = cms.untracked.string('cleanPatJets'),
|
77 |
|
|
JetCollectionPU = cms.untracked.string('cleanPatJetsNoPU'),
|
78 |
|
|
JetCollectionJPT = cms.untracked.string('cleanPatJetsJPT'),
|
79 |
|
|
JetCollectionPF = cms.untracked.string('cleanPatJetsPF'),
|
80 |
|
|
)
|
81 |
|
|
|
82 |
|
|
|
83 |
|
|
|
84 |
|
|
process.p = cms.Path(process.hltFilter *
|
85 |
|
|
( process.MuAllSpectra +
|
86 |
|
|
process.MuPreselSpectra +
|
87 |
|
|
process.JetLepton) )
|
88 |
|
|
|
89 |
|
|
|
90 |
|
|
|
91 |
|
|
|