How to use TEMPO2.ppt

tempo2tutorial georgehobbs 饺子 researchscientistcsiroaustraliatelescopenationalfacility internationalyoungscientistofntscgeorge hobbs csiro au content tempo2isasoftwarepackage itimplements pulsartiming algorithms section1 whattempo2hasbeenusedfor overviewoftempo2section2 howtousetempo2section3 howtodeveloptempo2majortempo2developers softwareandalgorithms g hobbs r edwards r manchester w coles x you m keith f jenet d yardley j verbiest csiro tempo2 whathastempo2beenusedforin2010 2011 searchingforgravitationalwavesignals e g yardleyetal 2011 yardleyetal 2010 abbottetal 2010 vanhaasterenetal 2010 usingpulsarsasnavigationalaids ruggiero2011 statisticalanalysisoftimingresiduals e g na x s etal 2011 hobbsetal 2010 studyingemissiongeometry e g noutsosetal 2011 searchingforgamma raypulsars rayetal 2010 determiningpulsarmasses demorestetal 2010 studyingpulsationsfrommain sequencestars ravietal 2010 testsofrelativity weisbergetal 2010 measuringjupiter smass championetal 2010 observationsofglitches chukwudeetal 2010 analysisofaccretingmillisecondx raypulsar patrunoetal 2010 relativisticspinprecession manchesteretal 2010 csiro tempo2 whatdoestempo2do usage1 startwithfileofpulsearrivaltimesstartwithmodelofapulsar sspin astrometricandorbitalparameters pulsartimingresidualsusage2 startwithpulsartimingmodel predictthepulsephaseandfrequencyatagivenobservatoryatagiventime csiro tempo2 aparameterfile par csiro tempo2 pulsarname astrometric rotational dispersionmeasureandorbitalparameters realisationofterrestrialtime solarsystemephemeris weightedfit fitforthisparameter anarrivaltimefile tim csiro tempo2 statesthatthisisatempo2formatfile filenameoridentifier observingfrequency mhz sitearrivaltime mjd uncertaintyonarrivaltime us telescopecode userdefinedflags howdoestempo2work thealgorithmsarebasedonthoseimplementedintempo1detailsinhobbs edwards manchester 2006 andedwards hobbs manchester 2006 csiro tempo2 conversionofsite arrival timetopulseemissiontime clockcorrections atmosphericdelays einsteindelay shapirodelay roemerdelay secularmotion e g radialvelocity dispersivedelay orbitalmotion details clockcorrection mustconvertfromtheobservatorytimestandardtoarealisationofterrestrialtimett usesetoftextfilescontainingthedifferencebetweentwotimestandards pks2gps clk gps2utc clk utc2tai clk tai2tt bipm2011 csiro tempo2 you llfindthesefilesin tempo2 clockdirectory details einsteindelay pulsearrivaltimes interrestrialtime mustbeconvertedtotheframeofthesolarsystembarycentre tempo2usesthetabulatedresultsofirwin fukushima 1999 csiro tempo2 details roemerdelay theroemerdelayisthevacuumlighttraveltimebetweenthepulsearrivingattheobservatoryandtheequivalentarrivaltimeatthessb csiro tempo2 r k ssb earth pulsar k unitvectorpointingatpulsar obtainedfromra dec pmra pmdec risvectorfromobservatorytossb usejplephemeris observatorycoordinatesandearthorientationparameters useobservatorycoordinatefiletoknowpositionofobservatory roemerdelay csiro gravitationalwavedetection pulsarinecliptic eclipticlatitude 90deg details shapirodelay timedelaycausedbythepassageofthepulsethroughlargegravitationalfields tempo2includesdelaycausedbysun 110us jupiter 180ns saturn 58ns neptune 12ns anduranus 10ns csiro tempo2 details dispersivedelay disthedispersionconstant dm cm 3pc 2 410 x10 16dfisthefrequencyoftheradiationatthesolarsystembarycentre csiro tempo2 details binarysystem forpulsarsinbinarysystems tempo2includesparametersdescribingtheorbitalmotion variousbinarymodelsexist suggestusing t2 modelthatcombinesmostearlierbinarymodels csiro tempo2 usingthepulsartimingmodel havepulseemissiontimeinthepulsarframe predictusingthepulsartimingmodelcanalsoincludesimplemodelofglitchevents csiro tempo2 phaseofpulsesequence pulsefrequency andtimederivatives pulseemissiontime timeatwhichdf dt n referencephase timingresiduals csiro tempo2 timingresidualfori thobservation pulsephase nearestintegertofi pulsefrequency example csiro tempo2 time measuredarrivaltimewithobservatoryclock measuredarrivaltimeinterrestrialtime dc measuredarrivaltimeinpulsarframe 1 3 4 2 5 predictedarrivaltimesinpulsarframe timingresidual pulsartimingresiduals ifpulsarmodelpredictstheobservationsperfectly andtheconversionfromtheobservatorytopulsarframeisperfect thenr 0 withinmeasurementuncertainty ifr 0thenthepulsarmodelis 1 notaccurateor 2 doesnotincludeaphysicalprocessthataffectsthemeasuredarrivaltimesor 3 thecorrectionfromtheobservatorytopulsarframeisnotcorrect csiro tempo2 pulsartimingresiduals incorrectf0 csiro tempo2 pmodel ptrue pulsartimingresiduals incorrectf1 csiro tempo2 pmodel ptrue pmodel ptrue pulsartimingresiduals incorrectposition csiro tempo2 pulsartimingresiduals incorrectpropermotion csiro tempo2 pulsartimingresiduals csiro tempo2 goodparameters donothaveaphaseconnectedsolution phasewraps inducedtimingresidualsfromismandfromgravitationalwaves csiro tempo2 50cm 600mhz 20cm 1400mhz 10cm 3000mhz inducedresidualsfromgravitationalwavebackground incorrectconversiontopulsarframe csiro tempo2 championetal 2010 inducedtimingresidualsfromerrorinmassofjupiterof5x10 10mo predictivemode tempo2canprovideapolynomialapproximationofthephaseandpulsefrequencyattheobservatoryforaspecifiedtimeinterval tempo1createspolycos tempo2cancreatepolycos butwerecommendpredictors predictorsmodelthephaseasafunctionoftimeusingchebyshevpolynomials tempo2providesa clibrary toreadpredictorstempo2 fmypar par pred sitenamemjd1mjd2freq1freq2ntimecoeffnfreqcoeffseg length csiro tempo2 observingfrequencyrange mhz timerange segmentlength s part2 usageoftempo2andexampleofpulsartimingtempo2iswritteninc c withalittlefortran andrunsonlinuxandmacos maintempo2website http www atnf csiro au research pulsar tempo2 see tutorials and documentation maintempo2downloadsite csiro tempo2 runningtempo2 needapulsarmodel usetheatnfpulsarcatalogue http www atnf csiro au research pulsar psrcat typeinpulsarnameandclick getephemeris getsomepulsearrivaltimes useparkesdataarchive seejonathan stalkand pat inpsrchive seewillem stalk tempo2 fpsrname parpsrname tim csiro tempo2 usingthe plk graphicalinterface tempo2 grplk fmypar parmytim tim csiro tempo2 whichparameterstofit re dothefit whattoplot theresiduals press h togethelp otherusefulgraphicalinterfaces spectrum csiro tempo2 otherusefulgraphicalinterfaces plotmany csiro tempo2 tempo2 grplotmany fpar1 partim1 tim fpar2 partim2 tim fpar3 parpar3 tim theglitchplugin csiro tempo2 yumeng inpreparation glitchevents pulsefrequency pulsefrequencywithmeanremoved pulsefrequencyderivative using output interfaces tempo2 outputgeneral s hellof0 f0 v f0 e nf1 f1 v f1 e n fmypar parmytim timtempo2 outputgeneral2 s freq freq sat sat bat bat res pre n fmypar parmytim timtempo2 outputmatrix fmypar parmytim tim csiro tempo2 part3 developingthesoftware youshouldneverneedtomodifythemaintempo2sourcecode pleasereportanybugstogeorge hobbs csiro au youshouldcreatenew plug inpackages thesecanbefortextoutput agraphicalinterface methodstoselectdataornewfittingalgorithms csiro tempo2 creatinganewtextoutputformat printoutthepulsarmodelparametersandthetimingresidualsinyourformat csiro tempo2 mustincludetempo2 h musthavethisfunction pulsarnumber parameterofinterest timederivative observationnumber observationparameter numberofobservations creatinganewgraphicalinterface csiro tempo2 mustincludethisfunction mustprocessthecommandlinearguments readparameterandobservationfile s formbarycentricarrivaltimes residualsandfit repeatwithpost fitparameters doyourplotting analysisetc mustdefinehowmanypulsarstoprocess outputformatversusgraphicalinterface outputformat tempo2runsasusualthengivescontroltoyourplugintodisplaytheoutput easy butnocontroloverfittingetc graphicalinterface tempo2providescontroltoyourpluginfromthebeginning yourpluginmustcallallthefunctionstoreaddatafiles dothefitetc butyouhavemorecontrol csiro tempo2 thetempo2toolkits toaiddevelopmentthefollowingtoolkitsareavailable t2toolkit h findminimumandmaximumvaluesinarray findrms mean range sortdataarray createrandomnumberstkspectrum h routinesforspectralestimationtkfit h fittingroutines e g svd choleskydecomposition tempo2util h calculatedotproducts turnshh mm ss dd mm ssintodegreesgwsim h gravitationalwavesimulationroutines note manypeopleuse numericalrecipes routinesinplugins thesearecopyrightandthesourcecodecannotbemadepublicallyavailable usethetempo2toolkitsinstead csiro tempo2 debuggingyourinstallation sometimesthingsdon twork ifitdoesn tcompile checkthatyourversionofgfortran gccandg arethesame checkthatpgplothasbeencorrectlyinstalled itcompilesandruns butdoesn tdoanything runyourcodewith tempo2 fmypar parmytim tim debug displayversion outfileandthensendyour par timandoutfiletogeorge hobbs csiro au oranytempo2developer yourplugindoesnotwork checkthatyouhaveinstalledthepluginin tempo2 plugins csiro tempo2 publishingresultsfromtempo2 forgeneralusage pleasecite hobbs edwards manchester 2006 fordetailsofthealgorithms edwards hobbs manchester 2006 forgravitationalwaveimplementation hobbsetal 2009 toenablecompletereproducibilityofresultsshouldincludeanappendixincludingoutputof displayversion important publishtheclockreference e g tt tai andsolarsystemephemerismodel e g de405 csiro tempo2 topicsnotcoveredinthispresentation simulatingdata usefake gwbkgrd plugins useof flags use select pass filteroptions atmosphericcorrectionsxiaopen