Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-04-06 12:27:53

0001 import FWCore.ParameterSet.Config as cms
0002 
0003 '''
0004 
0005 Configuration for Pi Zero producer plugins.
0006 
0007 Author: Evan K. Friis, UC Davis
0008 
0009 
0010 '''
0011 from RecoTauTag.RecoTau.PFRecoTauQualityCuts_cfi import PFTauQualityCuts
0012 
0013 # Produce a PiZero candidate for each photon - the "trivial" case
0014 allSinglePhotons = cms.PSet(
0015     name = cms.string("1"),
0016     plugin = cms.string("RecoTauPiZeroTrivialPlugin"),
0017     qualityCuts = PFTauQualityCuts,
0018 )
0019 
0020 # Produce a PiZero candidate for each possible photon pair
0021 combinatoricPhotonPairs = cms.PSet(
0022     name = cms.string("2"),
0023     plugin = cms.string("RecoTauPiZeroCombinatoricPlugin"),
0024     qualityCuts = PFTauQualityCuts,
0025     # Determine the maximum number of PiZeros to use. -1 for all
0026     maxInputGammas = cms.uint32(10),
0027     # Mass constraints taken care of during cleaning.
0028     minMass = cms.double(0.0),
0029     maxMass = cms.double(-1.0),
0030     choose = cms.uint32(2),
0031 )
0032 
0033 # Produce a "strips" of photons
0034 strips = cms.PSet(
0035     name = cms.string("s"),
0036     plugin = cms.string("RecoTauPiZeroStripPlugin"),
0037     qualityCuts = PFTauQualityCuts,
0038     # Clusterize photons and electrons (PF numbering)
0039     stripCandidatesParticleIds  = cms.vint32(2, 4),
0040     stripEtaAssociationDistance = cms.double(0.05),
0041     stripPhiAssociationDistance = cms.double(0.2),
0042     makeCombinatoricStrips = cms.bool(False),
0043     verbosity = cms.int32(0)
0044 )
0045 
0046 comboStrips = strips.clone(
0047     name = "cs",
0048     plugin = "RecoTauPiZeroStripPlugin",
0049     makeCombinatoricStrips = True,
0050     maxInputStrips = cms.int32(5),
0051     stripMassWhenCombining = cms.double(0.0), # assume photon like
0052 )
0053 
0054 # Produce a "strips" of photons
0055 # with no track quality cuts applied to PFElectrons
0056 modStrips = strips.clone(
0057     plugin = 'RecoTauPiZeroStripPlugin2',
0058     applyElecTrackQcuts = cms.bool(False),
0059     minGammaEtStripSeed = cms.double(1.0),
0060     minGammaEtStripAdd = cms.double(1.0),
0061     minStripEt = cms.double(1.0),
0062     updateStripAfterEachDaughter = cms.bool(False),
0063     maxStripBuildIterations = cms.int32(-1),
0064 )
0065 
0066 # Produce a "strips" of photons
0067 # with no track quality cuts applied to PFElectrons
0068 # and eta x phi size of strip increasing for low pT photons
0069 
0070 modStrips2 = strips.clone(
0071     plugin = 'RecoTauPiZeroStripPlugin3',
0072     qualityCuts = PFTauQualityCuts,
0073     applyElecTrackQcuts = cms.bool(False),
0074     stripEtaAssociationDistanceFunc = cms.PSet(
0075         function = cms.string("TMath::Min(0.15, TMath::Max(0.05, [0]*TMath::Power(pT, -[1])))"),
0076         par0 = cms.double(1.97077e-01),
0077         par1 = cms.double(6.58701e-01)
0078     ),
0079     stripPhiAssociationDistanceFunc = cms.PSet(
0080         function = cms.string("TMath::Min(0.3, TMath::Max(0.05, [0]*TMath::Power(pT, -[1])))"),
0081         par0 = cms.double(3.52476e-01),
0082         par1 = cms.double(7.07716e-01)
0083     ),
0084     makeCombinatoricStrips = False,
0085     minGammaEtStripSeed = cms.double(1.0),
0086     minGammaEtStripAdd = cms.double(1.0),
0087     minStripEt = cms.double(1.0),
0088     # CV: parametrization of strip size in eta and phi determined by Yuta Takahashi,
0089     #     chosen to contain 95% of photons from tau decays
0090     updateStripAfterEachDaughter = cms.bool(False),
0091     maxStripBuildIterations = cms.int32(-1),
0092 )