#include <MSLayersKeeperX0DetLayer.h>

Inheritance diagram for MSLayersKeeperX0DetLayer:

Public Member Functions
MSLayer	layer (const DetLayer *layer) const override

const MSLayersAtAngle &	layers (float cotTheta) const override

	MSLayersKeeperX0DetLayer ()

	~MSLayersKeeperX0DetLayer () override=default

Public Member Functions inherited from MSLayersKeeper
virtual	~MSLayersKeeper ()

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

Constructor & Destructor Documentation

◆ MSLayersKeeperX0DetLayer()

MSLayersKeeperX0DetLayer::MSLayersKeeperX0DetLayer ( )

Definition at line 13 of file MSLayersKeeperX0DetLayer.cc.

References GeomDetEnumerators::barrel, PixelRecoLineRZ::cotLine(), caHitNtupletGeneratorKernels::if(), slope, and SiPixelPI::zero.

                                                    {
   //  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();
 }