ATLAS Offline Software
Namespaces | Functions
testEgEfficiencyCorrFwd.cxx File Reference
#include "ElectronEfficiencyCorrection/AsgElectronEfficiencyCorrectionTool.h"
#include "SFHelpers.h"
#include <memory>
#include <TFile.h>
#include "xAODRootAccess/Init.h"
#include "xAODRootAccess/TEvent.h"
#include "xAODRootAccess/TStore.h"
#include "xAODEgamma/Electron.h"
#include "xAODEgamma/ElectronContainer.h"
#include "AsgMessaging/MessageCheck.h"
#include "AsgMessaging/MsgStream.h"
#include "AsgTools/StandaloneToolHandle.h"
#include "xAODRootAccess/tools/TFileAccessTracer.h"

Go to the source code of this file.

Namespaces

 asg
 

Functions

int main (int argc, char *argv[])
 

Function Documentation

◆ main()

int main ( int  argc,
char *  argv[] 
)

Definition at line 33 of file testEgEfficiencyCorrFwd.cxx.

34 {
35 
37  // The application's name:
38  const char* APP_NAME = argv[0];
39  using namespace asg::testEgEfficiencyCorrFwd;
40  ANA_CHECK_SET_TYPE(int);
41  MSG::Level mylevel = MSG::INFO;
42  setMsgLevel(mylevel);
43  msg().setName(APP_NAME);
44 
45  // Check if we received a file name:
46  if (argc < 2) {
47  ANA_MSG_ERROR(APP_NAME << "No file name received!");
49  APP_NAME
50  << " Usage: <<APP_NAME << [xAOD file name] [Num of events to use]");
51  return 1;
52  }
53 
54  // Initialise the application:
56 
57  // Open the input file:
58  const TString fileName = argv[1];
59  ANA_MSG_INFO("Opening file: " << fileName.Data());
60  std::unique_ptr<TFile> ifile(TFile::Open(fileName, "READ"));
61  ANA_CHECK(ifile.get());
62 
63  // Create a TEvent object:
65 
66  // Then the tools
67  std::vector<std::string> id_configFiles{
68  "ElectronEfficiencyCorrection/2012/offline/"
69  "efficiencySF.offline.FwdTight.2012.8TeV.rel17p2.GEO21.v02.root"
70  }; // we don't support keys for fwd electrons, yet. Our latest file is 2012,
71  // still
73  "AsgElectronEfficiencyCorrectionTool/ElEffCorrectionTool");
74  ANA_CHECK(
75  ElEffCorrectionTool.setProperty("CorrectionFileNameList", id_configFiles));
76  ANA_CHECK(ElEffCorrectionTool.setProperty("ForceDataType", 1));
77  ANA_CHECK(ElEffCorrectionTool.setProperty("OutputLevel", mylevel));
78  ANA_CHECK(ElEffCorrectionTool.setProperty("CorrelationModel", "FULL"));
79  ANA_CHECK(ElEffCorrectionTool.setProperty("UseRandomRunNumber", false));
80  ANA_CHECK(ElEffCorrectionTool.initialize());
81 
82  // Then open the file(s)
83  ANA_CHECK(event.readFrom(ifile.get()));
84  ANA_MSG_INFO("Number of available events to read in: "
85  << static_cast<long long int>(event.getEntries()));
86 
87  // Decide how many events to run over:
88  long long int entries = event.getEntries();
89  if (argc > 2) {
90  const long long int e = atoll(argv[2]);
91  if (e < entries) {
92  entries = e;
93  }
94  }
95  ANA_MSG_INFO("Number actual events to read in: " << entries);
96 
97  // Loop over the events:
98  for (long long int entry = 0; entry < entries; ++entry) {
99  event.getEntry(entry);
100  ANA_MSG_INFO(" \n ==> Event " << entry);
101 
102  const xAOD::ElectronContainer* electrons = nullptr;
103  ANA_CHECK(event.retrieve(electrons, "ForwardElectrons"));
104 
105  for (const xAOD::Electron* el : *electrons) {
106  if (el->pt() < 20000)
107  continue; // skip electrons outside of recommendations
108  if (fabs(el->caloCluster()->eta()) < 2.5)
109  continue; // skip electrons outside of recommendations
110  int index = ElEffCorrectionTool->systUncorrVariationIndex(*el);
111  /*
112  * Set up the systematic variations
113  */
114  bool isToys = false;
115  double nominalSF{};
116  double totalNeg{};
117  double totalPos{};
118  ANA_CHECK(
120  ElEffCorrectionTool, *el, nominalSF, totalPos, totalNeg, isToys) ==
121  0);
122 
123  ANA_MSG_INFO("===> electron : Pt = "
124  << el->pt() << " : eta = " << el->eta()
125  << " : Bin index = " << index << " : SF = " << nominalSF
126  << " + " << totalPos << " - " << totalNeg << " <===");
127  }
128  }
129 
130  ANA_MSG_INFO("===> DONE <===\n");
131  return 0;
132 }
AllowedVariables::e
e
Definition: AsgElectronSelectorTool.cxx:37
index
Definition: index.py:1
SFHelpers::result
int result(asg::StandaloneToolHandle< IAsgElectronEfficiencyCorrectionTool > &tool, const xAOD::Electron &el, double &nominalSF, double &totalPos, double &totalNeg, const bool isToys)
Definition: SFHelpers.cxx:20
xAOD::TFileAccessTracer::enableDataSubmission
static void enableDataSubmission(::Bool_t value)
Function for turning data submission on/off.
Definition: TFileAccessTracer.cxx:281
ANA_MSG_ERROR
#define ANA_MSG_ERROR(xmsg)
Macro printing error messages.
Definition: Control/AthToolSupport/AsgMessaging/AsgMessaging/MessageCheck.h:294
ANA_CHECK
#define ANA_CHECK(EXP)
check whether the given expression was successful
Definition: Control/AthToolSupport/AsgMessaging/AsgMessaging/MessageCheck.h:324
POOL::TEvent::readFrom
StatusCode readFrom(TFile *file)
Definition: PhysicsAnalysis/POOLRootAccess/src/TEvent.cxx:133
FortranAlgorithmOptions.fileName
fileName
Definition: FortranAlgorithmOptions.py:13
TrigConf::MSGTC::Level
Level
Definition: Trigger/TrigConfiguration/TrigConfBase/TrigConfBase/MsgStream.h:21
event
POOL::TEvent event(POOL::TEvent::kClassAccess)
POOL::TEvent::getEntries
long getEntries()
Definition: PhysicsAnalysis/POOLRootAccess/src/TEvent.cxx:124
asg::StandaloneToolHandle
an "initializing" ToolHandle for stand-alone applications
Definition: StandaloneToolHandle.h:44
LArCellNtuple.argv
argv
Definition: LArCellNtuple.py:152
ANA_MSG_INFO
#define ANA_MSG_INFO(xmsg)
Macro printing info messages.
Definition: Control/AthToolSupport/AsgMessaging/AsgMessaging/MessageCheck.h:290
APP_NAME
#define APP_NAME
Definition: BoostedXbbTag.cxx:23
plotIsoValidation.el
el
Definition: plotIsoValidation.py:197
DQHistogramMergeRegExp.argc
argc
Definition: DQHistogramMergeRegExp.py:20
DataVector
Derived DataVector<T>.
Definition: DataVector.h:794
GetAllXsec.entry
list entry
Definition: GetAllXsec.py:132
xAOD::Electron_v1
Definition: Electron_v1.h:34
ANA_CHECK_SET_TYPE
#define ANA_CHECK_SET_TYPE(TYPE)
set the type for ANA_CHECK to report failures
Definition: Control/AthToolSupport/AsgMessaging/AsgMessaging/MessageCheck.h:314
POOL::TEvent::retrieve
StatusCode retrieve(const T *&obj)
Definition: PhysicsAnalysis/POOLRootAccess/POOLRootAccess/TEvent.h:74
entries
double entries
Definition: listroot.cxx:49
LArCellNtuple.ifile
string ifile
Definition: LArCellNtuple.py:133
python.AutoConfigFlags.msg
msg
Definition: AutoConfigFlags.py:7
InDetDD::electrons
@ electrons
Definition: InDetDD_Defs.h:17
xAOD::TEvent
Tool for accessing xAOD files outside of Athena.
Definition: Control/xAODRootAccess/xAODRootAccess/TEvent.h:84
xAOD::Init
StatusCode Init(const char *appname)
Function initialising ROOT/PyROOT for using the ATLAS EDM.
Definition: Init.cxx:31