|
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
|
## adapted from produceCalibrationTree_template_cfg.py
import FWCore.ParameterSet.Config as cms
##from CalibTracker.SiStripCommon.shallowTree_test_template import * ## TODO get rid of this one
process = cms.Process('CALIB')
process.load('Configuration/StandardSequences/MagneticField_cff')
process.load('Configuration.Geometry.GeometryRecoDB_cff')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
from Configuration.AlCa.GlobalTag import GlobalTag
process.GlobalTag = GlobalTag(process.GlobalTag, "auto:run3_data_PromptAnalysis")
process.load('FWCore.MessageService.MessageLogger_cfi')
process.load('Configuration.StandardSequences.Services_cff')
process.maxEvents = cms.untracked.PSet(input=cms.untracked.int32(-1))
process.source = cms.Source("PoolSource",
fileNames = cms.untracked.vstring("/store/express/Run2023F/StreamExpress/ALCARECO/SiStripCalMinBias-Express-v1/000/373/710/00000/e2df2f78-b95a-4f33-ae22-add59aa2903f.root"))
process.options = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) )
process.MessageLogger.cerr.FwkReport.reportEvery = 10000
inTracks = cms.InputTag("ALCARECOSiStripCalMinBias")
process.load('CalibTracker.SiStripCommon.prescaleEvent_cfi')
process.load('CalibTracker.Configuration.Filter_Refit_cff')
## use CalibrationTracks (for clusters) and CalibrationTracksRefit (for tracks)
process.CalibrationTracks.src = inTracks
tracksForCalib = cms.InputTag("CalibrationTracksRefit")
process.prescaleEvent.prescale = 1
process.TkCalSeq = cms.Sequence(process.prescaleEvent*process.MeasurementTrackerEvent*process.trackFilterRefit)
process.load("PhysicsTools.NanoAOD.nano_cff")
process.load("PhysicsTools.NanoAOD.NanoAODEDMEventContent_cff")
## as a test: it should be possible to add tracks fully at configuration level (+ declaring the plugin)
from PhysicsTools.NanoAOD.common_cff import *
## this is equivalent to ShallowTrackProducer as configured for the gain calibration
process.tracksTable = cms.EDProducer("SimpleTrackFlatTableProducer",
src=tracksForCalib,
cut=cms.string(""),
name=cms.string("track"),
doc=cms.string("SiStripCalMinBias ALCARECO tracks"),
singleton=cms.bool(False),
extension=cms.bool(False),
variables=cms.PSet(
chi2ndof=Var("chi2()/ndof", float),
pt=Var("pt()", float),
hitsvalid=Var("numberOfValidHits()", int), ## unsigned?
phi=Var("phi()", float),
eta=Var("eta()", float),
)
)
process.load("CalibTracker.SiStripCommon.tkInstLumiTable_cfi")
process.tkInstLumiTable.extension = True
process.load("CalibTracker.SiStripCommon.siStripGainCalibTable_cfi")
process.siStripGainCalibTable.Tracks = tracksForCalib
process.nanoCTPath = cms.Path(process.TkCalSeq*
process.nanoMetadata*process.tkInstLumiTable
*process.tracksTable
*process.siStripGainCalibTable
)
process.out = cms.OutputModule("NanoAODOutputModule",
fileName=cms.untracked.string("CalibTreeMC_nano_G2.root"),
outputCommands=process.NANOAODEventContent.outputCommands
)
process.end = cms.EndPath(process.out)
|