Explorar el Código

added refgenome building methods

Avik Datta hace 3 años
padre
commit
3d989592b6
Se han modificado 1 ficheros con 65 adiciones y 1 borrados
  1. 65 1
      dags/dag17_create_transcriptome_ref.py

+ 65 - 1
dags/dag17_create_transcriptome_ref.py

@@ -7,6 +7,11 @@ from airflow.operators.python_operator import BranchPythonOperator
 from airflow.operators.bash_operator import BashOperator
 from igf_airflow.utils.dag17_create_transcriptome_ref_utils import download_gtf_file_func
 from igf_airflow.utils.dag17_create_transcriptome_ref_utils import create_star_index_func
+from igf_airflow.utils.dag17_create_transcriptome_ref_utils import create_rsem_index_func
+from igf_airflow.utils.dag17_create_transcriptome_ref_utils import create_reflat_index_func
+from igf_airflow.utils.dag17_create_transcriptome_ref_utils import create_ribosomal_interval_func
+from igf_airflow.utils.dag17_create_transcriptome_ref_utils import create_cellranger_ref_func
+
 
 args = {
     'owner': 'airflow',
@@ -53,4 +58,63 @@ with dag:
                 'star_options': [
                     '--sjdbOverhang', 149]
             },
-            python_callable=create_star_index_func)
+            python_callable=create_star_index_func)
+    ## TASK
+    create_rsem_index = \
+        PythonOperator(
+            task_id="create_rsem_index",
+            dag=dag,
+            queue='hpc_8G8t',
+            params={
+                'gtf_xcom_task': 'download_gtf_file',
+                'gtf_xcom_key': 'gtf_file',
+                'rsem_ref_xcom_key': 'rsem_ref',
+                'threads': 8
+            },
+            python_callable=create_rsem_index_func)
+    ## TASK
+    create_reflat_index = \
+        PythonOperator(
+            task_id='create_reflat_index',
+            dag=dag,
+            queue='hpc_4G',
+            params={
+                'gtf_xcom_task': 'download_gtf_file',
+                'gtf_xcom_key': 'gtf_file',
+                'refflat_ref_xcom_key': 'refflat_ref'
+            },
+            python_callable=create_reflat_index_func)
+    ## TASK
+    create_ribosomal_interval = \
+        PythonOperator(
+            task_id='create_ribosomal_interval',
+            dag=dag,
+            queue='hpc_4G',
+            params={
+                'gtf_xcom_task': 'download_gtf_file',
+                'gtf_xcom_key': 'gtf_file',
+                'ribosomal_ref_xcom_key': 'ribosomal_ref',
+                'skip_gtf_rows': 5
+            },
+            python_callable=create_ribosomal_interval_func)
+    ## TASK
+    create_cellranger_ref = \
+        PythonOperator(
+            task_id='create_cellranger_ref',
+            dag=dag,
+            queue='hpc_8G8t',
+            params={
+                'gtf_xcom_task': 'download_gtf_file',
+                'gtf_xcom_key': 'gtf_file',
+                'cellranger_ref_xcom_key': 'cellranger_ref',
+                'skip_gtf_rows': 5,
+                'threads': 8,
+                'memory': 8
+            },
+            python_callable=create_cellranger_ref_func)
+    ## PIPELINE
+    download_gtf_file >> create_star_index
+    download_gtf_file >> create_rsem_index
+    download_gtf_file >> create_reflat_index
+    download_gtf_file >> create_ribosomal_interval
+    download_gtf_file >> create_cellranger_ref