CMSSW Reference Manual

HICMeasurementEstimator::HICMeasurementEstimator	(	const GeometricSearchTracker *	theTracker0,
		const MagneticField *	mf,
		double	maxChi2,
		double	nSigma = `3.`
	)			`[inline, explicit]`

Construct with cuts on chi2 and nSigma.

The cut on Chi2 is used to define the acceptance of RecHits. The errors of the trajectory state are multiplied by nSigma to define acceptance of BoundPlane and maximalLocalDisplacement.

Definition at line 30 of file HICMeasurementEstimator.h.

References GeometricSearchTracker::barrelLayers(), bl, field, fneg, fpos, GeometricSearchTracker::negForwardLayers(), GeometricSearchTracker::posForwardLayers(), setHICDetMap(), and theTracker.

Referenced by clone().

00031                                                                :
00032   Chi2MeasurementEstimatorBase(maxChi2,nSigma)
00033     {
00034       theTracker = theTracker0;
00035       bl = theTracker->barrelLayers();
00036       fpos = theTracker->posForwardLayers();
00037       fneg = theTracker->negForwardLayers();
00038       field = mf;
00039       setHICDetMap();
00040     }

void HICMeasurementEstimator::chooseCuts ( int & i ) [virtual]

Definition at line 277 of file HICMeasurementEstimator.cc.

References theChi2Cut, theCutType, theNewCut, thePhiBound, thePhiBoundMean, thePhiCut, thePhiCutMean, thePhiWin, thePhiWinB, thePhiWinMean, theZBound, theZBoundMean, theZCut, theZCutMean, theZWin, theZWinB, and theZWinMean.

00278 {
00279 
00280       theChi2Cut = theNewCut;
00281 //      cout<<" Choose Chi2Cut "<<theChi2Cut<<endl;
00282       
00283    if( i == 1 )
00284    {
00285       thePhiBound = thePhiWin;
00286       theZBound = theZWin;
00287       thePhiBoundMean = thePhiWinMean;
00288       theZBoundMean = theZWinMean;
00289 //      cout<<" HICMeasurementEstimator::chooseCuts  "<<i<<" "<<thePhiBound<<" "<<theZBound<<endl;
00290    } 
00291       if( i == 2 )
00292    {
00293       thePhiBound = thePhiWinB;
00294       theZBound = theZWinB;
00295       thePhiBoundMean = thePhiWinMean;
00296       theZBoundMean = theZWinMean;
00297 //      cout<<" HICMeasurementEstimator::chooseCuts  "<<i<<" "<<thePhiBound<<" "<<theZBound<<endl;
00298    }  
00299       if( i == 3 )
00300    {
00301       thePhiBound = thePhiCut;
00302       theZBound = theZCut;
00303       thePhiBoundMean = thePhiCutMean;
00304       theZBoundMean = theZCutMean;
00305 //      cout<<" HICMeasurementEstimator::chooseCuts  "<<i<<" "<<thePhiBound<<" "<<theZBound<<endl;
00306    }
00307      
00308    theCutType = i;
00309 }

virtual HICMeasurementEstimator* HICMeasurementEstimator::clone ( void ) const [inline, virtual]

Implements MeasurementEstimator.

Definition at line 71 of file HICMeasurementEstimator.h.

References HICMeasurementEstimator().

00071                                                  {
00072     return new HICMeasurementEstimator(*this);
00073   }

bool HICMeasurementEstimator::estimate	(	const TrajectoryStateOnSurface &	ts,
		const BoundPlane &	plane
	)			const `[virtual]`

Returns true if the TrajectoryStateOnSurface is compatible with the BoundPlane, false otherwise.

The TrajectoryStateOnSurface must be on the plane.

Reimplemented from Chi2MeasurementEstimatorBase.

Definition at line 97 of file HICMeasurementEstimator.cc.

References BoundSurface::bounds(), GenMuonPlsPt100GeV_cfg::cout, lat::endl(), TrajectoryStateOnSurface::globalPosition(), Bounds::length(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), pi, GloballyPositioned< T >::position(), Bounds::width(), and PV3DBase< T, PVType, FrameType >::z().

00099 {
00100 
00101 //  cout<<" start estimate plane "<<endl;
00102   double pi = 4.*atan(1.);
00103   double twopi = 2.*pi;
00104   float theZError = plane.bounds().length() + 4.;
00105   float thePhiError = 2.*plane.bounds().width()/plane.position().perp(); 
00106 // Change 02.07.08
00107 //  float thePhiError = 4.*plane.bounds().width()/plane.position().perp();
00108 
00109 #ifdef DEBUG  
00110   cout<<" ======================================================================================== ";
00111   cout<<" Estimate detector::   tsos      :     detector   :   Error "<<endl;
00112   cout<<" R                 "<<ts.globalPosition().perp()<<" "<<plane.position().perp()<<" "<<theZError<<endl;
00113   cout<<" Phi               "<<ts.globalPosition().phi()<<" "<<plane.position().phi()<<" "<<thePhiError<<endl;
00114   cout<<" Z                 "<<ts.globalPosition().z()<<" "<<plane.position().z()<<" "<<theZError<<endl;
00115 #endif
00116 
00117   bool flag = false;
00118   if(fabs(ts.globalPosition().perp()-plane.position().perp())<theZError){
00119    if(fabs(ts.globalPosition().z()-plane.position().z())<theZError){
00120    float phi1 = ts.globalPosition().phi();
00121    float phi2 = plane.position().phi();
00122    if(phi1<0.) phi1 = twopi+phi1;
00123    if(phi2<0.) phi2 = twopi+phi2;
00124    float dfi = fabs(phi1-phi2);
00125    if(dfi>pi) dfi = twopi-dfi;
00126       if(dfi<thePhiError) flag = true;
00127    }
00128   }
00129 #ifdef DEBUG
00130     cout<<" Estimate = "<<flag<<endl;
00131 #endif
00132 
00133   return flag;
00134   
00135 }

template<unsigned int D>

std::pair<bool,double> HICMeasurementEstimator::estimate	(	const TrajectoryStateOnSurface &	ts,
		const TransientTrackingRecHit &	hit
	)			const `[inline, virtual]`

Returns pair( true, value) if the TrajectoryStateOnSurface is compatible with the RecHit, and pair( false, value) if it is not compatible.

The TrajectoryStateOnSurface must be on the same Surface as the RecHit. For an estimator where there is no value computed, e.g. fixed window estimator, only the first(bool) part is of interest.

Implements Chi2MeasurementEstimatorBase.

std::pair< bool, double > HICMeasurementEstimator::estimate	(	const TrajectoryStateOnSurface &	ts,
		const TransientTrackingRecHit &	hit
	)			const `[inline, virtual]`

Returns pair( true, value) if the TrajectoryStateOnSurface is compatible with the RecHit, and pair( false, value) if it is not compatible.

The TrajectoryStateOnSurface must be on the same Surface as the RecHit. For an estimator where there is no value computed, e.g. fixed window estimator, only the first(bool) part is of interest.

Implements Chi2MeasurementEstimatorBase.

Definition at line 11 of file HICMeasurementEstimator.cc.

References GenMuonPlsPt100GeV_cfg::cout, TrackingRecHit::dimension(), lat::endl(), Exception, and TrajectoryStateOnSurface::isValid().

00012                                                                                  {
00013     std::pair<bool,double> flag(false,0.);
00014     if(!tsos.isValid()) {
00015 #ifdef DEBUG
00016           std::cout<<" HICMeasurementEstimator::estimate::trajectory is not valid "<<std::endl;
00017 #endif
00018     return flag;
00019     }
00020 
00021     switch (aRecHit.dimension()) {
00022         case 1: return estimate<1>(tsos,aRecHit);
00023         case 2: return estimate<2>(tsos,aRecHit);
00024         case 3: return estimate<3>(tsos,aRecHit);
00025         case 4: return estimate<4>(tsos,aRecHit);
00026         case 5: return estimate<5>(tsos,aRecHit);
00027     }
00028     throw cms::Exception("RecHit of invalid size (not 1,2,3,4,5)");
00029 }

int HICMeasurementEstimator::getDetectorCode ( const DetLayer * a ) [virtual]

Definition at line 311 of file HICMeasurementEstimator.cc.

References GeomDetEnumerators::barrel, GeomDetEnumerators::endcap, DetLayer::location(), GloballyPositioned< T >::position(), edm::second(), GeometricSearchDet::surface(), theBackward, theBarrel, and theForward.

00312 {
00313      int layer = 0;
00314      if( a->location() == GeomDetEnumerators::barrel ) 
00315      {
00316            layer = (*theBarrel.find(a)).second;
00317      }     
00318      if( a->location() == GeomDetEnumerators::endcap ) 
00319      {
00320            if( a->surface().position().z() > 0. ) { layer = 100+(*theForward.find(a)).second;}
00321            if( a->surface().position().z() < 0. ) { layer = -100-(*theBackward.find(a)).second;}           
00322      }
00323      return layer;
00324 }

const MagneticField* HICMeasurementEstimator::getField ( ) [inline]

Definition at line 69 of file HICMeasurementEstimator.h.

References field.

00069 {return field;}

cms::HICConst* HICMeasurementEstimator::getHICConst ( ) [inline]

Definition at line 55 of file HICMeasurementEstimator.h.

References theHICConst.

00055 {return theHICConst;}

vector< double > HICMeasurementEstimator::setCuts	(	Trajectory &	traj,
		const DetLayer *	theCurrentLayer
	)			`[virtual]`

Definition at line 137 of file HICMeasurementEstimator.cc.

References a, GeomDetEnumerators::barrel, GenMuonPlsPt100GeV_cfg::cout, Trajectory::data(), GeomDetEnumerators::endcap, lat::endl(), first, prof2calltree::last, DetLayer::location(), Trajectory::measurements(), GeomDetEnumerators::PixelBarrel, GeomDetEnumerators::PixelEndcap, GloballyPositioned< T >::position(), edm::second(), DetLayer::subDetector(), GeometricSearchDet::surface(), theBackward, theBarrel, theForward, theLastLayer, theLayer, theLowMult, theNewCut, theNewCutB, thePhiCut, thePhiCutMean, thePhiWin, thePhiWinB, thePhiWinMean, thePhiWinMeanB, theTrajectorySize, theZCut, theZCutMean, theZWin, theZWinB, theZWinMean, and theZWinMeanB.

00138 {
00139      vector<double> theCuts;
00140      const DetLayer* a = traj.data().back().layer();
00141      const DetLayer* first = traj.data().front().layer();
00142      const DetLayer* last = traj.data().front().layer();
00143      thePhiWinMean = 0.;
00144      theZWinMean = 0.;
00145      thePhiWin = 0.;
00146      theZWin = 0.;
00147      theNewCut = 10.; // change 5->10 03.07.2008
00148      theNewCutB = 5.;
00149      
00150      thePhiWinMeanB = 0.002;
00151      theZWinMeanB = 0.;
00152      thePhiWinB = 0.008;
00153      theZWinB = 17.;
00154           
00155      theZCutMean = 0.;
00156      thePhiCutMean = 0.;
00157      thePhiCut = 0.;
00158      theZCut = 0.;
00159      
00160      theLayer = b;
00161      theLastLayer = a;
00162 
00163      
00164      theTrajectorySize = traj.data().size();
00165      
00166      
00167      if( theBarrel.size() == 0 || theForward.size() == 0 )
00168      {
00169 #ifdef DEBUG
00170        cout<<" HICMeasurementEstimator::setCuts:: no datector map "<<endl;
00171 #endif
00172         return theCuts;
00173      }
00174      
00175      if( a->location() == GeomDetEnumerators::barrel )
00176      {
00177        if( first->location() == GeomDetEnumerators::barrel ) 
00178        {
00179         thePhiWin = (*theHICConst).phiwinbar[(*theBarrel.find(first)).second][(*theBarrel.find(a)).second][(*theBarrel.find(b)).second];
00180         theZWin = (*theHICConst).zwinbar[(*theBarrel.find(first)).second][(*theBarrel.find(a)).second][(*theBarrel.find(b)).second];
00181         thePhiCut = (*theHICConst).phicutbar[(*theBarrel.find(first)).second][(*theBarrel.find(theLastLayer)).second][(*theBarrel.find(b)).second];
00182         theZCut = (*theHICConst).zcutbar[(*theBarrel.find(first)).second][(*theBarrel.find(theLastLayer)).second][(*theBarrel.find(b)).second];
00183 //        cout<<" Barrel first -Barrel cuts::layers "<<(*theBarrel.find(first)).second<<" "<<(*theBarrel.find(a)).second<<" "<<(*theBarrel.find(b)).second<<endl; 
00184 //        cout<<" Barrel first -Barrel cuts "<< thePhiWin<<" "<<theZWin <<" "<<thePhiCut <<" "<<theZCut<<endl;
00185        }
00186          else
00187          {
00188           if(first->surface().position().z() > 0. )
00189           {
00190           thePhiWin = (*theHICConst).phiwinfbb[(*theForward.find(first)).second][(*theBarrel.find(a)).second][(*theBarrel.find(b)).second];
00191           theZWin = (*theHICConst).zwinfbb[(*theForward.find(first)).second][(*theBarrel.find(a)).second][(*theBarrel.find(b)).second];
00192           thePhiCut = (*theHICConst).phicutfbb[(*theForward.find(first)).second][(*theBarrel.find(theLastLayer)).second][(*theBarrel.find(b)).second];
00193           theZCut = (*theHICConst).zcutfbb[(*theForward.find(first)).second][(*theBarrel.find(theLastLayer)).second][(*theBarrel.find(b)).second];
00194 //          cout<<" Endcap first positive -Barrel cuts::layers "<<(*theForward.find(first)).second<<" "<<(*theBarrel.find(a)).second<<" "<<(*theBarrel.find(b)).second<<endl;       
00195 //          cout<<" Endcap first positive -Barrel cuts "<< thePhiWin<<" "<<theZWin <<" "<<thePhiCut <<" "<<theZCut<<endl;
00196           } else {
00197           thePhiWin = (*theHICConst).phiwinfbb[(*theBackward.find(first)).second][(*theBarrel.find(a)).second][(*theBarrel.find(b)).second];
00198           theZWin = (*theHICConst).zwinfbb[(*theBackward.find(first)).second][(*theBarrel.find(a)).second][(*theBarrel.find(b)).second];
00199           thePhiCut = (*theHICConst).phicutfbb[(*theBackward.find(first)).second][(*theBarrel.find(theLastLayer)).second][(*theBarrel.find(b)).second];
00200           theZCut = (*theHICConst).zcutfbb[(*theBackward.find(first)).second][(*theBarrel.find(theLastLayer)).second][(*theBarrel.find(b)).second];
00201 //          cout<<" Endcap first negative -Barrel cuts::layers "<<(*theBackward.find(first)).second<<" "<<(*theBarrel.find(a)).second<<" "<<(*theBarrel.find(b)).second<<endl;
00202 //          cout<<" Endcap first negative -Barrel cuts "<< thePhiWin<<" "<<theZWin <<" "<<thePhiCut <<" "<<theZCut<<endl;
00203          }
00204          } //
00205 
00206 
00207         theCuts.push_back(thePhiWin); theCuts.push_back(theZWin);
00208         theCuts.push_back(thePhiCut); theCuts.push_back(theZCut);
00209 
00210         return theCuts;
00211      }
00212      if( a->location() == GeomDetEnumerators::endcap && b->location() == GeomDetEnumerators::endcap)
00213      {
00214         if( a->surface().position().z() > 0. )
00215         {
00216         thePhiWin = (*theHICConst).phiwinfrw[(*theForward.find(first)).second][(*theForward.find(a)).second][(*theForward.find(b)).second];
00217         theZWin = (*theHICConst).zwinfrw[(*theForward.find(first)).second][(*theForward.find(a)).second][(*theForward.find(b)).second];
00218         thePhiCut = (*theHICConst).phicutfrw[(*theForward.find(first)).second][(*theForward.find(theLastLayer)).second][(*theForward.find(b)).second];
00219         theZCut = (*theHICConst).zcutfrw[(*theForward.find(first)).second][(*theForward.find(theLastLayer)).second][(*theForward.find(b)).second];
00220         }
00221           else
00222           {
00223            thePhiWin = (*theHICConst).phiwinfrw[(*theBackward.find(first)).second][(*theBackward.find(a)).second][(*theBackward.find(b)).second];
00224            theZWin = (*theHICConst).zwinfrw[(*theBackward.find(first)).second][(*theBackward.find(a)).second][(*theBackward.find(b)).second];
00225            thePhiCut = (*theHICConst).phicutfrw[(*theBackward.find(first)).second][(*theBackward.find(theLastLayer)).second][(*theBackward.find(b)).second];
00226            theZCut = (*theHICConst).zcutfrw[(*theBackward.find(first)).second][(*theBackward.find(theLastLayer)).second][(*theBackward.find(b)).second];
00227           }
00228         
00229         if( theLowMult == 1 )
00230         {
00231         if( b->subDetector() == GeomDetEnumerators::PixelEndcap ) theNewCut = 20.;
00232         if( traj.measurements().size() == 1 ) theNewCut = 20.; 
00233         theNewCutB = 30.;
00234         }
00235         
00236         thePhiWinMeanB = 0.004;
00237         thePhiWinB = 0.04;
00238 
00239 
00240         theCuts.push_back(thePhiWin); theCuts.push_back(theZWin);
00241         theCuts.push_back(thePhiCut); theCuts.push_back(theZCut);
00242         
00243         return theCuts;
00244      }
00245      if( a->location() == GeomDetEnumerators::endcap && b->location() == GeomDetEnumerators::barrel )
00246      {
00247        
00248        if( a->surface().position().z() > 0. )
00249        {
00250         thePhiWin = (*theHICConst).phiwinbfrw[(*theForward.find(first)).second][(*theForward.find(a)).second][(*theBarrel.find(b)).second];
00251         theZWin = (*theHICConst).zwinbfrw[(*theForward.find(first)).second][(*theForward.find(a)).second][(*theBarrel.find(b)).second];
00252         thePhiCut = (*theHICConst).phicutbfrw[(*theForward.find(first)).second][(*theForward.find(theLastLayer)).second][(*theBarrel.find(b)).second];
00253         theZCut = (*theHICConst).zcutbfrw[(*theForward.find(first)).second][(*theForward.find(theLastLayer)).second][(*theBarrel.find(b)).second];
00254        }
00255           else
00256           {
00257             thePhiWin = (*theHICConst).phiwinbfrw[(*theBackward.find(first)).second][(*theBackward.find(a)).second][(*theBarrel.find(b)).second];
00258             theZWin = (*theHICConst).zwinbfrw[(*theBackward.find(first)).second][(*theBackward.find(a)).second][(*theBarrel.find(b)).second];
00259             thePhiCut = (*theHICConst).phicutbfrw[(*theBackward.find(first)).second][(*theBackward.find(theLastLayer)).second][(*theBarrel.find(b)).second];
00260             theZCut = (*theHICConst).zcutbfrw[(*theBackward.find(first)).second][(*theBackward.find(theLastLayer)).second][(*theBarrel.find(b)).second];
00261           }     
00262         
00263         if( b->subDetector() ==  GeomDetEnumerators::PixelBarrel) theNewCut = 20.;
00264 
00265         thePhiWinMeanB = 0.004;
00266         thePhiWinB = 0.016;
00267 
00268         theCuts.push_back(thePhiWin); theCuts.push_back(theZWin);
00269         theCuts.push_back(thePhiCut); theCuts.push_back(theZCut);
00270         
00271         return theCuts;
00272      }
00273 //     cout<<" HICMeasurementEstimator::setCuts::Error: unknown detector layer "<<endl;
00274      return theCuts;
00275 }

void HICMeasurementEstimator::setHICConst ( cms::HICConst * hh ) [inline]

Definition at line 57 of file HICMeasurementEstimator.h.

References theHICConst.

Referenced by cms::HITrackVertexMaker::produceTracks().

00057 {theHICConst = hh;}

void HICMeasurementEstimator::setHICDetMap ( ) [virtual]

Definition at line 326 of file HICMeasurementEstimator.cc.

References bl, GenMuonPlsPt100GeV_cfg::cout, lat::endl(), fneg, fpos, theBackward, theBarrel, and theForward.

Referenced by HICMeasurementEstimator().

00327 {
00328 
00329 #ifdef DEBUG
00330   std::cout<<" Set Detector Map... "<<std::endl;
00331 #endif
00332   
00333   int ila=0;
00334   for ( std::vector<BarrelDetLayer*>::const_iterator ilayer = bl.begin(); ilayer != bl.end(); ilayer++)
00335   {
00336      theBarrel[(*ilayer)]=ila;
00337      ila++;
00338   }
00339 //
00340 // The same for forward part.
00341 //
00342 
00343    int ilf1 = 0; 
00344    int ilf2 = 0;
00345    for ( vector<ForwardDetLayer*>::const_iterator ilayer = fpos.begin(); 
00346                                                     ilayer != fpos.end(); ilayer++)
00347    {
00348                 theForward[(*ilayer)] = ilf1;
00349                 ilf1++;
00350    }
00351    for ( vector<ForwardDetLayer*>::const_iterator ilayer = fneg.begin(); 
00352                                                     ilayer != fneg.end(); ilayer++)
00353    {
00354 #ifdef DEBUG
00355             cout<<" HICDetectorMap::negative layers "<<(**ilayer).position().z()<<" "<<ilf2<<endl;
00356 #endif
00357                    theBackward[(*ilayer)] = ilf2;
00358                    ilf2++;
00359    }
00360 
00361 }

virtual void HICMeasurementEstimator::setLastLayer ( DetLayer *& ll ) [inline, virtual]

Definition at line 59 of file HICMeasurementEstimator.h.

References theLastLayer.

00059 {theLastLayer = ll;};

virtual void HICMeasurementEstimator::setMult ( int aMult = 1 ) [inline, virtual]

Definition at line 65 of file HICMeasurementEstimator.h.

References theLowMult.

Referenced by cms::HITrackVertexMaker::produceTracks().

00065 {theLowMult=aMult;}

virtual void HICMeasurementEstimator::setSign ( int & i ) [inline, virtual]

Definition at line 67 of file HICMeasurementEstimator.h.

References theSign.

00067 {theSign=i;}

std::vector<BarrelDetLayer*> HICMeasurementEstimator::bl [private]

Definition at line 112 of file HICMeasurementEstimator.h.

Referenced by HICMeasurementEstimator(), and setHICDetMap().

const MagneticField* HICMeasurementEstimator::field [private]

Definition at line 116 of file HICMeasurementEstimator.h.

Referenced by getField(), and HICMeasurementEstimator().

std::vector<ForwardDetLayer*> HICMeasurementEstimator::fneg [private]

Definition at line 114 of file HICMeasurementEstimator.h.

Referenced by HICMeasurementEstimator(), and setHICDetMap().

std::vector<ForwardDetLayer*> HICMeasurementEstimator::fpos [private]

Definition at line 113 of file HICMeasurementEstimator.h.

Referenced by HICMeasurementEstimator(), and setHICDetMap().

std::map<const DetLayer*,int> HICMeasurementEstimator::theBackward [private]

Definition at line 81 of file HICMeasurementEstimator.h.

Referenced by getDetectorCode(), setCuts(), and setHICDetMap().

std::map<const DetLayer*,int> HICMeasurementEstimator::theBarrel [private]

Definition at line 79 of file HICMeasurementEstimator.h.

Referenced by getDetectorCode(), setCuts(), and setHICDetMap().

double HICMeasurementEstimator::theChi2Cut [private]

Definition at line 105 of file HICMeasurementEstimator.h.

Referenced by chooseCuts().

int HICMeasurementEstimator::theCutType [private]

Definition at line 110 of file HICMeasurementEstimator.h.

Referenced by chooseCuts().

const DetLayer* HICMeasurementEstimator::theFirstLayer [private]

Definition at line 89 of file HICMeasurementEstimator.h.

std::map<const DetLayer*,int> HICMeasurementEstimator::theForward [private]

Definition at line 80 of file HICMeasurementEstimator.h.

Referenced by getDetectorCode(), setCuts(), and setHICDetMap().

cms::HICConst* HICMeasurementEstimator::theHICConst [private]

Definition at line 76 of file HICMeasurementEstimator.h.

Referenced by getHICConst(), and setHICConst().

const DetLayer* HICMeasurementEstimator::theLastLayer [private]

Definition at line 87 of file HICMeasurementEstimator.h.

Referenced by setCuts(), and setLastLayer().

const DetLayer* HICMeasurementEstimator::theLayer [private]

Definition at line 88 of file HICMeasurementEstimator.h.

Referenced by setCuts().

int HICMeasurementEstimator::theLowMult [private]

Definition at line 109 of file HICMeasurementEstimator.h.

Referenced by setCuts(), and setMult().

double HICMeasurementEstimator::theMaxChi2 [private]

Reimplemented from Chi2MeasurementEstimatorBase.

Definition at line 77 of file HICMeasurementEstimator.h.

double HICMeasurementEstimator::theNewCut [private]

Definition at line 106 of file HICMeasurementEstimator.h.

Referenced by chooseCuts(), and setCuts().

double HICMeasurementEstimator::theNewCutB [private]

Definition at line 107 of file HICMeasurementEstimator.h.

Referenced by setCuts().

int HICMeasurementEstimator::theNSigma [private]

Reimplemented from Chi2MeasurementEstimatorBase.

Definition at line 78 of file HICMeasurementEstimator.h.

double HICMeasurementEstimator::thePhiBound [private]

Definition at line 83 of file HICMeasurementEstimator.h.

Referenced by chooseCuts().

double HICMeasurementEstimator::thePhiBoundMean [private]

Definition at line 85 of file HICMeasurementEstimator.h.

Referenced by chooseCuts().

double HICMeasurementEstimator::thePhiCut [private]

Definition at line 101 of file HICMeasurementEstimator.h.

Referenced by chooseCuts(), and setCuts().

double HICMeasurementEstimator::thePhiCutMean [private]

Definition at line 103 of file HICMeasurementEstimator.h.

Referenced by chooseCuts(), and setCuts().

double HICMeasurementEstimator::thePhiWin [private]

Definition at line 91 of file HICMeasurementEstimator.h.

Referenced by chooseCuts(), and setCuts().

double HICMeasurementEstimator::thePhiWinB [private]

Definition at line 96 of file HICMeasurementEstimator.h.

Referenced by chooseCuts(), and setCuts().

double HICMeasurementEstimator::thePhiWinMean [private]

Definition at line 93 of file HICMeasurementEstimator.h.

Referenced by chooseCuts(), and setCuts().

double HICMeasurementEstimator::thePhiWinMeanB [private]

Definition at line 98 of file HICMeasurementEstimator.h.

Referenced by setCuts().

int HICMeasurementEstimator::theSign [private]

Definition at line 115 of file HICMeasurementEstimator.h.

Referenced by setSign().

const GeometricSearchTracker* HICMeasurementEstimator::theTracker [private]

Definition at line 111 of file HICMeasurementEstimator.h.

Referenced by HICMeasurementEstimator().

int HICMeasurementEstimator::theTrajectorySize [private]

Definition at line 90 of file HICMeasurementEstimator.h.

Referenced by setCuts().

double HICMeasurementEstimator::theZBound [private]

Definition at line 84 of file HICMeasurementEstimator.h.

Referenced by chooseCuts().

double HICMeasurementEstimator::theZBoundMean [private]

Definition at line 86 of file HICMeasurementEstimator.h.

Referenced by chooseCuts().

double HICMeasurementEstimator::theZCut [private]

Definition at line 102 of file HICMeasurementEstimator.h.

Referenced by chooseCuts(), and setCuts().

double HICMeasurementEstimator::theZCutMean [private]

Definition at line 104 of file HICMeasurementEstimator.h.

Referenced by chooseCuts(), and setCuts().

double HICMeasurementEstimator::theZWin [private]

Definition at line 92 of file HICMeasurementEstimator.h.

Referenced by chooseCuts(), and setCuts().

double HICMeasurementEstimator::theZWinB [private]

Definition at line 97 of file HICMeasurementEstimator.h.

Referenced by chooseCuts(), and setCuts().

double HICMeasurementEstimator::theZWinMean [private]

Definition at line 94 of file HICMeasurementEstimator.h.

Referenced by chooseCuts(), and setCuts().

double HICMeasurementEstimator::theZWinMeanB [private]

Definition at line 99 of file HICMeasurementEstimator.h.

Referenced by setCuts().

HICMeasurementEstimator Class Reference

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation


Public Member Functions
virtual void	chooseCuts (int &i)
virtual HICMeasurementEstimator *	clone () const
bool	estimate (const TrajectoryStateOnSurface &ts, const BoundPlane &plane) const
	Returns true if the TrajectoryStateOnSurface is compatible with the BoundPlane, false otherwise.
template<unsigned int D>
std::pair< bool, double >	estimate (const TrajectoryStateOnSurface &, const TransientTrackingRecHit &) const
	Returns pair( true, value) if the TrajectoryStateOnSurface is compatible with the RecHit, and pair( false, value) if it is not compatible.
virtual std::pair< bool, double >	estimate (const TrajectoryStateOnSurface &, const TransientTrackingRecHit &) const
	Returns pair( true, value) if the TrajectoryStateOnSurface is compatible with the RecHit, and pair( false, value) if it is not compatible.
virtual int	getDetectorCode (const DetLayer *a)
const MagneticField *	getField ()
cms::HICConst *	getHICConst ()
	HICMeasurementEstimator (const GeometricSearchTracker theTracker0, const MagneticField mf, double maxChi2, double nSigma=3.)
	Construct with cuts on chi2 and nSigma.
virtual std::vector< double >	setCuts (Trajectory &traj, const DetLayer *theCurrentLayer)
void	setHICConst (cms::HICConst *hh)
virtual void	setHICDetMap ()
virtual void	setLastLayer (DetLayer *&ll)
virtual void	setMult (int aMult=1)
virtual void	setSign (int &i)
Private Attributes
std::vector< BarrelDetLayer * >	bl
const MagneticField *	field
std::vector< ForwardDetLayer * >	fneg
std::vector< ForwardDetLayer * >	fpos
std::map< const DetLayer *, int >	theBackward
std::map< const DetLayer *, int >	theBarrel
double	theChi2Cut
int	theCutType
const DetLayer *	theFirstLayer
std::map< const DetLayer *, int >	theForward
cms::HICConst *	theHICConst
const DetLayer *	theLastLayer
const DetLayer *	theLayer
int	theLowMult
double	theMaxChi2
double	theNewCut
double	theNewCutB
int	theNSigma
double	thePhiBound
double	thePhiBoundMean
double	thePhiCut
double	thePhiCutMean
double	thePhiWin
double	thePhiWinB
double	thePhiWinMean
double	thePhiWinMeanB
int	theSign
const GeometricSearchTracker *	theTracker
int	theTrajectorySize
double	theZBound
double	theZBoundMean
double	theZCut
double	theZCutMean
double	theZWin
double	theZWinB
double	theZWinMean
double	theZWinMeanB