Add DecFiles for DM searches

528c8707 · Saul Lopez Solino · 4a9013a7 · 528c8707 · 528c8707 · 528c8707
Commit 528c8707 authored Sep 30, 2021 by Saul Lopez Solino
--- a/dkfiles/B0_PsiDMLambda1520,Lambda1520_pK=FullGenEvtCut,mPsiDM=1500MeV.dec
+++ b/dkfiles/B0_PsiDMLambda1520,Lambda1520_pK=FullGenEvtCut,mPsiDM=1500MeV.dec
+# EventType: 11502001
+#
+# Descriptor: [B0 -> (Lambda(1520)0 -> p+ K-) H_30 ]cc
+#
+# NickName: B0_PsiDMLambda1520,Lambda1520_pK=FullGenEvtCut,mPsiDM=1500MeV
+#
+# Cuts: None
+# FullEventCuts: LoKi::FullGenEventCut/B0toDM
+#
+# Documentation:
+#   Decay a B0 to a Lambda(1520)0 -> pK and a redefined H_30 for our need, acting the latter as stable Dark Matter candidate. The mass of the H_30 is 1.5 GeV.
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# CPUTime: < 1 min
+# Responsible: Saul Lopez 
+# Email: saul.lopez.solino@cern.ch
+# Date: 20210930
+#
+#
+# InsertPythonCode:
+# from Configurables import LHCb__ParticlePropertySvc, LoKi__GenCutTool
+# LHCb__ParticlePropertySvc().Particles = [
+# "H_30     89       36      0.0     1.500000        1.000000e+16    A0      36      0.00"
+# ]
+# ## Generator level cuts:
+# from Configurables import LoKi__FullGenEventCut
+# Generation().addTool( LoKi__FullGenEventCut, "B0toDM" )
+# tracksInAcc = Generation().B0toDM
+# tracksInAcc.Code = " count ( isGoodB ) > 0 "
+# ### - HepMC::IteratorRange::descendants   4
+# tracksInAcc.Preambulo += [ "from GaudiKernel.SystemOfUnits import GeV, mrad"
+#                          , "inAcc = in_range(1.9, GETA, 5.0)" 
+#                          , "isGoodKaon     = ( ( GPT > 0.50*GeV ) & inAcc & ( 'K+' == GABSID ) )"
+#                          , "isGoodP        = ( ( GPT > 0.50*GeV ) & inAcc & ( 'p+' == GABSID ) )"
+#                          , "isGoodL0       = ( ( 'Lambda(1520)0' == GABSID ) & ( GNINTREE( isGoodKaon, 1 ) > 0 ) & ( GNINTREE( isGoodP, 1 ) > 0 ) )"
+#                          , "isGoodB        = ( ( 'B0' == GABSID ) & ( GNINTREE( isGoodL0, 1 ) > 0 ) )" ]
+# EndInsertPythonCode
+#
+Alias      MyLambda(1520)0       Lambda(1520)0
+Alias      Myanti-Lambda(1520)0  anti-Lambda(1520)0
+ChargeConj MyLambda(1520)0       Myanti-Lambda(1520)0
+#
+Alias  MyH_30     A0
+Alias  Myanti-H_30    A0
+ChargeConj MyH_30   Myanti-H_30
+#
+Decay B0sig
+    1.000   MyLambda(1520)0   MyH_30    PHSP;
+Enddecay
+CDecay anti-B0sig
+#
+Decay MyLambda(1520)0
+    1.000       p+          K-              PHSP;
+Enddecay
+CDecay Myanti-Lambda(1520)0
+#
+End
--- a/dkfiles/B0_PsiDMLambda1520,Lambda1520_pK=FullGenEvtCut,mPsiDM=2000MeV.dec
+++ b/dkfiles/B0_PsiDMLambda1520,Lambda1520_pK=FullGenEvtCut,mPsiDM=2000MeV.dec
+# EventType: 11502002
+#
+# Descriptor: [B0 -> (Lambda(1520)0 -> p+ K-) H_30 ]cc
+#
+# NickName: B0_PsiDMLambda1520,Lambda1520_pK=FullGenEvtCut,mPsiDM=2000MeV
+#
+# Cuts: None
+# FullEventCuts: LoKi::FullGenEventCut/B0toDM
+#
+# Documentation:
+#   Decay a B0 to a Lambda(1520)0 -> pK and a redefined H_30 for our need, acting the latter as stable Dark Matter candidate. The mass of the H_30 is 2.00 GeV.
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# CPUTime: < 1 min
+# Responsible: Saul Lopez 
+# Email: saul.lopez.solino@cern.ch
+# Date: 20210930
+#
+#
+# InsertPythonCode:
+# from Configurables import LHCb__ParticlePropertySvc, LoKi__GenCutTool
+# LHCb__ParticlePropertySvc().Particles = [
+# "H_30     89       36      0.0     2.000000        1.000000e+16    A0      36      0.00"
+# ]
+# ## Generator level cuts:
+# from Configurables import LoKi__FullGenEventCut
+# Generation().addTool( LoKi__FullGenEventCut, "B0toDM" )
+# tracksInAcc = Generation().B0toDM
+# tracksInAcc.Code = " count ( isGoodB ) > 0 "
+# ### - HepMC::IteratorRange::descendants   4
+# tracksInAcc.Preambulo += [ "from GaudiKernel.SystemOfUnits import GeV, mrad"
+#                          , "inAcc = in_range(1.9, GETA, 5.0)" 
+#                          , "isGoodKaon     = ( ( GPT > 0.50*GeV ) & inAcc & ( 'K+' == GABSID ) )"
+#                          , "isGoodP        = ( ( GPT > 0.50*GeV ) & inAcc & ( 'p+' == GABSID ) )"
+#                          , "isGoodL0       = ( ( 'Lambda(1520)0' == GABSID ) & ( GNINTREE( isGoodKaon, 1 ) > 0 ) & ( GNINTREE( isGoodP, 1 ) > 0 ) )"
+#                          , "isGoodB        = ( ( 'B0' == GABSID ) & ( GNINTREE( isGoodL0, 1 ) > 0 ) )" ]
+# EndInsertPythonCode
+#
+Alias      MyLambda(1520)0       Lambda(1520)0
+Alias      Myanti-Lambda(1520)0  anti-Lambda(1520)0
+ChargeConj MyLambda(1520)0       Myanti-Lambda(1520)0
+#
+Alias  MyH_30     A0
+Alias  Myanti-H_30    A0
+ChargeConj MyH_30   Myanti-H_30
+#
+Decay B0sig
+    1.000   MyLambda(1520)0   MyH_30    PHSP;
+Enddecay
+CDecay anti-B0sig
+#
+Decay MyLambda(1520)0
+    1.000       p+          K-              PHSP;
+Enddecay
+CDecay Myanti-Lambda(1520)0
+#
+End
--- a/dkfiles/B0_PsiDMLambda1520,Lambda1520_pK=FullGenEvtCut,mPsiDM=2400MeV.dec
+++ b/dkfiles/B0_PsiDMLambda1520,Lambda1520_pK=FullGenEvtCut,mPsiDM=2400MeV.dec
+# EventType: 11502003
+#
+# Descriptor: [B0 -> (Lambda(1520)0 -> p+ K-) H_30 ]cc
+#
+# NickName: B0_PsiDMLambda1520,Lambda1520_pK=FullGenEvtCut,mPsiDM=2400MeV
+#
+# Cuts: None
+# FullEventCuts: LoKi::FullGenEventCut/B0toDM
+#
+# Documentation:
+#   Decay a B0 to a Lambda(1520)0 -> pK and a redefined H_30 for our need, acting the latter as stable Dark Matter candidate. The mass of the H_30 is 2.4 GeV.
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# CPUTime: < 1 min
+# Responsible: Saul Lopez 
+# Email: saul.lopez.solino@cern.ch
+# Date: 20210930
+#
+#
+# InsertPythonCode:
+# from Configurables import LHCb__ParticlePropertySvc, LoKi__GenCutTool
+# LHCb__ParticlePropertySvc().Particles = [
+# "H_30     89       36      0.0     2.400000        1.000000e+16    A0      36      0.00"
+# ]
+# ## Generator level cuts:
+# from Configurables import LoKi__FullGenEventCut
+# Generation().addTool( LoKi__FullGenEventCut, "B0toDM" )
+# tracksInAcc = Generation().B0toDM
+# tracksInAcc.Code = " count ( isGoodB ) > 0 "
+# ### - HepMC::IteratorRange::descendants   4
+# tracksInAcc.Preambulo += [ "from GaudiKernel.SystemOfUnits import GeV, mrad"
+#                          , "inAcc = in_range(1.9, GETA, 5.0)" 
+#                          , "isGoodKaon     = ( ( GPT > 0.50*GeV ) & inAcc & ( 'K+' == GABSID ) )"
+#                          , "isGoodP        = ( ( GPT > 0.50*GeV ) & inAcc & ( 'p+' == GABSID ) )"
+#                          , "isGoodL0       = ( ( 'Lambda(1520)0' == GABSID ) & ( GNINTREE( isGoodKaon, 1 ) > 0 ) & ( GNINTREE( isGoodP, 1 ) > 0 ) )"
+#                          , "isGoodB        = ( ( 'B0' == GABSID ) & ( GNINTREE( isGoodL0, 1 ) > 0 ) )" ]
+# EndInsertPythonCode
+#
+Alias      MyLambda(1520)0       Lambda(1520)0
+Alias      Myanti-Lambda(1520)0  anti-Lambda(1520)0
+ChargeConj MyLambda(1520)0       Myanti-Lambda(1520)0
+#
+Alias  MyH_30     A0
+Alias  Myanti-H_30    A0
+ChargeConj MyH_30   Myanti-H_30
+#
+Decay B0sig
+    1.000   MyLambda(1520)0   MyH_30    PHSP;
+Enddecay
+CDecay anti-B0sig
+#
+Decay MyLambda(1520)0
+    1.000       p+          K-              PHSP;
+Enddecay
+CDecay Myanti-Lambda(1520)0
+#
+End
--- a/dkfiles/B0_PsiDMLambda1520,Lambda1520_pK=FullGenEvtCut,mPsiDM=3000MeV.dec
+++ b/dkfiles/B0_PsiDMLambda1520,Lambda1520_pK=FullGenEvtCut,mPsiDM=3000MeV.dec
+# EventType: 11502004
+#
+# Descriptor: [B0 -> (Lambda(1520)0 -> p+ K-) H_30 ]cc
+#
+# NickName: B0_PsiDMLambda1520,Lambda1520_pK=FullGenEvtCut,mPsiDM=3000MeV
+#
+# Cuts: None
+# FullEventCuts: LoKi::FullGenEventCut/B0toDM
+#
+# Documentation:
+#   Decay a B0 to a Lambda(1520)0 -> pK and a redefined H_30 for our need, acting the latter as stable Dark Matter candidate. The mass of the H_30 is 3.0 GeV.
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# CPUTime: < 1 min
+# Responsible: Saul Lopez 
+# Email: saul.lopez.solino@cern.ch
+# Date: 20210930
+#
+#
+# InsertPythonCode:
+# from Configurables import LHCb__ParticlePropertySvc, LoKi__GenCutTool
+# LHCb__ParticlePropertySvc().Particles = [
+# "H_30     89       36      0.0     3.00000        1.000000e+16    A0      36      0.00"
+# ]
+# ## Generator level cuts:
+# from Configurables import LoKi__FullGenEventCut
+# Generation().addTool( LoKi__FullGenEventCut, "B0toDM" )
+# tracksInAcc = Generation().B0toDM
+# tracksInAcc.Code = " count ( isGoodB ) > 0 "
+# ### - HepMC::IteratorRange::descendants   4
+# tracksInAcc.Preambulo += [ "from GaudiKernel.SystemOfUnits import GeV, mrad"
+#                          , "inAcc = in_range(1.9, GETA, 5.0)" 
+#                          , "isGoodKaon     = ( ( GPT > 0.50*GeV ) & inAcc & ( 'K+' == GABSID ) )"
+#                          , "isGoodP        = ( ( GPT > 0.50*GeV ) & inAcc & ( 'p+' == GABSID ) )"
+#                          , "isGoodL0       = ( ( 'Lambda(1520)0' == GABSID ) & ( GNINTREE( isGoodKaon, 1 ) > 0 ) & ( GNINTREE( isGoodP, 1 ) > 0 ) )"
+#                          , "isGoodB        = ( ( 'B0' == GABSID ) & ( GNINTREE( isGoodL0, 1 ) > 0 ) )" ]
+# EndInsertPythonCode
+#
+Alias      MyLambda(1520)0       Lambda(1520)0
+Alias      Myanti-Lambda(1520)0  anti-Lambda(1520)0
+ChargeConj MyLambda(1520)0       Myanti-Lambda(1520)0
+#
+Alias  MyH_30     A0
+Alias  Myanti-H_30    A0
+ChargeConj MyH_30   Myanti-H_30
+#
+Decay B0sig
+    1.000   MyLambda(1520)0   MyH_30    PHSP;
+Enddecay
+CDecay anti-B0sig
+#
+Decay MyLambda(1520)0
+    1.000       p+          K-              PHSP;
+Enddecay
+CDecay Myanti-Lambda(1520)0
+#
+End
--- a/dkfiles/B0_PsiDMLambda1520,Lambda1520_pK=FullGenEvtCut,mPsiDM=3500MeV.dec
+++ b/dkfiles/B0_PsiDMLambda1520,Lambda1520_pK=FullGenEvtCut,mPsiDM=3500MeV.dec
+# EventType: 11502005
+#
+# Descriptor: [B0 -> (Lambda(1520)0 -> p+ K-) H_30 ]cc
+#
+# NickName: B0_PsiDMLambda1520,Lambda1520_pK=FullGenEvtCut,mPsiDM=3500MeV
+#
+# Cuts: None
+# FullEventCuts: LoKi::FullGenEventCut/B0toDM
+#
+# Documentation:
+#   Decay a B0 to a Lambda(1520)0 -> pK and a redefined H_30 for our need, acting the latter as stable Dark Matter candidate. The mass of the H_30 is 3.5 GeV.
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# CPUTime: < 1 min
+# Responsible: Saul Lopez 
+# Email: saul.lopez.solino@cern.ch
+# Date: 20210930
+#
+#
+# InsertPythonCode:
+# from Configurables import LHCb__ParticlePropertySvc, LoKi__GenCutTool
+# LHCb__ParticlePropertySvc().Particles = [
+# "H_30     89       36      0.0     3.500000        1.000000e+16    A0      36      0.00"
+# ]
+# ## Generator level cuts:
+# from Configurables import LoKi__FullGenEventCut
+# Generation().addTool( LoKi__FullGenEventCut, "B0toDM" )
+# tracksInAcc = Generation().B0toDM
+# tracksInAcc.Code = " count ( isGoodB ) > 0 "
+# ### - HepMC::IteratorRange::descendants   4
+# tracksInAcc.Preambulo += [ "from GaudiKernel.SystemOfUnits import GeV, mrad"
+#                          , "inAcc = in_range(1.9, GETA, 5.0)" 
+#                          , "isGoodKaon     = ( ( GPT > 0.50*GeV ) & inAcc & ( 'K+' == GABSID ) )"
+#                          , "isGoodP        = ( ( GPT > 0.50*GeV ) & inAcc & ( 'p+' == GABSID ) )"
+#                          , "isGoodL0       = ( ( 'Lambda(1520)0' == GABSID ) & ( GNINTREE( isGoodKaon, 1 ) > 0 ) & ( GNINTREE( isGoodP, 1 ) > 0 ) )"
+#                          , "isGoodB        = ( ( 'B0' == GABSID ) & ( GNINTREE( isGoodL0, 1 ) > 0 ) )" ]
+# EndInsertPythonCode
+#
+Alias      MyLambda(1520)0       Lambda(1520)0
+Alias      Myanti-Lambda(1520)0  anti-Lambda(1520)0
+ChargeConj MyLambda(1520)0       Myanti-Lambda(1520)0
+#
+Alias  MyH_30     A0
+Alias  Myanti-H_30    A0
+ChargeConj MyH_30   Myanti-H_30
+#
+Decay B0sig
+    1.000   MyLambda(1520)0   MyH_30    PHSP;
+Enddecay
+CDecay anti-B0sig
+#
+Decay MyLambda(1520)0
+    1.000       p+          K-              PHSP;
+Enddecay
+CDecay Myanti-Lambda(1520)0
+#
+End
--- a/dkfiles/B0_PsiDMLambda1520,Lambda1520_pK=FullGenEvtCut,mPsiDM=940MeV.dec
+++ b/dkfiles/B0_PsiDMLambda1520,Lambda1520_pK=FullGenEvtCut,mPsiDM=940MeV.dec
+# EventType: 11502000
+#
+# Descriptor: [B0 -> (Lambda(1520)0 -> p+ K-) H_30 ]cc
+#
+# NickName: B0_PsiDMLambda1520,Lambda1520_pK=FullGenEvtCut,mPsiDM=940MeV
+#
+# Cuts: None
+# FullEventCuts: LoKi::FullGenEventCut/B0toDM
+#
+# Documentation:
+#   Decay a B0 to a Lambda(1520)0 -> pK and a redefined H_30 for our need, acting the latter as stable Dark Matter candidate. The mass of the H_30 is 0.94 GeV.
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# CPUTime: < 1 min
+# Responsible: Saul Lopez 
+# Email: saul.lopez.solino@cern.ch
+# Date: 20210930
+#
+#
+# InsertPythonCode:
+# from Configurables import LHCb__ParticlePropertySvc, LoKi__GenCutTool
+# LHCb__ParticlePropertySvc().Particles = [
+# "H_30     89       36      0.0     0.940000        1.000000e+16    A0      36      0.00"
+# ]
+# ## Generator level cuts:
+# from Configurables import LoKi__FullGenEventCut
+# Generation().addTool( LoKi__FullGenEventCut, "B0toDM" )
+# tracksInAcc = Generation().B0toDM
+# tracksInAcc.Code = " count ( isGoodB ) > 0 "
+# ### - HepMC::IteratorRange::descendants   4
+# tracksInAcc.Preambulo += [ "from GaudiKernel.SystemOfUnits import GeV, mrad"
+#                          , "inAcc = in_range(1.9, GETA, 5.0)" 
+#                          , "isGoodKaon     = ( ( GPT > 0.50*GeV ) & inAcc & ( 'K+' == GABSID ) )"
+#                          , "isGoodP        = ( ( GPT > 0.50*GeV ) & inAcc & ( 'p+' == GABSID ) )"
+#                          , "isGoodL0       = ( ( 'Lambda(1520)0' == GABSID ) & ( GNINTREE( isGoodKaon, 1 ) > 0 ) & ( GNINTREE( isGoodP, 1 ) > 0 ) )"
+#                          , "isGoodB        = ( ( 'B0' == GABSID ) & ( GNINTREE( isGoodL0, 1 ) > 0 ) )" ]
+# EndInsertPythonCode
+#
+Alias      MyLambda(1520)0       Lambda(1520)0
+Alias      Myanti-Lambda(1520)0  anti-Lambda(1520)0
+ChargeConj MyLambda(1520)0       Myanti-Lambda(1520)0
+#
+Alias  MyH_30     A0
+Alias  Myanti-H_30    A0
+ChargeConj MyH_30   Myanti-H_30
+#
+Decay B0sig
+    1.000   MyLambda(1520)0   MyH_30    PHSP;
+Enddecay
+CDecay anti-B0sig
+#
+Decay MyLambda(1520)0
+    1.000       p+          K-              PHSP;
+Enddecay
+CDecay Myanti-Lambda(1520)0
+#
+End
--- a/dkfiles/Bplus_PsiDMLambdac2595,Lambdac2595_pipiLambdacplus,Lambdacplus_pKpi=FullGenEvtCut,mPsiDM=1500MeV.dec
+++ b/dkfiles/Bplus_PsiDMLambdac2595,Lambdac2595_pipiLambdacplus,Lambdacplus_pKpi=FullGenEvtCut,mPsiDM=1500MeV.dec
+# EventType: 12565011
+#
+# Descriptor: [B+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> p+ K- pi+) pi+ pi-) H_30 ]cc
+#
+# NickName: Bplus_PsiDMLambdac2595,Lambdac2595_pipiLambdacplus,Lambdacplus_pKpi=FullGenEvtCut,mPsiDM=1500MeV
+#
+# Cuts: None
+# FullEventCuts: LoKi::FullGenEventCut/BplustoDM
+#
+# Documentation:
+#   Decay a B+ to a Lambda_c(2595)+ -> pipiLambda_c (Lambda_c -> pKpi) and a redefined H_30 for our need, acting the latter as stable Dark Matter candidate of 1.50 GeV.
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# CPUTime: < 1 min
+# Responsible: Saul Lopez
+# Email: saul.lopez.solino@cern.ch
+# Date: 20210930
+#
+#
+# InsertPythonCode:
+# from Configurables import LHCb__ParticlePropertySvc, LoKi__GenCutTool
+# LHCb__ParticlePropertySvc().Particles = [
+# "H_30     89       36      0.0     1.500000        1.000000e+16    A0      36      0.00"
+# ]
+# ## Generator level cuts:
+# from Configurables import LoKi__FullGenEventCut
+# Generation().addTool( LoKi__FullGenEventCut, "BplustoDM" )
+# tracksInAcc = Generation().BplustoDM
+# tracksInAcc.Code = " count ( isGoodB ) > 0 "
+# ### - HepMC::IteratorRange::descendants   4
+# tracksInAcc.Preambulo += [ "from GaudiKernel.SystemOfUnits import GeV, mrad"
+#                          , "inAcc = in_range(1.9, GETA, 5.0)" 
+#                          , "isGoodKaon     = ( ( GPT > 0.25*GeV ) & inAcc & ( 'K+' == GABSID ) )"
+#                          , "isGoodPiC      = ( ( GPT > 0.25*GeV ) & inAcc & ( 'pi+' == GABSID ) )"
+#                          , "isGoodPi2595   = ( ( GPT > 0.10*GeV ) & inAcc & ( 'pi+' == GABSID ) )"
+#                          , "isGoodP        = ( ( GPT > 0.25*GeV ) & inAcc & ( 'p+' == GABSID ) )"
+#                          , "isGoodLc       = ( ( 'Lambda_c+' == GABSID ) & ( GNINTREE( isGoodKaon, 1 ) > 0 ) & ( GNINTREE( isGoodP, 1 ) > 0 ) & ( GNINTREE( isGoodPiC, 1 ) > 0 ) )"
+#                          , "isGoodLc2595   = ( ( 'Lambda_c(2595)+' == GABSID ) & ( GNINTREE( isGoodLc, 1 ) > 0) & ( GNINTREE ( isGoodPi2595, 1 ) > 1 ) )"
+#                          , "isGoodB        = ( ( 'B+' == GABSID ) & ( GNINTREE( isGoodLc2595, 1 ) > 0 ) )" ]
+# EndInsertPythonCode
+#
+Alias      MyLambda_c(2593)+       Lambda_c(2593)+
+Alias      Myanti-Lambda_c(2593)-       anti-Lambda_c(2593)-
+ChargeConj MyLambda_c(2593)+       Myanti-Lambda_c(2593)-
+#
+Alias       MyLambda_c+        Lambda_c+
+Alias  Myanti-Lambda_c-        anti-Lambda_c-
+ChargeConj      MyLambda_c+    Myanti-Lambda_c-
+#
+Alias  MyH_30     A0
+Alias  Myanti-H_30    A0
+ChargeConj MyH_30   Myanti-H_30
+#
+Decay B+sig
+    1.000   MyLambda_c(2593)+   MyH_30    PHSP;
+Enddecay
+CDecay B-sig
+#
+Decay MyLambda_c(2593)+
+    1.000      MyLambda_c+ pi+ pi-             PHSP;
+Enddecay
+CDecay Myanti-Lambda_c(2593)-
+#
+Decay MyLambda_c+
+    1.000       p+  K-  pi+      PHSP;
+Enddecay
+CDecay Myanti-Lambda_c-
+#
+End
--- a/dkfiles/Bplus_PsiDMLambdac2595,Lambdac2595_pipiLambdacplus,Lambdacplus_pKpi=FullGenEvtCut,mPsiDM=2000MeV.dec
+++ b/dkfiles/Bplus_PsiDMLambdac2595,Lambdac2595_pipiLambdacplus,Lambdacplus_pKpi=FullGenEvtCut,mPsiDM=2000MeV.dec
+# EventType: 12565012
+#
+# Descriptor: [B+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> p+ K- pi+) pi+ pi-) H_30 ]cc
+#
+# NickName: Bplus_PsiDMLambdac2595,Lambdac2595_pipiLambdacplus,Lambdacplus_pKpi=FullGenEvtCut,mPsiDM=2000MeV
+#
+# Cuts: None
+# FullEventCuts: LoKi::FullGenEventCut/BplustoDM
+#
+# Documentation:
+#   Decay a B+ to a Lambda_c(2595)+ -> pipiLambda_c (Lambda_c -> pKpi) and a redefined H_30 for our need, acting the latter as stable Dark Matter candidate of 2.0 GeV.
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# CPUTime: < 1 min
+# Responsible: Saul Lopez
+# Email: saul.lopez.solino@cern.ch
+# Date: 20210930
+#
+#
+# InsertPythonCode:
+# from Configurables import LHCb__ParticlePropertySvc, LoKi__GenCutTool
+# LHCb__ParticlePropertySvc().Particles = [
+# "H_30     89       36      0.0     2.000000        1.000000e+16    A0      36      0.00"
+# ]
+# ## Generator level cuts:
+# from Configurables import LoKi__FullGenEventCut
+# Generation().addTool( LoKi__FullGenEventCut, "BplustoDM" )
+# tracksInAcc = Generation().BplustoDM
+# tracksInAcc.Code = " count ( isGoodB ) > 0 "
+# ### - HepMC::IteratorRange::descendants   4
+# tracksInAcc.Preambulo += [ "from GaudiKernel.SystemOfUnits import GeV, mrad"
+#                          , "inAcc = in_range(1.9, GETA, 5.0)" 
+#                          , "isGoodKaon     = ( ( GPT > 0.25*GeV ) & inAcc & ( 'K+' == GABSID ) )"
+#                          , "isGoodPiC      = ( ( GPT > 0.25*GeV ) & inAcc & ( 'pi+' == GABSID ) )"
+#                          , "isGoodPi2595   = ( ( GPT > 0.10*GeV ) & inAcc & ( 'pi+' == GABSID ) )"
+#                          , "isGoodP        = ( ( GPT > 0.25*GeV ) & inAcc & ( 'p+' == GABSID ) )"
+#                          , "isGoodLc       = ( ( 'Lambda_c+' == GABSID ) & ( GNINTREE( isGoodKaon, 1 ) > 0 ) & ( GNINTREE( isGoodP, 1 ) > 0 ) & ( GNINTREE( isGoodPiC, 1 ) > 0 ) )"
+#                          , "isGoodLc2595   = ( ( 'Lambda_c(2595)+' == GABSID ) & ( GNINTREE( isGoodLc, 1 ) > 0) & ( GNINTREE ( isGoodPi2595, 1 ) > 1 ) )"
+#                          , "isGoodB        = ( ( 'B+' == GABSID ) & ( GNINTREE( isGoodLc2595, 1 ) > 0 ) )" ]
+# EndInsertPythonCode
+#
+Alias      MyLambda_c(2593)+       Lambda_c(2593)+
+Alias      Myanti-Lambda_c(2593)-       anti-Lambda_c(2593)-
+ChargeConj MyLambda_c(2593)+       Myanti-Lambda_c(2593)-
+#
+Alias       MyLambda_c+        Lambda_c+
+Alias  Myanti-Lambda_c-        anti-Lambda_c-
+ChargeConj      MyLambda_c+    Myanti-Lambda_c-
+#
+Alias  MyH_30     A0
+Alias  Myanti-H_30    A0
+ChargeConj MyH_30   Myanti-H_30
+#
+Decay B+sig
+    1.000   MyLambda_c(2593)+   MyH_30    PHSP;
+Enddecay
+CDecay B-sig
+#
+Decay MyLambda_c(2593)+
+    1.000      MyLambda_c+ pi+ pi-             PHSP;
+Enddecay
+CDecay Myanti-Lambda_c(2593)-
+#
+Decay MyLambda_c+
+    1.000       p+  K-  pi+      PHSP;
+Enddecay
+CDecay Myanti-Lambda_c-
+#
+End
--- a/dkfiles/Bplus_PsiDMLambdac2595,Lambdac2595_pipiLambdacplus,Lambdacplus_pKpi=FullGenEvtCut,mPsiDM=2400MeV.dec
+++ b/dkfiles/Bplus_PsiDMLambdac2595,Lambdac2595_pipiLambdacplus,Lambdacplus_pKpi=FullGenEvtCut,mPsiDM=2400MeV.dec
+# EventType: 12565013
+#
+# Descriptor: [B+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> p+ K- pi+) pi+ pi-) H_30 ]cc
+#
+# NickName: Bplus_PsiDMLambdac2595,Lambdac2595_pipiLambdacplus,Lambdacplus_pKpi=FullGenEvtCut,mPsiDM=2400MeV
+#
+# Cuts: None
+# FullEventCuts: LoKi::FullGenEventCut/BplustoDM
+#
+# Documentation:
+#   Decay a B+ to a Lambda_c(2595)+ -> pipiLambda_c (Lambda_c -> pKpi) and a redefined H_30 for our need, acting the latter as stable Dark Matter candidate of 2.4 GeV.
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# CPUTime: < 1 min
+# Responsible: Saul Lopez
+# Email: saul.lopez.solino@cern.ch
+# Date: 20210930
+#
+#
+# InsertPythonCode:
+# from Configurables import LHCb__ParticlePropertySvc, LoKi__GenCutTool
+# LHCb__ParticlePropertySvc().Particles = [
+# "H_30     89       36      0.0     2.400000        1.000000e+16    A0      36      0.00"
+# ]
+# ## Generator level cuts:
+# from Configurables import LoKi__FullGenEventCut
+# Generation().addTool( LoKi__FullGenEventCut, "BplustoDM" )
+# tracksInAcc = Generation().BplustoDM
+# tracksInAcc.Code = " count ( isGoodB ) > 0 "
+# ### - HepMC::IteratorRange::descendants   4
+# tracksInAcc.Preambulo += [ "from GaudiKernel.SystemOfUnits import GeV, mrad"
+#                          , "inAcc = in_range(1.9, GETA, 5.0)" 
+#                          , "isGoodKaon     = ( ( GPT > 0.25*GeV ) & inAcc & ( 'K+' == GABSID ) )"
+#                          , "isGoodPiC      = ( ( GPT > 0.25*GeV ) & inAcc & ( 'pi+' == GABSID ) )"
+#                          , "isGoodPi2595   = ( ( GPT > 0.10*GeV ) & inAcc & ( 'pi+' == GABSID ) )"
+#                          , "isGoodP        = ( ( GPT > 0.25*GeV ) & inAcc & ( 'p+' == GABSID ) )"
+#                          , "isGoodLc       = ( ( 'Lambda_c+' == GABSID ) & ( GNINTREE( isGoodKaon, 1 ) > 0 ) & ( GNINTREE( isGoodP, 1 ) > 0 ) & ( GNINTREE( isGoodPiC, 1 ) > 0 ) )"
+#                          , "isGoodLc2595   = ( ( 'Lambda_c(2595)+' == GABSID ) & ( GNINTREE( isGoodLc, 1 ) > 0) & ( GNINTREE ( isGoodPi2595, 1 ) > 1 ) )"
+#                          , "isGoodB        = ( ( 'B+' == GABSID ) & ( GNINTREE( isGoodLc2595, 1 ) > 0 ) )" ]
+# EndInsertPythonCode
+#
+Alias      MyLambda_c(2593)+       Lambda_c(2593)+
+Alias      Myanti-Lambda_c(2593)-       anti-Lambda_c(2593)-
+ChargeConj MyLambda_c(2593)+       Myanti-Lambda_c(2593)-
+#
+Alias       MyLambda_c+        Lambda_c+
+Alias  Myanti-Lambda_c-        anti-Lambda_c-
+ChargeConj      MyLambda_c+    Myanti-Lambda_c-
+#
+Alias  MyH_30     A0
+Alias  Myanti-H_30    A0
+ChargeConj MyH_30   Myanti-H_30
+#
+Decay B+sig
+    1.000   MyLambda_c(2593)+   MyH_30    PHSP;
+Enddecay
+CDecay B-sig
+#
+Decay MyLambda_c(2593)+
+    1.000      MyLambda_c+ pi+ pi-             PHSP;
+Enddecay
+CDecay Myanti-Lambda_c(2593)-
+#
+Decay MyLambda_c+
+    1.000       p+  K-  pi+      PHSP;
+Enddecay
+CDecay Myanti-Lambda_c-
+#
+End
--- a/dkfiles/Bplus_PsiDMLambdac2595,Lambdac2595_pipiLambdacplus,Lambdacplus_pKpi=FullGenEvtCut,mPsiDM=940MeV.dec
+++ b/dkfiles/Bplus_PsiDMLambdac2595,Lambdac2595_pipiLambdacplus,Lambdacplus_pKpi=FullGenEvtCut,mPsiDM=940MeV.dec
+# EventType: 12565010
+#
+# Descriptor: [B+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> p+ K- pi+) pi+ pi-) H_30 ]cc
+#
+# NickName: Bplus_PsiDMLambdac2595,Lambdac2595_pipiLambdacplus,Lambdacplus_pKpi=FullGenEvtCut,mPsiDM=940MeV
+#
+# Cuts: None
+# FullEventCuts: LoKi::FullGenEventCut/BplustoDM
+#
+# Documentation:
+#   Decay a B+ to a Lambda_c(2595)+ -> pipiLambda_c (Lambda_c -> pKpi) and a redefined H_30 for our need, acting the latter as stable Dark Matter candidate of 0.94 GeV.
+# EndDocumentation