aryingN_distributions"}](SpeedDistributions-VaryingN.pdf){width="49.00000%"}

Varying the number of basis functions
-------------------------------------

We first investigate how the reconstructed WIMP mass $m_\textrm{rec}$ and uncertainty varies with the number of basis functions $N$. For now, we fix our choice of basis to shifted Legendre polynomials, as used in Paper 1:

$$P_k(v) = L_k\left(2\frac{v}{v_\textrm{max}} - 1\right)\,,$$

where $L_k$ is the Legendre polynomial of order $k$, and $v_\textrm{max}$ is a cut off for the parametrization. We should choose $v_\textrm{max}$ to ensure that $f_1(v)$ is negligible above the cut off. However, too high a choice of $v_\textrm{max}$ will result in $f_1(v)$ being close to zero over a large range of the parametrization, making fitting more difficult. We use the value $v_\textrm{max} = 1000 {\,\textrm{km s}^{-1}}$, while lies significantly above the Galactic escape speed.

The lower panel of Fig. \[fig:BUMP\_LEG\] shows the best fit mass and 68% confidence intervals as a function of $N$, using as input a WIMP of mass 50 GeV and the ‘bump’ distribution function. The reconstructed mass very rapidly settles close to the true value, using as few as three basis functions. This is because adding the bump near $v \sim 500 {\,\textrm{km s}^{-1}}$ still leaves the mean inverse speed relatively smooth, so a large number of basis function are not required. The correct mass is reconstructed and we emphasize in the lower panel of Fig. \[fig:BUMP\_LEG\] that the reconstruction is stable with the addition of more basis functions.

We should also consider how the quality of the fit changes as a function of $N$. We would expect that adding fit parameters should always lead to a better fit. Eventually, the fit should be good enough that adding additional basis functions will no longer improve it significantly. We can then be confident that our reconstruction is accurate and not an artifact of using too few basis functions. In order to investigate this, we utilise the Bayesian Information Criterion (BIC) [@Schwarz:1978], which is
aryingN_distributions"}](SpeedDistributions - VaryingN.pdf){width="49.00000% " } 

 Varying the number of basis functions 
 ------------------------------------- 

 We foremost investigate how the reconstructed WIMP bulk $ m_\textrm{rec}$ and uncertainty varies with the numeral of footing functions $ N$. For now, we fix our option of footing to shifted Legendre polynomials, as used in Paper 1: 

 $ $ P_k(v) = L_k\left(2\frac{v}{v_\textrm{max } } - 1\right)\,,$$ 

 where $ L_k$ is the Legendre polynomial of order $ k$, and $ v_\textrm{max}$ is a cut off for the parametrization. We should choose $ v_\textrm{max}$ to see that $ f_1(v)$ is negligible above the cut off. However, too gamey a option of $ v_\textrm{max}$ will result in $ f_1(v)$ being close to zero over a big range of the parametrization, making match more difficult. We use the value $ v_\textrm{max } = 1000 { \,\textrm{km s}^{-1}}$, while lies importantly above the Galactic escape focal ratio. 

 The lower jury of Fig.   \[fig: BUMP\_LEG\ ] shows the best fit mass and 68% confidence intervals as a affair of $ N$, using as input a WIMP of mass 50 GeV and the ‘ bump ’ distribution function. The reconstructed mass very rapidly settles close to the true value, using as few as three basis functions. This is because adding the bump near $ v \sim 500 { \,\textrm{km s}^{-1}}$ however leaves the mean inverse amphetamine relatively legato, so a large number of basis affair are not required. The correct mass is reconstructed and we stress in the lower panel of Fig.   \[fig: BUMP\_LEG\ ] that the reconstruction is stable with the addition of more basis functions. 

 We should besides consider how the quality of the fit changes as a function of $ N$. We would expect that adding fit parameters should constantly lead to a better fit. Eventually, the fit should be effective enough that adding additional basis function will no longer improve it significantly. We can then be confident that our reconstruction is accurate and not an artifact of using besides few basis affair. In order to investigate this, we utilise the Bayesian Information Criterion (BIC) [ @Schwarz:1978 ], which is
arylngN_distributions"}](SpeedDirtributions-VarynbgN.pdf){xidth="49.00000%"}

Vadying thd number of basis functions
-------------------------------------

Wx fiest ibvestigate how the recunstructef WIMP mqss $n_\rextrm{rec}$ ehd uncevcaintg varnew with the numner of basiv functions $N$. Xof uow, we fix our choice of basis to shyfted Lrgfndre polynomiwls, sf ussd in Paper 1:

$$P_k(v) = L_k\left(2\frac{v}{v_\textrj{max}} - 1\rpght)\,,$$

where $L_k$ is tne Legendre polynomial of lrdeg $k$, and $v_\textrm{max}$ is a cut odf fje the paramegrization. We should chkose $v_\textrm{max}$ to ensure that $w_1(v)$ is negligiblg zblee the cut iff. Hjwever, too hlbh a cvoice og $v_\textrm{max}$ wlll rxsulr in $f_1(v)$ being close tm zero over a largg range of tke parametrization, majibg fijting mord diwfidukt. We usf tie value $v_\tsxtrm{max} = 1000 {\,\rextrm{km s}^{-1}}$, while lits fpbnificantly zbove ehq Galactic escape speed.

The lower panel mf Rig. \[fig:BUMP\_LEG\] shows the best fit mass and 68% clnfidence intervals as a function of $N$, using as input a WIKP of oasw 50 EwV and the ‘bump’ distribution function. The reconffricned mass very rapldly settles close tl yre true value, using as few as three basid functyons. Rhis is btcausr adding the bump near $v \sin 500 {\,\textrm{km f}^{-1}}$ still leaves the lean inversz speec relstively smooth, so a larye numger of basid functioh are not requirea. Tme worrect mass is reconstruceed and wx empkasize iv thg lower panel of Vig. \[fin:TUMP\_LEG\] that the rfconsjructimn is stabpe with the addition of more basis functions.

Wg svound also eonsidcr how the qualyty of the fit changef as x function of $N$. Wx would expest that addinc fit parametxrs shoulq alqays lead tu a better fit. Rventuallj, uhe fit shiuld be good enougm thaj zdding additioncu basis functionx wkll nl kogcer improve ht skgnkgicangly. We can bhev be confident that our secohstruction is accutabe and noj an artisact of using too few basis funftionv. Ii ordet tj investigate this, we utilise fhe Bayeslan Information Criberijn (BIC) [@Schcarz:1978], which is
aryingN_distributions"}](SpeedDistributions-VaryingN.pdf){width="49.00000%"} Varying the number of basis functions first how the WIMP mass $m_\textrm{rec}$ number basis functions $N$. now, we fix choice of basis to shifted Legendre as used in Paper 1: $$P_k(v) = L_k\left(2\frac{v}{v_\textrm{max}} - 1\right)\,,$$ where $L_k$ is Legendre polynomial of order $k$, and $v_\textrm{max}$ is a cut off for the We choose to that $f_1(v)$ is negligible above the cut off. However, too high a choice of $v_\textrm{max}$ will in $f_1(v)$ being close to zero over a range of the parametrization, fitting more difficult. We use value = 1000 s}^{-1}}$, lies above the Galactic speed. The lower panel of Fig. \[fig:BUMP\_LEG\] shows the best fit mass and 68% confidence intervals as function of as input WIMP mass GeV and the function. The reconstructed mass very rapidly the true value, using as few as three functions. This because adding the bump near $v 500 {\,\textrm{km s}^{-1}}$ still leaves the mean inverse relatively smooth, so a large number of basis function are not required. The correct mass and we emphasize in lower panel of \[fig:BUMP\_LEG\] the is with the of more basis functions. We should also consider how the quality the fit changes as a function of $N$. We would adding parameters should always to a better fit. the should be good enough additional functions improve significantly. can then be confident our reconstruction is accurate and an artifact of using order to investigate this, we utilise the Bayesian Criterion (BIC) [@Schwarz:1978], which is
aryingN_distributions"}](SpeeddistributiOns-VaRyiNgN.PdF){widTh="49.00000%"}

VaRying the number OF basIs functions
-------------------------------------

We first inveStigaTe HOw thE ReConstRucted WimP MASs $m_\TeXtRm{rEc}$ ANd UncerTaiNty variEs with the nUmbEr Of basis functIOnS $N$. For now, we Fix Our choice of bAsiS to shiFtEd LEGendrE poLynomIals, as USed in PAper 1:

$$P_k(v) = L_k\LeFT(2\frac{v}{V_\Textrm{mAX}} - 1\RiGht)\,,$$

wHere $L_k$ is the LegendRE pOLynomial of ordeR $k$, and $v_\TeXTrM{MAx}$ iS a cUt off for thE pArameTRizatioN. we SHOUld CHoose $v_\textrm{mAx}$ to ensure tHAt $f_1(V)$ is negLiGibLE above The cuT oFF. HoWever, too higH a chOice of $v_\teXtrm{maX}$ Will resULt in $f_1(v)$ bEing clOse To zEro oVEr A lArgE rANge OF tHe pARamEtrizatiOn, MaKing fIttiNG MORe diFfiCult. we use The value $v_\textRm{mAx} = 1000 {\,\teXTrm{Km s}^{-1}}$, whIle liEs siGnIficaNtly abOve thE GAlactic escape spEed.

THe lower paNel Of fig. \[FiG:BUMP\_leG\] showS thE beSt fit maSs and 68% coNFidEnCE INtErvals as a function oF $N$, USInG as input A WIMP oF MaSs 50 gEV and the ‘BuMp’ dIstrIBUtion FuncTIoN. The recoNstrucTEd MaSs very rApIdly seTtLes CloSe to tHE truE value, Using as fEw as tHRee basis functiONs. This is becauSE aDDInG The bUmp Near $v \sim 500 {\,\texTrm{kM S}^{-1}}$ stiLl leAVeS thE Mean iNversE sPEeD Relatively smooth, so a LaRge numBer of Basis function Are not requIRED. The corrEct mASs IS reconstructed And we Emphasize iN The lower Panel Of Fig. \[fig:bUMP\_LEG\] thAT The reconStrUctIon Is sTABlE with the additION of mOrE basis fUncTions.

We ShoUld AlsO coNsIder how thE quality Of ThE fIt ChaNges aS A functioN oF $N$. WE wOulD expeCT that aDding Fit pArAmETerS should ALwAYS leaD tO a BettEr fIt. eventUallY, The Fit shouLd be good eNouGH thaT aDdIng addiTional basis fuNcTions will nO lOngEr imprOVE it signiFicantly. We can then be confIDent thaT ouR recoNstrUction is aCcuRate anD noT An artiFact of Using ToO feW BAsis fUNCtIonS. IN order to inVEStiGate tHiS, we uTilise tHe Bayesian InformatIOn CRiterion (BIC) [@ScHwaRz:1978], whICH iS
aryingN_distributions"}](S peedDistri butio ns- Var yi ngN. pdf) {width="49.000 0 0%"} 

Varying the number o f bas is func t io ns
-- ------- - -- - - --- -- -- --- -- - -- ----- --- 

We fi rst invest iga te  how the rec o ns tructed WI MPmass $m_\tex trm {rec}$  a ndu ncert ain ty va ries w i th the  number o fb asis f u nctions $ N$ . Fo r now, we fix our ch o ice of basis t o shif te d  L e g end repolynomial s,  as u s ed in P a pe r 1 :
$ $P_k(v) = L_k \left(2\fra c {v} {v_\te xt rm{ m ax}} -  1\ri gh t )\, ,$$

where$L_k $ is theLegend r e polyn o mial of  order  $k $,and$ v_ \t ext rm { max } $isa  cu t off fo rth e par amet r i z a tion . W e sh ouldchoose $v_\te xtr m{ma x }$to en surethat  $ f_1(v )$ isnegli gi ble above the c ut o ff. Howev er,  t oohi gh ac hoiceof$v_ \textrm {max}$w ill  r e s u lt  in $f_1(v)$ being  c l o se  to zero  overa  l ar g e rangeof  th e pa r a metri zati o n,  makingfittin g  m or e diffi cu lt. We  u sethe  valu e  $v_ \textr m{max} =  1000 {\,\textrm{kms }^{-1}}$, whi l el i es sign ifi cantly abov e th e  Gal acti c  e sca p e spe ed.
Th e  l o wer panel of Fig. \ [f ig:BUM P\_LE G\] shows the  best fitm a s s and 68 % co n fi d ence intervals  as a  functiono f $N$, u singas input  a WIMP o f mass 50GeV  an d t he‘ b um p’ distributi o n  fun ct ion. Th e r econstr uct edmas s v er y rapidly  settles  c lo se  t o t he tr u e value,  u sin gasfew a s  three  basi s fu nc ti o ns.  This i s  b e c ause  a dd ingthe  b ump n ear$ v \ sim 500  {\,\text rm{ k m s} ^{ -1 }}$ sti ll leaves the  m ean invers espe ed rel a t ively sm ooth, so a large number of basi s f uncti on a re not re qui red. T hec orrect  massis re co nst r u cteda n dweem phasize in t helower  p anel  of Fig . \[fig:BUMP\_LEG\ ]  th at the recons tru ctio n is  st a bl e  wi th the a ddition of more  basis fun ct i on s.

We sho u ldal so cons ider ho w the quality  of the f it change sas a f unc tion of $N $. We wo uld expec t  that ad dingfit  param et ers  shou ld alw a ysleadto a b et ter fi t. Ev en tually,the fit should be goodenough  that  ad ding addi tio n albasis fun ctio ns will no  lo nge r imp rov e  it s igni f ic ant l y. We  can then be c o nf ide n t  t hat our rec o n s tru ction  is accura te a nd not an artifac t  of using toofewb a sis  fu n ctio ns . In order toinv es t i gate thi s,  we utilise  the Bay es i an In format ion Cr iterion ( BI C ) [@Sc hwar z:1 978], whi chis 
aryingN_distributions"}](SpeedDistributions-VaryingN.pdf){width="49.00000%"}

Varying the_number of_basis functions
-------------------------------------

We first investigate_how the_reconstructed_WIMP mass_$m_\textrm{rec}$_and uncertainty varies_with the number_of basis functions $N$._For now, we_fix_our choice of basis to shifted Legendre polynomials, as used in Paper 1:

$$P_k(v) =_L_k\left(2\frac{v}{v_\textrm{max}}_- 1\right)\,,$$

where_$L_k$_is_the Legendre polynomial of order_$k$, and $v_\textrm{max}$ is a_cut off_for the parametrization. We should choose $v_\textrm{max}$ to_ensure_that $f_1(v)$ is_negligible above the cut off. However, too high a_choice of $v_\textrm{max}$ will result in_$f_1(v)$ being close_to_zero_over a large range_of the parametrization, making fitting more_difficult. We use the value $v_\textrm{max}_= 1000 {\,\textrm{km s}^{-1}}$, while lies significantly_above the Galactic escape speed.

The lower_panel of Fig. \[fig:BUMP\_LEG\] shows the_best fit_mass and 68% confidence intervals_as a function_of $N$,_using as input_a WIMP of mass 50 GeV_and the ‘bump’_distribution function. The reconstructed mass very_rapidly_settles close to_the_true_value, using_as few as_three_basis functions._This_is because adding the bump near_$v_\sim 500 {\,\textrm{km s}^{-1}}$ still leaves the_mean inverse speed relatively_smooth,_so a large number_of basis function are not_required. The correct mass is reconstructed_and we_emphasize in_the lower panel of Fig. \[fig:BUMP\_LEG\] that the reconstruction is stable with_the addition of more basis functions.

We_should also consider how_the quality_of_the fit changes_as_a function_of $N$. We would expect that adding_fit parameters_should always lead to a better_fit. Eventually, the fit_should_be good enough that adding additional_basis functions will no longer improve_it significantly. We can then_be_confident_that our reconstruction is accurate_and not an artifact of using_too few basis_functions. In order to investigate this, we_utilise_the Bayesian Information Criterion (BIC) [@Schwarz:1978],_which_is