// Silly macro that uses  FilterRootFiles to filter 
// root events

gROOT->Reset();

#include <iostream.h>
#include "FilterRootFiles.C"

void FilterEvents_BeamTest_CR_CERN322()
{

  // Software release to be used

  //   char* softrelease = "EngineeringModel-v6r070329p16";
  char* softrelease = "BeamtestRelease-v3r0907p2";
  
  //char* cuts = {"TkrNumTracks > 0 &&  GltTower > -1 && CalEnergyRaw > 10 && CalCsIRLn > 4"};
  char* cuts = "TMath::Abs((acos(Tkr1ZDir)*(180./3.14159))-180.) < 5 && CalEnergyRaw > 10";




  // merit tuple 



  char* inmerit = {"/nfs/farm/g/glast/u36/Integration/rootData/700000322/v3r0907p2/recon/recon-v1r030603p2_700000322_merit_merit.root"};
  // Cuts to be applied
  



// Root files to read in:

  char* indigi = {"/nfs/farm/g/glast/u36/Integration/rootData/700000322/v3r0907p2/digi/digitization-latte-v1r030603p2_700000322_digi_DIGI.root"};
  char* inrecon = {"/nfs/farm/g/glast/u36/Integration/rootData/700000322/v3r0907p2/recon/recon-v1r030603p2_700000322_recon_RECON.root"};
  
  char* inmc = {""};



  // Root files to write out

  char* outdigi = {"/nfs/farm/g/glast/u33/dpaneque/FilteredRootFiles/digitization-latte-v1r030603p2_700000322_digi_DIGI_5deg.root"};
  char* outrecon = {"/nfs/farm/g/glast/u33/dpaneque/FilteredRootFiles/recon-v1r030603p2_700000322_recon_RECON_5deg.root"};
  char* outmc = {""};



 // Filter events


 cout << "Executing FilterRootFiles " << endl;
 
 FilterRootFiles(softrelease, inmerit, cuts, 
		 indigi, inrecon, inmc, 
		 outdigi, outrecon, outmc);
 


}