Changeset 03d2c3f in git for validation

Timestamp:

Apr 13, 2017, 1:30:01 PM (8 years ago)

Author:

Pavel Demin <pavel-demin@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

54ba29e

Parents:

65dc4e7

Message:

move validation scripts to validation directory

Location:

validation

Files:

• FlatGunLHEventProducer.py (added)
• Validation.cpp (moved) (moved from examples/Validation.cpp )
• validation.sh (moved) (moved from examples/validation.sh ) (6 diffs)

validation/validation.sh

@@ -1 +1 @@
 #!/bin/sh
-##############################################################################################
+################################################################################
 #
 # This code produces at set of validation plots for a given detector card. 
@@ -19 +19 @@
 # This said, 500k events should be ok for most cases. 
 #
-############################################################################################
+################################################################################
 
 EXPECTED_ARGS=2
@@ -32 +32 @@
 
 cardbase=$(basename $1)
-carddir=$(dirname $1)
 nEvents=$2
-validationCard=$carddir/validation_$cardbase
 output=validation_${cardbase%.*}.root
 mainoutputdir=report_${cardbase%.*}
-outputroot=report_${cardbase%.*}/root
+outputrootdir=report_${cardbase%.*}/root
 cardlabel=${cardbase%.*}
 version=$(cat VERSION)
 outpdf=$mainoutputdir/${output%.*}.pdf
-qcdPgLocation="/home/fynu/mselvaggi/storage/DelphesValidationSamples"
+cardsdir=validation/cards
+samplesdir=validation/samples
+validationcard=$cardsdir/validation_$cardbase
 
-mkdir -p $outputroot
+mkdir -p $cardsdir
+mkdir -p $samplesdir
+mkdir -p $outputrootdir
 mkdir -p $mainoutputdir/www/fig
 
-sed 's/delphes_card_CMS.tcl/'$cardbase'/g' cards/validation_card.tcl  > $validationCard
-sed -i "1i set MaxEvents ${nEvents}" "$validationCard"
+sed "s/delphes_card_CMS.tcl/$cardbase/g" cards/validation_card.tcl > $validationcard
+sed -i "1i set MaxEvents $nEvents" $validationcard
 
 function runParticleGun {
   name=$1
   pid=$2
-  cmnd="examples/Pythia8/configParticleGun_$name.cmnd"
-  rootfile="particleGun_${name}_${cardlabel}.root"
-  sed '/Main:spareMode1/s/=[[:space:]]*[0-9]*/= '$pid'/' examples/Pythia8/configParticleGun.cmnd > examples/Pythia8/tmp.cmnd
-  sed '/Main:numberOfEvents/s/=[[:space:]]*[0-9]*/= '$nEvents'/' examples/Pythia8/tmp.cmnd  > $cmnd 
-  ./DelphesPythia8 $validationCard $cmnd $outputroot/$rootfile
-  
+  cmnd=$cardsdir/configParticleGun_$name.cmnd
+  outputroot=particleGun_${name}_${cardlabel}.root
+  sed "/Main:spareMode1/s/=[[:space:]]*[0-9]*/= $pid/; /Main:numberOfEvents/s/=[[:space:]]*[0-9]*/= $nEvents/" examples/Pythia8/configParticleGun.cmnd > $cmnd
+  ./DelphesPythia8 $validationcard $cmnd $outputrootdir/$outputroot
 }
 
@@ -64 +64 @@
   name=$1
   pid=$2
-  inputroot="${qcdPgLocation}/${pid}.root"
-  rootfile="particleGun_${name}_${cardlabel}.root"
-  ./DelphesROOT $validationCard $outputroot/$rootfile $inputroot
-  echo "./DelphesROOT $validationCard $outputroot/$rootfile $inputroot"
+  cmnd=$cardsdir/configLHE_$pid.cmnd
+  inputroot=$samplesdir/$pid.root
+  outputroot=particleGun_${name}_${cardlabel}.root
+
+  if [ ! -f $inputroot ]
+  then
+    python FlatGunLHEventProducer.py --pdg $pid --ptmin 1 --ptmax 50000 --etamin -6 --etamax 6 --size $nEvents --seed 1 --output $samplesdir/events_$pid.lhe --log --ecm 100000
+
+    cp examples/Pythia8/configLHE.cmnd $cmnd
+    echo "Beams:LHEF = $samplesdir/events_$pid.lhe" >> $cmnd
+    echo "Main:numberOfEvents = $nEvents" >> $cmnd
+
+    ./DelphesPythia8 cards/gen_card.tcl $cmnd $inputroot
+  fi
+
+  ./DelphesROOT $validationcard $outputrootdir/$outputroot $inputroot
 }
 
@@ -85 +97 @@
 
 
-./Validation $outputroot/particleGun_pion_$cardlabel.root $outputroot/particleGun_electron_$cardlabel.root $outputroot/particleGun_muon_$cardlabel.root $outputroot/particleGun_photon_$cardlabel.root $outputroot/particleGun_neutron_$cardlabel.root $outputroot/particleGun_jet_$cardlabel.root $outputroot/particleGun_bjet_$cardlabel.root $outputroot/particleGun_cjet_$cardlabel.root $outputroot/particleGun_taujet_$cardlabel.root $mainoutputdir/$output $version
+./Validation $outputrootdir/particleGun_pion_$cardlabel.root $outputrootdir/particleGun_electron_$cardlabel.root $outputrootdir/particleGun_muon_$cardlabel.root $outputrootdir/particleGun_photon_$cardlabel.root $outputrootdir/particleGun_neutron_$cardlabel.root $outputrootdir/particleGun_jet_$cardlabel.root $outputrootdir/particleGun_bjet_$cardlabel.root $outputrootdir/particleGun_cjet_$cardlabel.root $outputrootdir/particleGun_taujet_$cardlabel.root $mainoutputdir/$output $version
 
 
@@ -98 +110 @@
 mv HCal.png $mainoutputdir/www/fig/img_hcal.png
 
-
-
-