add jet track particle thinning in FTAG1

PhysicsAnalysis/DerivationFramework/DerivationFrameworkFlavourTag/python/FTAG1.py

@@ -23,7 +23,32 @@ def FTAG1KernelCfg(ConfigFlags, name='FTAG1Kernel', **kwargs):
 
     # thinning tools
     thinningTools = []
-
+    from DerivationFrameworkInDet.InDetToolsConfig import JetTrackParticleThinningCfg
+
+    #from DerivationFrameworkInDet.InDetToolsConfig import InDetTrackSelectionToolWrapperCfg
+    #DFCommonTrackSelection = acc.getPrimaryAndMerge(InDetTrackSelectionToolWrapperCfg(
+    #    ConfigFlags,
+    #    name           = "DFCommonTrackSelectionLoose",
+    #    ContainerName  = "InDetTrackParticles",
+    #    DecorationName = "DFJETM1Loose"))
+    #DFCommonTrackSelection.TrackSelectionTool.CutLevel = "Loose" 
+    #acc.addEventAlgo(CompFactory.DerivationFramework.CommonAugmentation("JETM1CommonKernel", AugmentationTools = [DFCommonTrackSelection]))
+
+    #FTAG1_thinning_expression = "InDetTrackParticles.DFJETM1Loose && ( abs(InDetTrackParticles.d0) < 5.0*mm ) && ( abs(DFCommonInDetTrackZ0AtPV*sin(InDetTrackParticles.theta)) < 5.0*mm )"
+    FTAG1_thinning_expression = "( abs(InDetTrackParticles.d0) < 5.0*mm ) && ( abs(DFCommonInDetTrackZ0AtPV*sin(InDetTrackParticles.theta)) < 5.0*mm )"
+    FTAG1Akt4PFlowJetTPThinningTool = acc.getPrimaryAndMerge(JetTrackParticleThinningCfg(
+        ConfigFlags,
+        name                         = "FTAG1Akt4PFlowJetTPThinningTool",
+        StreamName                   = kwargs['StreamName'],
+        JetKey                       = "AntiKt4EMPFlowJets",
+        SelectionString              = "",
+        #SelectionString              = "AntiKt4EMPFlowJets.pt > 18*GeV",
+        TrackSelectionString         = FTAG1_thinning_expression,
+        InDetTrackParticlesKey       = "InDetTrackParticles"))
+
+    thinningTools = [
+            FTAG1Akt4PFlowJetTPThinningTool,
+            ]
     # Finally the kernel itself
     DerivationKernel = CompFactory.DerivationFramework.DerivationKernel
     acc.addEventAlgo(DerivationKernel(name, ThinningTools = thinningTools))