Changeset e40b9cf in git
- Timestamp:
- Oct 30, 2017, 12:56:01 AM (7 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- 9458496b
- Parents:
- cde9f31
- Location:
- modules
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/TrackCountingBTagging.cc
rcde9f31 re40b9cf 72 72 fNtracks = GetInt("Ntracks", 3); 73 73 74 fUse3D = GetBool("Use3D", false); 75 74 76 // import input array(s) 75 77 … … 95 97 Candidate *jet, *track; 96 98 97 Double_t jpx, jpy ;99 Double_t jpx, jpy, jpz; 98 100 Double_t dr, tpt; 99 Double_t xd, yd, d0, dd0, ip, sip;101 Double_t xd, yd, zd, d0, dd0, dz, ddz, sip; 100 102 101 103 Int_t sign; … … 110 112 jpx = jetMomentum.Px(); 111 113 jpy = jetMomentum.Py(); 114 jpz = jetMomentum.Pz(); 112 115 113 116 // loop over all input tracks … … 122 125 xd = track->Xd; 123 126 yd = track->Yd; 127 zd = track->Zd; 124 128 d0 = TMath::Abs(track->D0); 125 dd0 = track->ErrorD0; 129 dd0 = TMath::Abs(track->ErrorD0); 130 dz = TMath::Abs(track->DZ); 131 ddz = TMath::Abs(track->ErrorDZ); 126 132 127 133 if(tpt < fPtMin) continue; … … 129 135 if(d0 > fIPmax) continue; 130 136 131 sign = (jpx*xd + jpy*yd > 0.0) ? 1 : -1; 132 133 ip = sign*d0; 134 sip = ip / TMath::Abs(dd0); 137 if(fUse3D){ 138 sign = (jpx*xd + jpy*yd + jpz*zd > 0.0) ? 1 : -1; 139 //add transvers and longitudinal significances in quadrature 140 sip = sign * TMath::Sqrt( TMath::Power(d0 / dd0, 2) + TMath::Power(dz / ddz, 2) ); 141 } 142 else { 143 sign = (jpx*xd + jpy*yd > 0.0) ? 1 : -1; 144 sip = sign * d0 / TMath::Abs(dd0); 145 } 135 146 136 147 if(sip > fSigMin) count++; -
modules/TrackCountingBTagging.h
rcde9f31 re40b9cf 54 54 Double_t fSigMin; 55 55 Int_t fNtracks; 56 Bool_t fUse3D; 56 57 57 58 TIterator *fItTrackInputArray; //!
Note:
See TracChangeset
for help on using the changeset viewer.