// File: hrichdisplaycal.cc
//
// Author: Walter Schoen <W.Schoen@physik.tu-muenchen.de>
// Last update by Walter Schoen: 99/09/14 17:23:38
// Modified by Thomas Eberl: 99/09/02
//
#include "hrichdisplaycal.h"
#include "hevent.h"
#include "hspectrometer.h"
#include "hrichdetector.h"
#include "hcategory.h"
//#include "hcategory.h"
#include "hmatrixcatiter.h"
#include "hlocation.h"
#include "hrichraw.h"
#include "hrichcal.h"
//#include "hrichcalibrater.h"
#include "hdebug.h"
#include "richdef.h"
#include "hrichgeometrypar.h"
#include "hdetgeompar.h"
#include "hruntimedb.h"

 Int_t HRichDisplayCal::producer(void){     
sharedFile = TMapFile::Create("RICHonline.map","RECREATE", 1000000,"shared memory file");

 counter = new TVector(1);
 (*counter)(0)=0;
    for (int i=0;i<sectorNum;i++){
         for (int j=0;j<2;j++){
             richHisto[i][j] = 0;
        }
    }

    char hname1[40];
    char hname2[40];
    

    for (int i =0; i< sectorNum ; i++){

         sprintf(hname1,"singleEvt_%d",i + 1);
         sprintf(hname2,"accumEvt_%d",i + 1);
         
	
         richHisto[i][0] = new TH2F(hname1,"single event",92,1,92,90,1,90);
         richHisto[i][1] = new TH2F(hname2,"accumulated event",92,1,92,90,1,90);
        
    }
  
   sharedFile->Add(counter,"counter");
   sharedFile->Update();
   
    return 0;
}

 HRichDisplayCal::HRichDisplayCal(Text_t *name, Text_t *title) : 
 HReconstructor(name,title) 
{ 

   pRichCal = 0;
}
 HRichDisplayCal::HRichDisplayCal()
{

   pRichCal = 0;
}
 HRichDisplayCal::~HRichDisplayCal(void) {
  
  
   for ( int i =0; i<sectorNum; i++){

  delete richHisto[i][0]; 
  delete richHisto[i][1]; 

   } 
}

 Bool_t HRichDisplayCal::init() {
  
    HRichDetector *pRichDet = 
	(HRichDetector*)gHades->getSetup()->getDetector("Rich");
    HRuntimeDb* rtdb=gHades->getRuntimeDb();
    HRichGeometryPar *pGeomPar = (HRichGeometryPar*)rtdb->
                                            getContainer("RichGeometryParameters");
      
  
    pRichCal=gHades->getCurrentEvent()->getCategory(catRichCal);
    sectorNum= pGeomPar->getSectorsNr(); 
    ncounter=0;
    producer();
    
  
  
    if (!pRichCal) {
      pRichCal=pRichDet->buildCategory(catRichCal);

      if (!pRichCal) return kFALSE;
      else gHades->getCurrentEvent()
                         ->addCategory(catRichCal, pRichCal, "Rich");
    }
    return kTRUE;
} 
 Bool_t HRichDisplayCal::finalize(void) {
    return kTRUE;
}
 Bool_t HRichDisplayCal::display(HLocation &refLoc){
#if DEBUG_LEVEL>2
  gDebuger->enterFunc("HRichDisplay::display");
  gDebuger->message("Raw cat points to %p",fRawCat);
#endif


  HRichCal *pCal=0;

  pCal=(HRichCal *)pRichCal->getObject(refLoc);
  
  if(pCal) {
    int sec = pCal->getSector();
    float charge = pCal->getCharge();
    // if (charge != 0) cout<<" **************** charge from cal:"<<charge<<endl;
    int row = pCal->getRow();
    int col = pCal->getCol();
    //   if (row==10 && col==45) cout<<" hrichdisplaycal  charge : "<<charge<<"  row : "<<row<<"  col  :"<<col<<endl;
      ///for (int k =0; k< sectorNum ; k++){//add logic for sectors here !!!
    richHisto[sec][0]->Fill(col,row,charge);
    richHisto[sec][1]->Fill(col,row,charge);
    //cout <<"hrichdisplay: " <<sec<<" "<<col<<" "<< row <<" "<< charge<<  endl; 
	    // }
      
      // sharedFile->Update();
  }
#if DEBUG_LEVEL>2
  gDebuger->leaveFunc("HRichDisplay::display");
#endif
  return kTRUE;

}
 Int_t HRichDisplayCal::execute() {
    HMatrixCatIter *pIter;
    HRichCal *pCal=NULL;
    HLocation loc;    
    pIter = (HMatrixCatIter*)getCalCat()->MakeIterator(); 
    Int_t n =0;
    for (int k =0; k< sectorNum ; k++){ 
	richHisto[k][0] -> Reset();
    } 
    while((pCal = (HRichCal *)pIter->Next()))
    {
	loc = pIter->getLocation();
	display(loc);
	n++;
    }
    ncounter++;
    (*counter)(0)=ncounter;
    cout << "event: " <<ncounter<< endl;
    sharedFile->Update();
    return kTRUE;
}
ClassImp(HRichDisplayCal)