void TreeAddWeightJuliet(char *outputfilename="PhysicsWeightedE-1_9string.MergedReduced.root") {
gROOT->Reset();
    
  //load shared libraries...not sure yet which of these are absolutely
  //critical... so load them all for now

  Load();

  char *input1="MCEventsPhysics1-rawfadc.tree.root";  


	SimulationCheckerTree *evt = new SimulationCheckerTree();
        TChain ch("JulietTree");
        ch.SetBranchAddress("branch",&evt);

	ch.Add(input1);
	ch.Add(input2);

        ch.Merge(outputfilename);

}

//////////////////////////////////////////////////
void Load()
{
gSystem->Load("libicetray");
gSystem->Load("libroot-icetray");
gSystem->Load("libdataclasses");
gSystem->Load("libphys-services");
gSystem->Load("libjuliet-interface.so");
//gSystem->Load("libsaved-weighting.so");
gSystem->Load("libweighting-module");
gSystem->Load("libfirst-guess");
}

void Style()
{
  gStyle->SetOptFit(0000);
  gStyle->SetPalette(1);
  gStyle->SetLabelSize(0.05,"x");
  gStyle->SetLabelSize(0.05,"y");
  gStyle->SetLabelSize(0.05,"z");
  gStyle->SetNdivisions(1005,"x");
  gStyle->SetNdivisions(1005,"y");
  gStyle->SetNdivisions(1005,"z");
  gStyle->SetOptTitle(1);
  gStyle->SetOptStat(0);
  gStyle->SetTitleW(0.98);
  gStyle->SetTitleH(0.08);
  gStyle->SetTitleBorderSize(0.0);

//margin
  Float_t rightmargin = 0.175;
  Float_t leftmargin = 0.175;
  Float_t topmargin = 0.15;
  Float_t bottommargin = 0.15;
  

    gPad->SetRightMargin(rightmargin);
    gPad->SetLeftMargin(leftmargin);
    gPad->SetTopMargin(topmargin);
    gPad->SetBottomMargin(bottommargin);
  
 }