{
// #include<iostream.h>

cout << " Executing Electroweak Analysis Macro ... " << endl;

// Clean up :
gROOT->Reset();

// Open input file :
TFile zeeFile("zee_anal.root");

// List the contents :
// zeeFile->ls();

// Go to the directory with the same name as the analysis module :
zeeFile->cd("HighPtElectrons");

// List the contents :
// zeeFile->ls();

// Find the various histograms and assign pointers :
TH1F *emClusterPt = (TH1F*) gROOT.FindObject("EM cluster pT   (GeV)");
TH1F *matchingTrackPt = (TH1F*) gROOT.FindObject("Matching Track pT   (GeV)");
TH1F *eOverP = (TH1F*) gROOT.FindObject("E over p for cluster, matching track");
TH1F *iMass = (TH1F*) gROOT.FindObject("Invariant Mass (2 EM Clusters, GeV)");

// Make a canvas to plot the histograms on :
c1 = new TCanvas("c1","Electroweak Analysis Variables",0,0,700,500);
c1->SetFillColor(18);

// Display fit results :
gStyle->SetOptFit(1);

// Make the pads for the individual histograms.
// The coordinates are (x,y) of the bottom-left-hand and upper-right-hand
// corners, respectively :
pad1 = new TPad("pad1"," ",0.01,0.01,0.99,0.99,21);
pad1->Draw();

// Activate pads and plot histogram :
pad1->cd();
pad1->SetFillColor(19);
iMass->Fit("gaus");
}