INT histo_booking(void){
   ... ...
  
  /*defination of histograms , UTILITY_ID_BASE is 100000*/
  
  /* chamber resolutions */
  HBOOK1( UTILITY_ID_BASE+20, "Resolution ch1 X", 80, -10, 10, 0.f);
  HBOOK1( UTILITY_ID_BASE+21, "Resolution ch1 Y", 80, -10, 10, 0.f);
  HBOOK1( UTILITY_ID_BASE+22, "Resolution ch1 Z", 80, -10, 10, 0.f);
  HBOOK2( UTILITY_ID_BASE+26, "ch1, delta Y-phi", 18,0.,360.,80,-10.,10., 0.f);
  HBOOK2( UTILITY_ID_BASE+27, "ch1, delta X-phi", 18,0.,360.,80,-10.,10., 0.f);
  HBOOK2( UTILITY_ID_BASE+28, "ch1, deltaphi-phi", 18,0.,360.,80,-5.,5., 0.f);

  HBOOK1( UTILITY_ID_BASE+30, "Resolution ch2 X", 80, -10, 10, 0.f);
  HBOOK1( UTILITY_ID_BASE+31, "Resolution ch2 Y", 80, -10, 10, 0.f);
  HBOOK1( UTILITY_ID_BASE+32, "Resolution ch2 Z", 80, -10, 10, 0.f);
  HBOOK2( UTILITY_ID_BASE+36, "ch2, delta Y-phi", 18,0.,360.,80,-10.,10., 0.f);
  HBOOK2( UTILITY_ID_BASE+37, "ch2, delta X-phi", 18,0.,360.,80,-10.,10., 0.f);
  HBOOK2( UTILITY_ID_BASE+38, "ch2, deltaphi-phi", 18,0.,360.,80,-5.,5., 0.f);
  ... ...
}

INT histo_filling(EVENT_HEADER *pheader, void *pevent) {

  ... ...

    /*  resolution of chambers */
  if ((trig[0] & TRIG_COSMIC_WED) &&
      mwpc->wc1_n_hits == 2       &&
      mwpc->wc2_n_hits == 2       &&   
            mwpc->wc1_theta[0]>-1.      &&
      mwpc->wc2_theta[0]>-1.      &&
      asum->csi_hit == 2          &&
      asum->csi_sum > 150.
      )
    {       

      if ((dif(wc1_c_p1, wc1_p1)+dif(wc1_c_p2, wc1_p2))<(dif(wc1_c_p1, wc1_p2)+dif(wc1_c_p2, wc1_p1)))
	{
	  HF1(UTILITY_ID_BASE+20, wc1_c_p1[0]-wc1_p1[0], 1.f);
	  HF1(UTILITY_ID_BASE+20, wc1_c_p2[0]-wc1_p2[0], 1.f);
	  HF1(UTILITY_ID_BASE+21, wc1_c_p1[1]-wc1_p1[1], 1.f);
	  HF1(UTILITY_ID_BASE+21, wc1_c_p2[1]-wc1_p2[1], 1.f);
	  HF1(UTILITY_ID_BASE+22, wc1_c_p1[2]-wc1_p1[2], 1.f);
	  HF1(UTILITY_ID_BASE+22, wc1_c_p2[2]-wc1_p2[2], 1.f);
	  HF2(UTILITY_ID_BASE+26, mwpc->wc1_phi[0],wc1_c_p1[1]-wc1_p1[1], 1.f);
	  HF2(UTILITY_ID_BASE+26, mwpc->wc1_phi[1],wc1_c_p2[1]-wc1_p2[1], 1.f);
	  HF2(UTILITY_ID_BASE+27, mwpc->wc1_phi[0],wc1_c_p1[0]-wc1_p1[0], 1.f);
	  HF2(UTILITY_ID_BASE+27, mwpc->wc1_phi[1],wc1_c_p2[0]-wc1_p2[0], 1.f);
	  wc1_c_phi1 = atan2(wc1_c_p1[1], wc1_c_p1[0])*180/PI;
	  if(wc1_c_phi1<0.0) wc1_c_phi1 += 360.;
	  phi_diff = wc1_c_phi1 - mwpc->wc1_phi[0];
	  if(phi_diff>180.) phi_diff -= 360.;
	  HF2(UTILITY_ID_BASE+28, mwpc->wc1_phi[0], phi_diff, 1.f);
	  wc1_c_phi2 = atan2(wc1_c_p2[1], wc1_c_p2[0])*180/PI;
	  if(wc1_c_phi2<0.0) wc1_c_phi2 += 360.;
	  phi_diff = wc1_c_phi2 - mwpc->wc1_phi[1];
	  if(phi_diff>180.) phi_diff -= 360.;
	  HF2(UTILITY_ID_BASE+28, mwpc->wc1_phi[1], phi_diff, 1.f);
	}
      else 
	{
	  HF1(UTILITY_ID_BASE+20, wc1_c_p1[0]-wc1_p2[0], 1.f);
	  HF1(UTILITY_ID_BASE+20, wc1_c_p2[0]-wc1_p1[0], 1.f);
	  HF1(UTILITY_ID_BASE+21, wc1_c_p1[1]-wc1_p2[1], 1.f);
	  HF1(UTILITY_ID_BASE+21, wc1_c_p2[1]-wc1_p1[1], 1.f);
	  HF1(UTILITY_ID_BASE+22, wc1_c_p1[2]-wc1_p2[2], 1.f);
	  HF1(UTILITY_ID_BASE+22, wc1_c_p2[2]-wc1_p1[2], 1.f);
	  HF2(UTILITY_ID_BASE+26, mwpc->wc1_phi[1],wc1_c_p1[1]-wc1_p2[1], 1.f);
	  HF2(UTILITY_ID_BASE+26, mwpc->wc1_phi[0],wc1_c_p2[1]-wc1_p1[1], 1.f);
	  HF2(UTILITY_ID_BASE+27, mwpc->wc1_phi[1],wc1_c_p1[0]-wc1_p2[0], 1.f);
	  HF2(UTILITY_ID_BASE+27, mwpc->wc1_phi[0],wc1_c_p2[0]-wc1_p1[0], 1.f);
	  wc1_c_phi1 = atan2(wc1_c_p1[1], wc1_c_p1[0])*180/PI;
	  if(wc1_c_phi1<0.0) wc1_c_phi1 += 360.;
	  phi_diff = wc1_c_phi1 - mwpc->wc1_phi[1];
	  if(phi_diff>180.) phi_diff -= 360.;
	  HF2(UTILITY_ID_BASE+28, mwpc->wc1_phi[1], phi_diff, 1.f);
	  wc1_c_phi2 = atan2(wc1_c_p2[1], wc1_c_p2[0])*180/PI;
	  if(wc1_c_phi2<0.0) wc1_c_phi2 += 360.;
	  phi_diff = wc1_c_phi2 - mwpc->wc1_phi[0];
	  if(phi_diff>180.) phi_diff -= 360.;
	  HF2(UTILITY_ID_BASE+28, mwpc->wc1_phi[0], phi_diff, 1.f);
	}
      if ((dif(wc2_c_p1, wc2_p1)+dif(wc2_c_p2, wc2_p2))<(dif(wc2_c_p1, wc2_p2)+dif(wc2_c_p2, wc2_p1)))
	{
	  HF1(UTILITY_ID_BASE+30, wc2_c_p1[0]-wc2_p1[0], 1.f);
	  HF1(UTILITY_ID_BASE+30, wc2_c_p2[0]-wc2_p2[0], 1.f);
	  HF1(UTILITY_ID_BASE+31, wc2_c_p1[1]-wc2_p1[1], 1.f);
	  HF1(UTILITY_ID_BASE+31, wc2_c_p2[1]-wc2_p2[1], 1.f);
       	  HF1(UTILITY_ID_BASE+32, wc2_c_p1[2]-wc2_p1[2], 1.f);
	  HF1(UTILITY_ID_BASE+32, wc2_c_p2[2]-wc2_p2[2], 1.f);
	  HF2(UTILITY_ID_BASE+36, mwpc->wc2_phi[0],wc2_c_p1[1]-wc2_p1[1], 1.f);
	  HF2(UTILITY_ID_BASE+36, mwpc->wc2_phi[1],wc2_c_p2[1]-wc2_p2[1], 1.f);
	  HF2(UTILITY_ID_BASE+37, mwpc->wc2_phi[0],wc2_c_p1[0]-wc2_p1[0], 1.f);
	  HF2(UTILITY_ID_BASE+37, mwpc->wc2_phi[1],wc2_c_p2[0]-wc2_p2[0], 1.f);
	  wc2_c_phi1 = atan2(wc2_c_p1[1], wc2_c_p1[0])*180/PI;
	  if(wc2_c_phi1<0.0) wc2_c_phi1 += 360.;
	  phi_diff = wc2_c_phi1 - mwpc->wc2_phi[0];
	  if(phi_diff>180.) phi_diff -= 360.;
	  HF2(UTILITY_ID_BASE+38, mwpc->wc2_phi[0], phi_diff, 1.f);
	  wc2_c_phi2 = atan2(wc2_c_p2[1], wc2_c_p2[0])*180/PI;
	  if(wc2_c_phi2<0.0) wc2_c_phi2 += 360.;
	  phi_diff = wc2_c_phi2 - mwpc->wc2_phi[1];
	  if(phi_diff>180.) phi_diff -= 360.;
	  HF2(UTILITY_ID_BASE+38, mwpc->wc2_phi[1], phi_diff, 1.f);
	}
      else 
	{
	  HF1(UTILITY_ID_BASE+30, wc2_c_p1[0]-wc2_p2[0], 1.f);
	  HF1(UTILITY_ID_BASE+30, wc2_c_p2[0]-wc2_p1[0], 1.f);
	  HF1(UTILITY_ID_BASE+31, wc2_c_p1[1]-wc2_p2[1], 1.f);
	  HF1(UTILITY_ID_BASE+31, wc2_c_p2[1]-wc2_p1[1], 1.f);
	  HF1(UTILITY_ID_BASE+32, wc2_c_p1[2]-wc2_p2[2], 1.f);
	  HF1(UTILITY_ID_BASE+32, wc2_c_p2[2]-wc2_p1[2], 1.f);
	  HF2(UTILITY_ID_BASE+36, mwpc->wc2_phi[1],wc2_c_p1[1]-wc2_p2[1], 1.f);
	  HF2(UTILITY_ID_BASE+36, mwpc->wc2_phi[0],wc2_c_p2[1]-wc2_p1[1], 1.f);
	  HF2(UTILITY_ID_BASE+37, mwpc->wc2_phi[1],wc2_c_p1[0]-wc2_p2[0], 1.f);
	  HF2(UTILITY_ID_BASE+37, mwpc->wc2_phi[0],wc2_c_p2[0]-wc2_p1[0], 1.f);
	  wc2_c_phi1 = atan2(wc2_c_p1[1], wc2_c_p1[0])*180/PI;
	  if(wc2_c_phi1<0.0) wc2_c_phi1 += 360.;
	  phi_diff = wc2_c_phi1 - mwpc->wc2_phi[1];
	  if(phi_diff>180.) phi_diff -= 360.;
	  HF2(UTILITY_ID_BASE+38, mwpc->wc2_phi[1], phi_diff, 1.f);
	  wc2_c_phi2 = atan2(wc2_c_p2[1], wc2_c_p2[0])*180/PI;
	  if(wc2_c_phi2<0.0) wc2_c_phi2 += 360.;
	  phi_diff = wc2_c_phi2 - mwpc->wc2_phi[0];
	  if(phi_diff>180.) phi_diff -= 360.;
	  HF2(UTILITY_ID_BASE+38, mwpc->wc2_phi[0], phi_diff, 1.f);

	}
     }
 ... ...
}