PndEmcHitProducer.h

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                 PndEmcHitProducer header file                 -----
// -----                 Created 14/08/06  by S.Spataro                -----
// -------------------------------------------------------------------------
#ifndef PNDEMCHITPRODUCER_H
#define PNDEMCHITPRODUCER_H

#include <PndPersistencyTask.h>
#include <map>
#include <string>
#include <limits>

#include "TH1F.h"
#include "TFile.h"
#include "TVector3.h"

#include "PndEmcHit.h"
#include "PndEmcStructure.h"
#include "PndEmcMapper.h"


class TClonesArray;
class TObjectArray;

class PndEmcDigiPar;
class PndEmcGeoPar;
class PndEmcDigiNonuniformityPar;

using std::cout;
using std::endl;
using std::map;

class PndEmcHitProducer : public PndPersistencyTask
{

 public:
  typedef std::map<Int_t, Float_t> mapper;

  PndEmcHitProducer();

  PndEmcHitProducer(Bool_t val);

  ~PndEmcHitProducer();


  virtual InitStatus Init();


  virtual void Exec(Option_t* opt);

  PndEmcHit* AddHit(Int_t trackID, Int_t detID, Float_t energy, Float_t time,std::vector <Int_t> &mctruth, FairMultiLinkedData entering, FairMultiLinkedData exiting);

  //PndEmcHit* AddHit(Int_t trackID, Int_t detID, Float_t energy, Float_t time,std::vector<PndEmcPoint*> pointList);
  // not implemented
  //void CreateStructure();

        void SetStorageOfData(Bool_t val); // Method to specify whether hits are stored or not.
        void SetNonuniformityFile(const char * filename){fNonuniformityFile = filename;};
  void SetDayOne(bool d=true){fDayOne=d;};

        void FinishTask();
 private:
        Int_t fUse_nonuniformity;
        TString fNonuniformityFile;
        TClonesArray* fPointArray;
        TClonesArray* fMCTrackArray;

        TClonesArray* fHitArray;
        //TClonesArray* fMcTrackArray;

        TObjArray *fVolumeArray;

        Int_t fMapVersion;
        Float_t fEnergyThreshold;

        mapper emcX;
        mapper emcY;
        mapper emcZ;

        PndEmcStructure* fEmcStr;
        PndEmcMapper* fMapper;

        PndEmcDigiPar*    fDigiPar;      
        PndEmcGeoPar*     fGeoPar;       
        PndEmcDigiNonuniformityPar *fNonuniformityPar;

        virtual void SetParContainers();

        Bool_t fStoreHits; // Flag which specify whether hits are stored or not

        void cleansortmclist( std::vector <Int_t> &newlist,TClonesArray* mcTrackArray);

        PndEmcHitProducer(const  PndEmcHitProducer& L);
        PndEmcHitProducer& operator= (const  PndEmcHitProducer&) {return *this;}

        ClassDef(PndEmcHitProducer,1);

        map<Int_t, Float_t> fTrackEnergy;
        map<Int_t, Float_t> fTrackTime;  //time of first point
        map<Int_t, std::vector <Int_t> > fTrackMcTruth;  //McTruth
        map<Int_t, std::vector <Int_t> > fPointMatch; //DetId , PointIds with same DetId
        map<Int_t, FairMultiLinkedData> fTrackEntering; //DetId, link to tracks entering same DetId
        map<Int_t, FairMultiLinkedData> fTrackExiting; //DetId, link to track exiting same DetId

  bool fDayOne;
  bool AcceptDayOne(PndEmcPoint* p);

};
#endif