#include <FastHFFibre.h>

Public Member Functions
double	attLength (double lambda)

	FastHFFibre (std::string &name, const DDCompactView &cpv, double cFibre)

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

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

	~FastHFFibre ()

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 19 of file FastHFFibre.h.

Constructor & Destructor Documentation

FastHFFibre::FastHFFibre	(	std::string &	name,
		const DDCompactView &	cpv,
		double	cFibre
	)

Definition at line 18 of file FastHFFibre.cc.

References DDFilteredView::addFilter(), attL, cFibre, DDSpecificsFilter::equals, edm::hlt::Exception, DDFilteredView::firstChild(), getDDDArray(), gpar, i, lambLim, longFL, DDFilteredView::mergedSpecifics(), mergeVDriftHistosByStation::name, nBinAtt, nBinR, radius, DDSpecificsFilter::setCriteria(), shortFL, AlCaHLTBitMon_QueryRunRegistry::string, and relativeConstraints::value.

                                                                                    { 
 
   cFibre = cLight;
   edm::LogInfo("FastCalorimetry") << "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,DDSpecificsFilter::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("FastCalorimetry") << "HFFibre: " << nBinAtt << " attL ";
     for (int it=0; it<nBinAtt; it++)  
       edm::LogInfo("FastCalorimetry") << "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("FastCalorimetry") << "HFFibre: Limits on lambda " << lambLim[0]
                         << " and " << lambLim[1];
 
     // Fibre Lengths
     nb       = 0;
     longFL   = getDDDArray("LongFL",sv,nb);
     edm::LogInfo("FastCalorimetry") << "HFFibre: " << nb << " Long Fibre Length";
     for (int it=0; it<nb; it++) 
       edm::LogInfo("FastCalorimetry") << "HFFibre: longFL[" << it << "] = " 
                           << longFL[it]/cm << " cm";
     nb       = 0;
     shortFL   = getDDDArray("ShortFL",sv,nb);
     edm::LogInfo("FastCalorimetry") << "HFFibre: " << nb << " Short Fibre Length";
     for (int it=0; it<nb; it++) 
       edm::LogInfo("FastCalorimetry") << "HFFibre: shortFL[" << it << "] = " 
                           << shortFL[it]/cm << " cm";
   } else {
     edm::LogError("FastCalorimetry") << "HFFibre: cannot get filtered "
                          << " view for " << attribute << " matching "
                          << name;
     throw cms::Exception("Unknown", "HFFibre")
       << "cannot match " << attribute << " to " << name <<"\n";
   }
 
   // Now geometry parameters
   attribute = "ReadOutName";
   value     = name;
   DDSpecificsFilter filter2;
   DDValue           ddv2(attribute,value,0);
   filter2.setCriteria(ddv2,DDSpecificsFilter::equals);
   DDFilteredView fv2(cpv);
   fv2.addFilter(filter2);
   dodet     = fv2.firstChild();
   if (dodet) {
     DDsvalues_type sv(fv2.mergedSpecifics());
 
     // Special Geometry parameters
     int nb    = -1;
     gpar      = getDDDArray("gparHF",sv,nb);
     edm::LogInfo("FastCalorimetry") << "HFFibre: " << nb <<" gpar (cm)";
     for (int i=0; i<nb; i++) 
       edm::LogInfo("FastCalorimetry") << "HFFibre: gpar[" << i << "] = "
                           << gpar[i]/cm << " cm";
     
     nBinR     = -1;
     radius    = getDDDArray("rTable",sv,nBinR);
     edm::LogInfo("FastCalorimetry") << "HFFibre: " << nBinR <<" rTable (cm)";
     for (int i=0; i<nBinR; i++) 
       edm::LogInfo("FastCalorimetry") << "HFFibre: radius[" << i << "] = "
                           << radius[i]/cm << " cm";
   } else {
     edm::LogError("FastCalorimetry") << "HFFibre: cannot get filtered "
                          << " view for " << attribute << " matching "
                          << name;
     throw cms::Exception("Unknown", "HFFibre")
       << "cannot match " << attribute << " to " << name <<"\n";
   }
 }

FastHFFibre::~FastHFFibre ( )

Definition at line 108 of file FastHFFibre.cc.

108 {}

Member Function Documentation

double FastHFFibre::attLength ( double lambda )

Definition at line 110 of file FastHFFibre.cc.

References attL, i, j, lambLim, and nBinAtt.

Referenced by FastHFShowerLibrary::getHits().

                                            {
 
   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("FastCalorimetry") << "HFFibre::attLength for Lambda " << lambda
                           << " index " << i  << " " << j << " Att. Length " 
                       << att;
 #endif
   return att;
 }

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

protected

Definition at line 177 of file FastHFFibre.cc.

References DDfetch(), DDValue::doubles(), edm::hlt::Exception, LogDebug, and relativeConstraints::value.

Referenced by FastHFFibre().

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

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

Definition at line 128 of file FastHFFibre.cc.

References cFibre, cond::rpcobgas::time, and zShift().

Referenced by FastHFShowerLibrary::getHits().

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

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

Definition at line 140 of file FastHFFibre.cc.

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

Referenced by FastHFShowerLibrary::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("FastCalorimetry") << "HFFibre::zShift for point " << point
                << " (R = " << hR/cm << " cm, Index = " << ieta 
                << ", depth = " << depth << ", Fibre Length = " 
                << length/cm       << " cm = " << zFibre/cm  
                << " cm)";
 #endif
   return zFibre;
 }