gROOT->Reset();
#include <iostream.h>

void clusterSubtract(){

TFile *f = new TFile("jethist_pt20_Cent10_low.root");
gStyle->SetNdivisions(505,"Y");//as per Horner's suggestion
gStyle->SetNdivisions(505,"X");
gStyle->SetPadBorderMode(0);
gStyle->SetPadBorderSize(0);
gStyle->SetPadColor(kWhite);
gStyle->SetPadTickX(1);
gStyle->SetPadTickY(1);
gROOT->SetStyle("Plain");

TCanvas* c1 = new TCanvas("DeltaPhiMix","DeltaPhiMix",0,0,650,500);

Double_t scale_data;
for(Int_t i = 0; i < 51; i++)
scale_data += dRSeedAssocClustersData->GetBinContent(i+50);

Double_t scale_mix;
for(Int_t j = 0; j < 51; j++)
scale_mix += dRSeedAssocClusters_Mix->GetBinContent(j+50);

dRSeedAssocClusters_Mix->Scale(scale_data/scale_mix);
dRSeedAssocClusters_Mix->Draw();

dRSeedAssocClustersData->SetLineColor(2);
dRSeedAssocClustersData->Draw("same");

Double_t data;
for(Int_t k = 0; k < 29; k++)
data += dRSeedAssocClustersData->GetBinContent(k+1);

Double_t mix;
for(Int_t m = 0; m < 29; m++)
mix += dRSeedAssocClusters_Mix->GetBinContent(m+1);

cout << "fraction of background clusters " << mix/data << endl;
cout << "fraction of real clusters " << 1 - mix/data << endl;
}