class HRichCorrelator: public HReconstructor

cout<<"we are in init() now !"<<endl;

cout<<"in reinit function"<<endl;
 called, when program enters execute() for the first time
      HSpecGeomPar* pSpecGeomPar = (HSpecGeomPar*)getSpecGeomPar();
      cout<<"in reinit: pSpecGeomPar: "<<pSpecGeomPar<<endl;

 various counters
 COUNTER INITIALIZATION

 initialize parameters from RTDB containers
 different cuts to be applied during hit matching

        if (nCounterProcessedNbEvents != 0 &&
  	nCounterProcessedNbEvents%50000==0) HRichUtilFunc::saveHistos(pFileOut,pHistArray);
         if (nCounterProcessedNbEvents != 0 &&
  	   nCounterProcessedNbEvents%10000==0)
         {
  	   cout<<"********** "<<nCounterProcessedNbEvents
  	       <<" EVTS PROCESSED **********"<<endl;
         }

 which kicktracks don't have a match in hitmatchs and vice versa

 remove rich hit obj from HLinearCategory
 afterwards delete empty slot

cout<<"in finalize"<<endl;
      pProfRMTHTtheta=pHistRMTRingHouTraTheta->
  	ProfileX("prof_RMT_HTtheta",0,nNbBinsForHouTra);
      pHistArray->Add(pProfRMTHTtheta);
      pProfRMSHTtheta=pHistRMSRingHouTraTheta->
  	ProfileX("prof_RMS_HTtheta",0,nNbBinsForHouTra);
      pHistArray->Add(pProfRMSHTtheta);
      pProfRMTPMtheta=pHistRMTRingPatMatTheta->
  	ProfileX("prof_RMT_PMtheta",0,nNbBinsForPatMat);
      pHistArray->Add(pProfRMTPMtheta);
      pProfRMHTtheta=pHistRMRingHouTraTheta->
  	ProfileX("prof_RM_HTtheta",0,nNbBinsForHouTra);
      pHistArray->Add(pProfRMHTtheta);
      pProfRMPMtheta=pHistRMRingPatMatTheta->
  	ProfileX("prof_RM_PMtheta",0,nNbBinsForPatMat);
      pHistArray->Add(pProfRMPMtheta);
      pProfRMSPMtheta=pHistRMSRingPatMatTheta->
  	ProfileX("prof_RMS_PMtheta",0,nNbBinsForPatMat);
      pHistArray->Add(pProfRMSPMtheta);
      pProfALLHTtheta=pHistALLRingHouTraTheta->
  	ProfileX("prof_ALL_HTtheta",0,nNbBinsForHouTra);
      pHistArray->Add(pProfALLHTtheta);
      pProfALLPMtheta=pHistALLRingPatMatTheta->
  	ProfileX("prof_ALL_PMtheta",0,nNbBinsForPatMat);
      pHistArray->Add(pProfALLPMtheta);

 HHitMatchHeader is meant to be a bookkeeping obj
 storing info on an event-by-event basis

 check if index i is already stored in iarr
 if not then store it
 iarr is initialized to -2

count arr slots that are different from -2 (ini value)

 iterate over HHitMatch (correlation or tracklet) objs
 in this evt, categorize and count them. Finally fill
 HHitMatchHeader obj for this evt.

create a new HHitMatch object, input sector

 copy info from RICH hit into HHitMatch
 copy info from RICH hit fit into HHitMatch

 copy info from KickTrack into HHitMatch
 retrieve objects from which kick track was built and
 use them to fill HHitMatch

 copy info from MDC Seg into HHitMatch

 copy info from Tof hit to HHitMatch

 copy info from ShowerHitTof to HHitMatch

    cout<<"in fillHitmatch rich - kick"<<endl;

 check whether this correlation object already holds info
 if not fill in detector specific info

 fill histos for all rich hits w/o selection or correlation

 use HHitMatch objs to fill histos under certain conditions
 e.g. only for corr objs that contain a RMT

 shift MDC Seg phi according to the sector

 multiply MdcSeg polar angles by a const factor
 to account for shifted target
 E.g. target was 18mm more downstream than seg finder assumes
 (info W. Koenig, apr01 gen1)
 correction factor can be set in param container, only rough corr !!!
 not necessary if alignment was done for DST production

 apply rich ring quality check
cout<<"in isGoodRing()"<<endl;

 corr RICH - MDC

 is RICH corr with this already MATCHed HIT (MDC part)
 not used !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

 corr RICH TOF

 information by J. Otwinowski
 check whether there was a shower indicating a lepton

 corr RICH SHOWERHITTOF

 corr MDCSEG - SHOWERHITTOF

 corr MDCSEG - TOF

 FIXME: ring counters are unchecked and maybe not up-to-date

gStyle->SetLineWidth(3);