From b071d44faccf05869d864e2b21158d4773fc0641 Mon Sep 17 00:00:00 2001 From: Maurizio Martinelli <maurizio.martinelli@cern.ch> Date: Thu, 5 Jan 2017 10:58:52 +0100 Subject: [PATCH 1/4] Moved some RICH information variables to default since they are useful to determine how the PID info is obtained --- Phys/DecayTreeTuple/src/TupleToolPid.cpp | 34 +++++++++++++----------- 1 file changed, 18 insertions(+), 16 deletions(-) diff --git a/Phys/DecayTreeTuple/src/TupleToolPid.cpp b/Phys/DecayTreeTuple/src/TupleToolPid.cpp index 1af69de3a..67aa67967 100755 --- a/Phys/DecayTreeTuple/src/TupleToolPid.cpp +++ b/Phys/DecayTreeTuple/src/TupleToolPid.cpp @@ -81,6 +81,24 @@ StatusCode TupleToolPid::fill( const Particle* // RICH const RichPID * richPID = proto->richPID(); test &= tuple->column( prefix+"_hasRich", richPID != NULL ); + // Moved some RICH information variables to default since they are useful to determine + // how the PID information has been obtained. + test &= tuple->column( prefix+"_UsedRichAerogel", + richPID ? richPID->usedAerogel() : false ); + test &= tuple->column( prefix+"_UsedRich1Gas", + richPID ? richPID->usedRich1Gas() : false ); + test &= tuple->column( prefix+"_UsedRich2Gas", + richPID ? richPID->usedRich2Gas() : false ); + test &= tuple->column( prefix+"_RichAboveElThres", + richPID ? richPID->electronHypoAboveThres() : false ); + test &= tuple->column( prefix+"_RichAboveMuThres", + richPID ? richPID->muonHypoAboveThres() : false ); + test &= tuple->column( prefix+"_RichAbovePiThres", + richPID ? richPID->pionHypoAboveThres() : false ); + test &= tuple->column( prefix+"_RichAboveKaThres", + richPID ? richPID->kaonHypoAboveThres() : false ); + test &= tuple->column( prefix+"_RichAbovePrThres", + richPID ? richPID->protonHypoAboveThres() : false ); // CALO test &= tuple->column( prefix+"_hasCalo", !proto->calo().empty() ); @@ -90,22 +108,6 @@ StatusCode TupleToolPid::fill( const Particle* { // RICH variables - test &= tuple->column( prefix+"_UsedRichAerogel", - richPID ? richPID->usedAerogel() : false ); - test &= tuple->column( prefix+"_UsedRich1Gas", - richPID ? richPID->usedRich1Gas() : false ); - test &= tuple->column( prefix+"_UsedRich2Gas", - richPID ? richPID->usedRich2Gas() : false ); - test &= tuple->column( prefix+"_RichAboveElThres", - richPID ? richPID->electronHypoAboveThres() : false ); - test &= tuple->column( prefix+"_RichAboveMuThres", - richPID ? richPID->muonHypoAboveThres() : false ); - test &= tuple->column( prefix+"_RichAbovePiThres", - richPID ? richPID->pionHypoAboveThres() : false ); - test &= tuple->column( prefix+"_RichAboveKaThres", - richPID ? richPID->kaonHypoAboveThres() : false ); - test &= tuple->column( prefix+"_RichAbovePrThres", - richPID ? richPID->protonHypoAboveThres() : false ); test &= tuple->column( prefix+"_RichDLLe", proto->info(LHCb::ProtoParticle::RichDLLe,-1000) ); test &= tuple->column( prefix+"_RichDLLmu", -- GitLab From c921e118bfd575beb8a910867b63486d1a1f196c Mon Sep 17 00:00:00 2001 From: ibelyaev <Ivan.Belyaev@cern.ch> Date: Wed, 11 Jan 2017 14:38:09 +0100 Subject: [PATCH 2/4] Phys/DaVinciTrackScaling: add scripts for momentum scaling for 2015&2016 --- Phys/DaVinciTrackScaling/doc/release.notes | 3 + .../DaVinciTrackScaling/momentum_scale.py | 215 +++++++----------- .../DaVinciTrackScaling/scaling_2011.py | 103 +++++++++ .../DaVinciTrackScaling/scaling_2015.py | 105 +++++++++ .../DaVinciTrackScaling/scaling_2016.py | 146 ++++++++++++ 5 files changed, 435 insertions(+), 137 deletions(-) create mode 100644 Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/scaling_2011.py create mode 100644 Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/scaling_2015.py create mode 100644 Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/scaling_2016.py diff --git a/Phys/DaVinciTrackScaling/doc/release.notes b/Phys/DaVinciTrackScaling/doc/release.notes index ca787c92b..fb0c96053 100644 --- a/Phys/DaVinciTrackScaling/doc/release.notes +++ b/Phys/DaVinciTrackScaling/doc/release.notes @@ -4,6 +4,9 @@ ! Purpose : Track scaling and smearing !----------------------------------------------------------------------------- +! 2017-01-11 - Vanya Belyaev + - prepare scripts for momenentum scaling for -2015&2016 + !=================== DaVinciTrackScaling v1r6 2016-02-04 ===================== ! 2015-12-05 - Gerhard Raven diff --git a/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/momentum_scale.py b/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/momentum_scale.py index 0135bcfaa..72722cc65 100644 --- a/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/momentum_scale.py +++ b/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/momentum_scale.py @@ -1,8 +1,6 @@ #!/usr/bin/env python # ============================================================================= -# $Id:$ -# ============================================================================= -# @file DaVinciTools/momentum_scale.py +# @file momentum_scale.py # # Helper script to convert input data (histograms + offsets) to XML-format # @@ -21,10 +19,6 @@ # @date 2013-04-20 # @see TrackScaleState # -# $Revision:$ -# Last modification $Date:$ -# by $Author:$ -# # ============================================================================ """ Helper script to convert input data (histograms + offsets) to XML-format @@ -46,76 +40,24 @@ Input data are specified close to the start of the script __author__ = "Vanya BELYAEV Ivan.Belyaev@itep.ru" __date__ = "2013-04-20" __version__ = "$Revision:$" -# ============================================================================ -print 100*'*' -print __doc__ -print 100*'*' +__all__ = ( + 'logger' , + 'buildXML' , + 'offsets_to_histo' + ) # ============================================================================ import ROOT from Ostap.PyRoUts import h1_axis,VE - - -# ============================================================================= -## get the data -# ============================================================================= -# -year = '2k+11' -reco = 'Reco14' -# -directory = '/afs/cern.ch/user/m/mneedham/public/2011-repro-Calib/' -filename = 'BFitResults-2011repro.root' -the_file = ROOT.TFile( directory + filename ) -print 100*'*' -print 'HISTOS are read from ' , directory + filename -print the_file.ls() -print 100*'*' - -delta = 1.0e-4 -print 'global DELTA is defined as %s ' % delta - -# -# get histograms -# -idp_plus = the_file.Get('idp-plus' ) -idp_minus = the_file.Get('idp-minus') - -# ============================================================================= -## CalScale.h: -## if (run >= 87219 && run <= 89777 ) alpha =0.0; -## else if (run >= 89333 && run <= 90256 ) alpha = 6.36699e-05 ; -## else if (run >= 90257 && run <= 90763 ) alpha = 7.11719e-05 ; -## else if (run >= 91556 && run <= 93282) alpha = 1.53674e-05 ; -## else if (run >= 93398 && run <= 94386 ) alpha = 0.000144135 ; -## else if (run >= 95948 && run <= 97028 ) alpha = 0.000214408 ; -## else if (run >= 97114 && run <= 98882 ) alpha = 3.41493e-05 ; -## else if (run >= 98900 && run <= 101862 ) alpha = 0.000137622; -## else if (run >= 101891 && run <= 102452 ) alpha = 3.73981e-05 ; -## else if (run >= 102499 && run <= 102907 ) alpha = 0.000169023 ; -## else if (run >= 103049 && run <= 103687 ) alpha = 3.20303e-05 ; -## else if (run >= 103954 && run <= 104414 ) alpha = 0.00017584 ; - -offsets = [ - ## start stop offset - ( 87219 , 89777 , 0.0 ) , - ( 89333 , 90256 , 6.36699e-05 ) , - ( 90257 , 90763 , 7.11719e-05 ) , - ( 91556 , 93282 , 1.53674e-05 ) , - ( 93398 , 94386 , 0.000144135 ) , - ( 95948 , 97028 , 0.000214408 ) , - ( 97114 , 98882 , 3.41493e-05 ) , - ( 98900 , 101862 , 0.000137622 ) , - ( 101891 , 102452 , 3.73981e-05 ) , - ( 102499 , 102907 , 0.000169023 ) , - ( 103049 , 103687 , 3.20303e-05 ) , - ( 103954 , 104414 , 0.00017584 ) , - ] - +# ============================================================================ +from Ostap.Logger import getLogger +logger = getLogger('MomentumScale') +# ============================================================================ # ============================================================================ ## convert run offsets into the histogram -# @param offset the list of run-dependes offsets: [ (start,stop,offset) ] -# @param title (optional) the histogram title -# @code +# @param offset the list of run-dependes offsets: [ (start,stop,offset) ] +# @param title (optional) the histogram title +# @code # offsets_2011 = [ # ## start stop offset # ( 89333 , 90256 , 9.97629e-05 ) , @@ -129,8 +71,7 @@ offsets = [ # @date 2013-04-20 # def offsets_to_histo ( offsets , title = '' ) : - """ - Helper function to convert the offset list into histogram + """Helper function to convert the offset list into histogram offsets_2011 = [ ## start stop offset @@ -150,7 +91,7 @@ def offsets_to_histo ( offsets , title = '' ) : ticks.add ( i[0] ) ticks.add ( i[1] + 1 ) ## take care about equality - ## create the histogram with appropriate biunings + ## create the histogram with appropriate binnings h = h1_axis ( list ( ticks ) , double = True , title = title ) for i in offsets : @@ -161,72 +102,72 @@ def offsets_to_histo ( offsets , title = '' ) : return h +line1 = '<?xml version="1.0" encoding="ISO-8859-1"?>' +line2 = '<!DOCTYPE DDDB SYSTEM "../../../DTD/structure.dtd">' # ============================================================================= # Build XML-document # ============================================================================= - -## -line1 = '<?xml version="1.0" encoding="ISO-8859-1"?>' -line2 = '<!DOCTYPE DDDB SYSTEM "../../../DTD/structure.dtd">' -## - -## import lxml.etree -import lxml.builder - -E = lxml.builder.ElementMaker() -DDDB = E.DDDB -CONDITION = E.condition -PARAM = E.param - -xml_delta = PARAM ( - str ( delta ), - name = 'Delta' , - type = 'double' , - comment = 'Global delta scale %s %s' % ( reco , year ) - ) - -xml_idplus = PARAM ( - idp_plus.toString() , - name = "IdpPlus" , - type = "Histo2D" , - comment = "Calibration for IdpPlus %s %s" % ( reco , year ) - ) - -xml_idpminus = PARAM ( - idp_minus.toString () , - name = "IdpMinus" , - type = "Histo2D" , - comment = "Calibration for IdpMinus %s %s" % ( reco , year ) - ) - -xml_offsets = PARAM ( - offsets_to_histo ( offsets , 'Run-dependent offsets' ).toString() , - name = "Offsets" , - type = "Histo1D" , - comment = "Run-dependent offsets %s %s" % ( reco , year ) - ) - -document = DDDB( - CONDITION ( - xml_delta , - xml_idplus , - xml_idpminus , - xml_offsets , - ## comment = "Momentum scale %s %s" % ( reco , year ) , - name = 'MomentumScale' - ), - ) - - -xml_name = 'MomentumScale.xml' -xml_file = open ( xml_name ,'w') - -print >> xml_file, line1+'\n'+line2+'\n\n\n'+lxml.etree.tostring( document , pretty_print=True) - -xml_file.close() -print 100*'*' -print 'XML file "%s" is generated in CWD ' % xml_name -print 100*'*' +def buildXML ( reco , + year , + idp_plus , + idp_minus , + offsets , + delta ) : + + import lxml.builder + + E = lxml.builder.ElementMaker() + DDDB = E.DDDB + CONDITION = E.condition + PARAM = E.param + + xml_delta = PARAM ( + str ( delta ), + name = 'Delta' , + type = 'double' , + comment = 'Global delta scale %s %s' % ( reco , year ) + ) + + xml_idplus = PARAM ( + idp_plus.toString() , + name = "IdpPlus" , + type = "Histo2D" , + comment = "Calibration for IdpPlus %s %s" % ( reco , year ) + ) + + xml_idpminus = PARAM ( + idp_minus.toString () , + name = "IdpMinus" , + type = "Histo2D" , + comment = "Calibration for IdpMinus %s %s" % ( reco , year ) + ) + + xml_offsets = PARAM ( + offsets_to_histo ( offsets , 'Run-dependent offsets' ).toString() , + name = "Offsets" , + type = "Histo1D" , + comment = "Run-dependent offsets %s %s" % ( reco , year ) + ) + + document = DDDB( + CONDITION ( + xml_delta , + xml_idplus , + xml_idpminus , + xml_offsets , + ## comment = "Momentum scale %s %s" % ( reco , year ) , + name = 'MomentumScale' + ), + ) + + xml_name = 'MomentumScale.xml' + xml_file = open ( xml_name ,'w') + + print >> xml_file, line1+'\n'+line2+'\n\n\n'+lxml.etree.tostring( document , pretty_print=True) + + xml_file.close() + logger.info ( 'XML file "%s" is generated in CWD ' % xml_name ) + # ============================================================================= diff --git a/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/scaling_2011.py b/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/scaling_2011.py new file mode 100644 index 000000000..cb69d43b8 --- /dev/null +++ b/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/scaling_2011.py @@ -0,0 +1,103 @@ +#!/usr/bin/env python +# ============================================================================= +# @file scaling_2011.py +# +# Momentum scaling for 2011 +# +# Input data are specified close to the start of the script: +# +# - input ROOT file with historams +# the file should have two 2D-histograms, named as 'idp-plus' and 'idp-minus' +# +# - global delta +# - list of run-dependent offsets as list of triplets: [ (start,stop,offset) ], +# where (start,stop) represent the run range (both edges are inclusive) +# +# As the output xml-file MomentumScale.xml is generated in cwd +# +# @author Matt NEEDHAM +# @author Vanya BELYAEV Ivan.Belyaev@itep.ru +# @date 2013-04-20 +# @see TrackScaleState +# ============================================================================ +"""Momentum sclaing for 2011 +- input ROOT file with historams +the file should have two 2D-histograms, named as 'idp-plus' and 'idp-minus' + +- global delta +- list of run-dependent offsets as list of triplets: [ (start,stop,offset) ], +where (start,stop) represent the run range (both edges are inclusive) + +As the output xml-file MomentumScale.xml is generated in cwd + +""" +# ============================================================================ +__author__ = "Vanya BELYAEV Ivan.Belyaev@itep.ru" +__date__ = "2013-04-20" +__version__ = "$Revision:$" +# ============================================================================ +import ROOT, Ostap.PyRoUts +# ============================================================================ +from momentum_scale import buildXML, logger +# ============================================================================ +logger.info(100*'*') +logger.info(__doc__) +logger.info(100*'*') +# ============================================================================ + +# ============================================================================ +## the histograms +with ROOT.TFile( 'BFitResults-2011repro.root','READ') as the_file : + idp_plus = the_file.Get('idp-plus' ).clone() + idp_minus = the_file.Get('idp-minus').clone() + logger.info ( 'HISTOS are read from %s' % the_file ) + +# ============================================================================= +## global delta +delta = 1.0e-4 +logger.info ( 'Global delta: %s' % delta) + +# ============================================================================= +## CalScale.h: +## if (run >= 87219 && run <= 89777 ) alpha =0.0; +## else if (run >= 89333 && run <= 90256 ) alpha = 6.36699e-05 ; +## else if (run >= 90257 && run <= 90763 ) alpha = 7.11719e-05 ; +## else if (run >= 91556 && run <= 93282) alpha = 1.53674e-05 ; +## else if (run >= 93398 && run <= 94386 ) alpha = 0.000144135 ; +## else if (run >= 95948 && run <= 97028 ) alpha = 0.000214408 ; +## else if (run >= 97114 && run <= 98882 ) alpha = 3.41493e-05 ; +## else if (run >= 98900 && run <= 101862 ) alpha = 0.000137622; +## else if (run >= 101891 && run <= 102452 ) alpha = 3.73981e-05 ; +## else if (run >= 102499 && run <= 102907 ) alpha = 0.000169023 ; +## else if (run >= 103049 && run <= 103687 ) alpha = 3.20303e-05 ; +## else if (run >= 103954 && run <= 104414 ) alpha = 0.00017584 ; + +offsets = [ + ## start stop offset + ( 87219 , 89777 , 0.0 ) , + ( 89333 , 90256 , 6.36699e-05 ) , + ( 90257 , 90763 , 7.11719e-05 ) , + ( 91556 , 93282 , 1.53674e-05 ) , + ( 93398 , 94386 , 0.000144135 ) , + ( 95948 , 97028 , 0.000214408 ) , + ( 97114 , 98882 , 3.41493e-05 ) , + ( 98900 , 101862 , 0.000137622 ) , + ( 101891 , 102452 , 3.73981e-05 ) , + ( 102499 , 102907 , 0.000169023 ) , + ( 103049 , 103687 , 3.20303e-05 ) , + ( 103954 , 104414 , 0.00017584 ) , + ] + +# ============================================================================= +# Build XML-document +# ============================================================================= +buildXML ( year = '2k+11' , + reco = 'Reco14' , + idp_plus = idp_plus , + idp_minus = idp_minus , + offsets = offsets , + delta = delta ) + +# ============================================================================= +# The END +# ============================================================================= diff --git a/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/scaling_2015.py b/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/scaling_2015.py new file mode 100644 index 000000000..9168fb4c2 --- /dev/null +++ b/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/scaling_2015.py @@ -0,0 +1,105 @@ +#!/usr/bin/env python +# ============================================================================= +# @file scaling_2015.py +# +# Momentum scaling for 2015 +# +# @author Matt NEEDHAM +# @author Vanya BELYAEV Ivan.Belyaev@itep.ru +# @date 2016-01-11 +# ============================================================================ +"""Momentum scaling fro 2015 + +Input data are specified close to the start of the script + +- input ROOT file with historams +the file should have two 2D-histograms, named as 'idp-plus' and 'idp-minus' + +- global delta +- list of run-dependent offsets as list of triplets: [ (start,stop,offset) ], +where (start,stop) represent the run range (both edges are inclusive) + +As the output xml-file MomentumScale.xml is generated in cwd + +""" +# ============================================================================ +__author__ = "Vanya BELYAEV Ivan.Belyaev@itep.ru" +__date__ = "2016-01-11" +__version__ = "" +# ============================================================================ +import ROOT, Ostap.PyRoUts +# ============================================================================ +from momentum_scale import buildXML, logger +# ============================================================================ +logger.info(100*'*') +logger.info(__doc__) +logger.info(100*'*') +# ============================================================================ + +# ============================================================================ +## the histograms +with ROOT.TFile('BFitResults-fine-2015.root','READ') as the_file : + idp_plus = the_file.Get('idp-plus' ).clone() + idp_minus = the_file.Get('idp-minus').clone() + logger.info ( 'HISTOS are read from %s' % the_file ) + +# ============================================================================= +## global delta +delta = -0.165063e-3 +logger.info ( 'Global delta: %s' % delta) + +## //down +## if (run >= 162200 && run <= 163200 ) alpha = 0.0; +## //gap +## else if (run >= 163250 && run <= 164000 ) alpha = 0.0406160e-3 ; +## else if (run >= 164001 && run <= 164500 ) alpha = -0.0405840e-3; +## else if (run >= 164501 && run <= 164900 ) alpha = -0.0751180e-3; +## else if (run >= 164901 && run <= 165200 ) alpha = -0.0706750e-3; +## //gap +## else if (run >= 165201 && run <= 165400 ) alpha = -0.0708270e-3 ; +## else if (run >= 165401 && run <= 165600 ) alpha = -0.0670200e-3 ; +## // up +## else if (run >= 166250 && run <= 166400 ) alpha = 0.0940130e-3 ; +## else if (run >= 166401 && run <= 166550 ) alpha = 0.0539590e-3 ; +## // gap +## else if (run >= 166600 && run <= 166700 ) alpha = 0.0157780e-3 ; +## else if (run >= 166700 && run <= 166800 ) alpha = 0.0421130e-3 ; +## // gap +## else if (run >= 166850 && run <= 167000 ) alpha = -0.0553680e-3 ; +## else if (run >= 167000 && run <= 167150 ) alpha = -0.0373790e-3; + +## Run-dependent offsets +offsets = [ + ## start stop offset + ( 163250 , 164000 , 0.0406160e-3 ) , + ( 164001 , 164500 , -0.0405840e-3 ) , + ( 164501 , 164900 , -0.0751180e-3 ) , + ( 164901 , 165200 , -0.0706750e-3 ) , + ## + ( 165201 , 165400 , -0.0708270e-3 ) , + ( 165401 , 165600 , -0.0670200e-3 ) , + ## + ( 166250 , 166400 , 0.0940130e-3 ) , + ( 166401 , 166550 , 0.0539590e-3 ) , + ## + ( 166600 , 166700 , 0.0157780e-3 ) , + ( 166700 , 166800 , 0.0421130e-3 ) , + ## + ( 166850 , 167000 , -0.0553680e-3 ) , + ( 167000 , 167150 , -0.0373790e-3 ) + ] + +# ============================================================================= +# Build XML-document +# ============================================================================= +buildXML ( year = '2k+15' , + reco = 'Reco16' , + idp_plus = idp_plus , + idp_minus = idp_minus , + offsets = offsets , + delta = delta ) + +# ============================================================================ +# The END +# ============================================================================ + diff --git a/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/scaling_2016.py b/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/scaling_2016.py new file mode 100644 index 000000000..2a16ef462 --- /dev/null +++ b/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/scaling_2016.py @@ -0,0 +1,146 @@ +#!/usr/bin/env python +# ============================================================================= +# @file scaling_2016.py +# +# Momentum scaling for 2016 +# +# @author Matt NEEDHAM +# @author Vanya BELYAEV Ivan.Belyaev@itep.ru +# @date 2016-01-11 +# ============================================================================ +"""Momentum scaling fro 2016 + +Input data are specified close to the start of the script + +- input ROOT file with historams +the file should have two 2D-histograms, named as 'idp-plus' and 'idp-minus' + +- global delta +- list of run-dependent offsets as list of triplets: [ (start,stop,offset) ], +where (start,stop) represent the run range (both edges are inclusive) + +As the output xml-file MomentumScale.xml is generated in cwd + +""" +# ============================================================================ +__author__ = "Vanya BELYAEV Ivan.Belyaev@itep.ru" +__date__ = "2016-01-11" +__version__ = "" +# ============================================================================ +import ROOT, Ostap.PyRoUts +# ============================================================================ +from momentum_scale import buildXML, logger +# ============================================================================ +logger.info(100*'*') +logger.info(__doc__) +logger.info(100*'*') +# ============================================================================ + +# ============================================================================ +## the histograms +with ROOT.TFile('BFitResults-fine-2016.root','READ') as the_file : + idp_plus = the_file.Get('idp-plus' ).clone() + idp_minus = the_file.Get('idp-minus').clone() + logger.info ( 'HISTOS are read from %s' % the_file ) + +# ============================================================================= +## global delta +delta = -0.19e-3 +logger.info ( 'Global delta: %s' % delta) + +## //down +## if (run >=174500 && run <= 175100 ) alpha = 0; +## else if (run >= 175101 && run <= 175900 ) alpha =0.0302171e-3 ; +## else if (run >= 175901 && run <= 176100 ) alpha =0.0558559e-3 ; +## else if (run >= 176101 && run <= 176330 ) alpha = 0.0874364e-3; +## //up +## else if (run >= 176657 && run <= 177000) alpha = 0.130490e-3; +## else if (run >= 177001 && run <= 177250) alpha = 0.146611e-3; +## else if (run >= 177251 && run <= 177350) alpha = 0.166163e-3; +## else if (run >= 177351 && run <= 177500) alpha = 0.153974e-3; +## // +## else if (run >= 177501 && run <= 177600) alpha = 0.239742e-3; +## else if (run >= 177601 && run <= 177800) alpha = 0.299172e-3; +## // +## else if (run >= 178001 && run <= 178250) alpha = 0.238894e-3; +## else if (run >= 178251 && run <= 178500) alpha = 0.221683e-3; +## else if (run >= 178501 && run <= 179001) alpha = 0.232447e-3; +## else if (run >= 179001 && run <= 179500) alpha = 0.251556e-3; +## else if (run >= 179501 && run <= 179999) alpha = 0.246533e-3; +## else if (run >= 180000 && run <= 180300) alpha = 0.270993e-3; +## else if (run >= 180301 && run <= 180640) alpha = 0.255366e-3; +## // down +## else if (run >= 181200 && run <= 181700 ) alpha = 0.0686362e-3 ; +## else if (run >= 181700 && run <= 182300 ) alpha =0.0945660e-3; +## else if (run >= 182301 && run <= 182500 ) alpha =0.0683000e-3 ; +## else if (run >= 182501 && run <= 182800 ) alpha = 0.0877048e-3; +## else if (run >= 182801 && run <= 183000 ) alpha = 0.126342e-3; +## else if (run >= 183001 && run <= 183200 ) alpha = 0.119248e-3; +## //up +## else if (run >= 183800 && run <= 184200 ) alpha =0.222877e-3 ; +## else if (run >= 184201 && run <= 184657 ) alpha = 0.211130e-3; +## // down +## else if (run >= 184700 && run <= 185000 ) alpha =0.2111e-3 ; +## else if (run >= 185000 && run <= 185100 ) alpha =0.05129e-3 ; +## else if (run >= 185101 && run <= 185200 ) alpha = -0.00853040e-3 ; +## else if (run >= 185201 && run <= 185300 ) alpha =0.0154273e-3 ; +## else if (run >= 185301 && run <= 185400 ) alpha =0.0351537e-3 ; +## else if (run >= 185401 && run <= 185600 ) alpha = 0.0677423e-3 ; + + +## Run-dependent offsets +offsets = [ + ## start stop offset + ## //down + ( 175101 , 175900 , 0.0302171e-3 ) , + ( 175901 , 176100 , 0.0558559e-3 ) , + ( 176101 , 176330 , 0.0874364e-3 ) , + ## //up + ( 176657 , 177000 , 0.130490e-3 ) , + ( 177001 , 177250 , 0.146611e-3 ) , + ( 177251 , 177350 , 0.166163e-3 ) , + ( 177351 , 177500 , 0.153974e-3 ) , + ## // + ( 177501 , 177600 , 0.239742e-3 ) , + ( 177601 , 177800 , 0.299172e-3 ) , + ## // + ( 178001 , 178250 , 0.238894e-3 ) , + ( 178251 , 178500 , 0.221683e-3 ) , + ( 178501 , 179001 , 0.232447e-3 ) , + ( 179001 , 179500 , 0.251556e-3 ) , + ( 179501 , 179999 , 0.246533e-3 ) , + ( 180000 , 180300 , 0.270993e-3 ) , + ( 180301 , 180640 , 0.255366e-3 ) , + ## // down + ( 181200 , 181700 , 0.0686362e-3 ) , + ( 181700 , 182300 , 0.0945660e-3 ) , + ( 182301 , 182500 , 0.0683000e-3 ) , + ( 182501 , 182800 , 0.0877048e-3 ) , + ( 182801 , 183000 , 0.126342e-3 ) , + ( 183001 , 183200 , 0.119248e-3 ) , + ## //up + ( 183800 , 184200 , 0.222877e-3 ) , + ( 184201 , 184657 , 0.211130e-3 ) , + ## // down + ( 184700 , 185000 , 0.2111e-3 ) , + ( 185000 , 185100 , 0.05129e-3 ) , + ( 185101 , 185200 , -0.00853040e-3 ) , + ( 185201 , 185300 , 0.0154273e-3 ) , + ( 185301 , 185400 , 0.0351537e-3 ) , + ( 185401 , 185600 , 0.0677423e-3 ) + ] + +# ============================================================================= +# Build XML-document +# ============================================================================= +buildXML ( year = '2k+16' , + reco = 'Reco16' , + idp_plus = idp_plus , + idp_minus = idp_minus , + offsets = offsets , + delta = delta ) + +# ============================================================================ +# The END +# ============================================================================ + -- GitLab From 39124fda6ee74eaab4cfc8be5eb271ad9cc4aa20 Mon Sep 17 00:00:00 2001 From: ibelyaev <Ivan.Belyaev@cern.ch> Date: Wed, 11 Jan 2017 14:49:47 +0100 Subject: [PATCH 3/4] Phys/DaVinciTrackScaling: - add original ROOT files for momentum scaling --- .../BFitResults-2011repro.root | Bin 0 -> 11567 bytes .../BFitResults-fine-2015.root | Bin 0 -> 7771 bytes .../BFitResults-fine-2016.root | Bin 0 -> 5219 bytes 3 files changed, 0 insertions(+), 0 deletions(-) create mode 100644 Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/BFitResults-2011repro.root create mode 100644 Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/BFitResults-fine-2015.root create mode 100644 Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/BFitResults-fine-2016.root diff --git a/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/BFitResults-2011repro.root b/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/BFitResults-2011repro.root new file mode 100644 index 0000000000000000000000000000000000000000..e5f5804570bb46ffce4c267efa85dd6d3e39f35f GIT binary patch literal 11567 zcmbW7bC6}-(x*?kyKLLGZKJzvv&*)PF59+kS9RH4wr%sP?mP3|iTP*d-V?E6XT;gL zbDecoJUiDf)6Ux30RR|x0RRAo000dg06;zeX{&t>pwHn2`Y!{+0{}o+0RV9rlMSW6 zBWx$3Lw!Ro@4Wc@_W$eC3i$INbBQlG5CFNqPyUPs0Dviqnpqg(2#T6H$Q#)^S~%F# z7@Juc(KtBJI2hSG(As^@4&K_)!N$=6UtUg55#RVfza0SfC+gqG03h$rpkANjZ*%~F zO6C9clR%mO%pB7HuRq<NnFHbe@ds5DXB38Tw4t%FaJ2v5S9R$xuOI+_NPr8_0IH7Z zZ=>5)?^!GG3gcURPI`S_bYE}s({qpf1~9@D_pgn?B55!*Ne`GO>dm=O?x74Zw-ZT= z#L#-8LaW~%hl7w=Mr(G#{anxnjwvlIoi%Ff?I)Z&vz!;qw}hN%AGkK2J`UQrOu%y? z@C|i`NxsawB91_JLw=bG#s{r-819^JsAzi_z#)B-kBNg_8&zexrCQ5{H{@*u5y}8U z4}>`MP^=FXS;wE5Uz|sT(B}iv^XZmBP_{lE7d9wknm)J|$~g78ATA*dH(Foh1X7f4 z55}IvYGkygAD=qAqk3zZ@tavQB=V8W8~J``J9V7*Ock;x4Tj}Yi73+HUL_LP5k5~x zNlLlxQi)y$v`OZyu@DN1<&%ggL2|)D;c9mu61V+8D`wg|w=SaC&%2;RhaC*nkFuWL zIVUZ`s4=UuB{MU-pfc{(!mH=P*K}>2`@3ys?6Q?DRjLAe(QBM*;K!4_OJf*J<{*Rd z0q13@3*0sx18e!9`8I2)LO##jk`I&M#V-!mqdXV?a(OC<E4s*;GzR?Z2kCD^<gR5@ zLl2$WSIOgGd~dsBF;#<n{Yay71KuazAbec6`mg?DMtEHBXskS07q^jdH5_Z_H*hZv zG^G#;QFvV!zi}V$H*p+-pO2Osk9vc4s<v~)2u(mb+^dr}I&8U{^l#L4-mS0FyS5+k z%-&?rw>dgT>r^XmbGJkI7lJe2bk?~)jQuy{#M->Lcor76c{x3ciY^qeFFL2efp6e* z1SGm=6MTd0bYD>g)<%C`FrSGtle(;ecOE?99tLLv9s0m$OpdmL*X(j=>6h-XEu~(7 zH%{&hIJ6!N?`tK84$2E<jI*;cUiaO$n%ENZIRtQT-lcZ_!uZhc3WDE~p1Ta3p?HKO zx?Wzl)L<3cPEQ7MSv020HF7wJf!nL8AVMQwtDdWAAj#bP$)KV6BRaXQLee0&<fkL1 zPYU^got<+V3)qX4!5}VITu>}OttOjaK74PcbUSyF`Z&UUz0eZ#Rt9|j4GD(|O_<_) zuSRGGjm*`E>fW;nC5CuSwJbSUzP;F(q|&b0t*%R-!ijSb+l)Z`oJipNPCn|%hh<nw zwDdkQZ3O+4`?5q-Q#@bUEe!5q{lmejL$6sI(50L@)UXAOc2X!`L>y<2qk^fmM*!b7 zKb{Mi3G%^dP!3)}Iix0h{a}%+x8h##+c<r}2gu?$y)JOA@A_c-fERwpX7QT1*J}BU z&dTGu-0gy3q1EhG{9s!2fD^Qdk;1?w+`9KLwQR|our0AvTJ(OACxh-8v}p&FGF*|R z=9u}%zS8HvT7Kb|<l<YXoDliEVgCV6%$qn40D%7Z%Z#AE%wYJ#3^4hBnDOP48J1>N zj`sf;pVXiS1qA)2hO)awoDwphPM6adCsL686&%P<bV7ny^knhsbSlUdU_>-yKM5iE zF&}*&iSQBy8*WIZX=)-+a?v~=U;$L3s;iI<bSU^$alQT@-2s*i++mE}m1D^JEmcv= zjts9%a<lH|uB#gy4<0&`6d(_n^cc`DVCnVm09R;TV)R=AMxg5(6~(S<;%Wzxsr6{F z&P#l{<pW!t8_u)U%}uv*Fxdp?O>S<{lepeA#Q=!PYCh~Qz$IBAWIk8?=tZ(i2+r<e zzt*aVHI6v|8*xUI8RXTRq+bLFBpw3h)yN9QIT82j$Cs;ic_dM(TYpY`5hTg0&Aa|Z zG=qntElTaFXahbKy(5pueO!>PEWMq(P25s1t0U?LSRCfL{Y~0~!H+K{&k~gy?{bKZ z;0}6@AMW}y0Y-mBo{jA>K=cOMKE{z30b@zgc9xUDY`;XSsk8qMYmftSB&~+!V}N0y zkI&#&`do`n_{k-sH%qCrzRE^{i9y%P)$xK!+CpUUVA`J8LBb~G=~tb0?}sVq_SwM_ zoVd*S$*v<WD21GzaFsi)!A>R!(n*|m*ut!*$Y`X;>{O7>U%KyDXjT^zOgjp&{mDC1 z4uaah#;B^i7|r%EMj9B2HVD?X(KOtC(R5<|vSu6qQW>e!>q%8n8|A%hI@%`ZT%i|i ze&wEm`a|u-#QtYfU142)&0*f@=q73@y^;~YI+>@S!|-9{F(s_ZX#7MGSACvc?cINH zg}K*pb+KkFvdosba&Y$9k<ptA&vM}^o*RE%B`wObK=F+8sP%~R==sd^Sod=9U8K0B z6zxalTYY(^lv6uF90t1QX0E72e;XCqj}cL4zVCZ!?+<1B0j~t>RzJ}%E(d{Z=zTcp z%`!-zSwiC4<qi&yT^j3dH5CSUB1H05eo?!i?O0jo7kB|V)i-FhbLlZ^OEp@MfWGF{ zUztr<KMf)h;W{&Q<~T5R3T~x#dEB2=wfJ!}fYWOl(@BHXAJE`IcGN8|Q)neM|0?IF z?<W(bd?|&Mr;VX~|E84uG6qs<s82)M99%GrtxH+l^r(t#hp2aG$3E24e|Dv-N645O ztf~y<O0f;qBVgX5D)vUDxOH2Qv*%ly{p%#<u{Icq!&n>H4&ezp<mh{()pit@)S#T^ zp+_icZ2!P}nNs-dhXz)ihKGRynDzZ4O)!&Fsa{&XvEr;+PziMh?gZMg=j06^C<6;{ zD#=1~DJ$!&(8}NiqSko0EpA~wf?OG^vj}d@PN6IEhUz=3LcjK_Oyn1C>NlF_ua~c9 z{6GA?OJeNBy3zDv!PN;FA!p?$rDHnbUXDzJq8IZn-IcT{5+*l}UoDlQux5Di*P}`! zD{^7HLcEXdYct2btrpnK&O~8czmJt1FIAnsc>t3Uj>X$Kwdcfs0F_<EcYl(`>@R6T z|B{C34{5%(|3jLuhSttjfAf=nTc50<0R@Eqi#3+gzm1UjK3)>nEuakMQ)k<BWmJ|E z>6T92aEXM9>|;F4@=$1OlKpJ0o&2>cQAX!e^Dg&hW`t*-dq<e7X7A7qS%8E6ofgDV zMTH1qg5V9;IVV@LMmI*&bq?t>Z_;mGCr0hh_a3^tb-=8^*&Ka>=kqoV06^7;!Tm7w zM#ptVbcquy?N@<Qi;?sS#6dRo4`1uslAOG3zVe9*i1UK5Ln)Yb$m9--VDw;7M8pU9 ztS)T=RmWc2EX}-q^`@uV@-{8jYZJ%6cc_gZEWBz2lmBLZy{7Nc{ypPrhk>Bs70S}P z)Urgs0htJ_;pa0p3vLt@_pt{>gLLs6bT;RBfp+V(l&eeGNVt039t*EMMj)U%4XZDe zs04#wBDIjxs-EJhCOl%1?nbt{nXWnk37ccJs1z#y*J1pQ^g2&4?}=YaJJZXmr}Oz< z`&kte-?+<(_Ik5H8Qsh38@el>hh1OUcWxf4%moz|q?KXJvoE?P*7<M8n|h~SWWmvD zi)~d$b6B%%TO>efj?fTptTPhGbV=T*kNf^0>>iGnqbM(_Hn(RG%HHmWi0^S~(+L_` zp~D?ImQU|Ph~TyB9>$KVzn@oTM^n;rdEXsz?ozX6F7mBxllYob-U+LnuAVwmIBqU| zES`+Y%ed3IJB}nY@XB48>sh9K@y;<jO@!pOF)c4iC&jqrPT$X8PbNk;+o~=Ts@^U$ zE6fYs<<wtBEnZ}_HZ~i-BRWxseNTY<sny|RHkvD!S-n&5x#mVXvCH6^aYuQ}EzLL4 z9D1DMn(CVDidk`>KlytzV#Ab7`z7NcjXBNtj?E>MTmrNWl0d$(^@C*CyWB=<y46Nc zb*jQq*{0nQSrb|ptJ%y`^_pc!yl2)-6(_CLZnFPIw3-FAtgY>lgE+DB$D($G<@orS zKNN01ln_iC1&H=oMVHVaAJl#<Va6p>O^%biS3%#4^wYAKT|2=vI!LpUBVAPCwDvb% zL(m20rf+S$u2C!`In^!<qUkD*F?EUpkbNQuzOIk~_AQ#$MGgje7!Q6&r9BvA`q9ed zp((-P0IPDT0LJIiTH&9y<d2on+Hkb$D#MQV&FUCbl6`w5B02o@)ix6RUw{;>A&&y+ z7v`W%v$Jl;GevIDDj%o|1F2o>PHY<@-Zxyv_8%uhQ+8avb@4pCLl!Sh2HsNZc3jl1 zrAVA4VWZL@qGuPXG7r^nbVxqb_qMvelt%GYnix>Me(yx7tC3omwT7NAetjKsv0T6_ zNHU%HmbaEvl82|639q#~zY=&0{*o40zvv2%0xPD7vmdLM!e07vJl)Ql{jq{R<@kc4 z{mSv3sPD0@wW$twBX8fIX4~-#uk{YyF}I4M*FHs}{xRc}y)-3X-K|IL$mHsX(|j8B z3twDpi|{9CRQ`e{;xA}e{(uIS<R8#L{sqn7q~`xHKT(4L3W)d@YRpY%Rgn2ym40TZ zlI3L?&k=+*(;a_f&OE6@@sQRq^oK5AY-P99inG*mtQ&z06E7)d38U9wyOlzKQ2&0v zQEK(A%#br=JDhSx@~Bh>32CkTVD+QJ?3GRSH9K?0{ha5(In#0?6!5@m^99`jD&2uq zAQ6BL!um{L2Y_~}qJ1fJm~rzyQt_>)$^n-`COFC*aQWrg+R`&&R}Ne#2y}@CopYlV zN{a9Vsv`0|kPkro4oC+y7Rhs%iwgCh;b7TLMzAbYgLK9jqpBHVsgO#2CN=kwi5Mm_ z8Mc)D^&5{O9Z{uChB<GBL|=bnCE%jq%B8p{zR&2JDGY=0TL$HJdj8H$7hSv!51Sn= zFlt3d0*aQb#4A&Zg#vg7<x}dgP&@^jL_P99d-iM{!avb{Q|s^g{hTj+!+a-B^nEFD zf!!_-0V-4N#q+v#T0y<Z{^%@akZ;SK%EUU$*dd+ka(2xxz_!Rt&E)PyX0Smns7l?V zWx<r-GTA1$Mrs=@vrFUhaHc^s-#Z%yCdR^|?5+JZ`2fy>&&6Y=8pY|9^EuyhEt2Cg zx_pS~-fnS()kAou#l3NJhn5xn=ZA>g7~kA=n~#h<e+z2i-iv6JOEoR8t))BhzzG|! zx0~0SN%2#ZyWGXg&BXlNOZDj<>D+aeGu;=9P)#wKh4*8_sMUG>YNz3daD^z3u$7?_ zb~TpGH~;Id8!jU6n?Nm!4y8`_fRl)d9O?p<$E3koCr+<~(n(CYQI?6+rdYTTle0bq zrK0nshAQQbL9_Q*FZG?S#wNfDlFf8B@%GV(u+WlIRaEiiPHt5=l%I$IZd2nZ_nswL zH!WayHJHOrI}S|kI(>rvcTEA!h5T$SXpm6Ovn`3vWpH_GOCe}ged9Y$1~(3u{`kuf zh|?a+KI6r`48AD$HmXOC!;Du!cSEViaC=K7fQ%Fnl(r<L)QI^My`wef7BhoGe4@Gw z2kL4-Tv4zrp=vA`<~M*6q4t$Jf97zE#QCkl)hn>1di_B*7@sJn!zv^L09zgfYR6gj z%PzW>a_E?-6a3VHHvXXaz;ADi4kqy>90a$@WP}+G320U-P+B+{h=B)k2yZ!Oq6AM? zZRa%v)~&?$?W0d!y{kchwAoQOx%$%3{-?FaDfJKBqEKlAj2-U2HEoaLd0Y9b<n*N{ zbQ7lMFnE@tj-!T;etHdS_K0={Q%l_@j?3!Vekl4R*!yCV?45^XDeyK|mR%No`NAub z`(Yd22npxBXQC1}hjy&=1Z&`RMW?RP0>^}bogKF&J;lA>W?nKq0ku*?I8E4i#m-F_ z3z%{~BMLhomPWZzVWnaZ45i0Yjy%e5XDnCfN}fj)F##&$kMl0p$8XTzNSgINu_FC{ zpOg?kg{J?du6!{w{Fl7)kMW5V#D7jo-r~PMk+Sj7x9V(VC~+&+3?YJ^ry5^lQ)OQQ z1Erpp)*mL;C%(f%wcrk8(m;*qhp?@zA-Bla(Sr*8QH)Qug$0RH&}V#07~;l4ZX)aX zeA`0kJn^<(X04MY?B77NZa;aD?(WrbzyV1Q=cmzKud$hPgKYTK$AkVbcnBss17}sf zp}N%0+4NNqb!o=0{4jH%@!YC4Wm6y9+_yv#pc$-km`+jgmC$>$g0>7c$Y;}a6No-j zVlOP$bnJHWXUYY!#>o=UhD4p8_O7ay6miaWKt~}%+4#KL0ipEIc{7qM#ahD{Q8%Km z2@#S}Be;`G_r<RT12b|YXaE?Du;*)$@N3nKMjgd;;;KOhhDNY%7dxazb&?jq7i00# zeUZ3WJWlt1uOsKMzBlpmOu;=LI#y8PT-Q|2Nz5F%c1p~;pvrW$_u<>|`R*r~z~%?J zvZR0qZz)x{@J5UD8>vQ=No(4plF=qBH)C%stCITVNSvWm=L4FEAdN;(7B@I9-e##p z-gX>1tC6xN96an8fK!&ob@la0!lLt0Q74`(O>e&ssEw+z-sa)S$?%Sqhok`N<i%hO zn-+)pgf}HL#fN}$f!R~%92ILh)iL?2D+unANyWPN8P=<I<4AuVQ#*~Q!vXPn^QqTl zx)RhcTg$4Jd*{dNH~(#miaK;|#lG>CaF&rM1=zCYh7B!tS=w**WR)9c9orwAnG<3! zvsdZO^Hd)4+-nePlcJya2(PPBzv$1nSGb+U_~45JL(khR5RmAfO$ZIDz0!>+zU`}} z8{vaZ3Upn?%;4Vo$+U?G2=0+NkmRTm=-lPV&@0{f?Q&H$PC{`SAuKyq^$2ZH*j3rk z3EZ~4Vhv5J78w;=@oO059+hp}cpXiZKf+}U`7b|6HrLwC*V@Y331p=K(oJ<qrP!@E zL&2IK%lEL73sfhnSkQs41`+I1jtI?btQ1g5Ds|(ey5ENEQCreXWuAc*1m5asg+KWL zw+4%b(R?LL9lxT$mLzmq9$1dBt<fWLfa%(K74@y<;_S->WMZ&~*Tzhfr?g?RHn|~T zZ3^;i5wy&9l^Y})&Dyb=D1=09^}){{PtV?12E79s9mis?RB{t&9`a-+9J~qJxY8(7 zt#H)c$Z(JQQ8vwjh@{6L<vrh*;R-+OrO7>p&cR7{MNJzELXq0oHJO9n&ZKb`%hn&{ z<NXY=c;l!DBS9eLD7^hCtGp7DMpS&l(ssuYv+8dzTYr}QM*8JGlJ#D+TS-lnt&%Z; z)z-2a`-@%wgrQ{D%Sx<=p11LZ@q3CV(yJ12Z+xlsm^Ku!VipVU&=D%zR`zl5d0r_< zuwq7^6W&u*Dwn>y_tU{f533(0003P7hb3PC|2_+$4*kQDuYXwbH$(Zi^+^*K;6EOp zYONBVX`BGEw@#Nw0}SV*2{MYqx9YgLg2IdL9%(@-G=gw^!Jr@6kSvtoHJi$WB9`Du z?1zLQLK?`G_2LLLs+^$M{rv^3kY+N8<P8_+>r%UQ-&^}tnm0Ngye>BO+RFTJ_<B^} zp#=lxccw}Ju)bp6#Jv0}&ewvaChV?9$Kv*4>9nH+xMy52D<6Vw<PPXu#@$2Fz^Z;$ z>JRFyDda?;c?V)~MEv_i@IaTW<#k<jIw6CwXt_l{>~~0gm*R@D2%V1ikd&iiEZX`D z4>IH?r7?$((8~FY)_M&i^~>>3hP(`UC)j!IIyZhu6pMgGujv<?ui#MTZnU)ziL2wY zVhl0>FvnS1wrh(4u^P+?>sAW)CTtir`ZtZLuP7+iTaYOxbzjX9pj(0k-BuLa+rG6{ z6f!3BtTZjY>|*nV49eevZTDwxCM0#OcCQZ0yM7#H2GM#5Yz1EjLiCKiCBc;#f?VW> zn`EA~`I=48k#kR1+3@8C?<9Xrx>Z~Uo%OLfzRs_>|Lp9&Q~k)x9HTp!KIZQ8;t&%E zMap>Ea<iO8>9+cENI=L?rK0feB|?4C<+=X6bIb2~+FkEuHSb&1tp8iUkgKw$V=!u+ z>|h~EGhI5gZeK~3XCCR^zN3<FQ3DaF3!fHCz0*fz;$?}aw;6w7WMJo{<*I*rnY*%$ zQ`LkzWhI(+J2tR6)6O>aMiNIrlg;c9PbBb0*|oFvM9#6f^9tg@XoWIr%h?_~Kq={K z*+m5PHk&zbNEvO_9D*<}nn&M3!^oSmG+L-Zv+Z<Pn2nC(JP%vcwPreS<jdZK*wX~9 z<MPk6<&e$0za^E=J<EOgchB<wNvr=9S^$P%iqdBG4$uk?c1C)ZMs^Zb#@2|!){goX zuJ{&aR^~>ApBH<F&kY=l0666S+`zvz%q_4N9C^8RuN&&uWA8%puA~ad)d<asDzY$m z$3#P=N46<RX-*j0xP&6>xrD~YVx@kEl39i2!{Lih5n9G%i@@ib&H{ri=f3X&AeWKJ zXL7s3fX}iT1vcw_Vst8Ik<XkbGPEYjorkr8SU07{>@QX?GVLdMs@%<{>ON%EM&QE+ z$>Q_Cjn-fJ;$H<1HH-(Y>6#nGTwY)>=|f<1T0N%xM?9Z6=E$c%+TX5VxQ)M^usk_+ z+L+HgjKjR2quyC~`Kr4JGlY~<BvUX|Kc2xH_dG$~ZF^9DozD_0Q)XC=N|pWI#m$$s zz~d?;tiaw#-JuW%MVAf!EUpKM?!bA@DAXm1Igy2WC>A@VmUgsSHRVub*pg$wl_)11 z$9}(Ms_~o>1E3GF5{-8SUl5tG>7-_NLoAw{;L=;djPw=s%BH2~7v5hAIi#$V8?nKW zWz{*r*LpA7v#lT9al|Y2Mp>W^DqyI5l=&JUaNB}h7x3`(Qwf7Ijajjo!XaG_7-m^T z^Fa1@E4@84X#wS=+5DQAa0mNxVkZ^a>ZPbDJdId>#ipA-b9xe}EpcdaLs2?CyQYmG z4u}Hnk$`E?Z+@MJS44!ff$6WvOg2RF(CtLf+%x<@&fwQ^dM?wpG;T*&lipyKj1Y9$ zy!JM=DgHdg-THxgNfTtjZ#klM#Y~3wX|vi(H^C##v*8}$c=eCn1&+huRGE{|$M6OO ztK=GxM4KA*ttaSoTm?2FdgV8k2k2vOaWUKPZ<ft`C-I}V;oeG`!gZUx9sbXoyl>Sh zFVbk3#8mt5EHjt4wILG*X6N>4y>A_?=MG42)(}Ed1kum}z{R7sy0~v{`7XDky?J>s z`uI6Q)H!si@0R36SsouflSts`I@DWGiPK4`Es9Cg8VpI0nmA$64HrwIeZtHZWy*na z5OW>+vFFVsWkHl<wp*N)G*Uw9QVi<yTZImuI5NpB*m;h_6#m8VJMDs)nb(3!k~IhV zrKss`%9XJG?u3g2mWrbDwC=`3g!c<E5fB|DJGRw5tPx_Y-qT+NmmkyE-hEbt_cjEK zP#5Xg;J|)_ch7&@$`q)V<Y8Z#Im)uGmk@TphMLET5W*Oox>`kYDOI<cW`r{Oj&tSi zY4-!Op^3(Elzh`0k3znaEZ-v^D13sYXNK?n+>_33XV7XYe$*QAc6wZ6jT$t^8dazR z@W4-(IF7G7b|#roh<i(Wb0;p7w1Jo!Nn$41JCpRjHE?2zt~sg8IpjgX3@7i|NuAfr zvjcLaoG$JPLuv+Sr@Xdl6#x||?dRgy?iHEl(STElmG{{ovvoMmT@5vf5K0wLM6M>A z(dAZgXi}$8<?T~5Dg*W=_rcWlgonbzPf^lXUM*L3PR#2189(Q#B3e9b68cu>`f&0Q zj>=U0&Z2lgaho=t2SS13`dW}2=ELk?OsL46RH#1)J^p>8b~0yHiLGAJ?W^@|%%Vn{ z#ZoH~^HS$QP|>s%UnCH=rBq|H{&OC+?UT5~>VWcXhyl0IhTCr!G#TH_vAp7m<tJZb zH|9ac5KYDA;4&Hb2~{q>%={ucx(eL1ayhiPjtd-FOyet+^^-fTYwXM6gZF^^Ef0mq z)LKN(#@n4kCHM1G+9qnk`*iK=(kWSP=U+A!!-G8-CGJGnBI6OplLP$>mpa0d-`<7M zSCGTM;NHDH&ShN%7r%kc8FRRHq*OU_)EIEI$YPe#K=M#DAs&hyVCGgPtE0IHzWV1e zAxv^(Q^c?n@bE=rx39mFks6ab?*}?hVZ&W|=g>a9WbAvH7k68bi%H^2%D<r;J){`M z@t9se_U|ysis<;hT_t+f)e4y%g5f?AUps4E#^ZQ0hb~m@S{v&ajIlwaH&QtoWdNBi zV`js5npR$#c_Rjd=6<~1FT+jayA^msD)o?Bwx-Bfv~n5EK3-jE=_zDQb4r@onWKgo zFV*TzaeJ$tJKp)g-sCUZ=|j(vfoEHZQ5KTN1{IkT0_8~<->p1ltE+9Wr#}eRmD4Fw zaU(6zn7EQ)Yvs~zW&02ZNQin+fHA}O$*gW5Uc3pFTh<F}o+2V_$7+NM-`cThgsro8 z@HoCe=wgd<K46y$-rC`rqjuyKcjrEEQm*znY8gk=@=uceB9qOC*w?GrTIrr&%&bb| z)VPvQUziemJSwMk1+WUtenR+XiJ=3%8P@cvHY}eblW=D4Hy9Bq>(E|4(ir>*L5uGp zO4F{A=wP9F<w73vK=Tt@iyUnx-EC;>1q}+Y-wD4g<8j?_3C;RkFR9)xgk&ij74|6a zz!gA<b@8awkL+!Z70M~vtm+zi?BA+sc}<Iv8q1!TRHkR?)OiTQoLE3S&B9s)3n1rP zT(wKOZoMhZ*eIve#<QN+p^c{^o$Aaz<tCM9A99s_;}E{Zj%2f~r&)R>=4IS`6>+!} zYqkQaBc_+SU>l<Xz8RSWy-9IFFU$8Ol?!6nLxf#YDCDzV@kQa~WaMRKWj!A5GfPPZ zn_XDy@Xrc2Q*TV<!4;z6wK3CrMq}qAq9|A9Jz|pJz}$ruzurX%q+5zx^9!~^g7J!D z_7`t~E@2+>dGmt;(*=8B!@+&Ldj+p|G)jlB`J&9IFVUhjWM6B`IoGS~0Hge~hf;8{ z7Wy<!qH{3OX4?>wrddHG6CzKIEldl<`G=o5jm)Kfll0jf=z@CM=D>6%!bN|&%%Gov z!~`t9ySzQ+?~$v=(yElbQW7NNhSNZTgWqlO8M6F1D!!#Bv_4j&E9pE3PYG&QE3j7; zYEI{=>7GkDANK>OoSrj#cZdOHm=luGe2K;n%!1#mj<XcCoeqyk$W&>Frl?BsvXke_ z8Ktgo*!$6X8jMQZPr(s1&kd}%j$?yv!p5?*q<|lmZ15JyaDT8rB)4<hs!!xtMW2Zm zu*P1gjOd*=9jK&?zOdz?%W2d9jAQ2<H}|Lj>p}uLGZAN6R1;DasFRopK3yuMGHcvi z-I@HQbGYe&p_kC#c$amDfJFu0dX#wkL-*;l(i8hvBCcnaA%qwZqyZIWT}93Kl7iAo zqsK%i^sfiEZ21Jw?xV4U$ApleDiGP}vKg+DgtYhOqa76_mjwq_P;oMGIKaFX^hVMt zmtV83xa<7`p>1QG4vL98lS$Zo0)$|e=Oqp9mVn~@cZ_qVdmgFKzT$5CeW4wzI_ric zIWC7LHGb>C+8HOFlF3X1IETqXf&(#|Qu{1>4W!HlG^$8<^oeq-ghY}nQ!d_p&}^T% zgs!O{R|I=gOv!3o@>*IV6E(2GnwODNn_oj)+ayxIKLn(<$BFDO<WLwgV?Y~_EHQJ} z!ko$X4A;w2Hj|dC?JwtmU#K?qtK+C<z^PK8hBZo2s=W}RJ5^6KWNcfU+}#dKp7zId z`+cVr!0s(<P|d8vjBGU894TCQHD^$>jy8txAccHZkO1R}Ah`)+$vfR_KFADS&^Eze zh3lUf3JcfhU0TJjmlj@NJ*cADpt#Ow_Yv~r(9}rmq?n`%hr=7VFBMG{*p;%DuhErt zV;Y4Od_<QYo2TdzGmSzsaV)#9=e?hMx(dCK<0*)e<9_RplZ&dMJsXVeG0tR%xXl;m zN#7Oc18)V;>B3cUGgKJ4f<kvLFXO`O@#%q6z2X$r)cqnV$b$=sT-W@1VC<3+bB9h7 z!Fz-Cc}AOiA({(LJ%3Xfs<4S4NRo<_1E0paGF>Nuk;}QXCS7uVR3kAN|EhcK&y;pt zb|nj(SM&Yh>sv%5y5(F(_fM{g!N!{s`AGNgvD%tR+UZjF`a{Me)sTBnV%7P+K5&}5 zMks3OY$-oTdNC_4r3@NoBd=nK66Jw8SoDFLn}8P4Sw)Yahr?nTRZmctk^|eC7X~FX z?{EMOGm9k=lhuS01G6uR#2K*zmPslvCa$hZg<4eLg{#6A*}51e8>VnC!d#-I{m}K| zOfI%j<Cc_@!$|EeFS&cP!BoK`b_iSDQXPIq&WRCL;^ru>`GzrfSWZ9`qPWY3#8Nqu z&`zPvXHzf~!p8f)oJ`IrG!6Rw&|sc~FA!L@S-mHk?;5T&C5|!4tcf4~Fi#s_(FN$& zBZPTovclZj`AP)9vBKhEO*h_6pb@+;!8`zO&3kUhYeIpIC)7W<dVtivEb2NsU=8cT ziWK|H%PkpCCuK?#XDqBIOjw^22yh8KN$saL_`oOvy={#-cSflsLtt9N@Lx1su$Vva zD(YVI;+3?Yn$P)@8Q~PoYHkJHDARCjba52P9eY$&acza+;(OfHsQ0)j%QF!iXXe35 zX0XBxEo@781U{N+lU@xQ@)4J&JN6QwNF6s9SV0f{MBg+5a*+VwN&`dNERH~r4Sey~ zK<QswoV^6grG}RyBY;ol%&}v}Z4Gz}OgGoFgyj>8=k)vT`MbY1E09HFUUKabJ*3Gn zZ<|Jd?(4HKB994^2;{pe)gym?5trhj;FZNPE#uuvQgoetU*39aI6r>`W4*4RrU-1) zy)F2*W&Y`wDPPDO^<$E<LNjgN?JNZEn}W@^J_zqfI2XiM+cpURp&?|GR_UaqbCmrK z^($Bv>gaoz+V#BBV8rL!-}WHX_*aN!LeS$6?vwj#vtVBnaPg`E{hY)PErf8BIQ{9O z)5{Q51J5<*-YGjE@(G2*ic*3){A6?ZQB1i|tNz-A8VgBr?ORPyi!EEDoD&zLJLOwQ zdM;-ajlgBM*ln+M{>Ia04m}!fz$!=%7Fu`2fLAE9zo|fT>M|~5Wi9WjFCX4P68_QB zP0aNk9+HIhNK-Zaj6aq4gAB0!H#H$6r(N$E>ODuXcO}JUlgmWT#BMpsEIHz*nY66y zAO7sE8+8Ra1HZ)!`jkVmAZNTN<pN16uUb-$kivJKDKDh9*2G*y36Zh9aL{7VPTDJL z;c;|nzwXaKgtf2}ILLZrJ3#Kwoa+$|P`gGOhCQLhUi%>vfXf!J$q7OUqJ;gR9@4NC zV(2ZCl__ZA%Ic{|T$4ZXQhV_`3CJUw{CNbM%e~<MVkFAZ5@iT2m2eXH9E75e7H95= zX5JeMxg{Z$s-A(_h$eg$4&tFX;aD-CeqF;4x9co!e;`bRLs@*1k0I=A_?+*Aw`wBi z&TDaKRJW_>AN<?@d^FBtcV{iJx~k?y|5;?ynI&zLC=@E@HfH@K!1tF>xkogieA{#z z1H?B$;21Q@gMg87oU+5}Wpoa5Kzkf;6{i&*XWfCXxP$Se-gJtLdt>&>VN_w?!gR55 z^u58H(XY)j<AkLsNYjXqKmo;lh|3R}3|CMKMP5TLKJ}EH_-}XJ)Te5O!#HB+Qq&;L z8ws0EJ~hmpCP#0Bs_EYxr168FJEs|P4i_=N1Jp-_imSFuSz9~Vgod;c-u7<WU_}Z1 z-J-I2Cu5_SClQ^JT@sNA3tD2i?iD8B30>+rj%as2RQnGYpHT7aq5hmcg#}A~qZh5V z3J7c=4`usQ0}1x6>FOee!O9xHrN$}?N-U{4J&;R2CO*i6NQ6h-&sppbX%G}#hliQy zYuta09jSBp@WxV>R#<t>-Nq_XL5|*4jW++KeI*+rPoXg6>qIT@*Px?lJuA8qZgcVD z_GRPBI}MA4{06TWre@?-qfRxHAANB-9P+`YgZTZpJUL9Umo&{}@b)^Cll<;cN+?aY ztG>hWwE=7SKAhpj7U*0nXFE=U?gt)WAmj{it(<Sc1NHm?JOx_yKnJmzZ?^ikv>v}3 z0hBVcD<Mykmf*$|;UAN{Q|7ih$`XCEG)Gd$U5a-1?HBl%+g2X&B!cwEoiNuRUg@>) zt{0jjbq?IP6h|)gO7uhqk7XMTY6qsp(#2!FR~u&mJ3@6XWASh046kn7R=77o*ceT7 z#&F}>+znsikl1TG0^_KCev4IGa&T(|IMu$^DXVfjN8)Tv1jFd!@Y7|eFMFCkd(V-p zuY9(NdjF9lf$0BAJAYeVKbuYeN<05W%j>5j_-D1`Kj}sNy!}6GL;H02PrFGzUH<dB r@Tb#%R`vNWq}|W;Qv0*~0LXtY|KGmmPi5`DS^iI9j)u;eB_8lU!pz)= literal 0 HcmV?d00001 diff --git a/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/BFitResults-fine-2015.root b/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/BFitResults-fine-2015.root new file mode 100644 index 0000000000000000000000000000000000000000..dfc93c83d60a443c3b39c08d67464b46ed2ffbd0 GIT binary patch literal 7771 zcmbW6bx<AKvWGXW!QBGEg1f`U6P(~7xVy81OK^902@(<<Hty~ag1dWwU>kzG9nO8{ zzPkUscivR3nXa1YS!-5()jhv%X9ovY0N}V8005W*02DgV0SmRYp#ugweBu8xumk`A zbU6S_ER{9vSQGYcDxob*X=nQ>^zQ%ZXa@Xwk&Tuh7z9v$elv77002S#jg_q#xx^bQ zS5-3?kgcl=i@BA(8H=kci>sN7E2}efbyyA{S4WU5rK+;BI;HtPe;k1DXV&ML0Wf~h zNqwQ?d3FGREbxE(B-qG5D@QbZ_GkOEau|v~{_yHDoKi?2M;1q0kjsBxKPVup!2%2s z0q!vM>WVVzDtJJyi~bG*m02a*#JrHNrP0)T5{4`ij#k2icCQL8ZHcPlW#EPsdy_7? zCX3_W6`Ubijqb~0vlN7(nq&7{M*g5;u4W*V-COngNux!NRm^oZF0ytD8F#Z-C4-@5 zg@@<g>&)L%d5sUySTDqs2IPI?Cb37(jd`O`>`M%J7@!Fqr(LICh0SK%mm<3MkBS2j z1v`6cR;UFU&Yf4pAY-)>k^$B-tY`V?rtuY*R+=YcEwqaKh1E(ZEYog$J5aY3-Zyo% z4lQxwtTyht_cYVRpl*}U^TbT>IAo{fmoWo2Gt+<_8q+a1nn$lJ$^sC%ih1BD9e%RS z)u|Hz$TQO->}HbgWQ$(;sx>1}eDa_|tz&i<{K?pB`jx!quprsp??%eh>7J$-=h0H} z0_;uUDc~k6d?=py_ELt-&dS@zu=K^tmsum4<knTuh7w!u4tWCKGatj;?02uT?-Kfb zXIq@Rgq+W3c*y37y*K4^#xBx2#3#<8Fm)#O`xoX`epc5LS}aEYL~KJ6_Pgu<!BP{X zmacqyx8}Orfi%r`lild;6H?(n+He)2iRMp^^;wh=;RKcG@galZ*OQmDuU5&2S|2aR zVBQFw%$}JlU5@NVRi;mZx>+Lkj#{<^wLTGHT+OF1`b`U$no-w(q2NOj3B$f<c~swT z3S2lVxMR5}%2iDkbYXg8_?j6g7p)v7rBXLIqnN+n6FOpxXB3BMPc%HBgtAKZ2P?-Y z7b94Tj8;{e&GOvntD5MVK7sYC4@FfqOo0FZZ16LR#?L5n{6P`!%iky>Ls7J|vIn{R zZ9rj^hXok_3r0&dajcyCjM+elTcW#)C6h?Feg-F6z?J%BCJutbo6_Q5F-7c1tbJ(- zsrw<>Z=n(IOQjXcYdECwcs4N;B4tRK%A+MPi!a{=)C&F1xuV661W_Cxk9uF6@U`A( z`aco5!@E!Tr_S2Hrw0VA7le+&VIlVxo`}S&ou{=?NdcCBQlm<Hi8i9#CVw2OD<s44 zulWNb-3PwPbvT2y^L82Z*KKgEfQm9!FbEw#Q;07}M>5%44!cxDW6>WY@zl*ma*JWk zbY~j^G9wn~C-TTN>9&-_-QwSKN^oLoWG=mirPO<SxiAxPJ<K;6*{`-;>zrcw1~9+r zYgC>9y{Y<E-MGHsQeGpfG0%gPIY)qsB#G=mSMVEGid&LojQf}bnXhIuB?Y=6?^ev| zat8%{F|du)dnM8hyc#<en_YTdxI!K*(i6(zX0}OEFBoiS?(n76&Y)$?tM*C%%4$NT z)m@2#|LIr4$#lp-+1ifhi?toa_O})J)4XdhfOEUV+;>kY6B&FCt|INtjYr$I)f-uu zW5S4gKN`%d0rSYfJC&h*u;YBLDA$Wehg_L)d<`wR$}BZjPx(hfCkSmWYJ)o6&M1=^ zg0`v=kTy3dvFo_Y;r%*r;58+vDJy9WMMKFx6|IbEnUSg3fsSK^r>4jNSYplYHMpCt zb<pejVhyQ~df?7cZmOPn6X(W61q6!D6CgBT$86uZu!aDid|`%Tg>&=dHHkHat_&MT z_fXi5@E`B)8jsfIay>9URzV>Qc!tpQ8A7f<5F)bv4IzrDgPZ+xMEth}#qbR*!1ON+ zm1UB|@%(<DG#uK8zP=&a$QOH|tU{=@DrT7dDz*<10UJ5F+jwEo2v^#n8c*9R0<Su< z>LI^AAf6^@@^#Sl<jeP8q_y)|L+B62TGKVmSH$L5kuG-IZZg;IcYm+)!5*@Tz4P}Y z_a>PTYXevY07Oj1{X@=*a!BVEJu5G1!RA3$fC-@zNA60v*s8DGk&{u*X|QG-Fmb`` z(a}ppR6qQxJA`18fbc*uLp4I3rGbn1n6~qj`MV{d2(hs*-%;z*bp^{_M3?B!?{1Gu zni1bdar{CV2$~9^>n9%HBHIExY5xJSEpR`x1^;9u_78}7H3Cc-*Z?ypo%cH5^#QJ@ zl#ER}DtEPvZ|=)IzM^;GQ<WuUUD^s(aLo$sq~Xpn_lBEpY`b0(v<tTD3OW|N$UVhu zLy5Xs-}+RNd0LUfHpgQBOW$Ji<6~rBlLWg@v(NebkA+{&Q5)o(zu$`rb_l=LFP@xz zx#IWWuc~qTEvH)MEgOjcLk?c0Wa`8^6NFirb?b`!=d#Jfr)%9eNJ_*$+AcYS)NUQo z4+5LNW|#>8zbU0e8Hra2V;)+@xFH?&gUtKNe3X%?^+Unt%^x;dAc_6ER|2=mCh55` zj5NN3pyUwD-EDg7Z>AebE_iWYy*~}&R}4lYm_<4VnlM{uyf(*|pUYq-0g|VLFhPK% zdZ!xn^l1DJAGTAcOh5||q|;D@-O6>H0<<Jg^p1;m%S+VQxp2x#T;LI2h*;Xq_OBVn zM>GOzOK8A@O@1cO{Fy+WKLjE}x%^i~MtvsmIS~F26Uv|#EWrFP2IVc)RIswW?P{?_ z<-bs}F3M>YWYUcyh_~RlmH-JOkg@OLDWk3W7!g;(r24GIv2;G$2i@hfplY#fN+5{S z1iY6e0MjCrpD`A|@A=a8XKKEl^>Ex9AKe&jIvw+YEkL58q09<lFK;ml2FwFsJO?EJ zplZqN6>^r7p%{}dWi0uJltjy4<LY81;wHK3BBx@U=aGgtVGZFAt@PCTP~tB?xpE>T zts)!wi)j;P>L6WI_sUf={1k1=L;jJ&I3_l>ojJBR={Y>(q7-AG#!Goh9Mx)P3=M{) zTa80oL0uY79O_;$aXMm`=S2}1BXc2y<t4AX1qkf00QgQF)w)myX$Y~BwaZ>V7~xa< zPPIld=UY(`FbEtB2q;_33+*)Gvgw5JdHqbutv?<3{gpZvHZYr?MBa4@UgMW!dUuSF zjaU4`;P|YTPRI@WDt^G8@*p^ojde4c39}N&mec>KyzhIn9gBPaudKq)L?>U@D{dwn z1UZxUe0jVDmSREIIzFQVIUY^9JLMas;D@3)8A4NoomOvP#aeT#2#E2C;<goYT>nl? zRj=JF5}}hZY*L!GP&CBE?<F)c`ceN%)K8k@k98To#eF}En;_O<WUIx{LSazGcfgZh ziAe8+XCbFuHabk-FZ32Bs~$3OhwENyRK)DtaNkURwQ<iU{q{<>o;Pad48J@l;DMu! z%X0@)RWjkM@Px~5!U(jXGd;0`OEs-5HRIGA0>3|X(#0@xaZts3d0=dcY6qFw?`{(H zOJhCl=<fiiL^*wiLJ9iM+>9puH<ZX$rhok>{%t^k6#6?iLxZ9UR_@Bf>};FV0S7oT z6h~NEZ`q6j{!q45rv~9uYH&fcMpeBfE*#q{rh;7tPAN83568ZBG+M!=IQO{0SJJ&W z4CS=&$-58fY(9@LM%kFFt?9StUL6fW4)Ay@{{Dax7_tvG8$p0D|5#uRK<sr=sj?=G z15a`53o42An>ks$lgLix!P7dl=@{oU7%LT+Gn@!3mW3{)uHsC^aJ+ZT=+NIpFk_0= z^+bpX{nlRLZ8(N3<&8Q6lHh1|Qyw}rAad+2SJQvLTgHxd3UyT=GP5)$oA${1(MYa~ zuAvVwy}|Z5EiE9wDcqIbVnB6{k+y;Nr~fIMB0r4g-wQ1G{uaOJQGF<^1Md_@-cUOX zPEPYqk+*{LxZL?i1R8<y+`+ny!4o;zewmy|#>274_)U<B?BVO<mbxnEU1tQ~3B{>T z!(=zV4l=lS?`x?1saNX>=5j{psZe!SRwA!3aZh4-c~&#Od-S;%glNEXao>?6;~ZkE z51csrxv}NdI*Cul?R_mmeDm#g=6Y#qj@aS?9~uZtC}l8{ql$ygLL1UyUm8T(N;_+5 zQsnln?RGh+;tlEbi8_C$$=%Y@j@fq91IeKiippfxtkD^Sp9UP&sC~bHs8!7LG;&N~ zi>BuGb`<aU5&>gGPbg$4)(98!s#_XOT@{!XN!)tS%N$6;zakEtXUz@?A4`iJVa2qU z02m@e4hpPK%yoG>1~R5viYIzn$GMdAA4E&c)qg{w^ZU<S{Nh>1J!=z`KhQk~!hc&( z>?Hrr#eeA&ht(*7+Hl_#BNpA#4Mxjp^gW3<RvS_W6|64gnO^mIQ*P;iOd|9l9BVF| zdYL#_{9x*>ZfOM7A3E>jyVYZ`#G_6=yF9&VD*VKr{mJXG{rAf0xX>wzr663p|2|LH zFHfA0DE^Y`=doi*mN6w0W`b;=+4RTkrJ}4tWVY*iotELT`zsui*28SlkgQ~3*fE^+ zZ|sbN-k&b4x)L-lqpAQL^mv3~@Oz}J-6nIsrJyE=^)DX)fKjV2RYLk*SnO@MuVfaH z+kOyvR=D)PmQ37As@Tfvx7)0Fw}a5RPEb&Js_vRZz|BT7p*=eL%`zHL8ot!InDPsn zq=QSBm88E6FcP1qZmyW0`yp2=jxbDMAeI@oO)Cj{^r(_4edZMAfY|!p#-t5yP3Pk- zENvK<+y9=8CDYRvacrI5)T46eYxuzf?x+MmttuY2U)=$BM}c!okCUtOQCh9)1F3%O z4{=5L2V4GA9@1M=o+CE5QLBQ$k&EvAnj7#<?ZZHiO{iCUD8Kbj6Gv8q9%q;;B;)&n zn)93PYQNZY*rx{Ubc9<e7x~L*!kAKKBsviBNRh;(W_DraeG$sy?0)SR<^A{qEmYs0 zCe(M*_DuMt1$1oEGbuxAUmuVTVu_`Dl$A}J%o2u+4W<MSA;281Hy7qlNYMjphiL%G zFYPwye%23=K8Uiz=eF`Koo&fmQ{A|^(Wp(KaxKWA!#)9an}Pe?V|QYe@oR6Nijld} zdp2hc1BT~%08|>&u0AJXu|J9UPn7^*ilDAw<>HF2=IU%_WM}3qYj5s=Bjo@xw)LR2 zwX(M{GlgDUT%jz&paWKEp;~$7!-@l$7QU#^^<ubPmY@?*{-7u~Jdqo$6Ol0a@t9u| z&Pw;t?uf=NV_}8ZD`SD!%bXXyB%yT(mwYDygQ9VX@(9z5mMn%0i$@i=k0!GNLm_sE zrux!DoUy}Y_ousg3)pRiAlzo9iL2J_=V?x?*Lc(3;%BcXIcY|$h%Uk;5jy<UTjV0{ zm!V*d+8(raT&&f{xR04qh5M{8FV`2fnu91b?tUwG+~Irch*fHBFE#y&R5P+Hzi?wy zjeQIlczduBDXTuwrJxf(=3tOVp_^N*uX`%r6;|32q;euB9FSu#EyBa6rTVG0D|WiH z>DS4Up9kmf0yj_VM`3-#7`#1y@scAp&Px~5^}8sswZvdnL4`e;fa2w{lYxq@MJjB? zd<uz`x=Su<7(*OAhj$yyiwK;KK5~xg1SJy;P1P2dR(WYJOVQA=RNPy|UGd8qHWza- zFr9b=WHJPpD|_}Ied?&O(3&j|@^~O!G<TXY+TFju-_?W-Q^S#y<!`B@N=db3m-x69 z*&Cz@e3tz(;*tZ(m7LK3p*ZqxC3abQt%7kbR<YniU=(c~FL7f#W#il0Mro+9E_^WK z)5MC^(onJOyOJqKG)m4}D$R!DJU`_QO20HQ+U>$a%Uaa)T#iU#ZB|k6yWb5qFtXer z00M*#W;8#dI-CavVl+(Z7gFfTOZ2oIG)k#6pQm|Q1b8m!XyqKO0UjuN>KJ}H_&wAf z9~Gw%da$1ICJa4E&~>mObgAQDea|JFx@B?Q4ew}L=$)I~AP{68>BY_W6Ch6v3z5IL z9&YNqsCu)~*+H%hTB@AaT&UHc^{Ff_p2)u$81`$eJvBGrmBj`>eu<C!B}BDz!N4qI z^{Hr=f44+@+R!R6kIi<O@he$hzDBqPv8KF+FjGX0i^Hnu?^w~+mDWA5tQuyo8bq90 zyL1gd1H-~<xH|!ww^iSF)3?5_Y(V`I`j#)8U}}S4@#zybh%}s|b`Y?xp|WklfDKQM zTj<dlx}*N7aI@$e!TJru4<#0hJgvu9#B}krhM2ouyDO=-J99Q>LZ<vnwKnY2*a~mo z+yaKaEq+J7?3`z3W8=li&cwd8B~nG7tUHhjhOqZ@DHI5i48bn0F_I?0ekm_R(%Ho) zJ`)r4L2vxZjXyxBw^n?5mD}NuiGJDQ|HypuWfnUfK4$|omi$e$iHV%R2a${d*ON2| zu8QWBzkgdwK2yeZ5fKXl)zUHs$L<zrLW%uVdtYg;V1KPWkG>qc<}KI-o&Bq$1@*>b zmtv#nktc=E<>S1Y`}@TGo;R;1Fh4X+xEhJo@Z99m!zoWhD6Qu10b?1K*S6M`w2@2S z6Sv@rhl;%lK!2C&Y)H{ICF6dXakl?X4js?lpXU;Vsr6Io-aNU(ePob)@CYOTba_;+ zYkjyMnIUh+sSs1;=_C@5p&%ux#94M`<=~vljw0V+U73&OqT3h)wv(-|FPo}~U-mhf z4APb)DVJADmokfX2}T#D)fXn@z_2cxXMDkH<fvQ~S*k;lSFnXZkeS!L5-mm15Gj0J z9E8U%ad19FeD~mm211`R!FeF{p+dQ<(!l8!xL?ojVX4zN@Rg=<5++YzM-36?^l<Jk z|3ceK+GT7x19ozB$nxIIg>bScA`Pg1Btcm%j#u)Ny+-D0knMTXGAtguGMikH^>~(# zhsdY+85Uv+@eeTLA={Xik;8-ZO2y~UnVV%vYE=#KFAdWXp!FYyDdN7RYc`%XnOnWZ z#$00_8@kCn8Iy?|9@k=?_=sf>zY(|mbsMJV<a@oquYtw(pb|Q4GO*NycBn{74eHBl zp`ucK9lbRJRO1n~iz|{5)0vNishia!`;!PnNs66irfFoyH4<5)0Vm%L6Oon*d^Fu{ zj_l3}n!5ut&}+Q>hhO+kYVOd<M3Ui9+`nj+*g^DLEprY<>FYIzA-B)hqtzeeF`6yZ zHb3yWeviU|_ezORRSAnx-pN*DCesj<l&kc5ft!N5z@@aQ!0q6mrI>(V6obHqB>#9L zPC}c5w|>WS*|++Ta@Rx#G0vVF$Cydt*o=y|BqX3W$0dm^e0rHg$?F8{cqZDtRMW|V ztjnq}#t+7}D{Xi)y;{k#$n@dFgj3;_!KC59dFc1R#Z%p~`5Z#>wJ`OOikIZsX#yW5 zhS(1Gw?{cZ_Xlb&d?#2rhRI`8e$7P_>v=OY5zd)cO*BFQH(EoTJhElmA2WT%1L}9* zP@_eVYSb6X?X&hwVogz(p61(sh$v$4HpS5(sdE!uWHDq#8{L^8*^TJby_I5Tv!Qxp z>R)Us08GH&^yjh*c8S4VIB$DP-#;)g0PDd|H&SC8TvlTmJnv-~?|E(9iNozGJ>CA9 z109TFET}ceUacz?I&Qw<#?=bnb3?O_^TJPXJqEFZZhgpk2%)|u>jIe4P9CkhIB!Gi zP|Ac^T^g^C^|=>gJ3&Auq*{o?ZtG=Qz8-&}jN}T|YaRZ(cPWL77mBOtJrSCRVe?dd z#kj{fV71Ve+1U+N*S97m7;fclJ$gh(IF&0g-ohi(dSU)^%OBc4N;Mo7@66a~zZ%iU z8pv`zk&DU}YxadrugfKJ)~Ub{Obw{wWvX6yz4jsJcS7H7C3&A^WkdCv8n9%_C=^|n zP}-v-$>Or9Gd9rE-2&(Bz^Ah2sTzxN6H;2yPUMw3k=fO^Dwxo%%;0E}b*&XgdWxE* zhnu5?@#RNb-lLS#&2D9Su~x0?QC+9Uw?f&2ESn>26XZ<g)WOv-Don93UjqLtz$!~( zC8VF5`spL7uPAY6s&Oh>L5DbE#aD?t+muTf3-C|*3@*?0?ZP`WgmIUG6pwXoKh3w~ zpFbz>XlaPlNMejqNSYnhTOB0gou!2<q=I!*hdQTQ6=?Q4kg^Y&x8|u%xk4}7A`+s_ zUxEm@Ta!%!Od(vt?%H{YNjj`UYn2-@x?`go9>Je=8Y_l5nG_7kK=Q`F6JHDy-CbFH z>FT%Z?mZ9O?fh!GZl%W`=hELZc>DWe@FEve(z1^bS6W7x!`aUp7LeXJFTo&$9vG)# zy5wjLziu&}H84oTv+k6S(bqtNb3wS4l-5C;Kc0`(_LNwQX>z}^k3QX{;)lItXNiFt z=uEsBdFzLF#bxFNsZsp#q<e;4`KylmhOqUfrMM4*PVDZf+M6b?aV?O_k43Q>#t4ZQ zj7!GwX^5te7?&8kx}h%6{0LHU{vJbW+yslv|LwC5<6*xNiIsc!WJRJ84sCkUa!x@0 zw|cm%9W$E3OhbN}dE%7!8_b&WX?2}hB-tS~TzH>*=_|8t%u`pRSIVuMYZht=zG?5y zTMqN_dMqvgY?jQ5q%~&w7ml^Z+X=HH4QJcw20n(=6sYKOD%h$>?KqQ;0R>pHgs~RR z5;aQ2m!0(onP$XR>d;G&lEz$bk4W&MWiJ!_eOZnKLL7T_YAtY(M6e+xRUpe6ENpD> zDL-qxOn5A7!92LC%MD`KgRz3j7S>j3{B9D=(zZi+rce;_vCzhJogt}do3}8KGgUk0 z!{%>`Z}aWSlkqy8k+{S9LWTXxxfyA~pu0Ng`C?ttqJJUFZObqlksOG5`HiVVXe8D* ztT}4kG)-9~T~(O^Ca@B$9|w1Ua!e#h7)lxr`iXng>eLFWFDoaGUI$iTmu>Tt^ds6| z?WiN!7YA?&Lg+_?M>^7vm$C0PRx<5_WM~B2jT*nZ)d-Qkym+ZTfDaBgd8DzwM~Q(E z)MGq*#TjZ#aq{Vy43Vr@K4A}O+wzmgpoN|Ds!6rlD-~n7lVaKxZ$T0^rphc7B86>Q ze=pTeb<3#l#0ll`V16L>Sei6?%Rra3?j45zqV`c1Igy-<9PecVt<G&0owK$3$gzAv zHDA2WoRk@h<#^48RLoUoPrVVw^vAG3nv7EEVDvm7o4in}{!Wozlhf2S;Z~7~GKmU$ zt2$5CW|gL^(k%GHpyrFGuL)I-u7Sbd{AAzMFtyrA^Tw@nklyevG)jK6!z@ip&Qs*} zH7Gk5s2i_u*HfZ&LoQ@6DPK#qeo*S%?^xeLe+O-;JquYF_GcyW+#ZBB=AM<r|D!z! zbwq?}<A2h(gFgPxdMi+ee_98Jy8QFPCDiGk6_%h5|FkFo-7oDwbp!zFzqkLn2MiT* R|DWxL>Ocw|H=ZQGe*tq80_*?) literal 0 HcmV?d00001 diff --git a/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/BFitResults-fine-2016.root b/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/BFitResults-fine-2016.root new file mode 100644 index 0000000000000000000000000000000000000000..13983fc6321c54f474ad940c9b77eb22c96c824b GIT binary patch literal 5219 zcmbVQbyQUC*BwHsp}Pi2C8b+(=n`Q-x=VVPAteNqQX1(-Ktg&5N$EyW7-AT@5hNsq z5Bb)5e}8`8@4f4ud(T?WJ<r+qu6>_p?G1HzhXDY)EdT((7672Ixeff=+~_tiZX*c$ ze-j)$005{Cz=)Y6bQupIT~DQU!6@x)-@hIGe;qA=zY{qd%93^gbna%p?FInc(^qzK zv13$Fc7i>#^Mbp;ym;-M-0XN^FkYCQ7mN>ji;lz{4)cJ+9z4_0(SKn7-#-qx_qXd^ zX8@q_wrS99+;s;4$o2lOSHk@H7dh^iyLy4Y$T2{F>#_A!1;MyZwjR75E^x1Z5@Su= zdQ89<T!1e|14KtPS)L;Jb4$w%^MzM*wFZuew33;N9}{dMX;EONQx22WsT<QIhV9_F z6}WgJ|M+R@6L5T@@OvdNwo=+^zP@_-Xl<pLy!F=TNPvZCu8AMg=q%{SKIka(G^h^! z5(M6tkHSfqd#XNu^z<i&zCvWscXdF)SHu+9b$K9)s6+AXqwwz045dGUxa3EyCtRf} zHb}yj24BFb72^}Js);bwjjk~SQQi{cQ+$%4%)yM7uEA*J8{d<K-7o@Uxl^a&`>F<Q z^LBKC-;n;adO8<y`lqaK#b^nSiAkn~S84n}Dm=t~fAC$6tS7LV5bC54r;Wks{%-Ju zgaRn7S5l*9h+J*QG=6U$&ICd|aj=2IiO7t5x!^1B!K3p<5vnxkn|}Wld+<cMo9A-D zYu~mVgKOPO3>Vo-x-sdtO&jRAi~n4^F;6&OaOWWk{Og5<+s4wlnuXyMBbtJikPy_X zuzAR56eO2>+B>zX8TewD(`CJtor+PO)~BYi=I6{$`pX-W;GNu)i<78Nnf($w?s7Lw z*SNBaNxQ6Y%0q|9K+>EuRN(o~l>R!0uf1%S4_|h*M?VjBV;L;j+{ADR$2Eg}-$cM? z(H4q|Zl3Qcczf*rjBAWp#tM|bi)QS;*XPvIyaJKPDrIm}KXvJsV41|ff({60O_%TJ zrG1-B>2@o4F(l5^9bs>j8Z6&2v0FS2{krmS;U=I}eLqqbsK1i<L~>Oc^TI2upK-ly z895~K=wp9!0^{k?IP#^y_$8EK!UAz#UWL$3VHQ($=CM7Z(68EyP~olCDa`6POR2fU zK1Q45{xB&o;iyN=cs8(IAoxz+`MHblI|XV+F$;?p2OMOh#((tZdiR{_jirx|-WUv? zYt}Hu^i?xNyZ?slLQi5SZ!n86)vW;l%s+Q*@3~{U&|kJ=1OH(=-e0!6I=R8U{*&GE zy$ci2^AF#3Npk;`eD;nUveAzqdf^-0f)pf{+*GSJseSbdo8HCsxn+3*b^uVQ2%4`r zlg|fYQQ-8?q#|Xfp)xOc7WuYvf{hfv(QkirF6H!`T|tTf-I^|Tl!+d>{2ho{;D?JB zDIL!Lq)M&<fW*SMNC*cM)Yv+-2&72EKb5x@q&>n8+vAFzf{+v;tLn@0DbK5klVP|a zB_xF#u_BQqTQUpTyO_}Z*Zgvp5O+~w!tw7tjU3y=f9UoN0WY+q2<57w{vA}dTJ5}3 z4=so|2L}bFPM<rABo{vC9S}W~iA@(@y%Ir+Dmi0#IA?Q=xg(;uB<wsAM2vn+Z6t~% z06?{Crt}330M$>5f^|lT7f4DLK1b7B6-V>>3m1u>S8@sg72XIbhR@a<=Iswm_}W{d zIaM5p#wfi(ZeZEHKJ?G#jtgn{2zzwc9BUj_cgZgv=S&AzS=W^cdl2aA)l~5ymQURK zvqqu)u@RReRaf^xch3s+JDb}<`H05MjfkDTHDDFtFWv=ZubM#QS>|-=Y4(-iN|!^2 z0}*YG9=btA%zVV`FsrV&LHgy`bmb}91*4jf1@<*zB5-G*<K~TIQ@7g3wDbssg7@?$ zsZm8#$yV)qf{PzJU1G<W_Vy+RNu0bJ7vykQ)UY2SIA?;zu%dDSIZ{zGq;qbiPaPjq zra31jp~ZyF#fHw1JY@7Em_dDr30BS*jvMIpgm@@WL5!I*aZcFP(44H9DQiYw0y>*U zEP}q*lMlcTq^1LIj0cN&G6F4L{v;1}G;V0Jl`+QDl=-w>)2rKPIIx_!6yvR}_ong6 z_2TeQ%&QA8EJoVFR3I=ivWpx8F8NF7EMiBH>AYuzIyc>y+9F8QFYz?tl2m>&F49w$ znfzhsb~5=1W^$Q{vWx|(=*IrC;m|dKaJdZ4-TrpF)Uz+#kzejF`Ty>E0NZ=|noeFY zpdJirXXR=KRdcg<Cj+~~tzG;cxH!2v+u7dUykNI~0tOJUz<#^(|1d_mGZ<3JO8uUV za?O(Tl+oDwkQ<dGj4+8#obcWiH^6c--FDq!bIq7R(FSD9&<5CxBIXnfx3N+#6fp20 z$9UHmJp?IoxR?~o@tXvyYZxH8ZGvY^e)3!$Uh6-7?LQ-YP*hRQsAi!$pP&0VaRFn^ z8?!|LH|k1UB+*^e$xk>a7tV1%2s1?>Aa&gY?c^j2*YVz&;KEHlrsKs~!<KLm<l<7N z<AU;~iCmT8>Rj_-jGmQa<&n?hXK~k|{VH2aF>3nbU798dqwW@YAk*ApbJKl|uE^4k zaNRw1>ChZEr6(d{hR<HLcEwGWHXrWI1^WqJ7WnwTxRy4zjHTEJkuTYKEO_i?yLj<d z4xRLxPf~M3HMDrXY_I>*$}BUfRz64pRevnRieX7+=KgevXZD_;N1(cgK2^y$XLF4M zu~S|eQz-$EMAx@f9!6QoxjdUoMC>Udp_(DVQ`NnB9azswLeS#4l_zr6#dcv3YxL#v z+eK~U5GxiVL;i|BJ{bHtyCl%(gPTQ~M3LIZVXqu`uHv})cdg;4sJMA0^e3+AIIRNX zu(#~>qO?uz51Lf!nv@)7-xdC`6Eh%*Glb~6^q#0z-nUB0n=8%668yC<9G-q`W3}Ez zL69|T=06|(p07n$Gx%!#wS|@A5|s=<YU_uAH@^F?urQ*=3G+gbsfI##+g1};pXXOv zfJ3PNjEP~+4jS+W)LqZ{(>?f4-R@5Ddul(vebL0h8wHMz$M?GQ$w<ECQhz<?g{?<* zG|%)*Pb^VM@(lNo=LbtLCPlu{IQl)*+<8>3jOy%Q)Pc`cAq{5gAnbuv#l_?Kr~O00 zt#$kM7NTmTi0h9D@rP2(Ye$?ss!p#ytc$Og$WK~2h2=eVndd5J=*@>jL1+y$Aky5? zwO;NEvX^nPt*F)wgqj|4kKUO)t5GSMGK0v$Y3NHL2z&e4^DhrupI0^Fe~eHOi=z6v zL^XT!iWE*CB~Uj2ScK@V+HjI$Gm;njbw;e|zbah*P(ihL%K2TJ*CEgFnuV4lf!&gL zy=xtn>asTNY$s(aK3C_=KS`>oqI?b*teE|Vd)$fSfBaaKEIX6*+=b>DaH4(-{Q8W) zPe`*sif#~d7R^PUi21Ry5XWSllJ-DOG7#8QE=<`^ZH6x1KgsO~AqF0Q4!P#p`#43K zj-9gvA5E!{ZRTdA3cSw%BKM~c$Es#y2?=R?pU<7~`vVQ{J?6Q2B7yZ4__#JdOM7o= zu4G@Gn~1qOzri`e3&>yY;lR3d-KEteyW<ZEJibQy_&!hC>{ezOCpK;#hgr$hik#*? z#L^j$)?Ucnkcs1*N3Seu8{w8dr~OPJA0fvQ3VfOhwFI?&RrNj2IM{ru4y14k5jn== zZhcj{fn;>QiV4?vJ$x1lKi;V{eX+e6lc8ZJs2N-B@A)KwNE586Em(Ho<POc{$J1zh zQI${NW!e;ra8+w)DEnHQFdw)#5pJYNSE-?$4z_#e^%^M6ZZ1tNK=h&v$weq?<)Kp@ zQ))t&SFm!ADl_j3>Kz!=7$g0=IGloCVe8i*?Zuw}0yuEmhU^c0AT!=YHH7So#MNSc zH*Y;;D@cjWQ<^c6AODTCpdYla@*{f>eV4W44}_<Od)CY4+%rL5Eu8*m-gI~i#R=NM zYQJ%X8r8a$KM%>rp{7zkWnCZSQ`{HRD#ARoJqro76Jc1IR8^-q>j5X^aM;^rDH>J} ziq8#kP!aTPM{1F)aI{!Y+DtE;lM<tOMh8zb_eND?hQ<te#=S}0u$SWJ%U3bF_r5hq z9QMz)hnH}WG9bX?MiEcm*Wxq%mika?ZeoVEz_%XOJ365ow*BEvow{5zyg6}?CRwYq z%r=c-w^ku*By{haWfIO@L7;)J^N#B;s+KQd8Nk|rkRigL34=8b)ffgc&=p~e!rJ}d zg)(RaUT=>*2BTZP8N2y_h}Be~k^NTS@8@^|6fD|e&$LOnG(24(KQb-h$+_B(W`sek z1zx4i1wLC_pNpxeMu@1K>GF4%;uVYpL>t!p=YwjtAFSJ`-j8<^CbQ;N*tKI8EqN2# zljD{AIBIg9PCH-^;c+1QWv;f97uS?eb5tBb>RQ^UJh@QCJIignXCtV|VliRa|108K z*zCS(88U~O5gn;N{E3M%J59n{Ves+x=IV$5{Ax?jOKguM$1-J<Ik@G+_+s7<wrFVP zNi&;N=&9jgr-)kFs&{7KSZKq#GAlteJ*1&feUq<yg5)b}>3+VOar6hym$qaOx_Tek zSzb#%f|0dxy7lNz({nKYV`pY%+mK>g37JI7<q#p)*Iu#YGr!tY9`^P3_hb5ra}3wI zye_M?eT@vTOz{70-AN`4Q<`iq5&$CbtR)Qx_-jnT5o7jCK0>Y7-KT7uWQ4(zi?R1R zI2H$?gZCOfXC2ABZ)bdGI*W8>cQ558t|?70Cn@S7T1}KteN!Vv=eoj_mTx9rsH%t} zd1NAf@$`M+?2*<&dUv$Jb|jLyx0rmF451hCd1`8j52j*MLgZ7~)@??!Lso@~eJMRW zX%-nWJ#XCR4Q|{nUi;x{#4>D7(w_y}Q-7N+*Aj%8UY|<?HTgs&nHpLx%3U+_2yL8m znb><?#n|gu)mVeplPURX6#+D*R`Tmeb~tYCj0^Mkbs$see^u7<wbr51hMrxv6Bkwg zWNJZ0m-zNmW(YK^!Rq93Zr)Bk#F}Qu2ofil-*^Bv&2Hr;MAw-5)zq7bp&oSOs^v&K zheMe?JA@iXq(L*t670iJrz{23n?7OI8*lm`S=!E2>r{e*4*C0^$`vkL-XCK)Abx6O z2>CCr7G4nC8}ll7@3$x%Y@oteRFtx22zgSgNHhXcwA*QL+Df80NQ=@;MVO`zc22fx zvTbzWWN)>sAer}tB97al6W`f0!Ks8>Q*1(Q&xE9Xjq;L`P51`URZFp^qa#axuZv8Y zJ`D+SYg#hEHLNd_2#07cP8>dV^|^lO`4zU_S#G=NWF{W()z>|6et9%-luNAW*h@{W zq$(`{4StCUNN++aa7qEg;vu$k9xt#L9mcZy2WUhVJ@bir8|lc7sL{!39qjpI`6O*O zNp-|FSEx<kWS4F*>71)05q=nyb~)xenBqjpF5s+I>-&xA0cq7?z3?e@>*?p>-ZLPr zum7{2G!4jaiBxeC&_X2FS%PGI3Q=I=yNA(Z3Zo~dtOW+%a9!xPSbFPbOkDAbA``Cd zK5aTD->8XCNmgX+>B;jsq4^aJSSM?CY=xPY;%rFT_s^Gj3^dZ}JG1Dr-_#0G6!ko; z$~v`AU3iD8e9=-nQ%6-{w2pKf5)<{CodGz{*?mxgOo`9z8jZD6XUAAhwR7})zo{+I zH5Jr!(FLzT=|^QGc(bHQW)6}drSkJo^8xN3a;SP>2~P4T?0kog;+@)YVn`70j>H>} z9+Nr;GMp!*XC>8e$669nQpCPEUxI2>9A5zvQQhTpX4y@I!hamuR%-op0>Rt1_TWIX z;Ei{ov+W{ha`UQaVVGd5QLORur9%a>U1uV}q%($m$Xu$hPbW7cO&Wete|tSzRCMT@ z$?{pT%)ZYkLwsDp-61s`7Zlm@cF{IX=SliA9S}xX6~a6oYYT6eMv^*$J_`Pme5cj3 z71LZzT^?AE(B)Ta3swxK*<9$Tr`wbV2uYqj9F`vLNZ*|&y@H@J-NIGbB-^c;zWLNj z(K8(}>GxA2qHM0&+^+CqF(l2n4p;;uTtIuTb{Xz76l)}I;H)~n@*8k)g)Z3C=&|Tp zW9=2Qqh3nVJ?5^;!lTh#We*8>)~WCK_8V<tWdcH6hIBMd3HZ6cOUd+tC}h@XCyS9r zT~%H5xRKrDJc|SR!gqLABe6y-!DJe2$Lly&y9ACs$?R^hBAWD$3}eeEReBA~lX<Km zm1@5B!K~Ty>u>6n54t*Zy1-U_k*wuv1DN&{!g#=d@TNSm+5;B$x*}Ljxt6=tRY^2{ zQGotbbf!tM!j-r*EhSG&ILM;xmqh(oL%W&w10UQ%PMb<}>We?xotqttE5N6>nxXN| z?sB|S5_fjjt<iL+B>sinbz6ved;WJm)W4AOw|$K6VB`M>_MI1aD}eqD>|6arVdEo` G4EP_Gk%Ha; literal 0 HcmV?d00001 -- GitLab From 5be984c2fb45778f3afead774da4e7a6b43dc9c4 Mon Sep 17 00:00:00 2001 From: ibelyaev <Ivan.Belyaev@cern.ch> Date: Wed, 11 Jan 2017 15:50:04 +0100 Subject: [PATCH 4/4] add original ROOT file for momentum scaling for 2012 --- .../BFitResults-fine-2012.root | Bin 0 -> 13686 bytes 1 file changed, 0 insertions(+), 0 deletions(-) create mode 100644 Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/BFitResults-fine-2012.root diff --git a/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/BFitResults-fine-2012.root b/Phys/DaVinciTrackScaling/python/DaVinciTrackScaling/BFitResults-fine-2012.root new file mode 100644 index 0000000000000000000000000000000000000000..29b3d84a9bd6ba7b0458eb24ed6def14826b55dc GIT binary patch literal 13686 zcmbW8bx>U2w&x)@!9pNF2pU|2ySuw<(BRTE+CXp(1P|6oaCdhI?(Wc7<L)~Aa_8MQ zHUG@q*HydD>N;obb$Xw!&sqDs_i}b{aDjmtafg9{F@u5OaD8n=ueQ!>gMDp2@c%N% zI$&Vlt-`<}PMpsg%Dmf1!b4QXDQP=-J^TOlX@dEuBO7V`<!TtEzdgT3gMmR%m$b4q zClZyka#1x00Bv0Wj22e*=8P^bj4tK?7bfS|yTft-x;O$|NL7`T)JZM=^UGln{)zfG zG7PNmYfzup_BT2V3>*0W^^*|l|G9Haz+Zpnf9@QX_#c0Gb!k>{B%mXsqb(5dUq@Ye zWM2dr049t(?1+}H{EP-p`w3oG+#Nm8p2d9X)2PblKO+kl%jV&>qpW&E!JPWWIJ9@f z=qh*?40W?<+v4-qfunnL7NS8yi%AgOV&UYb)>`W~x~OIw#3!9%oujQcY0!g*%#rp( z*5kWXz!*=)#X;+dhpPN4mA<O`+dx7$EE&~Mc1VYiJ`-6ib$p^(B2LnM@`iU+n-l>^ zQ8l)V7<7)BzabW+T&##XL8<EKD=zo>Z7GV9j*m`^5vdp0@ht>@CTS>%2*cAG5U`rx zRr(MQdlktwq+%Kiw2Cnl#_@<twtUk5tyu~j?bpuVU6zJ!y3FR}eE?p);e}g)8|bF} zjA@41@AHW>bzB|wYzANC*nONdt`x%*FLZ-8$OEM_W=c{lR2$0`da`7qqL7OvVajdr zoT+&5LFGf1(Grm*VCARBrvAb_<HHZqsfL3+#;C!Ck%N1>(uurWW`%{8*^jI7=*{-W z_2dIO7hZ!*di6>cik=e!L;D)4)v2Y&)5}w-RY`hyde@3ex(4P%Am1!Dj`}xD`a0M9 zLsrz8YPI?U?{+;jSthf5DK19_#}-Ljd`=&k*@{U$A4jh?Ff}_1uDfW7WS=hSCxb{y zklNzF0x}mqz8fk$khj|%81Y3B7dKgggMTJYgL|V^PVXYipIU#Ei*QJWLqu0P*+P7F zeYhl}<K4E|<<}6d_>`i2UUnq8J&hB5l}_iJ_NZrKg366B-yFmUcvXzEj{s(!-1UZ; z4p<bAVvpNuaFbbtZoy|dqlyGc+o4qJfMwg8?v>enH=k!@jGsmib4)g4f_1RJ9P%}O zQQX-{$b+BEy{{u~QEm@O4k`tDJ)+MIMXuJpt|RUk`arYLhJX5a+=Ns=ITyigR(tQT z30~X{Kd}=(qinvb$S4`RGqvtB=h*=)NA_9P$(>R5tmz*xc-zkOXlFPRH`vQ3k6j;v z9t1K9Pawbc;^0Ab7nD~-OVBpsL1(D~l77QOzz`@A$~F|!M)Z=Ep6iH)g4x#KR9N4t z%T2-EbOP~&G(Oh9okco_+rK(#$lO7f`t2x7RuuU5CjaV@v?ecwu=f}1dVX2j#mGT$ z=I^{bA0G(ecOytb=>sZ6Rc~o~JUiF%^Rg#0W5-BZuuk~upx$?Qt8w}nxTg0v<pki7 zr8&Xsu1&2d8L~Y5lPe86!Qyw~@>iu+Hd1aYyipty&)-vw$*+~u`x!XN@BI8(DXUdz zsCK;(Q<^hU8a%?kc(aX7hqB<MtlM&{bRb-byiB(ZHllvj3zA-laEWl-02klzyojbW zIolX}srU)4B!LsI%y^k_BBRwJYn7>B+0x-#rKy{$NsL<lw8C@o?=^fGe@TO{aZO)+ zK}wbHJno#!cClQaPOcXY_x9MA3FoX`%tc!Jv29RH7G6&P0bq5$uOcm4&26hU4!ukZ z^6V5~19WAo4q1K2ny}He%NH<9HoN>P+-WJEWBg$4qjQ6YRwcl8CoJ8Zgm)3*m@^JH zb}=b)ef*I>bT_FD^O^m5;S7=z_P*@riq9s}5ZoQQ7zLNKns{Waugc(5MQ5BMn_nkU z+}fFE*M0mJetXj65eM*lD%F_kPkqdhov(vJLu;%)GN{VD$Zq18kpS-Y>)TRS(?08L zkU=iL5&LcB8g^^Z{<Z<XUl)4Hr<X6@oefx7g2*xHg6ErwdvJV$6)#$}(E{;OvkxQF zQ@nE14VAptsey-YphB`|=ZzWnSyR(9QDaXd-k1zR!dDnEFZv|MA%f>Ob)`yTJ<$6n z1>gbuQ}Xu=lYLrh<+j}m$~bueGgfcnh775K;~V*QS*{sId2gqr4I-*mST*Y6C5q9g zmx|FB7`}x<rq>tC$X{%`{>3)SKiEc``8T$aU$JdxWe)`W+js?c90H8%zrby08!v&i zdayEc%?=6+B*|eHMe9)yKn`AKdD|cIE-aMA8nrU8ya;5-Qftg$4QCxMjP>okJVG{p zZ472E{Dxixs)?L6890<QclJ$!u$7fu<)4KWk0?(&)7zjWmkc*c%bn!pl^c7-Vz69- zQSmqkcikgRc3L)<hUDe!ls+N8#@1=ke>iA}yb@_vy6R)epRq8>Ip3%ZQM%54TCHT2 zf`lT~0sb|sTAiL}^)Gy%`Pj0x`Lz+^&uf=PSbo<_ia|ss#=Es-=kSqslHGbMBrs2Q z9r>hxnhTI!`1EuA?%-Q@X=$xmr94nFvm<N&Zg)hg6L%Q5J7t)!v=QCoi)7zWQIN1W zArV`?Y#S?M7<3`$-y$1aqUMZUQ>{#wC=NR10*TCIYtmImWCH9tzj?_-N;R21xw74j z?O)p0ZS+vuq?1>4yBd|}xhUNKq(^F5J1Hyh@Ht#EG#?=J{$LLpTT1mVDTNdds4v^h zZYEqU6*`v@x$KW0tBhTyNN}8#%YOo$`xHGqZ$_lL<xS+;)H1iS`M_~pj})e;PbZ~J zdGZtOwnynwgWP`^s4fk^tJ&z{EtMSevG1^yS#+!D|7QJix%0X2p>)!3%5=?hw=#e? zrnSPxTQ35V(NX!qHS{J?aFR21F&zIxB)=M$G3cY6koyj>l{%@8VQt}g#h>7E&%j7> z^dyBu?67Z1%-58r<=5KuuZ&pq2hPf&TQk!ZlF|%(9RdN7R4&ey2T{BDis1I+it#&2 z(!-AE+v^5*-m`9u2G{-F^LI-bKmJs<DLJ+~9FTk7XgHkP2$g+4gd+<f-bdn@&ER`b zeC)?-(A~#)aTvUVAq#o=#KU|jh9OS0)zY*r^}HS;@VFzJ-igU`TS?*&^61vk<hp^x zj0F^5*6w!T@9lfbx0_Yw?AyHE*Xk#DOt>LRtZhmH`uQ4Q6>NXJ#>r!J+D&Z3Q$PQp zoPLyXyA-)_88k$3?b9{ZHkD7OMWTZ}6#|}S5q#0Bn;mr1K)8!Y4Y{l13P^=Fx&jg_ z3B;V8r|szw=<FP<X?^9zShG*aOwcWT<d3d?3aCoVbJ)!u%IwyNi^5T|$qIjbWV{*n zj>*n)G~7#Rbg3yh)b2i2ak-Y=qq-Kxqr$!7a%b_Tfiq>#8G#F#f~=wXS0^fZn}NQV zyqS`{a-LW!<-1*|AuV!qm(}43A*}I`ABZ$d#>nFloPqzq<VZ)L1E~5jDsLXRIQlfd z@?+A)wjhvw*}<r7-?;)a{<9zMV-4lUJ%c2!er1M1je%l(Xdyh9r1ql2SLF&^@x{h! zCr#R#8#rHqr0P*&(kHH=6#(RQ$_qCMJ-hA*eIQQ5!ykBGZfo)~gQ|TX%J+gvRZf-b zu4}N`ajRU=3mOPY8aJPEaui$aQoH|2)ufx$_rZ5!u9qSAm!RYB5tr_s*q+|QB6B{G zBBoEPk;3m7wIB|xR`5;S3K)&AsJ3R$MKNW<Ix1K!{Ew63CB6Kwh<cUW9-W<y)8)Qi z=kmosy(DM?x4`Ke6rX1o59MiBUWv<)EhpV`LCS}#Sl_rB#ql_*g}4^e+v*w<Y!Y`r zmkDo>YU!z=BO~#QoDVgH&t>hlGgE2o39wME^frdORV>X#_AG9d3%i`u-EVsRQLTFV zA!S`7?DxpRfp}7gW|x=FdDE7={qlQg^2b~s?)o?;`_ngtBIlAs``aykKNQs-4p}g! zVH+)vjs;aBPV43L75vdvqc5uL0e9c=q}xNEoR-;ceur(XV(l;PkLo>`kz83Z*I^XQ zxCVW-PP`BncyYZ^L>e5XdHvT4m5ndH-=CqmQ5=;*>QSniCYw0h7sqSI;5AG?Jf9VD zE8}KPKb&1}ow>!t^Qs5E!hH5G%su|Xob4Z&qh<da<|t+kZuWn(@c*`6ah`$z<M9{g zn*TcwM?wRS8tF7vnJuK0sU&~3x#ZUx>DXo{shugvLGl0i@kAxCZWn22gh#3ZD>0&w z-az8H|B=Vx>oy$&`a(5HN;HW-xixvPP{%)6IJFN+$(@ob!3}7v5egcGrkw<x#F#j? zok~MesaqP9a4S_q;UQfl`ov@#8<x!O99&M2$H$HY>dP=T!QGIt%EcuBS_Y`A^hlZF zvM5$GP=f!M3*UjAEIX98pQaxLnSG76MYO?l7}B#fLzdy|(tqP<IpU0-nQ>-DE{UcJ z+#Gom`;Bmz^U6EOFu}V1N@YoUP5hpB@<%~Ik9PHUmrMYu@4Us6-|6>LuWbwD&$(EH z_cZ6VHjJ!_N7W6&*tPC;61Y-1c7a$i<Y`;w#+Ah}s-bCu(s+jSRsaQd@-)tdCtI*) z^6FSj_oGG_&}`wry@9upqqOw)y0mD9FPpsqP=+e>Y<YfY>7aInDR{lomtPbwYPXdy zwm!trY9W&DwO(^HwuU4XYwYAD^Wl?#rekLF=G$*ZjVw;%Q}*0!h$p)D_`4u4ZTq__ zf^p#Ps-EC7ZlE0E30#AIY<j(5mSod;tzQ1+a%H*4Y!2n*L-nXvaiqlbVau5L@_Yp2 zMUtXqnMUia;hdDV8opv7aT_=JkoR-YEgquH4SNe^J&wd-@dvxZyX#!&dK@7aQk27B zaP%NCdJDn50Qp@S<z~g&H43`j=X?s@h#x&{$Ajkpw=p{u*iA>!%mbe$X>O=chL*1W zI-rgXW3K*gWl{{8C?R80J8@WgHMQoN!1&p`+r+g=lli-OdGR!t4~SWoDNB%Nb3VfZ zWU=`$G6r<~wAm_lk{`Mc1&r!RUsaQ2RF5`X3cjA=CrF}%Zu2XG>Z(clDz}~|yg~=Z zJQ`mP$69nO`A5fM56sfn2p`x+JuNMYrH;AI;C41!q--t0kwbU4X`Az-H*T9vldw`4 zb38iY$JH0l<+)ye&KMR(LCKg0X-OW+hHdquk8=u1T|kkZ^$^0O=f#5{G5_K{*0OwX zzdPl{fU%^q@3X<f)i7j{VOcelvMb=hj7<pneIYbhFc0mPEbV-1wuWh5OQmwZti%KB z{4T;K(Ne^~?#M6d3{%9$as*ZU`*r;OJ6*6WNyHoTe(Ec@$-OHZG#2Ka5qC3LdExhQ zqG*+eBNM`)+0UW)`x%DKczeAmym?S~e^38C=fr(U(=y-yvc-pQ8F+!5imVSK$ZokD zS_<m%Grm(%s1s08tm6<uBK`_G9Q;L1eqE@xGQ~zxcG~CV1|_zKf?{X%^c2G3p#e0C z+xNXcnuTETc(2RUFKfrC=`7#y+$&?8QmtHcd^&+&R%T0hT4FHHE*^dsnkQ$Im{;^J z8h$h~gr~N?T70XJ!<G8Zm@ADXLq!NwN*z=49<|)`>a^<%rnR!eQ#?n0ktSQ!@Rrq# zPSv6kck{QpI|-$a+^{4yPLBAUEO4JDHqW{^#?2)Y0pn&CNtG3T8&;=kC2Zwi(&s)e z)orMZK8_cB)e0e>VmjO5UFRxFx)p)j%P(YBC$X)(2-;{<KXBnL)^X3c?>_%A=sVsp z5-7=68j!W5iPl$$87cONyXTm)(PZa6UXoAvGQ`hc?7$%O=V{6KDf`?_#fYx})3UG} zigMz=O}0|UGOlow8MITlIgPpF9e#LJw{}|Gs^vJ{)>rUIbbp<rF7|eydRpt4d1PiX zw_H!J(q-zQU7Y_VuML=D^9$8iMW{h&xOiwoay$QqsPV=t)c5xY67`W+U_`u;qi+G` z%5*+48^E%Y^1eB%t<z%(X|M7LxT8g;SI;Z&ar><DNxufl+Pf21`N61HmtwEhBR6Dp z(~jiUTLyDG`h>4NZp=6L1@4}_5ABue&3~y5`b%|=f2fZASB3hktf2m-`rjn||1e*< z?u`Hg`oHv*2<72406GD0K45y>1rQ8nqIYDAaJH+R`=E%$&<QeQK151jc(EGZQPpM= zr7tnd)hsAw?$$M5$FHb#&aGgF9}I~`!o)dT_8RlpxPhekPM9xWv_ioTcWHp*yK@0) z77o%^c_ku{5Ni95I!O9G#$zW0n|!r7>vlW)awJHC`L0hOTcf2;R1OrSc3l>6yM?Qm z8^W3;S&NgcMzIb`s*}VuK~*B<Qr^WY08=lEa`#^D4B|fUwIl=hS}+3%>2bb$_NFJe zlMh(yWh4C*Us0dcX;&{&F;8dboMq7^$+{BJwNcRtvb=)Pg*$><=w>m+3`7bPib^G` zi5CTNi2N*Hw*Psi8PSOl+aPOfAn&+U!$T&{j@*1PS2?dPM=K72S4>(hn`%<VA&2KU zR-n?Vzg#gJ%Cnz0N7IQ~0#A^thuQKN|4`I>`DXKCskv33e8S+oYK0DRDS{mB>6NVK z&4rBR)atI>S>tG1G9U4!fw!s^#~wD@7D(zkkc@gPhudi~uheSaVW$7UZ9{+{?uyuD z3cPuWPB*B8Brjk|PWQ}QeR+JU%ZO$&Li(*8+lG*oD>gxCAvd>inUlZFmW$NN?Y`}w zri%ZK1f592sJQQruQ0m8<|uq8w^m8)NXfuqz+uQ3(*k1(2}zc(df%eu<7c6Tr_h&? zC^i+ZI#YomZ!bhVyZdH1cLmpLW5{h)*{YXju9QcStxfM&ou}up)?I6yq`8){=V@CZ z{_D*h%Oh8o#}LEns<qoIf%S`%$&`zZyRT0pTgQ9glTQy3eos(9q#SGT8u%vN1#`GM z;5aLC!$LrAwwmaI72^Rw?^9>`iO+|2<?ZI;%>w4C<v@dr63dyd75hXB3llWsV$l5) ziCi=I4nJg$!t^$JY5&JjqmS(k5)RU-VEJ8|_sEqrYp~F*lht}qjsNaL4=7{pskR*| zEfOE%WL+qTp^{z&o_lPo^tDxTcYT)8lY=yn^WQ>xv)Zn$jpNq2ebM`x`T+IY0y?pJ z*@7=onlnJ|ws(rOJ&KU7JCPM3WN+0g^!F|Vjzw$-F5Zgpx^(n5J5V9aTNxUMlKzmo zyu9o2*zIv5`{fAaQzgw621@3;581*Ckc55M6@r^$T=13NoM6r1zoZp{lPgVf!Qv0% z6qaRH){<`dyWF!l_j2k;Y_^y)VS^^7v3jw1%A8{;xjo~K;2_62HGrv{L`^cAl6p7_ z5*QlDJ<)-{LG-+}dfphG(;!cwp1Dku#$7&pkw(}bOTt9Ap!*A|G$KFJubmbShu#-L zNA>aB`|f1fV8eK~hNM-m7s(seG{;=AW(~GXC9Q)y-!B7!=kfPxtIm7mv=+PV!(XyJ z10;uOOFuWE;0J#35tQSrCJPJ-F;Yj~fW+TZ8+TBN4KfsCbWJN7vXU5k^-mV=)xm8_ zMIu}t>}*l!&yL6o2I49xlE5cuwFQlvKAIZtk+%0r3yJy|O?|yxFQ?*2IeN@ogvC|j zIzt@DV5){)!I>k41Oc0NghM}uhjs?NTQM^0s<38bH40QY&v@Ep^I5tp)jY3H9yK5( zj*31<oTFJ2*wX4@4>)CBxJhx%@W}n*5$FF%!U{FL9ZUAw_WhiF^g7Q?75W9z##2wo znFhQ^6(M}lXi6eG54gY@DX*-*5F`8p96c02$%;6?YlaPd)&qdgYXXi7Gxj&=W}q+- z{k#@uVwAB_<`tKy3m$PGA?t})&GR$+8QdD)9fFv;QMWptez~_g-Xrl(Z#1D+-!GCt z#gl|cHqfN$z8Y%R=6##>Iqi!6321dX*TvC0<)OZy?eH<K_8F5}xP|T&-I4#dx`NO4 z?}QxL%IsgN_`i);Y^(m4x}y14T>&>Ays9g27O_mx*hNEoGDR2y*IDR8OJzefe7{X% z$!-0nQ}v=Q#H}qr&t!p@$o+wW8`+g1`B~*>Z~cegm329&1Af_%K`k4-Z+7Lm!kCtB zmSsm?ZXcgRmt3GL3zivoi#Mw(U&_!xsBel%19lL87Q#^Zg*<<|g{i8N1=FRg+#kNp z>VYGxLi1O?Oe?Hx!|+7?u60m}eOaIZqtc&uCmBG6d*;{P($GQyp94lEfB`owqzO4# zQedTFqb$F-SuR1jWIGzpe?P|;Vv+2}Aj75TG8#SOIn9m4(7Vug!0g_dIwz?<H?G6> zbFTlP-+kHgtRK|>_}lhvPz>t&1&&CYd#;EDPiZ2a&Lo|li%e(b)HS>?Idc8QXZ<%E zs<81nj!he~>6}U8Y$nkV1-83I4V4teq^7SC2;rY<)2G%^;N|ubxx0~~6_*{G=Li<r zp;;$F0??ihZS{hHeU^O7hLa%)5)K}skGKrLs0`2ZRRsXCi|*^(9Q0WAhQXIh|HXxJ zeQY}fMd$m(V^doAH(1xl$>xWRWYxpgi=Csj61y(^PD29TU@;@+k#3*>e!U^hSZQg` zZQglIOz<2a-dAJ6vT@D3pv*&#@5>M6+TV|@MX$=okYZ=CP-_53^%Adx4maFd!T3a3 zH>2E6?beYeIa&Scdil+vegL^p1&^X`grPga80gseEUxq>WLC;YZhqlbY{_Q7co4q) zv>|t12AaQ%(C0(3kkk;=-jxtCcyr_MIZ63)FK0{PF$u7NbV0p<_p+6iDm3h2Gy;tE zhTyfXi!*bF(WgIsjN(dv{Z89LF2Q^)!E6IjuxwAf7<vffmlKpgE$^Y8EXDJxVs0r> z8hl}z^NzbHYs%UUI$66$VCF6ZymnO#;qBKw{xi#dlD>`rO>`d3bz5q;wm|@solb?M z2EV4%A`1S{9`ldU0{ShkcfTPU70-`ey1+ij0-Pa2#XAr?^Oq;MUvW(i(i7&!YV0rT zt6dlKx(2PLi{X}Eu#q)6><#@wQc^K*Mp7b<DtqmEzOUQ?xvV_OzVV{O#(bZ|+!`(L z{e8+zujjlb>MPBc2h*!}GBOLiEYZXIg{8l|EP8ye@F_x0^y^9|T9epStaEJ;PDoOc z&9LgYtocruI9F`S1UfDuy5h1CIxUiuG;a(1DI`kzo2casj-81DrI`b1BX$eBknzrd zsl0;rrG7(+nFkHmFFwGIK_fddH?OiZ2XdSlTaq7a?3~5ciqp+t<O2J|AiQMqbDw5M zS{cX1i?QlL6I_3}TH9!mjMv*{oRTkao0@bdtD>y|V#{<r4LqMJTH&sjNKdhWr4se% z>pwqK!hK7mco?y~p**aBu8<X@t&(hw9aHW^==fG-cp%M38}6J;9lL5mzB>s?`Te0a zUY0ob&vX}Q7DK<H^%}J!QJSYBQ2Jt|R87q29=mNd+PALlh0nG_y-C{YGi^t?v2%-i zdiVBlOy?kfsqM1%o^e)zunzF(TeW7{h+#Z<iX3iuK4adM7PWGyw7jTBtOH0#WKfEI z`kXj-mw4eFc%jhFdT+svo)`BxpJjw=U}Pmy8L@e=4saqQ8EiH7a!=VjsRFt#^>IeA z4e}q&5`!W4PlG;=Mkms(`#koAa}s60=Mx^L>~tYe7VDGJrMB~-*3E{82JqAHTHiz8 zk<4TIo#J`x5wt1I^3mv%1NFrv>vQz<6^X~oJjR+^g4IPPE00i;hp)SBo75`F1Mw#5 zP(!1l$QFoOO0Dl?w9mKpnXnt{=P(QmJnTO#Mmzm)afGP--z-M?hsA%h=>N7}Y5WHc z<`nI9k@AnGg0#8{PTSOQ!Dp4rYBUsUmB8Wxg}p|soZ?JZKgKsGXbj!{ljISBoY+BT z81S+0YNaduBXjhUrE^r^buJIo6AQC`rC7|io}`@+IHV0K7=9h)gMu%-q0sx)8Eg9t zRGS2_^1FwnlF1rXIgS431zb(K&zkDmu?G4zHb#5XaaRMUppp!4MZhYA_#wV<hnr7> zrhdOh)qu{z7tdpL*6f{z%ggnXMDAcpQgNtimW%_Nn}N5b#c0*&@Kv{)5;x3^t(&u5 z5C4T(ZM}3(?MTSTaygLRUU$B48SANNH%a_P+-|woGeVqtf%+_;@KDKL%KwFsi<G}G z!8$FkAC8TYwVU5WAKR)H-&W?2%;cx6cq8#IceD%*o~ra3><W|SYzm8fR2VUc{L45( z5Tz#4Kbjq&2MRt>#JYnw_a%kutzQ;D6|Ef&ZQ6zXzV}3!G24lsk2_q->YJnTOtS4) zToh`wGU!W*&_iJj1>F4{5PDo!uTUl-Y0x-?F2#^I2Jk)xJgkSj1nvLGmEtv`IQxuI zwQ)JZR)-LHYie%bfgXXWnmql(yFZi2sX|OhkFClY+BZt!vxOPPy?&B#ha+cX971w< z8bE~o$lC6^BgOC73FkBaHdkmKvp37v>3GfP4iepx1&8-I9)sl4BT^pQmb)*sW60<o zEOa>TM_BV&==LkgbBDl~w`>u~!|lr(jLXb;UK%FoJ5UFYmXCQ=t$XUriyAew(b4nS zSnx$*PL@$w<>711khenqY=%)M*k~pG=1?#hygFEtH{!ihDJ_VxJ#A|A(cIu}dhS>g z34O57y~nF~XPoG)edkB+a?QYjW*Z={yA=eu{ssMHS+c=USaijVS)anmXXebz*!N<1 z+x#oKZ9z|_GyL8Dhjy=FrC1~mSwnwQY&0|CD&L<OO>XWHh7FvTnv-%6FSUv;mI<5G zAO(aZjQ#B)H1sfW{`-l;pi=h!P1?=WV4O!1eV%ScpLB24F}wk}bIzm4gE=0Cr{t~s z0=@HiK~v&}w3w{~b)g-e;frh)#7vgt8~BU;P(sTAk9(_Xp!Z%<%LAoi84>$dZ{pl~ zMpFH-G>RaX%;V_`0}yS1=i^Tv%yWKVs0UE{{N4`~<RCxVr^VM(lE)}0<k&|b6!wy^ zULz|o7I>4Mh9m@r=HtqJeE%Y|_Yim*1xSD9c^ZX^OXZcXwNp?pvsB>O_p9t((9Grk z_K#zFk~S42N@x?mUEarc8I=p^5$fgcytkBlLAYT*wfVa(;}rn>yMF;7{Lei4pUoOD zW(ex?Rsa_?H5X^|uXg6nGWHe@*y0XA6I%~bTPu4TbF<e2z~z+&3GZNb`Cn-;qi19L zX$~mJw{>1yvl@F9`uj{Sk4A^w3JSS7{mLn8ruspP#Y!i7VTHgsXNADpB2k}CvBgP^ z!W%}!|ASQoJflCUMCu2c@efzbNXkEvuyQ1P5MU>T>;gKhU_ZBG=hH5pr13TLRj*RI z5iEXZW*KV*w0n&T+_W6$iYPr(=h8~he9TbBAdL4sdQq-=o8^#9IV1#KB!(PmS#^-| ze&0VA!KAw}hbsRTe!1vnBXpO~njJbn0v-3}oTL&yuAuK*f&6sc#aTiJQGU^~RootN z?04Q_Uu}TMF=jI)iZun7BT^I%+qs1@W&}LM#MQVt8C%uj-i~D>-b;UdgYUw7%qrF{ z{&AE9Zx?A}+$wc%xp>?q-?%Z$lrK?9JdXQ%{hQu>atzFu&@wIBEBLaA90*9o@se9H zBOTb4&G88~=uz0f_<QJdoJRIggCJcu@XIRbM}W&&-_Yi_q`m_t^%v|iD>ON8{j*eT z<bd5N6bqD-YQkbxbxx}SM}u8%LZpv1%#EX^RY3b4M{!vla9l}qM5w1z4WoC`Xd|R{ zoK&|gw$A^6W75$OyC;U*WD)G#wQf9$0z{Lq%tXo}NQCoVxp*sHLnuY9Jlyf4$Yh=I z>v|S1Y6IKPbhPe^D4U53Bk%i&Z*e^C7+1*waU-7&MG!vFqcyt-3Lf6fq)Mcdr7opi zhl7RGd(y#By~j)mB&Je1UhOlGg>ne~TZrU@Dn9I5_6$*<dnaumJ*sKC17Dx<F-bYy z@!Cr8L$zY1bo<58BY6tK&7Q(E(#@Xb$g*YAIH}INoARH*`FP?k#InvAOE||Tbx!%I z<6cYn{uq8{$R@%DQE(T#n8o2z?ezS-`qYTnygk_WX?tsd&i7ZC$4_Wqzn}fVZZjtE zZcCQqAZUH(MB6u$g{<=cF5{9H@2<1fxUErIp{Lr5(yO2===y|Lwz*g}om-2Q^#|8p zbZ%}uP^bHyZc-X+sPntg*XKq{uTF{MMU`#c*=T2O{iT80-(XqPJm)CeGMK?+t#HtN z{D77C;lQ@rJ*OI)PMY)Mr1f5FnTqa(fUdg3C(DCISr5L>AJq*lTVua-AyLt{V3q3K z(IAs{{TB=ge6c=2G^EhB(PL;BR{|)0ay3wBHN9Fj&G7sCi;YIixdGX>H&EN93?<R$ z66z+<_k)m6M2%{Rboc3p-pVjSmcx#i0Md={7~Tk}@o=f^k38Z_rmKXa6eeGCcA~H; zoRHkHm(#}LZ}|F3UwQjM100yDc5Mv5yNEUo&NS1C&Ogde!no^fLCn&T{nf+Vocn_! zvc2kv%X)HQ+LbnTrzIM3NAN?~18A6yqbWMw%XaM><Eufz)!~{b&&n_y9rvWzg2JrD zoVATsmB$Q1Yvpu>4~db=-IjjSGt1N|Ztdso-&lh7jU<6!QQu{L-U8<H+py5^@9Y5u z6zcp5IVQbmBPwfEW0_OtTs5ekunw0o{d&#c7Meb@FSJ2}sVD4(qF_JS$?>+E-2Y~9 zS}2NbxKOxUpL#Cg?q$2YN{XyMMe|g}%;vR$dZp}5DEr-QW_J<k+wFJiSV2QJFE0x< zqvdDwA_i7(IH}RAZzCLT?wQ;%ZnprwqLYqoOdtjAFhRWnKF|WVIeUmPB2*wBuZ6Yd z{n!f2r_?GZN6CG|w^Ax>xeI`IPH2xx8k97?)RA`K67bQUK294Xw_?-J4A#|YEtA7@ zDQWXDY5J7tGarJ$fqE#H-2kA9mV7A+KiPBVL~>~UaXDQKH<iE|9!bmUp%j$`4PbwP zYvL2;DYHt_!;=M|okQK26Q@{~7*+L2V(+HVJx<`;`EE})o1%=7-_u#5cXf?|)h+_@ zEycN;#%Vl}_m?QGxP5^u76&wLI}PhWUWCaS4q-CTgg#_yl^glJ#LLr<+KeobOW6nf zB7cQN7q--S7q`K@v)`M1U3r4Pf$@`DvkY4!+2-Au9?uoa2m7!@J2I7fJ@^F*Xu>Nc zMZd_EK7z|a;UjyKXB;d8oUOb!jc4DEWVb0ZoyqIhr)?)Vx`Z_a2E2*Iq9Q({#2FWc z=W=qBYjMM<teU`fs+1ECm9xmD%T8)sRlV%p`%#!Ggs}N#GC5xvH#_N0xO$gJ65&Iu zmNEpR?10rY9>J4WP*7#xjQy<*`Wj+g?=@&;dLbMyyCqn#ODP{X09S#M(pza&s6_cY zKt?ax!VPd+m`Zkz)96g2>#2tG%n}9v702X4v?P4=eHV@NCv&3pIMj+d%)ZeDBDmGo z4%*Um$50cI@zwsqHR2S+cqNJKcFD&YIyUsWkG?bl-1kR~EgB4n`5y&UMx6p;&XVKu zfp-!l$-j6CGn!mneMZjTag<aOZ7K0p5Fx7ktoU<5w;gnoB-5>9K8U^8FlNY3uUO%x z0*x@N@NYRO5#7d`P~$f2o~X;Ij=Syuexa57B-96d|BKR3$20${A<i#6U;1B4lHJ-= z=<FV>ktj)y{b1H7k0p$f)6ST)%F0R=C1S<Fp5f%S(PN@Vefg>E+QTz^>j^)(gUR-e zI|7-nwdXP!p+gg-Zo;wDm^7dP@F5O=eCoH;gB$|!L38OG5yV4F0~#3m`6Wm9^|=ee z@!13?(H?T`vrtNiCfV3P@1F2my$KF&o;UO;Ux!#LqI>T?>U4!|qSUNw5_k|tT8z-R zqAwSsl-00zMy#;bs=e2@vh&B2YL;xlb>bv`8})~W@B5+5AKC;tso%wI79WNZ<$uj& z;{~T_ov6C<!UFXXipuz4nkXIEuG64N7JnT0I_fPZ8tt5~%~PM3pt`%4GD&lCSY_7S z^dir~u`5%^60KFf`W_2e`FLr@wV}ey<v!m{bgRP|29!3jjNs=Ib9X8U#T^@?TcV2P zK916}1_<IK9alpH%aoS2+>Tj1JT_|5#8ya`G!hm=`;1UR8rQ%^=1eDc;w<16<B!x% zE-7=j&ToDC-YhD`T2xdASNL0O=3p*Vn`=$lb|!~pM751~-jR3Fh1uB%9en20Bz2i$ z+}-XtUoh}O?8_te&M-rgf<-l@r>icl9GO=$U#bI*wxLPgcxI|5cz5g#C)~Eup&KEj zWozZSVbgOxsZ2>R_nnsQ+9Ll{{7DQ`(v086yy^z;wE6S=(n$T!f}*cX+L7caPANTR zpdj$%fWIrC;P5x=^3ml|B6f(1A9UEOcg14Onk%Dj&%XM*VYFH9k*8{gKp*jLSuhz? z%CQKL0g!tSuEUdy^1>P|vByoMJp0BcxC5X0eLA6i-1>|RGS2o(hezGeKw`8KF+~5g zZhRdjtf^TsMdczer6o>cXU>n-jGX|@lyZTczY6_GwR4~*ioSuWL}zEwkM~5oUR;+* zI~`M-3_H9|j?U_VAHS}>xAthmcJJzPK=!aFtHWP~PWV$-UadB``eRg``MPSi&2vC# ziIxmosNn~Mr^*I@;@((~ffN;2$7R4+mrXef)-jTarP-f9j1XIUEM2N>o79WB#FH%# z)x&PG&QE2Ga%=HW+&>}ZlGjC|DI!~P4%n5(lAdfs7^1gWs>8D=ol@VT@$!Mi=M93_ z^9yHT)zJcZ(MtT!{lBB5E18aZV>>N~`QR=aq*q9w(n5&5NX*7~<@_wA5+3kqZHvpq z=;lHui1d#{`IS>r#j;W-?}=8lQJ@y?=`mOMjFEz8Qm<3q=?B1cNXqe(TYsrz{GTMb zU?s$f;)!poROm`srw&xBZnruJmLn*}=K*TyR~o+5eiQQr51~G?QjwhditD{ZmZnR0 z>a>!rG}6pvvdp7BUMyy8=GsB$k7C-B!W~FV=T+ehqp9+>*k94->T=r_t)$$~-pSGU zMSf&~T{b`}!6f|c5^N?Romv8nupRC)Shm}wD7)wUu4{ED-fO8Hle4iWKo_gnzmYgp z{8Zi2LN`^03-{ihoKv|Xvf800!aLjUa{@V{ZVo|9L*EW(40H<x-0rdJoYb8FC(#({ zxSng+RmCSK$(2(Z&W~Nye~;e$9p^J9JN|%b0d{}^xjm1I7-~#%V|ni|k%U|cK8u)0 zM*s73(;@XKZ_<mVp`3{UMfQ(vUtfdgdN@g}$ho3y;=L0IfmNb8(g2`@>_d&iBfMA- z3Frffu#kNdv4op)))8;SqGE;8hc}$sL@VM|zr^0{XHKiuSzYE*Yb1=?ILJT=y`;xS zm9PC9Ef|gAG=zNxE=EY!ad-TB_29<6j+#dnTQPBewT!827+#z2cg9skP05rwcIy;l zgUeRZCa8gJb93pA5|~(R_u8e#0`!H+Simu5q+H2V#4Y*>iZ<~_cILQS_ZAXtjZt36 z&u~iMvMLwimM;9MKMX%J%uO#WBDbwMq6Im8+X?E1>z38+8x@?;zD+D-`E+HTk2eM# z&!VHX6`Ua%WU@q^1#Q+Q3u9d=j5fLUAhaI6DwVU~OgwmR@+fd+GC6{84lamva*Bw~ zCCi=D4lnus=xzT(Bjl^YiUcXCB!TW1sxtDB)CX?32~O2($C?by8HU?>O|=H9-<PRg zg3pAG%H2pWp_uO2k9N&6Fyv;a^#(;sa&FPin+%UU<=6wy;X0?Y8oiNqbR165Fbk<r zyM*9dYyJb*3#(K=OH^rQDU+<k7lU^Z)SD&~`G=>W+NRzsZe8Qfur#A;yQM{B)ug{{ zu!c~T!t8sh5-M$FrNKA)@U{>Y^DG;8i7U-Z)GvH)xOzd0o@pCCtHO1M4Fq2q`C*pf zvIy}yB71zo*aN<aG^ETEqn1?(E(Z&dY$uWK?Nw`R4+&AkuEslj{hH&Y^}KBZzmr)N zJb2ju(xfwJ6ZI%p;Mlj$2xRahNf0E7xX9(=+txsIa%s1e6N@O;NWjn|^BI(R6<=`^ z8+xG@dC`nf2gNPSL2e`vf-q&<TcEzUVaY|)ps36DMv7l@Cv-a-f#W=9Ws|hO<oH6u zaH_XNHfVP+`?e8iMbXk^9f(6d1k&Nm{EDmJ2sR1SN)MNHa1#vBV{E3JeVV-HnM0A1 z0)Lz$68C7^LmZ5Bb|x9XPo$fBdxcKl!$dW=&omoIiQBcA$X?0IVapu4froTkkY=MB zQM)DPN8Eq%b7wS6l2=`M?iF_S1D7&~f(=99Yta07gheMU(oHeaN)q)#e_u1>SEJ82 ze=aXNPOLdQrLghHk3YArpuk0kXy1`a<hcHJGsVaZgeRw0A4N&{K~%h3v53z@1Je@s zrkvNFgtr<BF4`Ri2V}FQUfYV)Pprv@lxM1pkRIU`3`KPkj?h(;;Ev<JybUNAc7wd= z@m?U$7qs-b`qZ%Yp=9m*>CD&m$MQujCut)(w$pbW!I)Zs<|g@mRB~tS{u0KxBKmDA zI#DTrhFqH%q-@lyVdvuQ75uZ6{2sbj|0z-w^@>>78%%j(1ExIua>g*Ak+2hfGOxrf z7P>yMlhYjeV)`y)XiE<?%y_x+P0s`+SJO1j@c<DaNR2e$O&20Ex?6ys3=J<g_{v6B z7ll&7a$zNx2BNnrhfcxByQntRA5$+XyM>24)>C)*jv~Qw@A-wZBsA}J|HL0I#c~?b zflAW3sMd{itX$Q+pKwRg8Yr`dDy<5H-ecXZ#p6akQ!i@v)0xx0S)4s|S0p1QlyANo z+c48DcV7A57FV0lVWXtFBoa&wQNO}wgBjh<|Mj`=oXshITN;x%>_cd%4+@X&Z7mQ8 z#gUbM+b@WaxH_esA(oTY{o%$x2f^Q2(GKzLD57~uo&;uX!WQwO33j2tgK)XHw+<z` z5?2m$*o*4G{6$;9x6tM@pG4%&f;-NB%v+T~(XCbh?#SesqlQdhK&d|gLHF8A;s~#l zn7OF=hbW`k_fa>)8!T$qOTDKD6Ad22gwAHm*r@{SUS~bz@`~&I*B=zCU^GcO>$Xo_ zh!+RWmGr(o4->cs@{w?U2g0<vm{>ux_FiBSeRP*!t4)0Wk#}Lg{8e}UR?fWEto&7X z{s)yaua1bXt9SoNNz3cy|JhdRSBL*}yWXqIf4&#w)#*Q5h5R2#<*(zV^Uv_Zp#Jyp X|1FPt71#bJ!~d$yak#p1#KZhAvz)UG literal 0 HcmV?d00001 -- GitLab