Changeset fa068d3 in git
- Timestamp:
- Sep 6, 2016, 11:21:11 AM (8 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- e9971a7
- Parents:
- b9ae4c3
- Location:
- examples
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
examples/Validation.cpp
rb9ae4c3 rfa068d3 551 551 552 552 TF1 *f2 = new TF1("f2", "gaus", f1->GetParameter(1) - 2*f1->GetParameter(2), f1->GetParameter(1) + 2*f1->GetParameter(2)); 553 hist->Fit("f2","RQ"); 553 hist->Fit("f2","RQ"); 554 554 555 555 Double_t sig = f2->GetParameter(2); … … 914 914 DrawAxis(mg_elFwd, leg_elFwd, 0.2); 915 915 916 C_el1->Print("delphes_validation.pdf(","pdf"); 917 C_el2->Print("delphes_validation.pdf","pdf"); 916 TString pdfOutput(outputFile); 917 pdfOutput.ReplaceAll(".root", ".pdf"); 918 919 C_el1->Print(pdfOutput+"(","pdf"); 920 C_el2->Print(pdfOutput,"pdf"); 918 921 919 922 gDirectory->cd(0); … … 1015 1018 DrawAxis(mg_muFwd, leg_muFwd, 0.3, 1); 1016 1019 1017 C_mu1->Print( "delphes_validation.pdf","pdf");1018 C_mu->Print( "delphes_validation.pdf","pdf");1020 C_mu1->Print(pdfOutput,"pdf"); 1021 C_mu->Print(pdfOutput,"pdf"); 1019 1022 1020 1023 gDirectory->cd(0); … … 1116 1119 DrawAxis(mg_phFwd, leg_phFwd, 0.1); 1117 1120 1118 C_ph1->Print( "delphes_validation.pdf","pdf");1119 C_ph->Print( "delphes_validation.pdf","pdf");1121 C_ph1->Print(pdfOutput,"pdf"); 1122 C_ph->Print(pdfOutput,"pdf"); 1120 1123 1121 1124 gDirectory->cd(0); … … 1267 1270 mg_met->GetYaxis()->SetTitle("#sigma(ME_{x}) [GeV]"); 1268 1271 1269 C_jet->Print( "delphes_validation.pdf","pdf");1270 C_btag1->Print( "delphes_validation.pdf","pdf");1271 C_tautag1->Print( "delphes_validation.pdf","pdf");1272 C_met->Print( "delphes_validation.pdf)","pdf");1272 C_jet->Print(pdfOutput,"pdf"); 1273 C_btag1->Print(pdfOutput,"pdf"); 1274 C_tautag1->Print(pdfOutput,"pdf"); 1275 C_met->Print(pdfOutput+")","pdf"); 1273 1276 1274 1277 TFile *fout = new TFile(outputFile,"recreate"); -
examples/validation.sh
rb9ae4c3 rfa068d3 1 ############################################################################################## 2 # 3 # This code produces at set of validation plots for a given detector card. 4 # 5 # In order to run this you need to compile Delphes with Pythia8 first, see: 6 # 7 # https://cp3.irmp.ucl.ac.be/projects/delphes/wiki/WorkBook/Pythia8 8 # 9 # After you (re-)compiled Delphes with Pythia8 you are ready to go, execute from Delphes main dir: 10 # 11 # ./examples/validation.sh [detector_card] [number_of_events] 12 # 13 # e.g. 14 # 15 # ./examples/validation.sh delphes_card_CMS.tcl 100000 16 # 17 # Note that the more events you specify, the more accurate the controls plots will be ... 18 # 19 ############################################################################################ 20 1 21 #! /bin/sh 22 23 EXPECTED_ARGS=2 24 E_BADARGS=65 25 26 if [ $# -ne $EXPECTED_ARGS ] 27 then 28 echo "Usage: ./examples/validation.sh [detector_card] [number_of_events]" 29 echo "for instance: ./examples/validation.sh delphes_card_CMS.tcl 10000" 30 exit $E_BADARGS 31 fi 32 33 card=$1 34 nEvents=$2 35 validationCard=cards/validation_$card 36 output=validation_${card%.*}.root 37 38 sed 's/delphes_card_CMS.tcl/'$card'/g' cards/validation_card.tcl > $validationCard 39 2 40 function runParticleGun { 3 41 name=$1 … … 5 43 cmnd="examples/Pythia8/configParticleGun_$name.cmnd" 6 44 sed '/Main:spareMode1/s/=[[:space:]]*[0-9]*/= '$pid'/' examples/Pythia8/configParticleGun.cmnd > $cmnd 7 ./DelphesPythia8 cards/validation_card.tcl $cmnd delphes_ParticleGun_$name.root 45 sed '/Main:numberOfEvents/s/=[[:space:]]*[0-9]*/= '$nEvents'/' examples/Pythia8/configParticleGun.cmnd > $cmnd 46 ./DelphesPythia8 $validationCard $cmnd delphes_ParticleGun_$name.root 8 47 } 9 48 … … 15 54 runParticleGun taujet 15 16 55 17 ./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.root56 ./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
Note:
See TracChangeset
for help on using the changeset viewer.