ParticleCombiner has issue in intermediate decay with gamma
Trying to reconstruct "[B+ -> (chi_c1(1P) -> (J/psi(1S) -> mu+ mu-) gamma) K+]cc"
with b_to_jpsix0, while information of chi_c1(1P) decay is missing in the hlt2-filtered dst file, the whole decay chain in this dst file is showed below:
Print_B.PrintDecayTreeTool INFO <-------------------------------------------- Particle -------------------------------------------->
Print_B.PrintDecayTreeTool INFO Name E M P Pt phi Vz P(C/K) PP(C/K)
Print_B.PrintDecayTreeTool INFO MeV MeV MeV MeV mrad mm
Print_B.PrintDecayTreeTool INFO B- 85569.28 6083.77 85352.74 7795.29 229.73 22.48 0/0 N/A
Print_B.PrintDecayTreeTool INFO +-->J/psi(1S) 60093.32 3090.15 60013.82 6806.93 247.12 22.48 1/0 N/A
Print_B.PrintDecayTreeTool INFO |+-->mu+ 24601.15 105.66 24600.93 3183.85 743.21 19.28 2/14 3/69
Print_B.PrintDecayTreeTool INFO |+-->mu- 35483.87 105.66 35483.72 4282.61 -114.58 19.43 2/7 3/39
Print_B.PrintDecayTreeTool INFO +-->gamma 3205.87 0.00 3205.87 685.75 -1261.22 0.00 4/3 5/19
Print_B.PrintDecayTreeTool INFO +-->K- 22270.16 493.68 22264.68 1091.78 773.22 14.22 6/33 3/33
Print_B.PrintDecayTreeTool INFO
Print_B.PrintDecayTreeTool INFO Used TES locations :-
Print_B.PrintDecayTreeTool INFO 1 = '/Event/HLT2/B2CC_Jpsi2MuMu_Filter/OutputParticles'
Print_B.PrintDecayTreeTool INFO 6 = '/Event/HLT2/FunctionalParticleMaker#1/Particles'
Print_B.PrintDecayTreeTool INFO 2 = '/Event/HLT2/FunctionalParticleMaker/Particles'
Print_B.PrintDecayTreeTool INFO 0 = '/Event/HLT2/Hlt2B2CC_BuToChicKPlus_Line/Particles'
Print_B.PrintDecayTreeTool INFO 4 = '/Event/HLT2/LHCb__Phys__ParticleMakers__PhotonMaker/Particles'
Print_B.PrintDecayTreeTool INFO 3 = '/Event/HLT2/Rec/ProtoP/Charged'
Print_B.PrintDecayTreeTool INFO 5 = '/Event/HLT2/Rec/ProtoP/Neutrals'
I've also done a test by making a "chi_c1(1P)" from selected jpsi (defined in basic_builder) and kaon (line defined in b_to_jpsix0.py and hlt2_b2cc.py), named as Hlt2B2CC_BuToChicKPlusTest_Line
(decay descriptor: "[B+ -> (chi_c1(1P) -> (J/psi(1S) -> mu+ mu-) K+) gamma]cc"
). In this case, the "chi_c1(1P)" information and the whole decay topology are stored correctly in the dst file, like this:
Print_Btest.PrintDecayTreeTool INFO <-------------------------------------------- Particle -------------------------------------------->
Print_Btest.PrintDecayTreeTool INFO Name E M P Pt phi Vz P(C/K) PP(C/K)
Print_Btest.PrintDecayTreeTool INFO MeV MeV MeV MeV mrad mm
Print_Btest.PrintDecayTreeTool INFO B+ 91531.53 6961.79 91266.40 7417.83 457.19 22.48 7/22 N/A
Print_Btest.PrintDecayTreeTool INFO +-->chi_c1(1P) 82363.41 4644.99 82232.33 7770.71 317.81 22.48 8/1 N/A
Print_Btest.PrintDecayTreeTool INFO |+-->J/psi(1S) 60093.32 3090.15 60013.82 6806.93 247.12 22.48 1/0 N/A
Print_Btest.PrintDecayTreeTool INFO ||+-->mu+ 24601.15 105.66 24600.93 3183.85 743.21 19.28 2/14 3/69
Print_Btest.PrintDecayTreeTool INFO ||+-->mu- 35483.87 105.66 35483.72 4282.61 -114.58 19.43 2/7 3/39
Print_Btest.PrintDecayTreeTool INFO |+-->K- 22270.16 493.68 22264.68 1091.78 773.22 14.22 6/33 3/33
Print_Btest.PrintDecayTreeTool INFO +-->gamma 9168.13 0.00 9168.13 1114.74 2279.59 0.00 4/29 5/50
Print_Btest.PrintDecayTreeTool INFO
Print_Btest.PrintDecayTreeTool INFO Used TES locations :-
Print_Btest.PrintDecayTreeTool INFO 8 = '/Event/HLT2/B2CC_Chic2JpsiG_Filter#1/OutputParticles'
Print_Btest.PrintDecayTreeTool INFO 1 = '/Event/HLT2/B2CC_Jpsi2MuMu_Filter/OutputParticles'
Print_Btest.PrintDecayTreeTool INFO 6 = '/Event/HLT2/FunctionalParticleMaker#1/Particles'
Print_Btest.PrintDecayTreeTool INFO 2 = '/Event/HLT2/FunctionalParticleMaker/Particles'
Print_Btest.PrintDecayTreeTool INFO 7 = '/Event/HLT2/Hlt2B2CC_BuToChicKPlusTest_Line/Particles'
Print_Btest.PrintDecayTreeTool INFO 4 = '/Event/HLT2/LHCb__Phys__ParticleMakers__PhotonMaker/Particles'
Print_Btest.PrintDecayTreeTool INFO 3 = '/Event/HLT2/Rec/ProtoP/Charged'
Print_Btest.PrintDecayTreeTool INFO 5 = '/Event/HLT2/Rec/ProtoP/Neutrals'
It seems that there's still a problem with ParticleCombiner when using it to reconstruct the intermediate state with gamma as one of the daughter particles.
(The Moore and DaVinci scripts I used to do the test and related log files (stdout.davincilog stdout.moorelog ) are attached, as well as the tck file of the output dst file Hlt2B2JPSIX0LINE_CHIC_MU_test.tck.json.)