Variables
	start_time = time.time()

	parser = ap.ArgumentParser()

	type

	help

	required

	default

	int

	nargs

	action

	float

	const

	str

	args = parser.parse_args()

	outbasedirectory = args.outbasedir

string	logging_format = "%(asctime)s : %(levelname)s : %(message)s"

	logging_filename = os.path.join(outbasedirectory, f"PiecewiseSearchLog_{j}.log")

	filename

	filemode

	level

	format

	loglevel

	LatticeTilingProgressLogLevel

	job_num = os.getenv("SLURM_ARRAY_TASK_ID")

	pr = cProfile.Profile()

	tbankcode = args.tbankcode

	j = args.j

	trefsegfrac = args.trefsegfrac

	rtnsum = args.rtnsum

	tempsperfile = args.tempsperfile

	SFTFiles = args.SFTFiles

	UseLocal = args.UseLocal

	BaseSeed = args.BaseSeed

	iterations = args.iterations

int	Random_Seed = BaseSeed * 10000000 + j * 1000

string	Random_Seed_Str = "{:0>9d}".format(Random_Seed)

string	mode = "search"

	tbank = cd.TBank()

	fmin

	fmax

	knotslist

	dur

	knots

	flags_bbox

	flags_intbox

	detectors = tbank.detectors

	num_det = len(detectors)

	detector_str_list = lal.CreateStringVector(detectors[0])

	noise_path = tbank.noise_path

	psd_path = tbank.psd_path

list	noise_curve_data = []

string	path_to_data = noise_path + detector + "asd.txt"

int	min_freq_strain = -1

int	max_freq_strain = -1

	this_line = line.split()

	frequency = float(this_line[0])

	strain = float(this_line[1])

tuple	this_det_noise = (min_freq_strain + max_freq_strain) / 2

	finputdata = cd.FInput()

	Tsft

	tstart

	h0

	cosi

	psi

	phi0

	dt_wf

	Alpha

	Delta

	noise_sqrt_Sh

	dfreq

	sourceDeltaT

	inject_data

string	append_string = "_BS-" + str(args.BaseSeed) + "_j-" + str(j)

	local_disk_directory = os.getenv("JOBFS")

tuple	tbankdirectory

	minimum_h0 = args.min_h0

	maximum_h0 = args.max_h0

	power_list = np.linspace(np.log10(minimum_h0), np.log10(maximum_h0), iterations)

list	h0_list = [0] + [10**i for i in power_list]

int	plus_one = 0

int	this_seed = Random_Seed + i

string	sumstr = ""

string	iteration_string = "_i-" + str(i)

string	tempsfilename = "Temps_For_" + tbank.toString() + iteration_string + sumstr + ".txt"

tuple	mismatchtempname

	seed = random.seed(this_seed)

	signalparams = lal.gsl_matrix(tbank.s * len(bf.knotslist), 1)

	SFTdirectory

	parentdirectory

	rand_seed

	tiling = lp.CreateLatticeTiling(tbank.s * len(bf.knotslist))

	metric = scmm.PreCompMetric(tbank.s)

	randparams = lal.CreateRandomParams(this_seed)

	locator = lp.CreateLatticeTilingLocator(tiling)

	signal_params_vec = lal.gsl_vector(tbank.s * len(bf.knotslist))

	nearest_point = lal.gsl_vector(tbank.s * len(bf.knotslist))

	nearest_index = lal.CreateUINT8Vector(tbank.s * len(bf.knotslist))

	f0 = signalparams[0]

	f1 = signalparams[1]

	f2 = signalparams[2]

string	commentline = "{:20s}".format("#FStat") + " "

int	knotnum = -1

	derivnum = i % tbank.s

string	columntitle = "f_" + str(knotnum) + "_" + str(derivnum)

int	SFTfpad = 10

int	SFTfmin = gom.gte(tbank.dur, tbank.fmin, tbank.nmax, tbank.kmax) - SFTfpad - 10

int	SFTfmax = gom.gte(0, tbank.fmax, tbank.nmin, tbank.kmin) + SFTfpad

	Fstat_mismatch = args.Fstat_mismatch

	out_directory

	fstat_hist

string	temp_count_directory = "Template_Count_And_Timing_For_" + tbank.toString()

string	temp_count_file = "Template_Count_And_Timing_" + Random_Seed_Str + ".txt"

	time_taken = time.time() - start_time

string	line = str(time_taken) + " " + str(template_count) + "\n"

string	Two_F_directory = "Largest_2Fs_For_" + tbank.toString()

string	Two_F_file = "Largest_2F_" + Random_Seed_Str + ".txt"

string	mismatch_directory = "Mismatches_For_" + tbank.toString()

string	mismatch_file = "Mismatch_" + Random_Seed_Str + ".txt"

	s = io.StringIO()

	ps = pstats.Stats(pr, stream=s).sort_stats("cumtime")

Variable Documentation

◆ start_time

lalpulsar_PiecewiseSearch.start_time = time.time()

Definition at line 48 of file lalpulsar_PiecewiseSearch.py.

◆ parser

lalpulsar_PiecewiseSearch.parser = ap.ArgumentParser()

Definition at line 53 of file lalpulsar_PiecewiseSearch.py.

◆ type

lalpulsar_PiecewiseSearch.type

Definition at line 62 of file lalpulsar_PiecewiseSearch.py.

◆ help

lalpulsar_PiecewiseSearch.help

Definition at line 63 of file lalpulsar_PiecewiseSearch.py.

◆ required

lalpulsar_PiecewiseSearch.required

Definition at line 64 of file lalpulsar_PiecewiseSearch.py.

◆ default

lalpulsar_PiecewiseSearch.default

Definition at line 72 of file lalpulsar_PiecewiseSearch.py.

◆ int

lalpulsar_PiecewiseSearch.int

Definition at line 76 of file lalpulsar_PiecewiseSearch.py.

◆ nargs

lalpulsar_PiecewiseSearch.nargs

Definition at line 87 of file lalpulsar_PiecewiseSearch.py.

◆ action

lalpulsar_PiecewiseSearch.action

Definition at line 99 of file lalpulsar_PiecewiseSearch.py.

◆ float

lalpulsar_PiecewiseSearch.float

Definition at line 122 of file lalpulsar_PiecewiseSearch.py.

◆ const

lalpulsar_PiecewiseSearch.const

Definition at line 159 of file lalpulsar_PiecewiseSearch.py.

◆ str

lalpulsar_PiecewiseSearch.str

Definition at line 363 of file lalpulsar_PiecewiseSearch.py.

◆ args

lalpulsar_PiecewiseSearch.args = parser.parse_args()

Definition at line 396 of file lalpulsar_PiecewiseSearch.py.

◆ outbasedirectory

lalpulsar_PiecewiseSearch.outbasedirectory = args.outbasedir

Definition at line 397 of file lalpulsar_PiecewiseSearch.py.

◆ logging_format

string lalpulsar_PiecewiseSearch.logging_format = "%(asctime)s : %(levelname)s : %(message)s"

Definition at line 402 of file lalpulsar_PiecewiseSearch.py.

◆ logging_filename

lalpulsar_PiecewiseSearch.logging_filename = os.path.join(outbasedirectory, f"PiecewiseSearchLog_{j}.log")

Definition at line 404 of file lalpulsar_PiecewiseSearch.py.

◆ filename

lalpulsar_PiecewiseSearch.filename

Definition at line 406 of file lalpulsar_PiecewiseSearch.py.

◆ filemode

lalpulsar_PiecewiseSearch.filemode

Definition at line 407 of file lalpulsar_PiecewiseSearch.py.

◆ level

lalpulsar_PiecewiseSearch.level

Definition at line 408 of file lalpulsar_PiecewiseSearch.py.

◆ format

lalpulsar_PiecewiseSearch.format

Definition at line 409 of file lalpulsar_PiecewiseSearch.py.

◆ loglevel

lalpulsar_PiecewiseSearch.loglevel

Definition at line 413 of file lalpulsar_PiecewiseSearch.py.

◆ LatticeTilingProgressLogLevel

lalpulsar_PiecewiseSearch.LatticeTilingProgressLogLevel

Definition at line 417 of file lalpulsar_PiecewiseSearch.py.

◆ job_num

lalpulsar_PiecewiseSearch.job_num = os.getenv("SLURM_ARRAY_TASK_ID")

Definition at line 419 of file lalpulsar_PiecewiseSearch.py.

◆ pr

lalpulsar_PiecewiseSearch.pr = cProfile.Profile()

Definition at line 425 of file lalpulsar_PiecewiseSearch.py.

◆ tbankcode

lalpulsar_PiecewiseSearch.tbankcode = args.tbankcode

Definition at line 434 of file lalpulsar_PiecewiseSearch.py.

◆ j

lalpulsar_PiecewiseSearch.j = args.j

Definition at line 440 of file lalpulsar_PiecewiseSearch.py.

◆ trefsegfrac

lalpulsar_PiecewiseSearch.trefsegfrac = args.trefsegfrac

Definition at line 441 of file lalpulsar_PiecewiseSearch.py.

◆ rtnsum

lalpulsar_PiecewiseSearch.rtnsum = args.rtnsum

Definition at line 442 of file lalpulsar_PiecewiseSearch.py.

◆ tempsperfile

lalpulsar_PiecewiseSearch.tempsperfile = args.tempsperfile

Definition at line 443 of file lalpulsar_PiecewiseSearch.py.

◆ SFTFiles

lalpulsar_PiecewiseSearch.SFTFiles = args.SFTFiles

Definition at line 444 of file lalpulsar_PiecewiseSearch.py.

◆ UseLocal

lalpulsar_PiecewiseSearch.UseLocal = args.UseLocal

Definition at line 445 of file lalpulsar_PiecewiseSearch.py.

◆ BaseSeed

lalpulsar_PiecewiseSearch.BaseSeed = args.BaseSeed

Definition at line 446 of file lalpulsar_PiecewiseSearch.py.

◆ iterations

lalpulsar_PiecewiseSearch.iterations = args.iterations

Definition at line 447 of file lalpulsar_PiecewiseSearch.py.

◆ Random_Seed

int lalpulsar_PiecewiseSearch.Random_Seed = BaseSeed * 10000000 + j * 1000

Definition at line 449 of file lalpulsar_PiecewiseSearch.py.

◆ Random_Seed_Str

string lalpulsar_PiecewiseSearch.Random_Seed_Str = "{:0>9d}".format(Random_Seed)

Definition at line 450 of file lalpulsar_PiecewiseSearch.py.

◆ mode

string lalpulsar_PiecewiseSearch.mode = "search"

Definition at line 455 of file lalpulsar_PiecewiseSearch.py.

◆ tbank

lalpulsar_PiecewiseSearch.tbank = cd.TBank()

Definition at line 473 of file lalpulsar_PiecewiseSearch.py.

◆ fmin

lalpulsar_PiecewiseSearch.fmin

Definition at line 495 of file lalpulsar_PiecewiseSearch.py.

◆ fmax

lalpulsar_PiecewiseSearch.fmax

Definition at line 496 of file lalpulsar_PiecewiseSearch.py.

◆ knotslist

lalpulsar_PiecewiseSearch.knotslist

Definition at line 506 of file lalpulsar_PiecewiseSearch.py.

◆ dur

lalpulsar_PiecewiseSearch.dur

Definition at line 509 of file lalpulsar_PiecewiseSearch.py.

◆ knots

lalpulsar_PiecewiseSearch.knots

Definition at line 515 of file lalpulsar_PiecewiseSearch.py.

◆ flags_bbox

lalpulsar_PiecewiseSearch.flags_bbox

Definition at line 529 of file lalpulsar_PiecewiseSearch.py.

◆ flags_intbox

lalpulsar_PiecewiseSearch.flags_intbox

Definition at line 532 of file lalpulsar_PiecewiseSearch.py.

◆ detectors

lalpulsar_PiecewiseSearch.detectors = tbank.detectors

Definition at line 538 of file lalpulsar_PiecewiseSearch.py.

◆ num_det

lalpulsar_PiecewiseSearch.num_det = len(detectors)

Definition at line 544 of file lalpulsar_PiecewiseSearch.py.

◆ detector_str_list

lalpulsar_PiecewiseSearch.detector_str_list = lal.CreateStringVector(detectors[0])

Definition at line 546 of file lalpulsar_PiecewiseSearch.py.

◆ noise_path

lalpulsar_PiecewiseSearch.noise_path = tbank.noise_path

Definition at line 552 of file lalpulsar_PiecewiseSearch.py.

◆ psd_path

lalpulsar_PiecewiseSearch.psd_path = tbank.psd_path

Definition at line 553 of file lalpulsar_PiecewiseSearch.py.

◆ noise_curve_data

list lalpulsar_PiecewiseSearch.noise_curve_data = []

Definition at line 555 of file lalpulsar_PiecewiseSearch.py.

◆ path_to_data

string lalpulsar_PiecewiseSearch.path_to_data = noise_path + detector + "asd.txt"

Definition at line 560 of file lalpulsar_PiecewiseSearch.py.

◆ min_freq_strain

lalpulsar_PiecewiseSearch.min_freq_strain = -1

Definition at line 565 of file lalpulsar_PiecewiseSearch.py.

◆ max_freq_strain

lalpulsar_PiecewiseSearch.max_freq_strain = -1

Definition at line 566 of file lalpulsar_PiecewiseSearch.py.

◆ this_line

lalpulsar_PiecewiseSearch.this_line = line.split()

Definition at line 574 of file lalpulsar_PiecewiseSearch.py.

◆ frequency

lalpulsar_PiecewiseSearch.frequency = float(this_line[0])

Definition at line 575 of file lalpulsar_PiecewiseSearch.py.

◆ strain

lalpulsar_PiecewiseSearch.strain = float(this_line[1])

Definition at line 576 of file lalpulsar_PiecewiseSearch.py.

◆ this_det_noise

tuple lalpulsar_PiecewiseSearch.this_det_noise = (min_freq_strain + max_freq_strain) / 2

Definition at line 585 of file lalpulsar_PiecewiseSearch.py.

◆ finputdata

lalpulsar_PiecewiseSearch.finputdata = cd.FInput()

Definition at line 597 of file lalpulsar_PiecewiseSearch.py.

◆ Tsft

lalpulsar_PiecewiseSearch.Tsft

Definition at line 598 of file lalpulsar_PiecewiseSearch.py.

◆ tstart

lalpulsar_PiecewiseSearch.tstart

Definition at line 599 of file lalpulsar_PiecewiseSearch.py.

◆ h0

lalpulsar_PiecewiseSearch.h0

Definition at line 601 of file lalpulsar_PiecewiseSearch.py.

◆ cosi

lalpulsar_PiecewiseSearch.cosi

Definition at line 602 of file lalpulsar_PiecewiseSearch.py.

◆ psi

lalpulsar_PiecewiseSearch.psi

Definition at line 603 of file lalpulsar_PiecewiseSearch.py.

◆ phi0

lalpulsar_PiecewiseSearch.phi0

Definition at line 604 of file lalpulsar_PiecewiseSearch.py.

◆ dt_wf

lalpulsar_PiecewiseSearch.dt_wf

Definition at line 605 of file lalpulsar_PiecewiseSearch.py.

◆ Alpha

lalpulsar_PiecewiseSearch.Alpha

Definition at line 606 of file lalpulsar_PiecewiseSearch.py.

◆ Delta

lalpulsar_PiecewiseSearch.Delta

Definition at line 607 of file lalpulsar_PiecewiseSearch.py.

◆ noise_sqrt_Sh

lalpulsar_PiecewiseSearch.noise_sqrt_Sh

Definition at line 609 of file lalpulsar_PiecewiseSearch.py.

◆ dfreq

lalpulsar_PiecewiseSearch.dfreq

Definition at line 612 of file lalpulsar_PiecewiseSearch.py.

◆ sourceDeltaT

lalpulsar_PiecewiseSearch.sourceDeltaT

Definition at line 613 of file lalpulsar_PiecewiseSearch.py.

◆ inject_data

lalpulsar_PiecewiseSearch.inject_data

Definition at line 614 of file lalpulsar_PiecewiseSearch.py.

◆ append_string

string lalpulsar_PiecewiseSearch.append_string = "_BS-" + str(args.BaseSeed) + "_j-" + str(j)

Definition at line 622 of file lalpulsar_PiecewiseSearch.py.

◆ local_disk_directory

lalpulsar_PiecewiseSearch.local_disk_directory = os.getenv("JOBFS")

Definition at line 625 of file lalpulsar_PiecewiseSearch.py.

◆ tbankdirectory

string lalpulsar_PiecewiseSearch.tbankdirectory

Initial value:

=  (
        os.path.join(local_disk_directory, tbank.toString()) + append_string
    )

Definition at line 626 of file lalpulsar_PiecewiseSearch.py.

◆ minimum_h0

lalpulsar_PiecewiseSearch.minimum_h0 = args.min_h0

Definition at line 640 of file lalpulsar_PiecewiseSearch.py.

◆ maximum_h0

lalpulsar_PiecewiseSearch.maximum_h0 = args.max_h0

Definition at line 641 of file lalpulsar_PiecewiseSearch.py.

◆ power_list

lalpulsar_PiecewiseSearch.power_list = np.linspace(np.log10(minimum_h0), np.log10(maximum_h0), iterations)

Definition at line 643 of file lalpulsar_PiecewiseSearch.py.

◆ h0_list

list lalpulsar_PiecewiseSearch.h0_list = [0] + [10**i for i in power_list]

Definition at line 645 of file lalpulsar_PiecewiseSearch.py.

◆ plus_one

int lalpulsar_PiecewiseSearch.plus_one = 0

Definition at line 647 of file lalpulsar_PiecewiseSearch.py.

◆ this_seed

int lalpulsar_PiecewiseSearch.this_seed = Random_Seed + i

Definition at line 667 of file lalpulsar_PiecewiseSearch.py.

◆ sumstr

string lalpulsar_PiecewiseSearch.sumstr = ""

Definition at line 670 of file lalpulsar_PiecewiseSearch.py.

◆ iteration_string

string lalpulsar_PiecewiseSearch.iteration_string = "_i-" + str(i)

Definition at line 674 of file lalpulsar_PiecewiseSearch.py.

◆ tempsfilename

string lalpulsar_PiecewiseSearch.tempsfilename = "Temps_For_" + tbank.toString() + iteration_string + sumstr + ".txt"

Definition at line 676 of file lalpulsar_PiecewiseSearch.py.

◆ mismatchtempname

tuple lalpulsar_PiecewiseSearch.mismatchtempname

Initial value:

=  (
        "Mismatch_Temps_For_" + tbank.toString() + iteration_string + sumstr + ".txt"
    )

Definition at line 677 of file lalpulsar_PiecewiseSearch.py.

◆ seed

lalpulsar_PiecewiseSearch.seed = random.seed(this_seed)

Definition at line 684 of file lalpulsar_PiecewiseSearch.py.

◆ signalparams

list lalpulsar_PiecewiseSearch.signalparams = lal.gsl_matrix(tbank.s * len(bf.knotslist), 1)

Definition at line 703 of file lalpulsar_PiecewiseSearch.py.

◆ SFTdirectory

lalpulsar_PiecewiseSearch.SFTdirectory

Initial value:

=  os.path.join(
            tbankdirectory,
            f"SFTs_h0-{finputdata.h0:.2e}_f0-{f0:.3f}_f1-{f1:.2e}_f2-{f2:.2e}_dur-{tbank.dur}_tstart-{tbank.tstart}/",
        )

Definition at line 703 of file lalpulsar_PiecewiseSearch.py.