27 double F[
nWF] = {1590,2090,2590};
36 sprintf(wf_name,
"RDE_%d_%1.1f",F[
n],TAU[n]);
38 wf.
name.ReplaceAll(
".",
"d");
40 for(
int m=0;
m<20;
m++) {
41 double cos_iota=gRandom->Uniform(-1,1);
42 double iota=acos(cos_iota)*180./
TMath::Pi();
47 wf.
par[0].name=
"freq"; wf.
par[0].value=F[
n];
48 wf.
par[1].name=
"tau"; wf.
par[1].value=TAU[
n];
49 wf.
par[2].name=
"iota"; wf.
par[2].value=iota;
65 vector<mdcpar>
par; par.resize(3);
66 par[0].name=
"entries";par[0].value=100000;
67 par[1].name=
"rho_min";par[1].value=1;
68 par[2].name=
"rho_max";par[2].value=100;
void SetInjRate(double inj_rate=MDC_INJ_RATE)
void SetInjLength(double inj_length=MDC_INJ_LENGTH)
mdcid AddWaveform(MDC_TYPE mdc_type, vector< mdcpar > par, TString uname="")
void SetSkyDistribution(MDC_DISTRIBUTION sky_distribution, vector< mdcpar > par, int seed=0, bool add=false)
void SetInjJitter(double inj_jitter=MDC_INJ_JITTER)
sprintf(tfres,"(1/%g)x(%g) (sec)x(Hz)", 2 *df, df)
void SetInjHrss(double inj_hrss=MDC_INJ_HRSS)
wavearray< double > GetRD(double frequency, double tau, double iota, bool polarization=0)