EcalTimeMapDigitizer.h

Go to the documentation of this file.

#ifndef EcalSimAlgos_EcalTimeMapDigitizer_h
#define EcalSimAlgos_EcalTimeMapDigitizer_h

#include "CalibFormats/CaloObjects/interface/CaloTSamples.h"
#include "SimDataFormats/CrossingFrame/interface/MixCollection.h"
#include "SimDataFormats/CaloHit/interface/PCaloHit.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
#include "CalibFormats/CaloObjects/interface/CaloTSamplesBase.h"
#include "DataFormats/EcalDigi/interface/EcalDigiCollections.h"

class CaloSubdetectorGeometry ;

class EcalTimeMapDigitizer
{
 public:

  struct time_average{

    static const unsigned short time_average_capacity=10;

    DetId id;
    float average_time[time_average_capacity];
    unsigned int nhits[time_average_capacity];
    float tot_energy[time_average_capacity];
    
    time_average(const DetId& myId):
      id(myId)
    {
      for (unsigned int i(0);i<time_average_capacity;++i)
    {
      average_time[i]=0;
      nhits[i]=0;
      tot_energy[i]=0;
    }

    };

    void 
    calculateAverage()
    {
      for (unsigned int i(0);i<time_average_capacity;++i)
    {
      if (nhits[i]>0)
        average_time[i]=average_time[i]/tot_energy[i];
      else
        average_time[i]=0;
    }
    };
    
    void 
    setZero()
    {
      for (unsigned int i(0);i<time_average_capacity;++i)
    {
      average_time[i]=0;
      nhits[i]=0;
      tot_energy[i]=0;
    }
    };

    bool
    zero()
    {
      for (unsigned int i(0);i<time_average_capacity;++i)
    {
      if(nhits[i]>0)
        return false;
    }
      return true;
    };
    
  };

  typedef time_average TimeSamples;

  typedef EcalTimeDigi Digi;

  typedef std::vector< unsigned int > VecInd ;

  explicit EcalTimeMapDigitizer(EcalSubdetector myDet);
  
  virtual ~EcalTimeMapDigitizer();
  
  void add(const std::vector<PCaloHit> & hits, int bunchCrossing);
  
  void setGeometry( const CaloSubdetectorGeometry* geometry ) ;

  void initializeMap();
  
  void run(EcalTimeDigiCollection& output  );

  void finalizeHits();
  
  inline void setTimeLayerId( const int& layerId ) { m_timeLayerId = layerId; };
  
  int getTimeLayerId() { return m_timeLayerId; };
  
  EcalSubdetector subdetector() { return m_subDet; };

  int  minBunch() const ;

  int  maxBunch() const  ;

  void blankOutUsedSamples() ; // blank out previously used elements
  

/*   const CaloVHitFilter* hitFilter() const ; */

 private:
  
  TimeSamples* findSignal( const DetId& detId );

  VecInd& index() ;
  
  const VecInd& index() const ;

  unsigned int samplesSize() const;

  unsigned int samplesSizeAll() const;

  const TimeSamples*  operator[]( unsigned int i ) const;

  TimeSamples*  operator[]( unsigned int i );

  TimeSamples*  vSam( unsigned int i );

  TimeSamples*  vSamAll( unsigned int i );

  const TimeSamples*  vSamAll( unsigned int i ) const;

  EcalSubdetector m_subDet;
  //time difference between bunches

  static const int BUNCHSPACE=25;

  static const float MIN_ENERGY_THRESHOLD; //50 KeV threshold to consider a valid hit in the timing detector

  static const int    m_minBunch=-4;
  static const int    m_maxBunch=5;
  
  int m_timeLayerId;
  
  const CaloSubdetectorGeometry* m_geometry ;
  
  double timeOfFlight( const DetId& detId , int layer) const; 

  std::vector<TimeSamples> m_vSam;

  VecInd m_index;  

};

#endif