|
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
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
|
#=================================
#inputs
globaltag = 'auto:phase1_2018_design' #APEs are copied from this GT (and IdealGeometry and TrackerTopology are used)
inputsqlitefile = None #if None, uses the GT alignment
alignmenttag = 'Alignments' #tag name for TrackerAlignmentRcd in the input file, also used for the output file
runnumberalignmentIOV = 1 #any run number in the iov that you want to start from
outputfilename = 'outputfile.db'
#misalignment amplitudes, -999 means no misalignment
#the commented numbers are the default magnitudes, which produce a maximum movement of around 600 microns
#see Alignment/TrackerAlignment/plugins/TrackerSystematicMisalignments.cc for definitions
#see also https://twiki.cern.ch/twiki/bin/viewauth/CMS/SystematicMisalignmentsofTracker
radialEpsilon = -999. # 5e-4
telescopeEpsilon = -999. # 5e-4
layerRotEpsilon = -999. # 9.43e-6 #cm^-1
bowingEpsilon = -999. # 6.77e-9 #cm^-2
zExpEpsilon = -999. # 2.02e-4
twistEpsilon = -999. # 2.04e-6 #cm^-1
ellipticalEpsilon = -999. # 5e-4
skewEpsilon = -999. # 5.5e-2 #cm
sagittaEpsilon = -999. # 5.0e-4
#phases for phi dependent misalignments
ellipticalDelta = 0.
skewDelta = 0.
sagittaDelta = 0.
#=================================
import FWCore.ParameterSet.Config as cms
process = cms.Process("TrackerSystematicMisalignments")
process.load("FWCore.MessageService.MessageLogger_cfi")
process.load("Configuration.Geometry.GeometryRecoDB_cff")
process.load("CondCore.CondDB.CondDB_cfi")
process.source = cms.Source("EmptySource",
firstRun=cms.untracked.uint32(runnumberalignmentIOV),
)
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(1)
)
# initial geom
# configure the database file - use survey one for default
process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
from Configuration.AlCa.GlobalTag import GlobalTag
process.GlobalTag = GlobalTag(process.GlobalTag, globaltag, '')
if inputsqlitefile is not None:
process.GlobalTag.toGet = cms.VPSet(
cms.PSet(
record = cms.string('TrackerAlignmentRcd'),
tag = cms.string(alignmenttag),
connect = cms.string('sqlite_file:'+inputsqlitefile),
),
)
# input
process.load("Alignment.TrackerAlignment.TrackerSystematicMisalignments_cfi")
process.TrackerSystematicMisalignments.fromDBGeom = True
#uncomment one or more of these to apply the misalignment(s)
process.TrackerSystematicMisalignments.radialEpsilon = radialEpsilon
process.TrackerSystematicMisalignments.telescopeEpsilon = telescopeEpsilon
process.TrackerSystematicMisalignments.layerRotEpsilon = layerRotEpsilon
process.TrackerSystematicMisalignments.bowingEpsilon = bowingEpsilon
process.TrackerSystematicMisalignments.zExpEpsilon = zExpEpsilon
process.TrackerSystematicMisalignments.twistEpsilon = twistEpsilon
process.TrackerSystematicMisalignments.ellipticalEpsilon = ellipticalEpsilon
process.TrackerSystematicMisalignments.skewEpsilon = skewEpsilon
process.TrackerSystematicMisalignments.sagittaEpsilon = sagittaEpsilon
#misalignment phases
process.TrackerSystematicMisalignments.ellipticalDelta = ellipticalDelta
process.TrackerSystematicMisalignments.skewDelta = skewDelta
process.TrackerSystematicMisalignments.sagittaDelta = sagittaDelta
# output
process.CondDB.connect = 'sqlite_file:'+outputfilename
process.PoolDBOutputService = cms.Service("PoolDBOutputService",
process.CondDB,
toPut = cms.VPSet(cms.PSet(record = cms.string('TrackerAlignmentRcd'),
tag = cms.string(alignmenttag),
),
cms.PSet(record = cms.string('TrackerAlignmentErrorExtendedRcd'),
tag = cms.string('AlignmentErrorsExtended'),
),
),
)
process.p = cms.Path( process.TrackerSystematicMisalignments )
|