void drawtwopacks(const int arm = 1,
		  const short SELGAP = 4,
		  const int runno = 247385)
{
  char infile[200];
  //  sprintf(infile,"onlmon_minbias_orig.root");
  sprintf(infile,"onlmon_minbias_mod.root");

  cout << " infile " << infile << endl;
  
  TFile *f = TFile::Open(infile);
  
  gStyle->SetTitleH(0.03);   
  gStyle->SetTitleW(0.5);

  const char *ARMSTR[2] = {"SOUTH", "NORTH"};
  const char *ORIENTSTR[2] = {"HORIZ", "VERT"};
  const char *ostr[2] = {"H", "V"};

  const int N_FILE_TYPES = 2;
  const char *file_type[2] = { "gif", "ps" };
  /*
const int N_FILE_TYPES = 1;
  const char *file_type[2] = { "ps" };
  */

  int colors[5] = {kBlack, kRed, 8, kBlue, 28};//{kBlack,kRed,kBlue,kGreen,kYellow};
  /*
  TGFrame* rootWin = (TGFrame*) gClient->GetRoot();
  int display_sizex = rootWin->GetDefaultWidth();
  int display_sizey = rootWin->GetDefaultHeight();
  */
  int display_sizex = 700;
  int display_sizey = 500;

  TCanvas* C1[2];
  C1[0]=new TCanvas("c1_0","",0,0,display_sizex,display_sizey);
  C1[1]=new TCanvas("c1_2","",0,0,display_sizex,display_sizey);
  C1[0]->SetFillColor(kWhite);
  C1[1]->SetFillColor(kWhite);
  C1[0]->Divide(3,2);
  C1[1]->Divide(3,2);
  TH1* h1[5];
  
  TLegend* legend = 0L;
  //TLegend* legend;
  char histname[128];
  char info[128];
  int muid_panel = 0;
  for(int iorient=0;iorient<2;iorient++){
    
    for(int ipanel=0;ipanel<6;ipanel++){
      muid_panel = ipanel;
      if (ipanel > 2) { // bottom is sort of reversed in MUID numbering
	muid_panel = 4 - ipanel + 4; // reverse around 4
      }
      
      for(int igap=0;igap<5;igap++){
	sprintf(histname,"MUI_twopack%d%d%d%d",arm, igap, muid_panel, iorient);
	h1[igap]=(TH1*)gROOT->FindObject(histname);
      }    
      int maxhits=0;
      int maxindex=0;
      for(int igap=0;igap<5;igap++){
	if(h1[igap]->GetMaximum()>maxhits){
	  maxhits=h1[igap]->GetMaximum();
	  maxindex=igap;
	}
      }
      
      C1[iorient]->cd(ipanel+1);
      //      C1[iorient]->GetPad(ipanel+1)->SetLogy();
      h1[maxindex]->SetStats(false);
      // change title for first one
      sprintf(info,"Arm %d Orient %d Panel %d", arm, iorient, muid_panel);
      h1[maxindex]->SetTitle(info);
      h1[maxindex]->Draw();
      for(int igap=0;igap<5;igap++){
	h1[igap]->SetLineColor(colors[igap]);
	h1[igap]->Draw("same");
      }
      
      if (ipanel == 1){
	legend = new TLegend(0.65,0.65,0.88,0.88,"Hits");
	legend->AddEntry(h1[0], "Gap 0");
	legend->AddEntry(h1[1], "Gap 1");
	legend->AddEntry(h1[2], "Gap 2");
	legend->AddEntry(h1[3], "Gap 3");
	legend->AddEntry(h1[4], "Gap 4");
	legend->SetTextSize(0.036);
	//legend->SetLineWidth(3);
	legend->DrawClone();
	
	// arm and orient info also
	Float_t x = 45 + 5*(1-iorient);
	Float_t yarm = maxhits*0.45 + 0.25;
	Float_t yorient = maxhits*0.35 - 0.25;
	
	armtex = new TLatex(x,yarm,ARMSTR[arm]);
	armtex->SetTextSize(0.05);
	armtex->SetLineWidth(2);
	armtex->Draw();
	
	orienttex = new TLatex(x,yorient,ORIENTSTR[iorient]);
	orienttex->SetTextSize(0.05);
	orienttex->SetLineWidth(2);
	orienttex->Draw();
      }
    }
    C1[iorient]->cd();
    Float_t yextra = 0.95;
    // add time
    TDatime dt;
    TText *tdate = new TText(0.99,yextra,dt.AsString());
    tdate->SetNDC();
    tdate->SetTextSize(0.02);
    tdate->SetTextAlign(31);
    tdate->Draw();
    
    // file info   
    TText *tfile = new TText(0.31,yextra,infile);
    tfile->SetNDC();
    tfile->SetTextSize(0.02);
    tfile->SetTextAlign(31);
    tfile->Draw();

    char text[64];

    C1[iorient]->Update();
    for (int ifile_type = 0; ifile_type<N_FILE_TYPES; ifile_type++) {

      sprintf (text, "plots/%sTwoPack_%d-%d-%d.%s",
	       ostr[iorient],arm, SELGAP, runno,file_type[ifile_type]);
      C1[iorient]->SaveAs(text);
    }

  }
}