|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
|
import FWCore.ParameterSet.Config as cms
process = cms.Process("LA")
process.load("Configuration.StandardSequences.Services_cff")
#process.load("Configuration.StandardSequences.GeometryDB_cff")
process.load('Configuration/StandardSequences/GeometryExtended_cff')
process.load('Configuration.StandardSequences.MagneticField_cff')
# process.load("Configuration.StandardSequences.FakeConditions_cff")
process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
# check for the correct tag on https://twiki.cern.ch/twiki/bin/view/CMS/SWGuideFrontierConditions
process.GlobalTag.globaltag = "GR09_PV7::All"
process.load("RecoTracker.Configuration.RecoTracker_cff")
from RecoVertex.BeamSpotProducer.BeamSpot_cff import *
process.offlineBeamSpot = offlineBeamSpot
process.load("RecoTracker/TrackProducer/TrackRefitters_cff")
# put here the tag of the tracks you want to use
# alcareco samples have special names for the tracks, in normal reco samples generalTracks can be used
process.TrackRefitter.src = "generalTracks"
# process.TrackRefitter.src = "ALCARECOTkAlZMuMu"
process.TrackRefitter.TrajectoryInEvent = True
process.load("RecoTracker.TransientTrackingRecHit.TransientTrackingRecHitBuilderWithoutRefit_cfi")
process.MessageLogger = cms.Service("MessageLogger",
cerr = cms.untracked.PSet(
enable = cms.untracked.bool(False)
),
cout = cms.untracked.PSet(
enable = cms.untracked.bool(True)
),
files = cms.untracked.PSet(
simul = cms.untracked.PSet(
threshold = cms.untracked.string('ERROR')
)
)
)
process.lorentzAngle = cms.EDAnalyzer("SiPixelLorentzAngle",
src = cms.string("TrackRefitter"),
fileName = cms.string("lorentzangle.root"),
fileNameFit = cms.string("lorentzFit.txt"),
binsDepth = cms.int32(50),
binsDrift = cms.int32(200),
ptMin = cms.double(3),
#in case of MC set this to true to save the simhits (does not work currently, Mixing Module needs to be included correctly)
simData = cms.bool(False),
normChi2Max = cms.double(2),
clustSizeYMin = cms.int32(4),
residualMax = cms.double(0.005),
clustChargeMax = cms.double(120000)
)
process.myout = cms.OutputModule("PoolOutputModule",
fileName = cms.untracked.string('LA_CMSSW.root')
)
process.p = cms.Path(process.offlineBeamSpot*process.TrackRefitter*process.lorentzAngle)
# uncomment this if you want to write out the new CMSSW root file (very large)
# process.outpath = cms.EndPath(process.myout)
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(100)
)
process.source = cms.Source("PoolSource",
#put here the sample you want to use
fileNames = cms.untracked.vstring(
#put your source file here
# ' '
),
# skipEvents = cms.untracked.uint32(100)
)
|