kkonganti@92: // Hold methods to print: kkonganti@92: // 1. Colored logo. kkonganti@92: // 2. Summary of parameters. kkonganti@92: // 3. Single dashed line. kkonganti@92: // 4. Double dashed line. kkonganti@92: // kkonganti@92: kkonganti@92: import groovy.json.JsonSlurper kkonganti@92: import nextflow.config.ConfigParser kkonganti@92: // import groovy.json.JsonOutput kkonganti@92: kkonganti@92: // ASCII logo kkonganti@92: def pipelineBanner() { kkonganti@92: kkonganti@92: def padding = (params.pad) ?: 30 kkonganti@92: Map fgcolors = getANSIColors() kkonganti@92: kkonganti@92: def banner = [ kkonganti@92: name: "${fgcolors.magenta}${workflow.manifest.name}${fgcolors.reset}", kkonganti@92: author: "${fgcolors.cyan}${workflow.manifest.author}${fgcolors.reset}", kkonganti@92: // workflow: "${fgcolors.magenta}${params.pipeline}${fgcolors.reset}", kkonganti@92: version: "${fgcolors.green}${workflow.manifest.version}${fgcolors.reset}", kkonganti@92: center: "${fgcolors.green}${params.center}${fgcolors.reset}", kkonganti@92: pad: padding kkonganti@92: ] kkonganti@92: kkonganti@92: manifest = addPadding(banner) kkonganti@92: kkonganti@92: return """${fgcolors.white}${dashedLine(type: '=')}${fgcolors.magenta} kkonganti@92: (o) kkonganti@92: ___ _ __ _ _ __ ___ ___ kkonganti@92: / __|| '_ \\ | || '_ \\ / _ \\/ __| kkonganti@92: | (__ | |_) || || |_) || __/\\__ \\ kkonganti@92: \\___|| .__/ |_|| .__/ \\___||___/ kkonganti@92: | | | | kkonganti@92: |_| |_|${fgcolors.reset} kkonganti@92: ${dashedLine()} kkonganti@92: ${fgcolors.blue}A collection of modular pipelines at CFSAN, FDA.${fgcolors.reset} kkonganti@92: ${dashedLine()} kkonganti@92: ${manifest} kkonganti@92: ${dashedLine(type: '=')} kkonganti@92: """.stripIndent() kkonganti@92: } kkonganti@92: kkonganti@92: // Add padding to keys so that kkonganti@92: // they indent nicely on the kkonganti@92: // terminal kkonganti@92: def addPadding(values) { kkonganti@92: kkonganti@92: def pad = (params.pad) ?: 30 kkonganti@92: values.pad = pad kkonganti@92: kkonganti@92: def padding = values.pad.toInteger() kkonganti@92: def nocapitalize = values.nocapitalize kkonganti@92: def stopnow = values.stopNow kkonganti@92: def help = values.help kkonganti@92: kkonganti@92: values.removeAll { kkonganti@92: k, v -> [ kkonganti@92: 'nocapitalize', kkonganti@92: 'pad', kkonganti@92: 'stopNow', kkonganti@92: 'help' kkonganti@92: ].contains(k) kkonganti@92: } kkonganti@92: kkonganti@92: values.keySet().each { k -> kkonganti@92: v = values[k] kkonganti@92: s = params.linewidth - (pad + 5) kkonganti@92: if (v.toString().size() > s && !stopnow) { kkonganti@92: def sen = '' kkonganti@92: v.toString().findAll(/.{1,${s}}\b(?:\W*|\s*)/).each { kkonganti@92: sen += ' '.multiply(padding + 2) + it + '\n' kkonganti@92: } kkonganti@92: values[k] = ( kkonganti@92: help ? sen.replaceAll(/^(\n|\s)*/, '') : sen.trim() kkonganti@92: ) kkonganti@92: } else { kkonganti@92: values[k] = (help ? v + "\n" : v) kkonganti@92: } kkonganti@92: k = k.replaceAll(/\./, '_') kkonganti@92: } kkonganti@92: kkonganti@92: return values.findResults { kkonganti@92: k, v -> nocapitalize ? kkonganti@92: k.padRight(padding) + ': ' + v : kkonganti@92: k.capitalize().padRight(padding) + ': ' + v kkonganti@92: }.join("\n") kkonganti@92: } kkonganti@92: kkonganti@92: // Method for error messages kkonganti@92: def stopNow(msg) { kkonganti@92: kkonganti@92: Map fgcolors = getANSIColors() kkonganti@92: Map errors = [:] kkonganti@92: kkonganti@92: if (msg == null) { kkonganti@92: msg = "Unknown error" kkonganti@92: } kkonganti@92: kkonganti@92: errors['stopNow'] = true kkonganti@92: errors["${params.cfsanpipename} - ${params.pipeline} - ERROR"] = """ kkonganti@92: ${fgcolors.reset}${dashedLine()} kkonganti@92: ${fgcolors.red}${msg}${fgcolors.reset} kkonganti@92: ${dashedLine()} kkonganti@92: """.stripIndent() kkonganti@92: // println dashedLine() // defaults to stdout kkonganti@92: // log.info addPadding(errors) // prints to stdout kkonganti@92: exit 1, "\n" + dashedLine() + kkonganti@92: "${fgcolors.red}\n" + addPadding(errors) kkonganti@92: } kkonganti@92: kkonganti@92: // Method to validate 4 required parameters kkonganti@92: // if input for entry point is FASTQ files kkonganti@92: def validateParamsForFASTQ() { kkonganti@92: switch (params) { kkonganti@92: case { params.metadata == null && params.input == null }: kkonganti@92: stopNow("Either metadata CSV file with 5 required columns\n" + kkonganti@92: "in order: sample, fq1, fq2, strandedness, single_end or \n" + kkonganti@92: "input directory of only FASTQ files (gzipped or unzipped) should be provided\n" + kkonganti@92: "using --metadata or --input options.\n" + kkonganti@92: "None of these two options were provided!") kkonganti@92: break kkonganti@92: case { params.metadata != null && params.input != null }: kkonganti@92: stopNow("Either metadata or input directory of FASTQ files\n" + kkonganti@92: "should be provided using --metadata or --input options.\n" + kkonganti@92: "Using both these options is not allowed!") kkonganti@92: break kkonganti@92: case { params.output == null }: kkonganti@92: stopNow("Please mention output directory to store all results " + kkonganti@92: "using --output option!") kkonganti@92: break kkonganti@92: } kkonganti@92: return 1 kkonganti@92: } kkonganti@92: kkonganti@92: // Method to print summary of parameters kkonganti@92: // before running kkonganti@92: def summaryOfParams() { kkonganti@92: kkonganti@92: def pipeline_specific_config = new ConfigParser().setIgnoreIncludes(true).parse( kkonganti@92: file("${params.workflowsconf}${params.fs}${params.pipeline}.config").text kkonganti@92: ) kkonganti@92: Map fgcolors = getANSIColors() kkonganti@92: Map globalparams = [:] kkonganti@92: Map localparams = params.subMap( kkonganti@92: pipeline_specific_config.params.keySet().toList() + params.logtheseparams kkonganti@92: ) kkonganti@92: kkonganti@92: if (localparams !instanceof Map) { kkonganti@92: stopNow("Need a Map of paramters. We got: " + localparams.getClass()) kkonganti@92: } kkonganti@92: kkonganti@92: if (localparams.size() != 0) { kkonganti@92: localparams['nocapitalize'] = true kkonganti@92: globalparams['nocapitalize'] = true kkonganti@92: globalparams['nextflow_version'] = "${nextflow.version}" kkonganti@92: globalparams['nextflow_build'] = "${nextflow.build}" kkonganti@92: globalparams['nextflow_timestamp'] = "${nextflow.timestamp}" kkonganti@92: globalparams['workflow_projectDir'] = "${workflow.projectDir}" kkonganti@92: globalparams['workflow_launchDir'] = "${workflow.launchDir}" kkonganti@92: globalparams['workflow_workDir'] = "${workflow.workDir}" kkonganti@92: globalparams['workflow_container'] = "${workflow.container}" kkonganti@92: globalparams['workflow_containerEngine'] = "${workflow.containerEngine}" kkonganti@92: globalparams['workflow_runName'] = "${workflow.runName}" kkonganti@92: globalparams['workflow_sessionId'] = "${workflow.sessionId}" kkonganti@92: globalparams['workflow_profile'] = "${workflow.profile}" kkonganti@92: globalparams['workflow_start'] = "${workflow.start}" kkonganti@92: globalparams['workflow_commandLine'] = "${workflow.commandLine}" kkonganti@92: return """${dashedLine()} kkonganti@92: Summary of the current workflow (${fgcolors.magenta}${params.pipeline}${fgcolors.reset}) parameters kkonganti@92: ${dashedLine()} kkonganti@92: ${addPadding(localparams)} kkonganti@92: ${dashedLine()} kkonganti@92: ${fgcolors.cyan}N E X T F L O W${fgcolors.reset} - ${fgcolors.magenta}${params.cfsanpipename}${fgcolors.reset} - Runtime metadata kkonganti@92: ${dashedLine()} kkonganti@92: ${addPadding(globalparams)} kkonganti@92: ${dashedLine()}""".stripIndent() kkonganti@92: } kkonganti@92: return 1 kkonganti@92: } kkonganti@92: kkonganti@92: // Method to display kkonganti@92: // Return dashed line either '-' kkonganti@92: // type or '=' type kkonganti@92: def dashedLine(Map defaults = [:]) { kkonganti@92: kkonganti@92: Map fgcolors = getANSIColors() kkonganti@92: def line = [color: 'white', type: '-'] kkonganti@92: kkonganti@92: if (!defaults.isEmpty()) { kkonganti@92: line.putAll(defaults) kkonganti@92: } kkonganti@92: kkonganti@92: return fgcolors."${line.color}" + kkonganti@92: "${line.type}".multiply(params.linewidth) + kkonganti@92: fgcolors.reset kkonganti@92: } kkonganti@92: kkonganti@92: // Return slurped keys parsed from JSON kkonganti@92: def slurpJson(file) { kkonganti@92: def slurped = null kkonganti@92: def jsonInst = new JsonSlurper() kkonganti@92: kkonganti@92: try { kkonganti@92: slurped = jsonInst.parse(new File ("${file}")) kkonganti@92: } kkonganti@92: catch (Exception e) { kkonganti@92: log.error 'Please check your JSON schema. Invalid JSON file: ' + file kkonganti@92: } kkonganti@92: kkonganti@92: // Declare globals for the nanofactory kkonganti@92: // workflow. kkonganti@92: return [keys: slurped.keySet().toList(), cparams: slurped] kkonganti@92: } kkonganti@92: kkonganti@92: // Default help text in a map if the entry point kkonganti@92: // to a pipeline is FASTQ files. kkonganti@92: def fastqEntryPointHelp() { kkonganti@92: kkonganti@92: Map helptext = [:] kkonganti@92: Map fgcolors = getANSIColors() kkonganti@92: kkonganti@92: helptext['Workflow'] = "${fgcolors.magenta}${params.pipeline}${fgcolors.reset}" kkonganti@92: helptext['Author'] = "${fgcolors.cyan}${params.workflow_built_by}${fgcolors.reset}" kkonganti@92: helptext['Version'] = "${fgcolors.green}${params.workflow_version}${fgcolors.reset}\n" kkonganti@92: helptext['Usage'] = "cpipes --pipeline ${params.pipeline} [options]\n" kkonganti@92: helptext['Required'] = "" kkonganti@92: helptext['--input'] = "Absolute path to directory containing FASTQ files. " + kkonganti@92: "The directory should contain only FASTQ files as all the " + kkonganti@92: "files within the mentioned directory will be read. " + kkonganti@92: "Ex: --input /path/to/fastq_pass" kkonganti@92: helptext['--output'] = "Absolute path to directory where all the pipeline " + kkonganti@92: "outputs should be stored. Ex: --output /path/to/output" kkonganti@92: helptext['Other options'] = "" kkonganti@92: helptext['--metadata'] = "Absolute path to metadata CSV file containing five " + kkonganti@92: "mandatory columns: sample,fq1,fq2,strandedness,single_end. The fq1 and fq2 " + kkonganti@92: "columns contain absolute paths to the FASTQ files. This option can be used in place " + kkonganti@92: "of --input option. This is rare. Ex: --metadata samplesheet.csv" kkonganti@92: helptext['--fq_suffix'] = "The suffix of FASTQ files (Unpaired reads or R1 reads or Long reads) if " + kkonganti@92: "an input directory is mentioned via --input option. Default: ${params.fq_suffix}" kkonganti@92: helptext['--fq2_suffix'] = "The suffix of FASTQ files (Paired-end reads or R2 reads) if an input directory is mentioned via " + kkonganti@92: "--input option. Default: ${params.fq2_suffix}" kkonganti@92: helptext['--fq_filter_by_len'] = "Remove FASTQ reads that are less than this many bases. " + kkonganti@92: "Default: ${params.fq_filter_by_len}" kkonganti@92: helptext['--fq_strandedness'] = "The strandedness of the sequencing run. This is mostly needed " + kkonganti@92: "if your sequencing run is RNA-SEQ. For most of the other runs, it is probably safe to use " + kkonganti@92: "unstranded for the option. Default: ${params.fq_strandedness}" kkonganti@92: helptext['--fq_single_end'] = "SINGLE-END information will be auto-detected but this option forces " + kkonganti@92: "PAIRED-END FASTQ files to be treated as SINGLE-END so only read 1 information is included in " + kkonganti@92: "auto-generated samplesheet. Default: ${params.fq_single_end}" kkonganti@92: helptext['--fq_filename_delim'] = "Delimiter by which the file name is split to obtain sample name. " + kkonganti@92: "Default: ${params.fq_filename_delim}" kkonganti@92: helptext['--fq_filename_delim_idx'] = "After splitting FASTQ file name by using the --fq_filename_delim option," + kkonganti@92: " all elements before this index (1-based) will be joined to create final sample name." + kkonganti@92: " Default: ${params.fq_filename_delim_idx}" kkonganti@92: kkonganti@92: return helptext kkonganti@92: } kkonganti@92: kkonganti@92: // Wrap help text with the following options kkonganti@92: def wrapUpHelp() { kkonganti@92: kkonganti@92: return [ kkonganti@92: 'Help options' : "", kkonganti@92: '--help': "Display this message.\n", kkonganti@92: 'help': true, kkonganti@92: 'nocapitalize': true kkonganti@92: ] kkonganti@92: } kkonganti@92: kkonganti@92: // Method to send email on workflow complete. kkonganti@92: def sendMail() { kkonganti@92: kkonganti@92: if (params.user_email == null) { kkonganti@92: return 1 kkonganti@92: } kkonganti@92: kkonganti@92: def pad = (params.pad) ?: 30 kkonganti@92: def contact_emails = [ kkonganti@92: stakeholder: (params.workflow_blueprint_by ?: 'Not defined'), kkonganti@92: author: (params.workflow_built_by ?: 'Not defined') kkonganti@92: ] kkonganti@92: def msg = """ kkonganti@92: ${pipelineBanner()} kkonganti@92: ${summaryOfParams()} kkonganti@92: ${params.cfsanpipename} - ${params.pipeline} kkonganti@92: ${dashedLine()} kkonganti@92: Please check the following directory for N E X T F L O W kkonganti@92: reports. You can view the HTML files directly by double clicking kkonganti@92: them on your workstation. kkonganti@92: ${dashedLine()} kkonganti@92: ${params.tracereportsdir} kkonganti@92: ${dashedLine()} kkonganti@92: Please send any bug reports to CFSAN Dev Team or the author or kkonganti@92: the stakeholder of the current pipeline. kkonganti@92: ${dashedLine()} kkonganti@92: Error messages (if any) kkonganti@92: ${dashedLine()} kkonganti@92: ${workflow.errorMessage} kkonganti@92: ${workflow.errorReport} kkonganti@92: ${dashedLine()} kkonganti@92: Contact emails kkonganti@92: ${dashedLine()} kkonganti@92: ${addPadding(contact_emails)} kkonganti@92: ${dashedLine()} kkonganti@92: Thank you for using ${params.cfsanpipename} - ${params.pipeline}! kkonganti@92: ${dashedLine()} kkonganti@92: """.stripIndent() kkonganti@92: kkonganti@92: def mail_cmd = [ kkonganti@92: 'sendmail', kkonganti@92: '-f', 'cfsan-hpc-noreply@fda.hhs.gov', kkonganti@92: '-F', 'cfsan-hpc-noreply', kkonganti@92: '-t', "${params.user_email}" kkonganti@92: ] kkonganti@92: kkonganti@92: def email_subject = "${params.cfsanpipename} - ${params.pipeline}" kkonganti@92: Map fgcolors = getANSIColors() kkonganti@92: kkonganti@92: if (workflow.success) { kkonganti@92: email_subject += ' completed successfully!' kkonganti@92: } kkonganti@92: else if (!workflow.success) { kkonganti@92: email_subject += ' has failed!' kkonganti@92: } kkonganti@92: kkonganti@92: try { kkonganti@92: ['env', 'bash'].execute() << """${mail_cmd.join(' ')} kkonganti@92: Subject: ${email_subject} kkonganti@92: Mime-Version: 1.0 kkonganti@92: Content-Type: text/html kkonganti@92:
kkonganti@92: ${msg.replaceAll(/\x1b\[[0-9;]*m/, '')}
kkonganti@92: 
kkonganti@92: """.stripIndent() kkonganti@92: } catch (all) { kkonganti@92: def warning_msg = "${fgcolors.yellow}${params.cfsanpipename} - ${params.pipeline} - WARNING" kkonganti@92: .padRight(pad) + ':' kkonganti@92: log.info """ kkonganti@92: ${dashedLine()} kkonganti@92: ${warning_msg} kkonganti@92: ${dashedLine()} kkonganti@92: Could not send mail with the sendmail command! kkonganti@92: ${dashedLine()} kkonganti@92: """.stripIndent() kkonganti@92: } kkonganti@92: return 1 kkonganti@92: } kkonganti@92: kkonganti@92: // Set ANSI colors for any and all kkonganti@92: // STDOUT or STDERR kkonganti@92: def getANSIColors() { kkonganti@92: kkonganti@92: Map fgcolors = [:] kkonganti@92: kkonganti@92: fgcolors['reset'] = "\033[0m" kkonganti@92: fgcolors['black'] = "\033[0;30m" kkonganti@92: fgcolors['red'] = "\033[0;31m" kkonganti@92: fgcolors['green'] = "\033[0;32m" kkonganti@92: fgcolors['yellow'] = "\033[0;33m" kkonganti@92: fgcolors['blue'] = "\033[0;34m" kkonganti@92: fgcolors['magenta'] = "\033[0;35m" kkonganti@92: fgcolors['cyan'] = "\033[0;36m" kkonganti@92: fgcolors['white'] = "\033[0;37m" kkonganti@92: kkonganti@92: return fgcolors kkonganti@92: }