run_HypSimpleAna_Geo2_Bg.C File Reference

Go to the source code of this file.

Functions
int	run_HypSimpleAna_Geo2 (TString base="Hit_hypIdReco", int nevts=0, TString outfile="")

Function Documentation

int run_HypSimpleAna_Geo2	(	TString	base = "Hit_hypIdReco",
		int	nevts = 0,
		TString	outfile = ""
	)

Definition at line 10 of file run_HypSimpleAna_Geo2_Bg.C.

References ctime, Double_t, PndHypSimpleAna::Finish(), gDebug, infile(), nEvents, outfile, parFile, printf(), rootlogon(), rtime, timer, and TString.

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

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

 gSystem->Load("librazhyp");
 gSystem->Load("libHyp");

 gSystem->Load("libRho");
 

  //if (outfile=="") outfile=infile+"_ana.root";

  //analysis
  
  
  //aida second step 0.5 T
   TString simfile = "../Alicia/Bup_hypFSG41TXm.root";//Bup_hypFSG41TXXb_Geo2.root";
       Int_t nEvents = 0;//50000;
       TString parFile="../Alicia/Bup_hypFS1TG4Xmparams.root";//Bup_hypFS1TG4XXb_Geo2params.root";
      TString recofile = "../Alicia/hypIdealrecoXm_Geo2.root";
      TString infile = "../Alicia/MicroSimpleG4Bup1TXm_Geo2.root";//MicroSimpleG4Bup1TXXbP8_Geo2.root";
      if (outfile=="")  outfile="AnaFSimpleBup1TXm_Bg.root";//AnaFSimpleP8Bup1TXXb.root";


 
  
 


  FairRunAna *fRunA= new FairRunAna();
   
 
   fRunA->SetInputFile(simfile.Data());
   fRunA->AddFriend(infile);
   fRunA->AddFriend(recofile);
   // fRunA->AddFile(simfile2.Data());
   // fRunA->AddFriend(infile);
   // fRunA->AddFriend(recofile);
  
   
   fRunA->SetOutputFile(outfile.Data());
  
  //append the analysis task
  //PndMicroWriter *tcandTask=new PndMicroWriter();
  //fRunAA->AddTask(tcandTask);
  
  PndHypSimpleAna *anaTask=new PndHypSimpleAna();
  fRunA->AddTask(anaTask);
  
  fRunA->Init();
  //if (nevts==0) nevts=10;
  fRunA->Run(0,nevts);
  anaTask->Finish(base.Data());
  
        
  timer.Stop();
  Double_t rtime = timer.RealTime();
  Double_t ctime = timer.CpuTime();
  printf("RealTime=%f seconds, CpuTime=%f seconds\n",rtime,ctime);
  return 0;
}