81 """A simple filter function for testing if we're running in MC
82 returns (bool, str) where the str contains an explanation of why the bool is False.
83 (probably worth re-allocating somehere else)"""
84 return flags.Input.isMC,
"Input file is not MC"
89 from JetMomentTools
import JetMomentToolsConfig
90 stdJetModifiers.update(
93 ClusterMoments = JetModifier(
"JetClusterMomentsTool",
"clsmoms", JetContainer = _jetname),
94 ECPSFrac = JetModifier(
"JetECPSFractionTool",
"ecpsfrac", JetContainer = _jetname),
95 Width = JetModifier(
"JetWidthTool",
"width", JetContainer = _jetname),
98 CaloEnergies = JetModifier(
"JetCaloEnergies",
"jetens",
99 prereqs=[
"mod:EMScaleMom"],
100 Calculations=[
"EMFrac",
"HECFrac",
"PSFrac",
"FracSamplingMax"], JetContainer = _jetname,
103 CaloEnergiesLargeR = JetModifier(
"JetCaloEnergies",
"jetenslargeR",
104 prereqs=[
"mod:EMScaleMom"],
105 Calculations=[
"EMFrac",
"HECFrac",
"PSFrac",
"EM3Frac",
"Tile0Frac",
"EffNClusts"], JetContainer = _jetname,
106 calcClusterBasedVars =
True,
110 CaloEnergiesClus = JetModifier(
"JetCaloEnergies",
"jetensclus",
111 prereqs=[
"mod:EMScaleMom"],
112 Calculations=[
"EMFrac",
"HECFrac",
"PSFrac",
"FracSamplingMax"], JetContainer = _jetname,
113 calcClusterBasedVars =
True),
115 CaloQuality = JetModifier(
"JetCaloQualityTool",
"caloqual",
117 Calculations = [
"LArQuality",
"N90Constituents",
"NegativeE",
"Timing",
"HECQuality",
"Centroid",
"AverageLArQF",
"BchCorrCell"],JetContainer = _jetname),
119 CaloQualityFE = JetModifier(
"JetCaloQualityToolFE",
"caloqualFE",
121 ThresholdCuts = [90],
122 Calculations = [
"LArQuality",
"NegativeE",
"Timing",
"HECQuality",
"Centroid",
"AverageLArQF",
"BchCorrCell"],JetContainer = _jetname),
124 N90 = JetModifier(
"JetCaloQualityTool",
"n90",
125 Calculations = [
"N90Constituents"],JetContainer = _jetname),
127 ConstitFourMom = JetModifier(
"JetConstitFourMomTool",
"constitfourmom_basename",
128 createfn=JetMomentToolsConfig.getConstitFourMomTool,),
129 EMScaleMom = JetModifier(
"JetEMScaleMomTool",
"emscalemom_basename",
130 createfn=JetMomentToolsConfig.getEMScaleMomTool,
131 JetContainer = _jetname),
133 JVF = JetModifier(
"JetVertexFractionTool",
"jvf",
134 createfn=JetMomentToolsConfig.getJVFTool,
136 JetContainer = _jetname),
137 JVT = JetModifier(
"JetVertexTaggerTool",
"jvt",
138 createfn=JetMomentToolsConfig.getJVTTool,
139 prereqs = [
"mod:JVF" ],JetContainer = _jetname),
140 NNJVT = JetModifier(
"JetVertexNNTagger",
"nnjvt",
141 createfn=JetMomentToolsConfig.getNNJvtTool,
142 prereqs = [
"mod:JVF" ],JetContainer = _jetname),
143 LArHVCorr = JetModifier(
"JetLArHVTool",
"larhvcorr",
144 prereqs = [
"mod:EMScaleMom"],JetContainer = _jetname),
145 OriginSetPV = JetModifier(
"JetOriginCorrectionTool",
"origin_setpv",
146 prereqs = [
"mod:JVF" ],JetContainer = _jetname, OnlyAssignPV=
True),
147 TrackMoments = JetModifier(
"JetTrackMomentsTool",
"trkmoms",
148 createfn=JetMomentToolsConfig.getTrackMomentsTool,
151 TrackSumMoments = JetModifier(
"JetTrackSumMomentsTool",
"trksummoms",
152 createfn=JetMomentToolsConfig.getTrackSumMomentsTool,
154 Charge = JetModifier(
"JetChargeTool",
"jetcharge",
155 prereqs = [
"ghost:Track" ]),
157 QGTagging = JetModifier(
"JetQGTaggerVariableTool",
"qgtagging",
158 createfn=JetMomentToolsConfig.getQGTaggingTool,
159 prereqs =
lambda _,jetdef :
161 ([
"mod:JetPtAssociation"]
if not isMC(jetdef._cflags)
else []),
162 JetContainer = _jetname),
164 fJVT = JetModifier(
"JetForwardPFlowJvtTool",
"fJVT",
165 createfn=JetMomentToolsConfig.getPFlowfJVTTool,
167 JetContainer = _jetname),
169 bJVT = JetModifier(
"JetBalancePFlowJvtTool",
"bJVT",
170 createfn=JetMomentToolsConfig.getPFlowbJVTTool,
172 JetContainer = _jetname),
174 ConstitFrac = JetModifier(
"JetConstituentFrac",
"constitFrac",
175 JetContainer = _jetname,
178 groomMRatio = JetModifier(
"JetGroomMRatio",
"groomMRatio",
179 JetContainer = _jetname,
182 JetPtAssociation = JetModifier(
"JetPtAssociationTool",
"jetPtAssociation",
184 createfn=JetMomentToolsConfig.getJetPtAssociationTool,
185 prereqs=[
"ghost:Truth"],
186 JetContainer = _jetname