ESProducers/python/ppsAssociationCuts_non_DB_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 from CalibPPS.ESProducers.ppsAssociationCutsESSource_cfi import *
0004
0005 p2016 = cms.PSet(
0006     validityRange=cms.EventRange("273725:min - 284044:max"),
0007     association_cuts_45=cms.PSet(
0008         xi_cut_mean =cms.string("0."),
0009         xi_cut_threshold=cms.string("0.010"),
0010
0011         ti_tr_min=cms.double(-1.5),
0012         ti_tr_max=cms.double(2.0)
0013     ),
0014     association_cuts_56=cms.PSet(
0015         xi_cut_mean= cms.string("0."),
0016         xi_cut_threshold=cms.string("0.015"),
0017
0018         ti_tr_min=cms.double(-1.5),
0019         ti_tr_max=cms.double(2.0)
0020     ),
0021 )
0022 ppsAssociationCutsESSource.configuration.append(p2016)
0023
0024 p2017 = cms.PSet(
0025     validityRange=cms.EventRange("297046:min - 307082:max"),
0026     association_cuts_45=cms.PSet(
0027         xi_cut_mean=cms.string("+6.0695e-5"),
0028         xi_cut_threshold=cms.string("5. * 0.00121"),
0029         ti_tr_min=cms.double(-1.5),
0030         ti_tr_max=cms.double(2.0)
0031     ),
0032     association_cuts_56=cms.PSet(
0033         y_cut_mean=cms.string("-0.022612"),
0034         y_cut_threshold=cms.string("5. * 0.14777"),
0035         xi_cut_mean=cms.string("+8.012857e-5"),
0036         xi_cut_threshold=cms.string("5. * 0.0020627"),
0037         ti_tr_min=cms.double(-1.5),
0038         ti_tr_max=cms.double(2.0)
0039     ),
0040 )
0041 ppsAssociationCutsESSource.configuration.append(p2017)
0042
0043 p2018 = cms.PSet(
0044     validityRange=cms.EventRange("314747:min - 325175:max"),
0045     association_cuts_45=cms.PSet(
0046         x_cut_mean=cms.string("-0.065194856"),
0047         x_cut_threshold=cms.string("4. * 0.16008188"),
0048         y_cut_mean=cms.string("+0.10973631"),
0049         y_cut_threshold=cms.string("4. * 0.1407986"),
0050         xi_cut_mean=cms.string("+3.113062e-5"),
0051         xi_cut_threshold=cms.string("4. * 0.0012403586"),
0052         ti_tr_min=cms.double(-1.5),
0053         ti_tr_max=cms.double(2.0)
0054     ),
0055     association_cuts_56=cms.PSet(
0056         x_cut_mean=cms.string("+0.073016431"),
0057         x_cut_threshold=cms.string("5. * 0.18126434"),
0058         y_cut_mean=cms.string("+0.064261029"),
0059         y_cut_threshold=cms.string("5. * 0.14990802"),
0060         xi_cut_mean=cms.string("-1.1852528e-5"),
0061         xi_cut_threshold=cms.string("5. * 0.002046409"),
0062         ti_tr_min=cms.double(-1.5),
0063         ti_tr_max=cms.double(2.0)
0064     ),
0065 )
0066 ppsAssociationCutsESSource.configuration.append(p2018)
0067
0068 p2021 = cms.PSet(
0069     validityRange=cms.EventRange("1234:1 - 1234:max"), # NB: a fake IOV, there are no LHC data from 2021 - keep it just for MC
0070     association_cuts_45=cms.PSet(
0071         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]) )"),
0072         x_cut_threshold = cms.string("4 * ( (0.091692+0.009316*[x_near]) + (-0.000727)*[y_near]*[y_near] )"),
0073         y_cut_mean = cms.string("- ( (0.672525+-0.619398/[x_near])*TMath::Erf([y_near]/(0.968499+0.597136*[x_near])) )"),
0074         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]) )"),
0075
0076         ti_tr_min=cms.double(-1.5),
0077         ti_tr_max=cms.double(2.0)
0078     ),
0079     association_cuts_56=cms.PSet(
0080         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]) )"),
0081         x_cut_threshold = cms.string("4 * ( (0.091692+0.009316*[x_near]) + (-0.000727)*[y_near]*[y_near] )"),
0082         y_cut_mean = cms.string("- ( (0.672525+-0.619398/[x_near])*TMath::Erf([y_near]/(0.968499+0.597136*[x_near])) )"),
0083         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]) )"),
0084
0085         ti_tr_min=cms.double(-1.5),
0086         ti_tr_max=cms.double(2.0)
0087     ),
0088 )
0089
0090 # NB: do not append the 2021 config - not used for any LHC data
0091
0092 p2022 = cms.PSet(
0093     validityRange=cms.EventRange("343890:min - 999999:max"),
0094     association_cuts_45=cms.PSet(
0095         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]) )"),
0096         x_cut_threshold = cms.string("4 * ( (0.091692+0.009316*[x_near]) + (-0.000727)*[y_near]*[y_near] )"),
0097         y_cut_mean = cms.string("- ( (0.672525+-0.619398/[x_near])*TMath::Erf([y_near]/(0.968499+0.597136*[x_near])) )"),
0098         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]) )"),
0099
0100         ti_tr_min=cms.double(-1.5),
0101         ti_tr_max=cms.double(2.0)
0102     ),
0103     association_cuts_56=cms.PSet(
0104         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]) )"),
0105         x_cut_threshold = cms.string("4 * ( (0.091692+0.009316*[x_near]) + (-0.000727)*[y_near]*[y_near] )"),
0106         y_cut_mean = cms.string("- ( (0.672525+-0.619398/[x_near])*TMath::Erf([y_near]/(0.968499+0.597136*[x_near])) )"),
0107         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]) )"),
0108
0109         ti_tr_min=cms.double(-1.5),
0110         ti_tr_max=cms.double(2.0)
0111     ),
0112 )
0113 ppsAssociationCutsESSource.configuration.append(p2022)
0114
0115 def use_single_infinite_iov_entry(ppsAssociationCutsESSource, iov):
0116     ppsAssociationCutsESSource.configuration = cms.VPSet()
0117     iov.validityRange = cms.EventRange("0:min - 999999:max")
0118     ppsAssociationCutsESSource.configuration.append(iov)