test
CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Types | Private Member Functions | Private Attributes
MuonSeedPtExtractor Class Reference

#include <MuonSeedPtExtractor.h>

Public Member Functions

 MuonSeedPtExtractor (const edm::ParameterSet &)
 Constructor with Parameter set and MuonServiceProxy. More...
 
virtual std::vector< double > pT_extract (MuonTransientTrackingRecHit::ConstMuonRecHitPointer firstHit, MuonTransientTrackingRecHit::ConstMuonRecHitPointer secondHit) const
 
void setBeamSpot (const GlobalVector &gv)
 
virtual ~MuonSeedPtExtractor ()
 Destructor. More...
 

Private Types

typedef std::map< std::string,
std::vector< double > > 		ParametersMap
 
typedef std::map< std::string,
std::vector< double > > 		ScalesMap
 

Private Member Functions

void fillParametersForCombo (const std::string &name, const edm::ParameterSet &pset)
 
void fillScalesForCombo (const std::string &name, const edm::ParameterSet &pset)
 
std::vector< double > getPt (const std::vector< double > &vPara, double eta, double dPhi) const
 
std::vector< double > getPt (const std::vector< double > &vPara, double eta, double dPhi, const std::string &combination, const DTChamberId &outerDetId) const
 
void init (const edm::ParameterSet &par)
 
int stationCode (MuonTransientTrackingRecHit::ConstMuonRecHitPointer hit) const
 

Private Attributes

bool scaleDT_
 
GlobalVector theBeamSpot
 
ParametersMap theParametersForCombo
 
ScalesMap theScalesForCombo
 

Detailed Description

Definition at line 15 of file MuonSeedPtExtractor.h.

Member Typedef Documentation

typedef std::map<std::string, std::vector<double> > MuonSeedPtExtractor::ParametersMap
private

Definition at line 45 of file MuonSeedPtExtractor.h.

typedef std::map<std::string, std::vector<double> > MuonSeedPtExtractor::ScalesMap
private

Definition at line 46 of file MuonSeedPtExtractor.h.

Constructor & Destructor Documentation

MuonSeedPtExtractor::MuonSeedPtExtractor ( const edm::ParameterSet &  par)

Constructor with Parameter set and MuonServiceProxy.

Definition at line 9 of file MuonSeedPtExtractor.cc.

References init().

10 : theBeamSpot(0.,0.,0.),
11  scaleDT_( par.getParameter<bool>("scaleDT") )
12 {
13  init(par);
14 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
MuonSeedPtExtractor::init
void init(const edm::ParameterSet &par)
Definition: MuonSeedPtExtractor.cc:17
MuonSeedPtExtractor::~MuonSeedPtExtractor ( )
virtual

Destructor.

Definition at line 113 of file MuonSeedPtExtractor.cc.

113  {
114 }

Member Function Documentation

void MuonSeedPtExtractor::fillParametersForCombo ( const std::string &  name,
const edm::ParameterSet &  pset 
)
private

Definition at line 117 of file MuonSeedPtExtractor.cc.

References edm::ParameterSet::getParameter(), mergeVDriftHistosByStation::name, and theParametersForCombo.

Referenced by init().

118 {
119  theParametersForCombo[name] = pset.getParameter<std::vector<double> >(name);
120 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
MuonSeedPtExtractor::theParametersForCombo
ParametersMap theParametersForCombo
Definition: MuonSeedPtExtractor.h:47
void MuonSeedPtExtractor::fillScalesForCombo ( const std::string &  name,
const edm::ParameterSet &  pset 
)
private

Definition at line 123 of file MuonSeedPtExtractor.cc.

References edm::ParameterSet::getParameter(), mergeVDriftHistosByStation::name, and theScalesForCombo.

Referenced by init().

124 {
125  theScalesForCombo[name] = pset.getParameter<std::vector<double> >(name);
126 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
std::vector< double > MuonSeedPtExtractor::getPt ( const std::vector< double > &  vPara,
double  eta,
double  dPhi 
) const
private

Definition at line 381 of file MuonSeedPtExtractor.cc.

References h.

Referenced by pT_extract().

381  {
382  //std::cout<<" eta = "<<eta<<" dPhi = "<<dPhi<<" vPara[0] = "<<vPara[0]<<" vPara[1] = "<<vPara[1]<<" vPara[2] = "<<vPara[2]<<std::endl;
383  double h = fabs(eta);
384  double estPt = ( vPara[0] + vPara[1]*h + vPara[2]*h*h ) / dPhi;
385  double estSPt = ( vPara[3] + vPara[4]*h + vPara[5]*h*h ) / dPhi;
386  // std::cout<<"estPt = "<<estPt<<std::endl;
387  std::vector<double> paraPt ;
388  paraPt.push_back( estPt );
389  paraPt.push_back( estSPt ) ;
390  return paraPt ;
391 }
h
FWCore Framework interface EventSetupRecordImplementation h
Helper function to determine trigger accepts.
Definition: L1TUtmAlgorithmRcd.h:4
dPhi
double dPhi(double phi1, double phi2)
Definition: JetUtil.h:30
std::vector< double > MuonSeedPtExtractor::getPt ( const std::vector< double > &  vPara,
double  eta,
double  dPhi,
const std::string &  combination,
const DTChamberId &  outerDetId 
) const
private

Definition at line 394 of file MuonSeedPtExtractor.cc.

References funct::abs(), h, theScalesForCombo, and DTChamberId::wheel().

395  {
396  double h = fabs(eta);
397  double estPt = ( vPara[0] + vPara[1]*h + vPara[2]*h*h ) / dPhi;
398  // changed by S.C.
399  double estSPt = ( vPara[3] + vPara[4]*h + vPara[5]*h*h ) / dPhi;
400 
401  // scale the pt and spt , changed by S.C.
402  int wheel = 0;
403  if(combination[0] == 'D') {
404  wheel = abs(outerDetId.wheel());
405  }
406 
407  std::ostringstream os;
408  os << combination << "_" << wheel << "_scale";
409 
410  ScalesMap::const_iterator scalesItr = theScalesForCombo.find(os.str());
411  if ( scalesItr != theScalesForCombo.end()) {
412  double t1 = scalesItr->second[0];
413  double scaleFactor = 1./( 1. + t1/( estPt + 10. ) ) ;
414 
415  estSPt = estSPt * scaleFactor ;
416  estPt = estPt * scaleFactor ;
417  }
418 
419  // std::cout<<"estPt = "<<estPt<<std::endl;
420  std::vector<double> paraPt ;
421  paraPt.push_back( estPt );
422  paraPt.push_back( estSPt ) ;
423  return paraPt ;
424 }
h
FWCore Framework interface EventSetupRecordImplementation h
Helper function to determine trigger accepts.
Definition: L1TUtmAlgorithmRcd.h:4
dPhi
double dPhi(double phi1, double phi2)
Definition: JetUtil.h:30
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
DTChamberId::wheel
int wheel() const
Return the wheel number.
Definition: DTChamberId.h:45
void MuonSeedPtExtractor::init ( const edm::ParameterSet &  par)
private

Definition at line 17 of file MuonSeedPtExtractor.cc.

References fillParametersForCombo(), and fillScalesForCombo().

Referenced by MuonSeedPtExtractor().

18 {
19  // load pT seed parameters
20  // DT combinations
21  fillParametersForCombo("DT_12", par);
22  fillParametersForCombo("DT_13", par);
23  fillParametersForCombo("DT_14", par);
24  fillParametersForCombo("DT_23", par);
25  fillParametersForCombo("DT_24", par);
26  fillParametersForCombo("DT_34", par);
27  // CSC combinations
28  fillParametersForCombo("CSC_01", par);
29  fillParametersForCombo("CSC_12", par);
30  fillParametersForCombo("CSC_02", par);
31  fillParametersForCombo("CSC_13", par);
32  fillParametersForCombo("CSC_03", par);
33  fillParametersForCombo("CSC_14", par);
34  fillParametersForCombo("CSC_23", par);
35  fillParametersForCombo("CSC_24", par);
36  fillParametersForCombo("CSC_34", par);
37 
38  // Overlap combinations
39  fillParametersForCombo("OL_1213", par);
40  fillParametersForCombo("OL_1222", par);
41  fillParametersForCombo("OL_1232", par);
42  fillParametersForCombo("OL_2213", par);
43  fillParametersForCombo("OL_2222", par);
44 
45  // Single segments (CSC)
46  fillParametersForCombo("SME_11", par);
47  fillParametersForCombo("SME_12", par);
48  fillParametersForCombo("SME_13", par);
49  fillParametersForCombo("SME_21", par);
50  fillParametersForCombo("SME_22", par);
51  fillParametersForCombo("SME_31", par);
52  fillParametersForCombo("SME_32", par);
53  fillParametersForCombo("SME_41", par);
54  fillParametersForCombo("SME_42", par);
55 
56  // Single segments (DT)
57  fillParametersForCombo("SMB_10", par);
58  fillParametersForCombo("SMB_11", par);
59  fillParametersForCombo("SMB_12", par);
60  fillParametersForCombo("SMB_20", par);
61  fillParametersForCombo("SMB_21", par);
62  fillParametersForCombo("SMB_22", par);
63  fillParametersForCombo("SMB_30", par);
64  fillParametersForCombo("SMB_31", par);
65  fillParametersForCombo("SMB_32", par);
66 
67  fillScalesForCombo("CSC_01_1_scale", par);
68  fillScalesForCombo("CSC_12_1_scale", par);
69  fillScalesForCombo("CSC_12_2_scale", par);
70  fillScalesForCombo("CSC_12_3_scale", par);
71  fillScalesForCombo("CSC_13_2_scale", par);
72  fillScalesForCombo("CSC_13_3_scale", par);
73  fillScalesForCombo("CSC_14_3_scale", par);
74  fillScalesForCombo("CSC_23_1_scale", par);
75  fillScalesForCombo("CSC_23_2_scale", par);
76  fillScalesForCombo("CSC_24_1_scale", par);
77  fillScalesForCombo("CSC_34_1_scale", par);
78  fillScalesForCombo("DT_12_1_scale", par);
79  fillScalesForCombo("DT_12_2_scale", par);
80  fillScalesForCombo("DT_13_1_scale", par);
81  fillScalesForCombo("DT_13_2_scale", par);
82  fillScalesForCombo("DT_14_1_scale", par);
83  fillScalesForCombo("DT_14_2_scale", par);
84  fillScalesForCombo("DT_23_1_scale", par);
85  fillScalesForCombo("DT_23_2_scale", par);
86  fillScalesForCombo("DT_24_1_scale", par);
87  fillScalesForCombo("DT_24_2_scale", par);
88  fillScalesForCombo("DT_34_1_scale", par);
89  fillScalesForCombo("DT_34_2_scale", par);
90  fillScalesForCombo("OL_1213_0_scale", par);
91  fillScalesForCombo("OL_1222_0_scale", par);
92  fillScalesForCombo("OL_1232_0_scale", par);
93  fillScalesForCombo("OL_2213_0_scale", par);
94  fillScalesForCombo("OL_2222_0_scale", par);
95  fillScalesForCombo("SMB_10_0_scale", par);
96  fillScalesForCombo("SMB_11_0_scale", par);
97  fillScalesForCombo("SMB_12_0_scale", par);
98  fillScalesForCombo("SMB_20_0_scale", par);
99  fillScalesForCombo("SMB_21_0_scale", par);
100  fillScalesForCombo("SMB_22_0_scale", par);
101  fillScalesForCombo("SMB_30_0_scale", par);
102  fillScalesForCombo("SMB_31_0_scale", par);
103  fillScalesForCombo("SMB_32_0_scale", par);
104  fillScalesForCombo("SME_11_0_scale", par);
105  fillScalesForCombo("SME_12_0_scale", par);
106  fillScalesForCombo("SME_13_0_scale", par);
107  fillScalesForCombo("SME_21_0_scale", par);
108  fillScalesForCombo("SME_22_0_scale", par);
109 
110 }
MuonSeedPtExtractor::fillParametersForCombo
void fillParametersForCombo(const std::string &name, const edm::ParameterSet &pset)
Definition: MuonSeedPtExtractor.cc:117
MuonSeedPtExtractor::fillScalesForCombo
void fillScalesForCombo(const std::string &name, const edm::ParameterSet &pset)
Definition: MuonSeedPtExtractor.cc:123
std::vector< double > MuonSeedPtExtractor::pT_extract ( MuonTransientTrackingRecHit::ConstMuonRecHitPointer  firstHit,
MuonTransientTrackingRecHit::ConstMuonRecHitPointer  secondHit 
) const
virtual

If firstHit and secondHit are the same, it calculates the pT from the phi bend of the single segment

Definition at line 129 of file MuonSeedPtExtractor.cc.

References funct::abs(), dPhi(), alignCSCRings::e, eta, PV3DBase< T, PVType, FrameType >::eta(), Exception, getPt(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), Pi, relativeConstraints::ring, CSCDetId::ring(), scaleDT_, jetcorrextractor::sign(), mathSSE::sqrt(), DTChamberId::station(), CSCDetId::station(), stationCode(), AlCaHLTBitMon_QueryRunRegistry::string, theBeamSpot, theParametersForCombo, DTChamberId::wheel(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by MuonDTSeedFromRecHits::computePtWithoutVtx(), MuonDTSeedFromRecHits::computePtWithVtx(), MuonCSCSeedFromRecHits::createDefaultEndcapSeed(), SETSeedFinder::estimateMomentum(), MuonOverlapSeedFromRecHits::makeSeed(), and MuonCSCSeedFromRecHits::makeSeed().

131 {
132  GlobalPoint innerPoint = firstHit->globalPosition() - theBeamSpot;
133  GlobalPoint outerPoint = secondHit->globalPosition() - theBeamSpot;
134  MuonTransientTrackingRecHit::ConstMuonRecHitPointer innerHit = firstHit;
135  MuonTransientTrackingRecHit::ConstMuonRecHitPointer outerHit = secondHit;
136 
137  // ways in which the hits could be in the wrong order
138  if( (outerHit->isDT() && innerHit->isCSC())
139  || (outerHit->isDT() && innerHit->isDT() && (innerPoint.perp() > outerPoint.perp()))
140  || (outerHit->isCSC() && innerHit->isCSC() && (fabs(innerPoint.z()) > fabs(outerPoint.z()))) )
141  {
142  innerHit = secondHit;
143  outerHit = firstHit;
144  innerPoint = innerHit->globalPosition() - theBeamSpot;
145  outerPoint = outerHit->globalPosition() - theBeamSpot;
146  }
147 
148  double phiInner = innerPoint.phi();
149  double phiOuter = outerPoint.phi();
150 
151  double etaInner = innerPoint.eta();
152  double etaOuter = outerPoint.eta();
153  //std::cout<<" inner pos = "<< innerPoint << " phi eta " << phiInner << " " << etaInner << std::endl;
154  //std::cout<<" outer pos = "<< outerPoint << " phi eta " << phiOuter << " " << etaOuter << std::endl;
155  //double thetaInner = firstHit->globalPosition().theta();
156  // if some of the segments is missing r-phi measurement then we should
157  // use only the 4D phi estimate (also use 4D eta estimate only)
158  // the direction is not so important (it will be corrected)
159  /*
160  bool firstOK = (4==allValidSets[iSet][firstMeasurement]->hit->dimension());
161  bool lastOK = (4==allValidSets[iSet][lastMeasurement]->hit->dimension());
162  if(!(firstOK * lastOK)){
163  if(!firstOK){
164  }
165  if(!firstOK){
166  }
167  }
168  */
169  double dPhi = phiInner - phiOuter;
170  if(dPhi < -TMath::Pi()){
171  dPhi += 2*TMath::Pi();
172  }
173  else if(dPhi > TMath::Pi()){
174  dPhi -= 2*TMath::Pi();
175  }
176  int sign = 1;
177  if ( dPhi< 0.) {
178  dPhi = -dPhi;
179  sign = -1;
180  }
181 
182  if (dPhi < 1.0e-6){
183  dPhi = 1.0e-6;
184  }
185  double eta = fabs(etaOuter);// what if it is 2D segment? use inner?
186 
187 
188  std::vector <int> stationCoded(2);
189  stationCoded[0] = stationCoded[1] = 999;
190 
191  DetId detId_inner = innerHit->hit()->geographicalId();
192  DetId detId_outer = outerHit->hit()->geographicalId();
193 
194  stationCoded[0] = stationCode(innerHit);
195  stationCoded[1] = stationCode(outerHit);
196 
197  std::ostringstream os0 ;
198  std::ostringstream os1 ;
199  os0 << abs(stationCoded[0]);
200  os1 << abs(stationCoded[1]);
201 
202  //std::cout<<" st1 = "<<stationCoded[0]<<" st2 = "<<stationCoded[1]<<std::endl;
203  //std::cout<<" detId_inner = "<<detId_inner.rawId()<<" detId_outer = "<< detId_outer.rawId()<<std::endl;
204  std::string combination = "0";
205  std::string init_combination = combination;
206  bool singleSegment = false;
207  //if(detId_first == detId_second){
208  if( stationCoded[0] == stationCoded[1]){
209  // single segment - DT or CSC
210  singleSegment = true;
211  //eta = innerPoint.eta();
212  GlobalVector gv = innerHit->globalDirection();
213  double gvPerpPos = gv.x()*gv.x() + gv.y()*gv.y();
214  if (gvPerpPos < 1e-32) gvPerpPos = 1e-32; gvPerpPos=sqrt(gvPerpPos);
215  // Psi is angle between the segment origin and segment direction
216  // Use dot product between two vectors to get Psi in global x-y plane
217  double cosDpsi = (gv.x()*innerPoint.x() + gv.y()*innerPoint.y());
218  if (cosDpsi!=0){
219  cosDpsi /= sqrt(innerPoint.x()*innerPoint.x() + innerPoint.y()*innerPoint.y());
220  cosDpsi /= gvPerpPos;
221  cosDpsi = cosDpsi > 1 ? 1 : cosDpsi;
222  cosDpsi = cosDpsi < -1 ? -1 : cosDpsi;
223  }
224 
225  double axb = ( innerPoint.x()*gv.y() ) - ( innerPoint.y()*gv.x() ) ;
226  sign = (axb < 0.) ? 1 : -1;
227 
228  double dpsi = fabs(acos(cosDpsi)) ;
229  if ( dpsi > TMath::Pi()/2.) {
230  dpsi = TMath::Pi() - dpsi;
231  }
232  if (fabs(dpsi) < 0.00005) {
233  dpsi = 0.00005;
234  }
235  dPhi = dpsi;
236 
237  if(innerHit->isDT())
238  {
239  DTChamberId dtCh(detId_inner);
240  std::ostringstream os;
241  os << dtCh.station() << abs(dtCh.wheel());
242  combination = "SMB_"+os.str();
243  }
244  else if(innerHit->isCSC())
245  {
246  CSCDetId cscId(detId_inner);
247  std::ostringstream os;
248  int ring = cscId.ring();
249  if(ring == 4) ring = 1;
250  os << cscId.station() << ring;
251  combination = "SME_"+os.str();
252  }
253  else
254  {
255  throw cms::Exception("MuonSeedPtExtractor") << "Bad hit DetId";
256  }
257  }
258  else{
259  if(stationCoded[0]<0){
260  if(stationCoded[1]<0){
261  // DT-DT
262  combination = "DT_" + os0.str() + os1.str();
263  }
264  else{
265  // DT-CSC
266  eta = fabs(etaInner);
267  if(-1==stationCoded[0]){
268  switch (stationCoded[1]){
269  case 1:
270  combination = "OL_1213";
271  break;
272  case 2:
273  combination = "OL_1222";
274  break;
275  case 3:
276  combination = "OL_1232";
277  break;
278  default:
279  // can not be
280  break;
281  }
282  }
283  else if (-2==stationCoded[0]){
284  if(1==stationCoded[1]){
285  combination = "OL_2213";
286  }
287  else{
288  // can not be (not coded?)
289  combination = "OL_2222";// in case
290  }
291  }
292  else{
293  // can not be
294  }
295  }
296  }
297  else{
298  if(stationCoded[1]<0){
299  // CSC-DT
300  // can not be
301  }
302  else{
303  // CSC-CSC
304  combination = "CSC_" + os0.str() + os1.str();
305  if("CSC_04" == combination){
306  combination = "CSC_14";
307  }
308  }
309  }
310  }
311 
312  std::vector<double> pTestimate(2);//
313  //std::cout<<" combination = "<<combination<<std::endl;
314  if(init_combination!=combination){
315  //std::cout<<" combination = "<<combination<<" eta = "<<eta<<" dPhi = "<<dPhi<<std::endl;
316  ParametersMap::const_iterator parametersItr = theParametersForCombo.find(combination);
317  if(parametersItr == theParametersForCombo.end()) {
318  // edm::LogWarning("RecoMuon|MuonSeedGenerator|MuonSeedPtExtractor") << "Cannot find parameters for combo " << combination;
319  edm::LogWarning("BadSegmentCombination") << "Cannot find parameters for combo " << combination;
320  pTestimate[0] = pTestimate[1] = 100;
321  // throw cms::Exception("MuonSeedPtEstimator") << "Cannot find parameters for combo " << combination;
322  } else {
323 
324  if(scaleDT_ && outerHit->isDT() )
325  {
326  pTestimate = getPt(parametersItr->second, eta, dPhi, combination, detId_outer);
327  }
328  else
329  {
330  pTestimate = getPt(parametersItr->second, eta, dPhi);
331  }
332 
333  if(singleSegment){
334  pTestimate[0] = fabs(pTestimate[0]);
335  pTestimate[1] = fabs(pTestimate[1]);
336  }
337  pTestimate[0] *= double(sign);
338  }
339  }
340  else{
341  // often a MB3 - ME1/3 seed
342  pTestimate[0] = pTestimate[1] = 100;
343  // hmm
344  }
345  // std::cout<<" pTestimate[0] = "<<pTestimate[0]<<" pTestimate[1] = "<<pTestimate[1]<<std::endl;
346  /*
347  MuonRecHitContainer_clusters[cluster][iHit+1]->isDT());
348  if(specialCase){
349  DTChamberId dtCh(detId);
350  DTChamberId dtCh_2(detId_2);
351  specialCase = (dtCh.station() == dtCh_2.station());
352  }
353  */
354  //return vPara;
355  return pTestimate;
356 }
Pi
const double Pi
Definition: CosmicMuonParameters.h:18
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:72
jetcorrextractor::sign
double sign(double x)
Definition: JetCorrExtractorT.h:32
PV3DBase::phi
Geom::Phi< T > phi() const
Definition: PV3DBase.h:69
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
dPhi
double dPhi(double phi1, double phi2)
Definition: JetUtil.h:30
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
MuonSeedPtExtractor::stationCode
int stationCode(MuonTransientTrackingRecHit::ConstMuonRecHitPointer hit) const
Definition: MuonSeedPtExtractor.cc:359
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
MuonSeedPtExtractor::theParametersForCombo
ParametersMap theParametersForCombo
Definition: MuonSeedPtExtractor.h:47
DetId
Definition: DetId.h:18
PV3DBase::eta
T eta() const
Definition: PV3DBase.h:76
MuonTransientTrackingRecHit::ConstMuonRecHitPointer
std::shared_ptr< MuonTransientTrackingRecHit const > ConstMuonRecHitPointer
Definition: MuonTransientTrackingRecHit.h:24
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
MuonSeedPtExtractor::getPt
std::vector< double > getPt(const std::vector< double > &vPara, double eta, double dPhi) const
Definition: MuonSeedPtExtractor.cc:381
void MuonSeedPtExtractor::setBeamSpot ( const GlobalVector &  gv)
inline

Definition at line 30 of file MuonSeedPtExtractor.h.

References theBeamSpot.

Referenced by MuonSeedVFinder::setBeamSpot().

30 {theBeamSpot = gv;}
int MuonSeedPtExtractor::stationCode ( MuonTransientTrackingRecHit::ConstMuonRecHitPointer  hit) const
private

Definition at line 359 of file MuonSeedPtExtractor.cc.

References mps_fire::result, CSCDetId::ring(), DTChamberId::station(), and CSCDetId::station().

Referenced by pT_extract().

360 {
361  DetId detId(hit->hit()->geographicalId());
362  int result = -999;
363  if(hit->isDT() ){
364  DTChamberId dtCh(detId);
365  //std::cout<<"first (DT) St/W/S = "<<dtCh.station()<<"/"<<dtCh.wheel()<<"/"<<dtCh.sector()<<"/"<<std::endl;
366  result = -1 * dtCh.station();
367  }
368  else if( hit->isCSC() ){
369  CSCDetId cscID(detId);
370  //std::cout<<"first (CSC) E/S/R/C = "<<cscID.endcap()<<"/"<<cscID.station()<<"/"<<cscID.ring()<<"/"<<cscID.chamber()<<std::endl;
371  result = cscID.station();
372  if(result == 1 && (1 == cscID.ring() || 4 == cscID.ring()) )
373  result = 0;
374  }
375  else if(hit->isRPC()){
376  }
377  return result;
378 }
mps_fire.result
tuple result
Definition: mps_fire.py:84
DetId
Definition: DetId.h:18

Member Data Documentation

bool MuonSeedPtExtractor::scaleDT_
private

Definition at line 50 of file MuonSeedPtExtractor.h.

Referenced by pT_extract().

GlobalVector MuonSeedPtExtractor::theBeamSpot
private

Definition at line 49 of file MuonSeedPtExtractor.h.

Referenced by pT_extract(), and setBeamSpot().

ParametersMap MuonSeedPtExtractor::theParametersForCombo
private

Definition at line 47 of file MuonSeedPtExtractor.h.

Referenced by fillParametersForCombo(), and pT_extract().

ScalesMap MuonSeedPtExtractor::theScalesForCombo
private

Definition at line 48 of file MuonSeedPtExtractor.h.

Referenced by fillScalesForCombo(), and getPt().

Generated for CMSSW Reference Manual by doxygen 1.8.5