Skip to content
Snippets Groups Projects

Associated Jpsi etac production (MC for Jpsi to pp)

Merged Raoul Henderson requested to merge associated_jpsi_etac_v3 into master
1 unresolved thread
Files
7
+ 68
0
 
# Processes the event with no selection, for acceptance efficiency
 
 
def fillMCTuple(MCTuple):
 
 
from Configurables import TupleToolMCTruth
 
MCTruth = TupleToolMCTruth("MCTruth")
 
MCTuple.addTool(MCTruth)
 
MCTruth.ToolList = ["MCTupleToolAngles"
 
, "MCTupleToolHierarchy"
 
, "MCTupleToolKinematic"
 
, "MCTupleToolReconstructed"
 
, "MCTupleToolPID"
 
, "MCTupleToolDecayType"
 
, "MCTupleToolPrompt"
 
]
 
 
from LoKiMC.functions import MCMOTHER, MCVFASPF, MCPRIMARY
 
from Configurables import LoKi__Hybrid__MCTupleTool
 
MCTuple.addTool( LoKi__Hybrid__MCTupleTool, name = "MCLoKiHybrid" )
 
MCTuple.ToolList += [ "LoKi::Hybrid::MCTupleTool/MCLoKiHybrid" ]
 
MCTuple.MCLoKiHybrid.Preambulo = [ "from LoKiCore.functions import switch" ]
 
MCTuple.MCLoKiHybrid.Variables = {
 
"TRUEID" : "MCID",
 
"TRUEY" : "MCY",
 
"TRUEETA" : "MCETA",
 
"TRUEM" : "MCM",
 
"MCVZ" : "MCVFASPF(MCVZ)",
 
"TRUEETA" : "MCETA",
 
"TRUEPHI" : "MCPHI",
 
"TRUETHETA" : "MCTHETA",
 
"TRUE_Tz" : "(3.3)*(MCMOTHER(MCVFASPF(MCVZ), 0)-MCVFASPF(MCVZ))*MCM/MCPZ",
 
"MC_ISPRIMARY" : "MCVFASPF(switch(MCPRIMARY,1,0))"
 
}
 
 
from Configurables import LoKi__Hybrid__EvtTupleTool
 
MCTuple.addTool( LoKi__Hybrid__EvtTupleTool, name = "ELoKiHybrid" )
 
MCTuple.ToolList += [ "LoKi::Hybrid::EvtTupleTool/ELoKiHybrid" ]
 
MCTuple.ELoKiHybrid.Preambulo = [ "from LoKiCore.basic import LHCb" ]
 
MCTuple.ELoKiHybrid.VOID_Variables = {
 
"nSPDHits" : "RECSUMMARY( LHCb.RecSummary.nSPDhits, -1, '/Event/Rec/Summary', False )"
 
}
 
 
return MCTuple
 
 
def build_algo_MC(parent, parent_name, child1, child2, child_name):
 
from Gaudi.Configuration import *
 
from LHCbKernel.Configuration import *
 
from DaVinci.Configuration import *
 
from DecayTreeTuple.Configuration import *
 
from Configurables import DaVinci
 
 
Decay = "{} -> ^{} ^{}".format(parent, child1, child2)
 
Branches = {
 
parent_name: "{} -> {} {}".format(parent, child1, child2),
 
"{}P".format(child_name): "{} -> ^{} {}".format(parent, child1, child2),
 
"{}M".format(child_name): "{} -> {} ^{}".format(parent, child1, child2)
 
}
 
 
from Configurables import MCDecayTreeTuple
 
MCTuple = MCDecayTreeTuple("MCDecayTreeTuple")
 
MCTuple.Decay = Decay
 
MCTuple.Branches = Branches
 
 
fillMCTuple(MCTuple)
 
 
DaVinci().UserAlgorithms += [MCTuple]
 
from Configurables import MessageSvc
 
MessageSvc().setWarning = [ 'RFileCnv' ]
 
\ No newline at end of file
Loading