![Logo](../../ATLAS-Logo-Square-Blue-RGB.png) |
ATLAS Offline Software
|
Go to the documentation of this file.
42 TString
fileName =
"/eos/atlas/atlascerngroupdisk/perf-jets/ReferenceFiles/mc16_13TeV.361028.Pythia8EvtGen_A14NNPDF23LO_jetjet_JZ8W.deriv.DAOD_FTAG1.e3569_s3126_r9364_r9315_p3260/DAOD_FTAG1.12133096._000074.pool.root.1";
48 Info(
APP_NAME,
"==============================================" );
52 Info(
APP_NAME,
" $> %s -n X | X = number of events you want to run on",
APP_NAME );
53 Info(
APP_NAME,
" $> %s -e X | X = specific number of the event to run on - for debugging",
APP_NAME );
55 Info(
APP_NAME,
"==============================================" );
59 Info(
APP_NAME,
"No arguments - using default file" );
71 if(
options.find(
"-f")!=std::string::npos){
72 for(
int ipos=0; ipos<
argc ; ipos++ ) {
73 if(std::string(
argv[ipos]).compare(
"-f")==0){
81 if(
options.find(
"-event")!=std::string::npos){
82 for(
int ipos=0; ipos<
argc ; ipos++ ) {
85 Info(
APP_NAME,
"Argument (-event) : Running only on event # %i", ievent );
91 if(
options.find(
"-n")!=std::string::npos){
92 for(
int ipos=0; ipos<
argc ; ipos++ ) {
101 if(
options.find(
"-v")!=std::string::npos){
103 Info(
APP_NAME,
"Argument (-v) : Setting verbose");
111 StatusCode::enableFailure();
117 TChain *
chain =
new TChain (
"CollectionTree",
"CollectionTree");
128 Long64_t
entries =
event.getEntries();
131 TFile*
outputFile = TFile::Open(
"output_JetQGTagger.root",
"recreate" );
133 TTree*
Tree =
new TTree(
"tree",
"test_tree" );
134 Tree->Branch(
"pass", &pass,
"pass/I" );
145 std::cout<<
"Initializing QG Tagger"<<std::endl;
149 ANA_CHECK( m_Tagger.
setProperty(
"ConfigFile",
"SmoothedWZTaggers/SmoothedContainedWTagger_AntiKt10LCTopoTrimmed_FixedSignalEfficiency50_MC15c_20161215.dat") );
159 event.getEntry(
entry );
163 if(
event.
retrieve( evtInfo,
"EventInfo" ) != StatusCode::SUCCESS){
166 if(ievent!=-1 &&
static_cast <int> (evtInfo->
eventNumber())!=ievent) {
172 if(
event.
retrieve( myJets,
"AntiKt10LCTopoTrimmedPtFrac5SmallR20Jets" ) != StatusCode::SUCCESS)
180 std::cout <<
"Testing W Tagger " << std::endl;
181 std::cout <<
"RunningTag : " <<
jet->auxdata<
bool>(
"Tagged") << std::endl;
184 pass =
jet->auxdata<
bool>(
"Tagged");
189 Info(
APP_NAME,
"===>>> done processing event #%i, run #%i %i events processed so far <<<===",
static_cast< int >( evtInfo->
eventNumber() ),
static_cast< int >( evtInfo->
runNumber() ),
static_cast< int >(
entry + 1 ) );
205 std::cout<<std::endl<<std::endl;
206 std::cout<<
"Smart Slimming Checker :"<<std::endl;
208 std::cout<<std::endl<<std::endl;
std::pair< int, int > compare(const AmgSymMatrix(N) &m1, const AmgSymMatrix(N) &m2, double precision=1e-9, bool relative=false)
compare two matrices, returns the indices of the first element that fails the condition,...
uint64_t eventNumber() const
The current event's event number.
ReadStats & stats()
Access the object belonging to the current thread.
@ kAthenaAccess
Access containers/objects like Athena does.
uint32_t runNumber() const
The current event's run number.
int main(int argc, char *argv[])
POOL::TEvent event(POOL::TEvent::kClassAccess)
static IOStats & instance()
Singleton object accessor.
void printSmartSlimmingBranchList(bool autoIncludeLinks=false) const
Print the accessed variables, formatted for smart slimming.
virtual StatusCode tag(const xAOD::Jet &jet, const xAOD::Vertex *pv) const override
Class describing the basic event information.
A relatively simple transient store for objects created in analysis.
StatusCode retrieve(const T *&obj)
int atoi(std::string_view str)
Helper functions to unpack numbers decoded in string into integers and doubles The strings are requir...
Error
The different types of error that can be flagged in the L1TopoRDO.
Tool for accessing xAOD files outside of Athena.
StatusCode Init(const char *appname)
Function initialising ROOT/PyROOT for using the ATLAS EDM.