[Swift-commit] r2852 - SwiftApps/SIDGrid/swift/projects/andric/peakfit_pilots/PK2

[noreply at svn.ci.uchicago.edu]

Author: andric
Date: 2009-04-09 16:35:01 -0500 (Thu, 09 Apr 2009)
New Revision: 2852

Added:
   SwiftApps/SIDGrid/swift/projects/andric/peakfit_pilots/PK2/PK2.swift
Log:
for doing peak analysis

Added: SwiftApps/SIDGrid/swift/projects/andric/peakfit_pilots/PK2/PK2.swift
===================================================================
--- SwiftApps/SIDGrid/swift/projects/andric/peakfit_pilots/PK2/PK2.swift	                        (rev 0)
+++ SwiftApps/SIDGrid/swift/projects/andric/peakfit_pilots/PK2/PK2.swift	2009-04-09 21:35:01 UTC (rev 2852)
@@ -0,0 +1,36 @@
+#----  doing peakfit on PK2 pilot subject
+type file{}
+
+type PeakResult{
+    file pdf;
+    file a;
+    file b;
+}
+
+app (PeakResult outObj) peakScript (file scriptFile, file tsfile, file preprocR, file peakfitR){
+    RInvoke @filename(scriptFile) @filename(tsfile) @filename(outObj);
+}
+
+
+file script<"peak_scripts/Shellpeak_PK2runbyrun.R">;
+file preprocR<"peak_scripts/preprocess.Enewsmooth.R">;
+file peakfitR<"peak_scripts/peakfit_2.1.R">;
+
+string declarelist[] = ["PK2"];
+foreach subject in declarelist{
+    int runs[] = [1:8];
+    foreach rn in runs{
+        string hemilist[] = ["lh","rh"];
+        foreach h in hemilist{
+            int inputRegionNo[] = [15,33,34,60];
+            foreach rr in inputRegionNo{
+                string project_dir=@strcat("/disks/ci-gpfs/fmri/cnari/swift/projects/andric/peakfit_pilots/",subject);
+                ## map inputs
+                file inFile<single_file_mapper; file=@strcat(project_dir,"/input/",subject,"_peakTS_data",h,".",rr,"run",rn)>;
+                ## map results
+                PeakResult outObj<simple_mapper; prefix=@strcat(subject,".",rr,"_",h,"run",rn,".result.")>;
+                (outObj)=peakScript(script,inFile, preprocR, peakfitR);
+            }
+        }
+    }
+}