ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/CmsHi/JetAnalysis/python/HiTrackJets_cff.py
Revision: 1.1
Committed: Wed Feb 13 18:19:26 2013 UTC (12 years, 2 months ago) by yilmaz
Content type: text/x-python
Branch: MAIN
Log Message:
TracksJets and rho calculation

File Contents

# User Rev Content
1 yilmaz 1.1 import FWCore.ParameterSet.Config as cms
2    
3     from RecoJets.JetProducers.ak5TrackJets_cfi import *
4    
5     from PhysicsTools.RecoAlgos.TrackWithVertexRefSelector_cfi import *
6     from RecoJets.JetProducers.TracksForJets_cff import *
7    
8     trackWithVertexRefSelector.src = cms.InputTag("hiTracks")
9     trackWithVertexRefSelector.ptErrorCut = cms.double(99999)
10     trackWithVertexRefSelector.d0Max = cms.double(99999)
11     trackWithVertexRefSelector.numberOfValidHits = cms.uint32(999999)
12     trackWithVertexRefSelector.useVtx = cms.bool(False)
13    
14     ak5TrackJets.doAreaFastjet = cms.bool(True)
15     ak5TrackJets.doRhoFastjet = cms.bool(True)
16     ak5TrackJets.doFastJetNonUniform = cms.bool(True)
17     ak5TrackJets.puCenters = cms.vdouble(-5,-4.5,-4,-3.5,-3,-2.5,-2,-1.5,-1,-0.5,0,0.5,1,1.5,2,2.5,3,3.5,4,4.5,5)
18     ak5TrackJets.puWidth = cms.double(0.5)
19     ak5TrackJets.nExclude = cms.uint32(0)
20     ak5TrackJets.Rho_EtaMax = cms.double(5.5)
21     ak5TrackJets.Ghost_EtaMax = cms.double(5.5)
22    
23     kt5TrackJets = ak5TrackJets.clone(jetAlgorithm = cms.string("Kt"))
24    
25     ak1TrackJets05 = ak5TrackJets.clone(rParam = cms.double(0.1),puWidth = 0.5)
26     ak2TrackJets05 = ak5TrackJets.clone(rParam = cms.double(0.2),puWidth = 0.5)
27     ak3TrackJets05 = ak5TrackJets.clone(rParam = cms.double(0.3),puWidth = 0.5)
28     ak4TrackJets05 = ak5TrackJets.clone(rParam = cms.double(0.4),puWidth = 0.5)
29     ak5TrackJets05 = ak5TrackJets.clone(rParam = cms.double(0.5),puWidth = 0.5)
30     ak6TrackJets05 = ak5TrackJets.clone(rParam = cms.double(0.6),puWidth = 0.5)
31    
32     ak1TrackJets1 = ak5TrackJets.clone(rParam = cms.double(0.1),puWidth = 1)
33     ak2TrackJets1 = ak5TrackJets.clone(rParam = cms.double(0.2),puWidth = 1)
34     ak3TrackJets1 = ak5TrackJets.clone(rParam = cms.double(0.3),puWidth = 1)
35     ak4TrackJets1 = ak5TrackJets.clone(rParam = cms.double(0.4),puWidth = 1)
36     ak5TrackJets1 = ak5TrackJets.clone(rParam = cms.double(0.5),puWidth = 1)
37     ak6TrackJets1 = ak5TrackJets.clone(rParam = cms.double(0.6),puWidth = 1)
38    
39     ak1TrackJets2 = ak5TrackJets.clone(rParam = cms.double(0.1),puWidth = 2)
40     ak2TrackJets2 = ak5TrackJets.clone(rParam = cms.double(0.2),puWidth = 2)
41     ak3TrackJets2 = ak5TrackJets.clone(rParam = cms.double(0.3),puWidth = 2)
42     ak4TrackJets2 = ak5TrackJets.clone(rParam = cms.double(0.4),puWidth = 2)
43     ak5TrackJets2 = ak5TrackJets.clone(rParam = cms.double(0.5),puWidth = 2)
44     ak6TrackJets2 = ak5TrackJets.clone(rParam = cms.double(0.6),puWidth = 2)
45    
46     kt1TrackJets05 = kt5TrackJets.clone(rParam = cms.double(0.1),puWidth = 0.5)
47     kt2TrackJets05 = kt5TrackJets.clone(rParam = cms.double(0.2),puWidth = 0.5)
48     kt3TrackJets05 = kt5TrackJets.clone(rParam = cms.double(0.3),puWidth = 0.5)
49     kt4TrackJets05 = kt5TrackJets.clone(rParam = cms.double(0.4),puWidth = 0.5)
50     kt5TrackJets05 = kt5TrackJets.clone(rParam = cms.double(0.5),puWidth = 0.5)
51     kt6TrackJets05 = kt5TrackJets.clone(rParam = cms.double(0.6),puWidth = 0.5)
52    
53     kt1TrackJets1 = kt5TrackJets.clone(rParam = cms.double(0.1),puWidth = 1)
54     kt2TrackJets1 = kt5TrackJets.clone(rParam = cms.double(0.2),puWidth = 1)
55     kt3TrackJets1 = kt5TrackJets.clone(rParam = cms.double(0.3),puWidth = 1)
56     kt4TrackJets1 = kt5TrackJets.clone(rParam = cms.double(0.4),puWidth = 1)
57     kt5TrackJets1 = kt5TrackJets.clone(rParam = cms.double(0.5),puWidth = 1)
58     kt6TrackJets1 = kt5TrackJets.clone(rParam = cms.double(0.6),puWidth = 1)
59    
60     kt1TrackJets2 = kt5TrackJets.clone(rParam = cms.double(0.1),puWidth = 2)
61     kt2TrackJets2 = kt5TrackJets.clone(rParam = cms.double(0.2),puWidth = 2)
62     kt3TrackJets2 = kt5TrackJets.clone(rParam = cms.double(0.3),puWidth = 2)
63     kt4TrackJets2 = kt5TrackJets.clone(rParam = cms.double(0.4),puWidth = 2)
64     kt5TrackJets2 = kt5TrackJets.clone(rParam = cms.double(0.5),puWidth = 2)
65     kt6TrackJets2 = kt5TrackJets.clone(rParam = cms.double(0.6),puWidth = 2)
66    
67    
68     recoTrackJets = cms.Sequence(trackWithVertexRefSelector*trackRefsForJets
69     *
70     ak1TrackJets05+ak2TrackJets05+ak3TrackJets05+ak4TrackJets05+ak5TrackJets05+ak6TrackJets05
71     +
72     ak1TrackJets1+ak2TrackJets1+ak3TrackJets1+ak4TrackJets1+ak5TrackJets1+ak6TrackJets1
73     +
74     ak1TrackJets2+ak2TrackJets2+ak3TrackJets2+ak4TrackJets2+ak5TrackJets2+ak6TrackJets2
75     +
76     kt1TrackJets05+kt2TrackJets05+kt3TrackJets05+kt4TrackJets05+kt5TrackJets05+kt6TrackJets05
77     +
78     kt1TrackJets1+kt2TrackJets1+kt3TrackJets1+kt4TrackJets1+kt5TrackJets1+kt6TrackJets1
79     +
80     kt1TrackJets2+kt2TrackJets2+kt3TrackJets2+kt4TrackJets2+kt5TrackJets2+kt6TrackJets2
81     )
82    
83    
84    
85    
86    
87    
88    
89    
90    
91    
92