LALInference  4.1.6.1-b72065a
lalinference_evolve_spins_and_append_samples Namespace Reference

Functions

def tilts_and_phi12_from_Cartesian_spins_and_L (chi1_v, chi2_v, Ln_v)
 

Variables

float dt = 0.1
 
 approx = lalsim.GetApproximantFromString("SpinTaylorT5")
 
 parser = argparse.ArgumentParser(description = 'Evolve the LALInference posterior samples of spins and append to HDF5 file of samples')
 
 help
 
 required
 
 type
 
 str
 
 default
 
 args = parser.parse_args()
 
 hdf_pos_file = args.sample_file
 
int v_final = 6.**-0.5
 
string label = '_isco'
 
 data = read_samples(hdf_pos_file)
 
dictionary metadata = {}
 
 run_identifier = extract_metadata(hdf_pos_file, metadata)
 
 q = np.atleast_1d(data['q'])
 
 eta = q/(1. + q)/(1. + q)
 
 m1
 
 m2
 
 chi1
 
 chi2
 
 tilt1
 
 tilt2
 
 phi12 = np.atleast_1d(data['phi12'])
 
 spinfreq_enum = np.array([lalsim.SimInspiralGetSpinFreqFromApproximant(int(lal_approx)) for lal_approx in data['LAL_APPROXIMANT']])
 
 f_start = np.where(np.array(spinfreq_enum == lalsim.SIM_INSPIRAL_SPINS_FLOW), data['flow'], data['f_ref'])
 
 tilt1_evol = np.zeros_like(m1)
 
 tilt2_evol = np.zeros_like(m1)
 
 phi12_evol = np.zeros_like(m1)
 
tuple mtot_s = (m1[i] + m2[i])*MTSUN_SI
 
int f_final = v_final*v_final*v_final/(mtot_s*np.pi)
 
 chi1x_v_data
 
 chi1y_v_data
 
 chi1z_v_data
 
 chi2x_v_data
 
 chi2y_v_data
 
 chi2z_v_data
 
 Lnx_v_data
 
 Lny_v_data
 
 Lnz_v_data
 
 deltaT
 
 MSUN_SI
 
 fStart
 
 fEnd
 
 s1x
 
 s1y
 
 s1z
 
 s2x
 
 s2y
 
 s2z
 
 lnhatx
 
 lnhaty
 
 lnhatz
 
 e1x
 
 e1y
 
 e1z
 
 lambda1
 
 lambda2
 
 quadparam1
 
 quadparam2
 
 spinO
 
 tideO
 
 phaseO
 
 lscorr
 
int idx_use = -1
 
 chi1_v = np.array([chi1x_v_data.data.data[idx_use], chi1y_v_data.data.data[idx_use], chi1z_v_data.data.data[idx_use]])
 
 chi2_v = np.array([chi2x_v_data.data.data[idx_use], chi2y_v_data.data.data[idx_use], chi2z_v_data.data.data[idx_use]])
 
 Ln_v = np.array([Lnx_v_data.data.data[idx_use], Lny_v_data.data.data[idx_use], Lnz_v_data.data.data[idx_use]])
 
 colnames = np.array(data.colnames)
 
 tilt1_id = np.where(colnames=='tilt1')[0][0]
 
 tilt2_id = np.where(colnames=='tilt2')[0][0]
 
 phi12_id = np.where(colnames=='phi12')[0][0]
 
 f_ref_id = np.where(colnames=='f_ref')[0][0]
 
 meta_tilt1 = tuple(data.columns.items())[tilt1_id][1].meta
 
 meta_tilt2 = tuple(data.columns.items())[tilt2_id][1].meta
 
 meta_phi12 = tuple(data.columns.items())[phi12_id][1].meta
 
 meta_f_ref = tuple(data.columns.items())[f_ref_id][1].meta
 
 tilts_evol = Table([Column(tilt1_evol, name='tilt1' + label, meta=meta_tilt1), Column(tilt2_evol, name='tilt2' + label, meta=meta_tilt2), Column(phi12_evol, name='phi12' + label, meta=meta_phi12)])
 
 data_joined = hstack([data, tilts_evol])
 
 vfinal_col = Table([Column(v_final*np.ones_like(tilt1_evol), name='vfinal', meta=meta_f_ref)])
 
 f = h5py.File(hdf_pos_file, 'r')
 
string path = '/lalinference/'+run_identifier+'/posterior_samples'
 
 level = f[path]
 
 arrt = level.attrs
 
 names = np.array([list(arrt.items())[i][0] for i in range(len(list(arrt.items())))])
 
int num_names = 0
 
 overwrite
 

Function Documentation

◆ tilts_and_phi12_from_Cartesian_spins_and_L()

def lalinference_evolve_spins_and_append_samples.tilts_and_phi12_from_Cartesian_spins_and_L (   chi1_v,
  chi2_v,
  Ln_v 
)

Variable Documentation

◆ dt

float lalinference_evolve_spins_and_append_samples.dt = 0.1

◆ approx

lalinference_evolve_spins_and_append_samples.approx = lalsim.GetApproximantFromString("SpinTaylorT5")

◆ parser

lalinference_evolve_spins_and_append_samples.parser = argparse.ArgumentParser(description = 'Evolve the LALInference posterior samples of spins and append to HDF5 file of samples')

◆ help

lalinference_evolve_spins_and_append_samples.help

◆ required

lalinference_evolve_spins_and_append_samples.required

◆ type

lalinference_evolve_spins_and_append_samples.type

◆ str

lalinference_evolve_spins_and_append_samples.str

◆ default

lalinference_evolve_spins_and_append_samples.default

◆ args

lalinference_evolve_spins_and_append_samples.args = parser.parse_args()

◆ hdf_pos_file

lalinference_evolve_spins_and_append_samples.hdf_pos_file = args.sample_file

◆ v_final

lalinference_evolve_spins_and_append_samples.v_final = 6.**-0.5

◆ label

string lalinference_evolve_spins_and_append_samples.label = '_isco'

◆ data

lalinference_evolve_spins_and_append_samples.data = read_samples(hdf_pos_file)

◆ metadata

lalinference_evolve_spins_and_append_samples.metadata = {}

◆ run_identifier

lalinference_evolve_spins_and_append_samples.run_identifier = extract_metadata(hdf_pos_file, metadata)

◆ q

lalinference_evolve_spins_and_append_samples.q = np.atleast_1d(data['q'])

◆ eta

lalinference_evolve_spins_and_append_samples.eta = q/(1. + q)/(1. + q)

◆ m1

lalinference_evolve_spins_and_append_samples.m1

◆ m2

lalinference_evolve_spins_and_append_samples.m2

◆ chi1

lalinference_evolve_spins_and_append_samples.chi1

◆ chi2

lalinference_evolve_spins_and_append_samples.chi2

◆ tilt1

lalinference_evolve_spins_and_append_samples.tilt1

◆ tilt2

lalinference_evolve_spins_and_append_samples.tilt2

◆ phi12

lalinference_evolve_spins_and_append_samples.phi12 = np.atleast_1d(data['phi12'])

◆ spinfreq_enum

lalinference_evolve_spins_and_append_samples.spinfreq_enum = np.array([lalsim.SimInspiralGetSpinFreqFromApproximant(int(lal_approx)) for lal_approx in data['LAL_APPROXIMANT']])

◆ f_start

lalinference_evolve_spins_and_append_samples.f_start = np.where(np.array(spinfreq_enum == lalsim.SIM_INSPIRAL_SPINS_FLOW), data['flow'], data['f_ref'])

◆ tilt1_evol

lalinference_evolve_spins_and_append_samples.tilt1_evol = np.zeros_like(m1)

◆ tilt2_evol

lalinference_evolve_spins_and_append_samples.tilt2_evol = np.zeros_like(m1)

◆ phi12_evol

lalinference_evolve_spins_and_append_samples.phi12_evol = np.zeros_like(m1)

◆ mtot_s

lalinference_evolve_spins_and_append_samples.mtot_s = (m1[i] + m2[i])*MTSUN_SI

◆ f_final

lalinference_evolve_spins_and_append_samples.f_final = v_final*v_final*v_final/(mtot_s*np.pi)

◆ chi1x_v_data

lalinference_evolve_spins_and_append_samples.chi1x_v_data

◆ chi1y_v_data

lalinference_evolve_spins_and_append_samples.chi1y_v_data

◆ chi1z_v_data

lalinference_evolve_spins_and_append_samples.chi1z_v_data

◆ chi2x_v_data

lalinference_evolve_spins_and_append_samples.chi2x_v_data

◆ chi2y_v_data

lalinference_evolve_spins_and_append_samples.chi2y_v_data

◆ chi2z_v_data

lalinference_evolve_spins_and_append_samples.chi2z_v_data

◆ Lnx_v_data

lalinference_evolve_spins_and_append_samples.Lnx_v_data

◆ Lny_v_data

lalinference_evolve_spins_and_append_samples.Lny_v_data

◆ Lnz_v_data

lalinference_evolve_spins_and_append_samples.Lnz_v_data

◆ deltaT

lalinference_evolve_spins_and_append_samples.deltaT

◆ MSUN_SI

lalinference_evolve_spins_and_append_samples.MSUN_SI

◆ fStart

lalinference_evolve_spins_and_append_samples.fStart

◆ fEnd

lalinference_evolve_spins_and_append_samples.fEnd

◆ s1x

lalinference_evolve_spins_and_append_samples.s1x

◆ s1y

lalinference_evolve_spins_and_append_samples.s1y

◆ s1z

lalinference_evolve_spins_and_append_samples.s1z

◆ s2x

lalinference_evolve_spins_and_append_samples.s2x

◆ s2y

lalinference_evolve_spins_and_append_samples.s2y

◆ s2z

lalinference_evolve_spins_and_append_samples.s2z

◆ lnhatx

lalinference_evolve_spins_and_append_samples.lnhatx

◆ lnhaty

lalinference_evolve_spins_and_append_samples.lnhaty

◆ lnhatz

lalinference_evolve_spins_and_append_samples.lnhatz

◆ e1x

lalinference_evolve_spins_and_append_samples.e1x

◆ e1y

lalinference_evolve_spins_and_append_samples.e1y

◆ e1z

lalinference_evolve_spins_and_append_samples.e1z

◆ lambda1

lalinference_evolve_spins_and_append_samples.lambda1

◆ lambda2

lalinference_evolve_spins_and_append_samples.lambda2

◆ quadparam1

lalinference_evolve_spins_and_append_samples.quadparam1

◆ quadparam2

lalinference_evolve_spins_and_append_samples.quadparam2

◆ spinO

lalinference_evolve_spins_and_append_samples.spinO

◆ tideO

lalinference_evolve_spins_and_append_samples.tideO

◆ phaseO

lalinference_evolve_spins_and_append_samples.phaseO

◆ lscorr

lalinference_evolve_spins_and_append_samples.lscorr

◆ idx_use

int lalinference_evolve_spins_and_append_samples.idx_use = -1

◆ chi1_v

lalinference_evolve_spins_and_append_samples.chi1_v = np.array([chi1x_v_data.data.data[idx_use], chi1y_v_data.data.data[idx_use], chi1z_v_data.data.data[idx_use]])

◆ chi2_v

lalinference_evolve_spins_and_append_samples.chi2_v = np.array([chi2x_v_data.data.data[idx_use], chi2y_v_data.data.data[idx_use], chi2z_v_data.data.data[idx_use]])

◆ Ln_v

lalinference_evolve_spins_and_append_samples.Ln_v = np.array([Lnx_v_data.data.data[idx_use], Lny_v_data.data.data[idx_use], Lnz_v_data.data.data[idx_use]])

◆ colnames

lalinference_evolve_spins_and_append_samples.colnames = np.array(data.colnames)

◆ tilt1_id

lalinference_evolve_spins_and_append_samples.tilt1_id = np.where(colnames=='tilt1')[0][0]

◆ tilt2_id

lalinference_evolve_spins_and_append_samples.tilt2_id = np.where(colnames=='tilt2')[0][0]

◆ phi12_id

lalinference_evolve_spins_and_append_samples.phi12_id = np.where(colnames=='phi12')[0][0]

◆ f_ref_id

lalinference_evolve_spins_and_append_samples.f_ref_id = np.where(colnames=='f_ref')[0][0]

◆ meta_tilt1

lalinference_evolve_spins_and_append_samples.meta_tilt1 = tuple(data.columns.items())[tilt1_id][1].meta

◆ meta_tilt2

lalinference_evolve_spins_and_append_samples.meta_tilt2 = tuple(data.columns.items())[tilt2_id][1].meta

◆ meta_phi12

lalinference_evolve_spins_and_append_samples.meta_phi12 = tuple(data.columns.items())[phi12_id][1].meta

◆ meta_f_ref

lalinference_evolve_spins_and_append_samples.meta_f_ref = tuple(data.columns.items())[f_ref_id][1].meta

◆ tilts_evol

lalinference_evolve_spins_and_append_samples.tilts_evol = Table([Column(tilt1_evol, name='tilt1' + label, meta=meta_tilt1), Column(tilt2_evol, name='tilt2' + label, meta=meta_tilt2), Column(phi12_evol, name='phi12' + label, meta=meta_phi12)])

◆ data_joined

lalinference_evolve_spins_and_append_samples.data_joined = hstack([data, tilts_evol])

◆ vfinal_col

lalinference_evolve_spins_and_append_samples.vfinal_col = Table([Column(v_final*np.ones_like(tilt1_evol), name='vfinal', meta=meta_f_ref)])

◆ f

lalinference_evolve_spins_and_append_samples.f = h5py.File(hdf_pos_file, 'r')

◆ path

lalinference_evolve_spins_and_append_samples.path = '/lalinference/'+run_identifier+'/posterior_samples'

◆ level

lalinference_evolve_spins_and_append_samples.level = f[path]

◆ arrt

lalinference_evolve_spins_and_append_samples.arrt = level.attrs

◆ names

lalinference_evolve_spins_and_append_samples.names = np.array([list(arrt.items())[i][0] for i in range(len(list(arrt.items())))])

◆ num_names

int lalinference_evolve_spins_and_append_samples.num_names = 0

◆ overwrite

lalinference_evolve_spins_and_append_samples.overwrite