Pamenend/bmesogenesis

Added 1 DaVinci script with a function for each case:

• Sb+ -> pi + (Lb -> DM + K + (D -> pi + pi + K)), on sim
• Sb- -> pi + (Lb -> DM + K + (D -> pi + pi + K)), on sim
• Sb+ (and Sb-) -> pi + (Lb -> K + (D -> pi + pi + K)), on data

and 3 python scripts calling the function for each case.

Different tuples are needed for Sb+ and Sb- since they are not charge conjugate of each other. Input files correspond to:

• 5 simulation samples for different masses of Dark Matter in the Sb+ channel, both mags.
• 4 simulation samples for different masses of Dark Matter in the Sbstar+ channel, both mags. In this case we only use 4 samples, since event type 16565023 didnt finish sim yet.
• 5 simulation samples for different masses of Dark Matter in the Sb- channel, both mags.
• 5 simulation samples for different masses of Dark Matter in the Sbstar- channel, both mags.
• Data from 2018 collisions and Stripping34r0p3 and Sb+ (or Sb-) CombineParticles, both mags.