realnewrun.cxx

Go to the documentation of this file.

#include <iostream>
#include <fstream>
#include <iomanip>
#include <cstring>
#include <string>
#include <vector>
#include <TFile.h>
#include <TTree.h>
#include <TString.h>
#include <TBranch.h>
#include <TLeaf.h>
#include <TH1F.h>
#include <cstdio>
 
using namespace std;
 
int main(int argc, char *argv[]){
  if(argc<5){
    cout << "please append two root files' names you want to compare " << endl;
    cout << "realnewrun.exe new.root treename old.root treename [output.root]" << endl;
    return 1;
  }
 
 
  TFile f1(argv[1]);
  TTree *t1 = (TTree*)f1.Get(argv[2]);
  TFile f2(argv[3]);
  TTree *t2 = (TTree*)f2.Get(argv[4]);
  TString str_output = "output.root";
  if(argc>5) str_output = argv[5];
  TFile f3(str_output,"recreate");
  TTree *t3 = new TTree("track","track");
  int newrun, newevt;
  t3->Branch("run", &newrun, "run/I");
  t3->Branch("event", &newevt, "event/I");
 
  int run, evt;
  t1->SetBranchAddress("run", &run);
  t1->SetBranchAddress("event", &evt);
  char str[255];
  for(int i=0; i<t1->GetEntries(); i++){
    t1->GetEntry(i);
    sprintf(str, "run==%d && event==%d", run, evt);
    if(t2->GetEntries(str)>0) cout << "str: " << str << "  run: " << run << "  event: " << evt << endl;
    else{
      newrun = run;
      newevt = evt;
      t3->Fill();
    }
  }
 
  f1.Close();
  f2.Close();
  f3.cd();
  t3->Write();
  f3.Close();
  return 0;
}