HFFibre Class Reference

#include <HFFibre.h>

Public Member Functions
double	attLength (double lambda)
	HFFibre (std::string &name, const DDCompactView &cpv, edm::ParameterSet const &p)
void	initRun (HcalDDDSimConstants *)
double	tShift (const G4ThreeVector &point, int depth, int fromEndAbs=0)
double	zShift (const G4ThreeVector &point, int depth, int fromEndAbs=0)
	~HFFibre ()

Protected Member Functions
std::vector< double >	getDDDArray (const std::string &, const DDsvalues_type &, int &)

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

Detailed Description

Definition at line 20 of file HFFibre.h.

Constructor & Destructor Documentation

HFFibre::HFFibre	(	std::string &	name,
		const DDCompactView &	cpv,
		edm::ParameterSet const &	p
	)

Definition at line 19 of file HFFibre.cc.

References DDFilteredView::addFilter(), attL, cFibre, equals, Exception, DDFilteredView::firstChild(), getDDDArray(), edm::ParameterSet::getParameter(), lambLim, longFL, DDFilteredView::mergedSpecifics(), dataset::name, nBinAtt, DDSpecificsFilter::setCriteria(), shortFL, and AlCaHLTBitMon_QueryRunRegistry::string.

                                    {

  edm::ParameterSet m_HF = p.getParameter<edm::ParameterSet>("HFShower");
  cFibre           = c_light*(m_HF.getParameter<double>("CFibre"));
  
  edm::LogInfo("HFShower") << "HFFibre:: Speed of light in fibre " << cFibre
               << " m/ns";

  std::string attribute = "Volume"; 
  std::string value     = "HF";
  DDSpecificsFilter filter1;
  DDValue           ddv1(attribute,value,0);
  filter1.setCriteria(ddv1,DDCompOp::equals);
  DDFilteredView fv1(cpv);
  fv1.addFilter(filter1);
  bool dodet = fv1.firstChild();

  if (dodet) {
    DDsvalues_type sv(fv1.mergedSpecifics());

    // Attenuation length
    nBinAtt      = -1;
    attL         = getDDDArray("attl",sv,nBinAtt);
    edm::LogInfo("HFShower") << "HFFibre: " << nBinAtt << " attL ";
    for (int it=0; it<nBinAtt; it++) 
      edm::LogInfo("HFShower") << "HFFibre: attL[" << it << "] = " 
                   << attL[it]*cm << "(1/cm)";

    // Limits on Lambda
    int nb   = 2;
    std::vector<double>  nvec = getDDDArray("lambLim",sv,nb);
    lambLim[0] = static_cast<int>(nvec[0]);
    lambLim[1] = static_cast<int>(nvec[1]);
    edm::LogInfo("HFShower") << "HFFibre: Limits on lambda " << lambLim[0]
                 << " and " << lambLim[1];

    // Fibre Lengths
    nb       = 0;
    longFL   = getDDDArray("LongFL",sv,nb);
    edm::LogInfo("HFShower") << "HFFibre: " << nb << " Long Fibre Length";
    for (int it=0; it<nb; it++) 
      edm::LogInfo("HFShower") << "HFFibre: longFL[" << it << "] = " 
                   << longFL[it]/cm << " cm";
    nb       = 0;
    shortFL   = getDDDArray("ShortFL",sv,nb);
    edm::LogInfo("HFShower") << "HFFibre: " << nb << " Short Fibre Length";
    for (int it=0; it<nb; it++) 
      edm::LogInfo("HFShower") << "HFFibre: shortFL[" << it << "] = " 
                   << shortFL[it]/cm << " cm";
  } else {
    edm::LogError("HFShower") << "HFFibre: cannot get filtered "
                  << " view for " << attribute << " matching "
                  << name;
    throw cms::Exception("Unknown", "HFFibre")
      << "cannot match " << attribute << " to " << name <<"\n";
  }
}

HFFibre::~HFFibre

(

)

Definition at line 78 of file HFFibre.cc.

{}

Member Function Documentation

double HFFibre::attLength

(

double

lambda

)

Definition at line 97 of file HFFibre.cc.

References attL, i, createfilelist::int, j, lambLim, and nBinAtt.

Referenced by HFShowerLibrary::fillHits(), HFShower::getHits(), and HFShowerParam::HFShowerParam().

                                       {

  int i = int(nBinAtt*(lambda - lambLim[0])/(lambLim[1]-lambLim[0]));

  int j =i;
  if (i >= nBinAtt) 
    j = nBinAtt-1;
  else if (i < 0)
    j = 0;
  double att = attL[j];
#ifdef DebugLog
  edm::LogInfo("HFShower") << "HFFibre::attLength for Lambda " << lambda
               << " index " << i  << " " << j << " Att. Length " 
               << att;
#endif
  return att;
}

std::vector< double > HFFibre::getDDDArray ( const std::string &  str, const DDsvalues_type &  sv, int &  nmin  )

protected

Definition at line 164 of file HFFibre.cc.

References DDfetch(), DDValue::doubles(), Exception, LogDebug, harvestTrackValidationPlots::str, and relativeConstraints::value.

Referenced by HFFibre().

                                 {

#ifdef DebugLog
  LogDebug("HFShower") << "HFFibre:getDDDArray called for " << str 
               << " with nMin " << nmin;
#endif
  DDValue value(str);
  if (DDfetch(&sv,value)) {
#ifdef DebugLog
    LogDebug("HFShower") << value;
#endif
    const std::vector<double> & fvec = value.doubles();
    int nval = fvec.size();
    if (nmin > 0) {
      if (nval < nmin) {
    edm::LogError("HFShower") << "HFFibre : # of " << str << " bins " 
                  << nval << " < " << nmin << " ==> illegal";
    throw cms::Exception("Unknown", "HFFibre")
      << "nval < nmin for array " << str <<"\n";
      }
    } else {
      if (nval < 1 && nmin != 0) {
    edm::LogError("HFShower") << "HFFibre : # of " << str << " bins " 
                  << nval << " < 1 ==> illegal (nmin=" 
                  << nmin << ")";
    throw cms::Exception("Unknown", "HFFibre")
      << "nval < 1 for array " << str <<"\n";
      }
    }
    nmin = nval;
    return fvec;
  } else {
    if (nmin != 0) {
      edm::LogError("HFShower") << "HFFibre : cannot get array " << str;
      throw cms::Exception("Unknown", "HFFibre")
    << "cannot get array " << str <<"\n";
    } else {
      std::vector<double> fvec;
      return fvec;
    }
  }
}

void HFFibre::initRun

(

HcalDDDSimConstants *

hcons

)

Definition at line 80 of file HFFibre.cc.

References HcalDDDSimConstants::getGparHF(), HcalDDDSimConstants::getRTableHF(), gpar, i, createfilelist::int, nBinR, and radius.

Referenced by HFShower::initRun(), HFShowerParam::initRun(), and HFShowerLibrary::initRun().

                                                {

  // Now geometry parameters
  gpar      = hcons->getGparHF();
  edm::LogInfo("HFShower") << "HFFibre: " << gpar.size() <<" gpar (cm)";
  for (unsigned int i=0; i<gpar.size(); i++)
    edm::LogInfo("HFShower") << "HFFibre: gpar[" << i << "] = "
                             << gpar[i]/cm << " cm";

  radius    = hcons->getRTableHF();
  nBinR     = (int)(radius.size());
  edm::LogInfo("HFShower") << "HFFibre: " << radius.size() <<" rTable (cm)";
  for (unsigned int i=0; i<radius.size(); i++)
    edm::LogInfo("HFShower") << "HFFibre: radius[" << i << "] = "
                             << radius[i]/cm << " cm";
}

double HFFibre::tShift	(	const G4ThreeVector &	point,
		int	depth,
		int	fromEndAbs = 0
	)

Definition at line 115 of file HFFibre.cc.

References cFibre, ntuplemaker::time, and zShift().

Referenced by HFShowerLibrary::fillHits(), HFShower::getHits(), and HFShowerParam::getHits().

                                                                            {

  double zFibre = zShift(point, depth, fromEndAbs);
  double time   = zFibre/cFibre;
#ifdef DebugLog
  edm::LogInfo("HFShower") << "HFFibre::tShift for point " << point
               << " ( depth = " << depth <<", traversed length = " 
               << zFibre/cm  << " cm) = " << time/ns << " ns";
#endif
  return time;
}

double HFFibre::zShift	(	const G4ThreeVector &	point,
		int	depth,
		int	fromEndAbs = 0
	)

Definition at line 127 of file HFFibre.cc.

References gpar, i, longFL, nBinR, radius, shortFL, mathSSE::sqrt(), and geometryCSVtoXML::zz.

Referenced by HFShowerLibrary::fillHits(), HFShowerParam::getHits(), and tShift().

                                                                            { // point is z-local

  double zFibre = 0;
  double hR     = sqrt((point.x())*(point.x())+(point.y())*(point.y()));
  int    ieta   = 0;
  double length = 250*cm;
  if (fromEndAbs < 0) {
    zFibre = 0.5*gpar[1] - point.z(); // Never, as fromEndAbs=0 (?)
  } else {
    // Defines the Radius bin by radial subdivision
    for (int i = nBinR-1; i > 0; --i) if (hR < radius[i]) ieta = nBinR - i - 1;
    // define the length of the fibre
    if (depth == 2) {
      if ((int)(shortFL.size()) > ieta) length = shortFL[ieta];
    } else {
      if ((int)(longFL.size())  > ieta) length = longFL[ieta];
    }
    zFibre = length;
    if (fromEndAbs > 0) {
      zFibre   -= gpar[1]; // Never, as fromEndAbs=0 (M.K. ?)
    } else  {
      double zz = 0.5*gpar[1] + point.z();
      zFibre   -= zz;
    }
    if (depth == 2) zFibre += gpar[0]; // here zFibre is reduced for Short
  }

#ifdef DebugLog
  edm::LogInfo("HFShower") << "HFFibre::zShift for point " << point
               << " (R = " << hR/cm << " cm, Index = " << ieta 
               << ", depth = " << depth << ", Fibre Length = " 
               << length/cm       << " cm = " << zFibre/cm  
               << " cm)";
#endif
  return zFibre;
}

Member Data Documentation

std::vector<double> HFFibre::attL

private

Definition at line 46 of file HFFibre.h.

Referenced by attLength(), and HFFibre().

double HFFibre::cFibre

private

Definition at line 43 of file HFFibre.h.

Referenced by HFFibre(), and tShift().

std::vector<double> HFFibre::gpar

private

Definition at line 44 of file HFFibre.h.

Referenced by initRun(), and zShift().

double HFFibre::lambLim[2]

private

Definition at line 48 of file HFFibre.h.

Referenced by attLength(), and HFFibre().

std::vector<double> HFFibre::longFL

private

Definition at line 45 of file HFFibre.h.

Referenced by HFFibre(), and zShift().

int HFFibre::nBinAtt

private

Definition at line 47 of file HFFibre.h.

Referenced by attLength(), and HFFibre().

int HFFibre::nBinR

private

Definition at line 47 of file HFFibre.h.

Referenced by initRun(), and zShift().

std::vector<double> HFFibre::radius

private

Definition at line 44 of file HFFibre.h.

Referenced by initRun(), and zShift().

std::vector<double> HFFibre::shortFL

private

Definition at line 45 of file HFFibre.h.

Referenced by HFFibre(), and zShift().