00001 void testMixer(char *file){
00002
00003 gSystem->Load("MyEvent/MyEvent.so");
00004 gSystem->Load("gamma/analysis/lib/AnaCuts.so");
00005 gSystem->Load("gamma/analysis/lib/EventMixer.so");
00006
00007 TFile *mFile=new TFile(file,"OPEN");
00008 TTree *myEventTree;
00009 myEventTree=(TTree*)mFile->Get("mEventTree");
00010 MyEvent *ev=new MyEvent();
00011 myEventTree->SetBranchAddress("branch",&ev);
00012
00013 EventMixer *mix=new EventMixer("pp05");
00014
00015 Int_t i=0;
00016 while(myEventTree->GetEntry(i) && i<10000){
00017
00018 MyEvent *e_add=new MyEvent(*ev);
00019 mix->addEvent(e_add);
00020
00021 i++;
00022 }
00023
00024 TCanvas *c=new TCanvas();
00025 TH2F *h=new TH2F(*(TH2F*)mix->getMinvMB());
00026 h->Draw("colz");
00027 c->SaveAs("mix.eps");
00028
00029 }
