"--output", dqVetoes["combined-veto-file"][category] ] + veto_cat_files) inspiralutils.make_external_call(concat_veto_call) # FIXME: remove this block when the veto files from # ligolw_segments_from_cats have contain a segment table with start_time_ns # and end_time_ns, and when the process table doesn't need the domain, jobid, # and is_online columns added compat_veto_call = " ".join([cp.get("condor", "ligolw_segments_compat"), dqVetoes["combined-veto-file"][category]]) inspiralutils.make_external_call(compat_veto_call) segFile = {} for ifo in ifos: segFile[ifo], dqVetoes[ifo] = inspiralutils.findSegmentsToAnalyze( cp, ifo, veto_categories, opts.generate_segments, opts.use_available_data, opts.generate_veto_segments) # set correct name as we're in the segments directory: segFile[ifo] = "../segments/" + segFile[ifo] cp.set("input", ifo.lower() + "-segments", segFile[ifo]) for key in dqVetoes[ifo].keys(): dqVetoes[ifo][key] = "../segments/" + dqVetoes[ifo][key] for key in dqVetoes["combined-veto-file"].keys(): dqVetoes["combined-veto-file"][key] = "../segments/" + dqVetoes["combined-veto-file"][key] os.chdir("..") hw_inj_dir = "hardware_injection_summary" if opts.run_hardware_inj: inspiralutils.mkdir(hw_inj_dir) inspiralutils.get_hwinj_segments(cp,ifos,hw_inj_dir)
from glue import segmentsUtils parser=OptionParser() opts,args = parser.parse_args() config=ConfigParser.ConfigParser() config.read(args[0]) if config.has_option('analysis','ifos'): ifos=ast.literal_eval(config.get('analysis','ifos')) else: ifos=['H1','L1']#,'V1'] segments={} if config.has_option('datafind','veto-categories'): veto_categories=ast.literal_eval(config.get('datafind','veto-categories')) else: veto_categories=[] for ifo in ifos: (segFileName,dqVetoes)=inspiralutils.findSegmentsToAnalyze(config, ifo, veto_categories, generate_segments=True, use_available_data=False, data_quality_vetoes=False) segfile=open(segFileName) segments[ifo]=segmentsUtils.fromsegwizard(segfile) print segments[ifo]
vetoCats = ast.literal_eval(cp.get('segments', 'veto-categories')) vetoDefFile = inspiralutils.downloadVetoDefFile(cp, True) inspiralutils.generate_veto_cat_files(cp, vetoDefFile, True) # Populate dictionary with veto filenames for each IFO and generate the veto segment files veto_dic = {} seg_dic = {} # Download science segment files for each IFO for IFO in ifos: (seg_dic[IFO], veto_dic[IFO]) = inspiralutils.findSegmentsToAnalyze( cp, IFO, vetoCats, generate_segments=True, use_available_data=False, data_quality_vetoes=True) # Feed science and veto segment files to make_injtimes and generate injection times and timeslides print 'TIGER: Generating GPS times for unvetoed injections' ensure_dir(os.path.join(basefolder, 'injtimes')) IFOdict = {} for ifo in ifos: IFOdict[ifo] = make_injtimes.IFO(ifo, seg_dic[ifo], veto_dic[ifo][4], minlen=max(psdlen, seglen)) os.chdir(curdir)
gps_end_time)]) else: # # --- Run segdb query # if cp.has_option('datafind','veto-categories'): veto_categories=ast.literal_eval(cp.get('datafind','veto-categories')) else: veto_categories=[] curdir=os.getcwd() os.chdir(segment_dir) (segFileName,dqVetoes)=inspiralutils.findSegmentsToAnalyze(cp, ifo, veto_categories, generate_segments=True, use_available_data=False, data_quality_vetoes=False) segfile=open(segFileName) segmentList[ifo]=segmentsUtils.fromsegwizard(segfile) segmentList[ifo].coalesce() segfile.close() if segmentList[ifo] == []: print >> sys.stderr, "No matching segments for %s"%ifo sys.exit() os.chdir(curdir) # --------------------------------------------------------------------
curdir = os.getcwd() os.chdir(segfolder) vetoCats = ast.literal_eval(cp.get('segments', 'veto-categories')) vetoDefFile = inspiralutils.downloadVetoDefFile(cp, True) inspiralutils.generate_veto_cat_files(cp, vetoDefFile, True) # Populate dictionary with veto filenames for each IFO and generate the veto segment files veto_dic = {} seg_dic = {} # Download science segment files for each IFO for IFO in ifos: (seg_dic[IFO], veto_dic[IFO]) = inspiralutils.findSegmentsToAnalyze(cp, IFO, vetoCats, generate_segments=True, use_available_data=False, data_quality_vetoes=True) # Feed science and veto segment files to make_injtimes and generate injection times and timeslides print 'TIGER: Generating GPS times for unvetoed injections' ensure_dir(os.path.join(basefolder, 'injtimes')) IFOdict={} for ifo in ifos: IFOdict[ifo]=make_injtimes.IFO(ifo, seg_dic[ifo], veto_dic[ifo][4], minlen=max(psdlen, seglen)) os.chdir(curdir) # Generate injtimes file if timeslides: print 'TIGER: Injection times will be generated with timeslides' timesdict = {}