Skip to content

Add additional D0 and anti-D0 decay modes to Prompt D0 anti-D0 spectroscopy searches

Paras Naik requested to merge promptDzDzb_BandQ_pnaik into master

Add additional D0 and anti-D0 decay modes to Prompt D0 anti-D0 spectroscopy searches. Also as more data is collected the quantum correlations coming from C=+/-1 resonances can be studied.

A line already exists that includes prompt D0 anti-D0 but only when both decay to the Cabibbo-favoured Kpi mode. This line also includes all sorts of other prompt double charm. So the decision was made to create new lines for the purpose of adding additional D0 and anti-D0 decay modes. All new lines will have the same other properties (run in HLT2, persistReco = True, etc...) as the existing double charm line.

Note that from the trigger perspective since e.g. D0 -> KK is self-conjugate it doesn't make sense to loop over anti-D0. Thus for these lines the decay descriptor is [psi(3770) -> D0 D0]cc instead of [psi(3770) -> D~0 D0]cc. I am not sure the cc bit is really necessary, but it exists in the existing doublecharm lines...

Additions so far:

  • Add ten new D0 anti-D0 lines based on topology of number of particles and whether the Ks is reconstructed LL or DD
    • Hlt2BandQ_DoubleCharm_D0ToHH: first D0 -> HH second D0 -> HH
      • [note, the D0 -> Kpi, D0 -> Kpi counted in the existing line would also be counted here... not sure if this will be a minor problem and/or if this decay dominates the rates]
    • Hlt2BandQ_DoubleCharm_D0ToKsLLHH_D0ToHH: first D0 -> KsHH (LL Ks) second D0 -> HH
    • Hlt2BandQ_DoubleCharm_D0ToKsDDHH_D0ToHH: first D0 -> KsHH (DD Ks) second D0 -> HH
    • Hlt2BandQ_DoubleCharm_D0ToKsLLHH: first D0 -> KsHH (LL Ks) second D0 -> KsHH (LL Ks)
    • Hlt2BandQ_DoubleCharm_D0ToKsLLHH_D0ToKsDDHH: first D0 -> KsHH (LL Ks) second D0 -> KsHH (DD Ks)
    • Hlt2BandQ_DoubleCharm_D0ToKsDDHH: first D0 -> KsHH (DD Ks) second D0 -> KsHH (DD Ks)
    • Hlt2BandQ_DoubleCharm_D0ToHHHH_D0ToHH: first D0 -> HHHH second D0 -> HH
    • Hlt2BandQ_DoubleCharm_D0ToHHHH_D0ToKsLLHH: first D0 -> HHHH second D0 -> KsHH (LL Ks)
    • Hlt2BandQ_DoubleCharm_D0ToHHHH_D0ToKsDDHH: first D0 -> HHHH second D0 -> KsHH (DD Ks)
    • Hlt2BandQ_DoubleCharm_D0ToHHHH: first D0 -> HHHH second D0 -> HHHH
  • Based on work done in B2OC, I would expect an error putting in the above lines in the case that the input D0 decay modes are different: ERROR NBodyCombiner<name> : Got contiguous children with the same PIDs, but the corresponding algorithm inputs were not the same. If this is intentional, please set AllowDiffInputsForSameIDChildren = True.
    • To get past this, I added the capability to use AllowDiffInputsForSameIDChildren in ParticleCombiner by adding an option allowDiffInputsForSameIDChildren (default = False) to _make_doublecharm and allowDiffInputsForSameIDChildren_dc (default = False) to make_doublecharm in doublecharm.py

Current work status:

  • Line writing
  • determine retentions
  • determine retentions with more events / DaVinci (running)

Retentions

  • New test with ~170k events and DaVinci
[pnaik@lxplus944 stack_Moore_promptDzDzb_BandQ_pnaik]$  cat hlt_eff_checker_hlt2_rate_example.log | grep 'LAZY_AND:' | grep 'WithOutput' | grep 'DoubleCharm'
    LAZY_AND: Hlt2BandQ_DoubleCharmOppositeSignDecisionWithOutput                                                       #=171490  Sum=26          Eff=|(0.01516123 +- 0.00297314)%|
    LAZY_AND: Hlt2BandQ_DoubleCharmSameSignDecisionWithOutput                                                           #=171490  Sum=12          Eff=|(0.006997493 +- 0.00201993)%|
    LAZY_AND: Hlt2BandQ_DoubleCharm_D0ToHHDecisionWithOutput                                                            #=171490  Sum=9           Eff=|(0.005248119 +- 0.00174933)%|
    LAZY_AND: Hlt2BandQ_DoubleCharm_D0ToHHHHDecisionWithOutput                                                          #=171490  Sum=3           Eff=|(0.001749373 +- 0.00100999)%|
    LAZY_AND: Hlt2BandQ_DoubleCharm_D0ToHHHH_D0ToHHDecisionWithOutput                                                   #=171490  Sum=9           Eff=|(0.005248119 +- 0.00174933)%|
    LAZY_AND: Hlt2BandQ_DoubleCharm_D0ToHHHH_D0ToKsDDHHDecisionWithOutput                                               #=171490  Sum=0           Eff=|( 0.000000 +- 0.00000 )%|
    LAZY_AND: Hlt2BandQ_DoubleCharm_D0ToHHHH_D0ToKsLLHHDecisionWithOutput                                               #=171490  Sum=15          Eff=|(0.008746866 +- 0.00225833)%|
    LAZY_AND: Hlt2BandQ_DoubleCharm_D0ToKsDDHHDecisionWithOutput                                                        #=171490  Sum=0           Eff=|( 0.000000 +- 0.00000 )%|
    LAZY_AND: Hlt2BandQ_DoubleCharm_D0ToKsDDHH_D0ToHHDecisionWithOutput                                                 #=171490  Sum=0           Eff=|( 0.000000 +- 0.00000 )%|
    LAZY_AND: Hlt2BandQ_DoubleCharm_D0ToKsLLHHDecisionWithOutput                                                        #=171490  Sum=22          Eff=|(0.01282874 +- 0.00273492)%|
    LAZY_AND: Hlt2BandQ_DoubleCharm_D0ToKsLLHH_D0ToHHDecisionWithOutput                                                 #=171490  Sum=32          Eff=|(0.01865998 +- 0.00329834)%|
    LAZY_AND: Hlt2BandQ_DoubleCharm_D0ToKsLLHH_D0ToKsDDHHDecisionWithOutput                                             #=171490  Sum=0           Eff=|( 0.000000 +- 0.00000 )%|
  • Old test with ~50k events
[pnaik@lxplus914 stack_Moore_promptDzDzb_BandQ_pnaik]$ Moore/run gaudirun.py option.py > bandq_all_includingQC.log &
[pnaik@lxplus914 stack_Moore_promptDzDzb_BandQ_pnaik]$ cat bandq_all_includingQC.log | grep 'LAZY_AND:' | grep 'WithOutput' | grep 'DoubleCharm'
   LAZY_AND: Hlt2BandQ_DoubleCharmOppositeSignDecisionWithOutput                                                       #=50541   Sum=7           Eff=|(0.01385014 +- 0.00523450)%|
   LAZY_AND: Hlt2BandQ_DoubleCharmSameSignDecisionWithOutput                                                           #=50541   Sum=5           Eff=|(0.009892958 +- 0.00442405)%|
   LAZY_AND: Hlt2BandQ_DoubleCharm_D0ToHHDecisionWithOutput                                                            #=50541   Sum=2           Eff=|(0.003957183 +- 0.00279810)%|
   LAZY_AND: Hlt2BandQ_DoubleCharm_D0ToHHHHDecisionWithOutput                                                          #=50541   Sum=1           Eff=|(0.001978592 +- 0.00197857)%|
   LAZY_AND: Hlt2BandQ_DoubleCharm_D0ToHHHH_D0ToHHDecisionWithOutput                                                   #=50541   Sum=2           Eff=|(0.003957183 +- 0.00279810)%|
   LAZY_AND: Hlt2BandQ_DoubleCharm_D0ToHHHH_D0ToKsDDHHDecisionWithOutput                                               #=50541   Sum=0           Eff=|( 0.000000 +- 0.00000 )%|
   LAZY_AND: Hlt2BandQ_DoubleCharm_D0ToHHHH_D0ToKsLLHHDecisionWithOutput                                               #=50541   Sum=3           Eff=|(0.005935775 +- 0.00342692)%|
   LAZY_AND: Hlt2BandQ_DoubleCharm_D0ToKsDDHHDecisionWithOutput                                                        #=50541   Sum=0           Eff=|( 0.000000 +- 0.00000 )%|
   LAZY_AND: Hlt2BandQ_DoubleCharm_D0ToKsDDHH_D0ToHHDecisionWithOutput                                                 #=50541   Sum=0           Eff=|( 0.000000 +- 0.00000 )%|
   LAZY_AND: Hlt2BandQ_DoubleCharm_D0ToKsLLHHDecisionWithOutput                                                        #=50541   Sum=4           Eff=|(0.007914367 +- 0.00395703)%|
   LAZY_AND: Hlt2BandQ_DoubleCharm_D0ToKsLLHH_D0ToHHDecisionWithOutput                                                 #=50541   Sum=8           Eff=|(0.01582873 +- 0.00559586)%|
   LAZY_AND: Hlt2BandQ_DoubleCharm_D0ToKsLLHH_D0ToKsDDHHDecisionWithOutput                                             #=50541   Sum=0           Eff=|( 0.000000 +- 0.00000 )%|
Edited by Paras Naik

Merge request reports