[Rivet-svn] r2231 - validation/scripts

[blackhole at projects.hepforge.org]

Author: hoeth
Date: Mon Feb  1 14:19:59 2010
New Revision: 2231

Log:
Script for merging output files from different pythia runs.

Added:
   validation/scripts/merge_runs_pythia6.sh   (contents, props changed)
Modified:
   validation/scripts/README

Modified: validation/scripts/README
==============================================================================
--- validation/scripts/README	Mon Feb  1 14:11:49 2010	(r2230)
+++ validation/scripts/README	Mon Feb  1 14:19:59 2010	(r2231)
@@ -1,6 +1,6 @@
-The scripts are templates wich XXX being the placeholder for a run
-number or run id. One possible way to create scripts that can directly
-be submitted to a batch system is:
+The scripts in the $GENERATOR directories are templates wich XXX being
+the placeholder for a run number or run id. One possible way to create
+scripts that can directly be submitted to a batch system is:
 
 for n in MC09 Perugia0 Perugia6 Perugiastar profCTEQ6ll profLOstar profMRSTMCal profQ2 profpT tuneA tuneDW
   do
@@ -10,3 +10,7 @@
         chmod +x ../scripts/${name/XXX/$n}
       done
   done
+
+The merge_runs*.sh scripts take all the aida files of the separate runs
+and merge them into a single out.aida file.
+

Added: validation/scripts/merge_runs_pythia6.sh
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ validation/scripts/merge_runs_pythia6.sh	Mon Feb  1 14:19:59 2010	(r2231)
@@ -0,0 +1,93 @@
+#! /bin/bash
+
+### This script merges the output files for various Pythia6 runs
+### into a single .aida file.
+###
+### Note that for many analyses the histogram merging scale is
+### not equal to the kinematic cut applied on generator level.
+
+### Clear the output file
+rm -f out.{dat,aida}
+
+### Get LEP-1 analyses
+rivet-rmgaps `rivet-config --datadir`/ALEPH_1996_S3486095.aida LEP_91GeV/Rivet.aida tmp.aida
+aida2flat tmp.aida >> out.dat
+rm tmp.aida
+
+### Get CDF_2000_S4155203
+aida2flat CDF_2000_S4155203/Rivet.aida >>  out.dat
+
+### Get CDF_2001_S4751469
+rivet-mergeruns CDF_2001_S4751469/Rivet.aida:1800:0 \
+                CDF_2001_S4751469_CKIN10/Rivet.aida:1800:18 \
+                CDF_2001_S4751469_CKIN20/Rivet.aida:1800:30 \
+                >> out.dat
+
+### Get CDF_2002_S4796047
+aida2flat -m d01-x01-y01 -m d03-x01-y01 CDF_2002_S4796047_630/Rivet.aida >> out.dat
+aida2flat -m d02-x01-y01 -m d04-x01-y01 CDF_2002_S4796047_1800/Rivet.aida >> out.dat
+
+### Get CDF_2004_S5839831
+rivet-mergeruns CDF_2004_S5839831_630/Rivet.aida:630:0 \
+                CDF_2004_S5839831_630_CKIN10/Rivet.aida:630:20 \
+                CDF_2004_S5839831_630_CKIN20/Rivet.aida:630:35 \
+                CDF_2004_S5839831_630_CKIN40/Rivet.aida:630:50 \
+                CDF_2004_S5839831_1800/Rivet.aida:1800:0 \
+                CDF_2004_S5839831_1800_CKIN25/Rivet.aida:1800:40 \
+                CDF_2004_S5839831_1800_CKIN60/Rivet.aida:1800:80 \
+                CDF_2004_S5839831_1800_CKIN100/Rivet.aida:1800:120 \
+                CDF_2004_S5839831_1800_CKIN140/Rivet.aida:1800:160 \
+                CDF_2004_S5839831_1800_CKIN180/Rivet.aida:1800:200 \
+                >> out.dat
+
+### Get CDF_2008_LEADINGJETS
+rivet-mergeruns CDF_2008_LEADINGJETS_CDF_2009_S8233977/Rivet.aida:1960:0 \
+                CDF_2008_LEADINGJETS_CKIN10/Rivet.aida:1960:30 \
+                CDF_2008_LEADINGJETS_CKIN20/Rivet.aida:1960:50 \
+                CDF_2008_LEADINGJETS_CKIN50/Rivet.aida:1960:80 \
+                CDF_2008_LEADINGJETS_CKIN100/Rivet.aida:1960:130 \
+                CDF_2008_LEADINGJETS_CKIN150/Rivet.aida:1960:180 \
+                > tmp.dat
+flat2aida tmp.dat
+aida2flat -m CDF_2008_LEADINGJETS tmp.aida >> out.dat
+rm tmp.{aida,dat}
+
+### Get CDF_2008_NOTE_9351
+aida2flat CDF_2008_NOTE_9351/Rivet.aida >> out.dat
+
+### Get CDF_2009_S8233977
+aida2flat -m CDF_2009_S8233977 CDF_2008_LEADINGJETS_CDF_2009_S8233977/Rivet.aida >> out.dat
+
+### Get D0_2004_S5992206
+rivet-mergeruns D0_2004_S5992206_CKIN50/Rivet.aida:1960:50 \
+                D0_2004_S5992206_CKIN75/Rivet.aida:1960:75 \
+                D0_2004_S5992206_CKIN100/Rivet.aida:1960:100 \
+                D0_2004_S5992206_CKIN150/Rivet.aida:1960:150 \
+                > tmp.dat
+flat2aida tmp.dat
+aida2flat -m D0_2004_S5992206 tmp.aida >> out.dat
+rm tmp.{aida,dat}
+
+### Get STAR_2006_S6500200
+aida2flat -m STAR_2006_S6500200 STAR_2006_S6500200_STAR_2006_S6870392/Rivet.aida >> out.dat
+
+### Get STAR_2006_S6870392
+rivet-mergeruns STAR_2006_S6500200_STAR_2006_S6870392/Rivet.aida:200:0 \
+                STAR_2009_UE_HELEN_STAR_2006_S6870392_CKIN3/Rivet.aida:200:3 \
+                >> tmp.dat
+flat2aida tmp.dat
+aida2flat -m STAR_2006_S6870392 tmp.aida >> out.dat
+rm tmp.{aida,dat}
+
+### Get STAR_2009_UE_HELEN
+rivet-mergeruns STAR_2009_UE_HELEN_STAR_2006_S6870392_CKIN3/Rivet.aida:200:5 \
+                STAR_2009_UE_HELEN_CKIN5/Rivet.aida:200:10 \
+                STAR_2009_UE_HELEN_CKIN10/Rivet.aida:200:15 \
+                STAR_2009_UE_HELEN_CKIN17/Rivet.aida:200:25 \
+                >> tmp.dat
+flat2aida tmp.dat
+aida2flat -m STAR_2009_UE_HELEN tmp.aida >> out.dat
+rm tmp.{aida,dat}
+
+flat2aida out.dat
+