Short vs. Long GammaRay Bursts: A Comprehensive Study of Energetics and Prompt GammaRay Correlations
Abstract
We present the results of a comprehensive study of the luminosity function, energetics, prompt gammaray correlations, and classification methodology of short–hard and long–soft GRBs (SGRBs & LGRBs), based on observational data in the largest catalog of GRBs available to this date: BATSE catalog of 2702 GRBs. We find that: 1. The leastbiased classification method of GRBs into short and long, solely based on prompt–emission properties, appears to be the ratio of the observed spectral peak energy to the observed duration () with the dividing line at . 2. Once data is carefully corrected for the effects of the detection threshold of gammaray instruments, the population distribution of SGRBs and LGRBs can be individually well described as multivariate lognormal distribution in the –dimensional space of the isotropic peak gammaray luminosity, total isotropic gammaray emission, the intrinsic spectral peak energy, and the intrinsic duration. 3. Relatively large fractions of SGRBs and LGRBs with moderatetolow spectral peak energies have been missed by BATSE detectors. 4. Relatively strong and highly significant intrinsic hardness–brightness and duration–brightness correlations likely exist in both populations of SGRBs and LGRBs, once data is corrected for selection effects. The strengths of these correlations are very similar in both populations, implying similar mechanisms at work in both GRB classes, leading to the emergence of these prompt gammaray correlations.
keywords:
GammaRays: Bursts – GammaRays: observations – Methods: statistical1 Introduction
The field of GammaRay Bursts (GRBs) has witnessed rapid growth over the past decades, in particular, following the launches of NASA Compton GammaRay Observatory (Meegan et al., 1992), Swift (Gehrels et al., 2004) & Fermi (Michelson, Atwood & Ritz, 2010) missions. Early hints to the existence of distinct populations of gammaray transients (e.g., Mazets et al., 1981; Norris et al., 1984), and at least two classes of shorthard (TypeI) & longsoft (TypeII) GRBs (Dezalay et al., 1992) have now been extensively corroborated and confirmed by the promptemission data from independent gammaray detector missions (e.g., Kouveliotou et al., 1993a; Gehrels, RamirezRuiz & Fox, 2009; Zhang et al., 2012) or followup observations of the afterglows or host galaxies (e.g., Zhang et al., 2009; Berger, 2011, 2014). Although the possibility of more than two classes of GRBs with distinct progenitors has been extensively discussed and considered (e.g., Horvath, 1998; Mukherjee et al., 1998; Hakkila et al., 2001; Balastegui, RuizLapuente & Canal, 2001; Hakkila et al., 2004b; Horvath et al., 2006; Gehrels et al., 2006; Chattopadhyay et al., 2007; Horvath et al., 2008; Virgili, Liang & Zhang, 2009; Gao, Lu & Zhang, 2010; Horvath et al., 2012; Levan et al., 2014; KÃ³bori et al., 2014), it has remained a matter of debate and speculation to this date (e.g., Hakkila et al., 2000c, a, b, 2003, 2004a; Shahmoradi, 2013b; Zhang et al., 2014; Levan et al., 2014).
Beginning with the influential work of Kouveliotou et al. (1993a), GRBs have been traditionally classified into two populations of Short and Long GRBs (SGRBs & LGRBs respectively) based on a sharp cutoff on the bimodal distribution of the observed duration () of prompt gammaray emission, generally set to . Nevertheless, the detector & energy dependence of the observed GRB duration (e.g., Fenimore et al., 1995; Nemiroff, 2000; Qin et al., 2013) has prompted many studies in search of alternative lessbiased methods of GRB classification, typically based on a combination of the prompt gammaray emission andor afterglowhost properties (e.g., Gehrels, RamirezRuiz & Fox, 2009; Zhang et al., 2009; Goldstein, Preece & Briggs, 2010; LÃ¼ et al., 2010; Shahmoradi & Nemiroff, 2011; Zhang, Chen & Huang, 2012; Zhang et al., 2012; Shahmoradi, 2013b; LÃ¼ et al., 2014) or using the promptemission spectral correlations in conjunction with the traditional method of classification (e.g., Qin & Chen, 2013).
Ideally, a phenomenological classification method of GRBs should be based on their intrinsic (i.e., restframe) properties, free from potential biases due to data analysis, detector specifications, observational selection effects and sample incompleteness. Such method is currently far from reach due to detectorinduced heterogeneity in available GRB catalogs (e.g., Qin et al., 2013) and complex selection effects in the detection, analysis, and redshift measurement processes (e.g., Hakkila et al., 2000c, 2003; Band & Preece, 2005; Nakar & Piran, 2005; Nava et al., 2008; Butler et al., 2007; Butler, Kocevski & Bloom, 2009; Butler, Bloom & Poznanski, 2010; Shahmoradi & Nemiroff, 2009, 2011; Shahmoradi, 2013b, a; Coward et al., 2013, c.f., Shahmoradi & Nemiroff 2011 for a comprehensive review of relevant literature).
Historically, the intrinsic properties of short GRBs have also been less studied compared to long GRBs, potentially due to lack of redshift information for the majority of SGRBs (e.g., Coward et al., 2012). This has led to an additional obstacle towards a quantitative GRB taxonomy. While the population properties of long GRBs and their progenitors have been extensively researched and fairly well constrained (e.g., Dainotti, Ostrowski & Willingale, 2011; Shahmoradi, 2013b, a; Lien et al., 2013; Dainotti et al., 2013; Lien et al., 2014; Howell et al., 2014; Pescalli et al., 2015; Littlejohns & Butler, 2014; Dainotti et al., 2014, c.f., Shahmoradi 2013 for a comprehensive review of the literature), the cosmic rate and the intrinsic promptemission properties of short GRBs are far less understood and investigated (e.g., Guetta, Piran & Waxman, 2005; Guetta & Piran, 2006; Nakar, GalYam & Fox, 2006; Salvaterra et al., 2008; Chapman, Priddey & Tanvir, 2009; Virgili et al., 2011; Czerny et al., 2011; Wanderman & Piran, 2014).
Motivated by the existing gap and uncertainties in the current knowledge of the intrinsic population properties of short GRBs and the lack of an efficient, quantitative, biasfree classification method for gammaray bursts into long and short subgroups, here we present a methodology and model to constrain the energetics, luminosity function and the joint distributions and correlations of the prompt gammaray emission of SGRBs. Despite lacking a complete knowledge of the true cosmic rate and redshift distribution, here we argue and show that short GRBs exhibit very similar promptemission correlations and population distribution to those of long GRBs as presented by Shahmoradi (2013b), both qualitatively and, under plausible cosmic rate assumptions, also quantitatively. This result, along with other recent works on the timeresolved and timeintegrated spectral properties of SGRBs (e.g., Ghirlanda et al., 2011; Shahmoradi & Nemiroff, 2011; Zhang, Chen & Huang, 2012; Tsutsui et al., 2013; Calderone et al., 2014) points towards the possibility of a unified mechanism responsible for the prompt gammaray emission of the two GRB classes, independently of the diverse progenitor candidates for the two GRB populations.
The presented work also paves the way towards a detectorindependent minimallybiased phenomenological classification method for GRBs solely based on the observed prompt gammaray data of individual events. Towards this, we focus our attention on the largest sample of uniformlydetected gammaray bursts to this date: the BATSE catalog of 2702 GRBs (Paciesas et al., 1999; Goldstein et al., 2013). In the following sections, we present an example of data mining on BATSE data that showcases the tremendous amount of useful, yet unexplored information buried in this seemingly archaic GRB catalog.
We devote Sec. 2.1 of this manuscript to presenting an elaborate method of classification for the observed sample of GammaRay Bursts into two subclasses of shorthard and longsoft GRBs. In Sec. 2.2, we elaborate on the construction of a GRB world model that is capable of describing BATSE SGRB data. This is followed by a discussion of the procedure for fitting the model to BATSE data in Sec. 2.4. The predictions of the model together with univariate and multivariate goodnessoffit tests are presented in Sec. 3, followed by a discussion of the similarities and differences in the population distributions of SGRBs & LGRBs in Sec. 4. The main findings of the presented GRB world model are summarized in Sec. 5.
2 GRB World Model
2.1 Sample Selection
Depending on the triggering criteria, some gammaray detectors facilitate the detection of one class of bursts over the others (e.g., Paciesas et al., 1999; Hakkila, 2003; Band, 2003, 2006; Band et al., 2008; Lien et al., 2014). For example, the specific detector sensitivity of the Burst Alert Telescope (BAT) onboard Swift satellite results in better detections of LGRBs over SGRBs (e.g., Band, 2003, 2006). Compared to BAT however, BATSE Large Area Detectors (LAD) had an increased relative sensitivity to short GRBs (e.g., Hakkila, 2003; Band, 2003). Several studies have already offered new methods for a global detectorindependent phenomenological classification of GRBs based on empirical relations that are believed to hold for only a specific class of GRBs or based on the prompt or afterglow emission data (e.g., Zhang et al., 2009; LÃ¼ et al., 2010; Qin & Chen, 2013). These methods however, either suffer from sample incompleteness or require information, such as redshift, that might not be available for the majority of GRBs (e.g., Coward et al., 2012, 2013).
Here in this work, we ensure a minimallybiased analysis of shortduration class of GRBs by following the classification approach of Shahmoradi (2013b), which to the extent of our knowledge provides the least biased classification of BATSE GRBs, solely based on prompt emission data
The classification method used here is based on fuzzy Cmeans clustering algorithm of Dunn (1973); Bezdek (1981): Each BATSE GRB is assigned a probability (i.e., class coefficient) of belonging to LGRB (versus SGRB) population, the value of which depends on the set of GRB variables used in classification. This can be any combinations of the peak flux () in BATSE detection energy range , in three different time scales: , & ; bolometric fluence (); the observerframe spectral peak energy (), for which we use estimates by Shahmoradi & Nemiroff (2010), and the observed duration ( [s]).
To ensure minimal contamination of the two GRB classes, the lightcurves of bursts among BATSE GRBs with SGRB class coefficients in the range of were also visually inspected in the four main energy channels of BATSE Large Area Detectors. This led to reclassification of some events as Soft Gamma Repeaters (SGRs) or from one class to another, thus reducing the size of the original SGRB sample to (Table 1).
Trigger  Trigger  Trigger  Trigger  Trigger  Trigger  Trigger  Trigger  Trigger  Trigger  Trigger  Trigger  Trigger  Trigger 

108  138  185  207  218  229  254  289  297  373  432  474  480  486 
491  508  512  537  547  551  555  568  575  603  677  729  734  788 
799  809  830  834  836  845  856  867  878  906  909  929  936  942 
974  1051  1073  1076  1088  1096  1097  1102  1112  1128  1129  1154  1211  1223 
1289  1308  1346  1359  1404  1435  1443  1453  1461  1481  1518  1546  1553  1566 
1588  1634  1635  1636  1637  1659  1662  1665  1679  1680  1683  1694  1719  1736 
1741  1747  1760  1851  1953  1968  2003  2037  2040  2041  2043  2044  2049  2056 
2068  2099  2103  2115  2117  2125  2126  2132  2142  2145  2146  2155  2159  2161 
2163  2167  2201  2205  2206  2217  2220  2265  2268  2273  2283  2288  2312  2326 
2327  2330  2332  2352  2353  2357  2358  2360  2365  2368  2372  2377  2382  2384 
2395  2401  2424  2434  2448  2449  2454  2485  2487  2502  2504  2512  2513  2523 
2529  2536  2564  2583  2585  2597  2599  2614  2615  2623  2632  2633  2649  2679 
2680  2690  2693  2701  2715  2728  2748  2755  2757  2760  2776  2788  2795  2799 
2800  2801  2810  2814  2821  2823  2828  2834  2844  2846  2849  2851  2860  2861 
2873  2879  2892  2894  2910  2918  2933  2952  2964  2966  2973  2975  2977  2978 
2987  2988  2995  3016  3027  3037  3038  3039  3043  3051  3066  3073  3078  3087 
3094  3113  3114  3118  3121  3137  3144  3146  3152  3155  3160  3164  3173  3215 
3218  3266  3278  3280  3282  3286  3293  3294  3297  3308  3323  3333  3335  3338 
3340  3342  3349  3359  3374  3379  3384  3437  3441  3476  3477  3487  3494  3502 
3510  3530  3545  3606  3611  3640  3642  3665  3668  3722  3728  3735  3737  3742 
3751  3770  3774  3782  3791  3799  3810  3866  3867  3868  3888  3889  3894  3895 
3902  3904  3910  3919  3921  3936  3939  3940  4327  4660  4744  4776  4807  4871 
4955  5079  5206  5212  5277  5339  5439  5448  5453  5456  5458  5459  5461  5467 
5469  5471  5485  5488  5491  5498  5499  5500  5501  5527  5528  5529  5533  5536 
5537  5546  5547  5556  5560  5562  5564  5576  5592  5599  5607  5619  5620  5633 
5638  5647  5650  5664  5724  5730  5733  5740  5770  5992  6091  6096  6105  6117 
6120  6123  6135  6136  6145  6153  6166  6178  6180  6182  6204  6205  6215  6216 
6219  6230  6237  6238  6251  6263  6265  6275  6281  6284  6292  6299  6301  6307 
6314  6331  6338  6341  6342  6343  6347  6354  6361  6368  6372  6376  6385  6386 
6398  6401  6411  6412  6427  6436  6439  6443  6445  6447  6452  6462  6469  6486 
6488  6497  6535  6540  6542  6543  6547  6562  6569  6571  6573  6579  6580  6586 
6591  6606  6634  6635  6638  6641  6643  6645  6659  6662  6671  6679  6682  6689 
6693  6697  6700  6710  6715  6718  6753  6757  6786  6787  6788  6800  6824  6866 
6867  6870  6904  6916  6931  7009  7060  7063  7078  7102  7106  7133  7142  7148 
7159  7173  7187  7227  7240  7281  7283  7287  7290  7292  7294  7297  7305  7329 
7344  7353  7359  7361  7366  7367  7375  7378  7427  7430  7440  7447  7449  7453 
7455  7456  7472  7495  7496  7508  7514  7526  7547  7554  7559  7581  7584  7595 
7599  7601  7602  7626  7663  7671  7706  7710  7734  7745  7753  7754  7775  7784 
7789  7793  7800  7805  7827  7830  7901  7912  7922  7939  7943  7952  7970  7979 
7980  7988  7995  7999  8018  8027  8035  8041  8047  8072  8076  8077  8079  8082 
8085  8089  8097  8104  8120  —  —  —  —  —  —  —  —  — 
2.2 Model Construction
Our primary goal in this work is to derive a multivariate statistical model that, subject to BATSE detection threshold, is capable of reproducing the observational data of BATSE SGRBs. Examples of multivariate treatment of GRB luminosity function and energetics are rare in studies of GammaRay Bursts. Conversely, many authors have focused primarily on the univariate distribution of the spectral parameters, most importantly on the luminosity function. This is particularly true for the short population of GRBs where the prompt, afterglow and redshift information of individual events are scarcely available. In Shahmoradi (2013b) we argued that an accurate modelling of the luminosity function of LGRBs requires at least two GRB observable incorporated in the model: the bolometric peak flux () and the observed peak energy (). The parameter is required, since most gammaray detectors are photon counters, a quantity that depends on not only but also of the burst. This leads to the requirement of using a bivariate distribution as the minimum acceptable model for long GRBs, for the purpose of constraining the luminosity function.
For the class of short GRBs, the duration distribution (e.g., ) of the population spans a wide range from milliseconds to tens of seconds. The wide duration distribution is particularly important in modeling BATSE Large Area Detectors, since short GRBs could be potentially triggered on any of three triggering timescales: ms, ms & ms. Therefore, a SGRB world model should minimally incorporate the joint trivariate distribution of , and an appropriate definition of the observed duration (e.g., ). The duration variable is required in order to correctly account for the detection threshold of BATSE LADs. In addition, the definition of the observerframe parameter bolometric peak flux () and the corresponding restframe parameter merits special attention in the study of SGRBs. Here we use the definition of peak flux, , for SGRBs taken from BATSE catalog data, in contrast to the two other common definitions: & . Unlike the case for LGRBs (e.g., Shahmoradi, 2013b), we show in the addendum to this article (Appendix A & B) that is the least biased measure of peak flux for the majority of BATSE SGRBs, also the best definition for an appropriate modelling of the BATSE detection threshold in case of short GRBs.
Hereafter in the text and figures, the two parameters and implicitly refer to a definition of peak flux and luminosity wherever used for SGRBs and to the commonly used definitions wherever used for LGRBs.
Following the arguments of Shahmoradi (2013b) for LGRBs, we propose the multivariate lognormal distribution as the simplest natural candidate model capable of describing BATSE SGRB data. The motivation behind this choice of model comes from the available observational data that closely resembles a joint multivariate lognormal distribution for four most widely studied temporal and spectral parameters of both GRB classes in the observerframe: , (bolometric fluence), , : Since most SGRBs are expected to originate from low redshifts , the convolution of these observerframe parameters with the redshift distribution results in negligible variation in the shape of the restframe joint distribution of the same SGRB parameters. Therefore, the redshiftconvoluted 4Dimensional (4D) restframe distribution can be well approximated as a linear translation of the observerframe parameters to the restframe parameter space, keeping the shape of the distribution almost intact (e.g., Balazs et al., 2003). This implies that the joint distribution of the intrinsic SGRB variables: the isotropic peak luminosity (), the total isotropic emission (), the restframe timeintegrated spectral peak energy () and the restframe duration () might be indeed well described as a multivariate lognormal distribution.
We model the process of SGRB observation as a nonhomogeneous Poisson process whose mean rate parameter – the cosmic SGRB differential rate, – is the product of the differential comoving SGRB rate density with a dimensional lognormal Probability Density Function (pdf), , of four SGRB variables: , , and , with location vector and the scale (i.e., covariance) matrix ,
(1)  
where is the comoving SGRB rate density, and the factor in the denominator accounts for cosmological time dilation. The comoving volume element per unit redshift, , is given by,
(2) 
with standing for the luminosity distance,
(3) 
assuming a flat CDM cosmology, with parameters set to , and (Jarosik et al., 2011) for consistency with the work of Shahmoradi (2013b). The parameters & stand for the speed of light and the Hubble constant respectively.
The 4dimensional lognormal distribution of Equation (1), , has an intimate connection to the multivariate Gaussian distribution in the logarithmic space of SGRB observable parameters (c.f., Appendix D in Shahmoradi, 2013b).
Finally, in order to obtain the observed rate () of SGRBs detected by BATSE LADs, the cosmic SGRB rate, , in Equation (1) must be convolved with an accurate model of BATSE trigger efficiency for short GRBs, ,
(4) 
In reality, the variable is a highly complex function of observational conditions and prompt emission characteristics, almost unique to each individual GRB. Nevertheless, we show in Appendix A that it can be approximated as a generic function of the burst’s redshift (), isotropic luminosity (), the restframe spectral peak energy () and the restframe duration ().
2.3 The SGRB Rate Density
The largest source of uncertainty in population studies of short GRBs originates from the lack of an accurate knowledge of their cosmic rate. Only a small fraction of heterogeneouslydetected SGRBs have measured redshifts to this date (e.g., Coward et al., 2012) and redshift completeness often limits studies to the brightest events (e.g., D’Avanzo et al., 2014). It is therefore, perceivable that the current observed redshift distribution of SGRBs is likely strongly biased and not representative of the entire population of SGRBs (e.g., Nakar, 2007; Coward et al., 2013). An alternative approach to empirical determination of the rate of SGRBs is through population synthesis simulations (e.g., Belczynski et al., 2006), based on the assumption of ‘compact binary mergers’ as the progenitor of the majority of SGRBs (e.g., Narayan, Paczynski & Piran, 1992; Eichler et al., 1989; Paczynski, 1986). In this scenario, the cosmic rate of SGRBs follows the Star Formation Rate (SFR) convolved with a distribution of the delay time between the formation of a binary system and its coalescence due to gravitational radiation.
There is currently no consensus on the statistical moments and shape of the delay time distribution, solely based on observations of individual events and their host galaxies. The median delays vary widely in the range of billion years depending on the assumptions involved in estimation methods or in the dominant binary formation channels considered (e.g., Belczynski et al., 2006; Hopman et al., 2006; Bogomazov, Lipunov & Tutukov, 2007; Berger et al., 2007; Zheng & RamirezRuiz, 2007; Berger, 2011; Hao & Yuan, 2013; Guelbenzu et al., 2014). Recent results from population synthesis simulations however, favor very short delay times of a few hundred million years with a long negligible tail towards several billion years (e.g., O’Shaughnessy, Belczynski & Kalogera, 2008; Belczynski et al., 2010).
The computational expenses and limitations imposed on this work strongly limit the number of possible scenarios that could be considered for the cosmic rate of short GRBs. Thus, in order to approximate the comoving rate density of SGRBs, we adopt the Star Formation Rate (SFR) of Hopkins & Beacom (2006) in the form of a piecewise powerlaw function,
(5) 
with parameters set to the bestfit values of an updated SFR fit by Li (2008). The SFR is then convolved with a lognormal model of the delay time distribution (e.g., Nakar, 2007),
(6) 
with parameters in units of billion years (Gyrs) estimated from the population synthesis simulation results of (e.g., Belczynski et al., 2010), such that the comoving rate density of SGRBs is calculated as,
(7) 
with the universe’s age at redshift given by,
(8) 
We also fit data with alternative cosmic rates of SGRBs, assuming that SGRBs follow SFR of Hopkins & Beacom (2006) or the convolution of SFR with a long merger delay time of lognormal form with parameters as suggested by Nakar (2007). For all redshift scenarios, we find that the resulting bestfit parameters are qualitatively the same, although some parameters may exhibit quantitative differences at significance level.
2.4 Model Fitting
Now, with a statistical model at hand for the observed rate of short GRBs (i.e., Equation 4), we proceed to obtain the best fit parameters of the model to BATSE short GRBs data. In principle, any model fitting must take into account the observational uncertainties and any prior knowledge from independent sources, which can be achieved via Bayesian multilevel methodology (e.g., Hobson et al., 2010). This is done by first constructing the likelihood function, taking into account the uncertainties in observational data (e.g., Loredo, 2004): Under the assumption of symmetric Gaussian uncertainties, as it is the case with BATSE catalog data, the full Poisson likelihood of data given the parameters of the SGRB world model in Eqns. 1 & 4 can be written as,
(9) 
in which is a factor that properly normalizes the cosmic rate of SGRBs () and the vector , standing for the th SGRB Observation in BATSE catalog, has the likelihood of having the true parameters in the restframe 4dimensional observation space () that can be described as a Gaussian probability density function with parameters obtained from BATSE catalog such that,
(10) 
In this sense, the term in Equation 2.4 acts as a Bayesian prior for . In the absence of this prior knowledge however, as it is the case with BATSE short GRBs, the Empirical Bayes approach provides an alternative solution, in which an ad hoc estimate of the model parameters based on the observed data – excluding uncertainties – serves as the prior for the same data – including uncertainties – at the second level of analysis (e.g., Hobson et al., 2010). The joint posterior of the unknown parameters of the model can be then written as,
(11) 
As for the choice of hyperprior, , we adopt the noninformative uniform prior for the mean vector . A variety of noninformative or weaklyinformative priors for the covariance matrix have been already proposed and considered in the literature, with Inverse Wishart familty of distributions among the most popular choices (c.f., John Barnard (2000) and references therein). Here, to avoid problems and complications associated with Inverse Wishart priors, we adopt a separation strategy (e.g., Browne, 2006) by decomposing the covariance matrix into a correlation matrix and a set of standard deviations. We then use uniform priors on the logtransformation of all standard deviations, also on the Fishertransformation (Fisher, 1915) of all correlation coefficients.
Due to the complex truncation imposed on SGRB data and the world model by BATSE detection threshold (c.f., Appendix A), maximization of the likelihood function of Eqn. 2.4 is analytically intractable. Calculation of the normalization factor by itself requires a multivariate integral over the 4dimensional space of SGRB variables at any given redshift. In addition, due to lack of redshift () information for BATSE SGRBs, the probability for observation of each SGRB given the model parameters must be marginalized over all possible redshifts. These numerical integrations make sampling from the posterior distribution of Eqn. 11 an extremely difficult task. Given the potential presence of unknown systematic biases in BATSE catalog data as discussed in Appendix B, (also by Koshut et al., 1996; Hakkila et al., 2000c, 2003; Shahmoradi, 2013b, Appendix C) and the high level of uncertainty in the redshift distribution of short GRBs, we take a bold but reasonable and highly simplifying step and drop data uncertainties in the calculation of the likelihood function (Eqn. 2.4) in order to bring the problem into the realm of current computational technologies. The joint posterior distribution of the model parameters is then obtained by iterative sampling using a variant of Markov Chain Monte Carlo (MCMC) techniques known as Adaptive MetropolisHastings (e.g., Haario, Saksman & Tamminen, 2001). To further the efficiency of MCMC sampling, we implement all algorithms in Fortran (Backus, 1978; Metcalf, Reid & Cohen, 2011) and approximate the numerical integration in the definition of the luminosity distance of Eqn. 3 by the analytical expressions of Wickramasinghe & Ukwatta (2010). This integration is encountered on the order of billion times during MCMC sampling from the posterior distribution (c.f., Appendix C in Shahmoradi (2013b) for further details of the MCMC sampling method).
3 Results & GoodnessofFit Tests
Parameter  SGRBs World Model  LGRBs World Model 
Redshift Parameters (Equation 5)  
Lognormal Merger Delay (Equation 6)  
0.1  –  
1.12  –  
Location Parameters  
Scale Parameters  
Correlation Coefficients  
BATSE Detection Efficiency (Eqn. A)  
Note.— The full Markov Chain sampling of the above parameters from the 16dimensional parameter space of the likelihood function are available for download at https://sites.google.com/site/amshportal/research/aca/inthenews/lgrbworldmodel for the LGRB world model and at https://bitbucket.org/AmirShahmoradi/grbworldmodel for SGRBs world model.
The resulting mean and standard deviations of the model parameters are tabulated in Table 2. It is notable that the parameters of the model exhibit strong covariance as illustrated in the upper triangle of the correlation matrix of Table 3. For comparison, the lower triangle of the table represents the correlation matrix of the same parameters for the LGRB world model based on BATSE LGRBs (Table in Shahmoradi (2013b)). All location parameters of the model appear to strongly correlate with each other, so do the scale parameters. The correlations among the four variables of the world models also weaken with increasing the location parameters. Therefore, given BATSE data, a higher cosmic rates of LGRBs and SGRBs at high redshifts generally implies weaker correlations among the prompt gammaray emission variables, in particular, the brightnesshardness type of relations. A comparison of the predicted redshift distribution of BATSE SGRBs with LGRBs is illustrated in the left plot of Figure 2.
The lack of redshift information combined with the relatively small sample size of BATSE GRBs strongly limit the variety of GoodnessofFit (GoF) tests that can be employed to assess the performance of the model. In addition to small sample size, marginialization of the likelihood function (Eqn. 2.4) over redshift variable, and the systematic biases in data close to detection threshold (c.f., Appendix B), have major contributions to the uncertainties of the bestfit parameters of the model (Table 2) and increase the risk of overfitting. Following Shahmoradi (2013b), the fitting results can be first visually inspected by superposing the univariate distribution of four variables of the SGRB world model on BATSE data, as illustrated in Figure 3. For comparison, the results of LGRB world model fit to BATSE LGRB data from Shahmoradi (2013b) are also shown by the redcolor lines and histograms.
At the second level of GoF tests, the joint bivariate distributions of pairs of GRB variables can be compared to model predictions as illustrated in Figures 4, 5, 6. The trace of BATSE trigger threshold on the populations of both SGRBs and LGRBs is visible in most bivariate distribution plots. Following Shahmoradi (2013b) and in order to ensure a good fit to the joint bivariate distributions of the four prompt emission parameters of SGRBs, we also compare the model predictions with observational data along the principal axes of the bivariate distributions, as depicted in the center & bottom plots of Figures 4, 5, 6. Although, statistically not a sufficient test for multivariate goodnessoffit of the model to observational data, this comparison can provide strong evidence in favor of or against a good fit, at much higher confidence than solely comparing the model predictions and data via marginal distributions, as illustrated in Figure 3.
A comparison of the bivariate model predictions with data along the principal axes immediately reveals the potential systematic biases that exist in BATSE catalog data close to detection threshold (e.g., centerright & bottomleft plots of Figure 4). We show in Appendix B that this bias has its origins primarily in the duration () and peak photon flux () measurements of BATSE GRBs.
Parameter 

















1.00  0.95  0.59  0.77  0.87  0.78  0.28  0.42  0.17  0.04  0.16  0.04  0.25  0.27  0.49  0.26  
0.96  1.00  0.66  0.83  0.88  0.86  0.33  0.48  0.33  0.05  0.27  0.11  0.33  0.35  0.47  0.25  
0.88  0.90  1.00  0.43  0.62  0.70  0.58  0.23  0.43  0.67  0.10  0.71  0.09  0.52  0.28  0.16  
0.34  0.43  0.39  1.00  0.71  0.76  0.22  0.67  0.40  0.06  0.61  0.05  0.64  0.33  0.35  0.17  
0.90  0.90  0.81  0.35  1.00  0.90  0.36  0.40  0.27  0.06  0.19  0.18  0.25  0.39  0.32  0.15  
0.85  0.89  0.83  0.52  0.94  1.00  0.48  0.51  0.58  0.25  0.41  0.33  0.41  0.56  0.26  0.12  
0.56  0.60  0.78  0.32  0.57  0.66  1.00  0.10  0.36  0.60  0.05  0.64  0.04  0.27  0.17  0.11  
0.13  0.14  0.13  0.20  0.13  0.18  0.10  1.00  0.23  0.05  0.42  0.09  0.55  0.20  0.16  0.08  
0.00  0.10  0.17  0.52  0.01  0.31  0.31  0.02  1.00  0.47  0.58  0.37  0.32  0.41  0.00  0.01  
0.45  0.50  0.74  0.32  0.48  0.60  0.83  0.08  0.43  1.00  0.01  0.92  0.13  0.36  0.00  0.03  
0.42  0.35  0.31  0.56  0.45  0.22  0.14  0.01  0.59  0.05  1.00  0.04  0.89  0.46  0.04  0.06  
0.49  0.54  0.76  0.33  0.54  0.65  0.84  0.09  0.36  0.96  0.10  1.00  0.13  0.54  0.01  0.01  
0.41  0.35  0.33  0.55  0.44  0.26  0.19  0.11  0.38  0.14  0.95  0.16  1.00  0.54  0.02  0.02  
0.01  0.05  0.01  0.59  0.05  0.20  0.04  0.13  0.35  0.04  0.61  0.15  0.66  1.00  0.00  0.04  
0.66  0.64  0.59  0.19  0.49  0.44  0.34  0.08  0.06  0.24  0.26  0.26  0.24  0.06  1.00  0.74  
0.42  0.41  0.38  0.11  0.29  0.26  0.21  0.05  0.05  0.15  0.16  0.16  0.15  0.05  0.78  1.00 
This method of scanning the model and data along the principal axes of the joint bivariate distributions can be generalized to trivariate and quadruvariate joint distributions. For brevity, however, only the bivariate tests are presented here.
4 Discussion
Throughout previous sections, we presented and explained an elaborate analysis to constrain the energetics, luminosity function and the prompt gammaray correlations of short GRBs, subject to detection threshold of gammaray detectors, here BATSE LADs. The methodology employed in this work is similar to the LGRB population study of Shahmoradi (2013b). We have argued and shown that the intrinsic population distributions of both long and short GRBs can be well described by a multivariate lognormal model in the 4dimensional prompt gammaray emission parameter space of the isotropic peak luminosity (), total isotropic emission (), the spectral peak energy (), and the prompt duration (). This was done by first employing a fuzzy clustering method to segregate the BATSE catalog GRBs into long and short classes (2.1), followed by a careful modelling of the effects of the BATSE detection threshold on the two GRB populations (Appendix A).
Ideally, if the population properties of both SGRBs and LGRBs can be well described by multivariate lognormal distribution as argued in Sec. 2.2, also by Shahmoradi (2013a) & Shahmoradi (2013b), then a multivariate lognormal mixture model ought to be used in order to simultaneously fit for the distributions of the entire BATSE catalog of LGRBs and SGRBs together. Therefore, the classification of the observed sample prior to model fitting as performed in Sec. 2.1 would be unnecessary and the members of the two GRB populations would be automatically determined by the best fit parameters of the mixture model. Nevertheless, our experimentation with mixture models generally led to either degenerate or very poor fitting results for BATSE GRB data. This is primarily due to the morphological differences in the lightcurves and spectra of short and long GRBs, which in turn result in challenging difficulties in modelling the BATSE triggering algorithm in a unified framework for both GRB classes in the multivariate mixture model.
In the following sections, we discuss the key findings of the presented work and compare the results with data and findings from other gammaray experiments and population studies of GRBs.
4.1 Luminosity Function & Isotropic Emission
Despite being an illdefined quantity, the isotropic peak luminosity () is among the most widely used and investigated parameters of the prompt emission of GRBs. The vague and conventional definition of stems from its dependence on the timescale used to define this quantity, commonly set to [sec] (). Although the definition of the peak luminosity is more or less independent of the duration (e.g., ) of the long class of GRBs, it is strongly durationdependent for the class of short GRBs (Figure 8). This is primarily due to the diverse range of prompt durations of SGRBs, with many of the detected bursts lasting only a fraction of a second. Alternative definitions of the peak luminosity have been already proposed in order to provide a universal definition of the peak luminosity, independent of the duration and type of the burst, such as the effective luminosity definition of Butler, Bloom & Poznanski (2010). Nevertheless, such global luminosity definitions can not be employed in the presented analysis due to the specific triggering algorithm of BATSE LADs, which is defined for distinct fixed timescales: , & . Therefore, in order to minimize the effects of GRB duration on the definition of peak luminosity, we have used the highest resolution timescale () for the definition of the isotropic peak luminosity of BATSE SGRBs. By contrast, the onesecond peak luminosity () is sufficiently accurate and almost independent of the prompt duration of almost all LGRBs, whether detected or undetected (c.f., Shahmoradi, 2013b).
Given BATSE SGRB data, we find a range of for the distribution of the peak luminosities of short GRBs. Depending on the burst duration from long to short, the peak luminosity of SGRBs can be on average to times larger than the conventional timescale definition of the peak luminosity, commonly used for LGRBs. Knowing that the most luminous bursts generally tend to be the longest, we obtain a conservative range of for short GRBs. For comparison, Shahmoradi (2013b) finds a range of for the population of Long GRBs.
It should be noted however, that the promptemission GRB quantities, in particular, the isotropic peak luminosity and the cosmic rates of GRBs are very sensitive to the detection threshold of the specific detector used to collect observational data, and to the systematic biases in measurements close to detection threshold (c.f., Appendix B). This has also been noted in an earlier study of the population properties of Swift LGRBs by Butler, Bloom & Poznanski (2010).
As for the population distribution of the total isotropic emission of short GRBs, we obtain a range of . For comparison, a range of was obtained for LGRBs by Shahmoradi (2013b). Therefore, it appears that LGRBs are on average more than one order of magnitude brighter than SGRBs. It is however, notable that the width of the distribution of is approximately the same in both populations of LGRBs and SGRBs. The same also holds for the width of the distribution of in the two GRB classes.
4.2 Spectral Peak Energy & Prompt Duration
The spectral peak energies of SGRBs have long been observed to be systematically higher than the typical peak energies of the class of LGRBs (e.g., Kouveliotou et al., 1993b). The extent and significance of the difference between the two populations however, has remained a matter of speculation due to sample incompleteness and unknown selection biases in observational data of both GRB classes. The methodology presented in this work enabled us for the first time to set stringent constraints on the potential underlying population distribution of the spectral peak energies of both GRB classes in both observer and rest frames.
For the distribution of the intrinsic spectral peak energies of short GRBs, the model predicts an approximate range of , with an intrinsic average peak energy of . This corresponds to an observerframe range of , with an observerframe average peak energy of . For comparison, Shahmoradi (2013b) finds approximate ranges of & for the intrinsic and observerframe spectral peak energy distributions of LGRBs, with population averages & respectively.
In contrast, the distributions of the observed spectral peak energies of BATSE catalog LGRBs & SGRBs as found by Shahmoradi & Nemiroff (2010) (Figure 13 therein) are and respectively, slightly larger than the inferred values for the underlying population of the two classes here in this work. Similarly, Nava et al. (2011) also find slightly larger average values for a sample of Fermi SGRBs and LGRBs. The discrepancy between the detected sample and the underlying population can be explained by sample incompleteness due to BATSE LAD & Fermi GBM detection thresholds. This is line with previous findings on the role of selection effects due to gammaray detectors in shaping the observed properties of the two GRB classes (e.g., Band & Preece, 2005; Shahmoradi & Nemiroff, 2009; Butler, Kocevski & Bloom, 2009; Shahmoradi & Nemiroff, 2011; Kocevski, 2012b).
As for the prompt duration of short GRBs, the model predicts a range of for the intrinsic duration distribution of SGRBs with a population average of . This corresponds to a dynamic range of for the duration distribution of short GRBs in the observer frame with a population average of . When compared to BATSE detected sample of SGRBs with an average , it is apparent that the majority of undetected SGRBs were likely among the shortest events in the population. This corroborates the early hints on the potential existence of very shortduration GRBs that could have been missed by BATSE Large Area Detectors (e.g., Figure (1) & (2) in Nemiroff et al., 1998). The minor excess predicted by the SGRB world model may also be partly attributed to the systematic biases and errors in BATSE data and the simplified model of BATSE trigger threshold for SGRBs as explained in Appendix A & B.
For comparison, the LGRB world model of Shahmoradi (2013b) predicts an approximate range of for the intrinsic duration distribution of LGRBs with a population average of , corresponding to an approximate dynamic range of in the observer frame with a population average of .
The presented LGRB and SGRB world models also predict that intrinsically shorter duration GRBs in both classes, although likely exist, have lower chances of detection by gammaray instruments such as BATSE LAD and Fermi GBM. These results corroborate the recent findings of Littlejohns et al. (2013) & Littlejohns & Butler (2014) based on a sample of Swift LGRBs and further extend it to the population of short GRBs as illustrated in right plot of Figure 2.
4.3 Temporal & Spectral Correlations
There is already an extensive body of literature on the potential existence of correlations among the prompt spectral and temporal parameters of LGRBs (c.f., Shahmoradi & Nemiroff, 2011, for a comprehensive review of the literature). Much of the efforts so far has mainly focused on brightnesshardness types of relations, such as the Amati relation () and the Yonetoku relation (). Currently, the general consensus is that these relations do exist with high significance (e.g., Butler, Bloom & Poznanski, 2010; Shahmoradi, 2013b), but with far less strength than the original findings of Amati et al. (2002) and Yonetoku et al. (2004). In addition to early hints on the potential existence of durationbrightness relations (e.g., Horvath et al., 2005), recently Butler, Bloom & Poznanski (2010) also found, through an elaborate and comprehensive analysis of Swift data, some tentative signatures of a significant positive correlation between the intrinsic duration and the total isotropic emission of LGRBs. Later, Shahmoradi (2013b) provided evidence in support of intrinsic durationbrightness correlations among LGRBs using independent methods and models applied to BATSE catalog data, and further showed that the brightnessduration relations (i.e., & ) are of comparable strength and significance to brightnesshardness relations (i.e., & ). This positive correlation is also evident in the results of the LGRB world model in Table 2. Unlike the work of Shahmoradi (2013b), here in this work we have applied an energy band correction of the form to the observed durations () of BATSE LGRBs (e.g., Gehrels et al., 2006; Butler, Bloom & Poznanski, 2010). This correction if applied, slightly relaxes the strength of the correlation between the intrinsic duration and brightness.
Contrary to LGRBs, correlations among the promptemission parameters of short GRBs have been far less investigated. Recently, Zhang, Chen & Huang (2012) and Tsutsui et al. (2013) used a small sample of SGRBs with measured redshift to argue for the potential existence of intrinsic hardnessbrightness correlations in the population of SGRBs similar to the class of LGRBs. Nevertheless, the strength and significance of these positive correlations in the underlying population of SGRBs could not be determined in their studies due to strong selection bias and sample incompleteness in observational data.
The results of the SGRB world model, as presented in Table 2, confirm the existence of intrinsic hardnessbrightness correlations among SGRBs. Specifically, the model predicts highly significant correlation strengths of & for & relations respectively. The two correlation strengths are very similar to the predictions of the LGRB world model for the same relations in the population of LGRBs ( & , respectively). This is illustrated in the plots of Figure 7.
The similarity of the two GRB classes in the joint distributions and correlations is not limited to only hardnessbrightness relations. Indeed, all four promptemission variables in both populations (, , , ) appear to be similarly related to each other in both populations, as shown in Table 2.
A potential correlation of the form,
(12) 
has also been derived and suggested by Putten (2008), in which stands for the beaminganglecorrected output energy from the burst. The presented analysis is consistent with the existence of such universal relation. The strength and significance of it however, cannot be determined solely based on BATSE observational data, as it requires a knowledge of redshift and the beaming angle of individual events. We also caution against the use of the hardnessbrightness correlations, such as the Amati relation, to infer the characteristics of the inner engines of the two GRB classes. Our predictions based on BATSE catalog GRBs is that the Amati relation in its current form – as presented by Putten et al. (2014) – is likely strongly affected by sample incompleteness. In addition, the significant overlap of the two GRB populations in hardnessbrightness plots, as illustrated in Figure 7, likely renders the Amati relation an ineffective tool for GRB classification.
5 Concluding Remarks
The primary goal of the presented analysis was to constrain the energetics, luminosity function, and the prompt gammaray correlations of shorthard class of GRBs, using the wealth of information that has remained untouched in the largest catalog of GRBs available to this date, the current BATSE GRB catalog. In the following lines we summarize the steps we have taken to constrain the population properties of shorthard and similarly, longsoft GRBs (SGRBs & LGRBs, respectively):

A sample of shorthard and longsoft bursts were first segregated and selected from the current BATSE catalog of 2130 GRBs, for which complete data were available, including the bolometric peak flux (), the bolometric fluence (), the observed spectral peak energy () and the observed duration (). The classification method is based on fuzzy clustering algorithms on the two promptemission variables & which are least affected by the detection threshold of gammaray detectors (Sec. 2.1). This methodology can be readily applied to other GRB catalogs, in particular Fermi GBM.

We propose that the intrinsic joint distribution of the four main promptemission parameters of SGRBs: the isotropic peak gammaray luminosity (), the total isotropic gammaray emission (), the intrinsic spectral peak energy () and the intrinsic duration () can be well described as a multivariate (4dimensional) lognormal distribution, once the observational data is corrected for effects of detection threshold and sample incompleteness (Sec. 2.2).

The bestfit parameters of the model are then found by maximizing the likelihood function of the model given BATSE SGRB data (Eqn. 2.4) subject to the Bayesian priors of Eqn. 11, with a SGRB rate density (Eqn. 7) that is the result of the convolution of Star Formation Rate (Eqn 5) with a lognormal binary merger delay distribution of Eqn. 6. The resulting bestfit parameters are summarized in Table 2.
We highlight, in the following lines, the main conclusions of the presented analysis and the important similarities and differences that we find in the prompt gammaray emission properties of the two classes of shorthard and longsoft GRBs.

Population Distribution.
The population distributions of LGRBs and SGRBs appears to be well described by two separate multivariate lognormal distributions in the 4dimensional parameter space of the isotropic peak gammaray luminosity (), total isotropic gammaray emission (), the intrinsic spectral peak energy () and the intrinsic duration (), once corrected for the effects of detector threshold and sample incompleteness. This is in line with previous findings of Shahmoradi (2013b) & Shahmoradi (2013a). 
GRB Classification.
According to the predictions of our GRB world, the most accurate and the quickest method of individual GRB classification – solely based on prompt emission properties – appears to be the observerframe ratio . We find that of all LGRBs have , and of all SGRBs have (c.f., Figure 6, center right).We caution against the use of other similar quantities, such as the ratio of the observed spectral peak energy to bolometric fluence as proposed by Goldstein, Preece & Briggs (2010). Although this ratio seems to be a good discriminator in the sample of detected GRBs, it is strongly affected by sample incompleteness and detector threshold effects. Figure 5 (center right) illustrates the effects of sample incompleteness on the observed distribution of this ratio.

Energetics & Luminosity Function.
The presented GRB world model predicts a range of for the isotropic peak luminosity of SGRBs. A translation of this range to peak luminosity using Eqn. 15 approximately corresponds to . This range is very close and similar to the predictions of the GRB world model for peak luminosity distribution of LGRBs with a range of .Also predicted by the model are the ranges of & for the total isotropic gammaray emission of SGRBs & LGRBs respectively. The two variables and are strongly correlated with each other in both GRB classes (c.f., Table 2).

Prompt Duration & Spectral Peak Energy.
The population distribution of the restframe spectral peak energies () of both SGRBs and LGRBs appears to be described well by lognormal distributions with population averages & , and ranges & respectively. In the observerframe, this corresponds approximately to average & with ranges & for the two SGRB and LGRB classes, respectively.While the underlying duration distribution of LGRBs () does not seem to be significantly affected by the detection threshold of BATSE Large Area Detectors, there is tentative evidence that very shortduration SGRBs had, in general, lower chances of detection by BATSE (Figure 3, bottom right). For the population of SGRBs, we find a range of intrinsic duration with a population average of . In contrast, for distribution of LGRBs we find a range of with population average (Sec. 4.2).

Prompt GammaRay Correlations.
All four prompt gammaray variables appear to be strongly and positively correlated with each other in both GRB classes, with the exception of the two variables & which tend to be weakly, yet positively, correlated with each other. The intrinsic hardness–brightness relations (e.g., the Amati & the Yonetoku relations) are confirmed but with much higher dispersions than originally reported for these relations (Figure 7). The presented GRB world model reveals startling similarities in the strengths of the corresponding hardness–brightness correlations in the two GRB classes. Specifically, the model predicts a Pearson’s correlation strength of for – relation and for – relation, similarly in both GRB classes.
In summary, we have presented a mathematical model with minimal free parameters that enables us, for the first time, to constrain the main characteristics of the prompt gammaray emission of short–hard and long–soft GRBs, jointly and simultaneously, while paying careful attention to selection biases and sample incompleteness due to gammaray detector thresholds. Our model predicts a high level of similarity in the joint population distribution of the promptemission properties of the two GRB classes, a finding that merits further investigation of the potential similarities in the prompt emission mechanisms of both GRB classes.
Appendix A BATSE Detection Threshold
An accurate modelling of the detection threshold of gammaray instruments is an integral part of any population study of GRBs. We have already argued in Sec. 2.2 that modeling the trigger efficiency of gammaray detectors solely based on a measure of peak photon/energy flux – as is generally done in most GRB population studies – can potentially lead to systematic biases in the derived quantities. Although, the detection efficiency of most gamma–ray detectors depends solely on the observed peak photon flux in a limited energy window, the quantity of interest that is most often modelled and studied is the bolometric peak flux (). This variable depends on the observed peak photon flux and the spectral peak energy () for the class of LGRBs (e.g., Shahmoradi, 2013b), also on the observed duration (e.g., ) of the burst for the class of SGRBs. The effect of GRB duration on the peak flux measurement is very well illustrated in the left plot of Figure 8, where we show that for GRBs with , the timescale used for the definition of the peak flux does indeed matter. This is particularly important in modelling the triggering algorithm of BATSE Large Area Detectors, when a short burst could be potentially detected on any of the three different peak flux timescales used in the triggering algorithm: , & . Therefore, we adopt the following approach to construct a minimallybiased model of BATSE trigger efficiency for the population study of shorthard bursts.
First, since only one definition (i.e., timescale) of the bolometric peak flux can be incorporated in the SGRB world model of Sec. 2.2, we use the least biased definition of peak flux for SGRBs – the timescale definition – in the GRB world model. Although, this definition is durationindependent for virtually all BATSE GRBs, it becomes an increasingly biased measure of the peak flux for very long duration GRBs () close to detection threshold. We then approximate the discrete timescale trigger efficiency of BATSE LADs with a sigmoidal function that increases monotonically with increasing duration of the burst, from to . In other words, we convert the peak flux used in our GRB world model to an effective triggering peak flux , for which the detection efficiency of BATSE becomes durationindependent.
To expand on this, consider an idealized GRB lightcurve containing only a single squareshaped pulse with an exact duration of and a signal strength that is times the required significance for its detection on a peak flux timescale. In contrast, if there were only one triggering timescale available on BATSE, the signal strength of this event would fall right on the detection threshold of BATSE LADs. Thus, a burst of peak flux would be equivalent to an effective peak flux,
(13) 
for the triggering algorithm of BATSE on a timescale.
In reality however, GRB lightcurves are far more diverse than a single square pulse. Thus in order to build a more realistic model of BATSE LAD triggering algorithm, we fit a complementary Error function of the mathematical form,
(14) 
to the logarithm of the ratio of to peak fluxes () as a function of the observed duration () of BATSE GRBs, as illustrated in the left plot of Figure 8. The resulting bestfit function for has the form,
(15)  
The effective triggering peak flux in the SGRB world model is then calculated using the following relation,
(16) 
Once is obtained, we follow the approach of Shahmoradi (2013b) to calculate the detection probability () of a given SGRB with an effective triggering peak flux ,
(17) 
where & are the detection threshold parameters that are found by fitting the SGRB world model to BATSE observational data (c.f., Table 2), and is the effective triggering peak flux in BATSE energy range of detection , –[keV], calculated from the peak flux () in BATSE detection energy range using Eqn. 16. The connection between the restframe GRB parameters, , and the peak flux is obtained by fitting a smoothly broken powerlaw known as the Band model (Band et al., 1993) of the mathematical form,
(18) 
to SGRBs differential photon spectra, such that,
(19) 
where is the luminosity distance of Eqn. 3. In order to bring the above calculations into the realm of current computational technologies, we simplify the integration limits in the denominator of Eqn. 19 to a redshiftindependent energy range and fix the low– & high– energy photon indices of the Band model (Eqn. 18) to their corresponding population averages & . Butler, Bloom & Poznanski (2010) show that these simplifications result in an uncertainty of dex in the estimated peak flux, which is negligible compared to the existing systematic biases in BATSE data (c.f., Appendix B and uncertainties in the spectral peak energy estimates of Shahmoradi & Nemiroff (2010) used in this work). The resulting bestfit model of BATSE detection efficiency as a function of for the class of shorthard bursts is illustrated and compared to the detection efficiency of longsoft bursts in the right plot of Figure 8.
Appendix B Systematic Biases in BATSE GRB Data
As argued by Hogg & Turner (1998), astronomical catalogs and surveys are prone to systematic biases in measurements, in particular, close to the detection threshold of the observational instruments. The BATSE catalog of GRBs is no exception to such biases, also noted by BATSE team and others (e.g., Nemiroff et al., 1994; Paciesas et al., 1999; Stern et al., 2001; Stern, Atteia & Hurley, 2002). Although, throughout this work we relied on BATSE catalog data in their original form, here we present the results of our search for the potential signatures of systematic biases in BATSE data, which will pave the way for more accurate and rigorous population studies of GRBs in future.
In order to identify the extent of systematic bias in BATSE GRB data close to detection threshold, we first calculate the average background photon counts for each individual BATSE GRB lightcurve, in each of the four main energy channels of BATSE Large Area Detectors. we then subtract the calculated average background gammaray photon counts from the corresponding BATSE GRB lightcurves. The backgroundsubtracted lightcurves are then used to calculate the peak photon fluxes of all GRBs in the sample, in three , , & timescale definitions. Although the calculated peak fluxes are already contaminated and biased by background noise, we assume they represent the ‘true peak fluxes’ of BATSE GRBs and use them as our reference to simulate and investigate the effects of background noise in the calculation of peak flux at very low signaltonoise ratios. To do so, we add synthetic background noise to the entire timebins of each of the backgroundsubtracted lightcurves. The noise count () for each timebin in a given GRB lightcurve is drawn from the Poisson distribution,
(20) 
The mean of the noise count () for each BATSE GRB lightcurve is set to the original mean background photon counts found in each of the original lightcurves. We then subtract the average background counts () from the newly obtained lightcurves and measure the new peak fluxes in three timescales , & . This procedure of adding and subtracting synthetic background noise is then repeated time for each GRB lightcurve to obtain a sample of peak flux measurements for each of BATSE catalog GRB. The mean of this sample is then compared to the original ‘true peak flux’ measured after the initial background subtraction. The results are illustrated in the left plot of Figure 9. It is evident that peak the flux estimates become increasingly biased with decreasing SignaltoNoise Ratio (SNR) in BATSE data. This result corroborates the findings of Hogg & Turner (1998), who argue that flux measurements in astronomical surveys tend to be overestimated at very low SNR.
We have also tested the BATSE fluence and spectral peak energy data for the potential existence of systematic biases at low SNR, for which we find no significant evidence. We identify, however, a systematic bias in the duration measurements of BATSE catalog GRBs, specifically, in very long measurements. This is illustrated in the right plot of Figure 9, where we show that the measures of BATSE LGRBs are likely systematically underestimated at very long durations, corresponding to durations . As argued by Koshut et al. (1996); Hakkila et al. (2000c); Kocevski (2012a), the definition of GRB duration seems to be prone to systematic underestimations for very long duration lowSNR GRBs, since it is more likely that the latetime weak signals in the lightcurves would vanish in the background noise. The effects of this duration bias are also evident in the predictions of the LGRB world model shown in Figures 5 & 6 (top–left & center–left plots in both figures).
Footnotes
 footnotemark:
 pagerange: Short vs. Long GammaRay Bursts: A Comprehensive Study of Energetics and Prompt GammaRay Correlations–References
 pubyear: 2014
 The authors have already investigated multiple classification methodologies, on a variety of GRB characteristics, that are documented in Shahmoradi & Nemiroff (2011) & Shahmoradi (2013b).
 All BATSE GRB data are taken from the current BATSE Catalog and the spectral peak energies () are taken from Shahmoradi & Nemiroff (