|
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
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
|
import FWCore.ParameterSet.Config as cms
from CalibPPS.ESProducers.ppsAssociationCutsESSource_cfi import *
p2016 = cms.PSet(
validityRange=cms.EventRange("273725:min - 284044:max"),
association_cuts_45=cms.PSet(
xi_cut_mean =cms.string("0."),
xi_cut_threshold=cms.string("0.010"),
ti_tr_min=cms.double(-1.5),
ti_tr_max=cms.double(2.0)
),
association_cuts_56=cms.PSet(
xi_cut_mean= cms.string("0."),
xi_cut_threshold=cms.string("0.015"),
ti_tr_min=cms.double(-1.5),
ti_tr_max=cms.double(2.0)
),
)
ppsAssociationCutsESSource.configuration.append(p2016)
p2017 = cms.PSet(
validityRange=cms.EventRange("297046:min - 307082:max"),
association_cuts_45=cms.PSet(
xi_cut_mean=cms.string("+6.0695e-5"),
xi_cut_threshold=cms.string("5. * 0.00121"),
ti_tr_min=cms.double(-1.5),
ti_tr_max=cms.double(2.0)
),
association_cuts_56=cms.PSet(
y_cut_mean=cms.string("-0.022612"),
y_cut_threshold=cms.string("5. * 0.14777"),
xi_cut_mean=cms.string("+8.012857e-5"),
xi_cut_threshold=cms.string("5. * 0.0020627"),
ti_tr_min=cms.double(-1.5),
ti_tr_max=cms.double(2.0)
),
)
ppsAssociationCutsESSource.configuration.append(p2017)
p2018 = cms.PSet(
validityRange=cms.EventRange("314747:min - 325175:max"),
association_cuts_45=cms.PSet(
x_cut_mean=cms.string("-0.065194856"),
x_cut_threshold=cms.string("4. * 0.16008188"),
y_cut_mean=cms.string("+0.10973631"),
y_cut_threshold=cms.string("4. * 0.1407986"),
xi_cut_mean=cms.string("+3.113062e-5"),
xi_cut_threshold=cms.string("4. * 0.0012403586"),
ti_tr_min=cms.double(-1.5),
ti_tr_max=cms.double(2.0)
),
association_cuts_56=cms.PSet(
x_cut_mean=cms.string("+0.073016431"),
x_cut_threshold=cms.string("5. * 0.18126434"),
y_cut_mean=cms.string("+0.064261029"),
y_cut_threshold=cms.string("5. * 0.14990802"),
xi_cut_mean=cms.string("-1.1852528e-5"),
xi_cut_threshold=cms.string("5. * 0.002046409"),
ti_tr_min=cms.double(-1.5),
ti_tr_max=cms.double(2.0)
),
)
ppsAssociationCutsESSource.configuration.append(p2018)
p2021 = cms.PSet(
validityRange=cms.EventRange("1234:1 - 1234:max"), # NB: a fake IOV, there are no LHC data from 2021 - keep it just for MC
association_cuts_45=cms.PSet(
x_cut_mean = cms.string("- ( (-0.530895+0.112595*[x_near]+-0.006785*[x_near]*[x_near]) - (0.046487+0.179333/[x_near])*abs([y_near]) )"),
x_cut_threshold = cms.string("4 * ( (0.091692+0.009316*[x_near]) + (-0.000727)*[y_near]*[y_near] )"),
y_cut_mean = cms.string("- ( (0.672525+-0.619398/[x_near])*TMath::Erf([y_near]/(0.968499+0.597136*[x_near])) )"),
y_cut_threshold = cms.string("4 * ( (0.015202+0.962944/[x_near]/[x_near]) + (0.009804+-7.710241/pow([x_near],5))*abs([y_near]) )"),
ti_tr_min=cms.double(-1.5),
ti_tr_max=cms.double(2.0)
),
association_cuts_56=cms.PSet(
x_cut_mean = cms.string("- ( (-0.530895+0.112595*[x_near]+-0.006785*[x_near]*[x_near]) - (0.046487+0.179333/[x_near])*abs([y_near]) )"),
x_cut_threshold = cms.string("4 * ( (0.091692+0.009316*[x_near]) + (-0.000727)*[y_near]*[y_near] )"),
y_cut_mean = cms.string("- ( (0.672525+-0.619398/[x_near])*TMath::Erf([y_near]/(0.968499+0.597136*[x_near])) )"),
y_cut_threshold = cms.string("4 * ( (0.015202+0.962944/[x_near]/[x_near]) + (0.009804+-7.710241/pow([x_near],5))*abs([y_near]) )"),
ti_tr_min=cms.double(-1.5),
ti_tr_max=cms.double(2.0)
),
)
# NB: do not append the 2021 config - not used for any LHC data
p2022 = cms.PSet(
validityRange=cms.EventRange("343890:min - 999999:max"),
association_cuts_45=cms.PSet(
x_cut_mean = cms.string("- ( (-0.530895+0.112595*[x_near]+-0.006785*[x_near]*[x_near]) - (0.046487+0.179333/[x_near])*abs([y_near]) )"),
x_cut_threshold = cms.string("4 * ( (0.091692+0.009316*[x_near]) + (-0.000727)*[y_near]*[y_near] )"),
y_cut_mean = cms.string("- ( (0.672525+-0.619398/[x_near])*TMath::Erf([y_near]/(0.968499+0.597136*[x_near])) )"),
y_cut_threshold = cms.string("4 * ( (0.015202+0.962944/[x_near]/[x_near]) + (0.009804+-7.710241/pow([x_near],5))*abs([y_near]) )"),
ti_tr_min=cms.double(-1.5),
ti_tr_max=cms.double(2.0)
),
association_cuts_56=cms.PSet(
x_cut_mean = cms.string("- ( (-0.530895+0.112595*[x_near]+-0.006785*[x_near]*[x_near]) - (0.046487+0.179333/[x_near])*abs([y_near]) )"),
x_cut_threshold = cms.string("4 * ( (0.091692+0.009316*[x_near]) + (-0.000727)*[y_near]*[y_near] )"),
y_cut_mean = cms.string("- ( (0.672525+-0.619398/[x_near])*TMath::Erf([y_near]/(0.968499+0.597136*[x_near])) )"),
y_cut_threshold = cms.string("4 * ( (0.015202+0.962944/[x_near]/[x_near]) + (0.009804+-7.710241/pow([x_near],5))*abs([y_near]) )"),
ti_tr_min=cms.double(-1.5),
ti_tr_max=cms.double(2.0)
),
)
ppsAssociationCutsESSource.configuration.append(p2022)
def use_single_infinite_iov_entry(ppsAssociationCutsESSource, iov):
ppsAssociationCutsESSource.configuration = cms.VPSet()
iov.validityRange = cms.EventRange("0:min - 999999:max")
ppsAssociationCutsESSource.configuration.append(iov)
|