CaloGeometryHelper.h

Go to the documentation of this file.

#ifndef FastSimulation_CaloGeometryTools_CaloGeometryHelper
#define FastSimulation_CaloGeometryTools_CaloGeometryHelper

#include "FastSimulation/CalorimeterProperties/interface/Calorimeter.h"
#include "Geometry/CaloTopology/interface/CaloDirection.h"
#include "FastSimulation/CaloGeometryTools/interface/BaseCrystal.h"

#include <vector>

class DetId;
class Crystal;

namespace edm { 
  class ParameterSet;
}

class CaloGeometryHelper:public Calorimeter
{

 public:

  typedef math::XYZVector XYZVector;
  typedef math::XYZVector XYZPoint;

  CaloGeometryHelper();
  CaloGeometryHelper(const edm::ParameterSet& fastCalo);
  ~CaloGeometryHelper();


  // more user friendly getClosestCell  
  DetId getClosestCell(const XYZPoint& point, bool ecal, bool central) const;

  // more user friendly getWindow
  void getWindow(const DetId& pivot,int s1,int s2,std::vector<DetId> &) const;
  
  
  double preshowerZPosition(int layer) const
  {
    return (layer==1) ? psLayer1Z_: psLayer2Z_ ; 
  }
  
  // the Crystal constructor
  void buildCrystal(const DetId& id,Crystal&) const;

  void initialize(double bField);  

  // get the <=8 neighbours
  const std::vector<DetId> & getNeighbours(const DetId& det) const ;

  inline double magneticField() const {return bfield_;}
  
  // temporary. Requires a missing geometry tool 
  bool borderCrossing(const DetId&, const DetId&) const ;

  bool move(DetId& cell, const CaloDirection& dir,bool fast=true) const;

  inline bool preshowerPresent() const {return preshowerPresent_;};

 private:
  void buildNeighbourArray();
  void buildCrystalArray();
  bool simplemove(DetId& cell, const CaloDirection& dir) const;
  bool diagonalmove(DetId& cell, const CaloDirection& dir) const;

 private:
    // Preshower layer positions
  double psLayer1Z_,psLayer2Z_;

  // array of neighbours the hashed index is used for the first vector
  std::vector<std::vector<DetId> > barrelNeighbours_;
  std::vector<std::vector<DetId> > endcapNeighbours_;

  std::vector<BaseCrystal> barrelCrystals_;
  std::vector<BaseCrystal> endcapCrystals_;

  bool neighbourmapcalculated_;
  
  //mag field at 0,0,0
  double bfield_;
  bool preshowerPresent_;
};
#endif