53
53
< li class ="toctree-l2 current "> < a class ="current reference internal " href ="# "> Crab Analysis</ a > < ul >
54
54
< li class ="toctree-l3 "> < a class ="reference internal " href ="#Search-the-best-frequency "> Search the best frequency</ a > </ li >
55
55
< li class ="toctree-l3 "> < a class ="reference internal " href ="#Bayrcentric-correction "> Bayrcentric correction</ a > </ li >
56
+ < li class ="toctree-l3 "> < a class ="reference internal " href ="#Play-with-the-profile "> Play with the profile</ a > </ li >
56
57
</ ul >
57
58
</ li >
58
59
< li class ="toctree-l2 "> < a class ="reference internal " href ="BinaryAnalysis.html "> Binary System Analysis</ a > </ li >
357
358
</ style >
358
359
< div class ="section " id ="Pulsar-Timing-Analysis ">
359
360
< h1 > Pulsar Timing Analysis< a class ="headerlink " href ="#Pulsar-Timing-Analysis " title ="Permalink to this headline "> </ a > </ h1 >
361
+ < p > Welcome to our step-by-step Jupyter notebook tutorial on pulsar timing analysis using the TAT-Pulsar Python package.</ p >
362
+ < p > Throughout this guide, we’ll provide you with hands-on examples of how to use the key features of TAT-Pulsar.</ p >
360
363
< div class ="section " id ="Search-the-best-frequency ">
361
364
< h2 > Search the best frequency< a class ="headerlink " href ="#Search-the-best-frequency " title ="Permalink to this headline "> </ a > </ h2 >
362
365
< div class ="nbinput nblast docutils container ">
@@ -368,7 +371,7 @@ <h2>Search the best frequency<a class="headerlink" href="#Search-the-best-freque
368
371
< span class ="kn "> import</ span > < span class ="nn "> wget</ span > < span class ="o "> ,</ span > < span class ="nn "> os</ span >
369
372
< span class ="kn "> from</ span > < span class ="nn "> astropy.io</ span > < span class ="kn "> import</ span > < span class ="n "> fits</ span >
370
373
371
- < span class ="n "> mpl</ span > < span class ="o "> .</ span > < span class ="n "> rcParams</ span > < span class ="p "> [</ span > < span class ="s1 "> 'figure.dpi'</ span > < span class ="p "> ]</ span > < span class ="o "> =</ span > < span class ="mi "> 150 </ span >
374
+ < span class ="n "> mpl</ span > < span class ="o "> .</ span > < span class ="n "> rcParams</ span > < span class ="p "> [</ span > < span class ="s1 "> 'figure.dpi'</ span > < span class ="p "> ]</ span > < span class ="o "> =</ span > < span class ="mi "> 250 </ span >
372
375
< span class ="n "> plt</ span > < span class ="o "> .</ span > < span class ="n "> style</ span > < span class ="o "> .</ span > < span class ="n "> use</ span > < span class ="p "> ([</ span > < span class ="s1 "> 'science'</ span > < span class ="p "> ,</ span > < span class ="s1 "> 'nature'</ span > < span class ="p "> ,</ span > < span class ="s1 "> 'no-latex'</ span > < span class ="p "> ])</ span >
373
376
</ pre > </ div >
374
377
</ div >
@@ -444,24 +447,24 @@ <h2>Bayrcentric correction<a class="headerlink" href="#Bayrcentric-correction" t
444
447
</ div >
445
448
< div class ="output_area docutils container ">
446
449
< div class ="highlight "> < pre >
447
- PSR_B 0531+21
448
- RA_hh 5
449
- RA_mm 34
450
- RA_ss 31.972
451
- DEC_hh 22
452
- DEC_mm 0
453
- DEC_ss 52.07
454
- MJD1 57966
455
- MJD2 57997
456
- t0geo 57981
457
- f0 29.6394
458
- f1 -3.68701e-10
459
- f2 9.29e-21
460
- RMS 0.6
461
- O J
462
- B DE200
463
- name 0531+21
464
- Notes NaN
450
+ PSR_B 0531+21
451
+ RA_hh 5
452
+ RA_mm 34
453
+ RA_ss 31.972
454
+ DEC_hh 22
455
+ DEC_mm 0
456
+ DEC_ss 52.07
457
+ MJD1 57966
458
+ MJD2 57997
459
+ t0geo 57981.0
460
+ f0 29.639378
461
+ f1 -0.0
462
+ f2 0.0
463
+ RMS 0.6
464
+ O J
465
+ B DE200
466
+ name 0531+21
467
+ Notes NaN
465
468
Name: 374, dtype: object
466
469
</ pre > </ div > </ div >
467
470
</ div >
@@ -556,6 +559,69 @@ <h2>Bayrcentric correction<a class="headerlink" href="#Bayrcentric-correction" t
556
559
</ div >
557
560
</ div >
558
561
</ div >
562
+ < div class ="section " id ="Play-with-the-profile ">
563
+ < h2 > Play with the profile< a class ="headerlink " href ="#Play-with-the-profile " title ="Permalink to this headline "> </ a > </ h2 >
564
+ < div class ="nbinput docutils container ">
565
+ < div class ="prompt highlight-none notranslate "> < div class ="highlight "> < pre > < span > </ span > [8]:
566
+ </ pre > </ div >
567
+ </ div >
568
+ < div class ="input_area highlight-ipython3 notranslate "> < div class ="highlight "> < pre > < span > </ span > < span class ="n "> profile</ span > < span class ="o "> .</ span > < span class ="n "> cycles</ span > < span class ="o "> =</ span > < span class ="mi "> 2</ span > < span class ="c1 "> # set the phase cycls of profile</ span >
569
+ < span class ="n "> plt</ span > < span class ="o "> .</ span > < span class ="n "> errorbar</ span > < span class ="p "> (</ span > < span class ="n "> profile</ span > < span class ="o "> .</ span > < span class ="n "> phase</ span > < span class ="p "> ,</ span >
570
+ < span class ="n "> profile</ span > < span class ="o "> .</ span > < span class ="n "> counts</ span > < span class ="p "> ,</ span >
571
+ < span class ="n "> yerr</ span > < span class ="o "> =</ span > < span class ="n "> profile</ span > < span class ="o "> .</ span > < span class ="n "> error</ span > < span class ="p "> ,</ span >
572
+ < span class ="n "> ds</ span > < span class ="o "> =</ span > < span class ="s1 "> 'steps-mid'</ span > < span class ="p "> )</ span >
573
+ < span class ="n "> plt</ span > < span class ="o "> .</ span > < span class ="n "> xlabel</ span > < span class ="p "> (</ span > < span class ="s2 "> "Phase"</ span > < span class ="p "> )</ span >
574
+ < span class ="n "> plt</ span > < span class ="o "> .</ span > < span class ="n "> ylabel</ span > < span class ="p "> (</ span > < span class ="s2 "> "Counts"</ span > < span class ="p "> )</ span >
575
+ </ pre > </ div >
576
+ </ div >
577
+ </ div >
578
+ < div class ="nboutput docutils container ">
579
+ < div class ="prompt highlight-none notranslate "> < div class ="highlight "> < pre > < span > </ span > [8]:
580
+ </ pre > </ div >
581
+ </ div >
582
+ < div class ="output_area docutils container ">
583
+ < div class ="highlight "> < pre >
584
+ Text(0, 0.5, 'Counts')
585
+ </ pre > </ div > </ div >
586
+ </ div >
587
+ < div class ="nboutput nblast docutils container ">
588
+ < div class ="prompt empty docutils container ">
589
+ </ div >
590
+ < div class ="output_area docutils container ">
591
+ < img alt ="../_images/notebooks_PulsarTimingAnalysis_16_1.png " src ="../_images/notebooks_PulsarTimingAnalysis_16_1.png " />
592
+ </ div >
593
+ </ div >
594
+ < div class ="nbinput docutils container ">
595
+ < div class ="prompt highlight-none notranslate "> < div class ="highlight "> < pre > < span > </ span > [20]:
596
+ </ pre > </ div >
597
+ </ div >
598
+ < div class ="input_area highlight-ipython3 notranslate "> < div class ="highlight "> < pre > < span > </ span > < span class ="n "> profile</ span > < span class ="o "> .</ span > < span class ="n "> norm</ span > < span class ="p "> (</ span > < span class ="n "> method</ span > < span class ="o "> =</ span > < span class ="mi "> 1</ span > < span class ="p "> )</ span > < span class ="c1 "> # Normalization</ span >
599
+ < span class ="n "> profile_new</ span > < span class ="o "> =</ span > < span class ="n "> profile</ span > < span class ="o "> .</ span > < span class ="n "> resample</ span > < span class ="p "> (</ span > < span class ="n "> sample_num</ span > < span class ="o "> =</ span > < span class ="mi "> 1000</ span > < span class ="p "> )</ span >
600
+
601
+ < span class ="n "> plt</ span > < span class ="o "> .</ span > < span class ="n "> imshow</ span > < span class ="p "> (</ span > < span class ="n "> profile_new</ span > < span class ="p "> ,</ span > < span class ="n "> aspect</ span > < span class ="o "> =</ span > < span class ="s1 "> 'auto'</ span > < span class ="p "> ,</ span >
602
+ < span class ="n "> extent</ span > < span class ="o "> =</ span > < span class ="p "> [</ span > < span class ="mi "> 0</ span > < span class ="p "> ,</ span > < span class ="mi "> 2</ span > < span class ="p "> ,</ span > < span class ="mi "> 0</ span > < span class ="p "> ,</ span > < span class ="mi "> 1000</ span > < span class ="p "> ])</ span >
603
+ < span class ="n "> plt</ span > < span class ="o "> .</ span > < span class ="n "> xlabel</ span > < span class ="p "> (</ span > < span class ="s2 "> "phase"</ span > < span class ="p "> )</ span >
604
+ < span class ="n "> plt</ span > < span class ="o "> .</ span > < span class ="n "> ylabel</ span > < span class ="p "> (</ span > < span class ="s2 "> "numbers of profiles"</ span > < span class ="p "> )</ span >
605
+ </ pre > </ div >
606
+ </ div >
607
+ </ div >
608
+ < div class ="nboutput docutils container ">
609
+ < div class ="prompt highlight-none notranslate "> < div class ="highlight "> < pre > < span > </ span > [20]:
610
+ </ pre > </ div >
611
+ </ div >
612
+ < div class ="output_area docutils container ">
613
+ < div class ="highlight "> < pre >
614
+ Text(0, 0.5, 'numbers of profiles')
615
+ </ pre > </ div > </ div >
616
+ </ div >
617
+ < div class ="nboutput nblast docutils container ">
618
+ < div class ="prompt empty docutils container ">
619
+ </ div >
620
+ < div class ="output_area docutils container ">
621
+ < img alt ="../_images/notebooks_PulsarTimingAnalysis_17_1.png " src ="../_images/notebooks_PulsarTimingAnalysis_17_1.png " />
622
+ </ div >
623
+ </ div >
624
+ </ div >
559
625
</ div >
560
626
561
627
@@ -570,7 +636,7 @@ <h2>Bayrcentric correction<a class="headerlink" href="#Bayrcentric-correction" t
570
636
571
637
< div role ="contentinfo ">
572
638
< p > © Copyright 2022, Youli.
573
- < span class ="lastupdated "> Last updated on 2023 Jun 06 at 12:09:33 UTC.
639
+ < span class ="lastupdated "> Last updated on 2023 Jun 07 at 15:30:44 UTC.
574
640
</ span > </ p >
575
641
</ div >
576
642
0 commit comments