{
  #include <string>
  gSystem->Load("libDelphes.so");
  gSystem->Load("libMinuit.so");
  gROOT->ProcessLine(".L Ttilde.C++");
  
  TChain *inputChain = new TChain("Delphes"); //sometimes the tree name is "CollectionTree"
  for (int ifile=1; ifile<21; ifile++){
    int imass = ifile; 
    stringstream filemass;
    filemass << imass;    
    string input = "../../MG5/forRequest33_2/output/Ts_ht/Ts_ht_unweighted_events_33TEV_L3_M1000_"+filemass.str()+"_NoPileUp.root";
    inputChain->Add(input.c_str());    
  }
  
  cout << "Chain has " << inputChain->GetEntries() << " entries" << endl;

  Ttilde(inputChain);
}