upBounds=cms.vdouble(2.000, 999.000),
                 values=cms.vdouble(0.0165),
                 uncertainties=cms.vdouble(0.0015)),
        cms.PSet(lowBounds=cms.vdouble(2.000, -999.000),
                 upBounds=cms.vdouble(6.000, 0.940),
                 values=cms.vdouble(0.0215),
                 uncertainties=cms.vdouble(0.0010)),
        cms.PSet(lowBounds=cms.vdouble(2.000, 0.940),
                 upBounds=cms.vdouble(6.000, 999.000),
                 values=cms.vdouble(0.0194),
                 uncertainties=cms.vdouble(0.0008)),
    ))

from flashgg.Systematics.ecalElfTools import getRunDependentScaleBins

photonScaleBinsData = getRunDependentScaleBins(
    'flashgg/Systematics/data/test_scales_800_plus_DCSOnly.dat')

photonScaleUncertBins = cms.PSet(
    variables=cms.vstring("abs(superCluster.eta)", "r9"),
    bins=cms.VPSet(
        ## only uncertainties here. scales loaded by SystematicsCustomize according to process.photonScaleBinsData
        cms.PSet(lowBounds=cms.vdouble(0.000, -999.000),
                 upBounds=cms.vdouble(1.000, 0.940),
                 values=cms.vdouble(0.),
                 uncertainties=cms.vdouble(0.00050)),
        cms.PSet(lowBounds=cms.vdouble(0.000, 0.940),
                 upBounds=cms.vdouble(1.000, 999.000),
                 values=cms.vdouble(0.),
                 uncertainties=cms.vdouble(0.00050)),
        cms.PSet(lowBounds=cms.vdouble(1.000, -999.000),
                 upBounds=cms.vdouble(1.500, 0.940),
photonSmearBins = cms.PSet(
    variables = cms.vstring("abs(superCluster.eta)","r9"),
    bins = cms.VPSet(
        cms.PSet(lowBounds=cms.vdouble(0.000,0.940),    upBounds=cms.vdouble(1.000,999.000), values=cms.vdouble(0.0083), uncertainties=cms.vdouble(0.0002)),
        cms.PSet(lowBounds=cms.vdouble(0.000,-999.000), upBounds=cms.vdouble(1.000,0.940),   values=cms.vdouble(0.0095), uncertainties=cms.vdouble(0.0002)),
        cms.PSet(lowBounds=cms.vdouble(1.000,0.940),    upBounds=cms.vdouble(1.500,999.000), values=cms.vdouble(0.0105), uncertainties=cms.vdouble(0.0008)),
        cms.PSet(lowBounds=cms.vdouble(1.000,-999.000), upBounds=cms.vdouble(1.500,0.940),   values=cms.vdouble(0.0159), uncertainties=cms.vdouble(0.0002)),
        cms.PSet(lowBounds=cms.vdouble(1.500,0.940),    upBounds=cms.vdouble(2.000,999.000), values=cms.vdouble(0.0182), uncertainties=cms.vdouble(0.0005)),
        cms.PSet(lowBounds=cms.vdouble(1.500,-999.000), upBounds=cms.vdouble(2.000,0.940),   values=cms.vdouble(0.0218), uncertainties=cms.vdouble(0.0003)),
        cms.PSet(lowBounds=cms.vdouble(2.000,0.940),    upBounds=cms.vdouble(6.000,999.000), values=cms.vdouble(0.0229), uncertainties=cms.vdouble(0.0003)),
        cms.PSet(lowBounds=cms.vdouble(2.000,-999.000), upBounds=cms.vdouble(6.000,0.940),   values=cms.vdouble(0.0247), uncertainties=cms.vdouble(0.0003)),
        ))

from flashgg.Systematics.ecalElfTools import getRunDependentScaleBins

photonScaleBinsData = getRunDependentScaleBins('flashgg/Systematics/data/80X_Golden15July_plusDCS_scales.dat')

photonScaleUncertBins = cms.PSet(
    variables = cms.vstring("abs(superCluster.eta)","r9"),
    bins = cms.VPSet(
        ## only uncertainties here. scales loaded by SystematicsCustomize according to process.photonScaleBinsData
            cms.PSet( lowBounds = cms.vdouble(0.000 , -999.000 ) , upBounds = cms.vdouble(1.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00050 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(0.000 , 0.940    ) , upBounds = cms.vdouble(1.000 , 999.000 ) , values = cms.vdouble( 0.  ) , uncertainties = cms.vdouble( 0.00050 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.000 , -999.000 ) , upBounds = cms.vdouble(1.500 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00120 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.000 , 0.940    ) , upBounds = cms.vdouble(1.500 , 999.000 ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00060 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.500 , -999.000 ) , upBounds = cms.vdouble(2.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00200 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.500 , 0.940    ) , upBounds = cms.vdouble(2.000 , 999.000 ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00300 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(2.000 , -999.000 ) , upBounds = cms.vdouble(6.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00200 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(2.000 , 0.940    ) , upBounds = cms.vdouble(6.000 , 999.000 ) , values = cms.vdouble( 0.  ) , uncertainties = cms.vdouble( 0.00300 )   ) ,
                    ))
Пример #3
0
photonSmearBins = cms.PSet(
    variables = cms.vstring("abs(superCluster.eta)","r9"),
    bins = cms.VPSet(
        cms.PSet(lowBounds=cms.vdouble(0.000,-999.000), upBounds=cms.vdouble(1.000,0.940),   values=cms.vdouble(0.0088), uncertainties=cms.vdouble(0.0003)),
        cms.PSet(lowBounds=cms.vdouble(0.000,0.940),    upBounds=cms.vdouble(1.000,999.000), values=cms.vdouble(0.0075), uncertainties=cms.vdouble(0.0003)),
        cms.PSet(lowBounds=cms.vdouble(1.000,-999.000), upBounds=cms.vdouble(1.500,0.940),   values=cms.vdouble(0.0148), uncertainties=cms.vdouble(0.0004)),
        cms.PSet(lowBounds=cms.vdouble(1.000,0.940),    upBounds=cms.vdouble(1.500,999.000), values=cms.vdouble(0.0116), uncertainties=cms.vdouble(0.0012)),
        cms.PSet(lowBounds=cms.vdouble(1.500,-999.000), upBounds=cms.vdouble(2.000,0.940),   values=cms.vdouble(0.0192), uncertainties=cms.vdouble(0.0006)),
        cms.PSet(lowBounds=cms.vdouble(1.500,0.940),    upBounds=cms.vdouble(2.000,999.000), values=cms.vdouble(0.0170), uncertainties=cms.vdouble(0.0009)),
        cms.PSet(lowBounds=cms.vdouble(2.000,-999.000), upBounds=cms.vdouble(6.000,0.940),   values=cms.vdouble(0.0220), uncertainties=cms.vdouble(0.0006)),
        cms.PSet(lowBounds=cms.vdouble(2.000,0.940),    upBounds=cms.vdouble(6.000,999.000), values=cms.vdouble(0.0200), uncertainties=cms.vdouble(0.0004)),
        ))

from flashgg.Systematics.ecalElfTools import getRunDependentScaleBins

photonScaleBinsData = getRunDependentScaleBins('flashgg/Systematics/data/Golden22June_scales.dat')

photonScaleUncertBins = cms.PSet(
    variables = cms.vstring("abs(superCluster.eta)","r9"),
    bins = cms.VPSet(
        ## only uncertainties here. scales loaded by SystematicsCustomize according to process.photonScaleBinsData
            cms.PSet( lowBounds = cms.vdouble(0.000 , -999.000 ) , upBounds = cms.vdouble(1.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00050 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(0.000 , 0.940    ) , upBounds = cms.vdouble(1.000 , 999.000 ) , values = cms.vdouble( 0.  ) , uncertainties = cms.vdouble( 0.00050 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.000 , -999.000 ) , upBounds = cms.vdouble(1.500 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00120 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.000 , 0.940    ) , upBounds = cms.vdouble(1.500 , 999.000 ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00060 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.500 , -999.000 ) , upBounds = cms.vdouble(2.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00200 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.500 , 0.940    ) , upBounds = cms.vdouble(2.000 , 999.000 ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00300 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(2.000 , -999.000 ) , upBounds = cms.vdouble(6.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00200 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(2.000 , 0.940    ) , upBounds = cms.vdouble(6.000 , 999.000 ) , values = cms.vdouble( 0.  ) , uncertainties = cms.vdouble( 0.00300 )   ) ,
                    ))
Пример #4
0
                 values=cms.vdouble(0.0078),
                 uncertainties=cms.vdouble(0.0046)),
        cms.PSet(lowBounds=cms.vdouble(2.000, -999.000),
                 upBounds=cms.vdouble(6.000, 0.940),
                 values=cms.vdouble(0.0200),
                 uncertainties=cms.vdouble(0.0024)),
        cms.PSet(lowBounds=cms.vdouble(2.000, 0.940),
                 upBounds=cms.vdouble(6.000, 999.000),
                 values=cms.vdouble(0.0193),
                 uncertainties=cms.vdouble(0.0013)),
    ))

from flashgg.Systematics.ecalElfTools import getRunDependentScaleBins

# photonScaleBinsData = getRunDependentScaleBins('EgammaAnalysis/ElectronTools/data/76X_16DecRereco_2015_scales.dat')
photonScaleBinsData = getRunDependentScaleBins(
    'flashgg/Systematics/data/test_2016B_scales.dat')

photonScaleUncertBins = cms.PSet(
    variables=cms.vstring("abs(superCluster.eta)", "r9"),
    bins=cms.VPSet(
        ## only uncertainties here. scales loaded by SystematicsCustomize according to process.photonScaleBinsData
        cms.PSet(lowBounds=cms.vdouble(0.000, -999.000),
                 upBounds=cms.vdouble(1.000, 0.940),
                 values=cms.vdouble(0.),
                 uncertainties=cms.vdouble(0.00050)),
        cms.PSet(lowBounds=cms.vdouble(0.000, 0.940),
                 upBounds=cms.vdouble(1.000, 999.000),
                 values=cms.vdouble(0.),
                 uncertainties=cms.vdouble(0.00050)),
        cms.PSet(lowBounds=cms.vdouble(1.000, -999.000),
                 upBounds=cms.vdouble(1.500, 0.940),
Пример #5
0
            # smearings for rereco 2015, photons,
            # see Fasanella et al., ECAL DPG 17/12/2015, https://indico.cern.ch/event/402608/contribution/1/attachments/1206871/1758949/DPG_reReco_corrections.pdf
        cms.PSet(lowBounds=cms.vdouble(0.000,-999.000), upBounds=cms.vdouble(1.000,0.940),   values=cms.vdouble(0.0111), uncertainties=cms.vdouble(0.0005)),
        cms.PSet(lowBounds=cms.vdouble(0.000,0.940),    upBounds=cms.vdouble(1.000,999.000), values=cms.vdouble(0.0096), uncertainties=cms.vdouble(0.0005)),
        cms.PSet(lowBounds=cms.vdouble(1.000,-999.000), upBounds=cms.vdouble(1.500,0.940),   values=cms.vdouble(0.0148), uncertainties=cms.vdouble(0.0009)),
        cms.PSet(lowBounds=cms.vdouble(1.000,0.940),    upBounds=cms.vdouble(1.500,999.000), values=cms.vdouble(0.0143), uncertainties=cms.vdouble(0.0045)),
        cms.PSet(lowBounds=cms.vdouble(1.500,-999.000), upBounds=cms.vdouble(2.000,0.940),   values=cms.vdouble(0.002), uncertainties=cms.vdouble(0.0056)),
        cms.PSet(lowBounds=cms.vdouble(1.500,0.940),    upBounds=cms.vdouble(2.000,999.000), values=cms.vdouble(0.0078), uncertainties=cms.vdouble(0.0046)),
        cms.PSet(lowBounds=cms.vdouble(2.000,-999.000), upBounds=cms.vdouble(6.000,0.940),   values=cms.vdouble(0.0200), uncertainties=cms.vdouble(0.0024)),
        cms.PSet(lowBounds=cms.vdouble(2.000,0.940),    upBounds=cms.vdouble(6.000,999.000), values=cms.vdouble(0.0193), uncertainties=cms.vdouble(0.0013)),
        ))

from flashgg.Systematics.ecalElfTools import getRunDependentScaleBins

# photonScaleBinsData = getRunDependentScaleBins('EgammaAnalysis/ElectronTools/data/76X_16DecRereco_2015_scales.dat')
photonScaleBinsData = getRunDependentScaleBins('flashgg/Systematics/data/test_2016B_scales.dat')

photonScaleUncertBins = cms.PSet(
    variables = cms.vstring("abs(superCluster.eta)","r9"),
    bins = cms.VPSet(
        ## only uncertainties here. scales loaded by SystematicsCustomize according to process.photonScaleBinsData
            cms.PSet( lowBounds = cms.vdouble(0.000 , -999.000 ) , upBounds = cms.vdouble(1.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00050 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(0.000 , 0.940    ) , upBounds = cms.vdouble(1.000 , 999.000 ) , values = cms.vdouble( 0.  ) , uncertainties = cms.vdouble( 0.00050 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.000 , -999.000 ) , upBounds = cms.vdouble(1.500 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00120 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.000 , 0.940    ) , upBounds = cms.vdouble(1.500 , 999.000 ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00060 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.500 , -999.000 ) , upBounds = cms.vdouble(2.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00200 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.500 , 0.940    ) , upBounds = cms.vdouble(2.000 , 999.000 ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00300 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(2.000 , -999.000 ) , upBounds = cms.vdouble(6.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00200 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(2.000 , 0.940    ) , upBounds = cms.vdouble(6.000 , 999.000 ) , values = cms.vdouble( 0.  ) , uncertainties = cms.vdouble( 0.00300 )   ) ,
                    ))
Пример #6
0
                 values=cms.vdouble(0.0221),
                 uncertainties=cms.vdouble(0.0006)),
        cms.PSet(lowBounds=cms.vdouble(2.000, -999.000),
                 upBounds=cms.vdouble(6.000, 0.940),
                 values=cms.vdouble(0.0268),
                 uncertainties=cms.vdouble(0.0008)),
        cms.PSet(lowBounds=cms.vdouble(2.000, 0.940),
                 upBounds=cms.vdouble(6.000, 999.000),
                 values=cms.vdouble(0.0230),
                 uncertainties=cms.vdouble(0.0010)),
    ))

from flashgg.Systematics.ecalElfTools import getRunDependentScaleBins

# photonScaleBinsData = getRunDependentScaleBins('EgammaAnalysis/ElectronTools/data/76X_16DecRereco_2015_scales.dat')
photonScaleBinsData = getRunDependentScaleBins(
    'flashgg/Systematics/data/76X_16DecRereco_2015_scales.dat')

photonScaleUncertBins = cms.PSet(
    variables=cms.vstring("abs(superCluster.eta)", "r9"),
    bins=cms.VPSet(
        ## only uncertainties here. scales loaded by SystematicsCustomize according to process.photonScaleBinsData
        cms.PSet(lowBounds=cms.vdouble(0.000, -999.000),
                 upBounds=cms.vdouble(1.000, 0.940),
                 values=cms.vdouble(0.),
                 uncertainties=cms.vdouble(0.00050)),
        cms.PSet(lowBounds=cms.vdouble(0.000, 0.940),
                 upBounds=cms.vdouble(1.000, 999.000),
                 values=cms.vdouble(0.),
                 uncertainties=cms.vdouble(0.00050)),
        cms.PSet(lowBounds=cms.vdouble(1.000, -999.000),
                 upBounds=cms.vdouble(1.500, 0.940),
Пример #7
0
                                        upBounds=cms.vdouble(2.000, 999.000),
                                        values=cms.vdouble(0.0170),
                                        uncertainties=cms.vdouble(0.0009)),
                               cms.PSet(lowBounds=cms.vdouble(2.000, -999.000),
                                        upBounds=cms.vdouble(6.000, 0.940),
                                        values=cms.vdouble(0.0220),
                                        uncertainties=cms.vdouble(0.0006)),
                               cms.PSet(lowBounds=cms.vdouble(2.000, 0.940),
                                        upBounds=cms.vdouble(6.000, 999.000),
                                        values=cms.vdouble(0.0200),
                                        uncertainties=cms.vdouble(0.0004)),
                           ))

from flashgg.Systematics.ecalElfTools import getRunDependentScaleBins

photonScaleBinsData = getRunDependentScaleBins(
    'flashgg/Systematics/data/Golden22June_scales.dat')

photonScaleUncertBins = cms.PSet(
    variables=cms.vstring("abs(superCluster.eta)", "r9"),
    bins=cms.VPSet(
        ## only uncertainties here. scales loaded by SystematicsCustomize according to process.photonScaleBinsData
        cms.PSet(lowBounds=cms.vdouble(0.000, -999.000),
                 upBounds=cms.vdouble(1.000, 0.940),
                 values=cms.vdouble(0.),
                 uncertainties=cms.vdouble(0.00050)),
        cms.PSet(lowBounds=cms.vdouble(0.000, 0.940),
                 upBounds=cms.vdouble(1.000, 999.000),
                 values=cms.vdouble(0.),
                 uncertainties=cms.vdouble(0.00050)),
        cms.PSet(lowBounds=cms.vdouble(1.000, -999.000),
                 upBounds=cms.vdouble(1.500, 0.940),
Пример #8
0
                 upBounds=cms.vdouble(2.000, 999.000),
                 values=cms.vdouble(0.0171),
                 uncertainties=cms.vdouble(0.0012)),
        cms.PSet(lowBounds=cms.vdouble(2.000, -999.000),
                 upBounds=cms.vdouble(6.000, 0.940),
                 values=cms.vdouble(0.0211),
                 uncertainties=cms.vdouble(0.0007)),
        cms.PSet(lowBounds=cms.vdouble(2.000, 0.940),
                 upBounds=cms.vdouble(6.000, 999.000),
                 values=cms.vdouble(0.0189),
                 uncertainties=cms.vdouble(0.0006)),
    ))

from flashgg.Systematics.ecalElfTools import getRunDependentScaleBins

photonScaleBinsData = getRunDependentScaleBins(
    'flashgg/Systematics/data/Golden10June_plus_DCS.dat')

photonScaleUncertBins = cms.PSet(
    variables=cms.vstring("abs(superCluster.eta)", "r9"),
    bins=cms.VPSet(
        ## only uncertainties here. scales loaded by SystematicsCustomize according to process.photonScaleBinsData
        cms.PSet(lowBounds=cms.vdouble(0.000, -999.000),
                 upBounds=cms.vdouble(1.000, 0.940),
                 values=cms.vdouble(0.),
                 uncertainties=cms.vdouble(0.00050)),
        cms.PSet(lowBounds=cms.vdouble(0.000, 0.940),
                 upBounds=cms.vdouble(1.000, 999.000),
                 values=cms.vdouble(0.),
                 uncertainties=cms.vdouble(0.00050)),
        cms.PSet(lowBounds=cms.vdouble(1.000, -999.000),
                 upBounds=cms.vdouble(1.500, 0.940),
                                        upBounds=cms.vdouble(2.000, 0.940),
                                        values=cms.vdouble(0.0218),
                                        uncertainties=cms.vdouble(0.0003)),
                               cms.PSet(lowBounds=cms.vdouble(2.000, 0.940),
                                        upBounds=cms.vdouble(6.000, 999.000),
                                        values=cms.vdouble(0.0229),
                                        uncertainties=cms.vdouble(0.0003)),
                               cms.PSet(lowBounds=cms.vdouble(2.000, -999.000),
                                        upBounds=cms.vdouble(6.000, 0.940),
                                        values=cms.vdouble(0.0247),
                                        uncertainties=cms.vdouble(0.0003)),
                           ))

from flashgg.Systematics.ecalElfTools import getRunDependentScaleBins

photonScaleBinsData = getRunDependentScaleBins(
    'flashgg/Systematics/data/80X_Golden15July_plusDCS_scales.dat')

photonScaleUncertBins = cms.PSet(
    variables=cms.vstring("abs(superCluster.eta)", "r9"),
    bins=cms.VPSet(
        ## only uncertainties here. scales loaded by SystematicsCustomize according to process.photonScaleBinsData
        cms.PSet(lowBounds=cms.vdouble(0.000, -999.000),
                 upBounds=cms.vdouble(1.000, 0.940),
                 values=cms.vdouble(0.),
                 uncertainties=cms.vdouble(0.00050)),
        cms.PSet(lowBounds=cms.vdouble(0.000, 0.940),
                 upBounds=cms.vdouble(1.000, 999.000),
                 values=cms.vdouble(0.),
                 uncertainties=cms.vdouble(0.00050)),
        cms.PSet(lowBounds=cms.vdouble(1.000, -999.000),
                 upBounds=cms.vdouble(1.500, 0.940),
Пример #10
0
photonSmearBins = cms.PSet(
    variables = cms.vstring("abs(superCluster.eta)","r9"),
    bins = cms.VPSet(
        cms.PSet(lowBounds=cms.vdouble(0.000,0.940),    upBounds=cms.vdouble(1.000,999.000), values=cms.vdouble(0.0084), uncertainties=cms.vdouble(0.0001)),
        cms.PSet(lowBounds=cms.vdouble(0.000,-999.000), upBounds=cms.vdouble(1.000,0.940),   values=cms.vdouble(0.0097), uncertainties=cms.vdouble(0.0001)),
        cms.PSet(lowBounds=cms.vdouble(1.000,0.940),    upBounds=cms.vdouble(1.500,999.000), values=cms.vdouble(0.0108), uncertainties=cms.vdouble(0.0008)),
        cms.PSet(lowBounds=cms.vdouble(1.000,-999.000), upBounds=cms.vdouble(1.500,0.940),   values=cms.vdouble(0.0162), uncertainties=cms.vdouble(0.0002)),
        cms.PSet(lowBounds=cms.vdouble(1.500,0.940),    upBounds=cms.vdouble(2.000,999.000), values=cms.vdouble(0.0182), uncertainties=cms.vdouble(0.0006)),
        cms.PSet(lowBounds=cms.vdouble(1.500,-999.000), upBounds=cms.vdouble(2.000,0.940),   values=cms.vdouble(0.0220), uncertainties=cms.vdouble(0.0003)),
        cms.PSet(lowBounds=cms.vdouble(2.000,0.940),    upBounds=cms.vdouble(6.000,999.000), values=cms.vdouble(0.0229), uncertainties=cms.vdouble(0.0003)),
        cms.PSet(lowBounds=cms.vdouble(2.000,-999.000), upBounds=cms.vdouble(6.000,0.940),   values=cms.vdouble(0.0251), uncertainties=cms.vdouble(0.0003)),
        ))

from flashgg.Systematics.ecalElfTools import getRunDependentScaleBins

photonScaleBinsData = getRunDependentScaleBins('flashgg/Systematics/data/80X_ichep_2016_pho_scales.dat')

photonScaleUncertBins = cms.PSet(
    variables = cms.vstring("abs(superCluster.eta)","r9"),
    bins = cms.VPSet(
        ## only uncertainties here. scales loaded by SystematicsCustomize according to process.photonScaleBinsData
            cms.PSet( lowBounds = cms.vdouble(0.000 , -999.000 ) , upBounds = cms.vdouble(1.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00050 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(0.000 , 0.940    ) , upBounds = cms.vdouble(1.000 , 999.000 ) , values = cms.vdouble( 0.  ) , uncertainties = cms.vdouble( 0.00050 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.000 , -999.000 ) , upBounds = cms.vdouble(1.500 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00120 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.000 , 0.940    ) , upBounds = cms.vdouble(1.500 , 999.000 ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00060 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.500 , -999.000 ) , upBounds = cms.vdouble(2.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00200 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.500 , 0.940    ) , upBounds = cms.vdouble(2.000 , 999.000 ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00300 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(2.000 , -999.000 ) , upBounds = cms.vdouble(6.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00200 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(2.000 , 0.940    ) , upBounds = cms.vdouble(6.000 , 999.000 ) , values = cms.vdouble( 0.  ) , uncertainties = cms.vdouble( 0.00300 )   ) ,
                    ))
Пример #11
0
    bins = cms.VPSet(
            # smearings for rereco 2015, photons,
            # see Fasanella et al., ECAL DPG 17/12/2015, https://indico.cern.ch/event/402608/contribution/1/attachments/1206871/1758949/DPG_reReco_corrections.pdf
        cms.PSet(lowBounds=cms.vdouble(0.000,-999.000), upBounds=cms.vdouble(1.000,0.940),   values=cms.vdouble(0.0099), uncertainties=cms.vdouble(0.0004)),
        cms.PSet(lowBounds=cms.vdouble(0.000,0.940),    upBounds=cms.vdouble(1.000,999.000), values=cms.vdouble(0.0086), uncertainties=cms.vdouble(0.0004)),
        cms.PSet(lowBounds=cms.vdouble(1.000,-999.000), upBounds=cms.vdouble(1.500,0.940),   values=cms.vdouble(0.0146), uncertainties=cms.vdouble(0.0006)),
        cms.PSet(lowBounds=cms.vdouble(1.000,0.940),    upBounds=cms.vdouble(1.500,999.000), values=cms.vdouble(0.0123), uncertainties=cms.vdouble(0.0019)),
        cms.PSet(lowBounds=cms.vdouble(1.500,-999.000), upBounds=cms.vdouble(2.000,0.940),   values=cms.vdouble(0.0167), uncertainties=cms.vdouble(0.0010)),
        cms.PSet(lowBounds=cms.vdouble(1.500,0.940),    upBounds=cms.vdouble(2.000,999.000), values=cms.vdouble(0.0165), uncertainties=cms.vdouble(0.0015)),
        cms.PSet(lowBounds=cms.vdouble(2.000,-999.000), upBounds=cms.vdouble(6.000,0.940),   values=cms.vdouble(0.0215), uncertainties=cms.vdouble(0.0010)),
        cms.PSet(lowBounds=cms.vdouble(2.000,0.940),    upBounds=cms.vdouble(6.000,999.000), values=cms.vdouble(0.0194), uncertainties=cms.vdouble(0.0008)),
        ))

from flashgg.Systematics.ecalElfTools import getRunDependentScaleBins

photonScaleBinsData = getRunDependentScaleBins('flashgg/Systematics/data/test_scales_800_plus_DCSOnly.dat')

photonScaleUncertBins = cms.PSet(
    variables = cms.vstring("abs(superCluster.eta)","r9"),
    bins = cms.VPSet(
        ## only uncertainties here. scales loaded by SystematicsCustomize according to process.photonScaleBinsData
            cms.PSet( lowBounds = cms.vdouble(0.000 , -999.000 ) , upBounds = cms.vdouble(1.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00050 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(0.000 , 0.940    ) , upBounds = cms.vdouble(1.000 , 999.000 ) , values = cms.vdouble( 0.  ) , uncertainties = cms.vdouble( 0.00050 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.000 , -999.000 ) , upBounds = cms.vdouble(1.500 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00120 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.000 , 0.940    ) , upBounds = cms.vdouble(1.500 , 999.000 ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00060 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.500 , -999.000 ) , upBounds = cms.vdouble(2.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00200 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.500 , 0.940    ) , upBounds = cms.vdouble(2.000 , 999.000 ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00300 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(2.000 , -999.000 ) , upBounds = cms.vdouble(6.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00200 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(2.000 , 0.940    ) , upBounds = cms.vdouble(6.000 , 999.000 ) , values = cms.vdouble( 0.  ) , uncertainties = cms.vdouble( 0.00300 )   ) ,
                    ))
            # smearings for rereco 2015, photons,
            # see Fasanella et al., ECAL DPG 17/12/2015, https://indico.cern.ch/event/402608/contribution/1/attachments/1206871/1758949/DPG_reReco_corrections.pdf
        cms.PSet(lowBounds=cms.vdouble(0.000,-999.000), upBounds=cms.vdouble(1.000,0.940),   values=cms.vdouble(0.0094), uncertainties=cms.vdouble(0.0005)),
        cms.PSet(lowBounds=cms.vdouble(0.000,0.940),    upBounds=cms.vdouble(1.000,999.000), values=cms.vdouble(0.0080), uncertainties=cms.vdouble(0.0012)),
        cms.PSet(lowBounds=cms.vdouble(1.000,-999.000), upBounds=cms.vdouble(1.500,0.940),   values=cms.vdouble(0.0183), uncertainties=cms.vdouble(0.0021)),
        cms.PSet(lowBounds=cms.vdouble(1.000,0.940),    upBounds=cms.vdouble(1.500,999.000), values=cms.vdouble(0.0115), uncertainties=cms.vdouble(0.0005)),
        cms.PSet(lowBounds=cms.vdouble(1.500,-999.000), upBounds=cms.vdouble(2.000,0.940),   values=cms.vdouble(0.0201), uncertainties=cms.vdouble(0.0016)),
        cms.PSet(lowBounds=cms.vdouble(1.500,0.940),    upBounds=cms.vdouble(2.000,999.000), values=cms.vdouble(0.0221), uncertainties=cms.vdouble(0.0006)),
        cms.PSet(lowBounds=cms.vdouble(2.000,-999.000), upBounds=cms.vdouble(6.000,0.940),   values=cms.vdouble(0.0268), uncertainties=cms.vdouble(0.0008)),
        cms.PSet(lowBounds=cms.vdouble(2.000,0.940),    upBounds=cms.vdouble(6.000,999.000), values=cms.vdouble(0.0230), uncertainties=cms.vdouble(0.0010)),
        ))

from flashgg.Systematics.ecalElfTools import getRunDependentScaleBins

# photonScaleBinsData = getRunDependentScaleBins('EgammaAnalysis/ElectronTools/data/76X_16DecRereco_2015_scales.dat')
photonScaleBinsData = getRunDependentScaleBins('flashgg/Systematics/data/76X_16DecRereco_2015_scales.dat')

photonScaleUncertBins = cms.PSet(
    variables = cms.vstring("abs(superCluster.eta)","r9"),
    bins = cms.VPSet(
        ## only uncertainties here. scales loaded by SystematicsCustomize according to process.photonScaleBinsData
            cms.PSet( lowBounds = cms.vdouble(0.000 , -999.000 ) , upBounds = cms.vdouble(1.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00050 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(0.000 , 0.940    ) , upBounds = cms.vdouble(1.000 , 999.000 ) , values = cms.vdouble( 0.  ) , uncertainties = cms.vdouble( 0.00050 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.000 , -999.000 ) , upBounds = cms.vdouble(1.500 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00120 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.000 , 0.940    ) , upBounds = cms.vdouble(1.500 , 999.000 ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00060 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.500 , -999.000 ) , upBounds = cms.vdouble(2.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00200 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.500 , 0.940    ) , upBounds = cms.vdouble(2.000 , 999.000 ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00300 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(2.000 , -999.000 ) , upBounds = cms.vdouble(6.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00200 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(2.000 , 0.940    ) , upBounds = cms.vdouble(6.000 , 999.000 ) , values = cms.vdouble( 0.  ) , uncertainties = cms.vdouble( 0.00300 )   ) ,
                    ))
Пример #13
0
                                        upBounds=cms.vdouble(2.000, 0.940),
                                        values=cms.vdouble(0.0220),
                                        uncertainties=cms.vdouble(0.0003)),
                               cms.PSet(lowBounds=cms.vdouble(2.000, 0.940),
                                        upBounds=cms.vdouble(6.000, 999.000),
                                        values=cms.vdouble(0.0229),
                                        uncertainties=cms.vdouble(0.0003)),
                               cms.PSet(lowBounds=cms.vdouble(2.000, -999.000),
                                        upBounds=cms.vdouble(6.000, 0.940),
                                        values=cms.vdouble(0.0251),
                                        uncertainties=cms.vdouble(0.0003)),
                           ))

from flashgg.Systematics.ecalElfTools import getRunDependentScaleBins

photonScaleBinsData = getRunDependentScaleBins(
    'flashgg/Systematics/data/80X_ichep_2016_pho_scales.dat')

photonScaleUncertBins = cms.PSet(
    variables=cms.vstring("abs(superCluster.eta)", "r9"),
    bins=cms.VPSet(
        ## only uncertainties here. scales loaded by SystematicsCustomize according to process.photonScaleBinsData
        cms.PSet(lowBounds=cms.vdouble(0.000, -999.000),
                 upBounds=cms.vdouble(1.000, 0.940),
                 values=cms.vdouble(0.),
                 uncertainties=cms.vdouble(0.00050)),
        cms.PSet(lowBounds=cms.vdouble(0.000, 0.940),
                 upBounds=cms.vdouble(1.000, 999.000),
                 values=cms.vdouble(0.),
                 uncertainties=cms.vdouble(0.00050)),
        cms.PSet(lowBounds=cms.vdouble(1.000, -999.000),
                 upBounds=cms.vdouble(1.500, 0.940),
Пример #14
0
    bins = cms.VPSet(
            # smearings for rereco 2015, photons,
            # see Fasanella et al., ECAL DPG 17/12/2015, https://indico.cern.ch/event/402608/contribution/1/attachments/1206871/1758949/DPG_reReco_corrections.pdf
        cms.PSet(lowBounds=cms.vdouble(0.000,-999.000), upBounds=cms.vdouble(1.000,0.940),   values=cms.vdouble(0.0087), uncertainties=cms.vdouble(0.0004)),
        cms.PSet(lowBounds=cms.vdouble(0.000,0.940),    upBounds=cms.vdouble(1.000,999.000), values=cms.vdouble(0.0075), uncertainties=cms.vdouble(0.0003)),
        cms.PSet(lowBounds=cms.vdouble(1.000,-999.000), upBounds=cms.vdouble(1.500,0.940),   values=cms.vdouble(0.0145), uncertainties=cms.vdouble(0.0005)),
        cms.PSet(lowBounds=cms.vdouble(1.000,0.940),    upBounds=cms.vdouble(1.500,999.000), values=cms.vdouble(0.0123), uncertainties=cms.vdouble(0.0014)),
        cms.PSet(lowBounds=cms.vdouble(1.500,-999.000), upBounds=cms.vdouble(2.000,0.940),   values=cms.vdouble(0.0184), uncertainties=cms.vdouble(0.0007)),
        cms.PSet(lowBounds=cms.vdouble(1.500,0.940),    upBounds=cms.vdouble(2.000,999.000), values=cms.vdouble(0.0171), uncertainties=cms.vdouble(0.0012)),
        cms.PSet(lowBounds=cms.vdouble(2.000,-999.000), upBounds=cms.vdouble(6.000,0.940),   values=cms.vdouble(0.0211), uncertainties=cms.vdouble(0.0007)),
        cms.PSet(lowBounds=cms.vdouble(2.000,0.940),    upBounds=cms.vdouble(6.000,999.000), values=cms.vdouble(0.0189), uncertainties=cms.vdouble(0.0006)),
        ))

from flashgg.Systematics.ecalElfTools import getRunDependentScaleBins

photonScaleBinsData = getRunDependentScaleBins('flashgg/Systematics/data/Golden10June_plus_DCS.dat')

photonScaleUncertBins = cms.PSet(
    variables = cms.vstring("abs(superCluster.eta)","r9"),
    bins = cms.VPSet(
        ## only uncertainties here. scales loaded by SystematicsCustomize according to process.photonScaleBinsData
            cms.PSet( lowBounds = cms.vdouble(0.000 , -999.000 ) , upBounds = cms.vdouble(1.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00050 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(0.000 , 0.940    ) , upBounds = cms.vdouble(1.000 , 999.000 ) , values = cms.vdouble( 0.  ) , uncertainties = cms.vdouble( 0.00050 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.000 , -999.000 ) , upBounds = cms.vdouble(1.500 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00120 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.000 , 0.940    ) , upBounds = cms.vdouble(1.500 , 999.000 ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00060 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.500 , -999.000 ) , upBounds = cms.vdouble(2.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00200 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(1.500 , 0.940    ) , upBounds = cms.vdouble(2.000 , 999.000 ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00300 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(2.000 , -999.000 ) , upBounds = cms.vdouble(6.000 , 0.940   ) , values = cms.vdouble( 0. ) ,  uncertainties = cms.vdouble( 0.00200 )   ) ,
            cms.PSet( lowBounds = cms.vdouble(2.000 , 0.940    ) , upBounds = cms.vdouble(6.000 , 999.000 ) , values = cms.vdouble( 0.  ) , uncertainties = cms.vdouble( 0.00300 )   ) ,
                    ))