############################################################################### ####### lalinference_pipe example initialisation file ######################### ####### For more information, see docs page at ######################### # https://www.lsc-group.phys.uwm.edu/ligovirgo/cbcnote/ParameterEstimationModelSelection/LALInference_HowTo ############################################################################### [analysis] # IFOs to use for this analysis ifos=['H1','L1'] # Analysis type for accounting for computing use # See https://ldas-gridmon.ligo.caltech.edu/accounting/condor_groups/determine_condor_account_group.html # !!! If omitted, jobs may not run on all LDG clusters !!! #1;2c Common choices are: # ligo.dev.o1.cbc.pe.lalinferenceonline - day-latency followup for O1 development # ligo.dev.o1.cbc.pe.lalinference - follow up for detection candidates for O1 development # ligo.dev.o1.cbc.noise.lalinference - noise modelling for O1 development # ligo.dev.o1.cbc.testgr.tiger - TIGER for O1 development # # Replace dev with sim or prod for simulation or production runs # Replace o1 with s6,o2,o3 for other runs # accounting_group=ligo.dev.o2.cbc.pe.lalinference #EHT: accounting_group_user=rory.smith accounting_group_user=ethrane # select the engine to use, lalinferencenest, lalinferencemcmc or lalinferencebambimpi engine=lalinferencenest # Number of parallel runs for each event. Using parallel runs increases the number of sample, giving a better exploration of the parameter space. # Default 2 for lalinferencemcmc, 1 for lalinference_bambi and 4 for lalinferencenest. nparallel=4 # Turn the coherence test? coherence-test=True # Upload results to the gravitational-wave events database? upload-to-gracedb=False # Add the lvem tag to the sky maps and extrinsic triangle plot to share with EM partners add-lvem-tag=True #chain seed (incremented by one for each parallel chain). If commented the cpu time will be used to initialize the chain. # randomseed=4321 # noise generator seed (incremented by eventnum for each event). If set to 0 the cpu time will be used to initialize the noise generator (that may make it problematic to reproduce results). Mandatory parameter when using fake-cache in the [lalinference] section. # dataseed=1234 pegasus.transfer.links=false ##################################################################################### [paths] # the following can be given with command line (see below) #webdir is the base output dir for results pages #EHT: webdir=/home/rory.smith/public_html/monitor/offline_pe/G211117/ethraneTest webdir=/home/ethrane/public_html/files/lalinference/online_pe roq_b_matrix_directory = /home/cbc/ROQ_data/IMRPhenomPv2/ #baseurl is the www address of the above path, e.g: #EHT: baseurl=https://ldas-jobs.ligo.caltech.edu/~rory.smith/monitor/offline_pe/ethraneTest baseurl=https://ldas-jobs.ligo.caltech.edu/~ethrane/files/lalinference/online_pe/ ##################################################################################### [input] # User-specified length of the psd. if not specified, will be automatically calculated from segment availability # psd-length=1024 # User-specified psd start time (used for manual over-ride, if not given pipeline determines this automatically) # psd-start-time= 968621415 # Maximum length to use for automatically-determined psdlength options max-psd-length=1024 # spacing between trigger times and start of psd estimation padding=16 # all will analyze all the events in the injetion table. Ranges and/or single events can also be specified. E.g. events=[0,4,5,10:12] will analyze the event 0,4,5,10,11,12. events=all # Overlap between segments when analysing a continuous span of time # segment-overlap = 32 # Can manually over-ride time limits here - must use these when analysing # a stretch of time, i.e. when no other input files are given #gps-start-time= #gps-end-time= analyse-all-time=False # Can manually specify input files here or over-ride on the command line #gps-time-file= #injection-file= #sngl-inspiral-file= #coinc-inspiral-file= #pipedown-db= #gstlal-db= # If this option is set "true" and a pipedown database is used as input, the pipeline will analyse the time slide events timeslides=false # Uncomment the following line to ignore science segments. Useful when using fake noise ignore-science-segments=True # Uncomment the following line to ignore the psd.xml.gz file when analyzing a gracedb event. ignore-gracedb-psd=True # Options to select events from SQLite databases (optional) # Select combined false alarm rate in Hz # min-cfar = 0 # max-cfar = 1e-6 # Option to specify text file for dumping time slide information #time-slide-dump=path/to/file ##################################################################################### [condor] # Replace with your own executable locations datafind=/home/gracedb.processor/users/vraymond/opt/ligo_data_find_sleep_wait.sh mergescript=/home/ethrane/local/bin/lalapps_nest2pos resultspage=/home/ethrane/local/bin/cbcBayesPostProc.py segfind=/home/ethrane/local/bin/ligolw_segment_query ligolw_print=/home/ethrane/local/bin/ligolw_print coherencetest=/home/ethrane/local/bin/lalapps_coherence_test lalinferencenest=/home/ethrane/local/bin/lalinference_nest computeroqweights = /home/ethrane/local/bin/lalapps_compute_roq_weights lalinferencemcmc=/home/ethrane/local/bin/lalinference_mcmc lalinferencebambi=/home/ethrane/local/bin/lalinference_bambi lalinferencedatadump=/home/ethrane/local/bin/lalinference_datadump # Bayesline code from arXiv:1410.3852 # Disabled by default, uncomment line below to use it # bayesline=/home/rory.smith/local/bin/BayesLine # Skyareas code from https://github.com/farr/skyarea # Disabled by default, uncomment line below to use it skyarea=/bin/run_sky_area.py # MPIrun only needed for MCMC, if not available use /bin/true mpirun=/ldcg/intel/2016u1/compilers_and_libraries_2016.1.150/linux/mpi/intel64/bin/mpirun # Wrapper script for MPI jobs mpiwrapper=/home/ethrane/local/bin/lalinference_mpi_wrapper # GraceDB only needed for online analysis, if not available use /bin/true gracedb=/usr/bin/gracedb # Scripts needed only for PP check pipeline ppanalysis=/home/ethrane/local/bin/cbcBayesPPAnalysis.py pos_to_sim_inspiral=/home/ethrane/local/bin/cbcBayesPosToSimInspiral.py # Queue information if needed to run on specific nodes #queue=Online_CBC_PE_Daily ##################################################################################### # Section used by the datafind jobs (not used with simulated noise) [datafind] # Specifies we look for data in the local filesystem (as opposed to remote data) url-type=file # Can specify veto categories here (default to using no vetoes) # veto-categories= # "type" of data to look for with ligo_data_find #types={'H1':'H1_HOFT_C00','L1':'L1_HOFT_C00','V1':'V1Online'}# Offline types={'H1':'H1_HOFT_C02','L1':'L1_HOFT_C02','V1':'V1Online'}# Offline ##################################################################################### [data] # Channels to use to read h(t) from frame files # S5 has LSC-STRAIN, S6 has LDAS-STRAIN. channels={'H1':'H1:DCS-CALIB_STRAIN_C02','L1':'L1:DCS-CALIB_STRAIN_C02','V1':'V1:FAKE_h_16384Hz_4R'} ##################################################################################### [lalinference] # fake-cache is used to specify the fake names for simulated data # For simulated LIGO noise use 'LALLIGO', Virgo noise: 'LALVirgo' # For simulated advanced detectors use 'LALSimAdLIGO' and 'LALSimAdVirgo' # If not enabled, the pipeline will look for real data #fake-cache={'H1':'LALLIGO','L1':'LALLIGO','V1':'LALVirgo'} # fLow can be used to specify the low frequency cutoff for each detector flow = {'H1':40,'L1':40,'V1':40} #fhigh = {'H1':1024,'L1':1024,'V1':1024} ##################################################################################### [engine] # All options in this section are passed to lalinference_nest, lalinference_mcmc, # and lalinference_bambi. Some useful ones are noted below. # Options passed to a sampler which does not know them will simply be ignored # REQUIRED SETTINGS: # Nlive specifies the number of live points for each job nlive=512 #nmcmc = 500 nsteps=-50000 # Sampling rate for data srate=2048 # Segment length to use for analysis (should be long enough for whole template seglen=32 # OPTIONAL SETTINGS: #psdFit= # Use lalinference's resume ability if the run is interrupted. Comment out to disable. resume= # lalinference_bambi automatically resumes, use this is you want to force a start from scratch #noresume= # approx can be used to manually specify an approximant # If this is not given, the code will use whatever was injected in the case of a software injection # Or TaylorF2threePointFivePN if no injection was given. approx=IMRPhenomPv2pseudoFourPN # Control the amplitude order (default: max available) amporder=0 # maxmcmc set the maximum chain length for the nested sampling sub-chains. Default 5000 # Auto determination is on, but the length cannot be longer than that. #maxmcmc=5000 # Number of independent sample to obtain before quitting for lalinference_mcmc neff=1000 #differential-buffer-limit=1000000 # Priors # For all parameters known to lalinference, the min and max default prior can be overwritten with #parname-min=MIN #parname-max=MAX #distance-max=300 distance-max=2000 #mass1-min=1.001398 #mass1-max=30 #mass2-min=1.001398 #mass2-max=4 a_spin1-min=0 a_spin1-max=0.8 a_spin2-min=0 a_spin2-max=0.8 # The starting point for the MCMC chain(s) can be specified with #parname=START # Parameters can be fixed to some value with #fix-parname=FIXVALUE #currently known parameters, together with default [min-max] are (radiants for angle, Mpc for distance, Msun for masses) #time Waveform time [trigtime-0.1-trigtime+0.1] #chirpmass Chirpmass [1,15.3] #eta Symmetric massratio (needs --use-eta) [0,0.0312] #q Asymmetric massratio (a.k.a. q=m2/m1 with m1>m2) [0.033,1] #phase Coalescence phase [0,2Pi] #costheta_jn Cosine of angle between J and line of sight [-1,1] #logdistance Log Distance [log(1),log(2000)] #rightascension Rightascension [0,2Pi] #declination Declination [-Pi/2,Pi/2] #polarisation Polarisation angle [0,Pi] #Spin Parameters: #a_spin1 Spin1 magnitude [-1,1] for aligned, [0,1] for precessing #a_spin2 Spin2 magnitude [-1,1] for aligned, [0,1] for precessing #tilt_spin1 Angle between spin1 and orbital angular momentum [0,Pi] #tilt_spin2 Angle between spin2 and orbital angular momentum [0, Pi] #phi_12 Difference between spins' azimuthal angles [0,2Pi] #phi_jl Difference between total and orbital angular momentum azimuthal angles [0,2Pi] #Equation of State parameters (requires --use-tidal or --use-tidalT): #lambda1 lambda1 [0,3000] #lambda2 lambda2 [0,3000] #lambdaT lambdaT [0,3000] #dLambdaT dLambdaT [-500,500] # Settings for allowed component masses in Solar Masses, with default values comp-max=80 comp-min=1.001398 # Allowed total masses in Msun (note, used together with component masses, mc,q,eta priors may lead to inconsistencies. Be careful!) #mtotal-max=35 #mtotal-min=2 # Setting time prior [secs] #dt=0.1 # The following three options control various marginalised likelihoods. Use at most one. # Analytically marginalise over phase (only for Newtonian amplitude orders) #margphi= # Semi-analytically marginalise over time #margtime= # Semi-analytically marginalise over time and phase (only for Newtonian amplitude orders) #margtimephi= # By default the code use spins if the choosen approximant supports spin. NOTE that this include TaylorF2, which will be run with aligned spins. # Several options, here below, allows the user to choose a particular setting: #Disable spin for waveforms which would normally have spins #disable-spin= # Only allow aligned spins #aligned-spin= # Only allow the heavier mass to spin (can be used together with aligned-spin) #singleSpin= # Print progress information throughout the run progress= # lalinference_bambi allows you to set a target efficiency and tolerance - these are default values #eff=0.1 tolerance=0.1 # Sample in eta instead than q=m2/m1 #use-eta= ##################################################################################### [mpi] # Settings when running with MPI for lalinference_mcmc or lalinference_bambi # number of CPUs desired and how much memory on each (MB) mpi_task_count=1 machine-count=8 machine-memory=2000 ##################################################################################### [resultspage] # Settings for the results pages (see cbcBayesPostProc.py --help for more) # Size of the side of the sky bins in degrees skyres=0.5 # Create 2D plots, which take a while to generate # plot-2d= # For lalinferencemcmc, select points for plotting only after the chains have # reached a log likelihood of at least max_log_likelihood_found-. # This deltaLogL value should be the number of signal parameters divided by 2. deltaLogL=5 # Send an email linking to the final page #email=albert.einstein@ligo.org ##################################################################################### [segfind] # URL of the segment database segment-url=https://dqsegdb5.phy.syr.edu ##################################################################################### [segments] # Names of science segments to use l1-analyze = L1:DMT-ANALYSIS_READY:1 h1-analyze = H1:DMT-ANALYSIS_READY:1 v1-analyze = V1:ITF_SCIENCEMODE:7 [ppanalysis] # This section only used with lalinference_pp_pipe # Webdir is the output path for the PP check page webdir=/home/albert.einstein/public_html/lalinference/review/PPcheck/ [skyarea] # Optional arguments: #maxpts=1000 # maximum number of posterior samples to use #trials=50 # Maximum number of trials to build sky posterior #noskyarea= # Turn off sky area computation #nside=512 # HEALPix resultion (default 512) #seed=1234 # Random seed to use