{
	gStyle->SetOptStat(0);
	c_inj = new TCanvas("Mass Distribution","Mass Distribution",600,400);
     c_inj->SetBorderMode(0);
     c_inj->SetBorderSize(2);
     c_inj->SetGridx();
     c_inj->SetGridy();
     c_inj->SetRightMargin(0.0562249);
     c_inj->SetBottomMargin(0.1165803);
	
	TH2F *inj_events = new TH2F("Mass Distribution","Mass Distribution",100,5,25,100,5,25);
	inj_events->GetXaxis()->SetTitle("Component mass 1(M_{#odot})");
	inj_events->GetYaxis()->SetTitle("Component mass 2(M_{#odot})");
    inj_events->GetXaxis()->SetTitleOffset(1.3);
    inj_events->GetYaxis()->SetTitleOffset(1.3);
    inj_events->GetXaxis()->CenterTitle(kTRUE);
    inj_events->GetYaxis()->CenterTitle(kTRUE);
    inj_events->GetXaxis()->SetTitleFont(42);
    inj_events->GetXaxis()->SetLabelFont(42);
    inj_events->GetYaxis()->SetTitleFont(42);
    inj_events->GetYaxis()->SetLabelFont(42);
    inj_events->SetLineWidth(4);
    inj_events->SetLineColor(2);
    inj_events->SetFillStyle(3017);    
    inj_events->SetFillColor(kRed);
	
	TRandom rnd;
	double rnum, m1, m2;
	for(int i=0; i<80000; i++) {
		rnum = rnd.Rndm();
		m1=5+20*rnum;
		rnum = rnd.Rndm();
		m2=5+20*rnum;
		if(m1>m2) inj_events->Fill(m1,m2);
	}
	inj_events->Draw("cf");
	
}