HFFibre.h

Go to the documentation of this file.

#ifndef SimG4CMS_HFFibre_h
#define SimG4CMS_HFFibre_h 1
// File: HFFibre.h
// Description: Calculates attenuation length

#include "DetectorDescription/Core/interface/DDsvalues.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "Geometry/HcalCommonData/interface/HcalDDDSimConstants.h"

#include "G4ThreeVector.hh"

#include <vector>
#include <string>

class DDCompactView;    

class HFFibre {
  
public:
  
  //Constructor and Destructor
  HFFibre(const std::string & name, const DDCompactView & cpv,
          edm::ParameterSet const & p);
  ~HFFibre() = default;

  void                initRun(const HcalDDDSimConstants*);
  double              attLength(double lambda);
  double              tShift(const G4ThreeVector& point, int depth, 
                             int fromEndAbs=0);
  double              zShift(const G4ThreeVector& point, int depth, 
                             int fromEndAbs=0);

protected:

  std::vector<double> getDDDArray(const std::string&, 
                                  const DDsvalues_type&, int&);

private:

  double                      cFibre;
  std::vector<double>         gpar, radius;
  std::vector<double>         shortFL, longFL;
  std::vector<double>         attL;
  int                         nBinR, nBinAtt;
  double                      lambLim[2];

};
#endif