コード例 #1
0
ファイル: JetToolSupport.py プロジェクト: dguest/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:
     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
コード例 #2
0
 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
コード例 #3
0
        # 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]