timing analysis , other than being an important diagnostic tool by itself ( see * ? ? ?
* for a review ) , is particularly powerful for dissecting the inner accretion regions in black hole binaries and agn when combined with spectral modeling .
good examples of this combination in binaries are studies of the variability of a specific spectral component ( e.g. the iron k@xmath3 line and the reflection component , * ? ? ?
* ; * ? ? ?
* ) , the study of the time lags between different energy bands and the comparison with the expected behavior of a given spectral component ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) , the study of the covariance of the signal at multiple energies and the variability spectrum ( e.g. * ? ? ?
* ; * ? ? ?
a review of reverberation lag measurements in binaries and agn can be found in @xcite .    the _ nuclear spectroscopic telescope array _ ( _ nustar _ ; @xcite ) mission has deployed the first hard x - ray ( above 10kev ) focusing satellite .
it features two telescopes , focusing x - rays between 3 and 79kev onto two identical focal planes ( usually called focal plane modules a and b , or fpma and fpmb ) .
it has a field of view ( fov ) of 12@xmath412and an angular resolution of 18fwhm ; ( 58hpd ) .
these features , together with good spectral resolution in the iron k band ( @xmath00.4kev@6kev ) , have motivated a large number of observations of x - ray binaries .
a particularly interesting result was the accurate modeling of the reflection component in several black hole ( bh ) binaries @xcite resulting in improved constraints on black hole spin in several systems .
such measurements have also been made by _
nustar_for active galactic nuclei ( agn ; e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
_ nustar _ , _ xmm _ and _ suzaku _ have recently measured reverberation lags in the agn mcg-5 - 23 - 16 in both the iron line and compton reflection hump together for the first time @xcite .    from the point of view of timing , _ nustar_has an advantage over other imaging satellites : the time resolution is 10@xmath5 , and so one can in principle study variability over the whole range of interesting frequencies in accreting systems without switching to an observing mode with decreased spectral resolution .
the satellite has in fact been used to perform timing analysis on selected targets with very good results : for example , it was used to study a number of accreting and rotation - powered pulsars , permitting the detection of variable cyclotron resonant scattering features @xcite and the discovery of strong hard time lags in the emission of another accreting pulsar ( gs 0834 - 430 ; @xcite ) . _ nustar_was also used to detect the pulsations from a magnetar near the galactic center @xcite and to carry out a detailed timing study of this serendipitous source @xcite .
these pulsars are mostly slow rotators , with pulse periods above 1s .
aperiodic timing has been measured in faint sources ( e.g. * ? ? ?
* ) , or in bright sources but only over the low - frequency part of the power spectrum ( e.g. * ? ? ?
* ) .    from the point of view of fast ( @xmath6hz ) timing studies of bright sources , however ,
s default observing mode requires care to analyze properly .
each x - ray incident on a focal plane produces a trigger , and the energy signal is read out immediately .
the read time is dependent on the number of triggered pixels , and so the dead time for each event is slightly variable ( at the few percent level ) with an average value of @xmath7ms .
this dead time , together with additional dead time produced on a separate cadence by housekeeping operations and vetoed events , can be accurately measured and its value over a second is stored in the housekeeping files of the observations .
in addition , the live time since the previous x - ray event is recorded in a column in the event lists ( see section  sec : deadtime for details ) . for bright sources where the fraction of events vetoed by the anti coincidence system is small , light curves can be accurately produced with any time binning ( see madsen et al 2014 , _ submitted _
, for application of this technique to the crab pulsar ) .
however , for timing analysis ( e.g. the production of power density spectra ) , the dead time introduces spurious correlations between event arrival times that are reflected in a distortion of the power spectrum , making it difficult to subtract the poisson noise level above @xmath8hz .
this effect is almost completely negligible for count rates up to @xmath9@xmath10 , but it is very pronounced in observations of galactic x - ray binaries during outbursts , with typical incident count rates @xmath11@xmath10over the full band .    in this paper , we present a technique to avoid spurious noise introduced by dead time , in particular the distortion of the white noise level in the power spectrum .
we show how a slight modification to the methods commonly used for quasi - periodic oscillation ( qpo ) detection and power spectral fitting enables good - quality results at all frequencies .
we use these techniques to perform a basic timing analysis of several _
nustar_observations of bh binaries ( bhbs ) , whose spectral analysis is presented in other papers ( grs1915@xmath2105 , @xcite ; cygx-1 , @xcite ; gx339@xmath14 , frst et al .
_ in prep .
we chose these bhbs as ideal objects with which to demonstrate the timing techniques , since their timing behavior has been extensively studied by other missions and instruments ( see section  sec : analysis for details ) .
the paper is organized as follows : in section  sec : deadtime we discuss details _ nustar _ s dead time , and how it is measured on - board the instrument ; in sections  sec : cpds and  [ sec : simu ] we describe how to use the cross power density spectrum to obtain a proxy of the pds ; in sections  sec : red and  [ sec : analysis ] we use these techniques to analyze the data from two black hole x - ray binaries , cyg x-1 and gx  339 - 4 ; and in section  sec : concl we provide conclusions . in the appendix we present results from a timing analysis of grs  1915 + 105 , which was observed during the on - orbit commissioning period .
this observation suffered from spurious features associated with the cadence of instrument housekeeping functions and so requires a unique analysis not relevant to _ nustar _  science phase observations .
dead time in the _ nustar _
instrument is produced when the focal plane module electronics are busy processing an event , when the shield veto prevents the focal plane from triggering ( or interrupts an event that is being processed ) , and for some regular instrument housekeeping functions ( the housekeeping produces less than 1 ms of dead time per second ) .
_  has two independent telescopes , and the associated focal plane modules have independent processors , so that the dead time in the two modules is uncorrelated .
all of the events processed by the instrument are found in the unfiltered event file .
the cleaned event file is a subset of the events in the unfiltered event file where gtis and additional screening have been applied to remove non - science quality events .    for the production of light curves and flux measurements the _ nustar
_  instrument has two operating modes , one for faint ( @xmath1250 mcrab ) sources , and one that allows highly accurate flux measurements even at very high count rates and on arbitrarily short timescales .
the fraction per second of dead time from events that are not processed , like shield vetoes , as well as that due to housekeeping functions is saved in the housekeeping files on a one - second cadence . for light curves binned in 1  s intervals ( or a integer multiple thereof ) multiplying by the live time stored in the housekeeping files correctly accounts for dead time .
in addition , the live time since the previous x - ray event is stored in the event lists in the  prior " column . in the absence of events vetoed by the active anti - coincidence shield
, this column would accurately reflect the true live time . in the operating mode
used for most observations taken since in - orbit checkout , the events vetoed by the anti - coincidence shield are not telemetered to the ground to minimize data volume . for
this mode adding the prior column does not yield the proper live time since vetoed events are not included .
however , for source count rates significantly greater than the veto rate the error is negligible ( for the crab it is only @xmath13% ; see madsen et al .
2014 , _ submitted _ ) . for sources of intermediate brightness , where the veto rate can not be ignored
, a fully tested and calibrated instrument mode exists that includes vetoed events in the data stream so that fluxes and light curves can be produced with arbitrary accuracy .    for timing analyses that involve power spectral analysis ,
dead time produces systematic effects even if it is perfectly measured .
dead time is classified as either paralyzable or non - paralyzable depending on whether a photon hitting the detector during dead time produces new dead time or not .
the _ nustar _
readout architecture produces non - paralyzable dead time .
its main effect is on recorded count rates : if we have an incident rate of photons @xmath14 , the detected count rate will be approximately @xmath15 ( see , e.g. , * ? ? ? * and references therein ) , where @xmath16is the dead time produced by each event , assuming that it is constant .
but dead time also alters the sensitivity to variable signals , acting as a frequency filter .
power density spectra ( pds ) , in particular , are deformed to a `` wavy '' shape that depends on the magnitude of dead time and on count rates ( see figure  fig : cpdsvspds ) .
power at frequencies slightly above @xmath17 is _ quenched _ , as there is a lack of events whose separation is less than @xmath18 , while there is a relatively higher rate above @xmath18 , and therefore the power at frequencies just below @xmath17 is slightly _ amplified_. these `` waves '' have nodes at @xmath17 and multiples thereof , where the power ( in @xcite normalization ) is equal to 2 , the value that it would have without dead time , and maxima and minima in between are given by the relative contribution of the quenching and amplification . for frequencies
@xmath19 , the main effect is a general deficiency of events , and the power has a decreasing level that approaches @xmath20 @xcite . assuming that @xmath16is the same for each event and that only source events contribute to either non - paralyzable or paralyzable dead time , this distortion can be modeled precisely ( see * ? ? ?
* ; * ? ? ?
* for an exhaustive treatment ) . also , some statistical properties of the pds hold in dead time - affected data . for example , the standard deviation associated with the bin @xmath21 of the pds is always equal to @xmath22 , where @xmath23 is the power in the bin @xmath21 and @xmath24 is the number of averaged pdss ( see figure  fig : cpdsvspds ) .    in _
nustar_@xmath16is not strictly fixed at the same value for all events , but varies by a few percent depending on the number of pixels that are triggered . for this reason , the models available in literature do not correctly describe the dead time effects for this satellite : the `` wavy '' behavior of the pds shifts slightly , and to fully account for this effect and produce a white - noise subtracted pds , a very precise modeling of the dead time would be required . since at high count rates the `` waves
'' can be very prominent , any real variability feature such a qpo can easily be `` hidden '' and difficult to detect .    as an additional complication , the models described above assume that dead time is produced completely by the recorded signal . in _ nustar _ , additional dead time comes from events outside the source extraction region , from vetoed events , and from all events discarded for other reasons during the cleaning process in the pipeline ( the step from unfiltered to cleaned event files ) . in the following ,
we neglect the effect of vetoed events , since their dead time ( @xmath25s ) has characteristic frequency @xmath26khz , much higher than science events , and their total contribution to dead time is small .
we instead present a method that permits construction of a proxy of a white - noise subtracted pds , regardless of the count rate and the ratio between source and background ( or spurious ) events .
_ nustar_has two completely independent focal plane modules ( each containing four detectors ) that are read out by separate microprocessors .
it is therefore possible in principle to obtain the same information given by a pds through the so - called cross power density spectrum ( cpds ; for more details see @xcite ) : instead of considering the pdss in the two individual focal planes : @xmath27 where @xmath28 indicates the fourier transform of the lightcurve detected by the focal plane @xmath21 and @xmath29 is the frequency , one multiplies the complex conjugate of one fourier transform with the other fourier transform : @xmath30 the cpds is often used in other contexts to obtain information on the correlation between the signal in two energy bands .
it is a complex quantity : its real part is also called the _ cospectrum _ and gives a measure of the signal that is _ in phase _ between the two channels ; its imaginary part , or _ quadrature spectrum _ , gives instead a measure of the _ off - phase _ signal .
therefore , in principle , it should be possible to eliminate all variability that is not related between the two light curves , including the effects of dead time , by just considering the cospectrum ( the real part of the cpds ) . in figure
fig : cpdsvspds we show the statistical properties of the cospectrum in the case of pure poisson noise . in both the dead time affected and in the zero - dead time cases the cospectrum mean value is 0 ( in figure  fig : cpdsvspds , it has been shifted to 2 for graphical reasons ) .
this is a big advantage , as this is independent of whether the dead time is constant or not ( since the distribution of dead time is also independent between the two detectors ) , and therefore it is not necessary to conduct complicated studies of the dead time distribution in order to obtain a white noise - subtracted cospectrum , as opposed to a pds where this procedure would be needed .
another important point is that the standard deviation @xmath31 of the cospectrum is linked to the standard deviations ( @xmath32 ) of the pdss of the two light curves . in figure  1 , we show that the ratio between @xmath31 and the geometric average of the standard deviations of the two pdss is close to @xmath33 regardless of the frequency .
this is very convenient , as it allows us to assign to the cospectrum bin @xmath34 an uncertainty of @xmath35 , where @xmath36 is the geometric average of the two pdss and @xmath24 is the number of averaged pdss .
the argument ( or angle ) of the cpds is called _ phase lag _ ; if we divide it by @xmath37 we obtain the _ time lag _ @xmath38 , that is a measure of the time shift between the two channels . in our case , since we are using two lightcurves in the same energy bands , we expect any time lag to be zero .
but time lags between different energy bands have been used as indicators of how a signal produced in an emitting region is reprocessed in other regions : for example , if the disk emission is comptonized in a corona , or when the signal from a region is reflected or propagated in another region , this time delay between the different emissions can be detected through time lags ( see , e.g. , @xcite and references therein ) . we briefly discuss time lags between different energy bands in the following sections , but for the most part our analysis will be performed on two channels in the same energy band , and the only important quantity for our treatment will be the cospectrum .    amplitude ( 8% ) at different frequencies .
the units in this plot , and in all following power spectra or cospectra if not stated otherwise , are power@xmath4frequency .
grey points show the standard , dead time - affected pds . the white noise level to subtract in the dead time - affected datasets was calculated between 10 and 20hz ( minimum between two qpos ) , while in the dead time - free case we subtracted the theoretical level ( 2 in leahy normalization ) . ]
in order to provide a fully consistent treatment of dead time effects in the pds and in the cospectrum , we ran a large number of simulations . in each simulation
we produced two event series containing variability , one for each fpm , and analyzed the data with cospectrum and pds before and after applying a dead time filter . by doing so
, we studied in detail the properties of the cospectrum and compared them to those of the pds . in the following paragraphs we will explain the procedure in more detail , and demonstrate that the cospectrum can be considered a very good proxy of a white noise - subtracted pds , albeit with some corrections to account for in the measured @xmath39 .
[ [ light - curve - generation ] ] light curve generation + + + + + + + + + + + + + + + + + + + + + +    we used the procedure by @xcite , introduced to astronomers by @xcite , to simulate light curves between two times @xmath40 and @xmath41 , from a number of model pds shapes containing qpos .
the sampling frequency of the light curves was at least four times higher than the maximum frequency of the variability components included in the simulation .
we normalized the light curves in order to have the desired mean count rate and total @xmath39variability ( 7 - 10% ) . in order to be able to later calculate pdss with a given maximum timescale @xmath42 ( see `` calculation of the pdss and the cospectrum '' below ) , we simulated light curves at least 10 times longer , following the prescriptions in @xcite to avoid aliasing .
[ [ event - list - generation ] ] event list generation + + + + + + + + + + + + + + + + + + + + +    from every light curve , we generated two event lists , corresponding to the signal from the two focal planes .
each event list was produced as follows : first of all , we calculated the number @xmath43 of event times to be generated as a random sample from a poisson distribution centered on the number of total photons expected ( summing up all expected lightcurve counts ) ; then , we generated @xmath43 events by using a _ monte carlo acceptance - rejection method_. this is a classical monte carlo technique ; a more general treatment can be found in most textbooks on monte carlo methods ( e.g. * ? ? ?
* ) . in our case
, we used the following procedure : ( 1 ) for every event we simulated an event time @xmath44 uniformly distributed between @xmath40 and @xmath41 , and an associated random amplitude ( `` probability '' ) value @xmath45 between 0 and the maximum of the light curve ; ( 2 ) we rejected all @xmath44 whose associated @xmath45 was higher than the light curve at @xmath44 ; to avoid possible spurious effects given by a stepwise model light curve , we used a _
cubic spline _ interpolation to approximate the light curve between bins ; ( 3 ) we sorted the event list by @xmath44 .
to simulate the effects of background ( spurious events filtered out by the pipeline , events recorded outside the source regions , etc . )
, we also produced two background series , at constant average flux , one for each simulated source light curve .
[ [ dead - time - filtering ] ] dead time filtering + + + + + + + + + + + + + + + + + + +    for each event list , we created a corresponding dead time - affected event list by applying a simple dead time filter : for each event , we eliminated all events in the 2.5ms after it .
source and background events contributed equally to dead time .
we used different versions of the dead time filter by varying slightly the dead time between events ( of @xmath46ms ) .
but the pernicious effect of variable dead time is mainly on white noise subtraction . in our case
( figure  fig : cpdsvspds ) , the cospectrum allows us to overcome this problem as its white noise is 0 , and we verified that the other effects are not significantly different in the constant and variable dead time cases . in the following , we will treat the case with constant dead time .    [
[ calculation - of - the - pdss - and - the - cospectrum ] ] calculation of the pdss and the cospectrum + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    we divided each pair of event lists in segments of length @xmath42 , calculated the pds in each of segments and the cpds from each pair of them .
we then averaged the pdss and cpdss from all segments .
the cpds white noise level is already 0 . for the pdss ,
used only in the ideal zero - dead time case , we subtracted the theoretical poisson level ( 2 in leahy normalization ) .
we then rebinned the pdss and cpdss , either by a fixed rebin factor or by averaging a larger number of bins at high frequencies , following approximately a geometric progression .
we finally multiplied the leahy - normalized pdss by @xmath47 ( where b is the background count rate and s the mean source count rate ) in order to obtain the squared @xmath39normalization often used in literature @xcite .
the cpds was instead multiplied by the factor @xmath48 , where bars indicate the _ geometric _ averages of the count rates in each of the two event lists used to calculate it , it s as if each fourier transform contributing to it were multiplied by @xmath49 .
since the contribution to the cpds is by two separate fourier transforms , we multiplied by @xmath50 , that is equivalent to the geometric averages described . ] .    finally , we calculated the cospectrum by taking the real part of the cpds .
as described above , we assigned to each final bin of the cospectrum @xmath34 an uncertainty calculated from the geometrical average of the pdss in the two channels , divided by @xmath51 , where @xmath24 is the number of averaged spectra and @xmath52 is the number of subsequent bins averaged to obtain @xmath34 .
[ [ fitting - procedure ] ] fitting procedure + + + + + + + + + + + + + + + + +    cospectra do not need poisson noise subtraction ; for pdss , we fitted a constant to the interval outside of the frequency range containing qpos . since in the following paragraphs , in all our examples
, we will be showing power spectra obtained by averaging more than 50 pdss , we are in the gaussian regime and a fit with standard @xmath53 minimization routines is appropriate to the precision we are interested in @xcite .
then , we fitted the qpos with a lorentzian profile in xspec @xcite .
errors were calculated through a monte carlo markov chain , as those intervals where the @xmath53 of the fit with no frozen parameters increased by 1 . according to @xcite , the significance of the detection of qpos is expected to be @xmath54 where @xmath55 is the @xmath39and @xmath56 is the equivalent width of the feature ( for a lorentzian , @xmath57 ) .
the exact proportionality factor depends on the definition of significance . in our case ,
the significance of qpos was defined as the ratio between the amplitude of the lorentzian and its error .
the significance calculated in this way gives a value @xmath58 times lower than obtained by using the excess power  la @xcite ( a factor 2 is expected due to the fact that the gaussian errors are calculated over @xmath59 , and the excess power over @xmath60 ; see * ? ? ?
* ; * ? ? ?
* ) , but the trend in equation  eq : sig holds provided that the variability is dominated by poisson noise .      of a qpo at different peak frequencies , measured with the various techniques and with and without dead time .
each point represents a simulated qpo with @xmath39=10% and fwhm=2hz .
a total of 281 simulations were used for this plot .
( bottom ) significance measured with each method .
the total pds has about twice the significance of the single - module pds in the no - dead time case , as expected , due to double number of photons .
the cpds in the no - dead time case is a factor @xmath61 higher than the single pds , and lower by the same amount than the total pds .
the dead time - affected cpds , instead , has a much lower level due to the lack of photons .
the decrease of significance does not depend on the frequency of the qpo , but only on count rate ( see figure  fig : cts ) . ]    .
a total of 118 simulations were used in this plot .
the line shows equation  eq : drms .
it is not a fit , and describes remarkably well the data . in the bottom panel ,
we plot the detection significance for all the cases .
the lines , again , are not fitted , they just show the theoretical prediction from equation  eq : sig .
all dead time - free cases are in good agreement with a linear increase with count rate below @xmath62 , above which some curvature appears due to the departure from the quasi - poissonian regime .
the dead time - affected case is in good agreement with equation  eq : sig if , instead of the incident count rate ( dashed line ) , one uses the _ observed _ count rate ( solid line ; see equation  eq : drms ) .
]    on the relative contribution of the source in a given energy range to the total count rate .
the 132 simulations that compose this plot show how the @xmath39drop is stronger if the source signal dominates the total signal , since the source signal contributes more to the total dead time .
the solid black line shows equation  eq : drmstot .
it is not a fit . ]
[ [ first - look ] ] first look + + + + + + + + + +    the simulation in figure  fig : cpdsvspds shows the comparison between the pds and the cospectrum with and without dead time for pure poisson noise . from these simulations
it is clear that the white noise - subtracted pds and the cospectrum are equivalent in the case with no dead time .
it is immediately evident that the most problematic effect of dead time , the modulation of the white noise level , disappears in the cospectrum . in the following paragraphs
, we investigate the frequency and count rate dependence of these quantities in more detail .
[ [ frequency - dependence ] ] frequency dependence + + + + + + + + + + + + + + + + + + + +    the general statistical properties of the pds and the cospectrum are also very similar , both in the dead time - affected and in the zero - dead time case .
figure  fig : cpdsvspds , panel d , shows that the variance of the cospectrum and the pds maintains a constant ratio equal to 2 in both the clean and the dead time - affected datasets .
this makes it easy to calculate the variance of the cospectrum values for the subsequent analysis , by simply using the known properties of the pds where the variance is just equal to the square of the power ( in leahy normalization ) .
contrary to what might be imagined , it is possible to detect variability even at frequencies that are affected the most by dead time , i.e. those above @xmath17 .
figure  fig : qpos shows that qpos at all frequencies are detectable , albeit with some modulation of the observed @xmath39 . to measure this change of @xmath39 , we simulated @xmath0500 lightcurves using the method above , each containing a single qpo with frequencies equally distributed between 5 and 1000hz , @xmath39=10% and fwhm=2hz .
as explained above , from every light curve we obtained two event lists in order to simulate the signals from the two detectors .
we produced for every pair of event lists the cospectrum , the two pdss , and a total pds including the counts from both detectors , both in the zero - dead time and the dead time - affected cases .
we then fitted the resulting spectra with a lorentzian model in xspec .
figure  fig : freq shows the change of the @xmath39measured with the pds and the cospectrum , with and without dead time .
the measured @xmath39 in the zero - dead time pdss agrees with the zero - dead time cospectrum , whereas the dead time - affected cospectrum yields a frequency - dependent deviation from the true @xmath39 , with deviation following the same trend as the variance ( see also figure  fig : cpdsvspds ) .
the detection significance does _ not _ depend on the frequency in any case , with or without dead time .
the decrease of the significance is instead driven by the observed count rate , as we discuss shortly . in the no - dead time case ,
the significance of the single pdss is about half that of the total pds because the significance is directly proportional to the intensity of the signal ( equation  eq : sig ) , and in the total pds one uses twice the number of photons .
the zero - dead time cospectrum yields instead a significance @xmath61 lower than the total pds , and higher by the same amount than the single - module pds .
this is just an effect of the factor @xmath33 between the standard deviations of the cospectrum and the single - module pds . from figure
fig : freq , it is clear that it is advantageous to use the total pds for low count rates where dead time is negligible , and the cospectrum otherwise , but with the formulas and the simulations shown above to account for the frequency - dependent distortion of the @xmath39amplitude .    in summary , the important point that figure  fig : freq makes is that qpos are still detectable at any frequency , even those heavily affected by dead time , albeit with a change of the measured @xmath39that must be taken into account .
[ [ count - rate - dependence ] ] count rate dependence + + + + + + + + + + + + + + + + + + + + +    we now investigate how the measured @xmath39is influenced by count rate .
figure  fig : cts and  [ fig : ctsback ] show the variation with count rate in the detected @xmath39of a qpo at 30hz , fwhm=2hz and @xmath39=7.5% , in two cases : increasing total count rate , and fixed total ( source + background ) count rate with variable source count rate .
for the first case , since @xmath63 and the higher - order corrections are not needed , we use van der klis ( 1989 ; equations 3.8 and 4.8 ) to obtain @xmath64 this relation is plotted with a dashed line in the top panel of figure  fig : cts and it is in remarkably good agreement with the simulated data .    in general , one would expect the significance of detection in the pds to be proportional to the incident count rate and to the square of the @xmath39(equation  eq : sig ) .
this condition holds if the qpo can be considered a small disturbance in an otherwise poissonian process , or @xmath65 @xcite . in the bottom panel of figure
fig : cts , we fit equation  eq : sig below 600@xmath10 , with a multiplicative constant due to the slightly different definition of 1-@xmath66 error that we use ( @xmath67 instead of the @xcite definition ) .
the best - fit multiplicative constant , @xmath68 , turns out to be consistent with the factor 2 expected from the fact that we are using gaussian fitting instead of excess power ( see section  sec : simutec ) .
the departure from the linear condition above @xmath69cts is evident .
indeed , it is expected that at count rates above @xmath70 , the significance starts departing from the linear trend .
the total pds is visibly more affected because its count rate is double than that of the single - module pds .
the significance of detection with the cospectrum is @xmath61 lower than that of the total pds , while it is @xmath61 higher than that of the pds from a single module .
the dead time - affected cospectrum , instead , has a large deviation from the linear trend .
this is just an effect of the diminishing count rate due to dead time .
in fact , what is plotted is the _ incident _ count rate .
if one converts it to the _ detected _ count rate , the linear relation between count rate and significance still holds ( solid line ) .
the second case ( figure  fig : ctsback ) clearly shows a linear decrease of the measured @xmath39as the source gains counts with respect to the background .
again , by using van der klis ( 1989 ; equations 3.8 and 4.8 ) , but this time putting the _ total _ count rate ( @xmath71 , where @xmath72 is the non - source count rate ) in the relation between incident and observed count rates , one obtains @xmath73\ ] ] this means that the measure of @xmath39we obtain in our data will be generally affected more if the source signal dominates the background , as is the case in most _
nustar_observations of bright sources . in the examples that we present in the next sections ,
we make use of monte carlo simulations similar to the ones above to estimate the change of @xmath39at the count rate of the sources we observe .
we now demonstrate the methods outlined above by analyzing the power density spectra for a number of bright galactic black hole binaries .
a summary of the observations used in this work is listed in table  tab : obs .
we preprocessed all _
nustar_observations with the nustardas pipeline included in heasoft 6.15 .
we produced cleaned event files and calculated good time intervals ( gtis ) .
we processed both the unfiltered and the clean event files ( files tagged with `` _ uf '' and `` _ cl '' respectively in the `` event_cl '' directory of _
nustar_data ) . by analyzing the unfiltered and cleaned event files , it is possible to estimate the ratio of `` good '' events to the total of dead time - producing events .
we calculated the ratio @xmath74 of `` good '' events over intervals of 100s .
we found that in all observations @xmath74 was almost constantly @xmath75 , with some drops due to increased solar activity .
we selected the intervals in which this ratio did nt go below a certain threshold ( ranging from 0.9 to 0.92 in different files ) , in order to eliminate possible spurious effects on the observed variability .
in the subsequent analysis , we will always consider an additional 10% `` background '' event rate due to unrecorded and spurious events .
we selected events from a region around each source depending on how broad the point - spread distribution was ( see the next subsections for details for each source ) .
we used the barycorr ftool to calculate event times at the solar system barycenter using the de200 ephemeris and source positions from simbad .
additional variability is introduced into the cross - power spectrum in some cases due to mis - matched gti between fpma and fpmb .
this mis - match occurs often when the satellite traverses the south atlantic anomaly ( saa ) because the automatic algorithm that turns off the focal plane modules to conserve telemetry kicks in at different thresholds for fpma and fpmb .
the parameters that rule this behavior are finely tuned in order to maximize the number of photons available for spectral analysis while avoiding most of the contamination @xcite , but for timing analysis this produces unwanted effects that have to be accounted for .
a signature of this variability is a strong time lag between the two detectors in the same energy band ( therefore , unrelated to the physics of the source ) , that disappears when discarding an additional interval at the gti borders or by using a more aggressive filtering during the pipeline process .
we decided to use a `` safe '' additional cut of 300s at the borders of all gtis during timing analysis , and verified that it solved the problem in all affected obsids .
from each barycentered event list ( section  sec : red ) we obtained light curves sampled at @xmath76s ( @xmath77s ) .
we extracted separate light curves for each detector , and calculated the common good time intervals ( gtis ) .
each light curve was then divided in segments ( the length was chosen in order to optimize the use of gtis ) with the start and stop times locked between the two detectors , and a fourier transform was calculated from each of them .
we then calculated the cpdss from each segment by multiplying the fourier transforms as explained in section  sec : deadtime , and obtained a single cospectrum from the real part of the average of all cpdss .    when needed , we plot the error on the best - fit model as a hatched region around the model line . to calculate these errors , we ran a monte carlo markov chain on xspec to simulate 10000 models around the best fit .
we chose all the models having @xmath78 , and we calculated at each energy the maximum deviation from the best fit coming from any of these models .
cygx-1is a bright persistent bhb , discovered during one of the earliest rocket flights at the dawn of x - ray astronomy @xcite . as such
, it is among the best studied x - ray sources in the galaxy .
its spectral evolution is typical for a bh , with the succession of low / hard , high / soft and intermediate states with correspondingly different timing patterns ( see * ? ? ?
* for a broad view on bhbs ) .
variability in all spectral states has been observed .
it shows two dominant lorentzian components in the low / hard state .
when the source softens their characteristic frequencies increase and an additional cutoff power law - like component appears . in the soft state
the latter dominates the power spectrum ( * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* and references therein ) .
time lags have been observed in all spectral states @xcite ; the variability of single spectral components was also investigated , showing that the reflection component shows variability at lower frequencies than the bulk of the emission in the low / hard state @xcite while the @xmath39level and shape of the power spectrum is comparable at all energies in the soft state @xcite .
a very thorough investigation of the relation between variability in different energy ranges and spectral states can be found in @xcite . for past studies of the high - energy variability of the source
see @xcite .    _
nustar_observed cygx-1 in its soft state ( for a full spectral analysis description , see * ? ? ?
the spectral analysis shows a strong reflection component , that might produce detectable phase lags .
phase lags have been observed with _
rxte_in all spectral states ( @xcite ; see @xcite for a rigorous treatment ) .
the pds in this state is often well described by a cutoff power law - like shape almost `` flat '' between @xmath79 and @xmath80hz ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
this is the case in our observations , as shown in figure  fig : cyg_pdsvse .
there is a clear change in @xmath39variability between the two observations , and a power law - like red noise component appears at the lowest fourier frequencies in only one observation , but for the rest the cospectrum is very well described by a cutoff power law .
as shown in figure  fig : cts , a drop of @xmath39is expected at higher count rates , and the second observation has indeed a higher count rate than the first .
but in figure  fig : cyg : rms we show that the expected @xmath39drop due to dead time is far lower than observed .
this means that the rms change is largely due to an intrinsic lower variability of the source in the second observation .
more interesting is the cospectrum in different energy ranges , shown in figure  fig : cyg_pdsvse .
there is a very strong change in the shape and normalization of the cospectrum with energy . in general , the @xmath39increases with energy , and it is about an order of magnitude higher above 10 kev than below 5 kev .
this is a clear indication that the soft component of the spectrum ( coming from the disk ) is varying less than the reprocessed component ( reflected and comptonized ) .
again , in order to evaluate the @xmath39change with energy , it is necessary to apply a correction for dead time .
in fact , channels with more counts will have a larger drop of @xmath39than channels with a few counts ( compare to figure  fig : ctsback ) .
following the procedure in section  sec : simu , we simulated the expected relative drop in @xmath39 in all bands . in figure  fig
: cyg : rms we show the measured @xmath39and the curves showing the expected drop due to count rate .
these curves were calculated with the method depicted in figure  fig : cts , assuming a total count rate @xmath81 times the source count rate , as the total background plus the rejected events account for about 10% of the contributions
to dead time . as can be seen
, the @xmath39change between bands is mostly intrinsic , with a drop of only @xmath82% expected for the channels with the highest count rates .
our cospectrum results are consistent with the pds results from other works ( e.g. * ? ? ?
although we do not find any significant time lags between any of the energy bands , this is mostly due to the fact that our observation is not long enough .
we estimate that @xmath83ks should be sufficient to detect time lags similar to those observed in the past in cygx-1 , in this spectral state ( e.g. * ? ? ?
* ) .
[ fig : cyg : rms ]      gx339@xmath14is another well known bhb . discovered in 1973 by _ oso-7 _
@xcite , it has been observed frequently since then .
it is a transient system , with long outbursts and periods of quiescence lasting up to several years ( * ? ? ?
* and references therein ) .
monitored during its outbursts by _ rxte _ , it showed one of the best examples of q - track pattern in the hardness - intensity diagram ( @xcite ; for the q - track pattern see @xcite and @xcite for a comparison with other sources ) .
like many accreting bhs , its low / hard state pds is dominated by two or more lorentzian components , whose frequencies are positively correlated to count rate . as the source reaches the intermediate and soft state ,
the broad lorentzians move further to higher frequencies and various kinds of qpos appear ( see , e.g. , * ? ? ?
* ; * ? ? ?
* ) , including high frequency qpos .
thanks to this wealth of data , it has often been used , like cygx-1 , to test new timing approaches . in its very high state ,
the signal at low energy ranges is generally coherent , while coherence is lost between low and high energies @xcite .
@xcite measured coherence and lags in the hard state of the source , showing that the two lorentzians are likely produced by independent , but internally coherent , processes .
@xcite showed that hard time lags in this source are more likely to be produced by propagation of modulations of the accretion flow in the disk rather than from reflection .    _
nustar_observed gx339@xmath14four times during the rise of the last outburst . during this outburst
, the source remained in a canonical low / hard state without reaching the threshold luminosity that would have marked the transition to the intermediate and soft states .
the spectral analysis will be published in a separate publication ( frst et al 2014 , _ in prep . _ ) , but in this paper we will concentrate on the timing properties of the source , analyzed with the methods described above .    the whole band and energy - dependent cospectra in the four different observations are shown in figure  fig : gx_4obs . the approximate flux of the source increased almost linearly from @xmath82 to @xmath84 in each detector between the first and the last observation . the general shape and the variation of the cospectrum is similar to what has been observed in the past for this source in the hard state ( e.g. * ? ? ?
* ; * ? ? ?
* ) and for other bhbs @xcite , with two zero - centered lorentzian components dominating the shape ( l@xmath85 and l@xmath86 in @xcite ) , plus occasionally a third lorentzian needed to account for an excess between them .
the cospectrum evolves with l@xmath85 consistently increasing its characteristic frequency in subsequent observations ( as luminosity increases ) , while l@xmath86 does not change as much .
this is consistent with what was previously reported by @xcite using the pds .
the energy dependence of the cospectrum is very weak in the first observation , and tends to increase in the subsequent observations , while the total @xmath39consistently decreases as the count rate increases .
in the second , third and fourth observation the high energy cospectrum has consistently a lower @xmath39than the low energy one .
again , these cospectra results are consistent with previous pds results ( e.g. * ? ? ?
* ) .    from figure
fig : cts , and the discussion about the much brighter cygx-1 , it is clear that at these frequencies and for count rates changing from 10 to 40@xmath10we do not expect a such a strong change in @xmath39 , which must therefore be intrinsic to the source .    a more detailed discussion will be presented in frst et al . , _ in prep._.
this paper presents a method to exploit _ nustar _ s two independent detectors to work around the spurious effects resulting from dead time , enabling standard aperiodic timing analysis of x - ray binaries to be applied . as extensively studied in literature ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) , dead time produces a spurious correlation between event times that strongly modifies the shape of the pds . for frequencies @xmath19 ,
the only visible effect in pdss is a slight decrease of the white noise level . in timing - oriented x - ray missions ,
dead time is very short ( _ rxte_/pca : 10@xmath87s , @xcite ; _ astrosat_/laxpc : 1035@xmath87s , @xcite ; _ loft _ :
@xmath88ns , @xcite ) and the `` interesting '' range of frequencies for accretion ( the highest being the dynamical timescales around neutron stars @xmath89hz ) is well below @xmath17 . in _
@xmath16is around 2.5ms , right in the middle of the frequency range commonly analyzed in bh and ns binaries .
moreover , in order to obtain better precision in the measure of photon energies , the dead time is variable and depends on the event grades . due to this , standard techniques for the modeling of the pds shape ( e.g. * ? ? ?
* ; * ? ? ?
* ) are not applicable .
but _ nustar_has two independent focal plane modules that are read out by separate processors , and the dead time is therefore independent . if , instead of taking the pds of each detector , we take the cospectrum ( see section  sec : cpds ) of the signals in the two detectors , one is able to cancel out most correlations produced by dead time and obtain a close proxy of a white noise - subtracted pds .
we described in sections  [ sec : cpds ] and  [ sec : simu ] how to calculate the cospectrum and estimate the remaining effects of dead time on the signal , namely the modulation of the variance of data bins , that can not be described using the standard @xcite prescriptions , and the equivalent change of the observed @xmath39at different frequencies and for different count rates .
this estimate requires monte carlo simulations tailored to the source that is being analyzed , as described in section  sec : simu .
we also applied this technique to the analysis of _ nustar_data of two bhbs , cygx-1and gx339@xmath14 .
cygx-1was observed in the soft state .
we report a strong change of @xmath39with energy , with a significant rise of the variability at higher energies .
these cospectrum results are consistent with earlier pds results , and extend the energy coverage with respect to earlier results from _ rxte _ ( * ? ? ?
* e.g. ) . comparing these results with the spectral modeling presented by @xcite ,
it is apparent that the low energies are dominated by the disk emission , while at higher energies it s the reprocessed component ( comptonized and reflected ) that dominates .
this is consistent with what is found in many bhbs , and in particular to what has already been reported in this source ( the  stable disk and unstable corona  picture , e.g. * ? ? ?
* ) .    in gx339@xmath14 , we observed the rising part of the outburst , in a standard low / hard state .
the cospectrum shows that , as the outburst goes on , the @xmath39generally decreases and the characteristic frequency of the first lorentzian component used to described the cospectrum consistently rises . also ,
this is consistent with what has been observed in the low / hard state of this and other sources using the pds @xcite .
we do not find evidence for time lags in the sources analyzed in this paper , due to the short observing time available , but the cross power density spectrum can be used to calculate time lags without major impact from dead time , if data are adequately filtered to avoid intervals with high background activity and in particular at the start and end of each good time interval .
the techniques presented in this paper can be used to perform most of the standard timing and spectral timing analysis used in literature , if the change of @xmath39with frequency and count rate is properly taken into account ( section  sec : simu ) . for slow timing ( @xmath90 ) , the dead time correction included in the official pipeline is still recommended , because it effectively corrects for all sources of dead time , including shielded events . for low count rates ( @xmath91 ) the standard pds of the sum of the two lightcurves ( what we called the `` total pds '' above ) is still advantageous , because the effects of dead time are negligible and the sensitivity of the total pds , expressed as significance of detection of a standard variability feature , is @xmath61 higher than the cospectrum .    as an additional note , this technique can be used in any satellite that has completely independent detectors , and of course by using multiple satellites .
_ rxte _ , with its independent pcus , would be a good candidate if not for the presence of so - called very large events ( vle ) , particle events depositing more than 100kev and occurring many times per second ( e.g. * ? ? ?
these events often involve more than one pcu , invalidating the assumption of independence .
the large area xenon propotional counter ( laxpc ) onboard the upcoming _ astrosat _ mission @xcite will provide three independent detectors .
provided that the fraction of vle - equivalent events is negligible , this technique will be applicable to laxpc data .
mb and db wish to acknowledge the support from the centre national dtudes spatiales ( cnes ) and the centre national de la recherche scientifique ( cnrs ) .
cs acknowledges funding by the german bmwi under dlr grant numbers 50  qr  0801 , 50  qr  0903 , and 50  oo  1111 .
thanks stfc for support ( grant reference st / j003697/1 ) .
this work was supported under nasa contract no .
nng08fd60c , and made use of data from the _ nustar _ mission , a project led by the california institute of technology , managed by the jet propulsion laboratory , and funded by the national aeronautics and space administration .
we thank the _ nustar _ operations , software and calibration teams for support with the execution and analysis of these observations .
this research has made use of the _ nustar _ data analysis software ( nustardas ) jointly developed by the asi science data center ( asdc , italy ) and the california institute of technology ( usa ) .
the timing analysis was executed with matteo bachetti s libraries and tools in python for _ nustar_timing ( malt pynt ) .
this code is available upon request .
most of the plots were produced with the veusz software by jeremy sanders .
the authors wish to thank chris done , denis leahy and tomaso belloni for very insightful discussions .
76 natexlab#1#1    agrawal , p.  c. 2006 ,
asr , 38 , 2989    arnaud , k.  a. 1996 , adass , 101 , 17    artigue , r. , barret , d. , lamb , f.  k. , lo , k.  h. , & miller , m.  c. 2013 , arxiv , 1303.0248 , 1303.0248    axelsson , m. , borgonovo , l. , & larsson , s. 2005 , a&a , 438 , 999    bachetti , m. et  al .
2013 , apj , 778 , 163    barret , d. , & vaughan , s. 2012 , apj , 746 , 131    belloni , t. , & hasinger , g. 1990 , a&a , 230 , 103    belloni , t. , homan , j. , casella , p. , van  der klis , m. , nespoli , e. , lewin , w. h.  g. , miller , j.  m. , & mendez , m. 2005 , a&a , 440 , 207    bendat , j.  s. , & piersol , a.  g. 2011 , random data : analysis and measurement procedures ( wiley )    bck , m. et  al .
2011 , a&a , 533 , 8    boutelier , m. 2009 , phd thesis , phd thesis , cesr    boutelier , m. , barret , d. , & miller , m.  c. 2009 , mnras , 399 , 1901    bowyer , s. , byram , e.  t. , chubb , t.  a. , & friedman , h. 1965 , science , 147 , 394    cabanac , c. , roques , j .-
, & jourdain , e. 2011 , apj , 739 , 58    casella , p. , belloni , t. , & stella , l. 2005 , apj , 629 , 403    cassatella , p. , uttley , p. , wilms , j. , & poutanen , j. 2012 , mnras , 422 , 2407    castro - tirado , a.  j. , brandt , s. , & lund , n. 1992 , iau circ . , 5590 , 2    churazov , e. , gilfanov , m. , & revnivtsev , m. 2001 , mnras , 321 , 759    cui , w. , zhang , s.  n. , jahoda , k. , focke , w. , swank , j.  h. , heindl , w.  a. , & rothschild , r.  e. 1997 , the transparent universe , 382 , 209    davies , r.  b. , & harte , d.  s. 1987 , biometrika , 74 , 95    dunn , r. j.  h. , fender , r.  p. , krding , e.  g. , belloni , t. , & cabanac , c. 2010 , mnras , 403 , 61    fabian , a.  c. et  al .
2009 , nat . ,
459 , 540    feroci , m. et  al .
2012 , experimental astronomy , 34 , 415    foster , r.  s. , waltman , e.  b. , tavani , m. , harmon , b.  a. , zhang , s.  n. , paciesas , w.  s. , & ghigo , f.  d. 1996 , astrophysical journal letters v.467 , 467 , l81    frst , f. et  al .
2013 , apj , 779 , 69    . 2014 , apj , 780 , 133    gandhi , p.
2010 , mnras , 407 , 2166    gentle , j.  e. 2003 , random number generation and monte carlo methods ( springer )    gilfanov , m. , churazov , e. , & revnivtsev , m. 2000 , mnras , 316 , 923    gleissner , t. , wilms , j. , pottschmidt , k. , uttley , p. , nowak , m.  a. , & staubert , r. 2004 , a&a , 414 , 1091    grinberg , v. et  al .
2014 , a&a , 565 , 1    harrison , f.  a. et  al .
2013 , apj , 770 , 103    jahoda , k. , markwardt , c.  b. , radeva , y. , rots , a.  h. , stark , m.  j. , swank , j.  h. , strohmayer , t.  e. , & zhang , w. 2006 , apj supp .
ser . , 163 , 401    jin , c. , done , c. , middleton , m. , & ward , m. 2013 , mnras , 436 , 3173    jones , t.  a. , levine , a.  m. , morgan , e.  h. , & rappaport , s. 2008 , apj , 677 , 1241    kaspi , v.  m. et  al .
2014 , apj , 786 , 84    krding , e. , & falcke , h. 2004 , a&a , 414 , 795    leahy , d.  a. , darbro , w. , elsner , r.  f. , weisskopf , m.  c. , kahn , s. , sutherland , p.  g. , & grindlay , j.  e. 1983 , apj , 266 , 160    lewin , w. h.  g. , van paradijs , j. , & van  der klis , m. 1988 , ssrv , 46 , 273    maccarone , t.  j. , & coppi , p.  s. 2003 , mnras , 338 , 189    marinucci , a. et  al .
2014 , mnras , 440 , 2347    markert , t.  h. , canizares , c.  r. , clark , g.  w. , lewin , w. h.  g. , schnopper , h.  w. , & sprott , g.  f. 1973 , apj , 184 , l67    mcclintock , j.  e. , & remillard , r.  a. 2006 , in compact stellar x - ray sources , ed .
j.  e. mcclintock & r.  a. remillard , 157213    miller , j.  m. et  al .
2013 , apjl , 775 , l45    miyamoto , s. , kimura , k. , kitamoto , s. , dotani , t. , & ebisawa , k. 1991 , apj , 383 , 784    miyamoto , s. , & kitamoto , s. 1989 , nat . , 342 , 773    miyasaka , h. et  al .
2013 , apj , 775 , 65    morgan , e.  h. , remillard , r.  a. , & greiner , j. 1997 , apj , 482 , 993    mori , k. et  al .
2013 , apjl , 770 , l23    natalucci , l. et  al .
2014 , apj , 780 , 63    nowak , m.  a. , & vaughan , b.  a. 1996 , mnras , 280 , 227    nowak , m.  a. , vaughan , b.  a. , wilms , j. , dove , j.  b. , & begelman , m.  c. 1999 , apj , 510 , 874    nowak , m.  a. , wilms , j. , & dove , j.  b. 1999 , apj , 517 , 355    pahari , m. , neilsen , j. , yadav , j.  s. , misra , r. , & uttley , p. 2013
, apj , 778 , 136    papadakis , i.  e. , nandra , k. , & kazanas , d. 2001 , apj , 554 , l133    paul , b. , & team , l. 2009 , in astrophysics with all - sky x - ray observations , 362    pottschmidt , k. , wilms , j. , nowak , m.  a. , heindl , w.  a. , smith , d.  m. , & staubert , r. 2000 , a&a , 357 , l17    reig , p. , belloni , t. , van  der klis , m. , mendez , m. , kylafis , n.  d. , & ford , e.  c. 2000 , apj , 541 , 883    remillard , r.  a. , muno , m.  p. , mcclintock , j.  e. , & orosz , j.  a. 2003 , aas , 7 , 648    revnivtsev , m. , gilfanov , m. , & churazov , e. 1999 , a&a , 347 , l23    . 2000 , a&a , 363 , 1013    risaliti , g. et  al .
2013 , nat . , 494 , 449    strohmayer , t.  e. 2001 , apj , 554 , l169    timmer , j. , & koenig , m. 1995 , a&a , 300 , 707    tomsick , j.  a. , & kaaret , p. 2001 , apj , 548 , 401    tomsick , j.  a. et  al .
2014 , apj , 780 , 78    uttley , p. , cackett , e.  m. , fabian , a.  c. , kara , e. , & wilkins , d.  r. 2014 , arxiv , 6575 , 1405.6575    uttley , p. , wilkinson , t. , cassatella , p. , wilms , j. , pottschmidt , k. , hanke , m. , & bck , m. 2011 , mnras let . , 414 , l60    van  der klis , m. 1989 , in timing neutron stars : proceedings of the nato advanced study institute on timing neutron stars held april 4 - 15 , 27    vaughan , b.  a. , & nowak , m.  a. 1997 , apjl , 474 , l43    vaughan , s. 2013 , arxiv , 1309.6435v1    vikhlinin , a. , churazov , e. , & gilfanov , m. 1994 , a&a , 287 , 73    walton , d.  j. et  al .
2014 , arxiv , 5620 , 1404.5620    weisskopf , m.  c. 1985 , in talk presented at the workshop _ time variability in x - ray and gamma - ray sources _ , taos , nm , usa    zhang , w. , jahoda , k. , swank , j.  h. , morgan , e.  h. , & giles , a.  b. 1995 , apj , 449 , 930    zoghbi , a. et  al .
2014 , apj , 789 , 56
grs1915@xmath2105is a very bright persistent microquasar , discovered in 1992 by _ granat _ @xcite which shows accretion states that depart from the usual behavior of galactic bhs , being associated with very high - eddington fractions .
the variability of this source has been extensively studied in the past , in particular for the presence of very strong low frequency ( @xmath120.1hz ) , medium - low frequency ( 0.510hz ) , medium - high frequency ( @xmath040 and 67hz ) and high frequency ( 113 and 165hz ) qpos ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) . in our observation ,
the medium - low qpo is the only one observed to very high significance .
its behavior has been studied in great detail ( e.g. * ? ? ? * ) and its properties are known up to very high energies thanks to _
rxte_/hexte observations @xcite .
this source is discussed as a reference for timing analysis on sources observed very early in the mission .
the first version of the flight software , in fact , introduced a large number of spurious modulations in the signal that we will discuss and show a way to work around .    during the observatory commissioning phase prior to commencement of science observations , _ nustar_observed the galactic bhb grs1915@xmath2105 in its _ plateau _ state @xcite .
the spectral analysis of that observation was presented by @xcite .
the version of the flight software that was installed at that time on _ nustar_executed housekeeping operations at regular intervals ( 1 or 8s ) , producing periodic dead time .
therefore , the raw power spectrum and the cpds of this observation are plagued by spurious features at 1hz and 0.125hz as well as their harmonics and sub - harmonics ( see figure  fig : grscpds ) . in this section
we show a technique to clean up the cospectrum and perform a basic analysis of the strong qpo that this source .    to `` clean '' the cospectrum from these spurious features
, we first produced a cospectrum on a time scale @xmath92 of 128s .
the maintenance operations were very regular , and so in this cospectrum the corresponding features were very sharp and limited to a single bin each .
we replaced the value of each of these bins with a random value calculated from the mean of the two closest bins and their standard deviation , as shown in figure  fig : grscpds .
the possible spurious effects of this imperfect statistical treatment of the simulated bin value are negligible due to the small frequency bin ( simulated data are only one bin every @xmath93hz ) and the rebinning applied in the following steps of the analysis . the real and imaginary part were treated independently .    from this cleaned cpds we calculated the cospectrum .
an example is shown in red in figure  fig : grscpds .
the cospectrum shows a very strong and variable qpo at @xmath94hz , with @xmath39=10(1)% ( from the best - fit lorentzian amplitude , and corrected for the expected @xmath95% drop at @xmath96 , equation  eq : drms ) .
this feature has previously been observed in this spectral state ( e.g. * ? ? ?
* ) . in figure
fig : grscpds we show that with these techniques it is possible in principle to also apply techniques of qpo tracking ( a dynamical cospectrum ) given a sufficient @xmath39 .
a hint of the harmonic of this qpo is also observable , both in the single and in the dynamical cospectra .
no other significant qpos are detected .
this source is known to show high - frequency qpos in _
rxte_observations @xcite .
these features are not observed in our dataset .
we calculate a rough upper limit of @xmath97% for additional qpos between 50 and 100hz .
this upper limit is higher than the observed @xmath39of qpos in this frequency range observed in the past ( e.g. @xmath94% in the full _ rxte_band in @xcite ) , and our non - detection might be largely due to the low statistics of this dataset .