L1TMuonBarrelKalmanAlgo Class Reference

#include <L1TMuonBarrelKalmanAlgo.h>

Classes
class	StubSorter
class	TrackSorter

Public Types
typedef ROOT::Math::SMatrix< double, 2, 2, ROOT::Math::MatRepSym< double, 2 > >	CovarianceMatrix2
typedef ROOT::Math::SMatrix< double, 1, 3 >	Matrix13
typedef ROOT::Math::SMatrix< double, 2, 3 >	Matrix23
typedef ROOT::Math::SMatrix< double, 3, 1 >	Matrix31
typedef ROOT::Math::SMatrix< double, 3, 2 >	Matrix32
typedef ROOT::Math::SMatrix< double, 3, 3 >	Matrix33
typedef ROOT::Math::SVector< double, 2 >	Vector2

Public Member Functions
void	addBMTFMuon (int, const L1MuKBMTrack &, std::unique_ptr< l1t::RegionalMuonCandBxCollection > &)
std::pair< bool, L1MuKBMTrack >	chain (const L1MuKBMTCombinedStubRef &, const L1MuKBMTCombinedStubRefVector &)
L1MuKBMTrackCollection	cleanAndSort (const L1MuKBMTrackCollection &, uint)
	L1TMuonBarrelKalmanAlgo (const edm::ParameterSet &settings)
void	resolveEtaUnit (L1MuKBMTrackCollection &)

Private Member Functions
int	correctedPhi (const L1MuKBMTCombinedStubRef &, int)
int	correctedPhiB (const L1MuKBMTCombinedStubRef &)
int	customBitmask (unsigned int, unsigned int, unsigned int, unsigned int)
int	encode (bool ownwheel, int sector, bool tag)
void	estimateChiSquare (L1MuKBMTrack &)
bool	getBit (int, int)
std::pair< bool, uint >	getByCode (const L1MuKBMTrackCollection &tracks, int mask)
int	hitPattern (const L1MuKBMTrack &)
std::pair< bool, uint >	match (const L1MuKBMTrack &, const L1MuKBMTCombinedStubRefVector &)
void	propagate (L1MuKBMTrack &)
int	rank (const L1MuKBMTrack &)
void	setFloatingPointValues (L1MuKBMTrack &, bool)
std::map< int, int >	trackAddress (const L1MuKBMTrack &, int &)
uint	twosCompToBits (int)
bool	update (L1MuKBMTrack &, const L1MuKBMTCombinedStubRef &, int)
void	updateEta (L1MuKBMTrack &, const L1MuKBMTCombinedStubRef &)
bool	updateLUT (L1MuKBMTrack &, const L1MuKBMTCombinedStubRef &, int)
bool	updateOffline (L1MuKBMTrack &, const L1MuKBMTCombinedStubRef &)
bool	updateOffline1D (L1MuKBMTrack &, const L1MuKBMTCombinedStubRef &)
void	vertexConstraint (L1MuKBMTrack &)
void	vertexConstraintLUT (L1MuKBMTrack &)
void	vertexConstraintOffline (L1MuKBMTrack &)
int	wrapAround (int, int)

Private Attributes
std::vector< double >	aPhi_
std::vector< double >	aPhiB_
std::vector< double >	aPhiBNLO_
std::vector< double >	bPhi_
std::vector< double >	bPhiB_
std::vector< double >	chiSquare_
std::vector< int >	chiSquareCut_
std::vector< int >	chiSquareCutCurv_
std::vector< int >	chiSquareCutPattern_
std::vector< int >	combos1_
std::vector< int >	combos2_
std::vector< int >	combos3_
std::vector< int >	combos4_
std::vector< double >	eLoss_
std::vector< double >	etaLUT0_
std::vector< double >	etaLUT1_
uint	globalChi2Cut_
std::vector< double >	initK2_
std::vector< double >	initK_
std::unique_ptr< L1TMuonBarrelKalmanLUTs >	lutService_
std::vector< double >	mScatteringPhi_
std::vector< double >	mScatteringPhiB_
double	pointResolutionPhi_
double	pointResolutionPhiB_
double	pointResolutionVertex_
bool	useOfflineAlgo_
bool	verbose_

Detailed Description

Definition at line 16 of file L1TMuonBarrelKalmanAlgo.h.

Member Typedef Documentation

typedef ROOT::Math::SMatrix<double,2,2,ROOT::Math::MatRepSym<double,2> > L1TMuonBarrelKalmanAlgo::CovarianceMatrix2

Definition at line 19 of file L1TMuonBarrelKalmanAlgo.h.

typedef ROOT::Math::SMatrix<double,1,3> L1TMuonBarrelKalmanAlgo::Matrix13

Definition at line 22 of file L1TMuonBarrelKalmanAlgo.h.

typedef ROOT::Math::SMatrix<double,2,3> L1TMuonBarrelKalmanAlgo::Matrix23

Definition at line 21 of file L1TMuonBarrelKalmanAlgo.h.

typedef ROOT::Math::SMatrix<double,3,1> L1TMuonBarrelKalmanAlgo::Matrix31

Definition at line 23 of file L1TMuonBarrelKalmanAlgo.h.

typedef ROOT::Math::SMatrix<double,3,2> L1TMuonBarrelKalmanAlgo::Matrix32

Definition at line 20 of file L1TMuonBarrelKalmanAlgo.h.

typedef ROOT::Math::SMatrix<double,3,3> L1TMuonBarrelKalmanAlgo::Matrix33

Definition at line 24 of file L1TMuonBarrelKalmanAlgo.h.

typedef ROOT::Math::SVector<double,2> L1TMuonBarrelKalmanAlgo::Vector2

Definition at line 18 of file L1TMuonBarrelKalmanAlgo.h.

Constructor & Destructor Documentation

L1TMuonBarrelKalmanAlgo::L1TMuonBarrelKalmanAlgo

(

const edm::ParameterSet &

settings

)

Definition at line 5 of file L1TMuonBarrelKalmanAlgo.cc.

                                                                               :
  verbose_(settings.getParameter<bool>("verbose")),
  lutService_(new L1TMuonBarrelKalmanLUTs(settings.getParameter<std::string>("lutFile"))),
  initK_(settings.getParameter<std::vector<double> >("initialK")),
  initK2_(settings.getParameter<std::vector<double> >("initialK2")),
  eLoss_(settings.getParameter<std::vector<double> >("eLoss")),
  aPhi_(settings.getParameter<std::vector<double> >("aPhi")),
  aPhiB_(settings.getParameter<std::vector<double> >("aPhiB")),
  aPhiBNLO_(settings.getParameter<std::vector<double> >("aPhiBNLO")),
  bPhi_(settings.getParameter<std::vector<double> >("bPhi")),
  bPhiB_(settings.getParameter<std::vector<double> >("bPhiB")),
  globalChi2Cut_(settings.getParameter<unsigned int>("globalChi2Cut")),
  chiSquare_(settings.getParameter<std::vector<double> >("chiSquare")),
  chiSquareCutPattern_(settings.getParameter<std::vector<int> >("chiSquareCutPattern")),
  chiSquareCutCurv_(settings.getParameter<std::vector<int> >("chiSquareCutCurvMax")),
  chiSquareCut_(settings.getParameter<std::vector<int> >("chiSquareCut")),
  combos4_(settings.getParameter<std::vector<int> >("combos4")),
  combos3_(settings.getParameter<std::vector<int> >("combos3")),
  combos2_(settings.getParameter<std::vector<int> >("combos2")),
  combos1_(settings.getParameter<std::vector<int> >("combos1")),
  useOfflineAlgo_(settings.getParameter<bool>("useOfflineAlgo")),
  mScatteringPhi_(settings.getParameter<std::vector<double> >("mScatteringPhi")),
  mScatteringPhiB_(settings.getParameter<std::vector<double> >("mScatteringPhiB")),
  pointResolutionPhi_(settings.getParameter<double>("pointResolutionPhi")),
  pointResolutionPhiB_(settings.getParameter<double>("pointResolutionPhiB")),
  pointResolutionVertex_(settings.getParameter<double>("pointResolutionVertex"))  

{



}

Member Function Documentation

void L1TMuonBarrelKalmanAlgo::addBMTFMuon	(	int	bx,
		const L1MuKBMTrack &	track,
		std::unique_ptr< l1t::RegionalMuonCandBxCollection > &	out
	)

Definition at line 48 of file L1TMuonBarrelKalmanAlgo.cc.

References funct::abs(), l1t::bmtf, L1MuKBMTrack::coarseEta(), L1MuKBMTrack::curvatureAtMuon(), L1MuKBMTrack::curvatureAtVertex(), L1MuKBMTrack::dxy(), PVValHelper::dxy, PVValHelper::eta, L1MuKBMTrack::fineEta(), objects.autophobj::float, L1MuKBMTrack::hasFineEta(), CaloTPGTranscoder_cfi::HF, createfilelist::int, M_PI, metsig::muon, phi, L1MuKBMTrack::phiAtMuon(), EnergyCorrector::pt, hiDetachedQuadStep_cff::pt2, jets_cff::quality, rank(), L1MuKBMTrack::sector(), l1t::RegionalMuonCand::setDataword(), l1t::RegionalMuonCand::setHwDXY(), l1t::RegionalMuonCand::setHwHF(), l1t::RegionalMuonCand::setHwPt2(), Validation_hcalonly_cfi::sign, trackAddress(), twosCompToBits(), and L1MuKBMTrack::wheel().

                                                                                                                                {

  int  K = abs(track.curvatureAtVertex());
  //calibration
  int sign,signValid;

  if (track.curvatureAtVertex()==0) {
    sign=0;
    signValid=0;
  }
  else if (track.curvatureAtVertex()>0) {
    sign=0;
    signValid=1;
  }
  else  {
    sign=1;
    signValid=1;
  }   

  if (K<22)
    K=22;

  if (K>2047)
    K=2047;

  float lsb=1.25/float(1<<13);
  int pt = int(2*(1.0/(lsb*float(K))));
  if (pt>511)
    pt=511;

  int  K2 = abs(track.curvatureAtMuon());
  if (K2<22)
    K2=22;

  if (K2>2047)
    K2=2047;
  int pt2 = int(1.0/(lsb*float(K2))); 
  if (pt2>254)
    pt2=254;
  int eta  = track.hasFineEta() ? track.fineEta() : track.coarseEta();
  
  int phi=track.phiAtMuon();
  phi=int((phi*M_PI/(6.0*2048.0)+15.0*M_PI/180.)/(0.625*M_PI/180.0));

  int processor=track.sector();
  int HF = track.hasFineEta();
  
  int quality=rank(track)/64;

  int dxy=abs(track.dxy())>>7;

  int trackAddr;
  std::map<int,int> addr = trackAddress(track,trackAddr);

  l1t::RegionalMuonCand muon(pt,phi,eta,sign,signValid,quality,processor,l1t::bmtf,addr);
  muon.setHwHF(HF);
  muon.setHwPt2(pt2);
  muon.setHwDXY(dxy);


  //nw the words!
  uint32_t word1=pt;
  word1=word1 | quality<<9;
  word1=word1 | (twosCompToBits(eta))<<13;
  word1=word1 | HF<<22;
  word1=word1 | (twosCompToBits(phi))<<23;

  uint32_t word2=sign;
  word2=word2 | signValid<<1;
  word2=word2 | dxy<<2;
  word2=word2 | trackAddr<<4;
  word2=word2 | (twosCompToBits(track.wheel()))<<20;
  word2=word2 | pt2<<23;
  muon.setDataword(word2,word1);
  out->push_back(bx,muon);
}

std::pair< bool, L1MuKBMTrack > L1TMuonBarrelKalmanAlgo::chain	(	const L1MuKBMTCombinedStubRef &	seed,
		const L1MuKBMTCombinedStubRefVector &	stubs
	)

Definition at line 619 of file L1TMuonBarrelKalmanAlgo.cc.

References funct::abs(), ALCARECOTkAlJpsiMuMu_cff::charge, cleanAndSort(), combos1_, combos2_, combos3_, combos4_, correctedPhi(), correctedPhiB(), estimateChiSquare(), objects.autophobj::float, getBit(), globalChi2Cut_, hitPattern(), mps_fire::i, initK2_, initK_, createfilelist::int, RecoTauDiscriminantConfiguration::mask, match(), pointResolutionPhi_, pointResolutionPhiB_, propagate(), rank(), resolveEtaUnit(), setFloatingPointValues(), HiIsolationCommonParameters_cff::track, parallelization::uint(), update(), verbose_, and vertexConstraint().

Referenced by L1TMuonBarrelKalmanRegionModule::process().

                                                                                                                                         {
  L1MuKBMTrackCollection pretracks; 
  std::vector<int> combinatorics;
  switch(seed->stNum()) {
  case 1:
    combinatorics=combos1_;
    break;
  case 2:
    combinatorics=combos2_;
    break;

  case 3:
    combinatorics=combos3_;
    break;

  case 4:
    combinatorics=combos4_;
    break;

  default:
    printf("Something really bad happend\n");
  }

  L1MuKBMTrack nullTrack(seed,correctedPhi(seed,seed->scNum()),correctedPhiB(seed));

  for( const auto& mask : combinatorics) {
    L1MuKBMTrack track(seed,correctedPhi(seed,seed->scNum()),correctedPhiB(seed));
    int phiB = correctedPhiB(seed);
    if (seed->quality()<4)
      phiB=0;
    int charge;
    if (phiB==0)
      charge = 0;
    else
      charge=phiB/abs(phiB);
    int initialK = int(initK_[seed->stNum()-1]*phiB/(1+initK2_[seed->stNum()-1]*charge*phiB));
    if (initialK>8192)
      initialK=8192;
    if (initialK<-8192)
      initialK=-8192;
    
    

    track.setCoordinates(seed->stNum(),initialK,correctedPhi(seed,seed->scNum()),phiB);
    track.setHitPattern(hitPattern(track));
    //Set eta coarse
    track.setCoarseEta(seed->coarseEta());
    //Set fine eta:
    if (seed->qeta1()>=0) {
      if (seed->qeta2()>=0) {
    track.setFineEta((seed->eta1()+seed->eta2())/2);
      }
      else {
    track.setFineEta(seed->eta1());
      }
    }

    //set covariance
    L1MuKBMTrack::CovarianceMatrix covariance;  


    float DK=512*512.;
    covariance(0,0)=DK;
    covariance(0,1)=0;
    covariance(0,2)=0;
    covariance(1,0)=0;
    covariance(1,1)=float(pointResolutionPhi_);
    covariance(1,2)=0;
    covariance(2,0)=0;
    covariance(2,1)=0;
    covariance(2,2)=float(pointResolutionPhiB_);
    track.setCovariance(covariance);
    //
    if (verbose_) {
      printf("New Kalman fit staring at step=%d, phi=%d,phiB=%d with curvature=%d\n",track.step(),track.positionAngle(),track.bendingAngle(),track.curvature());
      printf("BITMASK:");
      for (unsigned int i=0;i<4;++i)
    printf("%d",getBit(mask,i));
      printf("\n");
      printf("------------------------------------------------------\n");
      printf("------------------------------------------------------\n");
      printf("------------------------------------------------------\n");
      printf("stubs:\n");
      for (const auto& stub: stubs) 
    printf("station=%d phi=%d phiB=%d qual=%d \n",stub->stNum(),correctedPhi(stub,seed->scNum()),correctedPhiB(stub),stub->quality()); 
      printf("------------------------------------------------------\n");
      printf("------------------------------------------------------\n");

    }

    while(track.step()>0) {
      // muon station 1 
      if (track.step()==1) {
    track.setCoordinatesAtMuon(track.curvature(),track.positionAngle(),track.bendingAngle());
    setFloatingPointValues(track,false);
    if (verbose_) 
      printf ("Unconstrained PT  in Muon System: pt=%f\n",track.ptUnconstrained());
      }
      
      propagate(track);
      if (verbose_)
    printf("propagated Coordinates step:%d,phi=%d,phiB=%d,K=%d\n",track.step(),track.positionAngle(),track.bendingAngle(),track.curvature());

      if (track.step()>0) 
    if (getBit(mask,track.step()-1)) {
      std::pair<bool,uint> bestStub = match(track,stubs);
          if ((!bestStub.first) || (!update(track,stubs[bestStub.second],mask)))
        break;
      if (verbose_) {
        printf("updated Coordinates step:%d,phi=%d,phiB=%d,K=%d\n",track.step(),track.positionAngle(),track.bendingAngle(),track.curvature());
      }
    }
      if (track.step()==0) {
    track.setCoordinatesAtVertex(track.curvature(),track.positionAngle(),track.bendingAngle());
    if (verbose_)
      printf(" Coordinates before vertex constraint step:%d,phi=%d,dxy=%d,K=%d\n",track.step(),track.phiAtVertex(),track.dxy(),track.curvatureAtVertex());
    estimateChiSquare(track);
    if (uint(track.approxChi2())>globalChi2Cut_)
      break;
    vertexConstraint(track);
    if (verbose_) {
      printf(" Coordinates after vertex constraint step:%d,phi=%d,dxy=%d,K=%d  maximum local chi2=%d\n",track.step(),track.phiAtVertex(),track.dxy(),track.curvatureAtVertex(),track.approxChi2());
      printf("------------------------------------------------------\n");
      printf("------------------------------------------------------\n");
    }
    setFloatingPointValues(track,true);
    track.setRank(rank(track));
    if (verbose_)
      printf ("Floating point coordinates at vertex: pt=%f, eta=%f phi=%f\n",track.pt(),track.eta(),track.phi());
    pretracks.push_back(track);
      }
    }
  }

  //Resolve eta
  resolveEtaUnit(pretracks);
  //Now for all the pretracks we need only one 
  L1MuKBMTrackCollection cleaned = cleanAndSort(pretracks,1);

  if (!cleaned.empty())
    return std::make_pair(true,cleaned[0]);
  return std::make_pair(false,nullTrack);
}       

L1MuKBMTrackCollection L1TMuonBarrelKalmanAlgo::cleanAndSort	(	const L1MuKBMTrackCollection &	tracks,
		uint	keep
	)

Definition at line 829 of file L1TMuonBarrelKalmanAlgo.cc.

References funct::abs(), chiSquareCut_, chiSquareCutCurv_, chiSquareCutPattern_, mps_fire::i, createfilelist::int, MillePedeFileConverter_cfg::out, MCScenario_CRAFT1_22X::sorter(), parallelization::uint(), verbose_, and TriggerAnalyzer::veto.

Referenced by chain(), L1TMuonBarrelKalmanRegionModule::process(), and L1TMuonBarrelKalmanSectorProcessor::process().

                                                                                                           {
  L1MuKBMTrackCollection out;

  if (verbose_) 
    printf(" -----Preselected Kalman Tracks-----\n");


  L1MuKBMTrackCollection pretracks;

  for(const auto& track1 : tracks) {
    if (verbose_)
      printf("Preselected Kalman Track charge=%d pt=%f eta=%f phi=%f curvature=%d curvature STA =%d stubs=%d chi2=%d pts=%f %f\n",track1.charge(),track1.pt(),track1.eta(),track1.phi(),track1.curvatureAtVertex(),track1.curvatureAtMuon(),int(track1.stubs().size()),track1.approxChi2(),track1.pt(),track1.ptUnconstrained()); 

  //Clean up
  bool veto=false;
  for (uint i=0;i<chiSquareCutPattern_.size();++i) {
    if (track1.hitPattern()==chiSquareCutPattern_[i] && abs(track1.curvatureAtVertex())<chiSquareCutCurv_[i] && track1.approxChi2()>chiSquareCut_[i] && track1.curvature()*track1.dxy()<0) {
    veto=true;
    break; 
      }
  }
  if (!veto)
    pretracks.push_back(track1);
  }

  for(const auto& track1 : pretracks) {
    bool keep=true;
    for(const auto& track2 : pretracks) {
      if (track1==track2)
    continue;
      if (!track1.overlapTrack(track2))
    continue;
      if (track1.rank()<track2.rank())
    keep=false;
    }
    if (keep) 
      out.push_back(track1);
  }

  if (verbose_) {
  printf(" -----Algo Result Kalman Tracks-----\n");
  for (const auto& track1 :out)
    printf("Final Kalman Track charge=%d pt=%f eta=%f phi=%f curvature=%d curvature STA =%d stubs=%d chi2=%d pts=%f %f\n",track1.charge(),track1.pt(),track1.eta(),track1.phi(),track1.curvatureAtVertex(),track1.curvatureAtMuon(),int(track1.stubs().size()),track1.approxChi2(),track1.pt(),track1.ptUnconstrained()); 
  }



  TrackSorter sorter;
  if (!out.empty())
    std::sort(out.begin(),out.end(),sorter);


  L1MuKBMTrackCollection exported;
  for (uint i=0;i<out.size();++i)
    if (i<=keep)
      exported.push_back(out[i]);
  return exported;
}

int L1TMuonBarrelKalmanAlgo::correctedPhi ( const L1MuKBMTCombinedStubRef &  stub, int  sector  )

private

Definition at line 155 of file L1TMuonBarrelKalmanAlgo.cc.

Referenced by chain(), estimateChiSquare(), match(), updateLUT(), updateOffline(), and updateOffline1D().

                                                                                        {
  if (stub->scNum()==sector) {
    return stub->phi();
  }
  else if ((stub->scNum()==sector-1) || (stub->scNum()==11 && sector==0)) {
    return stub->phi()-2144;
  }
  else if ((stub->scNum()==sector+1) || (stub->scNum()==0 && sector==11)) {
    return stub->phi()+2144;
  }
  return stub->phi();
} 

int L1TMuonBarrelKalmanAlgo::correctedPhiB ( const L1MuKBMTCombinedStubRef &  stub )

private

Definition at line 149 of file L1TMuonBarrelKalmanAlgo.cc.

Referenced by chain(), estimateChiSquare(), updateLUT(), and updateOffline().

                                                                              {
  //Promote phiB to 12 bits
  return 8*stub->phiB();

}

int L1TMuonBarrelKalmanAlgo::customBitmask ( unsigned int  bit1, unsigned int  bit2, unsigned int  bit3, unsigned int  bit4  )

private

Definition at line 179 of file L1TMuonBarrelKalmanAlgo.cc.

Referenced by rank().

                                                                                                                   {
  return bit1*1+bit2*2+bit3*4+bit4*8;
}

int L1TMuonBarrelKalmanAlgo::encode ( bool  ownwheel, int  sector, bool  tag  )

private

Definition at line 889 of file L1TMuonBarrelKalmanAlgo.cc.

Referenced by trackAddress().

                                                                     {
  if (ownwheel) {
    if (sector==0) {
      if (tag)
    return 8;
      else 
    return 9;
    }
    else if (sector==1) {
      if (tag)
    return 10;
      else 
    return 11;

    }
    else {
      if (tag)
    return 12;
      else 
    return 13;
    }

  }
  else {
    if (sector==0) {
      if (tag)
    return 0;
      else 
    return 1;
    }
    else if (sector==1) {
      if (tag)
    return 2;
      else 
    return 3;

    }
    else {
      if (tag)
    return 4;
      else 
    return 5;
    }
  }
  return 15;
} 

void L1TMuonBarrelKalmanAlgo::estimateChiSquare ( L1MuKBMTrack &  track )

private

Definition at line 768 of file L1TMuonBarrelKalmanAlgo.cc.

References funct::abs(), chiSquare_, correctedPhi(), correctedPhiB(), L1MuKBMTrack::curvatureAtMuon(), delta, L1MuKBMTrack::phiAtMuon(), L1MuKBMTrack::phiBAtMuon(), L1MuKBMTrack::sector(), L1MuKBMTrack::setApproxChi2(), L1MuKBMTrack::stubs(), and parallelization::uint().

Referenced by chain().

                                                                   {
  //here we have a simplification of the algorithm for the sake of the emulator - rsult is identical
  // we apply cuts on the firmware as |u -u'|^2 < a+b *K^2 
  int K = track.curvatureAtMuon();

  int chi=0;

  for (const auto& stub: track.stubs()) {
    uint delta=abs(correctedPhi(stub,track.sector())-track.phiAtMuon()+correctedPhiB(stub)-track.phiBAtMuon()-chiSquare_[stub->stNum()-1]*K);
     chi=chi+delta;    
   }
  chi=chi/2;
  if (chi>511)
    chi=511;
   track.setApproxChi2(chi);
}

bool L1TMuonBarrelKalmanAlgo::getBit ( int  bitmask, int  pos  )

private

Definition at line 183 of file L1TMuonBarrelKalmanAlgo.cc.

Referenced by chain().

                                                         {
  return (bitmask & ( 1 << pos )) >> pos;
}

std::pair< bool, uint > L1TMuonBarrelKalmanAlgo::getByCode ( const L1MuKBMTrackCollection &  tracks, int  mask  )

private

Definition at line 40 of file L1TMuonBarrelKalmanAlgo.cc.

References hitPattern(), mps_fire::i, and parallelization::uint().

                                                                                                   {
  for (uint i=0;i<tracks.size();++i) {
    printf("Code=%d, track=%d\n",tracks[i].hitPattern(),mask); 
    if (tracks[i].hitPattern()==mask)
      return std::make_pair(true,i);
  }
  return std::make_pair(false,0);
}

int L1TMuonBarrelKalmanAlgo::hitPattern ( const L1MuKBMTrack &  track )

private

Definition at line 169 of file L1TMuonBarrelKalmanAlgo.cc.

References RecoTauDiscriminantConfiguration::mask, funct::pow(), and L1MuKBMTrack::stubs().

Referenced by chain(), getByCode(), rank(), updateLUT(), updateOffline(), and updateOffline1D().

                                                                  {
  unsigned int mask = 0;
  for (const auto& stub : track.stubs()) {
    mask = mask+round(pow(2,stub->stNum()-1));
  }
  return mask;
}

std::pair< bool, uint > L1TMuonBarrelKalmanAlgo::match ( const L1MuKBMTrack &  track, const L1MuKBMTCombinedStubRefVector &  stubs  )

private

Definition at line 127 of file L1TMuonBarrelKalmanAlgo.cc.

References funct::abs(), correctedPhi(), SoftLeptonByDistance_cfi::distance, runEdmFileComparison::found, N, L1MuKBMTrack::positionAngle(), L1MuKBMTrack::sector(), L1MuKBMTrack::step(), and parallelization::uint().

Referenced by chain().

                                                                                                                       {
  L1MuKBMTCombinedStubRefVector selected;

  bool found=false;
  uint best=0;
  int distance=100000;
  uint N=0;
  for (const auto& stub :stubs)  {
    N=N+1;
    if (stub->stNum()!=track.step()) 
      continue;
    if (abs(track.positionAngle()-correctedPhi(stub,track.sector()))<distance) {
      distance = abs(track.positionAngle()-correctedPhi(stub,track.sector()));
      best=N-1;
      found=true;
    }
  }
  return std::make_pair(found,best);
}

void L1TMuonBarrelKalmanAlgo::propagate ( L1MuKBMTrack &  track )

private

Definition at line 189 of file L1TMuonBarrelKalmanAlgo.cc.

References a, funct::abs(), aPhi_, aPhiB_, aPhiBNLO_, b, L1MuKBMTrack::bendingAngle(), bPhi_, bPhiB_, ALCARECOTkAlJpsiMuMu_cff::charge, L1MuKBMTrack::covariance(), L1MuKBMTrack::curvature(), eLoss_, mScatteringPhi_, mScatteringPhiB_, phi, L1MuKBMTrack::positionAngle(), L1MuKBMTrack::setCoordinates(), L1MuKBMTrack::setCovariance(), L1MuKBMTrack::step(), verbose_, and wrapAround().

Referenced by chain().

                                                           {
  int K = track.curvature();
  int phi = track.positionAngle();
  int phiB = track.bendingAngle();
  unsigned int step = track.step();

  int charge=1;
  if (K!=0) 
    charge = K/abs(K);

  //energy loss term only for MU->VERTEX
  //int offset=int(charge*eLoss_[step-1]*K*K);
  //  if (abs(offset)>4096)
  //      offset=4096*offset/abs(offset);
  int KNew =wrapAround(int(K/(1+charge*eLoss_[step-1]*K)),8192);

  //phi propagation
  int phiNew =wrapAround(phi+int(aPhi_[step-1]*K)-int(bPhi_[step-1]*phiB),8192);

  //phiB propagation
  int phiBNew = wrapAround(int(aPhiB_[step-1]*K)  +int(bPhiB_[step-1]*phiB),2048);
  
  //Only for the propagation to vertex we use the LUT for better precision and the full function
  if (step==1)
    phiBNew = wrapAround(int(aPhiB_[step-1]*K/(1+charge*aPhiBNLO_[step-1]*K))+int(bPhiB_[step-1]*phiB),2048);

  //Rest of the stuff  is for the offline version only 
  //where we want to check what is happening in the covariaznce matrix 
  
  //Create the transformation matrix
  double a[9];
  a[0] = 1.;
  a[1] = 0.0;
  a[2] = 0.0;
  a[3] = aPhi_[step-1];
  //  a[3] = 0.0;
  a[4] = 1.0;
  a[5] = -bPhi_[step-1];
  //a[6]=0.0;
  a[6] = aPhiB_[step-1];
  a[7] = 0.0;
  a[8] = bPhiB_[step-1];


  ROOT::Math::SMatrix<double,3> P(a,9);

  const std::vector<double>& covLine = track.covariance();
  L1MuKBMTrack::CovarianceMatrix cov(covLine.begin(),covLine.end());
  cov = ROOT::Math::Similarity(P,cov);

  
  //Add the multiple scattering
  double phiRMS = mScatteringPhi_[step-1]*K*K;
  double phiBRMS = mScatteringPhiB_[step-1]*K*K;

  std::vector<double> b(6);
  b[0] = 0;
  b[1] = 0;
  b[2] =phiRMS;
  b[3] =0;
  b[4] = 0;
  b[5] = phiBRMS;

  reco::Candidate::CovarianceMatrix MS(b.begin(),b.end());

  cov = cov+MS;

  if (verbose_) {
    printf("Covariance term for phiB = %f\n",cov(2,2));
    printf("Multiple scattering term for phiB = %f\n",MS(2,2));
  }
 


  track.setCovariance(cov);
  track.setCoordinates(step-1,KNew,phiNew,phiBNew);

}

int L1TMuonBarrelKalmanAlgo::rank ( const L1MuKBMTrack &  track )

private

Definition at line 786 of file L1TMuonBarrelKalmanAlgo.cc.

References L1MuKBMTrack::approxChi2(), customBitmask(), hitPattern(), and L1MuKBMTrack::stubs().

Referenced by addBMTFMuon(), and chain().

                                                           {
  //    int offset=0;
    if (hitPattern(track)==customBitmask(0,0,1,1))
      return 65;
    //    return offset+(track.stubs().size()*2+track.quality())*80-track.approxChi2();
    return 700+(track.stubs().size())*80-track.approxChi2();

}

void L1TMuonBarrelKalmanAlgo::resolveEtaUnit

(

L1MuKBMTrackCollection &

tracks

)

Definition at line 806 of file L1TMuonBarrelKalmanAlgo.cc.

References setFloatingPointValues(), HiIsolationCommonParameters_cff::track, and parallelization::uint().

Referenced by chain(), and L1TMuonBarrelKalmanRegionModule::process().

                                                                           {
  int bestFineEta=0;
  uint bestSegments=0;
  for (const auto& track : tracks) {
    if (track.stubs().size()>bestSegments && track.hasFineEta()) {
      bestFineEta = track.fineEta();
      bestSegments=track.stubs().size();
    }
  }
  
  for (auto & track :tracks) {
    if (bestSegments!=0 && (!track.hasFineEta())) {
      track.setFineEta(bestFineEta);
      setFloatingPointValues(track,true);
    }
  }



}

void L1TMuonBarrelKalmanAlgo::setFloatingPointValues ( L1MuKBMTrack &  track, bool  vertex  )

private

Definition at line 567 of file L1TMuonBarrelKalmanAlgo.cc.

References funct::abs(), L1MuKBMTrack::coarseEta(), L1MuKBMTrack::curvatureAtMuon(), L1MuKBMTrack::curvatureAtVertex(), MillePedeFileConverter_cfg::e, PVValHelper::eta, L1MuKBMTrack::fineEta(), objects.autophobj::float, L1MuKBMTrack::hasFineEta(), M_PI, phi, L1MuKBMTrack::phiAtVertex(), EnergyCorrector::pt, L1MuKBMTrack::sector(), reco::LeafCandidate::setCharge(), L1MuKBMTrack::setPtEtaPhi(), and L1MuKBMTrack::setPtUnconstrained().

Referenced by chain(), and resolveEtaUnit().

                                                                                    {
  int K,phiINT,etaINT;

  if (track.hasFineEta())
    etaINT=track.fineEta();
  else
    etaINT=track.coarseEta();


  double lsb = 1.25/float(1 << 13);
  double lsbEta = 0.010875;


  if (vertex) {
    K  = track.curvatureAtVertex();
    if (K==0)
      track.setCharge(1);
    else
      track.setCharge(K/abs(K));

    phiINT = track.phiAtVertex();
    double phi= track.sector()*M_PI/6.0+phiINT*M_PI/(6*2048.)-2*M_PI;
    double eta = etaINT*lsbEta;
    if (phi<-M_PI)
      phi=phi+2*M_PI;
    if (K==0)
      K=1;    


    float FK=fabs(K);
    FK = fabs(0.912*FK+(2.557e-5)*FK*FK-6);

    if (FK<51)
      FK=51;
    
    double pt = 1.0/(lsb*(FK));

    track.setPtEtaPhi(pt,eta,phi);
  }
  else {
    K=track.curvatureAtMuon();
    if (K==0)
      K=1;
    if (abs(K)<46)
      K=46*K/abs(K);
    double pt = 1.0/(lsb*abs(K));
    track.setPtUnconstrained(pt);
  }
}

std::map< int, int > L1TMuonBarrelKalmanAlgo::trackAddress ( const L1MuKBMTrack &  track, int &  word  )

private

Definition at line 939 of file L1TMuonBarrelKalmanAlgo.cc.

References funct::abs(), encode(), l1t::RegionalMuonCand::kNumBmtfSubAddr, l1t::RegionalMuonCand::kSegSelStat1, l1t::RegionalMuonCand::kSegSelStat2, l1t::RegionalMuonCand::kSegSelStat3, l1t::RegionalMuonCand::kSegSelStat4, l1t::RegionalMuonCand::kStat1, l1t::RegionalMuonCand::kStat2, l1t::RegionalMuonCand::kStat3, l1t::RegionalMuonCand::kStat4, l1t::RegionalMuonCand::kWheelNum, l1t::RegionalMuonCand::kWheelSide, MillePedeFileConverter_cfg::out, L1MuKBMTrack::sector(), L1MuKBMTrack::stubs(), and L1MuKBMTrack::wheel().

Referenced by addBMTFMuon().

                                                                                         {
  std::map<int,int> out;
  if (track.wheel()>=0)
    out[l1t::RegionalMuonCand::kWheelSide] = 0;
  else
    out[l1t::RegionalMuonCand::kWheelSide] = 1;

  out[l1t::RegionalMuonCand::kWheelNum] = abs(track.wheel());
  out[l1t::RegionalMuonCand::kStat1]=3;
  out[l1t::RegionalMuonCand::kStat2]=15;
  out[l1t::RegionalMuonCand::kStat3]=15;
  out[l1t::RegionalMuonCand::kStat4]=15;
  out[l1t::RegionalMuonCand::kSegSelStat1]=0;
  out[l1t::RegionalMuonCand::kSegSelStat2]=0;
  out[l1t::RegionalMuonCand::kSegSelStat3]=0;
  out[l1t::RegionalMuonCand::kSegSelStat4]=0;
  out[l1t::RegionalMuonCand::kNumBmtfSubAddr]=0;


  for (const auto stub: track.stubs()) {
    bool ownwheel = stub->whNum() == track.wheel();
    int sector=0;
    if ((stub->scNum()==track.sector()+1) || (stub->scNum()==0 && track.sector()==11))
      sector=+1;
    if ((stub->scNum()==track.sector()-1) || (stub->scNum()==11 && track.sector()==0))
      sector=-1;
    int addr = encode(ownwheel,sector,stub->tag());
   
    if (stub->stNum()==4) {
      addr=addr & 3;
      out[l1t::RegionalMuonCand::kStat1]=addr;
    }      
    if (stub->stNum()==3) {
      out[l1t::RegionalMuonCand::kStat2]=addr;    
    }
    if (stub->stNum()==2) {
      out[l1t::RegionalMuonCand::kStat3]=addr;
    }
    if (stub->stNum()==1) {
      out[l1t::RegionalMuonCand::kStat4]=addr;
    }
  }
    
  word=0;
  word = word | out[l1t::RegionalMuonCand::kStat1]<<12;
  word = word | out[l1t::RegionalMuonCand::kStat2]<<8;
  word = word | out[l1t::RegionalMuonCand::kStat3]<<4;
  word = word | out[l1t::RegionalMuonCand::kStat4];

  

  return out;
}

uint L1TMuonBarrelKalmanAlgo::twosCompToBits ( int  q )

private

Definition at line 995 of file L1TMuonBarrelKalmanAlgo.cc.

References lumiQueryAPI::q.

Referenced by addBMTFMuon().

                                                  {
  if (q>=0)
    return q;
  else 
    return (~q)+1; 


}

bool L1TMuonBarrelKalmanAlgo::update ( L1MuKBMTrack &  track, const L1MuKBMTCombinedStubRef &  stub, int  mask  )

private

Definition at line 270 of file L1TMuonBarrelKalmanAlgo.cc.

References updateEta(), updateLUT(), updateOffline(), updateOffline1D(), and useOfflineAlgo_.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), chain(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                                                                                     {
  updateEta(track,stub);
  if (useOfflineAlgo_) {
    if (mask==3 || mask ==5 || mask==9 ||mask==6|| mask==10 ||mask==12)
          return updateOffline(track,stub);
        else
      return updateOffline1D(track,stub);

  }
  else
    return updateLUT(track,stub,mask);

}

void L1TMuonBarrelKalmanAlgo::updateEta ( L1MuKBMTrack &  track, const L1MuKBMTCombinedStubRef &  stub  )

private

Definition at line 472 of file L1TMuonBarrelKalmanAlgo.cc.

References funct::abs(), L1MuKBMTrack::fineEta(), L1MuKBMTrack::hasFineEta(), L1MuKBMTrack::setFineEta(), and parallelization::uint().

Referenced by update().

                                                                                               {

    if (stub->qeta1()>=0) {
      if (track.hasFineEta()) {
    uint dist2=1000;
    uint dist1 = abs(track.fineEta()-stub->eta1());
    if (stub->qeta2()>=0)
      dist2 = abs(track.fineEta()-stub->eta2());
    if (dist1<dist2)
      track.setFineEta((stub->eta1()+track.fineEta())/2);
    else
      track.setFineEta((stub->eta2()+track.fineEta())/2);
      }else {
    if (stub->qeta2()>=0)
      track.setFineEta((stub->eta1()+stub->eta2())/2);
    else
      track.setFineEta(stub->eta1());
      }

    }
}

bool L1TMuonBarrelKalmanAlgo::updateLUT ( L1MuKBMTrack &  track, const L1MuKBMTCombinedStubRef &  stub, int  mask  )

private

Definition at line 428 of file L1TMuonBarrelKalmanAlgo.cc.

References funct::abs(), L1MuKBMTrack::addStub(), L1MuKBMTrack::bendingAngle(), L1MuKBMTrack::coarseEta(), correctedPhi(), correctedPhiB(), L1MuKBMTrack::curvature(), hitPattern(), lutService_, RecoTauDiscriminantConfiguration::mask, phi, L1MuKBMTrack::positionAngle(), L1MuKBMTrack::sector(), L1MuKBMTrack::setCoarseEta(), L1MuKBMTrack::setCoordinates(), L1MuKBMTrack::setHitPattern(), L1MuKBMTrack::setKalmanGain(), L1MuKBMTrack::step(), reco::btau::trackPhi, parallelization::uint(), and wrapAround().

Referenced by update().

                                                                                                        {
    int trackK = track.curvature();
    int trackPhi = track.positionAngle();
    int trackPhiB = track.bendingAngle();

    //Update eta
    track.setCoarseEta(int((track.coarseEta()+stub->coarseEta())/2.0));

    int phi  = correctedPhi(stub,track.sector());
    int phiB = correctedPhiB(stub);
    //    if (stub->quality()<6)
    //      phiB=trackPhiB;

    Vector2 residual;
    residual[0] = phi-trackPhi;
    residual[1] = phiB-trackPhiB;

    uint absK = abs(trackK);
    if (absK>4095)
      absK = 4095;
    std::vector<float> GAIN = lutService_->trackGain(track.step(),mask,absK/4);
    track.setKalmanGain(track.step(),abs(trackK),GAIN[0],GAIN[1],1,0,GAIN[2],GAIN[3]);

    //For the three stub stuff use only gains 0 and 4
    if (!(mask==3 || mask ==5 || mask==9 ||mask==6|| mask==10 ||mask==12))  {
      GAIN[1]=0.0;
      GAIN[3]=0.0;
    }
      



    int KNew  = wrapAround(trackK+int(GAIN[0]*residual(0)+GAIN[1]*residual(1)),8192);
    int phiNew  = wrapAround(trackPhi+residual(0),8192);
    int phiBNew = wrapAround(trackPhiB+int(GAIN[2]*residual(0)+GAIN[3]*residual(1)),2048);
    track.setCoordinates(track.step(),KNew,phiNew,phiBNew);
    track.addStub(stub);
    track.setHitPattern(hitPattern(track));
    return true;
}

bool L1TMuonBarrelKalmanAlgo::updateOffline ( L1MuKBMTrack &  track, const L1MuKBMTCombinedStubRef &  stub  )

private

Definition at line 284 of file L1TMuonBarrelKalmanAlgo.cc.

References funct::abs(), L1MuKBMTrack::addStub(), L1MuKBMTrack::bendingAngle(), EnergyCorrector::c, L1MuKBMTrack::coarseEta(), correctedPhi(), correctedPhiB(), L1MuKBMTrack::covariance(), L1MuKBMTrack::curvature(), g4SimHits_cfi::Gain, class-composition::H, hitPattern(), createfilelist::int, phi, pointResolutionPhi_, pointResolutionPhiB_, L1MuKBMTrack::positionAngle(), dttmaxenums::R, S(), L1MuKBMTrack::sector(), L1MuKBMTrack::setCoarseEta(), L1MuKBMTrack::setCoordinates(), L1MuKBMTrack::setCovariance(), L1MuKBMTrack::setHitPattern(), L1MuKBMTrack::setKalmanGain(), L1MuKBMTrack::setResidual(), L1MuKBMTrack::step(), reco::btau::trackPhi, verbose_, and wrapAround().

Referenced by update().

                                                                                                   {
    int trackK = track.curvature();
    int trackPhi = track.positionAngle();
    int trackPhiB = track.bendingAngle();

    int phi  = correctedPhi(stub,track.sector());
    int phiB = correctedPhiB(stub);



    //Update eta
    track.setCoarseEta(int((track.coarseEta()+stub->coarseEta())/2.0));

    Vector2 residual;
    residual[0] = phi-trackPhi;
    residual[1] = phiB-trackPhiB;

   
   

    //    if (stub->quality()<4)
    //  phiB=trackPhiB;

    Matrix23 H;
    H(0,0)=0.0;
    H(0,1)=1.0;
    H(0,2)=0.0;
    H(1,0)=0.0;
    H(1,1)=0.0;
    H(1,2)=1.0;

    
    CovarianceMatrix2 R;
    R(0,0) = pointResolutionPhi_;
    R(0,1) = 0.0;
    R(1,0) = 0.0;
    R(1,1) = pointResolutionPhiB_;

    const std::vector<double>& covLine = track.covariance();
    L1MuKBMTrack::CovarianceMatrix cov(covLine.begin(),covLine.end());


    CovarianceMatrix2 S = ROOT::Math::Similarity(H,cov)+R;
    if (!S.Invert())
      return false;
    Matrix32 Gain = cov*ROOT::Math::Transpose(H)*S;

    track.setKalmanGain(track.step(),abs(trackK),Gain(0,0),Gain(0,1),Gain(1,0),Gain(1,1),Gain(2,0),Gain(2,1));

    int KNew  = (trackK+int(Gain(0,0)*residual(0)+Gain(0,1)*residual(1)));
    if (abs(KNew)>8192)
      return false;
    
    int phiNew  = wrapAround(trackPhi+residual(0),8192);
    int phiBNew = wrapAround(trackPhiB+int(Gain(2,0)*residual(0)+Gain(2,1)*residual(1)),2048);
    
    track.setResidual(stub->stNum()-1,abs(phi-phiNew)+abs(phiB-phiBNew)/8);


    if (verbose_) {
      printf(" K = %d + %f * %f + %f * %f\n",trackK,Gain(0,0),residual(0),Gain(0,1),residual(1));
      printf(" phiB = %d + %f * %f + %f * %f\n",trackPhiB,Gain(2,0),residual(0),Gain(2,1),residual(1));
    }


    track.setCoordinates(track.step(),KNew,phiNew,phiBNew);
    Matrix33 covNew = cov - Gain*(H*cov);
    L1MuKBMTrack::CovarianceMatrix c;
 
    c(0,0)=covNew(0,0); 
    c(0,1)=covNew(0,1); 
    c(0,2)=covNew(0,2); 
    c(1,0)=covNew(1,0); 
    c(1,1)=covNew(1,1); 
    c(1,2)=covNew(1,2); 
    c(2,0)=covNew(2,0); 
    c(2,1)=covNew(2,1); 
    c(2,2)=covNew(2,2); 
    if (verbose_) {
      printf("Post Fit Covariance Matrix %f %f %f \n",cov(0,0),cov(1,1),cov(2,2));
      
    }

    track.setCovariance(c);
    track.addStub(stub);
    track.setHitPattern(hitPattern(track));

    return true;
}

bool L1TMuonBarrelKalmanAlgo::updateOffline1D ( L1MuKBMTrack &  track, const L1MuKBMTCombinedStubRef &  stub  )

private

Definition at line 375 of file L1TMuonBarrelKalmanAlgo.cc.

References funct::abs(), L1MuKBMTrack::addStub(), L1MuKBMTrack::bendingAngle(), EnergyCorrector::c, L1MuKBMTrack::coarseEta(), correctedPhi(), L1MuKBMTrack::covariance(), L1MuKBMTrack::curvature(), g4SimHits_cfi::Gain, class-composition::H, hitPattern(), phi, pointResolutionPhi_, L1MuKBMTrack::positionAngle(), S(), L1MuKBMTrack::sector(), L1MuKBMTrack::setCoarseEta(), L1MuKBMTrack::setCoordinates(), L1MuKBMTrack::setCovariance(), L1MuKBMTrack::setHitPattern(), L1MuKBMTrack::setKalmanGain(), L1MuKBMTrack::step(), reco::btau::trackPhi, and wrapAround().

Referenced by update().

                                                                                                     {
    int trackK = track.curvature();
    int trackPhi = track.positionAngle();
    int trackPhiB = track.bendingAngle();


    int phi  = correctedPhi(stub,track.sector());
    track.setCoarseEta(int((track.coarseEta()+stub->coarseEta())/2.0));

    double residual= phi-trackPhi;

    Matrix13 H;
    H(0,0)=0.0;
    H(0,1)=1.0;
    H(0,2)=0.0;
    

    const std::vector<double>& covLine = track.covariance();
    L1MuKBMTrack::CovarianceMatrix cov(covLine.begin(),covLine.end());

    double S = ROOT::Math::Similarity(H,cov)(0,0)+pointResolutionPhi_;

    if (S==0.0)
      return false;
    Matrix31 Gain = cov*ROOT::Math::Transpose(H)/S;

    track.setKalmanGain(track.step(),abs(trackK),Gain(0,0),0.0,Gain(1,0),0.0,Gain(2,0),0.0);

    int KNew  = wrapAround(trackK+int(Gain(0,0)*residual),8192);
    int phiNew  = wrapAround(trackPhi+residual,8192);
    int phiBNew = wrapAround(trackPhiB+int(Gain(2,0)*residual),2048);
    track.setCoordinates(track.step(),KNew,phiNew,phiBNew);
    Matrix33 covNew = cov - Gain*(H*cov);
    L1MuKBMTrack::CovarianceMatrix c;
 
    c(0,0)=covNew(0,0); 
    c(0,1)=covNew(0,1); 
    c(0,2)=covNew(0,2); 
    c(1,0)=covNew(1,0); 
    c(1,1)=covNew(1,1); 
    c(1,2)=covNew(1,2); 
    c(2,0)=covNew(2,0); 
    c(2,1)=covNew(2,1); 
    c(2,2)=covNew(2,2); 
    track.setCovariance(c);
    track.addStub(stub);
    track.setHitPattern(hitPattern(track));

    return true;
}

void L1TMuonBarrelKalmanAlgo::vertexConstraint ( L1MuKBMTrack &  track )

private

Definition at line 499 of file L1TMuonBarrelKalmanAlgo.cc.

References useOfflineAlgo_, vertexConstraintLUT(), and vertexConstraintOffline().

Referenced by chain().

                                                                  {
  if (useOfflineAlgo_)
    vertexConstraintOffline(track);
  else
    vertexConstraintLUT(track);
 
}

void L1TMuonBarrelKalmanAlgo::vertexConstraintLUT ( L1MuKBMTrack &  track )

private

Definition at line 551 of file L1TMuonBarrelKalmanAlgo.cc.

References funct::abs(), L1MuKBMTrack::curvature(), L1MuKBMTrack::dxy(), L1MuKBMTrack::hitPattern(), lutService_, L1MuKBMTrack::positionAngle(), L1MuKBMTrack::setCoordinatesAtVertex(), L1MuKBMTrack::setKalmanGain(), L1MuKBMTrack::step(), parallelization::uint(), and wrapAround().

Referenced by vertexConstraint().

                                                                     {
  double residual = -track.dxy();
  uint absK = abs(track.curvature());
  if (absK>2047)
    absK = 2047;

std::pair<float,float> GAIN = lutService_->vertexGain(track.hitPattern(),absK/2);
  track.setKalmanGain(track.step(),abs(track.curvature()),GAIN.first,GAIN.second,-1);

  int KNew = wrapAround(int(track.curvature()+GAIN.first*residual),8192);
  int phiNew = wrapAround(int(track.positionAngle()+GAIN.second*residual),8192);
  track.setCoordinatesAtVertex(KNew,phiNew,-residual);
}

void L1TMuonBarrelKalmanAlgo::vertexConstraintOffline ( L1MuKBMTrack &  track )

private

Definition at line 508 of file L1TMuonBarrelKalmanAlgo.cc.

References funct::abs(), EnergyCorrector::c, L1MuKBMTrack::covariance(), L1MuKBMTrack::curvature(), L1MuKBMTrack::dxy(), g4SimHits_cfi::Gain, class-composition::H, pointResolutionVertex_, L1MuKBMTrack::positionAngle(), S(), L1MuKBMTrack::setCoordinatesAtVertex(), L1MuKBMTrack::setCovariance(), L1MuKBMTrack::setKalmanGain(), L1MuKBMTrack::step(), verbose_, and wrapAround().

Referenced by vertexConstraint().

                                                                         {
  double residual = -track.dxy();
  Matrix13 H;
  H(0,0)=0;
  H(0,1)=0;
  H(0,2)=1;

  const std::vector<double>& covLine = track.covariance();
  L1MuKBMTrack::CovarianceMatrix cov(covLine.begin(),covLine.end());
  
  double S = (ROOT::Math::Similarity(H,cov))(0,0)+pointResolutionVertex_;
  S=1.0/S;
  Matrix31 Gain = cov*(ROOT::Math::Transpose(H))*S;
  track.setKalmanGain(track.step(),abs(track.curvature()),Gain(0,0),Gain(1,0),Gain(2,0));

  if (verbose_) {
    printf("sigma3=%f sigma6=%f\n",cov(0,3),cov(3,3));
    printf(" K = %d + %f * %f\n",track.curvature(),Gain(0,0),residual);
  }

  int KNew = wrapAround(int(track.curvature()+Gain(0,0)*residual),8192);
  int phiNew = wrapAround(int(track.positionAngle()+Gain(1,0)*residual),8192);
  int dxyNew = wrapAround(int(track.dxy()+Gain(2,0)*residual),8192);
  if (verbose_)
    printf("Post fit impact parameter=%d\n",dxyNew);
  track.setCoordinatesAtVertex(KNew,phiNew,-residual);
  Matrix33 covNew = cov - Gain*(H*cov);
  L1MuKBMTrack::CovarianceMatrix c;
  c(0,0)=covNew(0,0); 
  c(0,1)=covNew(0,1); 
  c(0,2)=covNew(0,2); 
  c(1,0)=covNew(1,0); 
  c(1,1)=covNew(1,1); 
  c(1,2)=covNew(1,2); 
  c(2,0)=covNew(2,0); 
  c(2,1)=covNew(2,1); 
  c(2,2)=covNew(2,2); 
  track.setCovariance(c);
  //  track.covariance = track.covariance - Gain*H*track.covariance;
}

int L1TMuonBarrelKalmanAlgo::wrapAround ( int  value, int  maximum  )

private

Definition at line 797 of file L1TMuonBarrelKalmanAlgo.cc.

References relativeConstraints::value.

Referenced by propagate(), updateLUT(), updateOffline(), updateOffline1D(), vertexConstraintLUT(), and vertexConstraintOffline().

                                                             {
  if (value>maximum-1)
    return value-2*maximum;
  if (value<-maximum)
    return value+2*maximum;
  return value;

}

Member Data Documentation

std::vector<double> L1TMuonBarrelKalmanAlgo::aPhi_

private

Definition at line 74 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by propagate().

std::vector<double> L1TMuonBarrelKalmanAlgo::aPhiB_

private

Definition at line 75 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by propagate().

std::vector<double> L1TMuonBarrelKalmanAlgo::aPhiBNLO_

private

Definition at line 76 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by propagate().

std::vector<double> L1TMuonBarrelKalmanAlgo::bPhi_

private

Definition at line 77 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by propagate().

std::vector<double> L1TMuonBarrelKalmanAlgo::bPhiB_

private

Definition at line 78 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by propagate().

std::vector<double> L1TMuonBarrelKalmanAlgo::chiSquare_

private

Definition at line 84 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by estimateChiSquare().

std::vector<int> L1TMuonBarrelKalmanAlgo::chiSquareCut_

private

Definition at line 87 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by cleanAndSort().

std::vector<int> L1TMuonBarrelKalmanAlgo::chiSquareCutCurv_

private

Definition at line 86 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by cleanAndSort().

std::vector<int> L1TMuonBarrelKalmanAlgo::chiSquareCutPattern_

private

Definition at line 85 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by cleanAndSort().

std::vector<int> L1TMuonBarrelKalmanAlgo::combos1_

private

Definition at line 94 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by chain().

std::vector<int> L1TMuonBarrelKalmanAlgo::combos2_

private

Definition at line 93 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by chain().

std::vector<int> L1TMuonBarrelKalmanAlgo::combos3_

private

Definition at line 92 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by chain().

std::vector<int> L1TMuonBarrelKalmanAlgo::combos4_

private

Definition at line 91 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by chain().

std::vector<double> L1TMuonBarrelKalmanAlgo::eLoss_

private

Definition at line 73 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by propagate().

std::vector<double> L1TMuonBarrelKalmanAlgo::etaLUT0_

private

Definition at line 79 of file L1TMuonBarrelKalmanAlgo.h.

std::vector<double> L1TMuonBarrelKalmanAlgo::etaLUT1_

private

Definition at line 80 of file L1TMuonBarrelKalmanAlgo.h.

uint L1TMuonBarrelKalmanAlgo::globalChi2Cut_

private

Definition at line 83 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by chain().

std::vector<double> L1TMuonBarrelKalmanAlgo::initK2_

private

Definition at line 70 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by chain().

std::vector<double> L1TMuonBarrelKalmanAlgo::initK_

private

Definition at line 69 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by chain().

std::unique_ptr<L1TMuonBarrelKalmanLUTs> L1TMuonBarrelKalmanAlgo::lutService_

private

Definition at line 65 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by updateLUT(), and vertexConstraintLUT().

std::vector<double> L1TMuonBarrelKalmanAlgo::mScatteringPhi_

private

Definition at line 101 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by propagate().

std::vector<double> L1TMuonBarrelKalmanAlgo::mScatteringPhiB_

private

Definition at line 102 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by propagate().

double L1TMuonBarrelKalmanAlgo::pointResolutionPhi_

private

Definition at line 104 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by chain(), updateOffline(), and updateOffline1D().

double L1TMuonBarrelKalmanAlgo::pointResolutionPhiB_

private

Definition at line 106 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by chain(), and updateOffline().

double L1TMuonBarrelKalmanAlgo::pointResolutionVertex_

private

Definition at line 108 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by vertexConstraintOffline().

bool L1TMuonBarrelKalmanAlgo::useOfflineAlgo_

private

Definition at line 100 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by update(), and vertexConstraint().

bool L1TMuonBarrelKalmanAlgo::verbose_

private

Definition at line 39 of file L1TMuonBarrelKalmanAlgo.h.

Referenced by chain(), cleanAndSort(), propagate(), updateOffline(), and vertexConstraintOffline().