脉冲星搜索软件 PRESTO 的大致流程-附代码

PRESTO is a large suite of pulsar search and analysis software developed primarily by Scott Ransom mostly from scratch, and released under the GPL (v2). It was primarily designed to efficiently search for binary millisecond pulsars from long observations of globular clusters (although it has since been used in several surveys with short integrations and to process a lot of X-ray data as well). It is written primarily in ANSI C, with many of the recent routines in Python. According to Steve Eikenberry, PRESTO stands for: PulsaR Exploration and Search TOolkit!

 

1，教程：

https://github.com/scottransom/presto

2，代码：

https://github.com/374365283/Pulsar_Search

source /root/presto.sh
export PRESTO=/usr/local/astrosoft/presto
export INPUT=GBT_Lband_PSR.fil
#readfile $INPUT 
time1=`date +%s`
rfifind -time 2.0 -o Lband $INPUT  
time2=`date +%s`
rfifind -time 1.0 -o Lband $INPUT 
time3=`date +%s`
rfifind -nocompute -time 1.0 -freqsig 6.0 -mask Lband_rfifind.mask -o Lband $INPUT  
time4=`date +%s`
prepdata -nobary -o Lband_topo_DM0.00 -dm 0.0 -mask Lband_rfifind.mask -numout 530000 $INPUT  
time5=`date +%s`
#exploredat Lband_topo_DM0.00.dat 
realfft Lband_topo_DM0.00.dat 
time6=`date +%s`
#explorefft Lband_topo_DM0.00.fft 
accelsearch -numharm 4 -zmax 0 Lband_topo_DM0.00.dat 
time7=`date +%s`
#less -S Lband_topo_DM0.00_ACCEL_0
#jed Lband.birds
cp Lband_rfifind.inf Lband.inf
time8=`date +%s`
makezaplist.py Lband.birds
time9=`date +%s`
#explorefft Lband_topo_DM0.00.fft
echo|prepfold -p 1.0 $INPUT
time10=`date +%s`
echo|DDplan.py -d 500.0 -n 96 -b 96 -t 0.000072 -f 1400.0 -s 32 -r 0.5
time11=`date +%s`
prepsubband -nsub 32 -lodm 0.0 -dmstep 2.0 -numdms 24 -numout 132500 -downsamp 4 -mask Lband_rfifind.mask -o Lband $INPUT 
time12=`date +%s`
cp $PRESTO/tests/dedisp.py .
time13=`date +%s`
python dedisp.py
time14=`date +%s`
mkdir subbands
mv *.sub* subbands/
rm -rf Lband*topo*
time15=`date +%s`
ls *.dat | xargs -n 1 realfft
time16=`date +%s`
ls *.fft | xargs -n 1 zapbirds -zap -zapfile Lband.zaplist -baryv -5.69726e-05
time17=`date +%s`
ls *dat | xargs -n 1 accelsearch -zmax 0
time18=`date +%s`
cp $PRESTO/examplescripts/ACCEL_sift.py .
time19=`date +%s`
python ACCEL_sift.py > cands.txt
time20=`date +%s`
echo|prepfold -accelcand 2 -accelfile Lband_DM62.00_ACCEL_0.cand Lband_DM62.00.dat
time21=`date +%s`
#ls subbands/
echo|prepfold -accelcand 2 -accelfile Lband_DM62.00_ACCEL_0.cand -dm 62 subbands/Lband_DM72.00.sub??
time22=`date +%s`
echo|prepfold -n 64 -nsub 96 -p 0.004621638 -dm 62.0 $INPUT
time23=`date +%s`
single_pulse_search.py *dat
time24=`date +%s`
duration=`echo "$time24-$time1" | bc -l`
echo "total time = $duration sec"
duration=`echo "$time2-$time1" | bc -l`
echo "rfifind1 = $duration sec"
duration=`echo "$time3-$time2" | bc -l`
echo "rfifind2 = $duration sec"
duration=`echo "$time4-$time3" | bc -l`
echo "rfifind3 = $duration sec"
duration=`echo "$time5-$time4" | bc -l`
echo "preparedata = $duration sec"
duration=`echo "$time6-$time5" | bc -l`
echo "realfft = $duration sec"
duration=`echo "$time7-$time6" | bc -l`
echo "accelsearch = $duration sec"
duration=`echo "$time9-$time8" | bc -l`
echo "makezaplist = $duration sec"
duration=`echo "$time10-$time9" | bc -l`
echo "prepfold = $duration sec"
duration=`echo "$time11-$time10" | bc -l`
echo "DDplan = $duration sec"
duration=`echo "$time12-$time11" | bc -l`
echo "prepsubband = $duration sec"
duration=`echo "$time14-$time13" | bc -l`
echo "dedisp = $duration sec"
duration=`echo "$time16-$time15" | bc -l`
echo "realfft = $duration sec"
duration=`echo "$time17-$time16" | bc -l`
echo "zapbirds = $duration sec"
duration=`echo "$time18-$time17" | bc -l`
echo "accelsearch = $duration sec"
duration=`echo "$time20-$time19" | bc -l`
echo "ACCEL_sift = $duration sec"
duration=`echo "$time21-$time20" | bc -l`
echo "prepfold = $duration sec"
duration=`echo "$time22-$time21" | bc -l`
echo "prepfold = $duration sec"
duration=`echo "$time23-$time22" | bc -l`
echo "prepfold = $duration sec"
duration=`echo "$time24-$time23" | bc -l`
echo "single_pulse_search = $duration sec"