improved way to replace input file in ED filter

Author: diemort

run_fpmc.sh

@@ -2,7 +2,7 @@
 cd fpmc/
 source setup_lxplus.sh
 ./fpmc-hepmc \
-	--cfg Datacards/dataQED_WW \
+	--cfg Datacards/dataDPE_Dijets \
 	--comenergy 13000 \
-	--fileout dataWW.hepmc \
-	--nevents 10
+	--fileout data_dijets.hepmc \
+	--nevents 100

run_sim.sh

@@ -1,18 +1,28 @@
 
 # Details here: https://twiki.cern.ch/twiki/bin/view/Sandbox/PpsCepStudies#FPMC
+
+# variable JOBNAME is been used for IO filenames from simulation and event sample used as input:
+#
+# $JOBNAME.hepmc: input event sample from FPMC
+# $JOBNAME.py: config file for cmsRun
+# $JOBNAME.root: output file from simulation
+
 CMSSW=CMSSW_8_0_1
 AREA=$PWD/$CMSSW/src/
-OUTPUTFILE=dataWW.py
-INPUTFILE=dataWW.hepmc
+JOBNAME=data_dijets
+NEVENTS=100
 
-cp fpmc/$INPUTFILE $AREA
+echo "Copying event sample"
+cp fpmc/$JOBNAME.hepmc $AREA
 
 cd $AREA
+echo "Setting up CMSSW environment"
 eval `scramv1 runtime -sh`
 
+echo "Creating configuration file from cmsDriver"
 cmsDriver.py \
 	readHepMC_cff.py \
-	-n 10 \
+	-n $((NEVENTS)) \
 	--fast \
 	--conditions auto:run2_mc \
 	--eventcontent AODSIM \
@@ -23,29 +33,34 @@ cmsDriver.py \
 	--era Run2_25ns \
 	--customise FastSimulation/PPSFastSim/customise_FastSimCTPPS_cff.customise \
 	--no_exec \
-	--fileout=dataWW.root \
-	--python_filename=$OUTPUTFILE
+	--fileout="$JOBNAME".root \
+	--python_filename="$JOBNAME".py
 
 # Replace input file:
 
-sed -i -e 's/FPMC_WW_Inclusive_13TeV/dataWW/g' Configuration/Generator/python/readHepMC_cff.py
+echo "Replacing input filename for hepMC format"
+sed -i 's/file:.*.hepmc/file:'"$JOBNAME"'.hepmc/g' Configuration/Generator/python/readHepMC_cff.py
 
 # modify the resulting driver (in this examples, 'readHepMC_cff_py_GEN_SIM_RECOBEFMIX_DIGI_RECO.py')
 # by adding the following line, with 'source' as the input HepMC sample:
 # process.VtxSmeared.src = 'source'
 
-sed '94iprocess.VtxSmeared.src = "'$INPUTFILE'"' $OUTPUTFILE > test.py
-mv test.py $OUTPUTFILE
+echo "Adding hepMC event sample for simulation"
+sed '94iprocess.VtxSmeared.src = '"\"$JOBNAME"'.hepmc\"' "$JOBNAME".py > test.py
+mv test.py $JOBNAME.py
 
 # filter has to be modified to include specific channels; for WW see:
 # https://raw.githubusercontent.com/uerj-cms-cep-studies/tmp/master/inputHepMC.cc
+echo "Fetching ED filter"
 mkdir inputHepMC
 cd inputHepMC
 mkedfltr inputHepMC
 cd ..
 
 # when ready, run:
+echo "Compiling CMSSW area"
 scram b -j 8
 # and:
-cmsRun dataWW.py
+echo "Running job"
+cmsRun $JOBNAME.py
 # or equivalent driver.