MSLayersKeeperX0AtEta Class Reference

#include <MSLayersKeeperX0AtEta.h>

Inheritance diagram for MSLayersKeeperX0AtEta:

Public Member Functions
const MSLayersAtAngle &	layers (float cotTheta) const override
	MSLayersKeeperX0AtEta (const GeometricSearchTracker &tracker, const MagneticField &bfield)
	~MSLayersKeeperX0AtEta () override
Public Member Functions inherited from MSLayersKeeper
virtual MSLayer	layer (const DetLayer *dl) const
virtual	~MSLayersKeeper ()

Private Member Functions
float	eta (int idxBin) const
int	idxBin (float eta) const

Static Private Member Functions
static void	setX0 (std::vector< MSLayer > &, float eta, const SumX0AtEtaDataProvider &)

Private Attributes
float	theDeltaEta
int	theHalfNBins
std::vector< MSLayersAtAngle >	theLayersData

Friends
class	MSLayersKeeperX0Averaged

Additional Inherited Members
Protected Types inherited from MSLayersKeeper
typedef MSLayer::DataX0	DataX0
Static Protected Member Functions inherited from MSLayersKeeper
static const DataX0 &	getDataX0 (const MSLayer &l)
static void	setDataX0 (MSLayer &l, const DataX0 &x0Data)

Detailed Description

Definition at line 13 of file MSLayersKeeperX0AtEta.h.

Constructor & Destructor Documentation

MSLayersKeeperX0AtEta::MSLayersKeeperX0AtEta	(	const GeometricSearchTracker &	tracker,
		const MagneticField &	bfield
	)

Definition at line 37 of file MSLayersKeeperX0AtEta.cc.

References newFWLiteAna::bin, HLT_FULL_cff::candidates, PVValHelper::eta, Puppi_cff::etaMax, MSLayersAtAngle::findLayer(), newFWLiteAna::found, LayerTriplets::layers(), makeLayoutFileForGui::layout, geometryCSVtoXML::line, SiStripPI::max, MultipleScatteringX0Data::maxEta(), MultipleScatteringX0Data::nBinsEta(), and SiPixelPI::zero.

                                                                                                               {
  const float BIG = 99999.;

  // set size from data file
  MultipleScatteringX0Data msX0data;
  theHalfNBins = msX0data.nBinsEta();
  float etaMax = msX0data.maxEta();
  theDeltaEta = (theHalfNBins != 0) ? etaMax / theHalfNBins : BIG;

  theLayersData = vector<MSLayersAtAngle>(max(2 * theHalfNBins, 1));
  MultipleScatteringGeometry layout(tracker);
  for (int idxbin = 0; idxbin < 2 * theHalfNBins; idxbin++) {
    float etaValue = eta(idxbin);
    float cotTheta = sinh(etaValue);

    vector<MSLayer> layers = layout.detLayers(etaValue, 0, bfield);
    vector<MSLayer> tmplay = layout.otherLayers(etaValue);
    layers.insert(layers.end(), tmplay.begin(), tmplay.end());
    sort(layers.begin(), layers.end());
    setX0(layers, etaValue, msX0data);
    theLayersData[idxbin] = MSLayersAtAngle(layers);
    PixelRecoPointRZ zero(0., 0.);
    PixelRecoLineRZ line(zero, cotTheta);
    vector<MSLayer>::iterator it;
    for (it = layers.begin(); it != layers.end(); it++) {
      float x0 = getDataX0(*it).x0;
      float sumX0D = theLayersData[idxbin].sumX0D(zero, it->crossing(line).first);
      setDataX0(*it, DataX0(x0, sumX0D, cotTheta));
      theLayersData[idxbin].update(*it);
    }
  }

  // add layers not seen from nominal vertex but crossed if
  // vertex seperated from nominal by less than 3 sigma
  for (int idxbin = 0; idxbin < 2 * theHalfNBins; idxbin++) {
    float etaValue = eta(idxbin);
    for (int isign = -1; isign <= 1; isign += 2) {
      float z = isign * 15.9;  //3 sigma from zero
      const MSLayersAtAngle& layersAtAngle = theLayersData[idxbin];
      vector<MSLayer> candidates = layout.detLayers(etaValue, z, bfield);
      vector<MSLayer>::iterator it;
      for (it = candidates.begin(); it != candidates.end(); it++) {
        if (layersAtAngle.findLayer(*it))
          continue;
        const MSLayer* found = nullptr;
        int bin = idxbin;
        while (!found) {
          bin--;
          if (bin < 0)
            break;
          found = theLayersData[bin].findLayer(*it);
        }
        bin = idxbin;
        while (!found) {
          bin++;
          if (bin > 2 * theHalfNBins - 1)
            break;
          found = theLayersData[bin].findLayer(*it);
        }
        if (found)
          theLayersData[idxbin].update(*found);
      }
    }
  }

  // cout << "LAYERS, size=: "<<theLayersData.size()<< endl;
  /*
  for (int idxbin = 0; idxbin <= theHalfNBins; idxbin+=25) {
    float etaValue = eta(idxbin);
    const MSLayersAtAngle & layers= theLayersData[idxbin];
    cout << "ETA: "<< etaValue <<" (bin:"<<idxbin<<") #layers:"
         <<layers.size()<<endl;
    layers.print();
  }
  for (int idxbin = 2*theHalfNBins-1; idxbin > theHalfNBins; idxbin-=25) {
    float etaValue = eta(idxbin);
    const MSLayersAtAngle & layers= theLayersData[idxbin];
    cout << "ETA: "<< etaValue <<" (bin:"<<idxbin<<") #layers:"
         <<layers.size()<<endl;
    layers.print();
  }
*/
}

MSLayersKeeperX0AtEta::~MSLayersKeeperX0AtEta ( )

overridedefault

Member Function Documentation

float MSLayersKeeperX0AtEta::eta ( int  idxBin ) const

private

Definition at line 25 of file MSLayersKeeperX0AtEta.cc.

Referenced by Particle.Particle::__str__(), Jet.Jet::jetID(), and Jet.Jet::puJetId().

{ return (idxBin + 0.5 - theHalfNBins) * theDeltaEta; }

int MSLayersKeeperX0AtEta::idxBin ( float  eta ) const

private

Definition at line 28 of file MSLayersKeeperX0AtEta.cc.

References funct::abs(), alignCSCRings::e, and SiStripPI::max.

                                                 {
  float ieta = eta / theDeltaEta;
  if (std::abs(ieta) >= theHalfNBins - 1.e-3)
    return (eta > 0) ? max(2 * theHalfNBins - 1, 0) : 0;
  else
    return int(ieta + theHalfNBins);
}

const MSLayersAtAngle & MSLayersKeeperX0AtEta::layers ( float  cotTheta ) const

overridevirtual

Implements MSLayersKeeper.

Definition at line 19 of file MSLayersKeeperX0AtEta.cc.

References PVValHelper::eta.

Referenced by MSLayersKeeperX0Averaged::MSLayersKeeperX0Averaged().

                                                                         {
  float eta = unsafe_asinhf(cotTheta);
  return theLayersData[idxBin(eta)];
}

void MSLayersKeeperX0AtEta::setX0 ( std::vector< MSLayer > &  layers, float  eta, const SumX0AtEtaDataProvider &  sumX0  )

staticprivate

Definition at line 122 of file MSLayersKeeperX0AtEta.cc.

References GeomDetEnumerators::barrel, geometryCSVtoXML::line, sqr(), mathSSE::sqrt(), and SumX0AtEtaDataProvider::sumX0atEta().

                                                                                                         {
  const float BIG = 99999.;
  float cotTheta = sinh(eta);
  float sinTheta = 1 / sqrt(1 + sqr(cotTheta));
  float cosTheta = cotTheta * sinTheta;

  float sumX0atAngleLast = 0.;
  vector<MSLayer>::iterator il;
  for (il = layers.begin(); il != layers.end(); il++) {
    PixelRecoLineRZ line(PixelRecoPointRZ(0., 0.), cotTheta);
    float rL = (*il).crossing(line).first.r();
    float rN = (il + 1 != layers.end()) ? (il + 1)->crossing(line).first.r() : BIG;
    float rBound = (rL + rN) / 2.;
    float sumX0atAngle = sumX0.sumX0atEta(eta, rBound);

    float dX0 = (il->face() == GeomDetEnumerators::barrel) ? (sumX0atAngle - sumX0atAngleLast) * sinTheta
                                                           : (sumX0atAngle - sumX0atAngleLast) * fabs(cosTheta);

    setDataX0(*il, DataX0(dX0, 0, cotTheta));
    sumX0atAngleLast = sumX0atAngle;
  }
}

Friends And Related Function Documentation

friend class MSLayersKeeperX0Averaged

friend

Definition at line 28 of file MSLayersKeeperX0AtEta.h.

Member Data Documentation

float MSLayersKeeperX0AtEta::theDeltaEta

private

Definition at line 26 of file MSLayersKeeperX0AtEta.h.

int MSLayersKeeperX0AtEta::theHalfNBins

private

Definition at line 25 of file MSLayersKeeperX0AtEta.h.

Referenced by MSLayersKeeperX0Averaged::MSLayersKeeperX0Averaged().

std::vector<MSLayersAtAngle> MSLayersKeeperX0AtEta::theLayersData

private

Definition at line 27 of file MSLayersKeeperX0AtEta.h.

Referenced by MSLayersKeeperX0Averaged::MSLayersKeeperX0Averaged().