Skip to content
Snippets Groups Projects
Commit 6d4b1706 authored by Carlos Eduardo Cocha Toapaxi's avatar Carlos Eduardo Cocha Toapaxi Committed by Linxuan Zhu
Browse files

New AP with 24c3 & 24c4 for dst_to_dee data (supersedes !1788)

parent dde8b9d8
No related branches found
No related tags found
No related merge requests found
Expand all Hide whitespace changes
Inline Side-by-side
dst_to_dee/DV.py 0 → 100644
+ 362
− 0
View file @ 6d4b1706
import Functors as F
from Functors.math import in_range
from GaudiKernel.SystemOfUnits import MeV
from Functors.math import log
from Hlt2Conf.algorithms_thor import ParticleFilter
from DaVinci import Options, make_config
from DaVinci.algorithms import create_lines_filter
from PyConf.reading import get_particles
from FunTuple import FunctorCollection
from PyConf.reading import get_particles, get_pvs
import FunTuple.functorcollections as FC
from FunTuple import FunTuple_Particles as Funtuple
from DecayTreeFitter import DecayTreeFitter
from Hlt2Conf.lines.charm.dst_to_dee_makers import (dst_BDT_functor)
from .tupling import (
make_DeltaM_variable,
make_basic_variables,
make_composite_variables,
make_composite_variables_3body,
make_composite_variables_4body,
make_hlt2_event_variables,
make_basic_dtf_variables,
make_composite_dtf_variables,
make_composite_dtf_variables_3body,
make_composite_dtf_variables_4body,
Hlt1_lines,
photon_isolation_variables
)
extra_brem_variables= (FunctorCollection(
{
"BREMHYPODELTAX": F.BREMHYPODELTAX,
"BREMHYPOENERGY": F.BREMHYPOENERGY,
"BREMHYPOID": F.BREMHYPOID,
"BREMHYPOMATCH_CHI2": F.BREMHYPOMATCH_CHI2,
"BREMPIDE": F.BREMPIDE,
"INBREM": F.INBREM,
"MASS_WITH_BREM": F.MASS_WITH_BREM,
"PT_WITH_BREM": F.PT_WITH_BREM,
"P_WITH_BREM": F.P_WITH_BREM,
"ECALPIDE": F.ECALPIDE,
"HCALPIDE": F.HCALPIDE,
"ELECTRONSHOWEREOP": F.ELECTRONSHOWEREOP,
"CLUSTERMATCH_CHI2": F.CLUSTERMATCH_CHI2,
"ELECTRONMATCH_CHI2": F.ELECTRONMATCH_CHI2,
"ELECTRONENERGY": F.ELECTRONENERGY,
"ELECTRONID": F.ELECTRONID,
"HCALEOP": F.HCALEOP,
"CALO_NEUTRAL_SHOWER_SHAPE": F.CALO_NEUTRAL_SHOWER_SHAPE,
"RICH_DLL_E": F.VALUE_OR(F.NaN)@F.RICH_DLL_E,
"TRACKPX": F.VALUE_OR(-1) @ F.TRACK_PX,
"TRACKPY": F.VALUE_OR(-1) @ F.TRACK_PY,
"TRACKPZ": F.VALUE_OR(-1) @ F.TRACK_PZ,
})
)
decay_descriptor = {
'dst_kpi_os' : {
"Dst0": "[ D*(2007)0 -> (D0 -> K- pi+) (gamma -> e+ e-)]CC",
"D0": "[ D*(2007)0 -> ^(D0 -> K- pi+) (gamma -> e+ e-)]CC",
"Kminus": "[ D*(2007)0 -> (D0 -> ^K- pi+) (gamma -> e+ e-)]CC",
"piplus": "[ D*(2007)0 -> (D0 -> K- ^pi+) (gamma -> e+ e-)]CC",
"gamma": "[ D*(2007)0 -> (D0 -> K- pi+) ^(gamma -> e+ e-)]CC",
"eplus": "[ D*(2007)0 -> (D0 -> K- pi+) (gamma -> ^e+ e-)]CC",
"eminus": "[ D*(2007)0 -> (D0 -> K- pi+) (gamma -> e+ ^e-)]CC",
},
'dst_kpi_ss' : {
"Dst0": "[ D*(2007)0 -> (D0 -> K- pi+) ([gamma -> e+ e+]CC)]CC",
"D0": "[ D*(2007)0 -> ^(D0 -> K- pi+) ([gamma -> e+ e+]CC)]CC",
"Kminus": "[ D*(2007)0 -> (D0 -> ^K- pi+) ([gamma -> e+ e+]CC)]CC",
"piplus": "[ D*(2007)0 -> (D0 -> K- ^pi+) ([gamma -> e+ e+]CC)]CC",
"gamma": "[ D*(2007)0 -> (D0 -> K- pi+) ^([gamma -> e+ e+]CC)]CC",
"eplus": "[ D*(2007)0 -> (D0 -> K- pi+) ([gamma -> ^e+ e+]CC)]CC",
"eminus": "[ D*(2007)0 -> (D0 -> K- pi+) ([gamma -> e+ ^e+]CC)]CC",
},
'dst_k3pi_os' : {
"Dst0": "[ D*(2007)0 -> (D0 -> K- pi- pi+ pi+) (gamma -> e+ e-)]CC",
"D0": "[ D*(2007)0 -> ^(D0 -> K- pi- pi+ pi+) (gamma -> e+ e-)]CC",
"Kminus": "[ D*(2007)0 -> (D0 -> ^K- pi- pi+ pi+) (gamma -> e+ e-)]CC",
"piminus": "[ D*(2007)0 -> (D0 -> K- ^pi- pi+ pi+) (gamma -> e+ e-)]CC",
"piplus1": "[ D*(2007)0 -> (D0 -> K- pi- ^pi+ pi+) (gamma -> e+ e-)]CC",
"piplus2": "[ D*(2007)0 -> (D0 -> K- pi- pi+ ^pi+) (gamma -> e+ e-)]CC",
"gamma": "[ D*(2007)0 -> (D0 -> K- pi- pi+ pi+) ^(gamma -> e+ e-)]CC",
"eplus": "[ D*(2007)0 -> (D0 -> K- pi- pi+ pi+) (gamma -> ^e+ e-)]CC",
"eminus": "[ D*(2007)0 -> (D0 -> K- pi- pi+ pi+) (gamma -> e+ ^e-)]CC",
},
'dst_k3pi_ss' : {
"Dst0": "[ D*(2007)0 -> (D0 -> K- pi- pi+ pi+) ([gamma -> e+ e+]CC)]CC",
"D0": "[ D*(2007)0 -> ^(D0 -> K- pi- pi+ pi+) ([gamma -> e+ e+]CC)]CC",
"Kminus": "[ D*(2007)0 -> (D0 -> ^K- pi- pi+ pi+) ([gamma -> e+ e+]CC)]CC",
"piminus": "[ D*(2007)0 -> (D0 -> K- ^pi- pi+ pi+) ([gamma -> e+ e+]CC)]CC",
"piplus1": "[ D*(2007)0 -> (D0 -> K- pi- ^pi+ pi+) ([gamma -> e+ e+]CC)]CC",
"piplus2": "[ D*(2007)0 -> (D0 -> K- pi- pi+ ^pi+) ([gamma -> e+ e+]CC)]CC",
"gamma": "[ D*(2007)0 -> (D0 -> K- pi- pi+ pi+) ^([gamma -> e+ e+]CC)]CC",
"eplus": "[ D*(2007)0 -> (D0 -> K- pi- pi+ pi+) ([gamma -> ^e+ e+]CC)]CC",
"eminus": "[ D*(2007)0 -> (D0 -> K- pi- pi+ pi+) ([gamma -> e+ ^e+]CC)]CC",
},
'dsstp_2kpi_os' : {
"Dstp": "[ D*_s+ -> (D_s+ -> K- K+ pi+) (gamma -> e+ e-)]CC",
"DpDs": "[ D*_s+ -> ^(D_s+ -> K- K+ pi+) (gamma -> e+ e-)]CC",
"Kminus": "[ D*_s+ -> (D_s+ -> ^K- K+ pi+) (gamma -> e+ e-)]CC",
"Kplus": "[ D*_s+ -> (D_s+ -> K- ^K+ pi+) (gamma -> e+ e-)]CC",
"piplus": "[ D*_s+ -> (D_s+ -> K- K+ ^pi+) (gamma -> e+ e-)]CC",
"gamma": "[ D*_s+ -> (D_s+ -> K- K+ pi+) ^(gamma -> e+ e-)]CC",
"eplus": "[ D*_s+ -> (D_s+ -> K- K+ pi+) (gamma -> ^e+ e-)]CC",
"eminus": "[ D*_s+ -> (D_s+ -> K- K+ pi+) (gamma -> e+ ^e-)]CC",
},
'dsstp_2kpi_ss' : {
"Dstp": "[ D*_s+ -> (D_s+ -> K- K+ pi+) ([gamma -> e+ e+]CC)]CC",
"DpDs": "[ D*_s+ -> ^(D_s+ -> K- K+ pi+) ([gamma -> e+ e+]CC)]CC",
"Kminus": "[ D*_s+ -> (D_s+ -> ^K- K+ pi+) ([gamma -> e+ e+]CC)]CC",
"Kplus": "[ D*_s+ -> (D_s+ -> K- ^K+ pi+) ([gamma -> e+ e+]CC)]CC",
"piplus": "[ D*_s+ -> (D_s+ -> K- K+ ^pi+) ([gamma -> e+ e+]CC)]CC",
"gamma": "[ D*_s+ -> (D_s+ -> K- K+ pi+) ^([gamma -> e+ e+]CC)]CC",
"eplus": "[ D*_s+ -> (D_s+ -> K- K+ pi+) ([gamma -> ^e+ e+]CC)]CC",
"eminus": "[ D*_s+ -> (D_s+ -> K- K+ pi+) ([gamma -> e+ ^e+]CC)]CC",
},
}
def MVA_dst_variables(line, pvs):
return (FunctorCollection( { "MVA": dst_BDT_functor(pvs, line=line)}))
#make_dtf_variables(options, input_data, ptype):
def make_dtf_variables(options, input_data, pvs, ptype="basic", mass_constrain = ["D*(2007)0", "D0" ], nparticles=2):
DTF = DecayTreeFitter(
name="DTF_{hash}",
input_particles=input_data,
)
DTF_DMass_BestPV = DecayTreeFitter(
name="DTF_DMass_BestPV_{hash}",
input_particles=input_data,
mass_constraints=mass_constrain,
input_pvs=pvs,
fit_all_pvs=False,
)
DTF_Mass_BestPV = DecayTreeFitter(
name="DTF_Mass_BestPV_{hash}",
input_particles=input_data,
mass_constraints=[mass_constrain[1]],
input_pvs=pvs,
fit_all_pvs=False,
)
DTF_BestPV = DecayTreeFitter(
name="DTF_BestPV_{hash}",
input_particles=input_data,
#mass_constraints=["D*(2007)0", "D0" ],
input_pvs=pvs,
fit_all_pvs=False,
)
DTF_DMass = DecayTreeFitter(
name="DTF_DMass_{hash}",
input_particles=input_data,
mass_constraints=mass_constrain,
#output_level=3,
)
DTF_Mass = DecayTreeFitter(
name="DTF_Mass_{hash}",
input_particles=input_data,
mass_constraints=[mass_constrain[1]],
#output_level=3,
)
if ptype == "basic":
dtf_variables = make_basic_dtf_variables
elif ptype == "composite":
if nparticles == 2:
dtf_variables = make_composite_dtf_variables
elif nparticles == 3:
dtf_variables = make_composite_dtf_variables_3body
elif nparticles == 4:
dtf_variables = make_composite_dtf_variables_4body
#dtf_vars = dtf_variables(options, pvs, input_data,
# DTF=DTF,
# pv_constraint=False,
# mass_constraint=False)
dtf_vars = dtf_variables(options, pvs, input_data,
DTF=DTF_BestPV,
pv_constraint=True,
mass_constraint=False)
#dtf_vars += dtf_variables(options, pvs, input_data,
# DTF=DTF_Mass,
# pv_constraint=False,
# mass_constraint=True, particle_name="D")
#dtf_vars += dtf_variables(options, pvs, input_data,
# DTF=DTF_Mass_BestPV,
# pv_constraint=True,
# mass_constraint=True, particle_name="D")
#dtf_vars += dtf_variables(options, pvs, input_data,
# DTF=DTF_DMass,
# pv_constraint=False,
# mass_constraint=True, particle_name="DstD")
dtf_vars += dtf_variables(options, pvs, input_data,
DTF=DTF_DMass_BestPV,
pv_constraint=True,
mass_constraint=True, particle_name="DstD")
return dtf_vars
CHILD_1 = lambda func: F.CHILD(1, func)
CHILD_2 = lambda func: F.CHILD(2, func)
def make_Dst0ToD0EmEp_D0ToKmPip_tuple(options, line, pvs , rec_summary, dd='dst_kpi_os'):
init_data=get_particles(f"/Event/HLT2/{line}/Particles")
extra_photons = get_particles(f"/Event/HLT2/{line}/{line[10:]}_prompt_Photons/Particles")
my_filter = create_lines_filter(f"LineFilter_{line}_{{hash}}",[line],)
input_data = ParticleFilter(init_data, F.FILTER(F.require_all(
in_range(1900. * MeV, F.MASS, 2100. * MeV),
in_range(1825. * MeV, CHILD_1(F.MASS), 1900. * MeV),
F.MAXDOCA < 0.2,
CHILD_1(F.MAXDOCA) < 0.2,
CHILD_2(F.MAXDOCA) < 0.2,)))
fields = decay_descriptor[dd]
extra_variables= photon_isolation_variables(input_data, extra_photons)
composite_tuple_variables = {
"Dst0" : make_composite_variables(options, pvs, input_data) + make_dtf_variables(options, input_data, pvs, "composite") +\
MVA_dst_variables("Hlt2Charm_Dst0ToD0EmEp_D0ToKmPip", pvs) + make_DeltaM_variable(options),
"D0" : make_composite_variables(options, pvs, input_data) + make_dtf_variables(options, input_data, pvs, "composite"),
"gamma" : make_composite_variables(options, pvs, input_data) + make_dtf_variables(options, input_data, pvs, "composite")+extra_variables
}
basic_tuple_variables = {
"Kminus" : make_basic_variables(options, pvs, input_data)+make_dtf_variables(options, input_data, pvs, "basic"),
"piplus" : make_basic_variables(options, pvs, input_data)+make_dtf_variables(options, input_data, pvs, "basic"),
"eplus" : make_basic_variables(options, pvs, input_data)+make_dtf_variables(options, input_data, pvs, "basic") + extra_brem_variables,
"eminus" : make_basic_variables(options, pvs, input_data)+make_dtf_variables(options, input_data, pvs, "basic") + extra_brem_variables,
}
my_tuple = Funtuple(
name=f"Tuple_{line}_{dd}",
tuple_name="DecayTree",
fields=fields,
variables= {**composite_tuple_variables, **basic_tuple_variables},
event_variables=make_hlt2_event_variables(options, pvs, rec_summary),
inputs=input_data,
store_multiple_cand_info=True
)
return [ my_filter, my_tuple]
def make_Dst0ToD0EmEp_D0ToKmPimPipPip_tuple(options, line, pvs, rec_summary, dd='dst_k3pi_os'):
init_data=get_particles(f"/Event/HLT2/{line}/Particles")
extra_photons = get_particles(f"/Event/HLT2/{line}/{line[10:]}_prompt_Photons/Particles")
my_filter = create_lines_filter(f"LineFilter_{line}_{{hash}}",[line],)
input_data = ParticleFilter(init_data, F.FILTER(F.require_all(
in_range(1900. * MeV, F.MASS, 2100. * MeV),
in_range(1825. * MeV, CHILD_1(F.MASS), 1900. * MeV),
F.MAXDOCA < 0.2,
CHILD_1(F.MAXDOCA) < 0.2,
CHILD_2(F.MAXDOCA) < 0.2,)))
fields = decay_descriptor[dd]
extra_variables= photon_isolation_variables(input_data, extra_photons)
composite_tuple_variables = {
"Dst0" : make_composite_variables(options, pvs, input_data) + make_dtf_variables(options, input_data, pvs, "composite") +\
MVA_dst_variables("Hlt2Charm_Dst0ToD0EmEp_D0ToKmPimPipPip", pvs) + make_DeltaM_variable(options),
"D0" : make_composite_variables_4body(options, pvs, input_data) + make_dtf_variables(options, input_data, pvs, "composite", nparticles=4),
"gamma" : make_composite_variables(options, pvs, input_data) + make_dtf_variables(options, input_data, pvs, "composite")+extra_variables
}
basic_tuple_variables = {
"Kminus" : make_basic_variables(options, pvs, input_data)+make_dtf_variables(options, input_data, pvs, "basic"),
"piminus" : make_basic_variables(options, pvs, input_data)+make_dtf_variables(options, input_data, pvs, "basic"),
"piplus1" : make_basic_variables(options, pvs, input_data)+make_dtf_variables(options, input_data, pvs, "basic"),
"piplus2" : make_basic_variables(options, pvs, input_data)+make_dtf_variables(options, input_data, pvs, "basic"),
"eplus" : make_basic_variables(options, pvs, input_data)+make_dtf_variables(options, input_data, pvs, "basic") + extra_brem_variables,
"eminus" : make_basic_variables(options, pvs, input_data)+make_dtf_variables(options, input_data, pvs, "basic") + extra_brem_variables,
}
my_tuple = Funtuple(
name=f"Tuple_{line}_{dd}",
tuple_name="DecayTree",
fields=fields,
variables= {**composite_tuple_variables, **basic_tuple_variables},
event_variables=make_hlt2_event_variables(options, pvs, rec_summary),
inputs=input_data,
store_multiple_cand_info=True
)
return [ my_filter, my_tuple]
def make_DstpToDpDspEmEp_DpDspToKmKpPip_tuple(options, line, pvs , rec_summary, dd='dsstp_2kpi_os'):
init_data=get_particles(f"/Event/HLT2/{line}/Particles")
extra_photons = get_particles(f"/Event/HLT2/{line}/{line[10:]}_prompt_Photons/Particles")
my_filter = create_lines_filter(f"LineFilter_{line}_{{hash}}",[line],)
input_data = ParticleFilter(init_data, F.FILTER(F.require_all(
in_range(1900. * MeV, F.MASS, 2300. * MeV),
in_range(1900. * MeV, CHILD_1(F.MASS), 2049. * MeV),
F.MAXDOCA < 0.2,
CHILD_1(F.MAXDOCA) < 0.2,
CHILD_2(F.MAXDOCA) < 0.2,)))
fields = decay_descriptor[dd]
extra_variables= photon_isolation_variables(input_data, extra_photons)
composite_tuple_variables = {
"Dstp" : make_composite_variables(options, pvs, input_data) + make_dtf_variables(options, input_data, pvs, "composite", mass_constrain = ["D*_s+", "D_s+"]) +\
MVA_dst_variables("Hlt2Charm_DstpToDpDspEmEp_DpDspToKmKpPip", pvs) + make_DeltaM_variable(options),
"DpDs" : make_composite_variables_3body(options, pvs, input_data) + make_dtf_variables(options, input_data, pvs, "composite", mass_constrain = ["D*_s+", "D_s+"], nparticles=3),
"gamma" : make_composite_variables(options, pvs, input_data) + make_dtf_variables(options, input_data, pvs, "composite", mass_constrain = ["D*_s+", "D_s+"])+extra_variables,
}
basic_tuple_variables = {
"Kminus" : make_basic_variables(options, pvs, input_data)+make_dtf_variables(options, input_data, pvs, "basic", mass_constrain = ["D*_s+", "D_s+"]),
"Kplus" : make_basic_variables(options, pvs, input_data)+make_dtf_variables(options, input_data, pvs, "basic", mass_constrain = ["D*_s+", "D_s+"]),
"piplus" : make_basic_variables(options, pvs, input_data)+make_dtf_variables(options, input_data, pvs, "basic", mass_constrain = ["D*_s+", "D_s+"]),
"eplus" : make_basic_variables(options, pvs, input_data)+make_dtf_variables(options, input_data, pvs, "basic", mass_constrain = ["D*_s+", "D_s+"])+ extra_brem_variables,
"eminus" : make_basic_variables(options, pvs, input_data)+make_dtf_variables(options, input_data, pvs, "basic", mass_constrain = ["D*_s+", "D_s+"])+ extra_brem_variables,
}
my_tuple = Funtuple(
name=f"Tuple_{line}_{dd}",
tuple_name="DecayTree",
fields=fields,
variables= {**composite_tuple_variables, **basic_tuple_variables},
event_variables=make_hlt2_event_variables(options, pvs, rec_summary),
inputs=input_data,
store_multiple_cand_info=True
)
return [ my_filter, my_tuple]
from DaVinci import Options, make_config
from PyConf.reading import get_pvs, get_rec_summary
def main(options: Options):
# get ODIN and DecReports location
pvs = get_pvs()
rec_summary = get_rec_summary()
tuples = {
'Hlt2Charm_Dst0ToD0EmEp_D0ToKmPip_OS': make_Dst0ToD0EmEp_D0ToKmPip_tuple(options, "Hlt2Charm_Dst0ToD0EmEp_D0ToKmPip_OS", pvs, rec_summary, 'dst_kpi_os'),
'Hlt2Charm_Dst0ToD0EmEp_D0ToKmPip_SS': make_Dst0ToD0EmEp_D0ToKmPip_tuple(options, "Hlt2Charm_Dst0ToD0EmEp_D0ToKmPip_SS", pvs, rec_summary, 'dst_kpi_ss'),
'Hlt2Charm_Dst0ToD0EmEp_D0ToKmPip_OS_MVA': make_Dst0ToD0EmEp_D0ToKmPip_tuple(options, "Hlt2Charm_Dst0ToD0EmEp_D0ToKmPip_OS_MVA", pvs, rec_summary, 'dst_kpi_os'),
'Hlt2Charm_Dst0ToD0EmEp_D0ToKmPip_SS_MVA': make_Dst0ToD0EmEp_D0ToKmPip_tuple(options, "Hlt2Charm_Dst0ToD0EmEp_D0ToKmPip_SS_MVA", pvs, rec_summary, 'dst_kpi_ss'),
'Hlt2Charm_Dst0ToD0EmEp_D0ToKmPimPipPip_OS': make_Dst0ToD0EmEp_D0ToKmPimPipPip_tuple(options, "Hlt2Charm_Dst0ToD0EmEp_D0ToKmPimPipPip_OS", pvs, rec_summary,'dst_k3pi_os' ),
'Hlt2Charm_Dst0ToD0EmEp_D0ToKmPimPipPip_SS': make_Dst0ToD0EmEp_D0ToKmPimPipPip_tuple(options, "Hlt2Charm_Dst0ToD0EmEp_D0ToKmPimPipPip_SS", pvs, rec_summary, 'dst_k3pi_ss'),
'Hlt2Charm_Dst0ToD0EmEp_D0ToKmPimPipPip_OS_MVA': make_Dst0ToD0EmEp_D0ToKmPimPipPip_tuple(options, "Hlt2Charm_Dst0ToD0EmEp_D0ToKmPimPipPip_OS_MVA", pvs, rec_summary, 'dst_k3pi_os'),
'Hlt2Charm_Dst0ToD0EmEp_D0ToKmPimPipPip_SS_MVA': make_Dst0ToD0EmEp_D0ToKmPimPipPip_tuple(options, "Hlt2Charm_Dst0ToD0EmEp_D0ToKmPimPipPip_SS_MVA", pvs, rec_summary, 'dst_k3pi_ss'),
'Hlt2Charm_DstpToDpDspEmEp_DpDspToKmKpPip_OS': make_DstpToDpDspEmEp_DpDspToKmKpPip_tuple(options, "Hlt2Charm_DstpToDpDspEmEp_DpDspToKmKpPip_OS", pvs, rec_summary, "dsstp_2kpi_os"),
'Hlt2Charm_DstpToDpDspEmEp_DpDspToKmKpPip_SS': make_DstpToDpDspEmEp_DpDspToKmKpPip_tuple(options, "Hlt2Charm_DstpToDpDspEmEp_DpDspToKmKpPip_SS", pvs, rec_summary, "dsstp_2kpi_ss"),
'Hlt2Charm_DstpToDpDspEmEp_DpDspToKmKpPip_OS_MVA': make_DstpToDpDspEmEp_DpDspToKmKpPip_tuple(options, "Hlt2Charm_DstpToDpDspEmEp_DpDspToKmKpPip_OS_MVA", pvs, rec_summary, "dsstp_2kpi_os"),
'Hlt2Charm_DstpToDpDspEmEp_DpDspToKmKpPip_SS_MVA': make_DstpToDpDspEmEp_DpDspToKmKpPip_tuple(options, "Hlt2Charm_DstpToDpDspEmEp_DpDspToKmKpPip_SS_MVA", pvs, rec_summary, "dsstp_2kpi_ss"),
}
config = make_config(options, tuples)
return config
dst_to_dee/README.md 0 → 100644
+ 5
− 0
View file @ 6d4b1706
# Analysis production for dst_to_dee charm module.
Creates tuples for Charm 2024-validation MC and data
## Analysis motivation
The module dst_to_dee contains several HLT2 lines to search for dark photons produced by a Charm decay.
\ No newline at end of file
dst_to_dee/info.yaml 0 → 100644
+ 42
− 0
View file @ 6d4b1706
defaults:
inform:
- carlos.eduardo.cocha.toapaxi@cern.ch
wg: Charm
{%- set datainfo = [
( 'MagUp', '24c2'),
( 'MagDown', '24c2'),
( 'MagUp', '24c3'),
( 'MagDown', '24c3'),
( 'MagUp', '24c4'),
( 'MagDown', '24c4'),
]%}
{%- for polarity, sprucen in datainfo %}
DATA_2024_{{ polarity }}_{{sprucen}}:
application: "DaVinci/v64r10@x86_64_v2-el9-clang16-opt"
input:
bk_query: "/LHCb/Collision24/Beam6800GeV-VeloClosed-{{polarity}}/Real Data/Sprucing{{sprucen}}/94000000/CHARM.DST"
dq_flags:
- UNCHECKED
- OK
keep_running: True
output: DVTUPLE.ROOT
options:
entrypoint: dst_to_dee.DV:main
extra_options:
input_raw_format: 0.5
input_type: ROOT
simulation: False
data_type: "Upgrade"
event_store: HiveWhiteBoard
geometry_version: 'run3/2024.Q1.2-v00.00'
conditions_version: master
input_process: "TurboPass"
input_stream: "charm"
compression:
algorithm: ZSTD
level: 1
max_buffer_size: 1048576
{%- endfor %}
dst_to_dee/tupling.py 0 → 100644
+ 1023
− 0
View file @ 6d4b1706
This diff is collapsed.
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment