void iso()
{
   TFile *f = new TFile("May26_2.root");

   TTree *t2 = (TTree*)f->Get("demo/AnaTree");

   TCanvas *c1 = new TCanvas("c1","c1",800,600);
/*   c1->Divide(2,2);

   TH1F *h2 = new TH1F("h2","hist_nonIsoMu",100,0,300);
   h2->GetXaxis()->SetTitle("selNonIsoMu_pt");
   h2->GetYaxis()->SetTitle("entries");
   c1->cd(1);
   t2->Draw("selNonIsoMu_pt >> h2","","Hist");
*/
   TH1F *h1 = new TH1F("h1","hist_nonIsoMu_cuts",100,0,300);
   h1->GetXaxis()->SetTitle("selNonIsoMu_cuts_pt");
   h1->GetYaxis()->SetTitle("entries");

 float selNonIsoMu_pt, selNonIsoMu_iso;
int selNonIsoMu_idTight, selNonIsoMu_ismuon;
t2->SetBranchAddress("selNonIsoMu_pt",&selNonIsoMu_pt);
t2->SetBranchAddress("selNonIsoMu_idTight",&selNonIsoMu_idTight);
t2->SetBranchAddress("selNonIsoMu_ismuon",&selNonIsoMu_ismuon);
t2->SetBranchAddress("selNonIsoMu_iso",&selNonIsoMu_iso);

if(selNonIsoMu_pt < 20 ) continue;
if(selNonIsoMu_idTight == 0) continue;
if(selNonIsoMu_ismuon == 0) continue;
if(selNonIsoMu_iso < 0.12) continue;
//c1->cd(2);
   t2->Draw("selNonIsoMu_pt >> h1","","Hist");

}