tutorials/analysis/ana_simple.C File Reference

Go to the source code of this file.

Functions
void	ana_simple (TString infile, TString cfgfile="analysis.cfg", int nevts=0, TString outfile="")

Function Documentation

void ana_simple	(	TString	infile,
		TString	cfgfile = "analysis.cfg",
		int	nevts = 0,
		TString	outfile = ""
	)

Definition at line 10 of file tutorials/analysis/ana_simple.C.

References ctime, Double_t, gDebug, PndSimpleAnalysis::Init(), outfile, printf(), rootlogon(), rtime, and timer.

{
  TStopwatch timer;
  timer.Start();
  gDebug=0;

  // Load basic libraries
  gROOT->LoadMacro("$VMCWORKDIR/gconfig/rootlogon.C");rootlogon();

/*  // Load this example libraries
  gSystem->Load("libGeoBase");
  gSystem->Load("libParBase");
  gSystem->Load("libBase");
  //gSystem->Load("libPndData");
  //gSystem->Load("libGen");
  gSystem->Load("libRho");
  gSystem->Load("libfastsim");*/
 
  if (outfile=="") outfile=infile+"_ana.root";

  //analysis
  
  PndEventReader evr(infile);
  if (nevts==0) nevts=evr.GetEntries();
  
  cout <<"\n####### Processing "<<nevts <<" events...\n"<<endl;
        
  
  FairRunAna *fRunA= new FairRunAna();
  fRunA->SetInputFile(infile);
  fRunA->SetOutputFile(outfile);
  
  //append the analysis task
  PndSimpleAnalysis *anaTask=new PndSimpleAnalysis(cfgfile.Data());
  fRunA->AddTask(anaTask);

  fRunA->Init();
  fRunA->Run(0,nevts);
  
        
  timer.Stop();
  Double_t rtime = timer.RealTime();
  Double_t ctime = timer.CpuTime();
  printf("RealTime=%f seconds, CpuTime=%f seconds\n",rtime,ctime);
}