Implement synthetic flare injection and recovery in a new module

[ekaterinailin]

What needs to be created or improved?

• [x] create a function that injects flares generated by a semi-empirical model into a real light curve, avoiding true flares and return a new FlareLightCurve.
• [x] create a wrapper that repeats the above process n times with a given flare frequency in a physically motivated grid and return recovered energies and recovery probabilities as a function of ED.

Can you provide an example?

def inject_flares(lc, flare_frequency=0.1):
    [...]
    return new_lc

def synthetic_flare_analysis(lc, n=1000):
    [...]
    return recovery_probability, ed_fraction

What is the goal / expected behaviour?

The module includes everything to characterize a given light curve with respect to its flare detection performance.