解析コード作成

// Infrastructure include(s):

#include "xAODRootAccess/Init.h"

#include "xAODRootAccess/TEvent.h"

// ASG status code check

#include <AsgTools/MessageCheck.h>

job.useXAOD ();

ANA_CHECK(xAOD::Init());

xAOD::TEvent* event = wk()->xaodEvent();

Info("initialize()", "Number of events = %lli", event->getEntries() );

PACKAGE_DEP = EventLoop xAODRootAccess AsgTools

#include "xAODRootAccess/Init.h"

#include "SampleHandler/SampleHandler.h"

#include "SampleHandler/ScanDir.h"

#include "SampleHandler/ToolsDiscovery.h"

#include "EventLoop/Job.h"

#include "EventLoop/DirectDriver.h"

#include "SampleHandler/DiskListLocal.h"

#include "MyNewPackage/MyAnalysisAlg.h"

#include <TSystem.h>

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

// Take the submit directory from the input if provided:

std::string submitDir = "submitDir"; if( argc > 1 ) submitDir = argv[ 1 ];

// Set up the job for xAOD access:

xAOD::Init().ignore();

// Construct the samples to run on:

SH::SampleHandler sh;

// use SampleHandler to scan all of the subdirectories of a directory for particular MC single file:

const char* inputFilePath = gSystem->ExpandPathName ("/gpfs/fs2001/nobe/data2086b/sample_forTutorial2016/data16_13TeV.00310247.physics_Main.merge.DAOD_HIGG2D4.f755_m1699_p2880/");

SH::ScanDir().filePattern("*pool.root*").scan(sh,inputFilePath);

// Set the name of the input TTree. It's always "CollectionTree" // for xAOD files.

sh.setMetaString( "nc_tree", "CollectionTree" );

// Print what we found:

sh.print();

// Create an EventLoop job:

EL::Job job; job.sampleHandler( sh );

job.options()->setDouble (EL::Job::optMaxEvents, 500); // maximum number of events to be analyzed. If you want to see all events, please set it to -1

// Add our analysis to the job:

MyAnalysisAlg* alg = new MyAnalysisAlg (); job.algsAdd( alg );

// Run the job using the local/direct driver:

EL::DirectDriver driver; driver.submit( job, submitDir );

return 0;

}