Skip to content

Commit

Permalink
Add an alternative model to the figure.
Browse files Browse the repository at this point in the history
  • Loading branch information
Fabian Kislat committed Jul 24, 2017
1 parent d172cba commit 1c00b54
Show file tree
Hide file tree
Showing 2 changed files with 20 additions and 7 deletions.
3 changes: 3 additions & 0 deletions crab-twopole-caustic-model-incl65deg.gpstmodel
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
phase = [ -0.2 -0.195 -0.19 -0.185 -0.18 -0.175 -0.17 -0.165 -0.16 -0.155 -0.15 -0.145 -0.14 -0.135 -0.13 -0.125 -0.12 -0.115 -0.11 -0.105 -0.1 -0.095 -0.09 -0.085 -0.08 -0.075 -0.07 -0.065 -0.06 -0.055 -0.05 -0.045 -0.04 -0.035 -0.03 -0.025 -0.02 -0.015 -0.01 -0.005 1.66533453694e-16 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 0.055 0.06 0.065 0.07 0.075 0.08 0.085 0.09 0.095 0.1 0.105 0.11 0.115 0.12 0.125 0.13 0.135 0.14 0.145 0.15 0.155 0.16 0.165 0.17 0.175 0.18 0.185 0.19 0.195 0.2 0.205 0.21 0.215 0.22 0.225 0.23 0.235 0.24 0.245 0.25 0.255 0.26 0.265 0.27 0.275 0.28 0.285 0.29 0.295 0.3 0.305 0.31 0.315 0.32 0.325 0.33 0.335 0.34 0.345 0.35 0.355 0.36 0.365 0.37 0.375 0.38 0.385 0.39 0.395 0.4 0.405 0.41 0.415 0.42 0.425 0.43 0.435 0.44 0.445 0.45 0.455 0.46 0.465 0.47 0.475 0.48 0.485 0.49 0.495 0.5 0.505 0.51 0.515 0.52 0.525 0.53 0.535 0.54 0.545 0.55 0.555 0.56 0.565 0.57 0.575 0.58 0.585 0.59 0.595 0.6 0.605 0.61 0.615 0.62 0.625 0.63 0.635 0.64 0.645 0.65 0.655 0.66 0.665 0.67 0.675 0.68 0.685 0.69 0.695 0.7 0.705 0.71 0.715 0.72 0.725 0.73 0.735 0.74 0.745 0.75 0.755 0.76 0.765 0.77 0.775 0.78 0.785 0.79 0.795 0.8 ]
fraction = [ 0.72 0.719872896904 0.719747872367 0.719625140355 0.719506712515 0.71940178966 0.719307189755 0.719232658412 0.719182242154 0.719152790624 0.718590476325 0.674205709648 0.671989734944 0.681968679307 0.701412403593 0.719384886603 0.722971283192 0.72745439391 0.729538338688 0.730603294737 0.731219381995 0.73172791931 0.732133367276 0.732391227926 0.7325065246 0.732963242212 0.733353120252 0.733681716607 0.733889085766 0.733853021739 0.733800299035 0.734107922832 0.734152711445 0.661886090624 0.494346117694 0.351151257539 0.326440156551 0.309265687283 0.293866160869 0.280451493676 0.269624596957 0.272583510179 0.297644529441 0.317208443774 0.329519980308 0.269623379451 0.285002319794 0.301900415434 0.316880953379 0.322565563696 0.328353013678 0.332874552354 0.334002043883 0.335299740092 0.336765219092 0.338395758857 0.337458218795 0.336508823227 0.335676534365 0.334972649128 0.334398267477 0.333954294209 0.33364143366 0.333055834663 0.331711390235 0.33050342587 0.329445371895 0.328556467035 0.327832213416 0.327274383641 0.326884422443 0.326334877438 0.325691481304 0.32522976333 0.324951036478 0.325002990603 0.32526945019 0.325747970473 0.326437680724 0.327337169177 0.328696232275 0.330670690477 0.331569494526 0.325054640584 0.319404386097 0.314627956894 0.310729059342 0.307705275268 0.30548339807 0.301412305248 0.29863045032 0.297092470199 0.296741031487 0.297508451807 0.300344751513 0.305875762033 0.312403650202 0.319444634119 0.326761510853 0.3350929832 0.344267127331 0.354129359462 0.364408175429 0.375241401929 0.385734077167 0.396022326493 0.406605515847 0.417605124628 0.428929432874 0.440018914791 0.450522982189 0.461046184412 0.471693565766 0.482692883253 0.493587883469 0.502964824754 0.510919023482 0.51913156439 0.527220488574 0.535163497507 0.541648564713 0.532267788131 0.517251289104 0.505748303378 0.505860133889 0.52993739332 0.623499828738 0.642287074893 0.640305374362 0.637422277897 0.633517508896 0.629629638956 0.627770557231 0.629080692119 0.633353775135 0.640174649879 0.648972387628 0.657573572742 0.665071578584 0.668250664912 0.671265340453 0.674283756806 0.67730589851 0.680328406201 0.683208591894 0.685384969251 0.68756614675 0.689752147716 0.691942995649 0.694153289148 0.696424397321 0.698687656971 0.700942854275 0.703189776651 0.705428212741 0.707657952397 0.709883639057 0.712106053924 0.714312807454 0.716358040003 0.716357209256 0.716356871956 0.716357028875 0.716357680671 0.716358150125 0.716355544194 0.716353316435 0.716351467423 0.716349997669 0.716348907625 0.716348197683 0.716327708753 0.716306307177 0.716284906602 0.716263506854 0.71624210812 0.7162207104 0.716199313693 0.716177917999 0.716155933184 0.716133257629 0.716110569367 0.716087868407 0.716065154757 0.716042428426 0.716019689422 0.715996937755 0.715974173432 0.715951396462 0.715928606855 0.715905804617 0.715882989759 0.715860162288 0.715837322213 0.715814469544 0.715791604288 0.715768726454 0.715745836051 0.715722933087 0.715700017572 0.715677089514 ]
angle = [ 121.0 120.9988334 120.997496889 120.99599473 120.994211543 120.991678895 120.988799802 120.985325493 120.981423531 120.97743984 120.956313027 119.04651945 118.290535425 117.987838611 118.384737777 119.30103523 119.529363401 120.067958867 120.451572101 120.72501826 120.874661812 121.037747384 121.21474189 121.405520763 121.602665014 121.821294161 122.069787066 122.354073143 122.674665238 123.017281643 123.342600981 123.771388862 124.219781939 125.104237048 130.125433114 133.607573239 134.447204347 135.913767511 137.49457804 139.081243038 140.379080194 139.426007728 135.899778486 135.868423817 135.428457818 130.513196886 129.217773066 128.048035767 127.1195426 126.592121694 126.027366167 125.532027512 125.142214485 124.757017315 124.377101232 124.003095883 123.676287718 123.342315427 123.00166629 122.658133516 122.312185825 121.964304013 121.614978665 121.265306391 120.913173405 120.555982408 120.181134916 119.7877376 119.390780338 118.990966497 118.589020527 118.180458861 117.764896503 117.346866146 116.9272059 116.476057452 116.025284283 115.576446735 115.130596854 114.6887657 114.261873404 113.859538664 113.402599191 112.525200538 111.613119092 110.672462884 109.710164432 108.733811974 107.745717014 106.489250656 105.229474544 103.980991785 102.757847031 101.572900119 100.424333127 99.3376152744 98.3500298424 97.3639375527 96.2859136592 95.3378805869 94.5169773348 93.8177082537 93.0807470045 92.3843638573 91.7814562312 91.2639210713 90.8460531567 90.3854181792 89.970488062 89.6342973363 89.3459876457 89.1350148492 88.9534538133 88.744941923 88.6128008388 88.4407625246 88.1839685107 87.922046423 87.7144617526 87.5556585634 87.6149356213 89.2873592629 92.0356726967 95.8966865399 101.173289519 108.079309618 116.565033023 117.723579482 118.001507698 118.33551205 118.74463611 119.187466537 119.706815039 120.233212823 120.727239603 121.156073776 121.496489267 121.756060343 121.962346057 121.919018714 121.867786663 121.81618279 121.764210984 121.711734329 121.659041254 121.619803735 121.580509363 121.541158374 121.501750998 121.462883444 121.426206569 121.388898606 121.350969277 121.312428235 121.273285059 121.233549258 121.193404299 121.152883242 121.111563979 121.072370191 121.073550396 121.074728933 121.075905496 121.077079778 121.078231535 121.079278644 121.080324838 121.08136992 121.082413692 121.083455955 121.08449651 121.085018628 121.085517091 121.086015566 121.086514056 121.087012558 121.087511071 121.088009595 121.08850813 121.088991394 121.089456727 121.089921673 121.090386233 121.090850407 121.091314194 121.091777593 121.092240605 121.092703227 121.093165461 121.093627306 121.094088761 121.094549825 121.095010499 121.095470782 121.095930673 121.096390172 121.096849278 121.097307991 121.097766311 121.098224237 121.098681768 ]
24 changes: 17 additions & 7 deletions dataCrabPhase.gpst
Original file line number Diff line number Diff line change
Expand Up @@ -9,18 +9,24 @@ GPST 14
renormalize = [ 1000., 100, 1000 ]

; Place a legend on the flux panel
legend = [ 0.15 0.60 0.55 0.86 ]
legend = [ 0.15 0.51 0.55 0.86 ]

; Uncomment this line in order to get rid of the version / data information on
; the plot
;clean = true

; Name the simulation. This will be the label of the model curve in the legend
name = "Crab Nebula"
name = "Outer gap model"

; Set the color of the model curve
color = kGreen + 1

; Put a label for the observed data on the legend
datalabel = "Obs. Time: 10 days"

; Set the line thickness of the data points
width = 2

; Observation time. Will be multiplied by the 'ontime' value. The amego.det
; file included above initializes that to 0.1. But you could change it below.
time = 10
Expand All @@ -47,13 +53,13 @@ fluxunit = cgse9

; Set the range of the flux axis
fluxmin = 3.1
fluxmax = 14
fluxmax = 15

; Set the maximum of the polarization fraction axis
fracmax = 0.99

; Set the range of the polarization angle axis
chimin = -25
chimin = -35
chimax = 65

; Shift the polarization angle so that it is centered around 0
Expand All @@ -62,15 +68,19 @@ deltachi = -90
; Include the detailed model file
@include "crab-outer-gap-model-incl65deg.gpstmodel"

; Set the color of the model curve
color = kGreen + 1

; Define the spectrum. In case of a phase-dependent simulation there is
; no polarization information here.
component = "Crab nebula"
energy = [ 100 9999 ]
flux = [ 8.155e-9 4.087e-9 ]

; Show an alternative model. Note that this model file does not have a
; flux parameter since the light curve is exactly the same as the one used
; for the simulation.
model = "Two-pole caustic model"
@include "crab-twopole-caustic-model-incl65deg.gpstmodel"
color = kBlue + 1

; Show MDP in each phase bin, show it as a thick dark red line
showmdp = true
color = kRed + 2
Expand Down

0 comments on commit 1c00b54

Please sign in to comment.