CastorTopology.h

Go to the documentation of this file.

#ifndef GEOMETRY_CALOTOPOLOGY_CASTORTOPOLOGY_H
#define GEOMETRY_CALOTOPOLOGY_CASTORTOPOLOGY_H 1

#include <vector>
#include "DataFormats/HcalDetId/interface/HcalCastorDetId.h"
#include "Geometry/CaloTopology/interface/CaloSubdetectorTopology.h"

class CastorTopology : public CaloSubdetectorTopology {
public:

  CastorTopology();
  void exclude(const HcalCastorDetId& id);
  void exclude(int zside); 
  void exclude(int zside, HcalCastorDetId::Section section);
  int exclude(int zside, HcalCastorDetId::Section section1, int isec1, int imod1, HcalCastorDetId::Section section2, int isec2, int imod2);

  virtual bool valid(const HcalCastorDetId& id) const;

  virtual bool validRaw(const HcalCastorDetId& id) const;

  virtual std::vector<DetId> incSector(const DetId& id) const;
  
  virtual std::vector<DetId> incModule(const DetId& id) const;

  //** I have to put this here since they inherit from CaloSubdetectorTopology
  virtual std::vector<DetId> east(const DetId& id) const;
  virtual std::vector<DetId> west(const DetId& id) const; 
  virtual std::vector<DetId> north(const DetId& id) const;
  virtual std::vector<DetId> south(const DetId& id) const;
  virtual std::vector<DetId> up(const DetId& id) const;
  virtual std::vector<DetId> down(const DetId& id) const;
  
  // how many channels (deph) for a given section
  int ncells(HcalCastorDetId::Section section) const;

  //return first and last cell of each section
  int firstCell(HcalCastorDetId::Section section)const;  
  int lastCell(HcalCastorDetId::Section section)const;
 
 private:
  
  std::vector<HcalCastorDetId> exclusionList_;
  
  bool excludeEM_, excludeHAD_, excludeZP_, excludeZN_;
  
  int firstEMModule_, lastEMModule_, firstHADModule_, lastHADModule_;
   
  bool isExcluded(const HcalCastorDetId& id) const;
  
  int firstEMModule() const {return firstEMModule_;}
  int firstHADModule() const {return firstHADModule_;}  
  int lastEMModule()  const {return lastEMModule_;}
  int lastHADModule() const {return lastHADModule_;}  

};


#endif