#include <CTPPSToFDetector.h>

Public Member Functions
void	AddHit (double x, double y, double tof)

void	clear ()

	CTPPSToFDetector (int ncellx, int ncelly, std::vector< double > &cellw, double cellh, double pitchx, double pitchy, double pos, int res)

	CTPPSToFDetector (int ncellx, int ncelly, double cellw, double cellh, double pitchx, double pitchy, double pos, int res)

int	findCellId (double x, double y)

bool	get_CellCenter (int cell_id, double &x, double &y)

int	get_CellId (int idx)

int	get_CellMultiplicity ()

int	get_nHits_ ()

std::vector< double >	get_ToF (int cell)

int	getADC (int cell, int hit)

double	getHeight ()

int	getMultiplicityByCell (int cellid)

double	getPosition ()

double	getWidth ()

virtual	~CTPPSToFDetector ()

Private Types
typedef std::map< int, std::vector< double > >	theToFInfo_t

Private Attributes
std::vector< std::pair< double, double > >	cellColumn_

double	cellH_

std::vector< std::pair< double, double > >	cellRow_

std::vector< double >	cellW_

double	cellWq_

double	detH_

int	detId_

double	detPosition_

double	detW_

int	fToFResolution_

std::map< int, std::vector< int > >	nADC_

int	nCellX_

int	nCellY_

int	nHits_

double	pitchX_

double	pitchY_

std::map< int, std::vector< double > >	theToFInfo

Detailed Description

Definition at line 8 of file CTPPSToFDetector.h.

Member Typedef Documentation

typedef std::map<int, std::vector<double> > CTPPSToFDetector::theToFInfo_t

private

Definition at line 84 of file CTPPSToFDetector.h.

Constructor & Destructor Documentation

CTPPSToFDetector::CTPPSToFDetector	(	int	ncellx,
		int	ncelly,
		std::vector< double > &	cellw,
		double	cellh,
		double	pitchx,
		double	pitchy,
		double	pos,
		int	res
	)

Definition at line 4 of file CTPPSToFDetector.cc.

References cellColumn_, cellH_, cellRow_, cellW_, detH_, detPosition_, detW_, mps_fire::i, nCellX_, nCellY_, pitchX_, testProducerWithPsetDescEmpty_cfi::x1, and testProducerWithPsetDescEmpty_cfi::x2.

     : nCellX_(ncellx),
       nCellY_(ncelly),
       cellW_(cellw),
       cellH_(cellh),
       pitchX_(pitchx),
       pitchY_(pitchy),
       fToFResolution_(res),
       detPosition_(pos) {
   // the vertical positions starts from the negative(bottom) to the positive(top) corner
   // vector index points to the row number from below
   cellRow_.push_back(std::pair<double, double>(-cellH_ * 0.5, cellH_ * 0.5));
   cellColumn_.reserve(nCellX_);  // vector index points to the column number
   for (int i = 0; i < nCellX_; i++) {
     double x1 = 0., x2 = 0.;
     if (i == 0) {
       detW_ = pitchX_;
       x1 = -(detPosition_ + detW_);
     } else
       x1 = -detPosition_ + detW_;  //detPosition_ - shift the limit of a column depending on the detector position
     x2 = x1 - cellW_.at(i);
     detW_ += (x2 - x1) - pitchX_;
     cellColumn_.push_back(std::pair<double, double>(x1, x2));
   }
   //diamond geometry
   detH_ = nCellY_ * cellH_;
   detW_ = -detW_ - 2 * pitchX_;
 };

CTPPSToFDetector::CTPPSToFDetector	(	int	ncellx,
		int	ncelly,
		double	cellw,
		double	cellh,
		double	pitchx,
		double	pitchy,
		double	pos,
		int	res
	)

Definition at line 34 of file CTPPSToFDetector.cc.

References cellColumn_, cellH_, cellRow_, cellWq_, detH_, detPosition_, detW_, mps_fire::i, nCellX_, nCellY_, pitchX_, pitchY_, testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, testProducerWithPsetDescEmpty_cfi::y1, and testProducerWithPsetDescEmpty_cfi::y2.

     : nCellX_(ncellx),
       nCellY_(ncelly),
       cellWq_(cellwq),
       cellH_(cellh),
       pitchX_(pitchx),
       pitchY_(pitchy),
       fToFResolution_(res),
       detPosition_(pos) {
   //
   detW_ = nCellX_ * cellWq_ + (nCellX_ - 1) * pitchX_;
   detH_ = nCellY_ * cellH_ + (nCellY_ - 1) * pitchY_;
   // the vertical positions starts from the negative(bottom) to the positive(top) corner
   cellRow_.reserve(nCellY_);  // vector index points to the row number from below
   for (int i = 0; i < nCellY_; i++) {
     double y1 = cellH_ * (i - nCellY_ * 0.5) + pitchY_ * (i - (nCellY_ - 1) * 0.5);
     double y2 = y1 + cellH_;
     cellRow_.push_back(std::pair<double, double>(y1, y2));
   }
   cellColumn_.reserve(nCellX_);  // vector index points to the column number
   for (int i = 0; i < nCellX_; i++) {
     double x1 = -(cellWq_ * i + pitchX_ * i);
     x1 -= detPosition_;  // shift the limit of a column depending on the detector position
     double x2 = x1 - cellWq_;
     cellColumn_.push_back(std::pair<double, double>(x1, x2));
   }
 };

virtual CTPPSToFDetector::~CTPPSToFDetector ( )

inlinevirtual

Definition at line 20 of file CTPPSToFDetector.h.

20 {};

Member Function Documentation

void CTPPSToFDetector::AddHit	(	double	x,
		double	y,
		double	tof
	)

Definition at line 62 of file CTPPSToFDetector.cc.

References findCellId(), fToFResolution_, mps_fire::i, nADC_, nHits_, and theToFInfo.

Referenced by getADC().

                                                             {
   int cellid = findCellId(x, y);
   if (cellid == 0)
     return;
   if (theToFInfo.find(cellid) == theToFInfo.end())
     theToFInfo[cellid];  // add empty cell
   std::vector<double>* tofs = &(theToFInfo.find(cellid)->second);
   int ntof = tofs->size();
   int i = 0;
   double oneOverRes = 1.0 / fToFResolution_;
   for (; i < ntof; i++) {
     if (fabs(tofs->at(i) - tof) * oneOverRes < 3) {
       tofs->at(i) = (tofs->at(i) + tof) / 2.;
       nADC_.at(cellid).at(i)++;
       return;
     }
   }
   tofs->push_back(tof);  // no other ToF inside resolution found
   nHits_++;
   nADC_[cellid].push_back(1);
 }

void CTPPSToFDetector::clear ( void )

inline

Definition at line 56 of file CTPPSToFDetector.h.

References detId_, nHits_, and theToFInfo.

Referenced by Vispa.Views.WidgetView.WidgetView::closeEvent(), Vispa.Views.BoxDecayView.BoxDecayView::closeEvent(), Vispa.Share.FindAlgorithm.FindAlgorithm::findUsingFindDialog(), Vispa.Views.LineDecayView.LineDecayView::setDataObjects(), Vispa.Views.WidgetView.WidgetView::setDataObjects(), Vispa.Views.TreeView.TreeView::updateContent(), Vispa.Views.TableView.TableView::updateContent(), Vispa.Views.BoxDecayView.BoxDecayView::updateContent(), and Vispa.Views.PropertyView.PropertyView::updateContent().

                {
     detId_ = 0;
     nHits_ = 0;
     theToFInfo.clear();
   };

int CTPPSToFDetector::findCellId	(	double	x,
		double	y
	)

Definition at line 83 of file CTPPSToFDetector.cc.

References cellColumn_, cellRow_, and mps_fire::i.

Referenced by AddHit(), CTPPSFastTrackingProducer::FastReco(), getPosition(), and CTPPSRecHitProducer::produce().

                                                    {
   int y_idx, x_idx;
   // first, get the row number
   unsigned int i;
   unsigned int start_idx = 0;
   unsigned int end_idx = cellRow_.size();
   for (i = 0; i < cellRow_.size(); i++) {
     if (y >= cellRow_.at(i).first && y <= cellRow_.at(i).second)
       break;
   }
   if (i >= cellRow_.size())
     return 0;
   y_idx = i + 1;
   start_idx = 0;
   end_idx = cellColumn_.size();
   for (i = start_idx; i < end_idx; i++) {
     if (x <= cellColumn_.at(i).first && x > cellColumn_.at(i).second)
       break;
   }
   if (i >= end_idx)
     return 0;
   x_idx = i + 1 - start_idx;
   return 100 * y_idx + x_idx;
 }

bool CTPPSToFDetector::get_CellCenter	(	int	cell_id,
		double &	x,
		double &	y
	)

Definition at line 107 of file CTPPSToFDetector.cc.

References cellColumn_, cellRow_, and createfilelist::int.

Referenced by getPosition(), and CTPPSRecHitProducer::produce().

                                                                        {
   if (cell_id == 0)
     return false;
   //if(!isValidCellId(cell_id)) return 0;
   unsigned int y_idx = int(cell_id * 0.01);
   unsigned int x_idx = cell_id - y_idx * 100;
   x = (cellColumn_.at(x_idx - 1).first + cellColumn_.at(x_idx - 1).second) / 2.0;
   y = (cellRow_.at(y_idx - 1).first + cellRow_.at(y_idx - 1).second) / 2.0;
   return true;
 }

int CTPPSToFDetector::get_CellId ( int idx )

inline

Definition at line 27 of file CTPPSToFDetector.h.

References theToFInfo.

                           {
     if (idx >= (int)theToFInfo.size())
       return 0;
     std::map<int, std::vector<double> >::const_iterator it = theToFInfo.begin();
     std::advance(it, idx);
     return it->first;
   }

int CTPPSToFDetector::get_CellMultiplicity ( )

inline

Definition at line 34 of file CTPPSToFDetector.h.

References theToFInfo.

34 { return (int)theToFInfo.size(); }; // return the number of cells with hit

CTPPSToFDetector::theToFInfo

std::map< int, std::vector< double > > theToFInfo

Definition: CTPPSToFDetector.h:82

int CTPPSToFDetector::get_nHits_ ( )

inline

Definition at line 41 of file CTPPSToFDetector.h.

References nHits_.

41 { return nHits_; }; // return the total hit multiplicity (full det)

CTPPSToFDetector::nHits_

int nHits_

Definition: CTPPSToFDetector.h:79

std::vector<double> CTPPSToFDetector::get_ToF ( int cell )

inline

Definition at line 42 of file CTPPSToFDetector.h.

References theToFInfo.

                                       {
     if (!theToFInfo.count(cell))
       return std::vector<double>();
     return theToFInfo.at(cell);
   };  // return the tof of the given cell

int CTPPSToFDetector::getADC	(	int	cell,
		int	hit
	)

inline

Definition at line 47 of file CTPPSToFDetector.h.

References AddHit(), and nADC_.

                                 {
     if (!nADC_.count(cell))
       return 0;
     if ((int)nADC_.at(cell).size() < hit)
       return 0;
     return nADC_.at(cell).at(hit);
   }

double CTPPSToFDetector::getHeight ( )

inline

Definition at line 22 of file CTPPSToFDetector.h.

References detH_.

22 { return detH_; };

CTPPSToFDetector::detH_

double detH_

Definition: CTPPSToFDetector.h:75

int CTPPSToFDetector::getMultiplicityByCell ( int cellid )

inline

Definition at line 35 of file CTPPSToFDetector.h.

References theToFInfo.

                                         {
     if (!theToFInfo.count(cellid))
       return 0;
     return (int)theToFInfo.at(cellid).size();
   };  // return the hit multiplicity of the given cell