MuonPathAnalyzerInChamber Class Reference

#include <MuonPathAnalyzerInChamber.h>

Inheritance diagram for MuonPathAnalyzerInChamber:

Public Member Functions
int	bxTolerance (void)
void	finish () override
bool	hasPosRF (int wh, int sec)
void	initialise (const edm::EventSetup &iEventSetup) override
int	minHits4Fit (void)
cmsdt::MP_QUALITY	minQuality (void)
	MuonPathAnalyzerInChamber (const edm::ParameterSet &pset, edm::ConsumesCollector &iC)
void	run (edm::Event &iEvent, const edm::EventSetup &iEventSetup, MuonPathPtrs &inMpath, MuonPathPtrs &outMPath) override
void	run (edm::Event &iEvent, const edm::EventSetup &iEventSetup, MuonPathPtrs &inMpath, std::vector< cmsdt::metaPrimitive > &metaPrimitives) override
void	setBxTolerance (int t)
void	setChiSquareThreshold (float ch2Thr)
void	setMinHits4Fit (int h)
void	setMinimumQuality (cmsdt::MP_QUALITY q)
	~MuonPathAnalyzerInChamber () override
Public Member Functions inherited from MuonPathAnalyzer
	MuonPathAnalyzer (const edm::ParameterSet &pset, edm::ConsumesCollector &iC)
virtual	~MuonPathAnalyzer ()

Public Attributes
const DTGeometry *	dtGeo_
edm::ESGetToken< DTGeometry, MuonGeometryRecord >	dtGeomH
std::map< int, float >	shiftinfo_

Private Member Functions
void	analyze (MuonPathPtr &inMPath, MuonPathPtrs &outMPaths)
void	buildLateralities (MuonPathPtr &mpath)
void	calculateFitParameters (MuonPathPtr &mpath, TLateralities lat, int present_layer[NLayers])
void	evaluateQuality (MuonPathPtr &mPath)
void	setCellLayout (MuonPathPtr &mpath)
void	setLateralitiesInMP (MuonPathPtr &mpath, TLateralities lat)
void	setWirePosAndTimeInMP (MuonPathPtr &mpath)

Private Attributes
int	bxTolerance_
int	cellLayout_ [NLayers]
double	chi2Th_
float	chiSquareThreshold_
bool	debug_
std::vector< TLateralities >	lateralities_
std::vector< cmsdt::LATQ_TYPE >	latQuality_
short	minHits4Fit_
cmsdt::MP_QUALITY	minQuality_
edm::FileInPath	shift_filename_
int	totalNumValLateralities_

Detailed Description

Definition at line 17 of file MuonPathAnalyzerInChamber.h.

Constructor & Destructor Documentation

MuonPathAnalyzerInChamber()

MuonPathAnalyzerInChamber::MuonPathAnalyzerInChamber	(	const edm::ParameterSet &	pset,
		edm::ConsumesCollector &	iC
	)

Definition at line 11 of file MuonPathAnalyzerInChamber.cc.

    : MuonPathAnalyzer(pset, iC),
      debug_(pset.getUntrackedParameter<bool>("debug")),
      chi2Th_(pset.getUntrackedParameter<double>("chi2Th")),
      shift_filename_(pset.getParameter<edm::FileInPath>("shift_filename")),
      bxTolerance_(30),
      minQuality_(LOWQGHOST),
      chiSquareThreshold_(50),
      minHits4Fit_(pset.getUntrackedParameter<int>("minHits4Fit")) {
  // Obtention of parameters
 
  if (debug_)
    LogDebug("MuonPathAnalyzerInChamber") << "MuonPathAnalyzer: constructor";
 
  setChiSquareThreshold(chi2Th_ * 100.);
 
  //shift
  int rawId;
  std::ifstream ifin3(shift_filename_.fullPath());
  double shift;
  if (ifin3.fail()) {
    throw cms::Exception("Missing Input File")
        << "MuonPathAnalyzerInChamber::MuonPathAnalyzerInChamber() -  Cannot find " << shift_filename_.fullPath();
  }
  while (ifin3.good()) {
    ifin3 >> rawId >> shift;
    shiftinfo_[rawId] = shift;
  }
 
  dtGeomH = iC.esConsumes<DTGeometry, MuonGeometryRecord, edm::Transition::BeginRun>();
}

References edm::BeginRun, chi2Th_, debug_, dtGeomH, edm::ConsumesCollector::esConsumes(), Exception, edm::FileInPath::fullPath(), LogDebug, setChiSquareThreshold(), edm::shift, shift_filename_, and shiftinfo_.

~MuonPathAnalyzerInChamber()

MuonPathAnalyzerInChamber::~MuonPathAnalyzerInChamber ( )

override

Definition at line 43 of file MuonPathAnalyzerInChamber.cc.

                                                      {
  if (debug_)
    LogDebug("MuonPathAnalyzerInChamber") << "MuonPathAnalyzer: destructor";
}

References debug_, and LogDebug.

Member Function Documentation

analyze()

void MuonPathAnalyzerInChamber::analyze ( MuonPathPtr &  inMPath, MuonPathPtrs &  outMPaths  )

private

Definition at line 80 of file MuonPathAnalyzerInChamber.cc.

                                                                                    {
  if (debug_)
    LogDebug("MuonPathAnalyzerInChamber") << "DTp2:analyze \t\t\t\t starts";
 
  // Clone the analyzed object
  if (debug_)
    LogDebug("MuonPathAnalyzerInChamber") << inMPath->nprimitives();
  auto mPath = std::make_shared<MuonPath>(inMPath);
 
  if (debug_) {
    LogDebug("MuonPathAnalyzerInChamber") << "DTp2::analyze, looking at mPath: ";
    for (int i = 0; i < mPath->nprimitives(); i++)
      LogDebug("MuonPathAnalyzerInChamber")
          << mPath->primitive(i)->layerId() << " , " << mPath->primitive(i)->superLayerId() << " , "
          << mPath->primitive(i)->channelId() << " , " << mPath->primitive(i)->laterality();
  }
 
  if (debug_)
    LogDebug("MuonPathAnalyzerInChamber") << "DTp2:analyze \t\t\t\t\t is Analyzable? ";
  if (!mPath->isAnalyzable())
    return;
  if (debug_)
    LogDebug("MuonPathAnalyzerInChamber") << "DTp2:analyze \t\t\t\t\t yes it is analyzable " << mPath->isAnalyzable();
 
  // first of all, get info from primitives, so we can reduce the number of latereralities:
  buildLateralities(mPath);
  setWirePosAndTimeInMP(mPath);
 
  std::shared_ptr<MuonPath> mpAux;
  int bestI = -1;
  float best_chi2 = 99999.;
  for (int i = 0; i < totalNumValLateralities_; i++) {  // LOOP for all lateralities:
    if (debug_)
      LogDebug("MuonPathAnalyzerInChamber") << "DTp2:analyze \t\t\t\t\t Start with combination " << i;
    int NTotalHits = NUM_LAYERS_2SL;
    float xwire[NUM_LAYERS_2SL];
    int present_layer[NUM_LAYERS_2SL];
    for (int ii = 0; ii < 8; ii++) {
      xwire[ii] = mPath->xWirePos(ii);
      if (xwire[ii] == 0) {
        present_layer[ii] = 0;
        NTotalHits--;
      } else {
        present_layer[ii] = 1;
      }
    }
 
    while (NTotalHits >= minHits4Fit_) {
      mPath->setChiSquare(0);
      calculateFitParameters(mPath, lateralities_[i], present_layer);
      if (mPath->chiSquare() != 0)
        break;
      NTotalHits--;
    }
    if (mPath->chiSquare() > chiSquareThreshold_)
      continue;
 
    evaluateQuality(mPath);
 
    if (mPath->quality() < minQuality_)
      continue;
 
    double z = 0;
    double jm_x = (mPath->horizPos());
    int selected_Id = 0;
    for (int i = 0; i < mPath->nprimitives(); i++) {
      if (mPath->primitive(i)->isValidTime()) {
        selected_Id = mPath->primitive(i)->cameraId();
        mPath->setRawId(selected_Id);
        break;
      }
    }
    DTLayerId thisLId(selected_Id);
    if (thisLId.station() >= 3)
      z = Z_SHIFT_MB4;
 
    DTSuperLayerId MuonPathSLId(thisLId.wheel(), thisLId.station(), thisLId.sector(), thisLId.superLayer());
    GlobalPoint jm_x_cmssw_global = dtGeo_->chamber(MuonPathSLId)->toGlobal(LocalPoint(jm_x, 0., z));
    //jm_x is already extrapolated to the middle of the SL
    int thisec = MuonPathSLId.sector();
    if (thisec == 13)
      thisec = 4;
    if (thisec == 14)
      thisec = 10;
    double phi = jm_x_cmssw_global.phi() - PHI_CONV * (thisec - 1);
    double psi = atan(mPath->tanPhi());
    mPath->setPhi(jm_x_cmssw_global.phi() - PHI_CONV * (thisec - 1));
    mPath->setPhiB(hasPosRF(MuonPathSLId.wheel(), MuonPathSLId.sector()) ? psi - phi : -psi - phi);
 
    if (mPath->chiSquare() < best_chi2 && mPath->chiSquare() > 0) {
      mpAux = std::make_shared<MuonPath>(mPath);
      bestI = i;
      best_chi2 = mPath->chiSquare();
    }
  }
  if (mpAux != nullptr) {
    outMPath.push_back(std::move(mpAux));
    if (debug_)
      LogDebug("MuonPathAnalyzerInChamber")
          << "DTp2:analize \t\t\t\t\t Laterality " << bestI << " is the one with smaller chi2";
  } else {
    if (debug_)
      LogDebug("MuonPathAnalyzerInChamber")
          << "DTp2:analize \t\t\t\t\t No Laterality found with chi2 smaller than threshold";
  }
  if (debug_)
    LogDebug("MuonPathAnalyzerInChamber") << "DTp2:analize \t\t\t\t\t Ended working with this set of lateralities";
}

References buildLateralities(), calculateFitParameters(), DTGeometry::chamber(), chiSquareThreshold_, debug_, dtGeo_, evaluateQuality(), hasPosRF(), mps_fire::i, cuy::ii, lateralities_, LogDebug, minHits4Fit_, minQuality_, eostools::move(), cmsdt::NUM_LAYERS_2SL, phi, cmsdt::PHI_CONV, DTChamberId::sector(), setWirePosAndTimeInMP(), DTChamberId::station(), DTSuperLayerId::superLayer(), GeomDet::toGlobal(), totalNumValLateralities_, DTChamberId::wheel(), z, and cmsdt::Z_SHIFT_MB4.

Referenced by run().

buildLateralities()

void MuonPathAnalyzerInChamber::buildLateralities ( MuonPathPtr &  mpath )

private

For a combination of up to 8 cells, build all the lateralities to be tested,

Definition at line 210 of file MuonPathAnalyzerInChamber.cc.

                                                                    {
  if (debug_)
    LogDebug("MuonPathAnalyzerInChamber") << "MPAnalyzer::buildLateralities << setLateralitiesFromPrims ";
  mpath->setLateralCombFromPrimitives();
 
  totalNumValLateralities_ = 0;
  lateralities_.clear();
  latQuality_.clear();
 
  /* We generate all the possible laterality combinations compatible with the built 
     group in the previous step*/
  lateralities_.push_back(TLateralities());
  for (int ilat = 0; ilat < NLayers; ilat++) {
    // Get value from input
    LATERAL_CASES lr = (mpath->lateralComb())[ilat];
    if (debug_)
      LogDebug("MuonPathAnalyzerInChamber") << "[DEBUG_] Input[" << ilat << "]: " << lr;
 
    // If left/right fill number
    if (lr != NONE) {
      if (debug_)
        LogDebug("MuonPathAnalyzerInChamber") << "[DEBUG_]   - Adding it to " << lateralities_.size() << " lists...";
      for (unsigned int iall = 0; iall < lateralities_.size(); iall++) {
        lateralities_[iall][ilat] = lr;
      }
    }
    // both possibilites
    else {
      // Get the number of possible options now
      auto ncurrentoptions = lateralities_.size();
 
      // Duplicate them
      if (debug_)
        LogDebug("MuonPathAnalyzerInChamber") << "[DEBUG_]   - Duplicating " << ncurrentoptions << " lists...";
      copy(lateralities_.begin(), lateralities_.end(), back_inserter(lateralities_));
      if (debug_)
        LogDebug("MuonPathAnalyzerInChamber") << "[DEBUG_]   - Now we have " << lateralities_.size() << " lists...";
 
      // Asign LEFT to first ncurrentoptions and RIGHT to the last
      for (unsigned int iall = 0; iall < ncurrentoptions; iall++) {
        lateralities_[iall][ilat] = LEFT;
        lateralities_[iall + ncurrentoptions][ilat] = RIGHT;
      }
    }  // else
  }    // Iterate over input array
 
  totalNumValLateralities_ = (int)lateralities_.size();
  if (totalNumValLateralities_ > 128) {
    // ADD PROTECTION!
    LogDebug("MuonPathAnalyzerInChamber") << "[WARNING]: TOO MANY LATERALITIES TO CHECK !!";
    LogDebug("MuonPathAnalyzerInChamber") << "[WARNING]: skipping this muon";
    lateralities_.clear();
    latQuality_.clear();
    totalNumValLateralities_ = 0;
  }
 
  // Dump values
  if (debug_) {
    for (unsigned int iall = 0; iall < lateralities_.size(); iall++) {
      LogDebug("MuonPathAnalyzerInChamber") << iall << " -> [";
      for (int ilat = 0; ilat < NLayers; ilat++) {
        if (ilat != 0)
          LogDebug("MuonPathAnalyzerInChamber") << ",";
        LogDebug("MuonPathAnalyzerInChamber") << lateralities_[iall][ilat];
      }
      LogDebug("MuonPathAnalyzerInChamber") << "]";
    }
  }
}

References filterCSVwithJSON::copy, debug_, createfilelist::int, lateralities_, latQuality_, cmsdt::LEFT, LogDebug, NONE, cmsdt::RIGHT, and totalNumValLateralities_.

Referenced by analyze().

bxTolerance()

int MuonPathAnalyzerInChamber::bxTolerance ( void  )

inline

Definition at line 45 of file MuonPathAnalyzerInChamber.h.

{ return bxTolerance_; };

References bxTolerance_.

calculateFitParameters()

void MuonPathAnalyzerInChamber::calculateFitParameters ( MuonPathPtr &  mpath, TLateralities  lat, int  present_layer[NLayers]  )

private

Definition at line 336 of file MuonPathAnalyzerInChamber.cc.

                                                                                          {
  // First prepare mpath for fit:
  float xwire[NUM_LAYERS_2SL], zwire[NUM_LAYERS_2SL], tTDCvdrift[NUM_LAYERS_2SL];
  double b[NUM_LAYERS_2SL];
  for (int i = 0; i < 8; i++) {
    xwire[i] = mpath->xWirePos(i);
    zwire[i] = mpath->zWirePos(i);
    tTDCvdrift[i] = mpath->tWireTDC(i);
    b[i] = 1;
  }
 
 
  // fill hit position
  float xhit[NUM_LAYERS_2SL];
  for (int lay = 0; lay < 8; lay++) {
    if (debug_)
      LogDebug("MuonPathAnalyzerInChamber") << "In fitPerLat " << lay << " xwire " << xwire[lay] << " zwire "
                                            << zwire[lay] << " tTDCvdrift " << tTDCvdrift[lay];
    xhit[lay] = xwire[lay] + (-1 + 2 * laterality[lay]) * 1000 * tTDCvdrift[lay];
    if (debug_)
      LogDebug("MuonPathAnalyzerInChamber") << "In fitPerLat " << lay << " xhit " << xhit[lay];
  }
 
  //Proceed with calculation of fit parameters
  double cbscal = 0.0;
  double zbscal = 0.0;
  double czscal = 0.0;
  double bbscal = 0.0;
  double zzscal = 0.0;
  double ccscal = 0.0;
 
  for (int lay = 0; lay < 8; lay++) {
    if (present_layer[lay] == 0)
      continue;
    if (debug_)
      LogDebug("MuonPathAnalyzerInChamber")
          << " For layer " << lay + 1 << " xwire[lay] " << xwire[lay] << " zwire " << zwire[lay] << " b " << b[lay];
    if (debug_)
      LogDebug("MuonPathAnalyzerInChamber") << " xhit[lat][lay] " << xhit[lay];
    cbscal = (-1 + 2 * laterality[lay]) * b[lay] + cbscal;
    zbscal = zwire[lay] * b[lay] + zbscal;  //it actually does not depend on laterality
    czscal = (-1 + 2 * laterality[lay]) * zwire[lay] + czscal;
 
    bbscal = b[lay] * b[lay] + bbscal;          //it actually does not depend on laterality
    zzscal = zwire[lay] * zwire[lay] + zzscal;  //it actually does not depend on laterality
    ccscal = (-1 + 2 * laterality[lay]) * (-1 + 2 * laterality[lay]) + ccscal;
  }
 
  double cz = 0.0;
  double cb = 0.0;
  double zb = 0.0;
  double zc = 0.0;
  double bc = 0.0;
  double bz = 0.0;
 
  cz = (cbscal * zbscal - czscal * bbscal) / (zzscal * bbscal - zbscal * zbscal);
  cb = (czscal * zbscal - cbscal * zzscal) / (zzscal * bbscal - zbscal * zbscal);
 
  zb = (czscal * cbscal - zbscal * ccscal) / (bbscal * ccscal - cbscal * cbscal);
  zc = (zbscal * cbscal - czscal * bbscal) / (bbscal * ccscal - cbscal * cbscal);
 
  bc = (zbscal * czscal - cbscal * zzscal) / (ccscal * zzscal - czscal * czscal);
  bz = (cbscal * czscal - zbscal * ccscal) / (ccscal * zzscal - czscal * czscal);
 
  double c_tilde[NUM_LAYERS_2SL];
  double z_tilde[NUM_LAYERS_2SL];
  double b_tilde[NUM_LAYERS_2SL];
 
  for (int lay = 0; lay < 8; lay++) {
    if (present_layer[lay] == 0)
      continue;
    if (debug_)
      LogDebug("MuonPathAnalyzerInChamber")
          << " For layer " << lay + 1 << " xwire[lay] " << xwire[lay] << " zwire " << zwire[lay] << " b " << b[lay];
    c_tilde[lay] = (-1 + 2 * laterality[lay]) + cz * zwire[lay] + cb * b[lay];
    z_tilde[lay] = zwire[lay] + zb * b[lay] + zc * (-1 + 2 * laterality[lay]);
    b_tilde[lay] = b[lay] + bc * (-1 + 2 * laterality[lay]) + bz * zwire[lay];
  }
 
  //Calculate results per lat
  double xctilde = 0.0;
  double xztilde = 0.0;
  double xbtilde = 0.0;
  double ctildectilde = 0.0;
  double ztildeztilde = 0.0;
  double btildebtilde = 0.0;
 
  double rect0vdrift = 0.0;
  double recslope = 0.0;
  double recpos = 0.0;
 
  for (int lay = 0; lay < 8; lay++) {
    if (present_layer[lay] == 0)
      continue;
    xctilde = xhit[lay] * c_tilde[lay] + xctilde;
    ctildectilde = c_tilde[lay] * c_tilde[lay] + ctildectilde;
    xztilde = xhit[lay] * z_tilde[lay] + xztilde;
    ztildeztilde = z_tilde[lay] * z_tilde[lay] + ztildeztilde;
    xbtilde = xhit[lay] * b_tilde[lay] + xbtilde;
    btildebtilde = b_tilde[lay] * b_tilde[lay] + btildebtilde;
  }
 
  //Results for t0vdrift (BX), slope and position per lat
  rect0vdrift = xctilde / ctildectilde;
  recslope = xztilde / ztildeztilde;
  recpos = xbtilde / btildebtilde;
  if (debug_) {
    LogDebug("MuonPathAnalyzerInChamber") << " In fitPerLat Reconstructed values per lat "
                                          << " rect0vdrift " << rect0vdrift;
    LogDebug("MuonPathAnalyzerInChamber")
        << "rect0 " << rect0vdrift / DRIFT_SPEED << " recBX " << rect0vdrift / DRIFT_SPEED / 25 << " recslope "
        << recslope << " recpos " << recpos;
  }
 
  //Get t*v and residuals per layer, and chi2 per laterality
  double rectdriftvdrift[NUM_LAYERS_2SL];
  double recres[NUM_LAYERS_2SL];
  double recchi2 = 0.0;
  int sign_tdriftvdrift = {0};
  int incell_tdriftvdrift = {0};
  int physical_slope = {0};
 
  // Select the worst hit in order to get rid of it
  double maxDif = -1;
  int maxInt = -1;
 
  for (int lay = 0; lay < 8; lay++) {
    if (present_layer[lay] == 0)
      continue;
    rectdriftvdrift[lay] = tTDCvdrift[lay] - rect0vdrift / 1000;
    if (debug_)
      LogDebug("MuonPathAnalyzerInChamber") << rectdriftvdrift[lay];
    recres[lay] = xhit[lay] - zwire[lay] * recslope - b[lay] * recpos - (-1 + 2 * laterality[lay]) * rect0vdrift;
    if ((present_layer[lay] == 1) && (rectdriftvdrift[lay] < -0.1)) {
      sign_tdriftvdrift = -1;
      if (-0.1 - rectdriftvdrift[lay] > maxDif) {
        maxDif = -0.1 - rectdriftvdrift[lay];
        maxInt = lay;
      }
    }
    if ((present_layer[lay] == 1) && (abs(rectdriftvdrift[lay]) > 21.1)) {
      incell_tdriftvdrift = -1;  //Changed to 2.11 to account for resolution effects
      if (rectdriftvdrift[lay] - 21.1 > maxDif) {
        maxDif = rectdriftvdrift[lay] - 21.1;
        maxInt = lay;
      }
    }
  }
 
  if (fabs(recslope / 10) > 1.3)
    physical_slope = -1;
 
  if (physical_slope == -1 && debug_)
    LogDebug("MuonPathAnalyzerInChamber") << "Combination with UNPHYSICAL slope ";
  if (sign_tdriftvdrift == -1 && debug_)
    LogDebug("MuonPathAnalyzerInChamber") << "Combination with negative tdrift-vdrift ";
  if (incell_tdriftvdrift == -1 && debug_)
    LogDebug("MuonPathAnalyzerInChamber") << "Combination with tdrift-vdrift larger than half cell ";
 
  for (int lay = 0; lay < 8; lay++) {
    if (present_layer[lay] == 0)
      continue;
    recchi2 = recres[lay] * recres[lay] + recchi2;
  }
  if (debug_)
    LogDebug("MuonPathAnalyzerInChamber")
        << "In fitPerLat Chi2 " << recchi2 << " with sign " << sign_tdriftvdrift << " within cell "
        << incell_tdriftvdrift << " physical_slope " << physical_slope;
 
  //LATERALITY IS NOT VALID
  if (true && maxInt != -1) {
    present_layer[maxInt] = 0;
    if (debug_)
      LogDebug("MuonPathAnalyzerInChamber") << "We get rid of hit in layer " << maxInt;
  }
 
  // LATERALITY IS VALID...
  if (!(sign_tdriftvdrift == -1) && !(incell_tdriftvdrift == -1) && !(physical_slope == -1)) {
    mpath->setBxTimeValue((rect0vdrift / DRIFT_SPEED) / 1000);
    mpath->setTanPhi(-1 * recslope / 10);
    mpath->setHorizPos(recpos / 10000);
    mpath->setChiSquare(recchi2 / 100000000);
    setLateralitiesInMP(mpath, laterality);
    if (debug_)
      LogDebug("MuonPathAnalyzerInChamber")
          << "In fitPerLat "
          << "t0 " << mpath->bxTimeValue() << " slope " << mpath->tanPhi() << " pos " << mpath->horizPos() << " chi2 "
          << mpath->chiSquare() << " rawId " << mpath->rawId();
  }
}

References funct::abs(), b, debug_, cmsdt::DRIFT_SPEED, mps_fire::i, LogDebug, cmsdt::NUM_LAYERS_2SL, and setLateralitiesInMP().

Referenced by analyze().

evaluateQuality()

void MuonPathAnalyzerInChamber::evaluateQuality ( MuonPathPtr &  mPath )

private

Definition at line 529 of file MuonPathAnalyzerInChamber.cc.

                                                                  {
  mPath->setQuality(NOPATH);
 
  if (mPath->nprimitivesUp() >= 4 && mPath->nprimitivesDown() >= 4) {
    mPath->setQuality(HIGHHIGHQ);
  } else if ((mPath->nprimitivesUp() == 4 && mPath->nprimitivesDown() == 3) ||
             (mPath->nprimitivesUp() == 3 && mPath->nprimitivesDown() == 4)) {
    mPath->setQuality(HIGHLOWQ);
  } else if ((mPath->nprimitivesUp() == 4 && mPath->nprimitivesDown() <= 2 && mPath->nprimitivesDown() > 0) ||
             (mPath->nprimitivesUp() <= 2 && mPath->nprimitivesUp() > 0 && mPath->nprimitivesDown() == 4)) {
    mPath->setQuality(CHIGHQ);
  } else if ((mPath->nprimitivesUp() == 3 && mPath->nprimitivesDown() == 3)) {
    mPath->setQuality(LOWLOWQ);
  } else if ((mPath->nprimitivesUp() == 3 && mPath->nprimitivesDown() <= 2 && mPath->nprimitivesDown() > 0) ||
             (mPath->nprimitivesUp() <= 2 && mPath->nprimitivesUp() > 0 && mPath->nprimitivesDown() == 3) ||
             (mPath->nprimitivesUp() == 2 && mPath->nprimitivesDown() == 2)) {
    mPath->setQuality(CLOWQ);
  } else if (mPath->nprimitivesUp() >= 4 || mPath->nprimitivesDown() >= 4) {
    mPath->setQuality(HIGHQ);
  } else if (mPath->nprimitivesUp() == 3 || mPath->nprimitivesDown() == 3) {
    mPath->setQuality(LOWQ);
  }
}

References cmsdt::CHIGHQ, cmsdt::CLOWQ, cmsdt::HIGHHIGHQ, cmsdt::HIGHLOWQ, cmsdt::HIGHQ, cmsdt::LOWLOWQ, cmsdt::LOWQ, and cmsdt::NOPATH.

Referenced by analyze().

finish()

void MuonPathAnalyzerInChamber::finish ( )

overridevirtual

Reimplemented from MuonPathAnalyzer.

Definition at line 72 of file MuonPathAnalyzerInChamber.cc.

                                       {
  if (debug_)
    LogDebug("MuonPathAnalyzerInChamber") << "MuonPathAnalyzer: finish";
};

References debug_, and LogDebug.

Referenced by progressbar.ProgressBar::__next__().

hasPosRF()

bool MuonPathAnalyzerInChamber::hasPosRF ( int  wh, int  sec  )

inline

Definition at line 49 of file MuonPathAnalyzerInChamber.h.

{ return wh > 0 || (wh == 0 && sec % 4 > 1); };

References fileinputsource_cfi::sec.

Referenced by analyze().

initialise()

void MuonPathAnalyzerInChamber::initialise ( const edm::EventSetup &  iEventSetup )

overridevirtual

Reimplemented from MuonPathAnalyzer.

Definition at line 51 of file MuonPathAnalyzerInChamber.cc.

                                                                           {
  if (debug_)
    LogDebug("MuonPathAnalyzerInChamber") << "MuonPathAnalyzerInChamber::initialiase";
 
  const MuonGeometryRecord &geom = iEventSetup.get<MuonGeometryRecord>();
  dtGeo_ = &geom.get(dtGeomH);
}

References debug_, dtGeo_, dtGeomH, relativeConstraints::geom, edm::EventSetup::get(), and LogDebug.

minHits4Fit()

int MuonPathAnalyzerInChamber::minHits4Fit ( void  )

inline

Definition at line 46 of file MuonPathAnalyzerInChamber.h.

{ return minHits4Fit_; };

References minHits4Fit_.

minQuality()

cmsdt::MP_QUALITY MuonPathAnalyzerInChamber::minQuality ( void  )

inline

Definition at line 47 of file MuonPathAnalyzerInChamber.h.

{ return minQuality_; };

References minQuality_.

run() [1/2]

void MuonPathAnalyzerInChamber::run ( edm::Event &  iEvent, const edm::EventSetup &  iEventSetup, MuonPathPtrs &  inMpath, MuonPathPtrs &  outMPath  )

overridevirtual

Implements MuonPathAnalyzer.

Definition at line 59 of file MuonPathAnalyzerInChamber.cc.

                                                                {
  if (debug_)
    LogDebug("MuonPathAnalyzerInChamber") << "MuonPathAnalyzerInChamber: run";
 
  // fit per SL (need to allow for multiple outputs for a single mpath)
  for (auto muonpath = muonpaths.begin(); muonpath != muonpaths.end(); ++muonpath) {
    analyze(*muonpath, outmuonpaths);
  }
}

References analyze(), debug_, and LogDebug.

run() [2/2]

void MuonPathAnalyzerInChamber::run ( edm::Event &  iEvent, const edm::EventSetup &  iEventSetup, MuonPathPtrs &  inMpath, std::vector< cmsdt::metaPrimitive > &  metaPrimitives  )

inlineoverridevirtual

Implements MuonPathAnalyzer.

Definition at line 25 of file MuonPathAnalyzerInChamber.h.

                                                                   {}

setBxTolerance()

void MuonPathAnalyzerInChamber::setBxTolerance ( int  t )

inline

Definition at line 37 of file MuonPathAnalyzerInChamber.h.

{ bxTolerance_ = t; };

References bxTolerance_, and submitPVValidationJobs::t.

setCellLayout()

void MuonPathAnalyzerInChamber::setCellLayout ( MuonPathPtr &  mpath )

private

Definition at line 189 of file MuonPathAnalyzerInChamber.cc.

                                                                {
  for (int i = 0; i <= mpath->nprimitives(); i++) {
    if (mpath->primitive(i)->isValidTime())
      cellLayout_[i] = mpath->primitive(i)->channelId();
    else
      cellLayout_[i] = -99;
  }
 
  // putting info back into the mpath:
  mpath->setCellHorizontalLayout(cellLayout_);
  for (int i = 0; i <= mpath->nprimitives(); i++) {
    if (cellLayout_[i] >= 0) {
      mpath->setBaseChannelId(cellLayout_[i]);
      break;
    }
  }
}

References cellLayout_, and mps_fire::i.

setChiSquareThreshold()

void MuonPathAnalyzerInChamber::setChiSquareThreshold ( float  ch2Thr )

inline

Definition at line 39 of file MuonPathAnalyzerInChamber.h.

{ chiSquareThreshold_ = ch2Thr; };

References chiSquareThreshold_.

Referenced by MuonPathAnalyzerInChamber().

setLateralitiesInMP()

void MuonPathAnalyzerInChamber::setLateralitiesInMP ( MuonPathPtr &  mpath, TLateralities  lat  )

private

Definition at line 279 of file MuonPathAnalyzerInChamber.cc.

                                                                                         {
  LATERAL_CASES tmp[NUM_LAYERS_2SL];
  for (int i = 0; i < 8; i++)
    tmp[i] = lat[i];
 
  mpath->setLateralComb(tmp);
}

References mps_fire::i, cmsdt::NUM_LAYERS_2SL, and createJobs::tmp.

Referenced by calculateFitParameters().

setMinHits4Fit()

void MuonPathAnalyzerInChamber::setMinHits4Fit ( int  h )

inline

Definition at line 38 of file MuonPathAnalyzerInChamber.h.

{ minHits4Fit_ = h; };

References h, and minHits4Fit_.

setMinimumQuality()

void MuonPathAnalyzerInChamber::setMinimumQuality ( cmsdt::MP_QUALITY  q )

inline

Definition at line 40 of file MuonPathAnalyzerInChamber.h.

                                            {
    if (minQuality_ >= cmsdt::LOWQGHOST)
      minQuality_ = q;
  };

References cmsdt::LOWQGHOST, minQuality_, and submitPVResolutionJobs::q.

setWirePosAndTimeInMP()

void MuonPathAnalyzerInChamber::setWirePosAndTimeInMP ( MuonPathPtr &  mpath )

private

Definition at line 286 of file MuonPathAnalyzerInChamber.cc.

                                                                        {
  int selected_Id = 0;
  for (int i = 0; i < mpath->nprimitives(); i++) {
    if (mpath->primitive(i)->isValidTime()) {
      selected_Id = mpath->primitive(i)->cameraId();
      mpath->setRawId(selected_Id);
      break;
    }
  }
  DTLayerId thisLId(selected_Id);
  DTChamberId chId(thisLId.wheel(), thisLId.station(), thisLId.sector());
  if (debug_)
    LogDebug("MuonPathAnalyzerInChamber")
        << "Id " << chId.rawId() << " Wh " << chId.wheel() << " St " << chId.station() << " Se " << chId.sector();
  mpath->setRawId(chId.rawId());
 
  DTSuperLayerId MuonPathSLId1(thisLId.wheel(), thisLId.station(), thisLId.sector(), 1);
  DTSuperLayerId MuonPathSLId3(thisLId.wheel(), thisLId.station(), thisLId.sector(), 3);
  DTWireId wireId1(MuonPathSLId1, 2, 1);
  DTWireId wireId3(MuonPathSLId3, 2, 1);
 
  if (debug_)
    LogDebug("MuonPathAnalyzerInChamber")
        << "shift1=" << shiftinfo_[wireId1.rawId()] << " shift3=" << shiftinfo_[wireId3.rawId()];
 
  float delta = 42000;                                                                      //um
  float zwire[NUM_LAYERS_2SL] = {-13.7, -12.4, -11.1, -9.8002, 9.79999, 11.1, 12.4, 13.7};  // mm
  for (int i = 0; i < mpath->nprimitives(); i++) {
    if (mpath->primitive(i)->isValidTime()) {
      if (i < 4)
        mpath->setXWirePos(10000 * shiftinfo_[wireId1.rawId()] +
                               (mpath->primitive(i)->channelId() + 0.5 * (double)((i + 1) % 2)) * delta,
                           i);
      if (i >= 4)
        mpath->setXWirePos(10000 * shiftinfo_[wireId3.rawId()] +
                               (mpath->primitive(i)->channelId() + 0.5 * (double)((i + 1) % 2)) * delta,
                           i);
      mpath->setZWirePos(zwire[i] * 1000, i);  // in um
      mpath->setTWireTDC(mpath->primitive(i)->tdcTimeStamp() * DRIFT_SPEED, i);
    } else {
      mpath->setXWirePos(0., i);
      mpath->setZWirePos(0., i);
      mpath->setTWireTDC(-1 * DRIFT_SPEED, i);
    }
    if (debug_)
      LogDebug("MuonPathAnalyzerInChamber") << mpath->primitive(i)->tdcTimeStamp() << " ";
  }
  if (debug_)
    LogDebug("MuonPathAnalyzerInChamber");
}

References debug_, dumpMFGeometry_cfg::delta, cmsdt::DRIFT_SPEED, mps_fire::i, LogDebug, cmsdt::NUM_LAYERS_2SL, DetId::rawId(), DTChamberId::sector(), shiftinfo_, DTChamberId::station(), and DTChamberId::wheel().

Referenced by analyze().

Member Data Documentation

bxTolerance_

int MuonPathAnalyzerInChamber::bxTolerance_

private

Definition at line 76 of file MuonPathAnalyzerInChamber.h.

Referenced by bxTolerance(), and setBxTolerance().

cellLayout_

int MuonPathAnalyzerInChamber::cellLayout_[NLayers]

private

Definition at line 80 of file MuonPathAnalyzerInChamber.h.

Referenced by setCellLayout().

chi2Th_

double MuonPathAnalyzerInChamber::chi2Th_

private

Definition at line 74 of file MuonPathAnalyzerInChamber.h.

Referenced by MuonPathAnalyzerInChamber().

chiSquareThreshold_

float MuonPathAnalyzerInChamber::chiSquareThreshold_

private

Definition at line 78 of file MuonPathAnalyzerInChamber.h.

Referenced by analyze(), and setChiSquareThreshold().

debug_

bool MuonPathAnalyzerInChamber::debug_

private

Definition at line 73 of file MuonPathAnalyzerInChamber.h.

Referenced by analyze(), buildLateralities(), calculateFitParameters(), finish(), initialise(), MuonPathAnalyzerInChamber(), run(), setWirePosAndTimeInMP(), and ~MuonPathAnalyzerInChamber().

dtGeo_

const DTGeometry* MuonPathAnalyzerInChamber::dtGeo_

Definition at line 49 of file MuonPathAnalyzerInChamber.h.

Referenced by analyze(), and initialise().

dtGeomH

edm::ESGetToken<DTGeometry, MuonGeometryRecord> MuonPathAnalyzerInChamber::dtGeomH

Definition at line 53 of file MuonPathAnalyzerInChamber.h.

Referenced by initialise(), and MuonPathAnalyzerInChamber().

lateralities_

std::vector<TLateralities> MuonPathAnalyzerInChamber::lateralities_

private

Definition at line 70 of file MuonPathAnalyzerInChamber.h.

Referenced by analyze(), and buildLateralities().

latQuality_

std::vector<cmsdt::LATQ_TYPE> MuonPathAnalyzerInChamber::latQuality_

private

Definition at line 71 of file MuonPathAnalyzerInChamber.h.

Referenced by buildLateralities().

minHits4Fit_

short MuonPathAnalyzerInChamber::minHits4Fit_

private

Definition at line 79 of file MuonPathAnalyzerInChamber.h.

Referenced by analyze(), minHits4Fit(), and setMinHits4Fit().

minQuality_

cmsdt::MP_QUALITY MuonPathAnalyzerInChamber::minQuality_

private

Definition at line 77 of file MuonPathAnalyzerInChamber.h.

Referenced by analyze(), minQuality(), and setMinimumQuality().

shift_filename_

edm::FileInPath MuonPathAnalyzerInChamber::shift_filename_

private

Definition at line 75 of file MuonPathAnalyzerInChamber.h.

Referenced by MuonPathAnalyzerInChamber().

shiftinfo_

std::map<int, float> MuonPathAnalyzerInChamber::shiftinfo_

Definition at line 56 of file MuonPathAnalyzerInChamber.h.

Referenced by MuonPathAnalyzerInChamber(), and setWirePosAndTimeInMP().

totalNumValLateralities_

int MuonPathAnalyzerInChamber::totalNumValLateralities_

private

Definition at line 69 of file MuonPathAnalyzerInChamber.h.

Referenced by analyze(), and buildLateralities().