// -------------------------------------------------------------------------
// -----                 CbmStsHitProducerIdel header file             -----
// -----                  Created 10/01/06  by V. Friese               -----
// -------------------------------------------------------------------------


/** CbmStsHitProducerIdeal.h
 *@author V.Friese <v.friese@gsi.de>
 **
 ** The ideal hit producer produces hits of type CbmStsMapsHit by copying
 ** the MCPoint position. The position error is set to 1 mum, much 
 ** smaller than can be obtained by any detector. Using the hits from 
 ** this HitProducer is thus equivalent to using MC information
 ** directly, but with the correct data interface.
 **/


#ifndef MVDHYBRIDHITPRODUCER_H
#define MVDHYBRIDHITPRODUCER_H 

#include "CbmTask.h"
#include "MvdDigiPar.h"
#include "MvdGeoMappingPar.h"
#include "MvdHit.h"
#include "MvdPoint.h"
#include "MvdPixel.h"
#include "MvdDigiPixel.h"
#include "CbmGeoVector.h"
#include "CbmGeoTransform.h"
#include "TVector3.h"
#include "TRandom.h"
#include "TGeoMatrix.h"
#include "TGeoBBox.h"

#include <string>
#include <vector>
 
class TClonesArray;

class MvdHybridHitProducer : public CbmTask
{
 public:

  /** Default constructor **/  
 MvdHybridHitProducer();
  
  MvdHybridHitProducer(Double_t lx, Double_t ly, Double_t treshold, Double_t noise);

  /** Destructor **/
  virtual ~MvdHybridHitProducer();


  /** Virtual method Init **/
  virtual void SetParContainers();
  virtual InitStatus Init();
  virtual InitStatus ReInit();

  /** Virtual method Exec **/
  virtual void Exec(Option_t* opt);


 private:
  
  TString fBranchName;
  /** Input array of MvdPoints **/
  TClonesArray* fPointArray;

  /** Output array of MvdHits **/
  TClonesArray* fHitArray;  
  TClonesArray* fPixelArray;
  TClonesArray* fFePixelArray;
  
  MvdDigiPar* fDigiPar;
  MvdGeoMappingPar* fGeoMappingPar;
  
  void Register();
  void Reset();  
  void ProduceHits();
  
  TGeoHMatrix GetTransformation (std::string detName);
  void GetLocalHitPoints(MvdPoint* myPoint, CbmGeoVector& myHitIn, CbmGeoVector& myHitOut);
  MvdHit CalcGlobalPoint(std::vector<MvdPixel> pixels);
  TVector3 GetSensorDimensions(std::string detName);  
  Double_t flx;	//pixel widh in x;
  Double_t fly;	//pixel width in y;
  Double_t fthreshold; //pixel threshold in electrons
  Double_t fnoise; //pixel noise in electrons
  Int_t pixelHits;
  ClassDef(MvdHybridHitProducer,5);

};

#endif