Next: , Previous: , Up: Continuous-Wave Functions   [Contents][Index]

1.2.4 SensitivityDepthStackSlide()

SensitivityDepthStackSlide() estimates the sensitivity of a continuous gravitational-wave search. Given the search parameters: 

octave> Nseg = 90;
octave> Tseg = 60 * 3600;
octave> Tdata = 12080 * 1800;
octave> avg2Fth = 6.109;
octave> sum2Fth = Nseg * avg2Fth;
octave> pFA = falseAlarm_chi2(sum2Fth, 4 * Nseg);
octave> misHistSC = createDeltaHist (0.7);

and the desired false dismissal probability, the search sensitivity is estimated as a sensitivity depth: 

octave> pDET = [0.95; 0.9; 0.85];
octave> Depths = SensitivityDepthStackSlide("Nseg", Nseg, "Tdata", Tdata, "misHist", misHistSC, "pFD", 1-pDET, "pFA", pFA, "detectors", "H1,L1")
Depths =

   38.683
   41.919
   44.495

The inverse function DetectionProbabilityStackSlide() instead computes the detection probabilities from the desired sensitivity depths: 

octave> detProbs = DetectionProbabilityStackSlide("Nseg", Nseg, "Tdata", Tdata, "misHist", misHistSC, "pFA", pFA, "detectors", "H1,L1", "Depth", Depths)
detProbs =

   0.95000
   0.90000
   0.85000