CMS 3D CMS Logo

List of all members | Public Types | Public Member Functions | Private Member Functions | Private Attributes
MuonSeedCreator Class Reference

#include <MuonSeedCreator.h>

Public Types

typedef MuonTransientTrackingRecHit::MuonRecHitContainer SegmentContainer
 

Public Member Functions

TrajectorySeed createSeed (int type, const SegmentContainer &seg, const std::vector< int > &layers, int NShower, int NShowerSeg)
 Create a seed from set of segments. More...
 
 MuonSeedCreator (const edm::ParameterSet &pset)
 Constructor. More...
 
void setBField (const MagneticField *theField)
 Cache Magnetic Field for current event. More...
 
 ~MuonSeedCreator ()
 Destructor. More...
 

Private Member Functions

void estimatePtCSC (const SegmentContainer &seg, const std::vector< int > &layers, double &pt, double &spt)
 Estimate transverse momentum of track from CSC measurements. More...
 
void estimatePtDT (const SegmentContainer &seg, const std::vector< int > &layers, double &pt, double &spt)
 Estimate transverse momentum of track from DT measurements. More...
 
void estimatePtOverlap (const SegmentContainer &seg, const std::vector< int > &layers, double &pt, double &spt)
 Estimate transverse momentum of track from CSC + DT measurements. More...
 
void estimatePtShowering (int &NShowers, int &NShowerSeg, double &pt, double &spt)
 Estimate transverse momentum of track from showering segment. More...
 
void estimatePtSingle (const SegmentContainer &seg, const std::vector< int > &layers, double &pt, double &spt)
 Estimate transverse momentum of track from single segment. More...
 
std::vector< double > getPt (const std::vector< double > &vParameters, double eta, double dPhi)
 Compute pt from parameters. More...
 
double scaledPhi (double dphi, double t1)
 Scale the dPhi from segment position. More...
 
void weightedPt (const std::vector< double > &ptEstimate, const std::vector< double > &sptEstimate, double &ptAvg, double &sptAvg)
 Compute weighted mean pt from different pt estimators. More...
 

Private Attributes

const MagneticFieldBField
 
std::vector< double > CSC01
 
std::vector< double > CSC01_1
 
std::vector< double > CSC02
 
std::vector< double > CSC03
 
std::vector< double > CSC12
 
std::vector< double > CSC12_1
 
std::vector< double > CSC12_2
 
std::vector< double > CSC12_3
 
std::vector< double > CSC13
 
std::vector< double > CSC13_2
 
std::vector< double > CSC13_3
 
std::vector< double > CSC14
 
std::vector< double > CSC14_3
 
std::vector< double > CSC23
 
std::vector< double > CSC23_1
 
std::vector< double > CSC23_2
 
std::vector< double > CSC24
 
std::vector< double > CSC24_1
 
std::vector< double > CSC34
 
std::vector< double > CSC34_1
 
bool debug
 
float defaultMomentum
 
std::vector< double > DT12
 
std::vector< double > DT12_1
 
std::vector< double > DT12_2
 
std::vector< double > DT13
 
std::vector< double > DT13_1
 
std::vector< double > DT13_2
 
std::vector< double > DT14
 
std::vector< double > DT14_1
 
std::vector< double > DT14_2
 
std::vector< double > DT23
 
std::vector< double > DT23_1
 
std::vector< double > DT23_2
 
std::vector< double > DT24
 
std::vector< double > DT24_1
 
std::vector< double > DT24_2
 
std::vector< double > DT34
 
std::vector< double > DT34_1
 
std::vector< double > DT34_2
 
std::vector< double > OL1213
 
std::vector< double > OL1222
 
std::vector< double > OL1232
 
std::vector< double > OL2213
 
std::vector< double > OL2222
 
std::vector< double > OL_1213
 
std::vector< double > OL_1222
 
std::vector< double > OL_1232
 
std::vector< double > OL_2213
 
std::vector< double > OL_2222
 
std::vector< double > SMB10
 
std::vector< double > SMB11
 
std::vector< double > SMB12
 
std::vector< double > SMB20
 
std::vector< double > SMB21
 
std::vector< double > SMB22
 
std::vector< double > SMB30
 
std::vector< double > SMB31
 
std::vector< double > SMB32
 
std::vector< double > SMB_10S
 
std::vector< double > SMB_11S
 
std::vector< double > SMB_12S
 
std::vector< double > SMB_20S
 
std::vector< double > SMB_21S
 
std::vector< double > SMB_22S
 
std::vector< double > SMB_30S
 
std::vector< double > SMB_31S
 
std::vector< double > SMB_32S
 
std::vector< double > SME11
 
std::vector< double > SME12
 
std::vector< double > SME13
 
std::vector< double > SME21
 
std::vector< double > SME22
 
std::vector< double > SME31
 
std::vector< double > SME32
 
std::vector< double > SME41
 
std::vector< double > SME_11S
 
std::vector< double > SME_12S
 
std::vector< double > SME_13S
 
std::vector< double > SME_21S
 
std::vector< double > SME_22S
 
double sysError
 
float theMaxMomentum
 
float theMinMomentum
 

Detailed Description

Creates seed from vector of segment

Determine pt of seed using various combination of segments from different layers (stations) Parameterization used to determine pt between layers i and j:

pt = [ c_0 + c_1 * (Delta phi_ij) + c_2 * (Delta phi_ij)^2 ] / eta

Author
Dominique Fortin - UCR

Definition at line 30 of file MuonSeedCreator.h.

Member Typedef Documentation

◆ SegmentContainer

Definition at line 32 of file MuonSeedCreator.h.

Constructor & Destructor Documentation

◆ MuonSeedCreator()

MuonSeedCreator::MuonSeedCreator ( const edm::ParameterSet pset)
explicit

Constructor.

See header file for a description of this class.

Author
: Shih-Chuan Kao, Dominique Fortin - UCR

Definition at line 35 of file MuonSeedCreator.cc.

35  {
36  theMinMomentum = pset.getParameter<double>("minimumSeedPt");
37  theMaxMomentum = pset.getParameter<double>("maximumSeedPt");
38  defaultMomentum = pset.getParameter<double>("defaultSeedPt");
39  debug = pset.getParameter<bool>("DebugMuonSeed");
40  sysError = pset.getParameter<double>("SeedPtSystematics");
41  // load seed PT parameters
42  DT12 = pset.getParameter<std::vector<double> >("DT_12");
43  DT13 = pset.getParameter<std::vector<double> >("DT_13");
44  DT14 = pset.getParameter<std::vector<double> >("DT_14");
45  DT23 = pset.getParameter<std::vector<double> >("DT_23");
46  DT24 = pset.getParameter<std::vector<double> >("DT_24");
47  DT34 = pset.getParameter<std::vector<double> >("DT_34");
48 
49  CSC01 = pset.getParameter<std::vector<double> >("CSC_01");
50  CSC12 = pset.getParameter<std::vector<double> >("CSC_12");
51  CSC02 = pset.getParameter<std::vector<double> >("CSC_02");
52  CSC13 = pset.getParameter<std::vector<double> >("CSC_13");
53  CSC03 = pset.getParameter<std::vector<double> >("CSC_03");
54  CSC14 = pset.getParameter<std::vector<double> >("CSC_14");
55  CSC23 = pset.getParameter<std::vector<double> >("CSC_23");
56  CSC24 = pset.getParameter<std::vector<double> >("CSC_24");
57  CSC34 = pset.getParameter<std::vector<double> >("CSC_34");
58 
59  OL1213 = pset.getParameter<std::vector<double> >("OL_1213");
60  OL1222 = pset.getParameter<std::vector<double> >("OL_1222");
61  OL1232 = pset.getParameter<std::vector<double> >("OL_1232");
62  OL1213 = pset.getParameter<std::vector<double> >("OL_1213");
63  OL2222 = pset.getParameter<std::vector<double> >("OL_1222");
64 
65  SME11 = pset.getParameter<std::vector<double> >("SME_11");
66  SME12 = pset.getParameter<std::vector<double> >("SME_12");
67  SME13 = pset.getParameter<std::vector<double> >("SME_13");
68  SME21 = pset.getParameter<std::vector<double> >("SME_21");
69  SME22 = pset.getParameter<std::vector<double> >("SME_22");
70  SME31 = pset.getParameter<std::vector<double> >("SME_31");
71  SME32 = pset.getParameter<std::vector<double> >("SME_32");
72  SME41 = pset.getParameter<std::vector<double> >("SME_41");
73 
74  SMB10 = pset.getParameter<std::vector<double> >("SMB_10");
75  SMB11 = pset.getParameter<std::vector<double> >("SMB_11");
76  SMB12 = pset.getParameter<std::vector<double> >("SMB_12");
77  SMB20 = pset.getParameter<std::vector<double> >("SMB_20");
78  SMB21 = pset.getParameter<std::vector<double> >("SMB_21");
79  SMB22 = pset.getParameter<std::vector<double> >("SMB_22");
80  SMB30 = pset.getParameter<std::vector<double> >("SMB_30");
81  SMB31 = pset.getParameter<std::vector<double> >("SMB_31");
82  SMB32 = pset.getParameter<std::vector<double> >("SMB_32");
83 
84  // Load dphi scale parameters
85  CSC01_1 = pset.getParameter<std::vector<double> >("CSC_01_1_scale");
86  CSC12_1 = pset.getParameter<std::vector<double> >("CSC_12_1_scale");
87  CSC12_2 = pset.getParameter<std::vector<double> >("CSC_12_2_scale");
88  CSC12_3 = pset.getParameter<std::vector<double> >("CSC_12_3_scale");
89  CSC13_2 = pset.getParameter<std::vector<double> >("CSC_13_2_scale");
90  CSC13_3 = pset.getParameter<std::vector<double> >("CSC_13_3_scale");
91  CSC14_3 = pset.getParameter<std::vector<double> >("CSC_14_3_scale");
92  CSC23_1 = pset.getParameter<std::vector<double> >("CSC_23_1_scale");
93  CSC23_2 = pset.getParameter<std::vector<double> >("CSC_23_2_scale");
94  CSC24_1 = pset.getParameter<std::vector<double> >("CSC_24_1_scale");
95  CSC34_1 = pset.getParameter<std::vector<double> >("CSC_34_1_scale");
96 
97  DT12_1 = pset.getParameter<std::vector<double> >("DT_12_1_scale");
98  DT12_2 = pset.getParameter<std::vector<double> >("DT_12_2_scale");
99  DT13_1 = pset.getParameter<std::vector<double> >("DT_13_1_scale");
100  DT13_2 = pset.getParameter<std::vector<double> >("DT_13_2_scale");
101  DT14_1 = pset.getParameter<std::vector<double> >("DT_14_1_scale");
102  DT14_2 = pset.getParameter<std::vector<double> >("DT_14_2_scale");
103  DT23_1 = pset.getParameter<std::vector<double> >("DT_23_1_scale");
104  DT23_2 = pset.getParameter<std::vector<double> >("DT_23_2_scale");
105  DT24_1 = pset.getParameter<std::vector<double> >("DT_24_1_scale");
106  DT24_2 = pset.getParameter<std::vector<double> >("DT_24_2_scale");
107  DT34_1 = pset.getParameter<std::vector<double> >("DT_34_1_scale");
108  DT34_2 = pset.getParameter<std::vector<double> >("DT_34_2_scale");
109 
110  OL_1213 = pset.getParameter<std::vector<double> >("OL_1213_0_scale");
111  OL_1222 = pset.getParameter<std::vector<double> >("OL_1222_0_scale");
112  OL_1232 = pset.getParameter<std::vector<double> >("OL_1232_0_scale");
113  OL_2213 = pset.getParameter<std::vector<double> >("OL_2213_0_scale");
114  OL_2222 = pset.getParameter<std::vector<double> >("OL_2222_0_scale");
115 
116  SMB_10S = pset.getParameter<std::vector<double> >("SMB_10_0_scale");
117  SMB_11S = pset.getParameter<std::vector<double> >("SMB_11_0_scale");
118  SMB_12S = pset.getParameter<std::vector<double> >("SMB_12_0_scale");
119  SMB_20S = pset.getParameter<std::vector<double> >("SMB_20_0_scale");
120  SMB_21S = pset.getParameter<std::vector<double> >("SMB_21_0_scale");
121  SMB_22S = pset.getParameter<std::vector<double> >("SMB_22_0_scale");
122  SMB_30S = pset.getParameter<std::vector<double> >("SMB_30_0_scale");
123  SMB_31S = pset.getParameter<std::vector<double> >("SMB_31_0_scale");
124  SMB_32S = pset.getParameter<std::vector<double> >("SMB_32_0_scale");
125 
126  SME_11S = pset.getParameter<std::vector<double> >("SME_11_0_scale");
127  SME_12S = pset.getParameter<std::vector<double> >("SME_12_0_scale");
128  SME_13S = pset.getParameter<std::vector<double> >("SME_13_0_scale");
129  SME_21S = pset.getParameter<std::vector<double> >("SME_21_0_scale");
130  SME_22S = pset.getParameter<std::vector<double> >("SME_22_0_scale");
131 }

References CSC01_1, CSC02, CSC03, CSC12, CSC12_1, CSC12_2, CSC12_3, CSC13, CSC13_2, CSC13_3, CSC14, CSC14_3, CSC23, CSC23_1, CSC23_2, CSC24, CSC24_1, CSC34, CSC34_1, debug, defaultMomentum, DT12, DT12_1, DT12_2, DT13, DT13_1, DT13_2, DT14, DT14_1, DT14_2, DT23, DT23_1, DT23_2, DT24, DT24_1, DT24_2, DT34, DT34_1, DT34_2, OL1213, OL1222, OL1232, OL2222, OL_1213, OL_1222, OL_1232, OL_2213, OL_2222, muonDTDigis_cfi::pset, SMB10, SMB11, SMB12, SMB20, SMB21, SMB22, SMB30, SMB31, SMB32, SMB_10S, SMB_11S, SMB_12S, SMB_20S, SMB_21S, SMB_22S, SMB_30S, SMB_31S, SMB_32S, SME11, SME12, SME13, SME21, SME22, SME31, SME32, SME41, SME_11S, SME_12S, SME_13S, SME_21S, SME_22S, sysError, theMaxMomentum, and theMinMomentum.

◆ ~MuonSeedCreator()

MuonSeedCreator::~MuonSeedCreator ( )

Destructor.

Definition at line 136 of file MuonSeedCreator.cc.

136 {}

Member Function Documentation

◆ createSeed()

TrajectorySeed MuonSeedCreator::createSeed ( int  type,
const SegmentContainer seg,
const std::vector< int > &  layers,
int  NShower,
int  NShowerSeg 
)

Create a seed from set of segments.

get the Global position

get the Global direction

scale the magnitude of total momentum

Trasfer into local direction

get the Global position

get the Global direction

count the energy loss - from parameterization

scale the magnitude of total momentum

Trasfer into local direction

Definition at line 145 of file MuonSeedCreator.cc.

146  {
147  // The index of the station closest to the IP
148  int last = 0;
149 
150  double ptmean = theMinMomentum;
151  double sptmean = theMinMomentum;
152 
153  AlgebraicVector t(4);
154  AlgebraicSymMatrix mat(5, 0);
155  LocalPoint segPos;
156 
157  // Compute the pt according to station types used;
158  if (type == 1)
159  estimatePtCSC(seg, layers, ptmean, sptmean);
160  if (type == 2)
161  estimatePtOverlap(seg, layers, ptmean, sptmean);
162  if (type == 3)
163  estimatePtDT(seg, layers, ptmean, sptmean);
164  if (type == 4)
165  estimatePtSingle(seg, layers, ptmean, sptmean);
166  // type 5 are the seeding for ME1/4
167  if (type == 5)
168  estimatePtCSC(seg, layers, ptmean, sptmean);
169 
170  // for certain clear showering case, set-up the minimum value
171  if (NShowers > 0)
172  estimatePtShowering(NShowers, NShowerSegments, ptmean, sptmean);
173  //if ( NShowers > 0 ) std::cout<<" Showering happened "<<NShowers<<" times w/ "<< NShowerSegments<<std::endl; ;
174 
175  // Minimal pt
176  double charge = 1.0;
177  if (ptmean < 0.)
178  charge = -1.0;
179  if ((charge * ptmean) < theMinMomentum) {
180  ptmean = theMinMomentum * charge;
181  sptmean = theMinMomentum;
182  } else if ((charge * ptmean) > theMaxMomentum) {
183  ptmean = theMaxMomentum * charge;
184  sptmean = theMaxMomentum * 0.25;
185  }
186 
188 
189  double p_err = 0.0;
190  // determine the seed layer
191  int best_seg = 0;
192  double chi2_dof = 9999.0;
193  unsigned int ini_seg = 0;
194  // avoid generating seed from 1st layer(ME1/1)
195  if (type == 5)
196  ini_seg = 1;
197  for (size_t i = ini_seg; i < seg.size(); i++) {
198  double dof = static_cast<double>(seg[i]->degreesOfFreedom());
199  if (chi2_dof < (seg[i]->chi2() / dof))
200  continue;
201  chi2_dof = seg[i]->chi2() / dof;
202  best_seg = static_cast<int>(i);
203  }
204 
205  if (type == 1 || type == 5 || type == 4) {
206  // Fill the LocalTrajectoryParameters
208  last = best_seg;
209  // last = 0;
210  GlobalVector mom = seg[last]->globalPosition() - GlobalPoint();
211  segPos = seg[last]->localPosition();
213 
214  GlobalVector polar(GlobalVector::Spherical(mom.theta(), seg[last]->globalDirection().phi(), 1.));
215 
217  polar *= fabs(ptmean) / polar.perp();
219  LocalVector segDirFromPos = seg[last]->det()->toLocal(polar);
220  int chargeI = static_cast<int>(charge);
221  LocalTrajectoryParameters param1(segPos, segDirFromPos, chargeI);
222  param = param1;
223  p_err = (sptmean * sptmean) / (polar.mag() * polar.mag() * ptmean * ptmean);
224  mat = seg[last]->parametersError().similarityT(seg[last]->projectionMatrix());
225  mat[0][0] = p_err;
226  if (type == 5) {
227  mat[0][0] = mat[0][0] / fabs(tan(mom.theta()));
228  mat[1][1] = mat[1][1] / fabs(tan(mom.theta()));
229  mat[3][3] = 2.25 * mat[3][3];
230  mat[4][4] = 2.25 * mat[4][4];
231  }
232  if (type == 4) {
233  mat[0][0] = mat[0][0] / fabs(tan(mom.theta()));
234  mat[1][1] = mat[1][1] / fabs(tan(mom.theta()));
235  mat[2][2] = 2.25 * mat[2][2];
236  mat[3][3] = 2.25 * mat[3][3];
237  mat[4][4] = 2.25 * mat[4][4];
238  }
239  double dh = fabs(seg[last]->globalPosition().eta()) - 1.6;
240  if (fabs(dh) < 0.1 && type == 1) {
241  mat[1][1] = 4. * mat[1][1];
242  mat[2][2] = 4. * mat[2][2];
243  mat[3][3] = 9. * mat[3][3];
244  mat[4][4] = 9. * mat[4][4];
245  }
246 
247  //if ( !highPt && type != 1 ) mat[1][1]= 2.25*mat[1][1];
248  //if ( highPt && type != 1 ) mat[3][3]= 2.25*mat[1][1];
249  //mat[2][2]= 3.*mat[2][2];
250  //mat[3][3]= 2.*mat[3][3];
251  //mat[4][4]= 2.*mat[4][4];
252  } else {
253  // Fill the LocalTrajectoryParameters
255  last = 0;
256  segPos = seg[last]->localPosition();
257  GlobalVector mom = seg[last]->globalPosition() - GlobalPoint();
259  GlobalVector polar(GlobalVector::Spherical(mom.theta(), seg[last]->globalDirection().phi(), 1.));
260  //GlobalVector polar(GlobalVector::Spherical(seg[last]->globalDirection().theta(),seg[last]->globalDirection().phi(),1.));
261 
263  //double ptRatio = 1.0 - (2.808/(fabs(ptmean) -1)) + (4.546/( (fabs(ptmean)-1)*(fabs(ptmean)-1)) );
264  //ptmean = ptmean*ptRatio ;
265 
267  polar *= fabs(ptmean) / polar.perp();
269  LocalVector segDirFromPos = seg[last]->det()->toLocal(polar);
270  int chargeI = static_cast<int>(charge);
271  LocalTrajectoryParameters param1(segPos, segDirFromPos, chargeI);
272  param = param1;
273  p_err = (sptmean * sptmean) / (polar.mag() * polar.mag() * ptmean * ptmean);
274  mat = seg[last]->parametersError().similarityT(seg[last]->projectionMatrix());
275  //mat[0][0]= 1.44 * p_err;
276  mat[0][0] = p_err;
277  }
278 
279  if (debug) {
280  GlobalPoint gp = seg[last]->globalPosition();
281  float Geta = gp.eta();
282 
283  std::cout << "Type " << type << " Nsegments " << layers.size() << " ";
284  std::cout << "pt " << ptmean << " errpt " << sptmean << " eta " << Geta << " charge " << charge << std::endl;
285  }
286 
287  // this perform H.T() * parErr * H, which is the projection of the
288  // the measurement error (rechit rf) to the state error (TSOS rf)
289  // Legend:
290  // H => is the 4x5 projection matrix
291  // parError the 4x4 parameter error matrix of the RecHit
292 
293  LocalTrajectoryError error(asSMatrix<5>(mat));
294 
295  // Create the TrajectoryStateOnSurface
296  TrajectoryStateOnSurface tsos(param, error, seg[last]->det()->surface(), &*BField);
297 
298  // Take the DetLayer on which relies the segment
299  DetId id = seg[last]->geographicalId();
300 
301  // Transform it in a TrajectoryStateOnSurface
302 
304 
306  for (unsigned l = 0; l < seg.size(); l++) {
307  container.push_back(seg[l]->hit()->clone());
308  //container.push_back(seg[l]->hit());
309  }
310 
311  TrajectorySeed theSeed(seedTSOS, container, alongMomentum);
312  return theSeed;
313 }

References alongMomentum, BField, ALCARECOTkAlJpsiMuMu_cff::charge, hltPixelTracks_cff::chi2, clone(), gather_cfg::cout, debug, cuy::dh, relativeConstraints::error, estimatePtCSC(), estimatePtDT(), estimatePtOverlap(), estimatePtShowering(), estimatePtSingle(), PVValHelper::eta, runTauDisplay::gp, mps_fire::i, cmsLHEtoEOSManager::l, dqmdumpme::last, hgcalTopologyTester_cfi::layers, PV3DBase< T, PVType, FrameType >::perp(), trajectoryStateTransform::persistentState(), edm::OwnVector< T, P >::push_back(), submitPVValidationJobs::t, funct::tan(), theMaxMomentum, theMinMomentum, and PV3DBase< T, PVType, FrameType >::theta().

Referenced by MuonSeedBuilder::build().

◆ estimatePtCSC()

void MuonSeedCreator::estimatePtCSC ( const SegmentContainer seg,
const std::vector< int > &  layers,
double &  pt,
double &  spt 
)
private

Estimate transverse momentum of track from CSC measurements.

Definition at line 324 of file MuonSeedCreator.cc.

327  {
328  unsigned size = seg.size();
329  if (size < 2)
330  return;
331 
332  // reverse the segment and layer container first for pure CSC case
333  //if ( layers[0] > layers[ layers.size()-1 ] ) {
334  // reverse( layers.begin(), layers.end() );
335  // reverse( seg.begin(), seg.end() );
336  //}
337 
338  std::vector<double> ptEstimate;
339  std::vector<double> sptEstimate;
340 
341  thePt = defaultMomentum;
342  theSpt = defaultMomentum;
343 
344  double pt = 0.;
345  double spt = 0.;
346  GlobalPoint segPos[2];
347 
348  int layer0 = layers[0];
349  segPos[0] = seg[0]->globalPosition();
350  float eta = fabs(segPos[0].eta());
351  //float corr = fabs( tan(segPos[0].theta()) );
352  // use pt from vertex information
353  /*
354  if ( layer0 == 0 ) {
355  SegmentContainer seg0;
356  seg0.push_back(seg[0]);
357  std::vector<int> lyr0(1,0);
358  estimatePtSingle( seg0, lyr0, thePt, theSpt);
359  ptEstimate.push_back( thePt );
360  sptEstimate.push_back( theSpt );
361  }
362  */
363 
364  //std::cout<<" estimate CSC "<<std::endl;
365 
366  unsigned idx1 = size;
367  if (size > 1) {
368  while (idx1 > 1) {
369  idx1--;
370  int layer1 = layers[idx1];
371  if (layer0 == layer1)
372  continue;
373  segPos[1] = seg[idx1]->globalPosition();
374 
375  double dphi = segPos[0].phi() - segPos[1].phi();
376  //double temp_dphi = dphi/corr;
377  double temp_dphi = dphi;
378 
379  double sign = 1.0;
380  if (temp_dphi < 0.) {
381  temp_dphi = -1.0 * temp_dphi;
382  sign = -1.0;
383  }
384 
385  // Ensure that delta phi is not too small to prevent pt from blowing up
386  if (temp_dphi < 0.0001) {
387  temp_dphi = 0.0001;
388  pt = theMaxMomentum;
389  spt = theMaxMomentum * 0.25;
390  ptEstimate.push_back(pt * sign);
391  sptEstimate.push_back(spt);
392  }
393  // ME1 is inner-most
394  if (layer0 == 0 && temp_dphi >= 0.0001) {
395  // ME1/2 is outer-most
396  if (layer1 == 1) {
397  //temp_dphi = scaledPhi(temp_dphi, CSC01_1[3] );
398  pt = getPt(CSC01, eta, temp_dphi)[0];
399  spt = getPt(CSC01, eta, temp_dphi)[1];
400  }
401  // ME2 is outer-most
402  else if (layer1 == 2) {
403  //temp_dphi = scaledPhi(temp_dphi, CSC12_3[3] );
404  pt = getPt(CSC02, eta, temp_dphi)[0];
405  spt = getPt(CSC02, eta, temp_dphi)[1];
406  }
407  // ME3 is outer-most
408  else if (layer1 == 3) {
409  //temp_dphi = scaledPhi(temp_dphi, CSC13_3[3] );
410  pt = getPt(CSC03, eta, temp_dphi)[0];
411  spt = getPt(CSC03, eta, temp_dphi)[1];
412  }
413  // ME4 is outer-most
414  else {
415  //temp_dphi = scaledPhi(temp_dphi, CSC14_3[3]);
416  pt = getPt(CSC14, eta, temp_dphi)[0];
417  spt = getPt(CSC14, eta, temp_dphi)[1];
418  }
419  ptEstimate.push_back(pt * sign);
420  sptEstimate.push_back(spt);
421  }
422 
423  // ME1/2,ME1/3 is inner-most
424  if (layer0 == 1) {
425  // ME2 is outer-most
426  if (layer1 == 2) {
427  //if ( eta <= 1.2 ) { temp_dphi = scaledPhi(temp_dphi, CSC12_1[3]); }
428  //if ( eta > 1.2 ) { temp_dphi = scaledPhi(temp_dphi, CSC12_2[3]); }
429  pt = getPt(CSC12, eta, temp_dphi)[0];
430  spt = getPt(CSC12, eta, temp_dphi)[1];
431  }
432  // ME3 is outer-most
433  else if (layer1 == 3) {
434  temp_dphi = scaledPhi(temp_dphi, CSC13_2[3]);
435  pt = getPt(CSC13, eta, temp_dphi)[0];
436  spt = getPt(CSC13, eta, temp_dphi)[1];
437  }
438  // ME4 is outer-most
439  else {
440  temp_dphi = scaledPhi(temp_dphi, CSC14_3[3]);
441  pt = getPt(CSC14, eta, temp_dphi)[0];
442  spt = getPt(CSC14, eta, temp_dphi)[1];
443  }
444  ptEstimate.push_back(pt * sign);
445  sptEstimate.push_back(spt);
446  }
447 
448  // ME2 is inner-most
449  if (layer0 == 2 && temp_dphi > 0.0001) {
450  // ME4 is outer-most
451  if (layer1 == 4) {
452  temp_dphi = scaledPhi(temp_dphi, CSC24_1[3]);
453  pt = getPt(CSC24, eta, temp_dphi)[0];
454  spt = getPt(CSC24, eta, temp_dphi)[1];
455  }
456  // ME3 is outer-most
457  else {
458  // if ME2-4 is availabe , discard ME2-3
459  if (eta <= 1.7) {
460  temp_dphi = scaledPhi(temp_dphi, CSC23_1[3]);
461  }
462  if (eta > 1.7) {
463  temp_dphi = scaledPhi(temp_dphi, CSC23_2[3]);
464  }
465  pt = getPt(CSC23, eta, temp_dphi)[0];
466  spt = getPt(CSC23, eta, temp_dphi)[1];
467  }
468  ptEstimate.push_back(pt * sign);
469  sptEstimate.push_back(spt);
470  }
471 
472  // ME3 is inner-most
473  if (layer0 == 3 && temp_dphi > 0.0001) {
474  temp_dphi = scaledPhi(temp_dphi, CSC34_1[3]);
475  pt = getPt(CSC34, eta, temp_dphi)[0];
476  spt = getPt(CSC34, eta, temp_dphi)[1];
477  ptEstimate.push_back(pt * sign);
478  sptEstimate.push_back(spt);
479  }
480  }
481  }
482 
483  // Compute weighted average if have more than one estimator
484  if (!ptEstimate.empty())
485  weightedPt(ptEstimate, sptEstimate, thePt, theSpt);
486 }

References CSC02, CSC03, CSC12, CSC13, CSC13_2, CSC14, CSC14_3, CSC23, CSC23_1, CSC23_2, CSC24, CSC24_1, CSC34, CSC34_1, defaultMomentum, PVValHelper::eta, getPt(), hgcalTopologyTester_cfi::layers, PV3DBase< T, PVType, FrameType >::phi(), DiDispStaMuonMonitor_cfi::pt, scaledPhi(), Validation_hcalonly_cfi::sign, findQualityFiles::size, theMaxMomentum, and weightedPt().

Referenced by createSeed().

◆ estimatePtDT()

void MuonSeedCreator::estimatePtDT ( const SegmentContainer seg,
const std::vector< int > &  layers,
double &  pt,
double &  spt 
)
private

Estimate transverse momentum of track from DT measurements.

Definition at line 494 of file MuonSeedCreator.cc.

497  {
498  unsigned size = seg.size();
499  if (size < 2)
500  return;
501 
502  std::vector<double> ptEstimate;
503  std::vector<double> sptEstimate;
504 
505  thePt = defaultMomentum;
506  theSpt = defaultMomentum;
507 
508  double pt = 0.;
509  double spt = 0.;
510  GlobalPoint segPos[2];
511 
512  int layer0 = layers[0];
513  segPos[0] = seg[0]->globalPosition();
514  float eta = fabs(segPos[0].eta());
515 
516  //std::cout<<" estimate DT "<<std::endl;
517  // inner-most layer
518  //for ( unsigned idx0 = 0; idx0 < size-1; ++idx0 ) {
519  // layer0 = layers[idx0];
520  // segPos[0] = seg[idx0]->globalPosition();
521  // outer-most layer
522  // for ( unsigned idx1 = idx0+1; idx1 < size; ++idx1 ) {
523  for (unsigned idx1 = 1; idx1 < size; ++idx1) {
524  int layer1 = layers[idx1];
525  segPos[1] = seg[idx1]->globalPosition();
526 
527  //eta = fabs(segPos[1].eta());
528  //if (layer1 == -4) eta = fabs(segPos[0].eta());
529 
530  double dphi = segPos[0].phi() - segPos[1].phi();
531  double temp_dphi = dphi;
532 
533  // Ensure that delta phi is not too small to prevent pt from blowing up
534 
535  double sign = 1.0;
536  if (temp_dphi < 0.) {
537  temp_dphi = -temp_dphi;
538  sign = -1.0;
539  }
540 
541  if (temp_dphi < 0.0001) {
542  temp_dphi = 0.0001;
543  pt = theMaxMomentum;
544  spt = theMaxMomentum * 0.25;
545  ptEstimate.push_back(pt * sign);
546  sptEstimate.push_back(spt);
547  }
548 
549  // MB1 is inner-most
550  if (layer0 == -1 && temp_dphi > 0.0001) {
551  // MB2 is outer-most
552  if (layer1 == -2) {
553  if (eta <= 0.7) {
554  temp_dphi = scaledPhi(temp_dphi, DT12_1[3]);
555  }
556  if (eta > 0.7) {
557  temp_dphi = scaledPhi(temp_dphi, DT12_2[3]);
558  }
559  pt = getPt(DT12, eta, temp_dphi)[0];
560  spt = getPt(DT12, eta, temp_dphi)[1];
561  }
562  // MB3 is outer-most
563  else if (layer1 == -3) {
564  if (eta <= 0.6) {
565  temp_dphi = scaledPhi(temp_dphi, DT13_1[3]);
566  }
567  if (eta > 0.6) {
568  temp_dphi = scaledPhi(temp_dphi, DT13_2[3]);
569  }
570  pt = getPt(DT13, eta, temp_dphi)[0];
571  spt = getPt(DT13, eta, temp_dphi)[1];
572  }
573  // MB4 is outer-most
574  else {
575  if (eta <= 0.52) {
576  temp_dphi = scaledPhi(temp_dphi, DT14_1[3]);
577  }
578  if (eta > 0.52) {
579  temp_dphi = scaledPhi(temp_dphi, DT14_2[3]);
580  }
581  pt = getPt(DT14, eta, temp_dphi)[0];
582  spt = getPt(DT14, eta, temp_dphi)[1];
583  }
584  ptEstimate.push_back(pt * sign);
585  sptEstimate.push_back(spt);
586  }
587 
588  // MB2 is inner-most
589  if (layer0 == -2 && temp_dphi > 0.0001) {
590  // MB3 is outer-most
591  if (layer1 == -3) {
592  if (eta <= 0.6) {
593  temp_dphi = scaledPhi(temp_dphi, DT23_1[3]);
594  }
595  if (eta > 0.6) {
596  temp_dphi = scaledPhi(temp_dphi, DT23_2[3]);
597  }
598  pt = getPt(DT23, eta, temp_dphi)[0];
599  spt = getPt(DT23, eta, temp_dphi)[1];
600  }
601  // MB4 is outer-most
602  else {
603  if (eta <= 0.52) {
604  temp_dphi = scaledPhi(temp_dphi, DT24_1[3]);
605  }
606  if (eta > 0.52) {
607  temp_dphi = scaledPhi(temp_dphi, DT24_2[3]);
608  }
609  pt = getPt(DT24, eta, temp_dphi)[0];
610  spt = getPt(DT24, eta, temp_dphi)[1];
611  }
612  ptEstimate.push_back(pt * sign);
613  sptEstimate.push_back(spt);
614  }
615 
616  // MB3 is inner-most -> only marginally useful to pick up the charge
617  if (layer0 == -3 && temp_dphi > 0.0001) {
618  // MB4 is outer-most
619 
620  if (eta <= 0.51) {
621  temp_dphi = scaledPhi(temp_dphi, DT34_1[3]);
622  }
623  if (eta > 0.51) {
624  temp_dphi = scaledPhi(temp_dphi, DT34_2[3]);
625  }
626  pt = getPt(DT34, eta, temp_dphi)[0];
627  spt = getPt(DT34, eta, temp_dphi)[1];
628  ptEstimate.push_back(pt * sign);
629  sptEstimate.push_back(spt);
630  }
631  }
632  //}
633 
634  // Compute weighted average if have more than one estimator
635  if (!ptEstimate.empty())
636  weightedPt(ptEstimate, sptEstimate, thePt, theSpt);
637 }

References defaultMomentum, DT12, DT12_1, DT12_2, DT13, DT13_1, DT13_2, DT14, DT14_1, DT14_2, DT23, DT23_1, DT23_2, DT24, DT24_1, DT24_2, DT34, DT34_1, DT34_2, PVValHelper::eta, getPt(), hgcalTopologyTester_cfi::layers, PV3DBase< T, PVType, FrameType >::phi(), DiDispStaMuonMonitor_cfi::pt, scaledPhi(), Validation_hcalonly_cfi::sign, findQualityFiles::size, theMaxMomentum, and weightedPt().

Referenced by createSeed(), and estimatePtOverlap().

◆ estimatePtOverlap()

void MuonSeedCreator::estimatePtOverlap ( const SegmentContainer seg,
const std::vector< int > &  layers,
double &  pt,
double &  spt 
)
private

Estimate transverse momentum of track from CSC + DT measurements.

Definition at line 643 of file MuonSeedCreator.cc.

646  {
647  int size = layers.size();
648 
649  thePt = defaultMomentum;
650  theSpt = defaultMomentum;
651 
652  SegmentContainer segCSC;
653  std::vector<int> layersCSC;
654  SegmentContainer segDT;
655  std::vector<int> layersDT;
656 
657  // DT layers are numbered as -4 to -1, whereas CSC layers go from 0 to 4:
658  for (unsigned j = 0; j < layers.size(); ++j) {
659  if (layers[j] > -1) {
660  segCSC.push_back(seg[j]);
661  layersCSC.push_back(layers[j]);
662  } else {
663  segDT.push_back(seg[j]);
664  layersDT.push_back(layers[j]);
665  }
666  }
667 
668  std::vector<double> ptEstimate;
669  std::vector<double> sptEstimate;
670 
671  GlobalPoint segPos[2];
672  int layer0 = layers[0];
673  segPos[0] = seg[0]->globalPosition();
674  float eta = fabs(segPos[0].eta());
675  //std::cout<<" estimate OL "<<std::endl;
676 
677  if (!segDT.empty() && !segCSC.empty()) {
678  int layer1 = layers[size - 1];
679  segPos[1] = seg[size - 1]->globalPosition();
680 
681  double dphi = segPos[0].phi() - segPos[1].phi();
682  double temp_dphi = dphi;
683 
684  // Ensure that delta phi is not too small to prevent pt from blowing up
685 
686  double sign = 1.0;
687  if (temp_dphi < 0.) {
688  temp_dphi = -temp_dphi;
689  sign = -1.0;
690  }
691 
692  if (temp_dphi < 0.0001) {
693  temp_dphi = 0.0001;
694  thePt = theMaxMomentum;
695  theSpt = theMaxMomentum * 0.25;
696  ptEstimate.push_back(thePt * sign);
697  sptEstimate.push_back(theSpt);
698  }
699 
700  // MB1 is inner-most
701  if (layer0 == -1 && temp_dphi > 0.0001) {
702  // ME1/3 is outer-most
703  if (layer1 == 1) {
704  temp_dphi = scaledPhi(temp_dphi, OL_1213[3]);
705  thePt = getPt(OL1213, eta, temp_dphi)[0];
706  theSpt = getPt(OL1213, eta, temp_dphi)[1];
707  }
708  // ME2 is outer-most
709  else if (layer1 == 2) {
710  temp_dphi = scaledPhi(temp_dphi, OL_1222[3]);
711  thePt = getPt(OL1222, eta, temp_dphi)[0];
712  theSpt = getPt(OL1222, eta, temp_dphi)[1];
713  }
714  // ME3 is outer-most
715  else {
716  temp_dphi = scaledPhi(temp_dphi, OL_1232[3]);
717  thePt = getPt(OL1232, eta, temp_dphi)[0];
718  theSpt = getPt(OL1232, eta, temp_dphi)[1];
719  }
720  ptEstimate.push_back(thePt * sign);
721  sptEstimate.push_back(theSpt);
722  }
723  // MB2 is inner-most
724  if (layer0 == -2 && temp_dphi > 0.0001) {
725  // ME1/3 is outer-most
726  if (layer1 == 1) {
727  temp_dphi = scaledPhi(temp_dphi, OL_2213[3]);
728  thePt = getPt(OL2213, eta, temp_dphi)[0];
729  theSpt = getPt(OL2213, eta, temp_dphi)[1];
730  ptEstimate.push_back(thePt * sign);
731  sptEstimate.push_back(theSpt);
732  }
733  // ME2 is outer-most
734  if (layer1 == 2) {
735  temp_dphi = scaledPhi(temp_dphi, OL_2222[3]);
736  thePt = getPt(OL2222, eta, temp_dphi)[0];
737  theSpt = getPt(OL2222, eta, temp_dphi)[1];
738  }
739  }
740  }
741 
742  if (segDT.size() > 1) {
743  estimatePtDT(segDT, layersDT, thePt, theSpt);
744  ptEstimate.push_back(thePt);
745  sptEstimate.push_back(theSpt);
746  }
747 
748  /*
749  // not useful ....and pt estimation is bad
750  if ( segCSC.size() > 1 ) {
751  // don't estimate pt without ME1 information
752  bool CSCLayer1=false;
753  for (unsigned i=0; i< layersCSC.size(); i++) {
754  if ( layersCSC[i]==0 || layersCSC[i]==1 ) CSCLayer1 = true;
755  }
756  if (CSCLayer1) {
757  estimatePtCSC(segCSC, layersCSC, thePt, theSpt);
758  ptEstimate.push_back(thePt);
759  sptEstimate.push_back(theSpt);
760  }
761  }
762  */
763 
764  // Compute weighted average if have more than one estimator
765  if (!ptEstimate.empty())
766  weightedPt(ptEstimate, sptEstimate, thePt, theSpt);
767 }

References defaultMomentum, estimatePtDT(), PVValHelper::eta, getPt(), dqmiolumiharvest::j, hgcalTopologyTester_cfi::layers, OL1213, OL1222, OL1232, OL2213, OL2222, OL_1213, OL_1222, OL_1232, OL_2213, OL_2222, PV3DBase< T, PVType, FrameType >::phi(), scaledPhi(), Validation_hcalonly_cfi::sign, findQualityFiles::size, theMaxMomentum, and weightedPt().

Referenced by createSeed().

◆ estimatePtShowering()

void MuonSeedCreator::estimatePtShowering ( int &  NShowers,
int &  NShowerSeg,
double &  pt,
double &  spt 
)
private

Estimate transverse momentum of track from showering segment.

Definition at line 909 of file MuonSeedCreator.cc.

909  {
910  if (NShowers > 2 && thePt < 300.) {
911  thePt = 800.;
912  theSpt = 200.;
913  }
914  if (NShowers == 2 && NShowerSegments > 11 && thePt < 150.) {
915  thePt = 280.;
916  theSpt = 70.;
917  }
918  if (NShowers == 2 && NShowerSegments <= 11 && thePt < 50.) {
919  thePt = 80.;
920  theSpt = 40.;
921  }
922  if (NShowers == 1 && NShowerSegments <= 5 && thePt < 10.) {
923  thePt = 16.;
924  theSpt = 8.;
925  }
926 }

Referenced by createSeed().

◆ estimatePtSingle()

void MuonSeedCreator::estimatePtSingle ( const SegmentContainer seg,
const std::vector< int > &  layers,
double &  pt,
double &  spt 
)
private

Estimate transverse momentum of track from single segment.

Definition at line 773 of file MuonSeedCreator.cc.

776  {
777  thePt = defaultMomentum;
778  theSpt = defaultMomentum;
779 
780  GlobalPoint segPos = seg[0]->globalPosition();
781  double eta = segPos.eta();
782  GlobalVector gv = seg[0]->globalDirection();
783 
784  // Psi is angle between the segment origin and segment direction
785  // Use dot product between two vectors to get Psi in global x-y plane
786  double cosDpsi = (gv.x() * segPos.x() + gv.y() * segPos.y());
787  cosDpsi /= sqrt(segPos.x() * segPos.x() + segPos.y() * segPos.y());
788  cosDpsi /= sqrt(gv.x() * gv.x() + gv.y() * gv.y());
789 
790  double axb = (segPos.x() * gv.y()) - (segPos.y() * gv.x());
791  double sign = (axb < 0.) ? 1.0 : -1.0;
792 
793  double dpsi = fabs(acos(cosDpsi));
794  if (dpsi > 1.570796) {
795  dpsi = 3.141592 - dpsi;
796  sign = -1. * sign;
797  }
798  if (fabs(dpsi) < 0.00005) {
799  dpsi = 0.00005;
800  }
801 
802  // the 1st layer
803  if (layers[0] == -1) {
804  // MB10
805  if (fabs(eta) < 0.3) {
806  dpsi = scaledPhi(dpsi, SMB_10S[3]);
807  thePt = getPt(SMB10, eta, dpsi)[0];
808  theSpt = getPt(SMB10, eta, dpsi)[1];
809  }
810  // MB11
811  if (fabs(eta) >= 0.3 && fabs(eta) < 0.82) {
812  dpsi = scaledPhi(dpsi, SMB_11S[3]);
813  thePt = getPt(SMB11, eta, dpsi)[0];
814  theSpt = getPt(SMB11, eta, dpsi)[1];
815  }
816  // MB12
817  if (fabs(eta) >= 0.82 && fabs(eta) < 1.2) {
818  dpsi = scaledPhi(dpsi, SMB_12S[3]);
819  thePt = getPt(SMB12, eta, dpsi)[0];
820  theSpt = getPt(SMB12, eta, dpsi)[1];
821  }
822  }
823  if (layers[0] == 1) {
824  // ME13
825  if (fabs(eta) > 0.92 && fabs(eta) < 1.16) {
826  dpsi = scaledPhi(dpsi, SME_13S[3]);
827  thePt = getPt(SME13, eta, dpsi)[0];
828  theSpt = getPt(SME13, eta, dpsi)[1];
829  }
830  // ME12
831  if (fabs(eta) >= 1.16 && fabs(eta) <= 1.6) {
832  dpsi = scaledPhi(dpsi, SME_12S[3]);
833  thePt = getPt(SME12, eta, dpsi)[0];
834  theSpt = getPt(SME12, eta, dpsi)[1];
835  }
836  }
837  if (layers[0] == 0) {
838  // ME11
839  if (fabs(eta) > 1.6) {
840  dpsi = scaledPhi(dpsi, SMB_11S[3]);
841  thePt = getPt(SME11, eta, dpsi)[0];
842  theSpt = getPt(SME11, eta, dpsi)[1];
843  }
844  }
845  // the 2nd layer
846  if (layers[0] == -2) {
847  // MB20
848  if (fabs(eta) < 0.25) {
849  dpsi = scaledPhi(dpsi, SMB_20S[3]);
850  thePt = getPt(SMB20, eta, dpsi)[0];
851  theSpt = getPt(SMB20, eta, dpsi)[1];
852  }
853  // MB21
854  if (fabs(eta) >= 0.25 && fabs(eta) < 0.72) {
855  dpsi = scaledPhi(dpsi, SMB_21S[3]);
856  thePt = getPt(SMB21, eta, dpsi)[0];
857  theSpt = getPt(SMB21, eta, dpsi)[1];
858  }
859  // MB22
860  if (fabs(eta) >= 0.72 && fabs(eta) < 1.04) {
861  dpsi = scaledPhi(dpsi, SMB_22S[3]);
862  thePt = getPt(SMB22, eta, dpsi)[0];
863  theSpt = getPt(SMB22, eta, dpsi)[1];
864  }
865  }
866  if (layers[0] == 2) {
867  // ME22
868  if (fabs(eta) > 0.95 && fabs(eta) <= 1.6) {
869  dpsi = scaledPhi(dpsi, SME_22S[3]);
870  thePt = getPt(SME22, eta, dpsi)[0];
871  theSpt = getPt(SME22, eta, dpsi)[1];
872  }
873  // ME21
874  if (fabs(eta) > 1.6 && fabs(eta) < 2.45) {
875  dpsi = scaledPhi(dpsi, SME_21S[3]);
876  thePt = getPt(SME21, eta, dpsi)[0];
877  theSpt = getPt(SME21, eta, dpsi)[1];
878  }
879  }
880 
881  // the 3rd layer
882  if (layers[0] == -3) {
883  // MB30
884  if (fabs(eta) <= 0.22) {
885  dpsi = scaledPhi(dpsi, SMB_30S[3]);
886  thePt = getPt(SMB30, eta, dpsi)[0];
887  theSpt = getPt(SMB30, eta, dpsi)[1];
888  }
889  // MB31
890  if (fabs(eta) > 0.22 && fabs(eta) <= 0.6) {
891  dpsi = scaledPhi(dpsi, SMB_31S[3]);
892  thePt = getPt(SMB31, eta, dpsi)[0];
893  theSpt = getPt(SMB31, eta, dpsi)[1];
894  }
895  // MB32
896  if (fabs(eta) > 0.6 && fabs(eta) < 0.95) {
897  dpsi = scaledPhi(dpsi, SMB_32S[3]);
898  thePt = getPt(SMB32, eta, dpsi)[0];
899  theSpt = getPt(SMB32, eta, dpsi)[1];
900  }
901  }
902  thePt = fabs(thePt) * sign;
903  theSpt = fabs(theSpt);
904 
905  return;
906 }

References defaultMomentum, PVValHelper::eta, PV3DBase< T, PVType, FrameType >::eta(), getPt(), hgcalTopologyTester_cfi::layers, scaledPhi(), Validation_hcalonly_cfi::sign, SMB10, SMB11, SMB12, SMB20, SMB21, SMB22, SMB30, SMB31, SMB32, SMB_10S, SMB_11S, SMB_12S, SMB_20S, SMB_21S, SMB_22S, SMB_30S, SMB_31S, SMB_32S, SME11, SME12, SME13, SME21, SME22, SME_12S, SME_13S, SME_21S, SME_22S, mathSSE::sqrt(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by createSeed().

◆ getPt()

std::vector< double > MuonSeedCreator::getPt ( const std::vector< double > &  vParameters,
double  eta,
double  dPhi 
)
private

Compute pt from parameters.

Definition at line 1001 of file MuonSeedCreator.cc.

1001  {
1002  double h = fabs(eta);
1003  double estPt = (vPara[0] + vPara[1] * h + vPara[2] * h * h) / dPhi;
1004  double estSPt = (vPara[3] + vPara[4] * h + vPara[5] * h * h) * estPt;
1005  std::vector<double> paraPt;
1006  paraPt.push_back(estPt);
1007  paraPt.push_back(estSPt);
1008 
1009  //std::cout<<" pt:"<<estPt<<" +/-"<< estSPt<<" h:"<<eta<<" df:"<<dPhi<<std::endl;
1010  return paraPt;
1011 }

References HLT_FULL_cff::dPhi, PVValHelper::eta, and h.

Referenced by estimatePtCSC(), estimatePtDT(), estimatePtOverlap(), and estimatePtSingle().

◆ scaledPhi()

double MuonSeedCreator::scaledPhi ( double  dphi,
double  t1 
)
private

Scale the dPhi from segment position.

Definition at line 1013 of file MuonSeedCreator.cc.

1013  {
1014  if (dphi != 0.) {
1015  double oPhi = 1. / dphi;
1016  dphi = dphi / (1. + t1 / (oPhi + 10.));
1017  return dphi;
1018 
1019  } else {
1020  return dphi;
1021  }
1022 }

References RandomServiceHelper::t1.

Referenced by estimatePtCSC(), estimatePtDT(), estimatePtOverlap(), and estimatePtSingle().

◆ setBField()

void MuonSeedCreator::setBField ( const MagneticField theField)
inline

Cache Magnetic Field for current event.

Definition at line 43 of file MuonSeedCreator.h.

43 { BField = theField; };

References BField.

Referenced by MuonSeedBuilder::build().

◆ weightedPt()

void MuonSeedCreator::weightedPt ( const std::vector< double > &  ptEstimate,
const std::vector< double > &  sptEstimate,
double &  ptAvg,
double &  sptAvg 
)
private

Compute weighted mean pt from different pt estimators.

Definition at line 934 of file MuonSeedCreator.cc.

937  {
938  int size = ptEstimate.size();
939 
940  // If only one element, by-pass computation below
941  if (size < 2) {
942  thePt = ptEstimate[0];
943  theSpt = sptEstimate[0];
944  return;
945  }
946 
947  double charge = 0.;
948  // If have more than one pt estimator, first look if any estimator is biased
949  // by comparing the charge of each estimator
950 
951  for (unsigned j = 0; j < ptEstimate.size(); j++) {
952  //std::cout<<" weighting pt: "<< ptEstimate[j] <<std::endl;
953  if (ptEstimate[j] < 0.) {
954  // To prevent from blowing up, add 0.1
955  // weight by relative error on pt
956  charge -= 1. * (ptEstimate[j] * ptEstimate[j]) / (sptEstimate[j] * sptEstimate[j]);
957  } else {
958  // weight by relative error on pt
959  charge += 1. * (ptEstimate[j] * ptEstimate[j]) / (sptEstimate[j] * sptEstimate[j]);
960  }
961  }
962 
963  // No need to normalize as we want to know only sign ( + or - )
964  if (charge < 0.) {
965  charge = -1.;
966  } else {
967  charge = 1.;
968  }
969 
970  //int n = 0;
971  double weightPtSum = 0.;
972  double sigmaSqr_sum = 0.;
973 
974  // Now, we want to compute average Pt using estimators with "correct" charge
975  // This is to remove biases
976  for (unsigned j = 0; j < ptEstimate.size(); ++j) {
977  //if ( (minpt_ratio < 0.5) && (fabs(ptEstimate[j]) < 5.0) ) continue;
978  //if ( ptEstimate[j] * charge > 0. ) {
979  //n++;
980  sigmaSqr_sum += 1.0 / (sptEstimate[j] * sptEstimate[j]);
981  weightPtSum += fabs(ptEstimate[j]) / (sptEstimate[j] * sptEstimate[j]);
982  //}
983  }
984  /*
985  if (n < 1) {
986  thePt = defaultMomentum*charge;
987  theSpt = defaultMomentum;
988  return;
989  }
990  */
991  // Compute weighted mean and error
992 
993  thePt = (charge * weightPtSum) / sigmaSqr_sum;
994  theSpt = sqrt(1.0 / sigmaSqr_sum);
995 
996  //std::cout<<" final weighting : "<< thePt <<" ~ "<< fabs( theSpt/thePt ) <<std::endl;
997 
998  return;
999 }

References ALCARECOTkAlJpsiMuMu_cff::charge, dqmiolumiharvest::j, findQualityFiles::size, and mathSSE::sqrt().

Referenced by estimatePtCSC(), estimatePtDT(), and estimatePtOverlap().

Member Data Documentation

◆ BField

const MagneticField* MuonSeedCreator::BField
private

Definition at line 89 of file MuonSeedCreator.h.

Referenced by createSeed(), and setBField().

◆ CSC01

std::vector<double> MuonSeedCreator::CSC01
private

Definition at line 99 of file MuonSeedCreator.h.

◆ CSC01_1

std::vector<double> MuonSeedCreator::CSC01_1
private

Definition at line 136 of file MuonSeedCreator.h.

Referenced by MuonSeedCreator().

◆ CSC02

std::vector<double> MuonSeedCreator::CSC02
private

Definition at line 101 of file MuonSeedCreator.h.

Referenced by estimatePtCSC(), and MuonSeedCreator().

◆ CSC03

std::vector<double> MuonSeedCreator::CSC03
private

Definition at line 103 of file MuonSeedCreator.h.

Referenced by estimatePtCSC(), and MuonSeedCreator().

◆ CSC12

std::vector<double> MuonSeedCreator::CSC12
private

Definition at line 100 of file MuonSeedCreator.h.

Referenced by estimatePtCSC(), and MuonSeedCreator().

◆ CSC12_1

std::vector<double> MuonSeedCreator::CSC12_1
private

Definition at line 137 of file MuonSeedCreator.h.

Referenced by MuonSeedCreator().

◆ CSC12_2

std::vector<double> MuonSeedCreator::CSC12_2
private

Definition at line 138 of file MuonSeedCreator.h.

Referenced by MuonSeedCreator().

◆ CSC12_3

std::vector<double> MuonSeedCreator::CSC12_3
private

Definition at line 139 of file MuonSeedCreator.h.

Referenced by MuonSeedCreator().

◆ CSC13

std::vector<double> MuonSeedCreator::CSC13
private

Definition at line 102 of file MuonSeedCreator.h.

Referenced by estimatePtCSC(), and MuonSeedCreator().

◆ CSC13_2

std::vector<double> MuonSeedCreator::CSC13_2
private

Definition at line 140 of file MuonSeedCreator.h.

Referenced by estimatePtCSC(), and MuonSeedCreator().

◆ CSC13_3

std::vector<double> MuonSeedCreator::CSC13_3
private

Definition at line 141 of file MuonSeedCreator.h.

Referenced by MuonSeedCreator().

◆ CSC14

std::vector<double> MuonSeedCreator::CSC14
private

Definition at line 104 of file MuonSeedCreator.h.

Referenced by estimatePtCSC(), and MuonSeedCreator().

◆ CSC14_3

std::vector<double> MuonSeedCreator::CSC14_3
private

Definition at line 142 of file MuonSeedCreator.h.

Referenced by estimatePtCSC(), and MuonSeedCreator().

◆ CSC23

std::vector<double> MuonSeedCreator::CSC23
private

Definition at line 105 of file MuonSeedCreator.h.

Referenced by estimatePtCSC(), and MuonSeedCreator().

◆ CSC23_1

std::vector<double> MuonSeedCreator::CSC23_1
private

Definition at line 143 of file MuonSeedCreator.h.

Referenced by estimatePtCSC(), and MuonSeedCreator().

◆ CSC23_2

std::vector<double> MuonSeedCreator::CSC23_2
private

Definition at line 144 of file MuonSeedCreator.h.

Referenced by estimatePtCSC(), and MuonSeedCreator().

◆ CSC24

std::vector<double> MuonSeedCreator::CSC24
private

Definition at line 106 of file MuonSeedCreator.h.

Referenced by estimatePtCSC(), and MuonSeedCreator().

◆ CSC24_1

std::vector<double> MuonSeedCreator::CSC24_1
private

Definition at line 145 of file MuonSeedCreator.h.

Referenced by estimatePtCSC(), and MuonSeedCreator().

◆ CSC34

std::vector<double> MuonSeedCreator::CSC34
private

Definition at line 107 of file MuonSeedCreator.h.

Referenced by estimatePtCSC(), and MuonSeedCreator().

◆ CSC34_1

std::vector<double> MuonSeedCreator::CSC34_1
private

Definition at line 146 of file MuonSeedCreator.h.

Referenced by estimatePtCSC(), and MuonSeedCreator().

◆ debug

bool MuonSeedCreator::debug
private

◆ defaultMomentum

float MuonSeedCreator::defaultMomentum
private

◆ DT12

std::vector<double> MuonSeedCreator::DT12
private

Definition at line 92 of file MuonSeedCreator.h.

Referenced by estimatePtDT(), and MuonSeedCreator().

◆ DT12_1

std::vector<double> MuonSeedCreator::DT12_1
private

Definition at line 148 of file MuonSeedCreator.h.

Referenced by estimatePtDT(), and MuonSeedCreator().

◆ DT12_2

std::vector<double> MuonSeedCreator::DT12_2
private

Definition at line 149 of file MuonSeedCreator.h.

Referenced by estimatePtDT(), and MuonSeedCreator().

◆ DT13

std::vector<double> MuonSeedCreator::DT13
private

Definition at line 93 of file MuonSeedCreator.h.

Referenced by estimatePtDT(), and MuonSeedCreator().

◆ DT13_1

std::vector<double> MuonSeedCreator::DT13_1
private

Definition at line 150 of file MuonSeedCreator.h.

Referenced by estimatePtDT(), and MuonSeedCreator().

◆ DT13_2

std::vector<double> MuonSeedCreator::DT13_2
private

Definition at line 151 of file MuonSeedCreator.h.

Referenced by estimatePtDT(), and MuonSeedCreator().

◆ DT14

std::vector<double> MuonSeedCreator::DT14
private

Definition at line 94 of file MuonSeedCreator.h.

Referenced by estimatePtDT(), and MuonSeedCreator().

◆ DT14_1

std::vector<double> MuonSeedCreator::DT14_1
private

Definition at line 152 of file MuonSeedCreator.h.

Referenced by estimatePtDT(), and MuonSeedCreator().

◆ DT14_2

std::vector<double> MuonSeedCreator::DT14_2
private

Definition at line 153 of file MuonSeedCreator.h.

Referenced by estimatePtDT(), and MuonSeedCreator().

◆ DT23

std::vector<double> MuonSeedCreator::DT23
private

Definition at line 95 of file MuonSeedCreator.h.

Referenced by estimatePtDT(), and MuonSeedCreator().

◆ DT23_1

std::vector<double> MuonSeedCreator::DT23_1
private

Definition at line 154 of file MuonSeedCreator.h.

Referenced by estimatePtDT(), and MuonSeedCreator().

◆ DT23_2

std::vector<double> MuonSeedCreator::DT23_2
private

Definition at line 155 of file MuonSeedCreator.h.

Referenced by estimatePtDT(), and MuonSeedCreator().

◆ DT24

std::vector<double> MuonSeedCreator::DT24
private

Definition at line 96 of file MuonSeedCreator.h.

Referenced by estimatePtDT(), and MuonSeedCreator().

◆ DT24_1

std::vector<double> MuonSeedCreator::DT24_1
private

Definition at line 156 of file MuonSeedCreator.h.

Referenced by estimatePtDT(), and MuonSeedCreator().

◆ DT24_2

std::vector<double> MuonSeedCreator::DT24_2
private

Definition at line 157 of file MuonSeedCreator.h.

Referenced by estimatePtDT(), and MuonSeedCreator().

◆ DT34

std::vector<double> MuonSeedCreator::DT34
private

Definition at line 97 of file MuonSeedCreator.h.

Referenced by estimatePtDT(), and MuonSeedCreator().

◆ DT34_1

std::vector<double> MuonSeedCreator::DT34_1
private

Definition at line 158 of file MuonSeedCreator.h.

Referenced by estimatePtDT(), and MuonSeedCreator().

◆ DT34_2

std::vector<double> MuonSeedCreator::DT34_2
private

Definition at line 159 of file MuonSeedCreator.h.

Referenced by estimatePtDT(), and MuonSeedCreator().

◆ OL1213

std::vector<double> MuonSeedCreator::OL1213
private

Definition at line 109 of file MuonSeedCreator.h.

Referenced by estimatePtOverlap(), and MuonSeedCreator().

◆ OL1222

std::vector<double> MuonSeedCreator::OL1222
private

Definition at line 110 of file MuonSeedCreator.h.

Referenced by estimatePtOverlap(), and MuonSeedCreator().

◆ OL1232

std::vector<double> MuonSeedCreator::OL1232
private

Definition at line 111 of file MuonSeedCreator.h.

Referenced by estimatePtOverlap(), and MuonSeedCreator().

◆ OL2213

std::vector<double> MuonSeedCreator::OL2213
private

Definition at line 112 of file MuonSeedCreator.h.

Referenced by estimatePtOverlap().

◆ OL2222

std::vector<double> MuonSeedCreator::OL2222
private

Definition at line 113 of file MuonSeedCreator.h.

Referenced by estimatePtOverlap(), and MuonSeedCreator().

◆ OL_1213

std::vector<double> MuonSeedCreator::OL_1213
private

Definition at line 161 of file MuonSeedCreator.h.

Referenced by estimatePtOverlap(), and MuonSeedCreator().

◆ OL_1222

std::vector<double> MuonSeedCreator::OL_1222
private

Definition at line 162 of file MuonSeedCreator.h.

Referenced by estimatePtOverlap(), and MuonSeedCreator().

◆ OL_1232

std::vector<double> MuonSeedCreator::OL_1232
private

Definition at line 163 of file MuonSeedCreator.h.

Referenced by estimatePtOverlap(), and MuonSeedCreator().

◆ OL_2213

std::vector<double> MuonSeedCreator::OL_2213
private

Definition at line 164 of file MuonSeedCreator.h.

Referenced by estimatePtOverlap(), and MuonSeedCreator().

◆ OL_2222

std::vector<double> MuonSeedCreator::OL_2222
private

Definition at line 165 of file MuonSeedCreator.h.

Referenced by estimatePtOverlap(), and MuonSeedCreator().

◆ SMB10

std::vector<double> MuonSeedCreator::SMB10
private

Definition at line 124 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SMB11

std::vector<double> MuonSeedCreator::SMB11
private

Definition at line 125 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SMB12

std::vector<double> MuonSeedCreator::SMB12
private

Definition at line 126 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SMB20

std::vector<double> MuonSeedCreator::SMB20
private

Definition at line 127 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SMB21

std::vector<double> MuonSeedCreator::SMB21
private

Definition at line 128 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SMB22

std::vector<double> MuonSeedCreator::SMB22
private

Definition at line 129 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SMB30

std::vector<double> MuonSeedCreator::SMB30
private

Definition at line 130 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SMB31

std::vector<double> MuonSeedCreator::SMB31
private

Definition at line 131 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SMB32

std::vector<double> MuonSeedCreator::SMB32
private

Definition at line 132 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SMB_10S

std::vector<double> MuonSeedCreator::SMB_10S
private

Definition at line 167 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SMB_11S

std::vector<double> MuonSeedCreator::SMB_11S
private

Definition at line 168 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SMB_12S

std::vector<double> MuonSeedCreator::SMB_12S
private

Definition at line 169 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SMB_20S

std::vector<double> MuonSeedCreator::SMB_20S
private

Definition at line 170 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SMB_21S

std::vector<double> MuonSeedCreator::SMB_21S
private

Definition at line 171 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SMB_22S

std::vector<double> MuonSeedCreator::SMB_22S
private

Definition at line 172 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SMB_30S

std::vector<double> MuonSeedCreator::SMB_30S
private

Definition at line 173 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SMB_31S

std::vector<double> MuonSeedCreator::SMB_31S
private

Definition at line 174 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SMB_32S

std::vector<double> MuonSeedCreator::SMB_32S
private

Definition at line 175 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SME11

std::vector<double> MuonSeedCreator::SME11
private

Definition at line 115 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SME12

std::vector<double> MuonSeedCreator::SME12
private

Definition at line 116 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SME13

std::vector<double> MuonSeedCreator::SME13
private

Definition at line 117 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SME21

std::vector<double> MuonSeedCreator::SME21
private

Definition at line 118 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SME22

std::vector<double> MuonSeedCreator::SME22
private

Definition at line 119 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SME31

std::vector<double> MuonSeedCreator::SME31
private

Definition at line 120 of file MuonSeedCreator.h.

Referenced by MuonSeedCreator().

◆ SME32

std::vector<double> MuonSeedCreator::SME32
private

Definition at line 121 of file MuonSeedCreator.h.

Referenced by MuonSeedCreator().

◆ SME41

std::vector<double> MuonSeedCreator::SME41
private

Definition at line 122 of file MuonSeedCreator.h.

Referenced by MuonSeedCreator().

◆ SME_11S

std::vector<double> MuonSeedCreator::SME_11S
private

Definition at line 177 of file MuonSeedCreator.h.

Referenced by MuonSeedCreator().

◆ SME_12S

std::vector<double> MuonSeedCreator::SME_12S
private

Definition at line 178 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SME_13S

std::vector<double> MuonSeedCreator::SME_13S
private

Definition at line 179 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SME_21S

std::vector<double> MuonSeedCreator::SME_21S
private

Definition at line 180 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ SME_22S

std::vector<double> MuonSeedCreator::SME_22S
private

Definition at line 181 of file MuonSeedCreator.h.

Referenced by estimatePtSingle(), and MuonSeedCreator().

◆ sysError

double MuonSeedCreator::sysError
private

Definition at line 83 of file MuonSeedCreator.h.

Referenced by MuonSeedCreator().

◆ theMaxMomentum

float MuonSeedCreator::theMaxMomentum
private

◆ theMinMomentum

float MuonSeedCreator::theMinMomentum
private

Definition at line 78 of file MuonSeedCreator.h.

Referenced by createSeed(), and MuonSeedCreator().

Vector3DBase
Definition: Vector3DBase.h:8
MuonSeedCreator::SMB_32S
std::vector< double > SMB_32S
Definition: MuonSeedCreator.h:175
DOFs::dof
dof
Definition: AlignPCLThresholdsWriter.cc:37
mps_fire.i
i
Definition: mps_fire.py:428
MuonSeedCreator::CSC34_1
std::vector< double > CSC34_1
Definition: MuonSeedCreator.h:146
MuonSeedCreator::CSC23
std::vector< double > CSC23
Definition: MuonSeedCreator.h:105
MuonSeedCreator::DT12
std::vector< double > DT12
Definition: MuonSeedCreator.h:92
MuonSeedCreator::SME_22S
std::vector< double > SME_22S
Definition: MuonSeedCreator.h:181
MuonSeedCreator::SMB30
std::vector< double > SMB30
Definition: MuonSeedCreator.h:130
PV3DBase::x
T x() const
Definition: PV3DBase.h:59
DiDispStaMuonMonitor_cfi.pt
pt
Definition: DiDispStaMuonMonitor_cfi.py:39
MuonSeedCreator::CSC12_3
std::vector< double > CSC12_3
Definition: MuonSeedCreator.h:139
h
FWCore Framework interface EventSetupRecordImplementation h
Helper function to determine trigger accepts.
Definition: L1TUtmAlgorithmRcd.h:4
MuonSeedCreator::debug
bool debug
Definition: MuonSeedCreator.h:86
PV3DBase::theta
Geom::Theta< T > theta() const
Definition: PV3DBase.h:72
trajectoryStateTransform::persistentState
PTrajectoryStateOnDet persistentState(const TrajectoryStateOnSurface &ts, unsigned int detid)
Definition: TrajectoryStateTransform.cc:14
MuonSeedCreator::theMinMomentum
float theMinMomentum
Definition: MuonSeedCreator.h:78
gather_cfg.cout
cout
Definition: gather_cfg.py:144
MuonSeedCreator::DT13_1
std::vector< double > DT13_1
Definition: MuonSeedCreator.h:150
MuonSeedCreator::estimatePtSingle
void estimatePtSingle(const SegmentContainer &seg, const std::vector< int > &layers, double &pt, double &spt)
Estimate transverse momentum of track from single segment.
Definition: MuonSeedCreator.cc:773
MuonSeedCreator::OL1213
std::vector< double > OL1213
Definition: MuonSeedCreator.h:109
MuonSeedCreator::DT13
std::vector< double > DT13
Definition: MuonSeedCreator.h:93
MuonSeedCreator::OL_1222
std::vector< double > OL_1222
Definition: MuonSeedCreator.h:162
Validation_hcalonly_cfi.sign
sign
Definition: Validation_hcalonly_cfi.py:32
Geom::Spherical2Cartesian
Definition: CoordinateSets.h:58
MuonSeedCreator::CSC12_2
std::vector< double > CSC12_2
Definition: MuonSeedCreator.h:138
hltPixelTracks_cff.chi2
chi2
Definition: hltPixelTracks_cff.py:25
HLT_FULL_cff.dPhi
dPhi
Definition: HLT_FULL_cff.py:13765
MuonSeedCreator::CSC01_1
std::vector< double > CSC01_1
Definition: MuonSeedCreator.h:136
relativeConstraints.error
error
Definition: relativeConstraints.py:53
MuonSeedCreator::DT23_1
std::vector< double > DT23_1
Definition: MuonSeedCreator.h:154
MuonSeedCreator::SME11
std::vector< double > SME11
Definition: MuonSeedCreator.h:115
LocalTrajectoryParameters
Definition: LocalTrajectoryParameters.h:25
DetId
Definition: DetId.h:17
TrajectoryStateOnSurface
Definition: TrajectoryStateOnSurface.h:16
dqmdumpme.last
last
Definition: dqmdumpme.py:56
MuonSeedCreator::estimatePtCSC
void estimatePtCSC(const SegmentContainer &seg, const std::vector< int > &layers, double &pt, double &spt)
Estimate transverse momentum of track from CSC measurements.
Definition: MuonSeedCreator.cc:324
RandomServiceHelper.t1
t1
Definition: RandomServiceHelper.py:256
MuonSeedCreator::SMB20
std::vector< double > SMB20
Definition: MuonSeedCreator.h:127
MuonSeedCreator::DT14_1
std::vector< double > DT14_1
Definition: MuonSeedCreator.h:152
AlgebraicVector
CLHEP::HepVector AlgebraicVector
Definition: AlgebraicObjects.h:13
clone
TEveGeoShape * clone(const TEveElement *element, TEveElement *parent)
Definition: eve_macros.cc:135
MuonSeedCreator::SMB11
std::vector< double > SMB11
Definition: MuonSeedCreator.h:125
PVValHelper::eta
Definition: PVValidationHelpers.h:69
MuonSeedCreator::defaultMomentum
float defaultMomentum
Definition: MuonSeedCreator.h:80
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
MuonSeedCreator::SMB_30S
std::vector< double > SMB_30S
Definition: MuonSeedCreator.h:173
MuonSeedCreator::CSC12
std::vector< double > CSC12
Definition: MuonSeedCreator.h:100
MuonSeedCreator::SME_13S
std::vector< double > SME_13S
Definition: MuonSeedCreator.h:179
h
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
MuonSeedCreator::SMB31
std::vector< double > SMB31
Definition: MuonSeedCreator.h:131
Point3DBase< float, LocalTag >
MuonSeedCreator::SME22
std::vector< double > SME22
Definition: MuonSeedCreator.h:119
MuonSeedCreator::DT34_1
std::vector< double > DT34_1
Definition: MuonSeedCreator.h:158
MuonSeedCreator::OL1232
std::vector< double > OL1232
Definition: MuonSeedCreator.h:111
MuonSeedCreator::CSC03
std::vector< double > CSC03
Definition: MuonSeedCreator.h:103
MuonSeedCreator::estimatePtShowering
void estimatePtShowering(int &NShowers, int &NShowerSeg, double &pt, double &spt)
Estimate transverse momentum of track from showering segment.
Definition: MuonSeedCreator.cc:909
MuonSeedCreator::OL_2222
std::vector< double > OL_2222
Definition: MuonSeedCreator.h:165
MuonSeedCreator::CSC13_3
std::vector< double > CSC13_3
Definition: MuonSeedCreator.h:141
MuonSeedCreator::SME13
std::vector< double > SME13
Definition: MuonSeedCreator.h:117
ALCARECOTkAlJpsiMuMu_cff.charge
charge
Definition: ALCARECOTkAlJpsiMuMu_cff.py:47
LocalTrajectoryError
Definition: LocalTrajectoryError.h:20
runTauDisplay.gp
gp
Definition: runTauDisplay.py:431
MuonSeedCreator::CSC24
std::vector< double > CSC24
Definition: MuonSeedCreator.h:106
MuonSeedCreator::DT24_2
std::vector< double > DT24_2
Definition: MuonSeedCreator.h:157
MuonSeedCreator::DT23_2
std::vector< double > DT23_2
Definition: MuonSeedCreator.h:155
MuonSeedCreator::SMB10
std::vector< double > SMB10
Definition: MuonSeedCreator.h:124
type
type
Definition: SiPixelVCal_PayloadInspector.cc:37
MuonSeedCreator::theMaxMomentum
float theMaxMomentum
Definition: MuonSeedCreator.h:79
MuonSeedCreator::CSC12_1
std::vector< double > CSC12_1
Definition: MuonSeedCreator.h:137
MuonSeedCreator::DT12_2
std::vector< double > DT12_2
Definition: MuonSeedCreator.h:149
MuonSeedCreator::sysError
double sysError
Definition: MuonSeedCreator.h:83
MuonSeedCreator::OL_2213
std::vector< double > OL_2213
Definition: MuonSeedCreator.h:164
PV3DBase::eta
T eta() const
Definition: PV3DBase.h:73
MuonSeedCreator::DT23
std::vector< double > DT23
Definition: MuonSeedCreator.h:95
PV3DBase::y
T y() const
Definition: PV3DBase.h:60
funct::tan
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
MuonSeedCreator::DT14_2
std::vector< double > DT14_2
Definition: MuonSeedCreator.h:153
MuonSeedCreator::SMB32
std::vector< double > SMB32
Definition: MuonSeedCreator.h:132
MuonSeedCreator::SMB12
std::vector< double > SMB12
Definition: MuonSeedCreator.h:126
CSC01
Definition: QTest.h:565
MuonSeedCreator::scaledPhi
double scaledPhi(double dphi, double t1)
Scale the dPhi from segment position.
Definition: MuonSeedCreator.cc:1013
MuonSeedCreator::CSC24_1
std::vector< double > CSC24_1
Definition: MuonSeedCreator.h:145
MuonSeedCreator::SMB_11S
std::vector< double > SMB_11S
Definition: MuonSeedCreator.h:168
MuonSeedCreator::SME32
std::vector< double > SME32
Definition: MuonSeedCreator.h:121
AlgebraicSymMatrix
CLHEP::HepSymMatrix AlgebraicSymMatrix
Definition: AlgebraicObjects.h:15
MuonSeedCreator::CSC14
std::vector< double > CSC14
Definition: MuonSeedCreator.h:104
MuonSeedCreator::SME_12S
std::vector< double > SME_12S
Definition: MuonSeedCreator.h:178
MuonSeedCreator::OL2222
std::vector< double > OL2222
Definition: MuonSeedCreator.h:113
MuonSeedCreator::SMB_10S
std::vector< double > SMB_10S
Definition: MuonSeedCreator.h:167
cmsLHEtoEOSManager.l
l
Definition: cmsLHEtoEOSManager.py:204
MuonSeedCreator::SMB22
std::vector< double > SMB22
Definition: MuonSeedCreator.h:129
MuonSeedCreator::OL_1213
std::vector< double > OL_1213
Definition: MuonSeedCreator.h:161
MuonSeedCreator::DT24_1
std::vector< double > DT24_1
Definition: MuonSeedCreator.h:156
MuonSeedCreator::SME31
std::vector< double > SME31
Definition: MuonSeedCreator.h:120
MuonSeedCreator::SMB_12S
std::vector< double > SMB_12S
Definition: MuonSeedCreator.h:169
MuonSeedCreator::DT24
std::vector< double > DT24
Definition: MuonSeedCreator.h:96
MuonSeedCreator::SME21
std::vector< double > SME21
Definition: MuonSeedCreator.h:118
MuonSeedCreator::CSC13
std::vector< double > CSC13
Definition: MuonSeedCreator.h:102
MuonSeedCreator::estimatePtOverlap
void estimatePtOverlap(const SegmentContainer &seg, const std::vector< int > &layers, double &pt, double &spt)
Estimate transverse momentum of track from CSC + DT measurements.
Definition: MuonSeedCreator.cc:643
MuonSeedCreator::BField
const MagneticField * BField
Definition: MuonSeedCreator.h:89
MuonSeedCreator::SMB_31S
std::vector< double > SMB_31S
Definition: MuonSeedCreator.h:174
MuonSeedCreator::OL1222
std::vector< double > OL1222
Definition: MuonSeedCreator.h:110
MuonSeedCreator::DT12_1
std::vector< double > DT12_1
Definition: MuonSeedCreator.h:148
MuonSeedCreator::CSC23_2
std::vector< double > CSC23_2
Definition: MuonSeedCreator.h:144
MuonSeedCreator::OL2213
std::vector< double > OL2213
Definition: MuonSeedCreator.h:112
TrajectorySeed
Definition: TrajectorySeed.h:18
MuonSeedCreator::DT13_2
std::vector< double > DT13_2
Definition: MuonSeedCreator.h:151
MuonSeedCreator::SMB_20S
std::vector< double > SMB_20S
Definition: MuonSeedCreator.h:170
MuonSeedCreator::SegmentContainer
MuonTransientTrackingRecHit::MuonRecHitContainer SegmentContainer
Definition: MuonSeedCreator.h:32
MuonSeedCreator::SME12
std::vector< double > SME12
Definition: MuonSeedCreator.h:116
MuonSeedCreator::getPt
std::vector< double > getPt(const std::vector< double > &vParameters, double eta, double dPhi)
Compute pt from parameters.
Definition: MuonSeedCreator.cc:1001
edm::OwnVector::push_back
void push_back(D *&d)
Definition: OwnVector.h:326
MuonSeedCreator::DT34
std::vector< double > DT34
Definition: MuonSeedCreator.h:97
MuonSeedCreator::SMB_21S
std::vector< double > SMB_21S
Definition: MuonSeedCreator.h:171
MuonSeedCreator::CSC02
std::vector< double > CSC02
Definition: MuonSeedCreator.h:101
MuonSeedCreator::CSC13_2
std::vector< double > CSC13_2
Definition: MuonSeedCreator.h:140
MuonSeedCreator::CSC23_1
std::vector< double > CSC23_1
Definition: MuonSeedCreator.h:143
MuonSeedCreator::SME_21S
std::vector< double > SME_21S
Definition: MuonSeedCreator.h:180
PTrajectoryStateOnDet
Definition: PTrajectoryStateOnDet.h:10
MuonSeedCreator::DT34_2
std::vector< double > DT34_2
Definition: MuonSeedCreator.h:159
dqmiolumiharvest.j
j
Definition: dqmiolumiharvest.py:66
MuonSeedCreator::DT14
std::vector< double > DT14
Definition: MuonSeedCreator.h:94
submitPVValidationJobs.t
string t
Definition: submitPVValidationJobs.py:644
MuonSeedCreator::estimatePtDT
void estimatePtDT(const SegmentContainer &seg, const std::vector< int > &layers, double &pt, double &spt)
Estimate transverse momentum of track from DT measurements.
Definition: MuonSeedCreator.cc:494
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:69
MuonSeedCreator::weightedPt
void weightedPt(const std::vector< double > &ptEstimate, const std::vector< double > &sptEstimate, double &ptAvg, double &sptAvg)
Compute weighted mean pt from different pt estimators.
Definition: MuonSeedCreator.cc:934
MuonSeedCreator::CSC34
std::vector< double > CSC34
Definition: MuonSeedCreator.h:107
MuonSeedCreator::SME_11S
std::vector< double > SME_11S
Definition: MuonSeedCreator.h:177
alongMomentum
Definition: PropagationDirection.h:4
hgcalTopologyTester_cfi.layers
layers
Definition: hgcalTopologyTester_cfi.py:8
PV3DBase::phi
Geom::Phi< T > phi() const
Definition: PV3DBase.h:66
cuy.dh
dh
Definition: cuy.py:355
hit
Definition: SiStripHitEffFromCalibTree.cc:88
MuonSeedCreator::SMB_22S
std::vector< double > SMB_22S
Definition: MuonSeedCreator.h:172
muonDTDigis_cfi.pset
pset
Definition: muonDTDigis_cfi.py:27
edm::OwnVector< TrackingRecHit >
MuonSeedCreator::SMB21
std::vector< double > SMB21
Definition: MuonSeedCreator.h:128
MuonSeedCreator::OL_1232
std::vector< double > OL_1232
Definition: MuonSeedCreator.h:163
MuonSeedCreator::SME41
std::vector< double > SME41
Definition: MuonSeedCreator.h:122
findQualityFiles.size
size
Write out results.
Definition: findQualityFiles.py:443
MuonSeedCreator::CSC14_3
std::vector< double > CSC14_3
Definition: MuonSeedCreator.h:142