MuRodBarrelLayer Class Reference

#include <MuRodBarrelLayer.h>

Inheritance diagram for MuRodBarrelLayer:

Public Member Functions
const std::vector< const GeomDet * > &	basicComponents () const override
std::vector< DetWithState >	compatibleDets (const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const override
const std::vector< const GeometricSearchDet * > &	components () const override
std::vector< DetGroup >	groupedCompatibleDets (const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const override
	MuRodBarrelLayer (std::vector< const DetRod * > &rods)
	Constructor, takes ownership of pointers. More...
virtual const std::vector< const DetRod * > &	rods () const
	Return the vector of rods. More...
SubDetector	subDetector () const override
	~MuRodBarrelLayer () override

Private Member Functions
float	xError (const TrajectoryStateOnSurface &tsos, const MeasurementEstimator &est) const

Private Attributes
bool	isOverlapping
std::vector< const GeomDet * >	theBasicComps
BaseBinFinder< double > *	theBinFinder
std::vector< const GeometricSearchDet * >	theComponents
std::vector< const DetRod * >	theRods

Detailed Description

A cylinder composed of rods. Represents barrel muon DT/RPC stations.

Author

N. Amapane - INFN Torino

Definition at line 17 of file MuRodBarrelLayer.h.

Constructor & Destructor Documentation

MuRodBarrelLayer::MuRodBarrelLayer

(

std::vector< const DetRod * > &

rods

)

Constructor, takes ownership of pointers.

Definition at line 26 of file MuRodBarrelLayer.cc.

References basicComponents(), popcon2dropbox::copy(), BarrelDetLayer::initialize(), isOverlapping, PhiBorderFinder::isPhiOverlapping(), PhiBorderFinder::isPhiPeriodic(), LogTrace, metname, precomputed_value_sort(), AlCaHLTBitMon_QueryRunRegistry::string, theBasicComps, theBinFinder, theComponents, and theRods.

                                                              :
  RodBarrelLayer(false),
  theRods(rods),
  theBinFinder(nullptr),
  isOverlapping(false)
{
  // Sort rods in phi
  precomputed_value_sort(theRods.begin(),theRods.end(), geomsort::ExtractPhi<DetRod,float>());

  theComponents.reserve(theRods.size());
  std::copy(theRods.begin(),theRods.end(),back_inserter(theComponents));

  const std::string metname = "Muon|RecoMuon|RecoMuonDetLayers|MuRodBarrelLayer";

  // Cache chamber pointers (the basic components_)
  for (vector<const DetRod*>::const_iterator it=rods.begin();
       it!=rods.end(); it++) {
    vector<const GeomDet*> tmp2 = (*it)->basicComponents();
    theBasicComps.insert(theBasicComps.end(),tmp2.begin(),tmp2.end());
  }

  // Initialize the binfinder
  PhiBorderFinder bf(theRods);
  isOverlapping = bf.isPhiOverlapping();

  if ( bf.isPhiPeriodic() ) { 
    theBinFinder = new PeriodicBinFinderInPhi<double>
    (theRods.front()->position().phi(),theRods.size());
  } else {
    theBinFinder = new GeneralBinFinderInPhi<double>(bf);
  }

  // Compute the layer's surface and bounds (from the components())
  BarrelDetLayer::initialize(); 
  
  LogTrace(metname) << "Constructing MuRodBarrelLayer: "
            << basicComponents().size() << " Dets " 
            << theRods.size() << " Rods "
            << " R: " << specificSurface().radius()
            << " Per.: " << bf.isPhiPeriodic()
            << " Overl.: " << isOverlapping;
}

MuRodBarrelLayer::~MuRodBarrelLayer ( )

override

Definition at line 70 of file MuRodBarrelLayer.cc.

References mps_fire::i, theBinFinder, and theRods.

                                    {
  delete theBinFinder;
  for (vector <const DetRod*>::iterator i = theRods.begin();
       i<theRods.end(); i++) {delete *i;}
}

Member Function Documentation

const std::vector<const GeomDet*>& MuRodBarrelLayer::basicComponents ( ) const

inlineoverride

Definition at line 27 of file MuRodBarrelLayer.h.

References compatibleDets(), components(), groupedCompatibleDets(), subDetector(), and theBasicComps.

Referenced by MuRodBarrelLayer().

{return theBasicComps;}

vector< GeometricSearchDet::DetWithState > MuRodBarrelLayer::compatibleDets ( const TrajectoryStateOnSurface &  startingState, const Propagator &  prop, const MeasurementEstimator &  est  ) const

override

Definition at line 78 of file MuRodBarrelLayer.cc.

References BaseBinFinder< T >::binIndex(), Surface::bounds(), GeometricSearchDet::compatibleDets(), funct::false, TrajectoryStateOnSurface::globalPosition(), createfilelist::int, isOverlapping, TrajectoryStateOnSurface::localPosition(), LogTrace, metname, GetRecoTauVFromDQM_MC_cff::next, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), GloballyPositioned< T >::position(), mps_fire::result, AlCaHLTBitMon_QueryRunRegistry::string, DetRod::surface(), theBinFinder, theRods, Bounds::width(), PV3DBase< T, PVType, FrameType >::x(), and xError().

Referenced by basicComponents().

                                                        {

  const std::string metname = "Muon|RecoMuon|RecoMuonDetLayers|MuRodBarrelLayer";
  vector<DetWithState> result; 

 
  LogTrace(metname) << "MuRodBarrelLayer::compatibleDets, Cyl R: " 
                    << specificSurface().radius()
                    << " TSOS at R= " << startingState.globalPosition().perp()
            << " phi= " << startingState.globalPosition().phi();

  pair<bool, TrajectoryStateOnSurface> compat =
    compatible(startingState, prop, est);
  if (!compat.first) {
    LogTrace(metname) << "     MuRodBarrelLayer::compatibleDets: not compatible"
                      << " (should not have been selected!)";
    return vector<DetWithState>();
  } 

  TrajectoryStateOnSurface& tsos = compat.second;

  LogTrace(metname) << "     MuRodBarrelLayer::compatibleDets, reached layer at: "
            << tsos.globalPosition()
            << " R = " << tsos.globalPosition().perp()
            << " phi = " << tsos.globalPosition().phi();

  int closest = theBinFinder->binIndex(tsos.globalPosition().phi());
  const DetRod* closestRod = theRods[closest];

  // Check the closest rod
  LogTrace(metname) << "     MuRodBarrelLayer::compatibleDets, closestRod: " << closest
                    << " phi : " << closestRod->surface().position().phi()
                    << " FTS phi: " << tsos.globalPosition().phi();
  
  result = closestRod->compatibleDets(tsos, prop, est);

  int nclosest = result.size(); // Debug counter

  bool checknext = false ;
  double dist;

  if (!result.empty()) { 
    // Check if the track go outside closest rod, then look for closest. 
    TrajectoryStateOnSurface& predictedState = result.front().second;
    float xErr = xError(predictedState, est);
    float halfWid = closestRod->surface().bounds().width()/2.;
    dist = predictedState.localPosition().x();

    // If the layer is overlapping, additionally reduce halfWid by 10%
    // to account for overlap.
    // FIXME: should we account for the real amount of overlap?
    if (isOverlapping) halfWid *= 0.9;

    if (fabs(dist) + xErr > halfWid) {
      checknext = true;
    }
  } else { // Rod is not compatible
    //FIXME: Usually next cannot be either. Implement proper logic.
    // (in general at least one rod should be when this method is called by
    // compatibleDets() which calls compatible())
    checknext = true;
    
    // Look for the next-to closest in phi.
    // Note Geom::Phi, subtraction is pi-border-safe
    if ( tsos.globalPosition().phi()-closestRod->surface().position().phi()>0.)
    {
      dist = -1.;
    } else {
      dist = +1.;
    }

    
    LogTrace(metname) << "     MuRodBarrelLayer::fastCompatibleDets, none on closest rod!";
  }

  if (checknext) {
    int next;
    if (dist<0.) next = closest+1;
    else next = closest-1;

    next = theBinFinder->binIndex(next); // Bin Periodicity
    const DetRod* nextRod = theRods[next];

  
    LogTrace(metname) << "     MuRodBarrelLayer::fastCompatibleDets, next-to closest"
                      << " rod: " << next << " dist " << dist
                      << " phi : " << nextRod->surface().position().phi()
                      << " FTS phi: " << tsos.globalPosition().phi();   
    
    vector<DetWithState> nextRodDets =
      nextRod->compatibleDets(tsos, prop, est);
    result.insert(result.end(), 
          nextRodDets.begin(), nextRodDets.end());
  }
  
  
   LogTrace(metname) << "     MuRodBarrelLayer::fastCompatibleDets: found: "
                     << result.size()
                     << " on closest: " << nclosest
                     << " # checked rods: " << 1 + int(checknext);
  
  return result;
}

const vector< const GeometricSearchDet * > & MuRodBarrelLayer::components ( ) const

override

Definition at line 201 of file MuRodBarrelLayer.cc.

References theComponents.

Referenced by basicComponents().

                                   {
  return theComponents;
}

vector< DetGroup > MuRodBarrelLayer::groupedCompatibleDets ( const TrajectoryStateOnSurface &  startingState, const Propagator &  prop, const MeasurementEstimator &  est  ) const

override

Definition at line 186 of file MuRodBarrelLayer.cc.

References gather_cfg::cout.

Referenced by basicComponents().

                                                            {
  // FIXME should return only 1 group 
  cout << "dummy implementation of MuRodBarrelLayer::groupedCompatibleDets()" << endl;
  return vector<DetGroup>();
}

virtual const std::vector<const DetRod*>& MuRodBarrelLayer::rods ( ) const

inlinevirtual

Return the vector of rods.

Definition at line 48 of file MuRodBarrelLayer.h.

References theRods, and xError().

{return theRods;}

GeomDetEnumerators::SubDetector MuRodBarrelLayer::subDetector ( ) const

override

Definition at line 196 of file MuRodBarrelLayer.cc.

References theBasicComps.

Referenced by basicComponents().

                                                                  {
  return theBasicComps.front()->subDetector();
}

float MuRodBarrelLayer::xError ( const TrajectoryStateOnSurface &  tsos, const MeasurementEstimator &  est  ) const

private

Definition at line 205 of file MuRodBarrelLayer.cc.

References TrajectoryStateOnSurface::hasError(), TrajectoryStateOnSurface::localError(), LocalTrajectoryError::positionError(), mathSSE::sqrt(), and LocalError::xx().

Referenced by compatibleDets(), and rods().

                                                          {
  const float nSigmas = 3.f;
  if (tsos.hasError()) {
    return nSigmas * sqrt(tsos.localError().positionError().xx());
  }
  else return nSigmas * 0.5;
}

Member Data Documentation

bool MuRodBarrelLayer::isOverlapping

private

Definition at line 60 of file MuRodBarrelLayer.h.

Referenced by compatibleDets(), and MuRodBarrelLayer().

std::vector<const GeomDet*> MuRodBarrelLayer::theBasicComps

private

Definition at line 58 of file MuRodBarrelLayer.h.

Referenced by basicComponents(), MuRodBarrelLayer(), and subDetector().

BaseBinFinder<double>* MuRodBarrelLayer::theBinFinder

private

Definition at line 59 of file MuRodBarrelLayer.h.

Referenced by compatibleDets(), MuRodBarrelLayer(), and ~MuRodBarrelLayer().

std::vector<const GeometricSearchDet*> MuRodBarrelLayer::theComponents

private

Definition at line 57 of file MuRodBarrelLayer.h.

Referenced by components(), and MuRodBarrelLayer().

std::vector<const DetRod*> MuRodBarrelLayer::theRods

private

Definition at line 56 of file MuRodBarrelLayer.h.

Referenced by compatibleDets(), MuRodBarrelLayer(), rods(), and ~MuRodBarrelLayer().