Changeset fa068d3 in git for examples

Timestamp:

Sep 6, 2016, 11:21:11 AM (8 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

e9971a7

Parents:

b9ae4c3

Message:

added card and number_of_events parameters, and some instructions

Location:

examples

Files:

• Validation.cpp (modified) (5 diffs)
• validation.sh (modified) (3 diffs)

examples/Validation.cpp

@@ -551 +551 @@
 
   TF1 *f2 = new TF1("f2", "gaus", f1->GetParameter(1) - 2*f1->GetParameter(2), f1->GetParameter(1) + 2*f1->GetParameter(2));
-  hist->Fit("f2","RQ");  
+  hist->Fit("f2","RQ");
 
   Double_t sig = f2->GetParameter(2);
@@ -914 +914 @@
   DrawAxis(mg_elFwd, leg_elFwd, 0.2);
 
-  C_el1->Print("delphes_validation.pdf(","pdf");
-  C_el2->Print("delphes_validation.pdf","pdf");
+  TString pdfOutput(outputFile);
+  pdfOutput.ReplaceAll(".root", ".pdf");
+
+  C_el1->Print(pdfOutput+"(","pdf");
+  C_el2->Print(pdfOutput,"pdf");
 
   gDirectory->cd(0);
@@ -1015 +1018 @@
   DrawAxis(mg_muFwd, leg_muFwd, 0.3, 1);
 
-  C_mu1->Print("delphes_validation.pdf","pdf");
-  C_mu->Print("delphes_validation.pdf","pdf");
+  C_mu1->Print(pdfOutput,"pdf");
+  C_mu->Print(pdfOutput,"pdf");
 
   gDirectory->cd(0);
@@ -1116 +1119 @@
   DrawAxis(mg_phFwd, leg_phFwd, 0.1);
 
-  C_ph1->Print("delphes_validation.pdf","pdf");
-  C_ph->Print("delphes_validation.pdf","pdf");
+  C_ph1->Print(pdfOutput,"pdf");
+  C_ph->Print(pdfOutput,"pdf");
 
   gDirectory->cd(0);
@@ -1267 +1270 @@
   mg_met->GetYaxis()->SetTitle("#sigma(ME_{x}) [GeV]");
 
-  C_jet->Print("delphes_validation.pdf","pdf");
-  C_btag1->Print("delphes_validation.pdf","pdf");
-  C_tautag1->Print("delphes_validation.pdf","pdf");
-  C_met->Print("delphes_validation.pdf)","pdf");
+  C_jet->Print(pdfOutput,"pdf");
+  C_btag1->Print(pdfOutput,"pdf");
+  C_tautag1->Print(pdfOutput,"pdf");
+  C_met->Print(pdfOutput+")","pdf");
 
   TFile *fout = new TFile(outputFile,"recreate");

examples/validation.sh

@@ +1 @@
+##############################################################################################
+#
+# This code produces at set of validation plots for a given detector card. 
+#
+# In order to run this you need to compile Delphes with Pythia8 first, see: 
+# 
+# https://cp3.irmp.ucl.ac.be/projects/delphes/wiki/WorkBook/Pythia8
+#
+# After you (re-)compiled Delphes with Pythia8 you are ready to go, execute from Delphes main dir:
+#
+# ./examples/validation.sh [detector_card] [number_of_events] 
+# 
+#  e.g.
+# 
+# ./examples/validation.sh delphes_card_CMS.tcl 100000
+#
+# Note that the more events you specify, the more accurate the controls plots will be ...
+#
+############################################################################################
+
 #! /bin/sh
+
+EXPECTED_ARGS=2
+E_BADARGS=65
+
+if [ $# -ne $EXPECTED_ARGS ]
+then
+  echo "Usage: ./examples/validation.sh [detector_card] [number_of_events]"
+  echo "for instance: ./examples/validation.sh delphes_card_CMS.tcl 10000"
+  exit $E_BADARGS
+fi
+
+card=$1
+nEvents=$2
+validationCard=cards/validation_$card
+output=validation_${card%.*}.root
+
+sed 's/delphes_card_CMS.tcl/'$card'/g' cards/validation_card.tcl  > $validationCard
+
 function runParticleGun {
   name=$1
@@ -5 +43 @@
   cmnd="examples/Pythia8/configParticleGun_$name.cmnd"
   sed '/Main:spareMode1/s/=[[:space:]]*[0-9]*/= '$pid'/' examples/Pythia8/configParticleGun.cmnd > $cmnd
-  ./DelphesPythia8 cards/validation_card.tcl $cmnd delphes_ParticleGun_$name.root
+  sed '/Main:numberOfEvents/s/=[[:space:]]*[0-9]*/= '$nEvents'/' examples/Pythia8/configParticleGun.cmnd  > $cmnd 
+  ./DelphesPythia8 $validationCard $cmnd delphes_ParticleGun_$name.root
 }
 
@@ -15 +54 @@
 runParticleGun taujet 15
 
-./Validation delphes_ParticleGun_electron.root delphes_ParticleGun_muon.root delphes_ParticleGun_photon.root delphes_ParticleGun_jet.root delphes_ParticleGun_bjet.root delphes_ParticleGun_taujet.root delphes_validation.root
+./Validation delphes_ParticleGun_electron.root delphes_ParticleGun_muon.root delphes_ParticleGun_photon.root delphes_ParticleGun_jet.root delphes_ParticleGun_bjet.root delphes_ParticleGun_taujet.root $output