MSLayersKeeperX0DetLayer Class Reference

#include <MSLayersKeeperX0DetLayer.h>

Inheritance diagram for MSLayersKeeperX0DetLayer:

Public Member Functions
virtual void	init (const edm::EventSetup &iSetup)
virtual MSLayer	layer (const DetLayer *layer) const
virtual const MSLayersAtAngle &	layers (float cotTheta) const
	MSLayersKeeperX0DetLayer ()
virtual	~MSLayersKeeperX0DetLayer ()
Public Member Functions inherited from MSLayersKeeper
virtual	~MSLayersKeeper ()

Private Attributes
bool	isInitialised
MSLayersAtAngle	theLayersData

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

Constructor & Destructor Documentation

MSLayersKeeperX0DetLayer::MSLayersKeeperX0DetLayer ( )

inline

Definition at line 9 of file MSLayersKeeperX0DetLayer.h.

: isInitialised(false) { }

virtual MSLayersKeeperX0DetLayer::~MSLayersKeeperX0DetLayer ( )

inlinevirtual

Definition at line 10 of file MSLayersKeeperX0DetLayer.h.

{ }

Member Function Documentation

void MSLayersKeeperX0DetLayer::init ( const edm::EventSetup &  iSetup )

virtual

Reimplemented from MSLayersKeeper.

Definition at line 13 of file MSLayersKeeperX0DetLayer.cc.

References GeomDetEnumerators::barrel, PixelRecoLineRZ::cotLine(), citk::if(), and slope.

{
  if (isInitialised) return;
  isInitialised = true;
  //  vector<MSLayer> allLayers = MSLayersKeeperX0DetLayerGeom().detLayers();
  //MP
  vector<MSLayer> allLayers;
  theLayersData = MSLayersAtAngle(allLayers);
  
  vector<MSLayer>::iterator it;
  PixelRecoPointRZ zero(0.,0.); 
  for (it = allLayers.begin(); it != allLayers.end(); it++) {
    PixelRecoPointRZ middle = it->face()== GeomDetEnumerators::barrel ?
        PixelRecoPointRZ(it->position(), it->range().mean())
      : PixelRecoPointRZ(it->range().mean(), it->position());

    float cotTheta = PixelRecoLineRZ(zero,middle).cotLine();
    float x0 =  getDataX0(*it).x0;

    DataX0 dataX0;
    if (it->face()== GeomDetEnumerators::barrel) {
      float sumX0D = theLayersData.sumX0D(zero, middle);
      dataX0 = DataX0(x0, sumX0D, cotTheta);      
    } else {
      float hrange= (it->range().max()-it->range().min())/2.;
      float cot1 = it->position()/(it->range().mean()-hrange/2);
      float cot2 = it->position()/(it->range().mean()+hrange/2);
      PixelRecoLineRZ line1(zero,cot1);
      PixelRecoLineRZ line2(zero,cot2);
      float sum1 = theLayersData.sumX0D(zero,it->crossing(line1).first);
      float sum2 = theLayersData.sumX0D(zero,it->crossing(line2).first);
      float slope = (sum2-sum1)/(1/cot2-1/cot1);
      float sumX0D = sum1 + slope*(1/cotTheta-1/cot1);
      dataX0 = DataX0(x0, sumX0D, cotTheta);
      dataX0.setForwardSumX0DSlope(slope);
    }
    setDataX0(*it, dataX0);
    theLayersData.update(*it);
  }
  //cout << "MSLayersKeeperX0DetLayer LAYERS: "<<endl;
  //theLayersData.print();
}

virtual MSLayer MSLayersKeeperX0DetLayer::layer ( const DetLayer *  layer ) const

inlinevirtual

Reimplemented from MSLayersKeeper.

Definition at line 12 of file MSLayersKeeperX0DetLayer.h.

Referenced by geometryXMLparser.DTAlignable::index(), and geometryXMLparser.CSCAlignable::index().

    {return *theLayersData.findLayer(MSLayer(layer)); }

virtual const MSLayersAtAngle& MSLayersKeeperX0DetLayer::layers ( float  cotTheta ) const

inlinevirtual

Implements MSLayersKeeper.

Definition at line 14 of file MSLayersKeeperX0DetLayer.h.

    {return theLayersData;}

Member Data Documentation

bool MSLayersKeeperX0DetLayer::isInitialised

private

Definition at line 18 of file MSLayersKeeperX0DetLayer.h.

MSLayersAtAngle MSLayersKeeperX0DetLayer::theLayersData

private

Definition at line 19 of file MSLayersKeeperX0DetLayer.h.