EcalTBHodoscopeRecInfoAlgo Class Reference

#include <EcalTBHodoscopeRecInfoAlgo.h>

Classes
class	BeamTrack
	Class to hold track information. More...

Public Member Functions
	EcalTBHodoscopeRecInfoAlgo ()
	EcalTBHodoscopeRecInfoAlgo (int fitMethod, const std::vector< double > &planeShift, const std::vector< double > &zPosition)
EcalTBHodoscopeRecInfo	reconstruct (const EcalTBHodoscopeRawInfo &hodoscopeRawInfo) const

Private Member Functions
void	clusterPos (float &x, float &xQuality, const int &ipl, const int &xclus, const int &wclus) const
void	fitHodo (float &x, float &xQuality, const int &ipl, const int &nclus, const std::vector< int > &xclus, const std::vector< int > &wclus) const
void	fitLine (float &x, float &xSlope, float &xQuality, const int &ipl1, const int &nclus1, const std::vector< int > &xclus1, const std::vector< int > &wclus1, const int &ipl2, const int &nclus2, const std::vector< int > &xclus2, const std::vector< int > &wclus2) const

Private Attributes
int	fitMethod_
EcalTBHodoscopeGeometry	myGeometry_
std::vector< double >	planeShift_
std::vector< double >	zPosition_

Detailed Description

Definition at line 11 of file EcalTBHodoscopeRecInfoAlgo.h.

Constructor & Destructor Documentation

EcalTBHodoscopeRecInfoAlgo() [1/2]

EcalTBHodoscopeRecInfoAlgo::EcalTBHodoscopeRecInfoAlgo

(

)

EcalTBHodoscopeRecInfoAlgo() [2/2]

EcalTBHodoscopeRecInfoAlgo::EcalTBHodoscopeRecInfoAlgo ( int  fitMethod, const std::vector< double > &  planeShift, const std::vector< double > &  zPosition  )

explicit

Definition at line 6 of file EcalTBHodoscopeRecInfoAlgo.cc.

    : fitMethod_(fitMethod), planeShift_(planeShift), zPosition_(zPosition), myGeometry_() {}

Member Function Documentation

clusterPos()

void EcalTBHodoscopeRecInfoAlgo::clusterPos ( float &  x, float &  xQuality, const int &  ipl, const int &  xclus, const int &  wclus  ) const

private

Definition at line 11 of file EcalTBHodoscopeRecInfoAlgo.cc.

References EcalTBHodoscopeGeometry::getFibreLp(), EcalTBHodoscopeGeometry::getFibreRp(), myGeometry_, planeShift_, ApeEstimator_cff::width, and x.

Referenced by fitHodo(), and fitLine().

                                                                                         {
  if (width == 1) {
    // Single fiber
    x = (myGeometry_.getFibreLp(ipl, xclus) + myGeometry_.getFibreRp(ipl, xclus)) / 2.0 - planeShift_[ipl];
    xQuality = (myGeometry_.getFibreRp(ipl, xclus) - myGeometry_.getFibreLp(ipl, xclus));
  } else if (width == 2) {
    // Two half overlapped fibers
    x = (myGeometry_.getFibreLp(ipl, xclus + 1) + myGeometry_.getFibreRp(ipl, xclus)) / 2.0 - planeShift_[ipl];
    xQuality = (myGeometry_.getFibreRp(ipl, xclus) - myGeometry_.getFibreLp(ipl, xclus + 1));
  } else {
    // More then two fibers case
    x = (myGeometry_.getFibreLp(ipl, xclus) + myGeometry_.getFibreRp(ipl, xclus + width - 1)) / 2.0 - planeShift_[ipl];
    xQuality = (myGeometry_.getFibreRp(ipl, xclus + width - 1) - myGeometry_.getFibreLp(ipl, xclus));
  }
}

fitHodo()

void EcalTBHodoscopeRecInfoAlgo::fitHodo ( float &  x, float &  xQuality, const int &  ipl, const int &  nclus, const std::vector< int > &  xclus, const std::vector< int > &  wclus  ) const

private

Definition at line 28 of file EcalTBHodoscopeRecInfoAlgo.cc.

References clusterPos(), and x.

Referenced by fitLine(), and reconstruct().

                                                                            {
  if (nclus == 1) {
    // Fill real x as mean position inside the cluster
    // Quality - width of cluster
    // To calculate sigma one can do sigma=sqrt(Quality**2/12.0)
    clusterPos(x, xQuality, ipl, xclus[0], wclus[0]);
  } else {
    xQuality = -10 - nclus;
  }
}

fitLine()

void EcalTBHodoscopeRecInfoAlgo::fitLine ( float &  x, float &  xSlope, float &  xQuality, const int &  ipl1, const int &  nclus1, const std::vector< int > &  xclus1, const std::vector< int > &  wclus1, const int &  ipl2, const int &  nclus2, const std::vector< int > &  xclus2, const std::vector< int > &  wclus2  ) const

private

Definition at line 44 of file EcalTBHodoscopeRecInfoAlgo.cc.

References clusterPos(), fitHodo(), testProducerWithPsetDescEmpty_cfi::i1, testProducerWithPsetDescEmpty_cfi::i2, tracks, x, testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, and zPosition_.

Referenced by reconstruct().

                                                                             {
  if (nclus1 == 0) {  // Fit with one plane
    fitHodo(x, xQuality, ipl2, nclus2, xclus2, wclus2);
    xSlope = 0.0;  //?? Should we put another number indicating that is not fitted
    return;
  }
  if (nclus2 == 0) {  // Fit with one plane
    fitHodo(x, xQuality, ipl1, nclus1, xclus1, wclus1);
    xSlope = 0.0;  //?? Should we put another number indicating that is not fitted
    return;
  }

  // We have clusters in both planes

  float x1, x2, xQ1, xQ2;
  float xs, x0, xq;

  std::list<BeamTrack> tracks;

  for (int i1 = 0; i1 < nclus1; i1++) {
    for (int i2 = 0; i2 < nclus2; i2++) {
      clusterPos(x1, xQ1, ipl1, xclus1[i1], wclus1[i1]);
      clusterPos(x2, xQ2, ipl2, xclus2[i2], wclus2[i2]);

      xs = (x2 - x1) / (zPosition_[ipl2] - zPosition_[ipl1]);               // slope
      x0 = ((x2 + x1) - xs * (zPosition_[ipl2] + zPosition_[ipl1])) / 2.0;  // x0
      xq = (xQ1 + xQ2) / 2.0;                                               // Quality, how i can do better ?
      tracks.push_back(BeamTrack(x0, xs, xq));
    }
  }

  // find track with minimal slope
  tracks.sort();

  // Return results
  x = tracks.begin()->x;
  xSlope = tracks.begin()->xS;
  xQuality = tracks.begin()->xQ;
}

reconstruct()

EcalTBHodoscopeRecInfo EcalTBHodoscopeRecInfoAlgo::reconstruct

(

const EcalTBHodoscopeRawInfo &

hodoscopeRawInfo

)

const

Fit 0

Fit 1

Fit 2

Definition at line 94 of file EcalTBHodoscopeRecInfoAlgo.cc.

References first, fitHodo(), fitLine(), fitMethod_, EcalTBHodoscopeGeometry::getNFibres(), EcalTBHodoscopeGeometry::getNPlanes(), dqmdumpme::last, myGeometry_, cms::cuda::nh, nhits, x, and y.

Referenced by EcalTBHodoscopeRecInfoProducer::produce().

                                                                                                                   {
  // Reset Hodo data
  float x, y = -100.0;
  float xSlope, ySlope = 0.0;
  float xQuality, yQuality = -100.0;

  int nclus[4];
  std::vector<int> xclus[4];
  std::vector<int> wclus[4];

  for (int ipl = 0; ipl < myGeometry_.getNPlanes(); ipl++) {
    int nhits = hodoscopeRawInfo[ipl].numberOfFiredHits();
    // Finding clusters
    nclus[ipl] = 0;
    if (nhits > 0) {
      int nh = nhits;
      int first = 0;
      int last = 0;
      while (nh > 0) {
        while (hodoscopeRawInfo[ipl][first] == 0)
          first++;  // start
        last = first + 1;
        nh--;
        do {
          while (last < myGeometry_.getNFibres() && hodoscopeRawInfo[ipl][last]) {
            last++;
            nh--;
          }                                                                              //end
          if (last + 1 < myGeometry_.getNFibres() && hodoscopeRawInfo[ipl][last + 1]) {  //Skip 1 fibre hole
            last += 2;
            nh--;
            //std::cout << "Skip fibre " << ipl << " " << first << " "<< last << std::endl;
          } else {
            break;
          }
        } while (nh > 0 && last < myGeometry_.getNFibres());
        wclus[ipl].push_back(last - first);
        xclus[ipl].push_back(first);  // Left edge !!!
        nclus[ipl]++;

        first = last + 1;
      }
    }
    //    printClusters(ipl);
  }

  // Straight line fit for one axis
  if (fitMethod_ == 0) {
    fitLine(x,
            xSlope,
            xQuality,
            0,
            nclus[0],
            xclus[0],
            wclus[0],  // X1
            2,
            nclus[2],
            xclus[2],
            wclus[2]);  // X2
    fitLine(y,
            ySlope,
            yQuality,
            1,
            nclus[1],
            xclus[1],
            wclus[1],  // Y1
            3,
            nclus[3],
            xclus[3],
            wclus[3]);  // Y2
  } else if (fitMethod_ == 1) {
    // x1 and y2 hodoscope
    fitHodo(x, xQuality, 0, nclus[0], xclus[0], wclus[0]);  // X1
    //   if ( xQuality[1] < 0.0 ) {
    //     printFibres(0);
    //     printClusters(0);
    //   }
    fitHodo(y, yQuality, 1, nclus[1], xclus[1], wclus[1]);  // Y1
    //   if ( yQuality[1] < 0.0 ) {
    //     printFibres(1);
    //     printClusters(1);
    //   }
  } else if (fitMethod_ == 2) {
    //x2 and y2 hodoscope
    fitHodo(x, xQuality, 2, nclus[2], xclus[2], wclus[2]);  // X2
    //   if ( xQuality[2] < 0.0 ) {
    //     printFibres(2);
    //     printClusters(2);
    //   }
    fitHodo(y, yQuality, 3, nclus[3], xclus[3], wclus[3]);  // Y2
    //   if ( yQuality[2] < 0.0 ) {
    //     printFibres(3);
    //     printClusters(3);
    //   }
  }

  return EcalTBHodoscopeRecInfo(x, y, xSlope, ySlope, xQuality, yQuality);
}

Member Data Documentation

fitMethod_

int EcalTBHodoscopeRecInfoAlgo::fitMethod_

private

Definition at line 62 of file EcalTBHodoscopeRecInfoAlgo.h.

Referenced by reconstruct().

myGeometry_

EcalTBHodoscopeGeometry EcalTBHodoscopeRecInfoAlgo::myGeometry_

private

Definition at line 68 of file EcalTBHodoscopeRecInfoAlgo.h.

Referenced by clusterPos(), and reconstruct().

planeShift_

std::vector<double> EcalTBHodoscopeRecInfoAlgo::planeShift_

private

Definition at line 64 of file EcalTBHodoscopeRecInfoAlgo.h.

Referenced by clusterPos().

zPosition_

std::vector<double> EcalTBHodoscopeRecInfoAlgo::zPosition_

private

Definition at line 65 of file EcalTBHodoscopeRecInfoAlgo.h.

Referenced by fitLine().