---

MSLayersKeeperX0AtEta Class Reference

#include <RecoTracker/TkMSParametrization/interface/MSLayersKeeperX0AtEta.h>

Inheritance diagram for MSLayersKeeperX0AtEta:

List of all members.

Public Member Functions
virtual void	init (const edm::EventSetup &iSetup)
virtual const MSLayersAtAngle &	layers (float cotTheta) const
	MSLayersKeeperX0AtEta ()
virtual	~MSLayersKeeperX0AtEta ()
Private Member Functions
float	eta (int idxBin) const
int	idxBin (float eta) const
void	setX0 (std::vector< MSLayer > &, float eta, const SumX0AtEtaDataProvider &) const
Private Attributes
bool	isInitialised
float	theDeltaEta
int	theHalfNBins
std::vector< MSLayersAtAngle >	theLayersData
Friends
class	MSLayersKeeperX0Averaged

---

Detailed Description

Definition at line 9 of file MSLayersKeeperX0AtEta.h.

---

Constructor & Destructor Documentation

MSLayersKeeperX0AtEta::MSLayersKeeperX0AtEta

(

)

[inline]

Definition at line 11 of file MSLayersKeeperX0AtEta.h.

: isInitialised(false) { }

virtual MSLayersKeeperX0AtEta::~MSLayersKeeperX0AtEta

(

)

[inline, virtual]

Definition at line 12 of file MSLayersKeeperX0AtEta.h.

{ }

---

Member Function Documentation

float MSLayersKeeperX0AtEta::eta

(

int

idxBin

)

const [private]

Definition at line 16 of file MSLayersKeeperX0AtEta.cc.

References theDeltaEta, and theHalfNBins.

Referenced by init(), and layers().

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

int MSLayersKeeperX0AtEta::idxBin

(

float

eta

)

const [private]

Definition at line 20 of file MSLayersKeeperX0AtEta.cc.

References e, int, max, theDeltaEta, and theHalfNBins.

Referenced by layers().

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

void MSLayersKeeperX0AtEta::init

(

const edm::EventSetup &

iSetup

)

[virtual]

Reimplemented from MSLayersKeeper.

Definition at line 30 of file MSLayersKeeperX0AtEta.cc.

References eta(), gen_hydjet_hard_muon_trigger_cfg::etaMax, MSLayersAtAngle::findLayer(), MSLayersKeeper::getDataX0(), isInitialised, it, layers(), layout, parsecf::pyparsing::line(), max, MultipleScatteringX0Data::maxEta(), MultipleScatteringX0Data::nBinsEta(), MSLayersKeeper::setDataX0(), setX0(), python::multivaluedict::sort(), theDeltaEta, theHalfNBins, theLayersData, MSLayer::DataX0::x0, and z.

Referenced by MSLayersKeeperX0Averaged::init().

{
  if (isInitialised) return;
  isInitialised = true;
  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(iSetup);
  for (int idxbin = 0; idxbin < 2*theHalfNBins; idxbin++) {
    float etaValue = eta(idxbin);
    float cotTheta = sinh(etaValue);

    vector<MSLayer> layers = layout.detLayers(etaValue,0,iSetup);
    vector<MSLayer> tmplay = layout.otherLayers(etaValue,iSetup);
    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,iSetup);
      vector<MSLayer>::iterator it;
      for (it = candidates.begin(); it != candidates.end(); it++) {
        if (layersAtAngle.findLayer(*it)) continue;
        const MSLayer * found = 0;
        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();
  }
*/
}

const MSLayersAtAngle & MSLayersKeeperX0AtEta::layers

(

float

cotTheta

)

const [virtual]

Implements MSLayersKeeper.

Definition at line 9 of file MSLayersKeeperX0AtEta.cc.

References eta(), idxBin(), and theLayersData.

Referenced by init(), and MSLayersKeeperX0Averaged::init().

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

void MSLayersKeeperX0AtEta::setX0	(	std::vector< MSLayer > &	,
		float	eta,
		const SumX0AtEtaDataProvider &
	)			const [private]

Referenced by init().

---

Friends And Related Function Documentation

friend class MSLayersKeeperX0Averaged [friend]

Definition at line 25 of file MSLayersKeeperX0AtEta.h.

---

Member Data Documentation

bool MSLayersKeeperX0AtEta::isInitialised [private]

Definition at line 22 of file MSLayersKeeperX0AtEta.h.

Referenced by init().

float MSLayersKeeperX0AtEta::theDeltaEta [private]

Definition at line 23 of file MSLayersKeeperX0AtEta.h.

Referenced by eta(), idxBin(), and init().

int MSLayersKeeperX0AtEta::theHalfNBins [private]

Definition at line 23 of file MSLayersKeeperX0AtEta.h.

Referenced by eta(), idxBin(), init(), and MSLayersKeeperX0Averaged::init().

std::vector<MSLayersAtAngle> MSLayersKeeperX0AtEta::theLayersData [private]

Definition at line 24 of file MSLayersKeeperX0AtEta.h.

Referenced by init(), MSLayersKeeperX0Averaged::init(), and layers().

---

The documentation for this class was generated from the following files:

• RecoTracker/TkMSParametrization/interface/MSLayersKeeperX0AtEta.h
• RecoTracker/TkMSParametrization/src/MSLayersKeeperX0AtEta.cc

---