EcalElectronicsMapper.h

Go to the documentation of this file.

#ifndef _EcalElectronicsMapper_HH_
#define _EcalElectronicsMapper_HH_

/*
 *\ Class EcalElectronicsMapper
 *
 * Mapper for the ECAL electronics 
 
 * \file EcalElectronicsMapper.h
 *
 * \author N. Almeida
 * \author G. Franzoni
 *
*/

#include <iostream>
#include <fstream>
#include <string>
#include <map>

#include "DCCRawDataDefinitions.h"
#include <DataFormats/EcalDetId/interface/EcalTrigTowerDetId.h>
#include <DataFormats/EcalDetId/interface/EcalElectronicsId.h>
#include <DataFormats/EcalDetId/interface/EcalScDetId.h>
#include <DataFormats/EcalDigi/interface/EcalDataFrame.h>
#include <DataFormats/EcalDigi/interface/EcalTriggerPrimitiveDigi.h>
#include <DataFormats/EcalDigi/interface/EcalPseudoStripInputDigi.h>
#include <DataFormats/EcalDigi/interface/EcalSrFlag.h>
#include <iostream>
#include <sstream>

#include <DataFormats/EcalDetId/interface/EcalTriggerElectronicsId.h>

class EcalElectronicsMapping;

class EcalElectronicsMapper {
public:
  EcalElectronicsMapper(unsigned int numbOfXtalTSamples, unsigned int numbOfTriggerTSamples);

  ~EcalElectronicsMapper();

  void setEcalElectronicsMapping(const EcalElectronicsMapping*);

  void deletePointers();
  void resetPointers();

  bool setActiveDCC(unsigned int dccId);

  bool setDCCMapFilePath(std::string);

  std::string getDCCMapFilePath() const { return pathToMapFile_; }

  //  bool readDCCMapFile();
  //  bool readDCCMapFile(std::string );

  bool makeMapFromVectors(std::vector<int>&, std::vector<int>&);

  const std::map<unsigned int, unsigned int>& getDCCMap() const { return myDCCMap_; }

  DetId* getDetIdPointer(unsigned int feChannel, unsigned int strip, unsigned int xtal) {
    return xtalDetIds_[smId_ - 1][feChannel - 1][strip - 1][xtal - 1];
  }

  EcalTrigTowerDetId* getTTDetIdPointer(unsigned int tccId, unsigned int tower) {
    return ttDetIds_[tccId - 1][tower - 1];
  }

  EcalElectronicsId* getTTEleIdPointer(unsigned int tccId, unsigned int tower) {
    return ttEleIds_[tccId - 1][tower - 1];
  }

  EcalTriggerPrimitiveDigi* getTPPointer(unsigned int tccId, unsigned int tower) {
    return ttTPIds_[tccId - 1][tower - 1];
  }

  EcalScDetId* getSCDetIdPointer(unsigned int smId, unsigned int feChannel) {
    return scDetIds_[smId - 1][feChannel - 1];
  }

  EcalElectronicsId* getSCElectronicsPointer(unsigned int smId, unsigned int feChannel) {
    return scEleIds_[smId - 1][feChannel - 1];
  }

  EcalPseudoStripInputDigi* getPSInputDigiPointer(unsigned int tccId, unsigned int towerId, unsigned int psId) {
    return psInput_[tccId - 1][towerId - 1][psId - 1];
  }

  EcalPseudoStripInputDigi* getPSInputDigiPointer(unsigned int tccId, unsigned int psCounter) {
    return getPSInputDigiPointer(tccId, tTandPs_[tccId - 1][psCounter - 1][0], tTandPs_[tccId - 1][psCounter - 1][1]);
  }

  // this getter method needs be clarified.
  // Changed by Ph.G. on July 1, 09: return a vector instead of a single
  // element. One SRF can be associated to two  supercrystals, because of
  // channel grouping.
  std::vector<EcalSrFlag*> getSrFlagPointer(unsigned int feChannel) { return srFlags_[smId_ - 1][feChannel - 1]; }

  std::vector<unsigned int>* getTccs(unsigned int smId) { return mapSmIdToTccIds_[smId]; }

  unsigned int getActiveDCC() { return dccId_; }

  unsigned int getActiveSM() { return smId_; }

  unsigned int numbXtalTSamples() { return numbXtalTSamples_; }

  unsigned int numbTriggerTSamples() { return numbTriggerTSamples_; }

  unsigned int getUnfilteredTowerBlockLength() { return unfilteredFEBlockLength_; }

  unsigned int getEBTCCBlockLength() { return ebTccBlockLength_; }

  unsigned int getEETCCBlockLength() { return eeTccBlockLength_; }

  unsigned int getSRPBlockLength() { return srpBlockLength_; }

  unsigned int getDCCId(unsigned int aSMId) const;

  unsigned int getSMId(unsigned int aDCCId) const;

  unsigned int getNumChannelsInDcc(unsigned int aDCCId) { return numChannelsInDcc_[aDCCId - 1]; }

  const EcalElectronicsMapping* mapping() { return mappingBuilder_; }

  bool isTCCExternal(unsigned int TCCId);

  friend std::ostream& operator<<(std::ostream& o, const EcalElectronicsMapper& aEcalElectronicsMapper);

  // Mantain this here as long as everything is moved to a general mapping
  enum SMGeom_t {
    kModules = 4,            // Number of modules per supermodule
    kTriggerTowers = 68,     // Number of trigger towers per supermodule
    kTowersInPhi = 4,        // Number of trigger towers in phi
    kTowersInEta = 17,       // Number of trigger towers in eta
    kCrystals = 1700,        // Number of crystals per supermodule
    kPns = 10,               // Number of PN laser monitoring diodes per supermodule
    kCrystalsInPhi = 20,     // Number of crystals in phi
    kCrystalsInEta = 85,     // Number of crystals in eta
    kCrystalsPerTower = 25,  // Number of crystals per trigger tower
    kCardsPerTower = 5,      // Number of VFE cards per trigger tower
    kChannelsPerCard = 5,    // Number of channels per VFE card
    TTMAPMASK = 100
  };

private:
  void fillMaps();
  unsigned int computeUnfilteredFEBlockLength();
  unsigned int computeEBTCCBlockLength();
  unsigned int computeEETCCBlockLength();

  std::string pathToMapFile_;

  unsigned int numbXtalTSamples_;

  unsigned int numbTriggerTSamples_;

  std::map<unsigned int, unsigned int> myDCCMap_;

  std::map<unsigned int, std::vector<unsigned int>*> mapSmIdToTccIds_;

  unsigned int dccId_;

  unsigned int smId_;

  unsigned int unfilteredFEBlockLength_;

  unsigned int srpBlockLength_;

  unsigned int ebTccBlockLength_, eeTccBlockLength_;

  static const unsigned int numChannelsInDcc_[NUMB_SM];

  // ARRAYS OF DetId
  DetId* xtalDetIds_[NUMB_SM][NUMB_FE][NUMB_STRIP][NUMB_XTAL];
  EcalScDetId* scDetIds_[NUMB_SM][NUMB_FE];
  EcalElectronicsId* scEleIds_[NUMB_SM][NUMB_FE];
  EcalTrigTowerDetId* ttDetIds_[NUMB_TCC][NUMB_FE];
  EcalElectronicsId* ttEleIds_[NUMB_TCC][NUMB_FE];
  EcalTriggerPrimitiveDigi* ttTPIds_[NUMB_TCC][NUMB_FE];
  std::vector<EcalSrFlag*> srFlags_[NUMB_SM][NUMB_FE];
  EcalPseudoStripInputDigi* psInput_[NUMB_TCC][TCC_EB_NUMBTTS][NUMB_STRIP];

  short tTandPs_[NUMB_TCC][5 * EcalTrigTowerDetId::kEBTowersPerSM][2];

  const EcalElectronicsMapping* mappingBuilder_;

  // functions and fields to work with 'ghost' VFEs:
public:
  // check, does the given [FED (dcc), CCU (tower), VFE (strip)] belongs
  // to the list of VFEs with 'ghost' channels
  bool isGhost(const int FED, const int CCU, const int VFE);

private:
  void setupGhostMap();
  std::map<int, std::map<int, std::map<int, bool> > > ghost_;
};

#endif