ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitHzz4l/ElectronStudies/src/runPhotonEnergyRegressionTreeWriter.cc
Revision: 1.1
Committed: Thu Jul 26 22:49:16 2012 UTC (12 years, 9 months ago) by anlevin
Content type: text/plain
Branch: MAIN
CVS Tags: compiled, HEAD
Log Message:
added stuff for photon energy regression studies

File Contents

# Content
1 //./bin/runPhotonEnergyRegressionTreeWriter.exe --catalogDir "/home/cmsprod/catalog" --book "t2mit/filefi/028" --dataset "s12-zllm50-2-v9" --fileset "0006" --output_filename "delete.root" --n_events_to_process 1000 --rho rho_kt_6_pf_jets
2
3 #include "PhotonEnergyRegressionTreeWriterMod.h"
4 #include "ElectronEnergyRegressionTreeWriterArgs.h"
5 #include "TChain.h"
6
7 int main(int argc, char** argv)
8 {
9 ElectronEnergyRegressionTreeWriterFlags ctrl;
10
11 parse_electron_energy_regression_tree_writer_args( argc, argv, ctrl );
12 if( ctrl.catalog_dir.empty() || ctrl.book.empty() || ctrl.dataset.empty()|| ctrl.fileset.empty()|| ctrl.output_filename.empty()|| ctrl.n_events_to_process == -2 || ctrl.rho.empty())
13 {
14 cerr << "usage: runPhotonEnergyRegressionTreeWriter.exe <flags> " << endl;
15 ctrl.dump();
16
17 return 1;
18 }
19 ctrl.dump();
20 assert(ctrl.rho == "rho" || ctrl.rho == "rho_kt_6_pf_jets");
21
22
23 mithep::PhotonEnergyRegressionTreeWriterMod *pertwMod = new mithep::PhotonEnergyRegressionTreeWriterMod;
24 mithep::Analysis *ana = new mithep::Analysis;
25 if(ctrl.n_events_to_process == -1)
26 ana->SetProcessNEvents(TChain::kBigNumber);
27 else
28 ana->SetProcessNEvents(ctrl.n_events_to_process);
29 ana->AddSuperModule(pertwMod);
30
31 if(ctrl.rho == "rho_kt_6_pf_jets")
32 pertwMod->rho_version = mithep::PhotonEnergyRegressionTreeWriterMod::e_rho_kt_6_pf_jets;
33 else if (ctrl.rho == "rho")
34 pertwMod->rho_version = mithep::PhotonEnergyRegressionTreeWriterMod::e_rho;
35 else
36 assert(0);
37
38 mithep::Catalog *c = new mithep::Catalog(TString(ctrl.catalog_dir));
39 mithep::Dataset *d = c->FindDataset(TString(ctrl.book),TString(ctrl.dataset),TString(ctrl.fileset));
40 ana->AddDataset(d);
41 //ana->AddFile("root://castorcms//castor/cern.ch/user/p/paus/filefi/028/r12a-dmu-pr-v1/FEF8A575-A785-E111-A21A-003048D37538.root");
42 pertwMod->fOutputName = ctrl.output_filename;
43 ana->Run(kFALSE);
44
45
46 }
47