Pulsar Analysis lomb-scargle periodogram

I would like to create a graph showing the relationship between frequency and power using the Lomb-Scargle periodogram method with the Crab Pulsar data observed by the Suzaku satellite. The data I will use is from the observation on 2010-02-23 01:01:31, which was downloaded from Heasarc (ae104001070hxd_0_pinno_cl.evt).I don’t know how to load this data using SciPy, so I would like to know. Since the Crab Pulsar has a period of about 33 ms, I expect a peak to appear around 29.7 Hz. I would appreciate your help in writing this code.By the way, the time derivative of the frequency (fdot) is approximately -3.7×10^{-15}.