def addJetTrimmer(self, output, rclus, ptfrac, input, modifiersin ="groomed",
pseudojetRetriever ="jpjretriever",
isTrigger =False, useTriggerStore =False, doArea =True):
from JetRec.JetRecConf import JetTrimmer
from JetRec.JetRecConf import JetRecTool
groomer = JetTrimmer(output + "Groomer")
groomer.RClus = rclus
groomer.PtFrac = ptfrac
if doArea:
groomer.JetBuilder = self.jetBuilderWithArea
else:
groomer.JetBuilder = self.jetBuilderWithoutArea
self += groomer
jetrec = JetRecTool(output)
jetrec.JetGroomer = groomer
jetrec.InputContainer = input
jetrec.OutputContainer = output
jetrec.JetModifiers = self.getModifiers(modifiersin)
jetrec.Trigger = isTrigger or useTriggerStore
jetrec.Timer = jetFlags.timeJetRecTool()
if pseudojetRetriever in self.tools:
jetrec.JetPseudojetRetriever = self.tools[pseudojetRetriever]
else:
print "Requested jet pseudojet retriever is not a registered tool: " \
+ pseudojetRetriever
raise KeyError
self += jetrec
if isTrigger:
self.trigjetrecs += [jetrec]
else:
self.jetrecs += [jetrec]
self.jetcons += [output]
return jetrec
def addJetTrimmer(self, output, rclus, ptfrac, input, modifiersin ="groomed",
pseudojetRetriever ="jpjretriever",
isTrigger =False, useTriggerStore =False, doArea =True):
from JetRec.JetRecConf import JetTrimmer
from JetRec.JetRecConf import JetRecTool
groomer = JetTrimmer(output + "Groomer")
groomer.RClus = rclus
groomer.PtFrac = ptfrac
if doArea:
groomer.JetBuilder = self.jetBuilderWithArea
else:
groomer.JetBuilder = self.jetBuilderWithoutArea
self += groomer
jetrec = JetRecTool(output)
jetrec.JetGroomer = groomer
jetrec.InputContainer = input
jetrec.OutputContainer = output
jetrec.JetModifiers = self.getModifiers(modifiersin)
jetrec.Trigger = isTrigger or useTriggerStore
jetrec.Timer = jetFlags.timeJetRecTool()
if pseudojetRetriever in self.tools:
jetrec.JetPseudojetRetriever = self.tools[pseudojetRetriever]
else:
jetlog.info( "Requested jet pseudojet retriever is not a registered tool: " \
+ pseudojetRetriever )
raise KeyError
self += jetrec
if isTrigger:
self.trigjetrecs += [jetrec]
else:
self.jetrecs += [jetrec]
self.jetcons += [output]
return jetrec
# jetrec.JetModifiers += [pull]
# jetrec.JetModifiers += [charge]
# jetrec.JetModifiers += [subjetmaker]
# jetrec.JetModifiers += [subjetfinder]
ToolSvc += jetrec
jetPseudojetRetriever = JetPseudojetRetriever(
"JetPseudojetRetriever_" + name)
ToolSvc += jetPseudojetRetriever
# JetRec tool for finding.
jetrec_trimm = JetRecTool("JetRecTool_" + name + "Trimm")
jetrec_trimm.JetGroomer = groomer
jetrec_trimm.InputContainer = "AntiKt10" + name + "Jets"
jetrec_trimm.OutputContainer = "AntiKt10" + name + "TrimmedJets"
jetrec_trimm.JetPseudojetRetriever = jetPseudojetRetriever
#jetrec_trimm.JetModifiers += [nsubjettiness]
#jetrec_trimm.JetModifiers += [nsubjettinessratios]
# jetrec_trimm.JetModifiers += [ktsplittingscale]
# jetrec_trimm.JetModifiers += [dipolarity]
# jetrec_trimm.JetModifiers += [angularity]
# jetrec_trimm.JetModifiers += [ktdr]
# jetrec_trimm.JetModifiers += [ktmassdrop]
# jetrec_trimm.JetModifiers += [planarflow]
# jetrec_trimm.JetModifiers += [centerofmassshapes]
#jetrec_trimm.JetModifiers += [energycorrelator]
#jetrec_trimm.JetModifiers += [energycorrelatorratios]
# jetrec_trimm.JetModifiers += [pull]
# jetrec_trimm.JetModifiers += [charge]
# jetrec_trimm.JetModifiers += [subjetmaker]
# jetrec_trimm.JetModifiers += [subjetfinder]