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 Attributes
PVFitter Class Reference

#include <PVFitter.h>

Public Member Functions

void compressStore ()
 reduce size of primary vertex cache by increasing quality limit More...
 
void dumpTxtFile ()
 
void FitPerBunchCrossing ()
 
reco::BeamSpot getBeamSpot ()
 
std::map< int, reco::BeamSpotgetBeamSpotMap ()
 
int * getFitLSRange ()
 
int getNPVs ()
 
const std::map< int, int > & getNPVsperBX ()
 
std::vector< BeamSpotFitPVDatagetpvStore ()
 
time_t * getRefTime ()
 
double getWidthX ()
 
double getWidthXerr ()
 
double getWidthY ()
 
double getWidthYerr ()
 
double getWidthZ ()
 
double getWidthZerr ()
 
void initialize (const edm::ParameterSet &iConfig, edm::ConsumesCollector &iColl)
 
bool IsFitPerBunchCrossing ()
 
 PVFitter ()
 
 PVFitter (const edm::ParameterSet &iConfig, edm::ConsumesCollector &&iColl)
 
 PVFitter (const edm::ParameterSet &iConfig, edm::ConsumesCollector &iColl)
 
double pvQuality (const reco::Vertex &pv) const
 vertex quality measure More...
 
double pvQuality (const BeamSpotFitPVData &pv) const
 vertex quality measure More...
 
void readEvent (const edm::Event &iEvent)
 
void resetAll ()
 
void resetLSRange ()
 
void resetRefTime ()
 
void resizepvStore (unsigned int rmSize)
 
bool runBXFitter ()
 
bool runFitter ()
 
void setFitLSRange (int ls0, int ls1)
 
void setRefTime (std::time_t t0, std::time_t t1)
 
void setTree (TTree *tree)
 
virtual ~PVFitter ()
 

Private Attributes

std::map< int, std::vector
< BeamSpotFitPVData > > 
bxMap_
 
bool debug_
 
bool do3DFit_
 
double dynamicQualityCut_
 
double errorScale_
 
std::ofstream fasciiFile
 
reco::BeamSpot fbeamspot
 
int fbeginLumiOfFit
 
std::map< int, reco::BeamSpotfbspotMap
 
int fendLumiOfFit
 
bool fFitPerBunchCrossing
 
std::time_t freftime [2]
 
TTree * ftree_
 
double fwidthX
 
double fwidthXerr
 
double fwidthY
 
double fwidthYerr
 
double fwidthZ
 
double fwidthZerr
 
TH2F * hPVx
 
TH2F * hPVy
 
unsigned int maxNrVertices_
 
double maxVtxNormChi2_
 
double maxVtxR_
 
double maxVtxZ_
 
unsigned int minNrVertices_
 
double minVtxNdf_
 
unsigned int minVtxTracks_
 
double minVtxWgt_
 
std::map< int, int > npvsmap_
 
std::string outputTxt_
 
std::vector< double > pvQualities_
 
std::vector< BeamSpotFitPVDatapvStore_
 
double sigmaCut_
 
BeamSpotTreeData theBeamSpotTreeData_
 
edm::EDGetTokenT
< reco::VertexCollection
vertexToken_
 

Detailed Description

Definition at line 43 of file PVFitter.h.

Constructor & Destructor Documentation

PVFitter::PVFitter ( )
inline

Definition at line 45 of file PVFitter.h.

45 {}
PVFitter::PVFitter ( const edm::ParameterSet iConfig,
edm::ConsumesCollector &&  iColl 
)

class: PVFitter.cc package: RecoVertex/BeamSpotProducer

author: Francisco Yumiceva, Fermilab (yumic.nosp@m.eva@.nosp@m.fnal..nosp@m.gov) Geng-Yuan Jeng, UC Riverside (Geng-.nosp@m.Yuan.nosp@m..Jeng.nosp@m.@cer.nosp@m.n.ch)


Definition at line 51 of file PVFitter.cc.

References initialize().

53  : ftree_(0)
54 {
55  initialize(iConfig, iColl);
56 }
void initialize(const edm::ParameterSet &iConfig, edm::ConsumesCollector &iColl)
Definition: PVFitter.cc:65
TTree * ftree_
Definition: PVFitter.h:168
PVFitter::PVFitter ( const edm::ParameterSet iConfig,
edm::ConsumesCollector iColl 
)

Definition at line 58 of file PVFitter.cc.

References initialize().

60  :ftree_(0)
61 {
62  initialize(iConfig, iColl);
63 }
void initialize(const edm::ParameterSet &iConfig, edm::ConsumesCollector &iColl)
Definition: PVFitter.cc:65
TTree * ftree_
Definition: PVFitter.h:168
PVFitter::~PVFitter ( )
virtual

Definition at line 96 of file PVFitter.cc.

96  {
97 
98 }

Member Function Documentation

void PVFitter::compressStore ( )

reduce size of primary vertex cache by increasing quality limit

Definition at line 523 of file PVFitter.cc.

References dynamicQualityCut_, i, pvQualities_, pvQuality(), pvStore_, and python.multivaluedict::sort().

Referenced by readEvent().

524 {
525  //
526  // fill vertex qualities
527  //
528  pvQualities_.resize(pvStore_.size());
529  for ( unsigned int i=0; i<pvStore_.size(); ++i ) pvQualities_[i] = pvQuality(pvStore_[i]);
530  sort(pvQualities_.begin(),pvQualities_.end());
531  //
532  // Set new quality cut to median. This cut will be used to reduce the
533  // number of vertices in the store and also apply to all new vertices
534  // until the next reset
535  //
537  //
538  // remove all vertices failing the cut from the store
539  // (to be moved to a more efficient memory management!)
540  //
541  unsigned int iwrite(0);
542  for ( unsigned int i=0; i<pvStore_.size(); ++i ) {
543  if ( pvQuality(pvStore_[i])>dynamicQualityCut_ ) continue;
544  if ( i!=iwrite ) pvStore_[iwrite] = pvStore_[i];
545  ++iwrite;
546  }
547  pvStore_.resize(iwrite);
548  edm::LogInfo("PVFitter") << "Reduced primary vertex store size to "
549  << pvStore_.size() << " ; new dynamic quality cut = "
550  << dynamicQualityCut_ << std::endl;
551 
552 }
std::vector< double > pvQualities_
Definition: PVFitter.h:205
int i
Definition: DBlmapReader.cc:9
double dynamicQualityCut_
Definition: PVFitter.h:204
std::vector< BeamSpotFitPVData > pvStore_
Definition: PVFitter.h:202
double pvQuality(const reco::Vertex &pv) const
vertex quality measure
Definition: PVFitter.cc:555
void PVFitter::dumpTxtFile ( )

Definition at line 480 of file PVFitter.cc.

480  {
481 /*
482  fasciiFile << "Runnumber " << frun << std::endl;
483  fasciiFile << "BeginTimeOfFit " << fbeginTimeOfFit << std::endl;
484  fasciiFile << "EndTimeOfFit " << fendTimeOfFit << std::endl;
485  fasciiFile << "LumiRange " << fbeginLumiOfFit << " - " << fendLumiOfFit << std::endl;
486  fasciiFile << "Type " << fbeamspot.type() << std::endl;
487  fasciiFile << "X0 " << fbeamspot.x0() << std::endl;
488  fasciiFile << "Y0 " << fbeamspot.y0() << std::endl;
489  fasciiFile << "Z0 " << fbeamspot.z0() << std::endl;
490  fasciiFile << "sigmaZ0 " << fbeamspot.sigmaZ() << std::endl;
491  fasciiFile << "dxdz " << fbeamspot.dxdz() << std::endl;
492  fasciiFile << "dydz " << fbeamspot.dydz() << std::endl;
493  if (inputBeamWidth_ > 0 ) {
494  fasciiFile << "BeamWidthX " << inputBeamWidth_ << std::endl;
495  fasciiFile << "BeamWidthY " << inputBeamWidth_ << std::endl;
496  } else {
497  fasciiFile << "BeamWidthX " << fbeamspot.BeamWidthX() << std::endl;
498  fasciiFile << "BeamWidthY " << fbeamspot.BeamWidthY() << std::endl;
499  }
500 
501  for (int i = 0; i<6; ++i) {
502  fasciiFile << "Cov("<<i<<",j) ";
503  for (int j=0; j<7; ++j) {
504  fasciiFile << fbeamspot.covariance(i,j) << " ";
505  }
506  fasciiFile << std::endl;
507  }
508  // beam width error
509  if (inputBeamWidth_ > 0 ) {
510  fasciiFile << "Cov(6,j) 0 0 0 0 0 0 " << "1e-4" << std::endl;
511  } else {
512  fasciiFile << "Cov(6,j) 0 0 0 0 0 0 " << fbeamspot.covariance(6,6) << std::endl;
513  }
514  fasciiFile << "EmittanceX " << fbeamspot.emittanceX() << std::endl;
515  fasciiFile << "EmittanceY " << fbeamspot.emittanceY() << std::endl;
516  fasciiFile << "BetaStar " << fbeamspot.betaStar() << std::endl;
517 
518 */
519 }
void PVFitter::FitPerBunchCrossing ( )
inline

Definition at line 63 of file PVFitter.h.

References fFitPerBunchCrossing.

63 { fFitPerBunchCrossing = true; }
bool fFitPerBunchCrossing
Definition: PVFitter.h:141
reco::BeamSpot PVFitter::getBeamSpot ( )
inline

Definition at line 91 of file PVFitter.h.

References fbeamspot.

Referenced by AlcaBeamMonitor::endLuminosityBlock(), and BeamFitter::runPVandTrkFitter().

91 { return fbeamspot; }
reco::BeamSpot fbeamspot
Definition: PVFitter.h:139
std::map<int, reco::BeamSpot> PVFitter::getBeamSpotMap ( )
inline

Definition at line 92 of file PVFitter.h.

References fbspotMap.

Referenced by BeamFitter::runPVandTrkFitter().

92 { return fbspotMap; }
std::map< int, reco::BeamSpot > fbspotMap
Definition: PVFitter.h:140
int* PVFitter::getFitLSRange ( )
inline

Definition at line 94 of file PVFitter.h.

References fbeginLumiOfFit, fendLumiOfFit, and tmp.

94  {
95  int *tmp=new int[2];
96  tmp[0] = fbeginLumiOfFit;
97  tmp[1] = fendLumiOfFit;
98  return tmp;
99  }
int fbeginLumiOfFit
Definition: PVFitter.h:178
int fendLumiOfFit
Definition: PVFitter.h:179
std::vector< std::vector< double > > tmp
Definition: MVATrainer.cc:100
int PVFitter::getNPVs ( )
inline

Definition at line 120 of file PVFitter.h.

References pvStore_.

Referenced by BeamFitter::getNPVs().

120 { return pvStore_.size(); }
std::vector< BeamSpotFitPVData > pvStore_
Definition: PVFitter.h:202
const std::map<int, int>& PVFitter::getNPVsperBX ( )
inline

Definition at line 122 of file PVFitter.h.

References bxMap_, python.multivaluedict::map(), npvsmap_, and findQualityFiles::size.

Referenced by BeamFitter::getNPVsperBX().

122  {
123 
124  npvsmap_.clear();
125 
126  for ( std::map<int,std::vector<BeamSpotFitPVData> >::const_iterator pvStore = bxMap_.begin();
127  pvStore!=bxMap_.end(); ++pvStore) {
128 
129  //std::cout << "bx " << pvStore->first << " NPVs = " << (pvStore->second).size() << std::endl;
130  npvsmap_[ pvStore->first ] = (pvStore->second).size();
131 
132  }
133  return npvsmap_;
134  }
std::map< int, std::vector< BeamSpotFitPVData > > bxMap_
Definition: PVFitter.h:203
std::map< int, int > npvsmap_
Definition: PVFitter.h:138
tuple size
Write out results.
std::vector<BeamSpotFitPVData> PVFitter::getpvStore ( )
inline

Definition at line 61 of file PVFitter.h.

References pvStore_.

Referenced by BeamFitter::getPVvectorSize().

61 { return pvStore_; }
std::vector< BeamSpotFitPVData > pvStore_
Definition: PVFitter.h:202
time_t* PVFitter::getRefTime ( )
inline

Definition at line 101 of file PVFitter.h.

References freftime.

101  {
102  time_t *tmptime=new time_t[2];
103  tmptime[0]=freftime[0];
104  tmptime[1]=freftime[1];
105  return tmptime;
106  }
std::time_t freftime[2]
Definition: PVFitter.h:164
double PVFitter::getWidthX ( )
inline

Definition at line 54 of file PVFitter.h.

References fwidthX.

54 { return fwidthX; }
double fwidthX
Definition: PVFitter.h:182
double PVFitter::getWidthXerr ( )
inline

Definition at line 57 of file PVFitter.h.

References fwidthXerr.

Referenced by BeamFitter::runPVandTrkFitter().

57 { return fwidthXerr; }
double fwidthXerr
Definition: PVFitter.h:185
double PVFitter::getWidthY ( )
inline

Definition at line 55 of file PVFitter.h.

References fwidthY.

55 { return fwidthY; }
double fwidthY
Definition: PVFitter.h:183
double PVFitter::getWidthYerr ( )
inline

Definition at line 58 of file PVFitter.h.

References fwidthYerr.

58 { return fwidthYerr; }
double fwidthYerr
Definition: PVFitter.h:186
double PVFitter::getWidthZ ( )
inline

Definition at line 56 of file PVFitter.h.

References fwidthZ.

56 { return fwidthZ; }
double fwidthZ
Definition: PVFitter.h:184
double PVFitter::getWidthZerr ( )
inline

Definition at line 59 of file PVFitter.h.

References fwidthZerr.

59 { return fwidthZerr; }
double fwidthZerr
Definition: PVFitter.h:187
void PVFitter::initialize ( const edm::ParameterSet iConfig,
edm::ConsumesCollector iColl 
)

Definition at line 65 of file PVFitter.cc.

References edm::ConsumesCollector::consumes(), debug_, do3DFit_, dynamicQualityCut_, errorScale_, fFitPerBunchCrossing, edm::ParameterSet::getParameter(), hPVx, hPVy, maxNrVertices_, maxVtxNormChi2_, maxVtxR_, maxVtxZ_, minNrVertices_, minVtxNdf_, minVtxTracks_, minVtxWgt_, sigmaCut_, and vertexToken_.

Referenced by PVFitter().

67 {
68  debug_ = iConfig.getParameter<edm::ParameterSet>("PVFitter").getUntrackedParameter<bool>("Debug");
70  iConfig.getParameter<edm::ParameterSet>("PVFitter")
71  .getUntrackedParameter<edm::InputTag>("VertexCollection",
72  edm::InputTag("offlinePrimaryVertices")));
73  do3DFit_ = iConfig.getParameter<edm::ParameterSet>("PVFitter").getUntrackedParameter<bool>("Apply3DFit");
74  //writeTxt_ = iConfig.getParameter<edm::ParameterSet>("PVFitter").getUntrackedParameter<bool>("WriteAscii");
75  //outputTxt_ = iConfig.getParameter<edm::ParameterSet>("PVFitter").getUntrackedParameter<std::string>("AsciiFileName");
76 
77  maxNrVertices_ = iConfig.getParameter<edm::ParameterSet>("PVFitter").getUntrackedParameter<unsigned int>("maxNrStoredVertices");
78  minNrVertices_ = iConfig.getParameter<edm::ParameterSet>("PVFitter").getUntrackedParameter<unsigned int>("minNrVerticesForFit");
79  minVtxNdf_ = iConfig.getParameter<edm::ParameterSet>("PVFitter").getUntrackedParameter<double>("minVertexNdf");
80  maxVtxNormChi2_ = iConfig.getParameter<edm::ParameterSet>("PVFitter").getUntrackedParameter<double>("maxVertexNormChi2");
81  minVtxTracks_ = iConfig.getParameter<edm::ParameterSet>("PVFitter").getUntrackedParameter<unsigned int>("minVertexNTracks");
82  minVtxWgt_ = iConfig.getParameter<edm::ParameterSet>("PVFitter").getUntrackedParameter<double>("minVertexMeanWeight");
83  maxVtxR_ = iConfig.getParameter<edm::ParameterSet>("PVFitter").getUntrackedParameter<double>("maxVertexR");
84  maxVtxZ_ = iConfig.getParameter<edm::ParameterSet>("PVFitter").getUntrackedParameter<double>("maxVertexZ");
85  errorScale_ = iConfig.getParameter<edm::ParameterSet>("PVFitter").getUntrackedParameter<double>("errorScale");
86  sigmaCut_ = iConfig.getParameter<edm::ParameterSet>("PVFitter").getUntrackedParameter<double>("nSigmaCut");
87  fFitPerBunchCrossing=iConfig.getParameter<edm::ParameterSet>("PVFitter").getUntrackedParameter<bool>("FitPerBunchCrossing");
88 
89  // preset quality cut to "infinite"
90  dynamicQualityCut_ = 1.e30;
91 
92  hPVx = new TH2F("hPVx","PVx vs PVz distribution",200,-maxVtxR_, maxVtxR_, 200, -maxVtxZ_, maxVtxZ_);
93  hPVy = new TH2F("hPVy","PVy vs PVz distribution",200,-maxVtxR_, maxVtxR_, 200, -maxVtxZ_, maxVtxZ_);
94 }
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
T getParameter(std::string const &) const
bool debug_
Definition: PVFitter.h:145
double maxVtxNormChi2_
Definition: PVFitter.h:154
double minVtxWgt_
Definition: PVFitter.h:156
double minVtxNdf_
Definition: PVFitter.h:153
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9
double errorScale_
Definition: PVFitter.h:159
bool do3DFit_
Definition: PVFitter.h:146
double dynamicQualityCut_
Definition: PVFitter.h:204
bool fFitPerBunchCrossing
Definition: PVFitter.h:141
TH2F * hPVy
Definition: PVFitter.h:166
edm::EDGetTokenT< reco::VertexCollection > vertexToken_
Definition: PVFitter.h:147
double maxVtxR_
Definition: PVFitter.h:157
double maxVtxZ_
Definition: PVFitter.h:158
unsigned int maxNrVertices_
Definition: PVFitter.h:151
unsigned int minVtxTracks_
Definition: PVFitter.h:155
TH2F * hPVx
Definition: PVFitter.h:166
unsigned int minNrVertices_
Definition: PVFitter.h:152
double sigmaCut_
Definition: PVFitter.h:160
bool PVFitter::IsFitPerBunchCrossing ( )
inline

Definition at line 93 of file PVFitter.h.

References fFitPerBunchCrossing.

Referenced by BeamFitter::dumpTxtFile(), and BeamFitter::runPVandTrkFitter().

93 { return fFitPerBunchCrossing; }
bool fFitPerBunchCrossing
Definition: PVFitter.h:141
double PVFitter::pvQuality ( const reco::Vertex pv) const

vertex quality measure

Definition at line 555 of file PVFitter.cc.

References reco::Vertex::covariance().

Referenced by compressStore(), and readEvent().

556 {
557  //
558  // determinant of the transverse part of the PV covariance matrix
559  //
560  return
561  pv.covariance(0,0)*pv.covariance(1,1)-
562  pv.covariance(0,1)*pv.covariance(0,1);
563 }
double covariance(int i, int j) const
(i, j)-th element of error matrix, i, j = 0, ... 2
Definition: Vertex.h:123
double PVFitter::pvQuality ( const BeamSpotFitPVData pv) const

vertex quality measure

Definition at line 566 of file PVFitter.cc.

References BeamSpotFitPVData::posCorr, and BeamSpotFitPVData::posError.

567 {
568  //
569  // determinant of the transverse part of the PV covariance matrix
570  //
571  double ex = pv.posError[0];
572  double ey = pv.posError[1];
573  return ex*ex*ey*ey*(1-pv.posCorr[0]*pv.posCorr[0]);
574 }
void PVFitter::readEvent ( const edm::Event iEvent)

Definition at line 101 of file PVFitter.cc.

References BeamSpotFitPVData::bunchCrossing, BeamSpotTreeData::bunchCrossing(), edm::EventBase::bunchCrossing(), rpc-layouts::bx, bxMap_, compressStore(), dynamicQualityCut_, fFitPerBunchCrossing, ftree_, edm::Event::getByToken(), hPVx, hPVy, edm::EventBase::id(), BeamSpotTreeData::lumi(), edm::EventBase::luminosityBlock(), maxNrVertices_, minVtxNdf_, minVtxWgt_, BeamSpotFitPVData::posCorr, BeamSpotFitPVData::posError, BeamSpotFitPVData::position, BeamSpotTreeData::pvData(), pvQuality(), pvStore_, BeamSpotTreeData::run(), edm::EventID::run(), theBeamSpotTreeData_, and vertexToken_.

Referenced by AlcaBeamMonitor::analyze(), and BeamFitter::readEvent().

102 {
103 
104 // frun = iEvent.id().run();
105 // const edm::TimeValue_t ftimestamp = iEvent.time().value();
106 // const std::time_t ftmptime = ftimestamp >> 32;
107 
108 // if (fbeginLumiOfFit == -1) freftime[0] = freftime[1] = ftmptime;
109 // if (freftime[0] == 0 || ftmptime < freftime[0]) freftime[0] = ftmptime;
110 // const char* fbeginTime = formatTime(freftime[0]);
111 // sprintf(fbeginTimeOfFit,"%s",fbeginTime);
112 
113 // if (freftime[1] == 0 || ftmptime > freftime[1]) freftime[1] = ftmptime;
114 // const char* fendTime = formatTime(freftime[1]);
115 // sprintf(fendTimeOfFit,"%s",fendTime);
116 
117 // flumi = iEvent.luminosityBlock();
118 // frunFit = frun;
119 
120 // if (fbeginLumiOfFit == -1 || fbeginLumiOfFit > flumi) fbeginLumiOfFit = flumi;
121 // if (fendLumiOfFit == -1 || fendLumiOfFit < flumi) fendLumiOfFit = flumi;
122 // std::cout << "flumi = " <<flumi<<"; fbeginLumiOfFit = " << fbeginLumiOfFit <<"; fendLumiOfFit = "<<fendLumiOfFit<<std::endl;
123 
124  //------ Primary Vertices
126  bool hasPVs = false;
127  //edm::View<reco::Vertex> vertices;
128  //const reco::VertexCollection & vertices = 0;
129 
130  if ( iEvent.getByToken(vertexToken_, PVCollection ) ) {
131  //pv = *PVCollection;
132  //vertices = *PVCollection;
133  hasPVs = true;
134  }
135  //------
136 
137  if ( hasPVs ) {
138 
139  for (reco::VertexCollection::const_iterator pv = PVCollection->begin(); pv != PVCollection->end(); ++pv ) {
140 
141 
142  //for ( size_t ipv=0; ipv != pv.size(); ++ipv ) {
143 
144  //--- vertex selection
145  if ( pv->isFake() || pv->tracksSize()==0 ) continue;
146  if ( pv->ndof() < minVtxNdf_ || (pv->ndof()+3.)/pv->tracksSize()<2*minVtxWgt_ ) continue;
147  //---
148 
149  hPVx->Fill( pv->x(), pv->z() );
150  hPVy->Fill( pv->y(), pv->z() );
151 
152  //
153  // 3D fit section
154  //
155  // apply additional quality cut
156  if ( pvQuality(*pv)>dynamicQualityCut_ ) continue;
157  // if store exceeds max. size: reduce size and apply new quality cut
158  if ( pvStore_.size()>=maxNrVertices_ ) {
159  compressStore();
160  if ( pvQuality(*pv)>dynamicQualityCut_ ) continue;
161  }
162  //
163  // copy PV to store
164  //
165  int bx = iEvent.bunchCrossing();
166  BeamSpotFitPVData pvData;
167  pvData.bunchCrossing = bx;
168  pvData.position[0] = pv->x();
169  pvData.position[1] = pv->y();
170  pvData.position[2] = pv->z();
171  pvData.posError[0] = pv->xError();
172  pvData.posError[1] = pv->yError();
173  pvData.posError[2] = pv->zError();
174  pvData.posCorr[0] = pv->covariance(0,1)/pv->xError()/pv->yError();
175  pvData.posCorr[1] = pv->covariance(0,2)/pv->xError()/pv->zError();
176  pvData.posCorr[2] = pv->covariance(1,2)/pv->yError()/pv->zError();
177  pvStore_.push_back(pvData);
178 
179  if(ftree_ != 0){
180  theBeamSpotTreeData_.run(iEvent.id().run());
184  ftree_->Fill();
185  }
186 
187  if (fFitPerBunchCrossing) bxMap_[bx].push_back(pvData);
188 
189  }
190 
191  }
192 
193 
194 
195 
196 }
RunNumber_t run() const
Definition: EventID.h:42
std::map< int, std::vector< BeamSpotFitPVData > > bxMap_
Definition: PVFitter.h:203
void run(unsigned int run)
void compressStore()
reduce size of primary vertex cache by increasing quality limit
Definition: PVFitter.cc:523
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:434
double minVtxWgt_
Definition: PVFitter.h:156
double minVtxNdf_
Definition: PVFitter.h:153
int bunchCrossing() const
Definition: EventBase.h:62
edm::LuminosityBlockNumber_t luminosityBlock() const
Definition: EventBase.h:59
double dynamicQualityCut_
Definition: PVFitter.h:204
void bunchCrossing(unsigned int bunchCrossing)
void lumi(unsigned int lumi)
std::vector< BeamSpotFitPVData > pvStore_
Definition: PVFitter.h:202
BeamSpotTreeData theBeamSpotTreeData_
Definition: PVFitter.h:207
bool fFitPerBunchCrossing
Definition: PVFitter.h:141
TH2F * hPVy
Definition: PVFitter.h:166
void pvData(const BeamSpotFitPVData &pvData)
TTree * ftree_
Definition: PVFitter.h:168
edm::EDGetTokenT< reco::VertexCollection > vertexToken_
Definition: PVFitter.h:147
double pvQuality(const reco::Vertex &pv) const
vertex quality measure
Definition: PVFitter.cc:555
unsigned int maxNrVertices_
Definition: PVFitter.h:151
TH2F * hPVx
Definition: PVFitter.h:166
string bx
Definition: rpc-layouts.py:13
edm::EventID id() const
Definition: EventBase.h:56
void PVFitter::resetAll ( )
inline

Definition at line 80 of file PVFitter.h.

References align::BeamSpot, bxMap_, dynamicQualityCut_, fbeamspot, fbspotMap, hPVx, hPVy, pvStore_, resetLSRange(), and resetRefTime().

Referenced by BeamFitter::BeamFitter(), AlcaBeamMonitor::endLuminosityBlock(), and BeamFitter::resetPVFitter().

80  {
81  resetLSRange();
82  resetRefTime();
83  pvStore_.clear();
84  bxMap_.clear();
85  dynamicQualityCut_ = 1.e30;
86  hPVx->Reset();
87  hPVy->Reset();
89  fbspotMap.clear();
90  };
std::map< int, std::vector< BeamSpotFitPVData > > bxMap_
Definition: PVFitter.h:203
double dynamicQualityCut_
Definition: PVFitter.h:204
reco::BeamSpot fbeamspot
Definition: PVFitter.h:139
std::vector< BeamSpotFitPVData > pvStore_
Definition: PVFitter.h:202
void resetLSRange()
Definition: PVFitter.h:66
TH2F * hPVy
Definition: PVFitter.h:166
std::map< int, reco::BeamSpot > fbspotMap
Definition: PVFitter.h:140
void resetRefTime()
Definition: PVFitter.h:67
TH2F * hPVx
Definition: PVFitter.h:166
void PVFitter::resetLSRange ( )
inline

Definition at line 66 of file PVFitter.h.

References fbeginLumiOfFit, and fendLumiOfFit.

Referenced by resetAll().

int fbeginLumiOfFit
Definition: PVFitter.h:178
int fendLumiOfFit
Definition: PVFitter.h:179
void PVFitter::resetRefTime ( )
inline

Definition at line 67 of file PVFitter.h.

References freftime.

Referenced by resetAll().

67 { freftime[0] = freftime[1] = 0; }
std::time_t freftime[2]
Definition: PVFitter.h:164
void PVFitter::resizepvStore ( unsigned int  rmSize)
inline

Definition at line 109 of file PVFitter.h.

References pvStore_.

Referenced by BeamFitter::resizePVvector().

109  {
110  pvStore_.erase(pvStore_.begin(), pvStore_.begin()+rmSize);
111  }
std::vector< BeamSpotFitPVData > pvStore_
Definition: PVFitter.h:202
bool PVFitter::runBXFitter ( )

Definition at line 203 of file PVFitter.cc.

References align::BeamSpot, bxMap_, errorScale_, fbeamspot, fbspotMap, fcn(), fwidthX, fwidthXerr, fwidthY, fwidthYerr, fwidthZ, fwidthZerr, python.multivaluedict::map(), makeMuonMisalignmentScenario::matrix, minNrVertices_, funct::pow(), reco::BeamSpot::setBeamWidthX(), reco::BeamSpot::setBeamWidthY(), FcnBeamSpotFitPV::setLimits(), reco::BeamSpot::setType(), sigmaCut_, findQualityFiles::size, and reco::BeamSpot::Tracker.

Referenced by BeamFitter::runPVandTrkFitter().

203  {
204 
205  using namespace ROOT::Minuit2;
206 
207  edm::LogInfo("PVFitter") << " Number of bunch crossings: " << bxMap_.size() << std::endl;
208 
209  bool fit_ok = true;
210 
211  for ( std::map<int,std::vector<BeamSpotFitPVData> >::const_iterator pvStore = bxMap_.begin();
212  pvStore!=bxMap_.end(); ++pvStore) {
213 
214  // first set null beam spot in case
215  // fit fails
216  fbspotMap[pvStore->first] = reco::BeamSpot();
217 
218  edm::LogInfo("PVFitter") << " Number of PVs collected for PVFitter: " << (pvStore->second).size() << " in bx: " << pvStore->first << std::endl;
219 
220  if ( (pvStore->second).size() <= minNrVertices_ ) {
221  edm::LogWarning("PVFitter") << " not enough PVs, continue" << std::endl;
222  fit_ok = false;
223  continue;
224  }
225 
226  //bool fit_ok = false;
227  edm::LogInfo("PVFitter") << "Calculating beam spot with PVs ..." << std::endl;
228 
229  //
230  // LL function and fitter
231  //
232  FcnBeamSpotFitPV* fcn = new FcnBeamSpotFitPV(pvStore->second);
233  //
234  // fit parameters: positions, widths, x-y correlations, tilts in xz and yz
235  //
236  MnUserParameters upar;
237  upar.Add("x", 0., 0.02, -10., 10.); // 0
238  upar.Add("y", 0., 0.02, -10., 10.); // 1
239  upar.Add("z", 0., 0.20, -30., 30.); // 2
240  upar.Add("ex", 0.015, 0.01, 0., 10.); // 3
241  upar.Add("corrxy", 0., 0.02, -1., 1.); // 4
242  upar.Add("ey", 0.015, 0.01, 0., 10.); // 5
243  upar.Add("dxdz", 0., 0.0002, -0.1, 0.1); // 6
244  upar.Add("dydz", 0., 0.0002, -0.1, 0.1); // 7
245  upar.Add("ez", 1., 0.1, 0., 30.); // 8
246  upar.Add("scale", errorScale_, errorScale_/10.,
247  errorScale_/2., errorScale_*2.); // 9
248  MnMigrad migrad(*fcn, upar);
249 
250  //
251  // first iteration without correlations
252  //
253  upar.Fix(4);
254  upar.Fix(6);
255  upar.Fix(7);
256  upar.Fix(9);
257  FunctionMinimum ierr = migrad();
258  if ( !ierr.IsValid() ) {
259  edm::LogInfo("PVFitter") << "3D beam spot fit failed in 1st iteration" << std::endl;
260  fit_ok = false;
261  continue;
262  }
263  //
264  // refit with harder selection on vertices
265  //
266  fcn->setLimits(upar.Value(0)-sigmaCut_*upar.Value(3),
267  upar.Value(0)+sigmaCut_*upar.Value(3),
268  upar.Value(1)-sigmaCut_*upar.Value(5),
269  upar.Value(1)+sigmaCut_*upar.Value(5),
270  upar.Value(2)-sigmaCut_*upar.Value(8),
271  upar.Value(2)+sigmaCut_*upar.Value(8));
272  ierr = migrad();
273  if ( !ierr.IsValid() ) {
274  edm::LogInfo("PVFitter") << "3D beam spot fit failed in 2nd iteration" << std::endl;
275  fit_ok = false;
276  continue;
277  }
278  //
279  // refit with correlations
280  //
281  upar.Release(4);
282  upar.Release(6);
283  upar.Release(7);
284  ierr = migrad();
285  if ( !ierr.IsValid() ) {
286  edm::LogInfo("PVFitter") << "3D beam spot fit failed in 3rd iteration" << std::endl;
287  fit_ok = false;
288  continue;
289  }
290  // refit with floating scale factor
291  // minuitx.ReleaseParameter(9);
292  // minuitx.Minimize();
293 
294  //minuitx.PrintResults(0,0);
295 
296  fwidthX = upar.Value(3);
297  fwidthY = upar.Value(5);
298  fwidthZ = upar.Value(8);
299  fwidthXerr = upar.Error(3);
300  fwidthYerr = upar.Error(5);
301  fwidthZerr = upar.Error(8);
302 
304  // need to get the full cov matrix
305  matrix(0,0) = pow( upar.Error(0), 2);
306  matrix(1,1) = pow( upar.Error(1), 2);
307  matrix(2,2) = pow( upar.Error(2), 2);
308  matrix(3,3) = fwidthZerr * fwidthZerr;
309  matrix(4,4) = pow( upar.Error(6), 2);
310  matrix(5,5) = pow( upar.Error(7), 2);
311  matrix(6,6) = fwidthXerr * fwidthXerr;
312 
314  upar.Value(1),
315  upar.Value(2) ),
316  fwidthZ,
317  upar.Value(6), upar.Value(7),
318  fwidthX,
319  matrix );
323 
324  fbspotMap[pvStore->first] = fbeamspot;
325  edm::LogInfo("PVFitter") << "3D PV fit done for this bunch crossing."<<std::endl;
326  //delete fcn;
327  fit_ok = fit_ok & true;
328  }
329 
330  return fit_ok;
331 }
math::Error< dimension >::type CovarianceMatrix
Definition: BeamSpot.h:31
std::map< int, std::vector< BeamSpotFitPVData > > bxMap_
Definition: PVFitter.h:203
math::XYZPoint Point
point in the space
Definition: BeamSpot.h:29
double fwidthY
Definition: PVFitter.h:183
double errorScale_
Definition: PVFitter.h:159
double fwidthYerr
Definition: PVFitter.h:186
double fwidthZ
Definition: PVFitter.h:184
double fwidthXerr
Definition: PVFitter.h:185
void setType(BeamType type)
set beam type
Definition: BeamSpot.h:131
reco::BeamSpot fbeamspot
Definition: PVFitter.h:139
void setBeamWidthY(double v)
Definition: BeamSpot.h:109
std::map< int, reco::BeamSpot > fbspotMap
Definition: PVFitter.h:140
double fwidthZerr
Definition: PVFitter.h:187
void fcn(int &, double *, double &, double *, int)
unsigned int minNrVertices_
Definition: PVFitter.h:152
void setLimits(float xmin, float xmax, float ymin, float ymax, float zmin, float zmax)
double sigmaCut_
Definition: PVFitter.h:160
void setBeamWidthX(double v)
Definition: BeamSpot.h:108
double fwidthX
Definition: PVFitter.h:182
tuple size
Write out results.
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
bool PVFitter::runFitter ( )

Definition at line 334 of file PVFitter.cc.

References align::BeamSpot, do3DFit_, errorScale_, fbeamspot, fcn(), fwidthX, fwidthXerr, fwidthY, fwidthYerr, fwidthZ, fwidthZerr, hPVx, hPVy, edm::isNotFinite(), makeMuonMisalignmentScenario::matrix, minNrVertices_, funct::pow(), pvStore_, reco::BeamSpot::setBeamWidthX(), reco::BeamSpot::setBeamWidthY(), FcnBeamSpotFitPV::setLimits(), reco::BeamSpot::setType(), sigmaCut_, and reco::BeamSpot::Tracker.

Referenced by AlcaBeamMonitor::endLuminosityBlock(), and BeamFitter::runPVandTrkFitter().

334  {
335 
336  using namespace ROOT::Minuit2;
337  edm::LogInfo("PVFitter") << " Number of PVs collected for PVFitter: " << pvStore_.size() << std::endl;
338 
339  if ( pvStore_.size() <= minNrVertices_ ) return false;
340 
341  //bool fit_ok = false;
342 
343  if ( ! do3DFit_ ) {
344  TH1F *h1PVx = (TH1F*) hPVx->ProjectionX("h1PVx", 0, -1, "e");
345  TH1F *h1PVy = (TH1F*) hPVy->ProjectionX("h1PVy", 0, -1, "e");
346  TH1F *h1PVz = (TH1F*) hPVx->ProjectionY("h1PVz", 0, -1, "e");
347 
348  h1PVx->Fit("gaus","QLM0");
349  h1PVy->Fit("gaus","QLM0");
350  h1PVz->Fit("gaus","QLM0");
351 
352  TF1 *gausx = h1PVx->GetFunction("gaus");
353  TF1 *gausy = h1PVy->GetFunction("gaus");
354  TF1 *gausz = h1PVz->GetFunction("gaus");
355 
356  fwidthX = gausx->GetParameter(2);
357  fwidthY = gausy->GetParameter(2);
358  fwidthZ = gausz->GetParameter(2);
359  fwidthXerr = gausx->GetParError(2);
360  fwidthYerr = gausy->GetParError(2);
361  fwidthZerr = gausz->GetParError(2);
362 
364  matrix(2,2) = gausz->GetParError(1) * gausz->GetParError(1);
365  matrix(3,3) = fwidthZerr * fwidthZerr;
366  matrix(6,6) = fwidthXerr * fwidthXerr;
367 
368  fbeamspot = reco::BeamSpot( reco::BeamSpot::Point(gausx->GetParameter(1),
369  gausy->GetParameter(1),
370  gausz->GetParameter(1) ),
371  fwidthZ,
372  0., 0.,
373  fwidthX,
374  matrix );
378 
379  }
380  else { // do 3D fit
381  //
382  // LL function and fitter
383  //
385  //
386  // fit parameters: positions, widths, x-y correlations, tilts in xz and yz
387  //
388  MnUserParameters upar;
389  upar.Add("x", 0., 0.02, -10., 10.); // 0
390  upar.Add("y", 0., 0.02, -10., 10.); // 1
391  upar.Add("z", 0., 0.20, -30., 30.); // 2
392  upar.Add("ex", 0.005, 0.0005, 0.0001, 0.05); // 3
393  upar.Add("corrxy", 0., 0.02, -1., 1.); // 4
394  upar.Add("ey", 0.005, 0.0005, 0.0001, 0.05); // 5
395  upar.Add("dxdz", 0., 0.0002, -0.1, 0.1); // 6
396  upar.Add("dydz", 0., 0.0002, -0.1, 0.1); // 7
397  upar.Add("ez", 1., 0.1, 1.0, 30.); // 8
398  upar.Add("scale", errorScale_, errorScale_/10.,
399  errorScale_/2., errorScale_*2.); // 9
400  MnMigrad migrad(*fcn, upar);
401  //
402  // first iteration without correlations
403  //
404  upar.Fix(4);
405  upar.Fix(6);
406  upar.Fix(7);
407  upar.Fix(9);
408  FunctionMinimum ierr = migrad();
409  if ( !ierr.IsValid() ) {
410  edm::LogWarning("PVFitter") << "3D beam spot fit failed in 1st iteration" << std::endl;
411  return false;
412  }
413  //
414  // refit with harder selection on vertices
415  //
416  fcn->setLimits(upar.Value(0)-sigmaCut_*upar.Value(3),
417  upar.Value(0)+sigmaCut_*upar.Value(3),
418  upar.Value(1)-sigmaCut_*upar.Value(5),
419  upar.Value(1)+sigmaCut_*upar.Value(5),
420  upar.Value(2)-sigmaCut_*upar.Value(8),
421  upar.Value(2)+sigmaCut_*upar.Value(8));
422  ierr = migrad();
423  if ( !ierr.IsValid() ) {
424  edm::LogWarning("PVFitter") << "3D beam spot fit failed in 2nd iteration" << std::endl;
425  return false;
426  }
427  //
428  // refit with correlations
429  //
430  upar.Release(4);
431  upar.Release(6);
432  upar.Release(7);
433  ierr = migrad();
434  if ( !ierr.IsValid() ) {
435  edm::LogWarning("PVFitter") << "3D beam spot fit failed in 3rd iteration" << std::endl;
436  return false;
437  }
438  // refit with floating scale factor
439  // minuitx.ReleaseParameter(9);
440  // minuitx.Minimize();
441 
442  //minuitx.PrintResults(0,0);
443 
444  fwidthX = upar.Value(3);
445  fwidthY = upar.Value(5);
446  fwidthZ = upar.Value(8);
447  fwidthXerr = upar.Error(3);
448  fwidthYerr = upar.Error(5);
449  fwidthZerr = upar.Error(8);
450 
451  // check errors on widths and sigmaZ for nan
452  if ( edm::isNotFinite(fwidthXerr) || edm::isNotFinite(fwidthYerr) || edm::isNotFinite(fwidthZerr) ) {
453  edm::LogWarning("PVFitter") << "3D beam spot fit returns nan in 3rd iteration" << std::endl;
454  return false;
455  }
456 
458  // need to get the full cov matrix
459  matrix(0,0) = pow( upar.Error(0), 2);
460  matrix(1,1) = pow( upar.Error(1), 2);
461  matrix(2,2) = pow( upar.Error(2), 2);
462  matrix(3,3) = fwidthZerr * fwidthZerr;
463  matrix(6,6) = fwidthXerr * fwidthXerr;
464 
466  upar.Value(1),
467  upar.Value(2) ),
468  fwidthZ,
469  upar.Value(6), upar.Value(7),
470  fwidthX,
471  matrix );
475  }
476 
477  return true; //FIXME: Need to add quality test for the fit results!
478 }
math::Error< dimension >::type CovarianceMatrix
Definition: BeamSpot.h:31
math::XYZPoint Point
point in the space
Definition: BeamSpot.h:29
double fwidthY
Definition: PVFitter.h:183
double errorScale_
Definition: PVFitter.h:159
bool do3DFit_
Definition: PVFitter.h:146
double fwidthYerr
Definition: PVFitter.h:186
double fwidthZ
Definition: PVFitter.h:184
double fwidthXerr
Definition: PVFitter.h:185
void setType(BeamType type)
set beam type
Definition: BeamSpot.h:131
reco::BeamSpot fbeamspot
Definition: PVFitter.h:139
std::vector< BeamSpotFitPVData > pvStore_
Definition: PVFitter.h:202
bool isNotFinite(T x)
Definition: isFinite.h:10
void setBeamWidthY(double v)
Definition: BeamSpot.h:109
TH2F * hPVy
Definition: PVFitter.h:166
double fwidthZerr
Definition: PVFitter.h:187
void fcn(int &, double *, double &, double *, int)
TH2F * hPVx
Definition: PVFitter.h:166
unsigned int minNrVertices_
Definition: PVFitter.h:152
void setLimits(float xmin, float xmax, float ymin, float ymax, float zmin, float zmax)
double sigmaCut_
Definition: PVFitter.h:160
void setBeamWidthX(double v)
Definition: BeamSpot.h:108
double fwidthX
Definition: PVFitter.h:182
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
void PVFitter::setFitLSRange ( int  ls0,
int  ls1 
)
inline

Definition at line 73 of file PVFitter.h.

References fbeginLumiOfFit, and fendLumiOfFit.

73  {
74  fbeginLumiOfFit = ls0;
75  fendLumiOfFit = ls1;
76  }
int fbeginLumiOfFit
Definition: PVFitter.h:178
int fendLumiOfFit
Definition: PVFitter.h:179
void PVFitter::setRefTime ( std::time_t  t0,
std::time_t  t1 
)
inline

Definition at line 69 of file PVFitter.h.

References freftime.

69  {
70  freftime[0] = t0;
71  freftime[1] = t1;
72  }
std::time_t freftime[2]
Definition: PVFitter.h:164
void PVFitter::setTree ( TTree *  tree)

Definition at line 198 of file PVFitter.cc.

References BeamSpotTreeData::branch(), ftree_, theBeamSpotTreeData_, and diffTreeTool::tree.

Referenced by BeamFitter::BeamFitter().

198  {
199  ftree_ = tree;
201 }
BeamSpotTreeData theBeamSpotTreeData_
Definition: PVFitter.h:207
TTree * ftree_
Definition: PVFitter.h:168
void branch(TTree *tree)

Member Data Documentation

std::map< int, std::vector<BeamSpotFitPVData> > PVFitter::bxMap_
private

Definition at line 203 of file PVFitter.h.

Referenced by getNPVsperBX(), readEvent(), resetAll(), and runBXFitter().

bool PVFitter::debug_
private

Definition at line 145 of file PVFitter.h.

Referenced by initialize().

bool PVFitter::do3DFit_
private

Definition at line 146 of file PVFitter.h.

Referenced by initialize(), and runFitter().

double PVFitter::dynamicQualityCut_
private

Definition at line 204 of file PVFitter.h.

Referenced by compressStore(), initialize(), readEvent(), and resetAll().

double PVFitter::errorScale_
private

Definition at line 159 of file PVFitter.h.

Referenced by initialize(), runBXFitter(), and runFitter().

std::ofstream PVFitter::fasciiFile
private

Definition at line 143 of file PVFitter.h.

reco::BeamSpot PVFitter::fbeamspot
private

Definition at line 139 of file PVFitter.h.

Referenced by getBeamSpot(), resetAll(), runBXFitter(), and runFitter().

int PVFitter::fbeginLumiOfFit
private

Definition at line 178 of file PVFitter.h.

Referenced by getFitLSRange(), resetLSRange(), and setFitLSRange().

std::map<int,reco::BeamSpot> PVFitter::fbspotMap
private

Definition at line 140 of file PVFitter.h.

Referenced by getBeamSpotMap(), resetAll(), and runBXFitter().

int PVFitter::fendLumiOfFit
private

Definition at line 179 of file PVFitter.h.

Referenced by getFitLSRange(), resetLSRange(), and setFitLSRange().

bool PVFitter::fFitPerBunchCrossing
private

Definition at line 141 of file PVFitter.h.

Referenced by FitPerBunchCrossing(), initialize(), IsFitPerBunchCrossing(), and readEvent().

std::time_t PVFitter::freftime[2]
private

Definition at line 164 of file PVFitter.h.

Referenced by getRefTime(), resetRefTime(), and setRefTime().

TTree* PVFitter::ftree_
private

Definition at line 168 of file PVFitter.h.

Referenced by readEvent(), and setTree().

double PVFitter::fwidthX
private

Definition at line 182 of file PVFitter.h.

Referenced by getWidthX(), runBXFitter(), and runFitter().

double PVFitter::fwidthXerr
private

Definition at line 185 of file PVFitter.h.

Referenced by getWidthXerr(), runBXFitter(), and runFitter().

double PVFitter::fwidthY
private

Definition at line 183 of file PVFitter.h.

Referenced by getWidthY(), runBXFitter(), and runFitter().

double PVFitter::fwidthYerr
private

Definition at line 186 of file PVFitter.h.

Referenced by getWidthYerr(), runBXFitter(), and runFitter().

double PVFitter::fwidthZ
private

Definition at line 184 of file PVFitter.h.

Referenced by getWidthZ(), runBXFitter(), and runFitter().

double PVFitter::fwidthZerr
private

Definition at line 187 of file PVFitter.h.

Referenced by getWidthZerr(), runBXFitter(), and runFitter().

TH2F* PVFitter::hPVx
private

Definition at line 166 of file PVFitter.h.

Referenced by initialize(), readEvent(), resetAll(), and runFitter().

TH2F* PVFitter::hPVy
private

Definition at line 166 of file PVFitter.h.

Referenced by initialize(), readEvent(), resetAll(), and runFitter().

unsigned int PVFitter::maxNrVertices_
private

Definition at line 151 of file PVFitter.h.

Referenced by initialize(), and readEvent().

double PVFitter::maxVtxNormChi2_
private

Definition at line 154 of file PVFitter.h.

Referenced by initialize().

double PVFitter::maxVtxR_
private

Definition at line 157 of file PVFitter.h.

Referenced by initialize().

double PVFitter::maxVtxZ_
private

Definition at line 158 of file PVFitter.h.

Referenced by initialize().

unsigned int PVFitter::minNrVertices_
private

Definition at line 152 of file PVFitter.h.

Referenced by initialize(), runBXFitter(), and runFitter().

double PVFitter::minVtxNdf_
private

Definition at line 153 of file PVFitter.h.

Referenced by initialize(), and readEvent().

unsigned int PVFitter::minVtxTracks_
private

Definition at line 155 of file PVFitter.h.

Referenced by initialize().

double PVFitter::minVtxWgt_
private

Definition at line 156 of file PVFitter.h.

Referenced by initialize(), and readEvent().

std::map<int, int> PVFitter::npvsmap_
private

Definition at line 138 of file PVFitter.h.

Referenced by getNPVsperBX().

std::string PVFitter::outputTxt_
private

Definition at line 149 of file PVFitter.h.

std::vector<double> PVFitter::pvQualities_
private

Definition at line 205 of file PVFitter.h.

Referenced by compressStore().

std::vector<BeamSpotFitPVData> PVFitter::pvStore_
private
double PVFitter::sigmaCut_
private

Definition at line 160 of file PVFitter.h.

Referenced by initialize(), runBXFitter(), and runFitter().

BeamSpotTreeData PVFitter::theBeamSpotTreeData_
private

Definition at line 207 of file PVFitter.h.

Referenced by readEvent(), and setTree().

edm::EDGetTokenT<reco::VertexCollection> PVFitter::vertexToken_
private

Definition at line 147 of file PVFitter.h.

Referenced by initialize(), and readEvent().