CMS 3D CMS Logo
List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes
BeamFitter Class Reference

#include <BeamFitter.h>

Public Member Functions

 BeamFitter ()
 
 BeamFitter (const edm::ParameterSet &iConfig, edm::ConsumesCollector &&iColl)
 
void dumpBWTxtFile (std::string &)
 
void dumpTxtFile (std::string &, bool)
 
reco::BeamSpot getBeamSpot ()
 
std::map< int, reco::BeamSpot > getBeamSpotMap ()
 
reco::BeamSpot getBeamWidth ()
 
std::vector< BSTrkParameters > getBSvector ()
 
TH1F * getCutFlow ()
 
std::pair< int, int > getFitLSRange ()
 
int getNPVs ()
 
const std::map< int, int > & getNPVsperBX ()
 
int getNTracks ()
 
std::size_t getPVvectorSize ()
 
std::pair< time_t, time_t > getRefTime ()
 
int getRunNumber ()
 
void readEvent (const edm::Event &iEvent)
 
void resetCutFlow ()
 
void resetLSRange ()
 
void resetPVFitter ()
 
void resetRefTime ()
 
void resetTotTrk ()
 
void resetTrkVector ()
 
void resizeBSvector (unsigned int nsize)
 
void resizePVvector (unsigned int npvsize)
 
void runAllFitter ()
 
bool runBeamWidthFitter ()
 
bool runFitter ()
 
bool runFitterNoTxt ()
 
bool runPVandTrkFitter ()
 
void setFitLSRange (int ls0, int ls1)
 
void SetPVInfo (const std::vector< float > &v1_)
 
void setRefTime (time_t t0, time_t t1)
 
void setRun (int run)
 
void subtractFromCutFlow (const TH1F *toSubtract)
 
void write2DB ()
 
virtual ~BeamFitter ()
 

Static Public Member Functions

static void fillDescription (edm::ParameterSetDescription &)
 

Private Member Functions

void updateBTime ()
 

Private Attributes

std::vector< reco::TrackBase::TrackAlgorithm > algorithm_
 
bool appendRunTxt_
 
edm::EDGetTokenT< reco::BeamSpot > beamSpotToken_
 
double convergence_
 
int countPass [9]
 
bool debug_
 
bool falgo
 
std::ofstream fasciiDIP
 
std::ofstream fasciiFile
 
reco::BeamSpot fbeamspot
 
reco::BeamSpot fbeamWidthFit
 
int fbeginLumiOfFit
 
char fbeginTimeOfFit [32]
 
std::map< int, reco::BeamSpot > fbspotPVMap
 
std::vector< BSTrkParameters > fBSvector
 
int fcharge
 
double fcov [7][7]
 
double fd0
 
double fd0bs
 
double fdxdz
 
double fdxdzErr
 
double fdydz
 
double fdydzErr
 
int fendLumiOfFit
 
char fendTimeOfFit [32]
 
double feta
 
bool ffilename_changed
 
TFile * file_
 
bool fitted_
 
int flumi
 
double fnormchi2
 
int fnPixelLayerMeas
 
int fnPXBLayerMeas
 
int fnPXFLayerMeas
 
int fnStripLayerMeas
 
int fnTECLayerMeas
 
int fnTIBLayerMeas
 
int fnTIDLayerMeas
 
int fnTOBLayerMeas
 
int fnTotLayerMeas
 
std::vector< float > ForDIPPV_
 
double fphi0
 
double fpt
 
TTree * fPVTree_
 
bool fpvValid
 
double fpvx
 
double fpvy
 
double fpvz
 
bool fquality
 
std::time_t freftime [2]
 
int frun
 
int frunFit
 
double fsigmad0
 
double fsigmaZ
 
double fsigmaz0
 
double fsigmaZErr
 
double ftheta
 
int ftotal_tracks
 
TTree * ftree_
 
TTree * ftreeFit_
 
double fvx
 
double fvy
 
double fwidthX
 
double fwidthXErr
 
double fwidthY
 
double fwidthYErr
 
double fx
 
double fxErr
 
double fy
 
double fyErr
 
double fz
 
double fz0
 
double fzErr
 
TH1F * h1cutFlow
 
TH1F * h1ntrks
 
TH1F * h1vz_event
 
TH1F * h1z
 
double inputBeamWidth_
 
bool isMuon_
 
int min_Ntrks_
 
PVFitter * MyPVFitter
 
std::string outputDIPTxt_
 
std::string outputfilename_
 
std::string outputTxt_
 
std::vector< reco::TrackBase::TrackQuality > quality_
 
bool saveBeamFit_
 
bool saveNtuple_
 
bool savePVVertices_
 
edm::EDGetTokenT< reco::TrackCollection > tracksToken_
 
std::vector< std::string > trk_Algorithm_
 
double trk_MaxEta_
 
double trk_MaxIP_
 
double trk_MaxNormChi2_
 
double trk_MaxZ_
 
int trk_MinNPixLayers_
 
int trk_MinNTotLayers_
 
double trk_MinpT_
 
std::vector< std::string > trk_Quality_
 
edm::EDGetTokenT< edm::View< reco::Vertex > > vertexToken_
 
bool writeDIPBadFit_
 
bool writeDIPTxt_
 
bool writeTxt_
 

Detailed Description

Definition at line 41 of file BeamFitter.h.

Constructor & Destructor Documentation

◆ BeamFitter() [1/2]

BeamFitter::BeamFitter ( )
inline

Definition at line 43 of file BeamFitter.h.

43 {}

◆ BeamFitter() [2/2]

BeamFitter::BeamFitter ( const edm::ParameterSet &  iConfig,
edm::ConsumesCollector &&  iColl 
)

Definition at line 39 of file BeamFitter.cc.

39  : fPVTree_(nullptr) {
40  debug_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<bool>("Debug");
41  tracksToken_ = iColl.consumes<reco::TrackCollection>(
42  iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<edm::InputTag>("TrackCollection"));
43  vertexToken_ = iColl.consumes<edm::View<reco::Vertex> >(
44  iConfig.getUntrackedParameter<edm::InputTag>("primaryVertex", edm::InputTag("offlinePrimaryVertices")));
45  beamSpotToken_ = iColl.consumes<reco::BeamSpot>(
46  iConfig.getUntrackedParameter<edm::InputTag>("beamSpot", edm::InputTag("offlineBeamSpot")));
47  writeTxt_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<bool>("WriteAscii");
48  outputTxt_ =
49  iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<std::string>("AsciiFileName");
50  appendRunTxt_ =
51  iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<bool>("AppendRunToFileName");
52  writeDIPTxt_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<bool>("WriteDIPAscii");
53  // Specify whether we want to write the DIP file even if the fit is failed.
54  writeDIPBadFit_ =
55  iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<bool>("WriteDIPOnBadFit", true);
56  outputDIPTxt_ =
57  iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<std::string>("DIPFileName");
58  saveNtuple_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<bool>("SaveNtuple");
59  saveBeamFit_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<bool>("SaveFitResults");
60  savePVVertices_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<bool>("SavePVVertices");
61  isMuon_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<bool>("IsMuonCollection");
62 
63  trk_MinpT_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<double>("MinimumPt");
64  trk_MaxEta_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<double>("MaximumEta");
65  trk_MaxIP_ =
66  iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<double>("MaximumImpactParameter");
67  trk_MaxZ_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<double>("MaximumZ");
68  trk_MinNTotLayers_ =
69  iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<int>("MinimumTotalLayers");
70  trk_MinNPixLayers_ =
71  iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<int>("MinimumPixelLayers");
72  trk_MaxNormChi2_ =
73  iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<double>("MaximumNormChi2");
74  trk_Algorithm_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter")
75  .getUntrackedParameter<std::vector<std::string> >("TrackAlgorithm");
76  trk_Quality_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter")
77  .getUntrackedParameter<std::vector<std::string> >("TrackQuality");
78  min_Ntrks_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<int>("MinimumInputTracks");
79  convergence_ =
80  iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<double>("FractionOfFittedTrks");
81  inputBeamWidth_ =
82  iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<double>("InputBeamWidth", -1.);
83 
84  for (unsigned int j = 0; j < trk_Algorithm_.size(); j++)
85  algorithm_.push_back(reco::TrackBase::algoByName(trk_Algorithm_[j]));
86  for (unsigned int j = 0; j < trk_Quality_.size(); j++)
87  quality_.push_back(reco::TrackBase::qualityByName(trk_Quality_[j]));
88 
89  if (saveNtuple_ || saveBeamFit_ || savePVVertices_) {
90  outputfilename_ =
91  iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<std::string>("OutputFileName");
92  file_ = TFile::Open(outputfilename_.c_str(), "RECREATE");
93  }
94  if (saveNtuple_) {
95  ftree_ = new TTree("mytree", "mytree");
96  ftree_->AutoSave();
97 
98  ftree_->Branch("pt", &fpt, "fpt/D");
99  ftree_->Branch("d0", &fd0, "fd0/D");
100  ftree_->Branch("d0bs", &fd0bs, "fd0bs/D");
101  ftree_->Branch("sigmad0", &fsigmad0, "fsigmad0/D");
102  ftree_->Branch("phi0", &fphi0, "fphi0/D");
103  ftree_->Branch("z0", &fz0, "fz0/D");
104  ftree_->Branch("sigmaz0", &fsigmaz0, "fsigmaz0/D");
105  ftree_->Branch("theta", &ftheta, "ftheta/D");
106  ftree_->Branch("eta", &feta, "feta/D");
107  ftree_->Branch("charge", &fcharge, "fcharge/I");
108  ftree_->Branch("normchi2", &fnormchi2, "fnormchi2/D");
109  ftree_->Branch("nTotLayerMeas", &fnTotLayerMeas, "fnTotLayerMeas/i");
110  ftree_->Branch("nStripLayerMeas", &fnStripLayerMeas, "fnStripLayerMeas/i");
111  ftree_->Branch("nPixelLayerMeas", &fnPixelLayerMeas, "fnPixelLayerMeas/i");
112  ftree_->Branch("nTIBLayerMeas", &fnTIBLayerMeas, "fnTIBLayerMeas/i");
113  ftree_->Branch("nTOBLayerMeas", &fnTOBLayerMeas, "fnTOBLayerMeas/i");
114  ftree_->Branch("nTIDLayerMeas", &fnTIDLayerMeas, "fnTIDLayerMeas/i");
115  ftree_->Branch("nTECLayerMeas", &fnTECLayerMeas, "fnTECLayerMeas/i");
116  ftree_->Branch("nPXBLayerMeas", &fnPXBLayerMeas, "fnPXBLayerMeas/i");
117  ftree_->Branch("nPXFLayerMeas", &fnPXFLayerMeas, "fnPXFLayerMeas/i");
118  ftree_->Branch("cov", &fcov, "fcov[7][7]/D");
119  ftree_->Branch("vx", &fvx, "fvx/D");
120  ftree_->Branch("vy", &fvy, "fvy/D");
121  ftree_->Branch("quality", &fquality, "fquality/O");
122  ftree_->Branch("algo", &falgo, "falgo/O");
123  ftree_->Branch("run", &frun, "frun/i");
124  ftree_->Branch("lumi", &flumi, "flumi/i");
125  ftree_->Branch("pvValid", &fpvValid, "fpvValid/O");
126  ftree_->Branch("pvx", &fpvx, "fpvx/D");
127  ftree_->Branch("pvy", &fpvy, "fpvy/D");
128  ftree_->Branch("pvz", &fpvz, "fpvz/D");
129  }
130  if (saveBeamFit_) {
131  ftreeFit_ = new TTree("fitResults", "fitResults");
132  ftreeFit_->AutoSave();
133  ftreeFit_->Branch("run", &frunFit, "frunFit/i");
134  ftreeFit_->Branch("beginLumi", &fbeginLumiOfFit, "fbeginLumiOfFit/i");
135  ftreeFit_->Branch("endLumi", &fendLumiOfFit, "fendLumiOfFit/i");
136  ftreeFit_->Branch("beginTime", fbeginTimeOfFit, "fbeginTimeOfFit/C");
137  ftreeFit_->Branch("endTime", fendTimeOfFit, "fendTimeOfFit/C");
138  ftreeFit_->Branch("x", &fx, "fx/D");
139  ftreeFit_->Branch("y", &fy, "fy/D");
140  ftreeFit_->Branch("z", &fz, "fz/D");
141  ftreeFit_->Branch("sigmaZ", &fsigmaZ, "fsigmaZ/D");
142  ftreeFit_->Branch("dxdz", &fdxdz, "fdxdz/D");
143  ftreeFit_->Branch("dydz", &fdydz, "fdydz/D");
144  ftreeFit_->Branch("xErr", &fxErr, "fxErr/D");
145  ftreeFit_->Branch("yErr", &fyErr, "fyErr/D");
146  ftreeFit_->Branch("zErr", &fzErr, "fzErr/D");
147  ftreeFit_->Branch("sigmaZErr", &fsigmaZErr, "fsigmaZErr/D");
148  ftreeFit_->Branch("dxdzErr", &fdxdzErr, "fdxdzErr/D");
149  ftreeFit_->Branch("dydzErr", &fdydzErr, "fdydzErr/D");
150  ftreeFit_->Branch("widthX", &fwidthX, "fwidthX/D");
151  ftreeFit_->Branch("widthY", &fwidthY, "fwidthY/D");
152  ftreeFit_->Branch("widthXErr", &fwidthXErr, "fwidthXErr/D");
153  ftreeFit_->Branch("widthYErr", &fwidthYErr, "fwidthYErr/D");
154  }
155 
156  fBSvector.clear();
157  ftotal_tracks = 0;
158  fnTotLayerMeas = fnPixelLayerMeas = fnStripLayerMeas = fnTIBLayerMeas = 0;
159  fnTIDLayerMeas = fnTOBLayerMeas = fnTECLayerMeas = fnPXBLayerMeas = fnPXFLayerMeas = 0;
160  frun = flumi = -1;
161  frunFit = fbeginLumiOfFit = fendLumiOfFit = -1;
162  fquality = falgo = true;
163  fpvValid = true;
164  fpvx = fpvy = fpvz = 0;
165  fitted_ = false;
166  resetRefTime();
167 
168  //debug histograms
169  h1ntrks = new TH1F("h1ntrks", "number of tracks per event", 50, 0, 50);
170  h1vz_event = new TH1F("h1vz_event", "track Vz", 50, -30, 30);
171  h1cutFlow = new TH1F("h1cutFlow", "Cut flow table of track selection", 9, 0, 9);
172  h1cutFlow->GetXaxis()->SetBinLabel(1, "No cut");
173  h1cutFlow->GetXaxis()->SetBinLabel(2, "Traker hits");
174  h1cutFlow->GetXaxis()->SetBinLabel(3, "Pixel hits");
175  h1cutFlow->GetXaxis()->SetBinLabel(4, "norm. #chi^{2}");
176  h1cutFlow->GetXaxis()->SetBinLabel(5, "algo");
177  h1cutFlow->GetXaxis()->SetBinLabel(6, "quality");
178  h1cutFlow->GetXaxis()->SetBinLabel(7, "d_{0}");
179  h1cutFlow->GetXaxis()->SetBinLabel(8, "z_{0}");
180  h1cutFlow->GetXaxis()->SetBinLabel(9, "p_{T}");
181  resetCutFlow();
182 
183  // Primary vertex fitter
184  MyPVFitter = new PVFitter(iConfig, iColl);
185  MyPVFitter->resetAll();
186  if (savePVVertices_) {
187  fPVTree_ = new TTree("PrimaryVertices", "PrimaryVertices");
188  MyPVFitter->setTree(fPVTree_);
189  }
190 
191  // check filename
192  ffilename_changed = false;
193  if (writeDIPTxt_)
194  fasciiDIP.open(outputDIPTxt_.c_str());
195 }

References reco::TrackBase::algoByName(), algorithm_, appendRunTxt_, beamSpotToken_, convergence_, debug_, falgo, fasciiDIP, fbeginLumiOfFit, fbeginTimeOfFit, fBSvector, fcharge, fcov, fd0, fd0bs, fdxdz, fdxdzErr, fdydz, fdydzErr, fendLumiOfFit, fendTimeOfFit, feta, ffilename_changed, file_, fitted_, flumi, fnormchi2, fnPixelLayerMeas, fnPXBLayerMeas, fnPXFLayerMeas, fnStripLayerMeas, fnTECLayerMeas, fnTIBLayerMeas, fnTIDLayerMeas, fnTOBLayerMeas, fnTotLayerMeas, fphi0, fpt, fPVTree_, fpvValid, fpvx, fpvy, fpvz, fquality, frun, frunFit, fsigmad0, fsigmaZ, fsigmaz0, fsigmaZErr, ftheta, ftotal_tracks, ftree_, ftreeFit_, fvx, fvy, fwidthX, fwidthXErr, fwidthY, fwidthYErr, fx, fxErr, fy, fyErr, fz, fz0, fzErr, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), h1cutFlow, h1ntrks, h1vz_event, inputBeamWidth_, HLT_FULL_cff::InputTag, isMuon_, dqmiolumiharvest::j, min_Ntrks_, MyPVFitter, outputDIPTxt_, outputfilename_, outputTxt_, AlcaBeamMonitor_cfi::PVFitter, quality_, reco::TrackBase::qualityByName(), PVFitter::resetAll(), resetCutFlow(), resetRefTime(), saveBeamFit_, saveNtuple_, savePVVertices_, PVFitter::setTree(), tracksToken_, trk_Algorithm_, trk_MaxEta_, trk_MaxIP_, trk_MaxNormChi2_, trk_MaxZ_, trk_MinNPixLayers_, trk_MinNTotLayers_, trk_MinpT_, trk_Quality_, vertexToken_, writeDIPBadFit_, writeDIPTxt_, and writeTxt_.

◆ ~BeamFitter()

BeamFitter::~BeamFitter ( )
virtual

Definition at line 197 of file BeamFitter.cc.

197  {
198  if (saveNtuple_) {
199  file_->cd();
200  if (fitted_ && h1z)
201  h1z->Write();
202  h1ntrks->Write();
203  h1vz_event->Write();
204  if (h1cutFlow)
205  h1cutFlow->Write();
206  ftree_->Write();
207  }
208  if (saveBeamFit_) {
209  file_->cd();
210  ftreeFit_->Write();
211  }
212  if (savePVVertices_) {
213  file_->cd();
214  fPVTree_->Write();
215  }
216 
217  if (saveNtuple_ || saveBeamFit_ || savePVVertices_) {
218  file_->Close();
219  delete file_;
220  }
221  delete MyPVFitter;
222 }

References file_, fitted_, fPVTree_, ftree_, ftreeFit_, h1cutFlow, h1ntrks, h1vz_event, h1z, MyPVFitter, saveBeamFit_, saveNtuple_, and savePVVertices_.

Member Function Documentation

◆ dumpBWTxtFile()

void BeamFitter::dumpBWTxtFile ( std::string &  fileName)

Definition at line 656 of file BeamFitter.cc.

656  {
657  std::ofstream outFile;
658  outFile.open(fileName.c_str(), std::ios::app);
659  outFile << "---------------------------------------------------------------------------------------------------------"
660  "----------------------------------------------------"
661  << std::endl;
662  outFile << "Beam width(in cm) from Log-likelihood fit (Here we assume a symmetric beam(SigmaX=SigmaY)!)" << std::endl;
663  outFile << " " << std::endl;
664  outFile << "BeamWidth = " << fbeamWidthFit.BeamWidthX() << " +/- " << fbeamWidthFit.BeamWidthXError() << std::endl;
665  outFile.close();
666 }

References reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthXError(), fbeamWidthFit, MillePedeFileConverter_cfg::fileName, and L1TdeCSCTF_cfi::outFile.

Referenced by runBeamWidthFitter().

◆ dumpTxtFile()

void BeamFitter::dumpTxtFile ( std::string &  fileName,
bool  append 
)

Definition at line 668 of file BeamFitter.cc.

668  {
669  std::ofstream outFile;
670 
671  std::string tmpname = outputTxt_;
672  char index[15];
673  if (appendRunTxt_ && writeTxt_ && !ffilename_changed) {
674  sprintf(index, "%s%i", "_Run", frun);
675  tmpname.insert(outputTxt_.length() - 4, index);
676  fileName = tmpname;
677  ffilename_changed = true;
678  }
679 
680  if (!append)
681  outFile.open(fileName.c_str());
682  else
683  outFile.open(fileName.c_str(), std::ios::app);
684 
685  if (MyPVFitter->IsFitPerBunchCrossing()) {
686  for (std::map<int, reco::BeamSpot>::const_iterator abspot = fbspotPVMap.begin(); abspot != fbspotPVMap.end();
687  ++abspot) {
688  reco::BeamSpot beamspottmp = abspot->second;
689  int bx = abspot->first;
690 
691  outFile << "Runnumber " << frun << " bx " << bx << std::endl;
692  outFile << "BeginTimeOfFit " << fbeginTimeOfFit << " " << freftime[0] << std::endl;
693  outFile << "EndTimeOfFit " << fendTimeOfFit << " " << freftime[1] << std::endl;
694  outFile << "LumiRange " << fbeginLumiOfFit << " - " << fendLumiOfFit << std::endl;
695  outFile << "Type " << beamspottmp.type() << std::endl;
696  outFile << "X0 " << beamspottmp.x0() << std::endl;
697  outFile << "Y0 " << beamspottmp.y0() << std::endl;
698  outFile << "Z0 " << beamspottmp.z0() << std::endl;
699  outFile << "sigmaZ0 " << beamspottmp.sigmaZ() << std::endl;
700  outFile << "dxdz " << beamspottmp.dxdz() << std::endl;
701  outFile << "dydz " << beamspottmp.dydz() << std::endl;
702  outFile << "BeamWidthX " << beamspottmp.BeamWidthX() << std::endl;
703  outFile << "BeamWidthY " << beamspottmp.BeamWidthY() << std::endl;
704  for (int i = 0; i < 6; ++i) {
705  outFile << "Cov(" << i << ",j) ";
706  for (int j = 0; j < 7; ++j) {
707  outFile << beamspottmp.covariance(i, j) << " ";
708  }
709  outFile << std::endl;
710  }
711  outFile << "Cov(6,j) 0 0 0 0 0 0 " << beamspottmp.covariance(6, 6) << std::endl;
712  //}
713  outFile << "EmittanceX " << beamspottmp.emittanceX() << std::endl;
714  outFile << "EmittanceY " << beamspottmp.emittanceY() << std::endl;
715  outFile << "BetaStar " << beamspottmp.betaStar() << std::endl;
716  }
717  } //if bx results needed
718  else {
719  beamspot::BeamSpotContainer currentBS;
720 
721  currentBS.beamspot = fbeamspot;
722  currentBS.run = frun;
723  std::copy(fbeginTimeOfFit, fbeginTimeOfFit + 32, currentBS.beginTimeOfFit);
724  std::copy(fendTimeOfFit, fendTimeOfFit + 32, currentBS.endTimeOfFit);
725  currentBS.beginLumiOfFit = fbeginLumiOfFit;
726  currentBS.endLumiOfFit = fendLumiOfFit;
727  std::copy(freftime, freftime + 2, currentBS.reftime);
728 
729  beamspot::dumpBeamSpotTxt(outFile, currentBS);
730 
731  //write here Pv info for DIP only: This added only if append is false, which happen for DIP only :)
732  if (!append) {
733  outFile << "events " << (int)ForDIPPV_[0] << std::endl;
734  outFile << "meanPV " << ForDIPPV_[1] << std::endl;
735  outFile << "meanErrPV " << ForDIPPV_[2] << std::endl;
736  outFile << "rmsPV " << ForDIPPV_[3] << std::endl;
737  outFile << "rmsErrPV " << ForDIPPV_[4] << std::endl;
738  outFile << "maxPV " << (int)ForDIPPV_[5] << std::endl;
739  outFile << "nPV " << (int)ForDIPPV_[6] << std::endl;
740  } //writeDIPPVInfo_
741  } //else end here
742 
743  outFile.close();
744 }

References mps_setup::append, appendRunTxt_, beamspot::BeamSpotContainer::beamspot, reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), beamspot::BeamSpotContainer::beginLumiOfFit, beamspot::BeamSpotContainer::beginTimeOfFit, reco::BeamSpot::betaStar(), l1GtPatternGenerator_cfi::bx, filterCSVwithJSON::copy, reco::BeamSpot::covariance(), beamspot::dumpBeamSpotTxt(), reco::BeamSpot::dxdz(), reco::BeamSpot::dydz(), reco::BeamSpot::emittanceX(), reco::BeamSpot::emittanceY(), beamspot::BeamSpotContainer::endLumiOfFit, beamspot::BeamSpotContainer::endTimeOfFit, fbeamspot, fbeginLumiOfFit, fbeginTimeOfFit, fbspotPVMap, fendLumiOfFit, fendTimeOfFit, ffilename_changed, MillePedeFileConverter_cfg::fileName, ForDIPPV_, freftime, frun, mps_fire::i, createfilelist::int, PVFitter::IsFitPerBunchCrossing(), dqmiolumiharvest::j, MyPVFitter, L1TdeCSCTF_cfi::outFile, outputTxt_, beamspot::BeamSpotContainer::reftime, beamspot::BeamSpotContainer::run, reco::BeamSpot::sigmaZ(), AlCaHLTBitMon_QueryRunRegistry::string, reco::BeamSpot::type(), writeTxt_, reco::BeamSpot::x0(), reco::BeamSpot::y0(), and reco::BeamSpot::z0().

Referenced by runFitter(), and runPVandTrkFitter().

◆ fillDescription()

void BeamFitter::fillDescription ( edm::ParameterSetDescription &  iDesc)
static

Definition at line 224 of file BeamFitter.cc.

224  {
225  edm::ParameterSetDescription beamFitter;
226 
227  beamFitter.addUntracked<bool>("Debug");
228  beamFitter.addUntracked<edm::InputTag>("TrackCollection");
229  iDesc.addUntracked<edm::InputTag>("primaryVertex", edm::InputTag("offlinePrimaryVertices"));
230  iDesc.addUntracked<edm::InputTag>("beamSpot", edm::InputTag("offlineBeamSpot"));
231  beamFitter.addUntracked<bool>("WriteAscii");
232  beamFitter.addUntracked<std::string>("AsciiFileName");
233  beamFitter.addUntracked<bool>("AppendRunToFileName");
234  beamFitter.addUntracked<bool>("WriteDIPAscii");
235  // Specify whether we want to write the DIP file even if the fit is failed.
236  beamFitter.addUntracked<bool>("WriteDIPOnBadFit", true);
237  beamFitter.addUntracked<std::string>("DIPFileName");
238  beamFitter.addUntracked<bool>("SaveNtuple");
239  beamFitter.addUntracked<bool>("SaveFitResults");
240  beamFitter.addUntracked<bool>("SavePVVertices");
241  beamFitter.addUntracked<bool>("IsMuonCollection");
242 
243  beamFitter.addUntracked<double>("MinimumPt");
244  beamFitter.addUntracked<double>("MaximumEta");
245  beamFitter.addUntracked<double>("MaximumImpactParameter");
246  beamFitter.addUntracked<double>("MaximumZ");
247  beamFitter.addUntracked<int>("MinimumTotalLayers");
248  beamFitter.addUntracked<int>("MinimumPixelLayers");
249  beamFitter.addUntracked<double>("MaximumNormChi2");
250  beamFitter.addUntracked<std::vector<std::string> >("TrackAlgorithm");
251  beamFitter.addUntracked<std::vector<std::string> >("TrackQuality");
252  beamFitter.addUntracked<int>("MinimumInputTracks");
253  beamFitter.addUntracked<double>("FractionOfFittedTrks");
254  beamFitter.addUntracked<double>("InputBeamWidth", -1.);
255 
256  beamFitter.addUntracked<std::string>("OutputFileName", "");
257 
258  iDesc.add<edm::ParameterSetDescription>("BeamFitter", beamFitter);
259 }

References edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addUntracked(), HLT_FULL_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by AlcaBeamMonitor::fillDescriptions().

◆ getBeamSpot()

reco::BeamSpot BeamFitter::getBeamSpot ( )
inline

Definition at line 93 of file BeamFitter.h.

93 { return fbeamspot; }

References fbeamspot.

Referenced by BeamSpotAnalyzer::endJob(), BeamSpotAnalyzer::endLuminosityBlock(), and AlcaBeamSpotProducer::endLuminosityBlockProduce().

◆ getBeamSpotMap()

std::map<int, reco::BeamSpot> BeamFitter::getBeamSpotMap ( )
inline

Definition at line 94 of file BeamFitter.h.

94 { return fbspotPVMap; }

References fbspotPVMap.

Referenced by BeamMonitorBx::FitAndFill().

◆ getBeamWidth()

reco::BeamSpot BeamFitter::getBeamWidth ( )
inline

Definition at line 56 of file BeamFitter.h.

56 { return fbeamWidthFit; }

References fbeamWidthFit.

Referenced by BeamSpotAnalyzer::endJob().

◆ getBSvector()

std::vector<BSTrkParameters> BeamFitter::getBSvector ( )
inline

Definition at line 95 of file BeamFitter.h.

95 { return fBSvector; }

References fBSvector.

Referenced by BeamSpotAnalyzer::endLuminosityBlock(), and AlcaBeamSpotProducer::endLuminosityBlockProduce().

◆ getCutFlow()

TH1F* BeamFitter::getCutFlow ( )
inline

Definition at line 96 of file BeamFitter.h.

96 { return h1cutFlow; }

References h1cutFlow.

◆ getFitLSRange()

std::pair<int, int> BeamFitter::getFitLSRange ( )
inline

Definition at line 114 of file BeamFitter.h.

114 { return std::make_pair(fbeginLumiOfFit, fendLumiOfFit); }

References fbeginLumiOfFit, and fendLumiOfFit.

Referenced by BeamSpotAnalyzer::endJob(), BeamSpotAnalyzer::endLuminosityBlock(), AlcaBeamSpotProducer::endLuminosityBlockProduce(), and BeamMonitorBx::FitAndFill().

◆ getNPVs()

int BeamFitter::getNPVs ( )
inline

Definition at line 122 of file BeamFitter.h.

122 { return MyPVFitter->getNPVs(); }

References PVFitter::getNPVs(), and MyPVFitter.

Referenced by BeamSpotAnalyzer::endLuminosityBlock().

◆ getNPVsperBX()

const std::map<int, int>& BeamFitter::getNPVsperBX ( )
inline

Definition at line 123 of file BeamFitter.h.

123 { return MyPVFitter->getNPVsperBX(); }

References PVFitter::getNPVsperBX(), and MyPVFitter.

Referenced by BeamMonitorBx::FitAndFill().

◆ getNTracks()

int BeamFitter::getNTracks ( )
inline

Definition at line 121 of file BeamFitter.h.

121 { return fBSvector.size(); }

References fBSvector.

Referenced by BeamSpotAnalyzer::endLuminosityBlock().

◆ getPVvectorSize()

std::size_t BeamFitter::getPVvectorSize ( )
inline

Definition at line 75 of file BeamFitter.h.

75 { return (MyPVFitter->getpvStore()).size(); }

References PVFitter::getpvStore(), MyPVFitter, and findQualityFiles::size.

◆ getRefTime()

std::pair<time_t, time_t> BeamFitter::getRefTime ( )
inline

Definition at line 70 of file BeamFitter.h.

70 { return std::make_pair(freftime[0], freftime[1]); }

References freftime.

◆ getRunNumber()

int BeamFitter::getRunNumber ( )
inline

Definition at line 112 of file BeamFitter.h.

112 { return frun; }

References frun.

◆ readEvent()

void BeamFitter::readEvent ( const edm::Event &  iEvent)

Definition at line 261 of file BeamFitter.cc.

261  {
262  frun = iEvent.id().run();
263  const edm::TimeValue_t ftimestamp = iEvent.time().value();
264  const std::time_t ftmptime = ftimestamp >> 32;
265 
266  if (fbeginLumiOfFit == -1)
267  freftime[0] = freftime[1] = ftmptime;
268  if (freftime[0] == 0 || ftmptime < freftime[0])
269  freftime[0] = ftmptime;
270  if (freftime[1] == 0 || ftmptime > freftime[1])
271  freftime[1] = ftmptime;
272  // Update the human-readable string versions of the time
273  updateBTime();
274 
275  flumi = iEvent.luminosityBlock();
276  frunFit = frun;
277  if (fbeginLumiOfFit == -1 || fbeginLumiOfFit > flumi)
278  fbeginLumiOfFit = flumi;
279  if (fendLumiOfFit == -1 || fendLumiOfFit < flumi)
280  fendLumiOfFit = flumi;
281 
282  edm::Handle<reco::TrackCollection> TrackCollection;
283  iEvent.getByToken(tracksToken_, TrackCollection);
284 
285  //------ Primary Vertices
286  edm::Handle<edm::View<reco::Vertex> > PVCollection;
287  bool hasPVs = false;
288  edm::View<reco::Vertex> pv;
289  if (iEvent.getByToken(vertexToken_, PVCollection)) {
290  pv = *PVCollection;
291  hasPVs = true;
292  }
293  //------
294 
295  //------ Beam Spot in current event
296  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
297  const reco::BeamSpot *refBS = nullptr;
298  if (iEvent.getByToken(beamSpotToken_, recoBeamSpotHandle))
299  refBS = recoBeamSpotHandle.product();
300  //-------
301 
302  const reco::TrackCollection *tracks = TrackCollection.product();
303 
304  double eventZ = 0;
305  double averageZ = 0;
306 
307  for (reco::TrackCollection::const_iterator track = tracks->begin(); track != tracks->end(); ++track) {
308  if (!isMuon_) {
309  const reco::HitPattern &trkHP = track->hitPattern();
310 
311  fnPixelLayerMeas = trkHP.pixelLayersWithMeasurement();
312  fnStripLayerMeas = trkHP.stripLayersWithMeasurement();
313  fnTotLayerMeas = trkHP.trackerLayersWithMeasurement();
314  fnPXBLayerMeas = trkHP.pixelBarrelLayersWithMeasurement();
315  fnPXFLayerMeas = trkHP.pixelEndcapLayersWithMeasurement();
316  fnTIBLayerMeas = trkHP.stripTIBLayersWithMeasurement();
317  fnTIDLayerMeas = trkHP.stripTIDLayersWithMeasurement();
318  fnTOBLayerMeas = trkHP.stripTOBLayersWithMeasurement();
319  fnTECLayerMeas = trkHP.stripTECLayersWithMeasurement();
320  } else {
321  fnTotLayerMeas = track->numberOfValidHits();
322  }
323 
324  fpt = track->pt();
325  feta = track->eta();
326  fphi0 = track->phi();
327  fcharge = track->charge();
328  fnormchi2 = track->normalizedChi2();
329  fd0 = track->d0();
330  if (refBS)
331  fd0bs = -1 * track->dxy(refBS->position());
332  else
333  fd0bs = 0.;
334 
335  fsigmad0 = track->d0Error();
336  fz0 = track->dz();
337  fsigmaz0 = track->dzError();
338  ftheta = track->theta();
339  fvx = track->vx();
340  fvy = track->vy();
341 
342  for (int i = 0; i < 5; ++i) {
343  for (int j = 0; j < 5; ++j) {
344  fcov[i][j] = track->covariance(i, j);
345  }
346  }
347 
348  fquality = true;
349  falgo = true;
350 
351  if (!isMuon_) {
352  if (!quality_.empty()) {
353  fquality = false;
354  for (unsigned int i = 0; i < quality_.size(); ++i) {
355  if (debug_)
356  edm::LogInfo("BeamFitter") << "quality_[" << i << "] = " << track->qualityName(quality_[i]) << std::endl;
357  if (track->quality(quality_[i])) {
358  fquality = true;
359  break;
360  }
361  }
362  }
363 
364  // Track algorithm
365 
366  if (!algorithm_.empty()) {
367  if (std::find(algorithm_.begin(), algorithm_.end(), track->algo()) == algorithm_.end())
368  falgo = false;
369  }
370  }
371 
372  // check if we have a valid PV
373  fpvValid = false;
374 
375  if (hasPVs) {
376  for (size_t ipv = 0; ipv != pv.size(); ++ipv) {
377  if (!pv[ipv].isFake())
378  fpvValid = true;
379 
380  if (ipv == 0 && !pv[0].isFake()) {
381  fpvx = pv[0].x();
382  fpvy = pv[0].y();
383  fpvz = pv[0].z();
384  } // fix this later
385  }
386  }
387 
388  if (saveNtuple_)
389  ftree_->Fill();
390  ftotal_tracks++;
391 
392  countPass[0] = ftotal_tracks;
393  // Track selection
394  if (fnTotLayerMeas >= trk_MinNTotLayers_) {
395  countPass[1] += 1;
396  if (fnPixelLayerMeas >= trk_MinNPixLayers_) {
397  countPass[2] += 1;
398  if (fnormchi2 < trk_MaxNormChi2_) {
399  countPass[3] += 1;
400  if (falgo) {
401  countPass[4] += 1;
402  if (fquality) {
403  countPass[5] += 1;
404  if (std::abs(fd0) < trk_MaxIP_) {
405  countPass[6] += 1;
406  if (std::abs(fz0) < trk_MaxZ_) {
407  countPass[7] += 1;
408  if (fpt > trk_MinpT_) {
409  countPass[8] += 1;
410  if (std::abs(feta) < trk_MaxEta_
411  //&& fpvValid
412  ) {
413  if (debug_) {
414  edm::LogInfo("BeamFitter") << "Selected track quality = " << track->qualityMask()
415  << "; track algorithm = " << track->algoName()
416  << "= TrackAlgorithm: " << track->algo() << std::endl;
417  }
418  BSTrkParameters BSTrk(fz0, fsigmaz0, fd0, fsigmad0, fphi0, fpt, 0., 0.);
419  BSTrk.setVx(fvx);
420  BSTrk.setVy(fvy);
421  fBSvector.push_back(BSTrk);
422  averageZ += fz0;
423  }
424  }
425  }
426  }
427  }
428  }
429  }
430  }
431  } // track selection
432 
433  } // tracks
434 
435  averageZ = averageZ / (float)(fBSvector.size());
436 
437  for (std::vector<BSTrkParameters>::const_iterator iparam = fBSvector.begin(); iparam != fBSvector.end(); ++iparam) {
438  eventZ += fabs(iparam->z0() - averageZ);
439  }
440 
441  h1ntrks->Fill(fBSvector.size());
442  h1vz_event->Fill(eventZ / (float)(fBSvector.size()));
443  for (unsigned int i = 0; i < sizeof(countPass) / sizeof(countPass[0]); i++)
444  h1cutFlow->SetBinContent(i + 1, countPass[i]);
445 
446  MyPVFitter->readEvent(iEvent);
447 }

References funct::abs(), algorithm_, beamSpotToken_, countPass, debug_, falgo, fbeginLumiOfFit, fBSvector, fcharge, fcov, fd0, fd0bs, fendLumiOfFit, feta, spr::find(), dqmMemoryStats::float, flumi, fnormchi2, fnPixelLayerMeas, fnPXBLayerMeas, fnPXFLayerMeas, fnStripLayerMeas, fnTECLayerMeas, fnTIBLayerMeas, fnTIDLayerMeas, fnTOBLayerMeas, fnTotLayerMeas, fphi0, fpt, fpvValid, fpvx, fpvy, fpvz, fquality, freftime, frun, frunFit, fsigmad0, fsigmaz0, ftheta, ftotal_tracks, ftree_, fvx, fvy, fz0, h1cutFlow, h1ntrks, h1vz_event, mps_fire::i, iEvent, isMuon_, dqmiolumiharvest::j, MyPVFitter, reco::HitPattern::pixelBarrelLayersWithMeasurement(), reco::HitPattern::pixelEndcapLayersWithMeasurement(), reco::HitPattern::pixelLayersWithMeasurement(), reco::BeamSpot::position(), edm::Handle< T >::product(), MetAnalyzer::pv(), L1TEGammaOffline_cfi::PVCollection, quality_, PVFitter::readEvent(), saveNtuple_, BSTrkParameters::setVx(), BSTrkParameters::setVy(), reco::HitPattern::stripLayersWithMeasurement(), reco::HitPattern::stripTECLayersWithMeasurement(), reco::HitPattern::stripTIBLayersWithMeasurement(), reco::HitPattern::stripTIDLayersWithMeasurement(), reco::HitPattern::stripTOBLayersWithMeasurement(), HLT_FULL_cff::track, reco::HitPattern::trackerLayersWithMeasurement(), PDWG_EXOHSCP_cff::tracks, tracksToken_, trk_MaxEta_, trk_MaxIP_, trk_MaxNormChi2_, trk_MaxZ_, trk_MinNPixLayers_, trk_MinNTotLayers_, trk_MinpT_, updateBTime(), and vertexToken_.

Referenced by BeamSpotAnalyzer::analyze(), BeamMonitorBx::analyze(), and AlcaBeamSpotProducer::produce().

◆ resetCutFlow()

void BeamFitter::resetCutFlow ( )
inline

Definition at line 104 of file BeamFitter.h.

104  {
105  h1cutFlow->Reset();
106  ftotal_tracks = 0;
107  for (unsigned int i = 0; i < sizeof(countPass) / sizeof(countPass[0]); i++)
108  countPass[i] = 0;
109  }

References countPass, ftotal_tracks, h1cutFlow, and mps_fire::i.

Referenced by AlcaBeamSpotProducer::AlcaBeamSpotProducer(), BeamFitter(), BeamSpotAnalyzer::BeamSpotAnalyzer(), BeamSpotAnalyzer::endLuminosityBlock(), AlcaBeamSpotProducer::endLuminosityBlockProduce(), and BeamMonitorBx::FitAndFill().

◆ resetLSRange()

void BeamFitter::resetLSRange ( )
inline

Definition at line 60 of file BeamFitter.h.

60 { fbeginLumiOfFit = fendLumiOfFit = -1; }

References fbeginLumiOfFit, and fendLumiOfFit.

Referenced by AlcaBeamSpotProducer::AlcaBeamSpotProducer(), BeamMonitorBx::BeamMonitorBx(), BeamSpotAnalyzer::BeamSpotAnalyzer(), BeamSpotAnalyzer::endLuminosityBlock(), AlcaBeamSpotProducer::endLuminosityBlockProduce(), and BeamMonitorBx::FitAndFill().

◆ resetPVFitter()

void BeamFitter::resetPVFitter ( )
inline

Definition at line 72 of file BeamFitter.h.

72 { MyPVFitter->resetAll(); }

References MyPVFitter, and PVFitter::resetAll().

Referenced by AlcaBeamSpotProducer::AlcaBeamSpotProducer(), BeamMonitorBx::BeamMonitorBx(), BeamSpotAnalyzer::BeamSpotAnalyzer(), BeamSpotAnalyzer::endLuminosityBlock(), AlcaBeamSpotProducer::endLuminosityBlockProduce(), and BeamMonitorBx::FitAndFill().

◆ resetRefTime()

void BeamFitter::resetRefTime ( )
inline

Definition at line 61 of file BeamFitter.h.

61 { freftime[0] = freftime[1] = 0; }

References freftime.

Referenced by AlcaBeamSpotProducer::AlcaBeamSpotProducer(), BeamFitter(), BeamMonitorBx::BeamMonitorBx(), BeamSpotAnalyzer::BeamSpotAnalyzer(), BeamSpotAnalyzer::endLuminosityBlock(), AlcaBeamSpotProducer::endLuminosityBlockProduce(), and BeamMonitorBx::FitAndFill().

◆ resetTotTrk()

void BeamFitter::resetTotTrk ( )
inline

Definition at line 59 of file BeamFitter.h.

59 { ftotal_tracks = 0; }

References ftotal_tracks.

◆ resetTrkVector()

void BeamFitter::resetTrkVector ( )
inline

Definition at line 58 of file BeamFitter.h.

58 { fBSvector.clear(); }

References fBSvector.

Referenced by AlcaBeamSpotProducer::AlcaBeamSpotProducer(), BeamMonitorBx::BeamMonitorBx(), BeamSpotAnalyzer::BeamSpotAnalyzer(), BeamSpotAnalyzer::endLuminosityBlock(), AlcaBeamSpotProducer::endLuminosityBlockProduce(), and BeamMonitorBx::FitAndFill().

◆ resizeBSvector()

void BeamFitter::resizeBSvector ( unsigned int  nsize)
inline

Definition at line 77 of file BeamFitter.h.

77 { fBSvector.erase(fBSvector.begin(), fBSvector.begin() + nsize); }

References fBSvector.

◆ resizePVvector()

void BeamFitter::resizePVvector ( unsigned int  npvsize)
inline

Definition at line 80 of file BeamFitter.h.

80 { MyPVFitter->resizepvStore(npvsize); }

References MyPVFitter, and PVFitter::resizepvStore().

◆ runAllFitter()

void BeamFitter::runAllFitter ( )

Definition at line 784 of file BeamFitter.cc.

784  {
785  if (!fBSvector.empty()) {
786  BSFitter *myalgo = new BSFitter(fBSvector);
787  fbeamspot = myalgo->Fit_d0phi();
788 
789  // iterative
790  if (debug_)
791  std::cout << " d0-phi Iterative:" << std::endl;
792  BSFitter *myitealgo = new BSFitter(fBSvector);
793  myitealgo->Setd0Cut_d0phi(4.0);
794  reco::BeamSpot beam_ite = myitealgo->Fit_ited0phi();
795  if (debug_) {
796  std::cout << beam_ite << std::endl;
797  std::cout << "\n Now run tests of the different fits\n";
798  }
799  // from here are just tests
800  std::string fit_type = "chi2";
801  myalgo->SetFitVariable(std::string("z"));
802  myalgo->SetFitType(std::string("chi2"));
803  reco::BeamSpot beam_fit_z_chi2 = myalgo->Fit();
804  if (debug_) {
805  std::cout << " z Chi2 Fit ONLY:" << std::endl;
806  std::cout << beam_fit_z_chi2 << std::endl;
807  }
808 
809  fit_type = "combined";
810  myalgo->SetFitVariable("z");
811  myalgo->SetFitType("combined");
812  reco::BeamSpot beam_fit_z_lh = myalgo->Fit();
813  if (debug_) {
814  std::cout << " z Log-Likelihood Fit ONLY:" << std::endl;
815  std::cout << beam_fit_z_lh << std::endl;
816  }
817 
818  myalgo->SetFitVariable("d");
819  myalgo->SetFitType("d0phi");
820  reco::BeamSpot beam_fit_dphi = myalgo->Fit();
821  if (debug_) {
822  std::cout << " d0-phi0 Fit ONLY:" << std::endl;
823  std::cout << beam_fit_dphi << std::endl;
824  }
825  /*
826  myalgo->SetFitVariable(std::string("d*z"));
827  myalgo->SetFitType(std::string("likelihood"));
828  reco::BeamSpot beam_fit_dz_lh = myalgo->Fit();
829  if (debug_){
830  std::cout << " Log-Likelihood Fit:" << std::endl;
831  std::cout << beam_fit_dz_lh << std::endl;
832  }
833 
834  myalgo->SetFitVariable(std::string("d*z"));
835  myalgo->SetFitType(std::string("resolution"));
836  reco::BeamSpot beam_fit_dresz_lh = myalgo->Fit();
837  if(debug_){
838  std::cout << " IP Resolution Fit" << std::endl;
839  std::cout << beam_fit_dresz_lh << std::endl;
840 
841  std::cout << "c0 = " << myalgo->GetResPar0() << " +- " << myalgo->GetResPar0Err() << std::endl;
842  std::cout << "c1 = " << myalgo->GetResPar1() << " +- " << myalgo->GetResPar1Err() << std::endl;
843  }
844  */
845  } else if (debug_)
846  std::cout << "No good track selected! No beam fit!" << std::endl;
847 }

References gather_cfg::cout, debug_, fbeamspot, fBSvector, BSFitter::Fit(), BSFitter::Fit_d0phi(), BSFitter::Fit_ited0phi(), BSFitter::Setd0Cut_d0phi(), BSFitter::SetFitType(), BSFitter::SetFitVariable(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by BeamSpotAnalyzer::endJob().

◆ runBeamWidthFitter()

bool BeamFitter::runBeamWidthFitter ( )

Definition at line 620 of file BeamFitter.cc.

620  {
621  bool widthfit_ok = false;
622  // default fit to extract beam spot info
623  if (fBSvector.size() > 1) {
624  edm::LogInfo("BeamFitter") << "Calculating beam spot positions('d0-phi0' method) and width using llh Fit"
625  << std::endl
626  << "We will use " << fBSvector.size() << " good tracks out of " << ftotal_tracks
627  << std::endl;
628 
629  BSFitter *myalgo = new BSFitter(fBSvector);
630  myalgo->SetMaximumZ(trk_MaxZ_);
631  myalgo->SetConvergence(convergence_);
632  myalgo->SetMinimumNTrks(min_Ntrks_);
633  if (inputBeamWidth_ > 0)
634  myalgo->SetInputBeamWidth(inputBeamWidth_);
635 
636  myalgo->SetFitVariable(std::string("d*z"));
637  myalgo->SetFitType(std::string("likelihood"));
638  fbeamWidthFit = myalgo->Fit();
639 
640  //Add to .txt file
641  if (writeTxt_)
642  dumpBWTxtFile(outputTxt_);
643 
644  delete myalgo;
645 
646  // not fake
647  if (fbeamspot.type() != 0)
648  widthfit_ok = true;
649  } else {
650  fbeamspot.setType(reco::BeamSpot::Fake);
651  edm::LogWarning("BeamFitter") << "Not enough good tracks selected! No beam fit!" << std::endl;
652  }
653  return widthfit_ok;
654 }

References convergence_, dumpBWTxtFile(), reco::BeamSpot::Fake, fbeamspot, fbeamWidthFit, fBSvector, BSFitter::Fit(), ftotal_tracks, inputBeamWidth_, min_Ntrks_, outputTxt_, BSFitter::SetConvergence(), BSFitter::SetFitType(), BSFitter::SetFitVariable(), BSFitter::SetInputBeamWidth(), BSFitter::SetMaximumZ(), BSFitter::SetMinimumNTrks(), reco::BeamSpot::setType(), AlCaHLTBitMon_QueryRunRegistry::string, trk_MaxZ_, reco::BeamSpot::type(), and writeTxt_.

Referenced by BeamSpotAnalyzer::endJob().

◆ runFitter()

bool BeamFitter::runFitter ( )

Definition at line 609 of file BeamFitter.cc.

609  {
610  bool fit_ok = runFitterNoTxt();
611 
612  if (writeTxt_)
613  dumpTxtFile(outputTxt_, true); // all reaults
614  if (writeDIPTxt_ && (fit_ok || writeDIPBadFit_)) {
615  dumpTxtFile(outputDIPTxt_, false); // for DQM/DIP
616  }
617  return fit_ok;
618 }

References dumpTxtFile(), outputDIPTxt_, outputTxt_, runFitterNoTxt(), writeDIPBadFit_, writeDIPTxt_, and writeTxt_.

◆ runFitterNoTxt()

bool BeamFitter::runFitterNoTxt ( )

Definition at line 546 of file BeamFitter.cc.

546  {
547  edm::LogInfo("BeamFitter") << " [BeamFitterDebugTime] freftime[0] = " << freftime[0]
548  << "; address = " << &freftime[0] << " = " << fbeginTimeOfFit << std::endl;
549  edm::LogInfo("BeamFitter") << " [BeamFitterDebugTime] freftime[1] = " << freftime[1]
550  << "; address = " << &freftime[1] << " = " << fendTimeOfFit << std::endl;
551 
552  if (fbeginLumiOfFit == -1 || fendLumiOfFit == -1) {
553  edm::LogWarning("BeamFitter") << "No event read! No Fitting!" << std::endl;
554  return false;
555  }
556 
557  bool fit_ok = false;
558  // default fit to extract beam spot info
559  if (fBSvector.size() > 1) {
560  edm::LogInfo("BeamFitter") << "Calculating beam spot..." << std::endl
561  << "We will use " << fBSvector.size() << " good tracks out of " << ftotal_tracks
562  << std::endl;
563 
564  BSFitter *myalgo = new BSFitter(fBSvector);
565  myalgo->SetMaximumZ(trk_MaxZ_);
566  myalgo->SetConvergence(convergence_);
567  myalgo->SetMinimumNTrks(min_Ntrks_);
568  if (inputBeamWidth_ > 0)
569  myalgo->SetInputBeamWidth(inputBeamWidth_);
570 
571  fbeamspot = myalgo->Fit();
572 
573  // retrieve histogram for Vz
574  h1z = (TH1F *)myalgo->GetVzHisto();
575 
576  delete myalgo;
577  if (fbeamspot.type() > 0) { // save all results except for Fake and Unknown (all 0.)
578  fit_ok = true;
579  if (saveBeamFit_) {
580  fx = fbeamspot.x0();
581  fy = fbeamspot.y0();
582  fz = fbeamspot.z0();
583  fsigmaZ = fbeamspot.sigmaZ();
584  fdxdz = fbeamspot.dxdz();
585  fdydz = fbeamspot.dydz();
586  fwidthX = fbeamspot.BeamWidthX();
587  fwidthY = fbeamspot.BeamWidthY();
588  fxErr = fbeamspot.x0Error();
589  fyErr = fbeamspot.y0Error();
590  fzErr = fbeamspot.z0Error();
591  fsigmaZErr = fbeamspot.sigmaZ0Error();
592  fdxdzErr = fbeamspot.dxdzError();
593  fdydzErr = fbeamspot.dydzError();
594  fwidthXErr = fbeamspot.BeamWidthXError();
595  fwidthYErr = fbeamspot.BeamWidthYError();
596  ftreeFit_->Fill();
597  }
598  }
599  } else { // tracks <= 1
600  reco::BeamSpot tmpbs;
601  fbeamspot = tmpbs;
602  fbeamspot.setType(reco::BeamSpot::Fake);
603  edm::LogInfo("BeamFitter") << "Not enough good tracks selected! No beam fit!" << std::endl;
604  }
605  fitted_ = true;
606  return fit_ok;
607 }

References reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthXError(), reco::BeamSpot::BeamWidthY(), reco::BeamSpot::BeamWidthYError(), convergence_, reco::BeamSpot::dxdz(), reco::BeamSpot::dxdzError(), reco::BeamSpot::dydz(), reco::BeamSpot::dydzError(), reco::BeamSpot::Fake, fbeamspot, fbeginLumiOfFit, fbeginTimeOfFit, fBSvector, fdxdz, fdxdzErr, fdydz, fdydzErr, fendLumiOfFit, fendTimeOfFit, BSFitter::Fit(), fitted_, freftime, fsigmaZ, fsigmaZErr, ftotal_tracks, ftreeFit_, fwidthX, fwidthXErr, fwidthY, fwidthYErr, fx, fxErr, fy, fyErr, fz, fzErr, BSFitter::GetVzHisto(), h1z, inputBeamWidth_, min_Ntrks_, saveBeamFit_, BSFitter::SetConvergence(), BSFitter::SetInputBeamWidth(), BSFitter::SetMaximumZ(), BSFitter::SetMinimumNTrks(), reco::BeamSpot::setType(), reco::BeamSpot::sigmaZ(), reco::BeamSpot::sigmaZ0Error(), trk_MaxZ_, reco::BeamSpot::type(), reco::BeamSpot::x0(), reco::BeamSpot::x0Error(), reco::BeamSpot::y0(), reco::BeamSpot::y0Error(), reco::BeamSpot::z0(), and reco::BeamSpot::z0Error().

Referenced by runFitter(), and runPVandTrkFitter().

◆ runPVandTrkFitter()

bool BeamFitter::runPVandTrkFitter ( )

Definition at line 449 of file BeamFitter.cc.

449  {
450  // run both PV and track fitters
451  bool fit_ok = false;
452  bool pv_fit_ok = false;
453  reco::BeamSpot bspotPV;
454  reco::BeamSpot bspotTrk;
455 
456  // First run PV fitter
457  if (MyPVFitter->IsFitPerBunchCrossing()) {
458  edm::LogInfo("BeamFitter") << " [BeamFitterBxDebugTime] freftime[0] = " << freftime[0]
459  << "; address = " << &freftime[0] << " = " << fbeginTimeOfFit << std::endl;
460  edm::LogInfo("BeamFitter") << " [BeamFitterBxDebugTime] freftime[1] = " << freftime[1]
461  << "; address = " << &freftime[1] << " = " << fendTimeOfFit << std::endl;
462 
463  if (MyPVFitter->runBXFitter()) {
464  fbspotPVMap = MyPVFitter->getBeamSpotMap();
465  pv_fit_ok = true;
466  }
467  if (writeTxt_)
468  dumpTxtFile(outputTxt_, true); // all reaults
469  if (writeDIPTxt_ && (pv_fit_ok || writeDIPBadFit_)) {
470  dumpTxtFile(outputDIPTxt_, false); // for DQM/DIP
471  }
472  return pv_fit_ok;
473  }
474 
475  if (MyPVFitter->runFitter()) {
476  bspotPV = MyPVFitter->getBeamSpot();
477 
478  // take beam width from PV fit and pass it to track fitter
479  // assume circular transverse beam width
480  inputBeamWidth_ = sqrt(pow(bspotPV.BeamWidthX(), 2) + pow(bspotPV.BeamWidthY(), 2)) / sqrt(2);
481  pv_fit_ok = true;
482 
483  } else {
484  // problems with PV fit
485  bspotPV.setType(reco::BeamSpot::Unknown);
486  bspotTrk.setType(reco::BeamSpot::Unknown); //propagate error to trk beam spot
487  }
488 
489  if (runFitterNoTxt()) {
490  bspotTrk = fbeamspot;
491  fit_ok = true;
492  } else {
493  // beamfit failed, flagged as empty beam spot
494  bspotTrk.setType(reco::BeamSpot::Fake);
495  fit_ok = false;
496  }
497 
498  // combined results into one single beam spot
499 
500  // to pass errors I have to create another beam spot object
501  reco::BeamSpot::CovarianceMatrix matrix;
502  for (int j = 0; j < 7; ++j) {
503  for (int k = j; k < 7; ++k) {
504  matrix(j, k) = bspotTrk.covariance(j, k);
505  }
506  }
507  // change beam width error to one from PV
508  if (pv_fit_ok && fit_ok) {
509  matrix(6, 6) = MyPVFitter->getWidthXerr() * MyPVFitter->getWidthXerr();
510 
511  // get Z and sigmaZ from PV fit
512  matrix(2, 2) = bspotPV.covariance(2, 2);
513  matrix(3, 3) = bspotPV.covariance(3, 3);
514  reco::BeamSpot tmpbs(reco::BeamSpot::Point(bspotTrk.x0(), bspotTrk.y0(), bspotPV.z0()),
515  bspotPV.sigmaZ(),
516  bspotTrk.dxdz(),
517  bspotTrk.dydz(),
518  bspotPV.BeamWidthX(),
519  matrix,
520  bspotPV.type());
521  tmpbs.setBeamWidthY(bspotPV.BeamWidthY());
522  // overwrite beam spot result
523  fbeamspot = tmpbs;
524  }
525  if (pv_fit_ok && fit_ok) {
526  fbeamspot.setType(bspotPV.type());
527  } else if (!pv_fit_ok && fit_ok) {
528  fbeamspot.setType(reco::BeamSpot::Unknown);
529  } else if (pv_fit_ok && !fit_ok) {
530  fbeamspot.setType(reco::BeamSpot::Unknown);
531  } else if (!pv_fit_ok && !fit_ok) {
532  fbeamspot.setType(reco::BeamSpot::Fake);
533  }
534 
535  if (writeTxt_)
536  dumpTxtFile(outputTxt_, true); // all reaults
537  if (writeDIPTxt_ && ((fit_ok && pv_fit_ok) || writeDIPBadFit_)) {
538  dumpTxtFile(outputDIPTxt_, false); // for DQM/DIP
539  for (size_t i = 0; i < 7; i++)
540  ForDIPPV_.push_back(0.0);
541  }
542 
543  return fit_ok && pv_fit_ok;
544 }

References reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), reco::BeamSpot::covariance(), dumpTxtFile(), reco::BeamSpot::dxdz(), reco::BeamSpot::dydz(), reco::BeamSpot::Fake, fbeamspot, fbeginTimeOfFit, fbspotPVMap, fendTimeOfFit, ForDIPPV_, freftime, PVFitter::getBeamSpot(), PVFitter::getBeamSpotMap(), PVFitter::getWidthXerr(), mps_fire::i, inputBeamWidth_, PVFitter::IsFitPerBunchCrossing(), dqmiolumiharvest::j, dqmdumpme::k, makeMuonMisalignmentScenario::matrix, MyPVFitter, outputDIPTxt_, outputTxt_, funct::pow(), PVFitter::runBXFitter(), PVFitter::runFitter(), runFitterNoTxt(), reco::BeamSpot::setBeamWidthY(), reco::BeamSpot::setType(), reco::BeamSpot::sigmaZ(), mathSSE::sqrt(), reco::BeamSpot::type(), reco::BeamSpot::Unknown, writeDIPBadFit_, writeDIPTxt_, writeTxt_, reco::BeamSpot::x0(), reco::BeamSpot::y0(), and reco::BeamSpot::z0().

Referenced by BeamSpotAnalyzer::endJob(), BeamSpotAnalyzer::endLuminosityBlock(), AlcaBeamSpotProducer::endLuminosityBlockProduce(), and BeamMonitorBx::FitAndFill().

◆ setFitLSRange()

void BeamFitter::setFitLSRange ( int  ls0,
int  ls1 
)
inline

Definition at line 115 of file BeamFitter.h.

115  {
116  fbeginLumiOfFit = ls0;
117  fendLumiOfFit = ls1;
118  }

References fbeginLumiOfFit, and fendLumiOfFit.

Referenced by BeamSpotAnalyzer::endLuminosityBlock(), and AlcaBeamSpotProducer::endLuminosityBlockProduce().

◆ SetPVInfo()

void BeamFitter::SetPVInfo ( const std::vector< float > &  v1_)
inline

Definition at line 83 of file BeamFitter.h.

83  {
84  ForDIPPV_.clear();
85  ForDIPPV_.assign(v1_.begin(), v1_.end());
86  }

References ForDIPPV_.

◆ setRefTime()

void BeamFitter::setRefTime ( time_t  t0,
time_t  t1 
)
inline

Definition at line 62 of file BeamFitter.h.

62  {
63  freftime[0] = t0;
64  freftime[1] = t1;
65  // Make sure the string representation of the time
66  // is up-to-date
67  updateBTime();
68  }

References freftime, FrontierCondition_GT_autoExpress_cfi::t0, RandomServiceHelper::t1, and updateBTime().

Referenced by BeamSpotAnalyzer::endLuminosityBlock(), and AlcaBeamSpotProducer::endLuminosityBlockProduce().

◆ setRun()

void BeamFitter::setRun ( int  run)
inline

Definition at line 119 of file BeamFitter.h.

119 { frun = run; }

References frun, and writedatasetfile::run.

Referenced by BeamSpotAnalyzer::endLuminosityBlock(), and AlcaBeamSpotProducer::endLuminosityBlockProduce().

◆ subtractFromCutFlow()

void BeamFitter::subtractFromCutFlow ( const TH1F *  toSubtract)
inline

Definition at line 97 of file BeamFitter.h.

97  {
98  h1cutFlow->Add(toSubtract, -1.0);
99  for (unsigned int i = 0; i < sizeof(countPass) / sizeof(countPass[0]); i++) {
100  countPass[i] = h1cutFlow->GetBinContent(i + 1);
101  }
102  }

References countPass, h1cutFlow, and mps_fire::i.

◆ updateBTime()

void BeamFitter::updateBTime ( )
private

class: BeamFitter.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 29 of file BeamFitter.cc.

29  {
30  char ts[] = "yyyy.mn.dd hh:mm:ss zzz ";
31  char *fbeginTime = ts;
32  strftime(fbeginTime, sizeof(ts), "%Y.%m.%d %H:%M:%S GMT", gmtime(&freftime[0]));
33  sprintf(fbeginTimeOfFit, "%s", fbeginTime);
34  char *fendTime = ts;
35  strftime(fendTime, sizeof(ts), "%Y.%m.%d %H:%M:%S GMT", gmtime(&freftime[1]));
36  sprintf(fendTimeOfFit, "%s", fendTime);
37 }

References fbeginTimeOfFit, fendTimeOfFit, and freftime.

Referenced by readEvent(), and setRefTime().

◆ write2DB()

void BeamFitter::write2DB ( )

Definition at line 746 of file BeamFitter.cc.

746  {
747  BeamSpotObjects *pBSObjects = new BeamSpotObjects();
748 
749  pBSObjects->SetPosition(fbeamspot.position().X(), fbeamspot.position().Y(), fbeamspot.position().Z());
750  //std::cout << " wrote: x= " << fbeamspot.position().X() << " y= "<< fbeamspot.position().Y() << " z= " << fbeamspot.position().Z() << std::endl;
751  pBSObjects->SetSigmaZ(fbeamspot.sigmaZ());
752  pBSObjects->Setdxdz(fbeamspot.dxdz());
753  pBSObjects->Setdydz(fbeamspot.dydz());
754  //if (inputBeamWidth_ > 0 ) {
755  // std::cout << " beam width value forced to be " << inputBeamWidth_ << std::endl;
756  // pBSObjects->SetBeamWidthX(inputBeamWidth_);
757  // pBSObjects->SetBeamWidthY(inputBeamWidth_);
758  //} else {
759  // need to fix this
760  //std::cout << " using default value, 15e-4, for beam width!!!"<<std::endl;
761  pBSObjects->SetBeamWidthX(fbeamspot.BeamWidthX());
762  pBSObjects->SetBeamWidthY(fbeamspot.BeamWidthY());
763  //}
764 
765  for (int i = 0; i < 7; ++i) {
766  for (int j = 0; j < 7; ++j) {
767  pBSObjects->SetCovariance(i, j, fbeamspot.covariance(i, j));
768  }
769  }
770  edm::Service<cond::service::PoolDBOutputService> poolDbService;
771  if (poolDbService.isAvailable()) {
772  std::cout << "poolDBService available" << std::endl;
773  if (poolDbService->isNewTagRequest("BeamSpotObjectsRcd")) {
774  std::cout << "new tag requested" << std::endl;
775  poolDbService->createNewIOV<BeamSpotObjects>(
776  pBSObjects, poolDbService->beginOfTime(), poolDbService->endOfTime(), "BeamSpotObjectsRcd");
777  } else {
778  std::cout << "no new tag requested" << std::endl;
779  poolDbService->appendSinceTime<BeamSpotObjects>(pBSObjects, poolDbService->currentTime(), "BeamSpotObjectsRcd");
780  }
781  }
782 }

References cond::service::PoolDBOutputService::appendSinceTime(), reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), cond::service::PoolDBOutputService::beginOfTime(), gather_cfg::cout, reco::BeamSpot::covariance(), cond::service::PoolDBOutputService::createNewIOV(), cond::service::PoolDBOutputService::currentTime(), reco::BeamSpot::dxdz(), reco::BeamSpot::dydz(), cond::service::PoolDBOutputService::endOfTime(), fbeamspot, mps_fire::i, edm::Service< T >::isAvailable(), cond::service::PoolDBOutputService::isNewTagRequest(), dqmiolumiharvest::j, reco::BeamSpot::position(), BeamSpotObjects::SetBeamWidthX(), BeamSpotObjects::SetBeamWidthY(), BeamSpotObjects::SetCovariance(), BeamSpotObjects::Setdxdz(), BeamSpotObjects::Setdydz(), BeamSpotObjects::SetPosition(), BeamSpotObjects::SetSigmaZ(), and reco::BeamSpot::sigmaZ().

Referenced by BeamSpotAnalyzer::endJob().

Member Data Documentation

◆ algorithm_

std::vector<reco::TrackBase::TrackAlgorithm> BeamFitter::algorithm_
private

Definition at line 156 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ appendRunTxt_

bool BeamFitter::appendRunTxt_
private

Definition at line 137 of file BeamFitter.h.

Referenced by BeamFitter(), and dumpTxtFile().

◆ beamSpotToken_

edm::EDGetTokenT<reco::BeamSpot> BeamFitter::beamSpotToken_
private

Definition at line 140 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ convergence_

double BeamFitter::convergence_
private

Definition at line 158 of file BeamFitter.h.

Referenced by BeamFitter(), runBeamWidthFitter(), and runFitterNoTxt().

◆ countPass

int BeamFitter::countPass[9]
private

Definition at line 236 of file BeamFitter.h.

Referenced by readEvent(), resetCutFlow(), and subtractFromCutFlow().

◆ debug_

bool BeamFitter::debug_
private

Definition at line 136 of file BeamFitter.h.

Referenced by BeamFitter(), readEvent(), and runAllFitter().

◆ falgo

bool BeamFitter::falgo
private

Definition at line 204 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fasciiDIP

std::ofstream BeamFitter::fasciiDIP
private

Definition at line 134 of file BeamFitter.h.

Referenced by BeamFitter().

◆ fasciiFile

std::ofstream BeamFitter::fasciiFile
private

Definition at line 133 of file BeamFitter.h.

◆ fbeamspot

reco::BeamSpot BeamFitter::fbeamspot
private

Definition at line 129 of file BeamFitter.h.

Referenced by dumpTxtFile(), getBeamSpot(), runAllFitter(), runBeamWidthFitter(), runFitterNoTxt(), runPVandTrkFitter(), and write2DB().

◆ fbeamWidthFit

reco::BeamSpot BeamFitter::fbeamWidthFit
private

Definition at line 130 of file BeamFitter.h.

Referenced by dumpBWTxtFile(), getBeamWidth(), and runBeamWidthFitter().

◆ fbeginLumiOfFit

int BeamFitter::fbeginLumiOfFit
private

Definition at line 212 of file BeamFitter.h.

Referenced by BeamFitter(), dumpTxtFile(), getFitLSRange(), readEvent(), resetLSRange(), runFitterNoTxt(), and setFitLSRange().

◆ fbeginTimeOfFit

char BeamFitter::fbeginTimeOfFit[32]
private

Definition at line 214 of file BeamFitter.h.

Referenced by BeamFitter(), dumpTxtFile(), runFitterNoTxt(), runPVandTrkFitter(), and updateBTime().

◆ fbspotPVMap

std::map<int, reco::BeamSpot> BeamFitter::fbspotPVMap
private

Definition at line 131 of file BeamFitter.h.

Referenced by dumpTxtFile(), getBeamSpotMap(), and runPVandTrkFitter().

◆ fBSvector

std::vector<BSTrkParameters> BeamFitter::fBSvector
private

Definition at line 128 of file BeamFitter.h.

Referenced by BeamFitter(), getBSvector(), getNTracks(), readEvent(), resetTrkVector(), resizeBSvector(), runAllFitter(), runBeamWidthFitter(), and runFitterNoTxt().

◆ fcharge

int BeamFitter::fcharge
private

Definition at line 179 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fcov

double BeamFitter::fcov[7][7]
private

Definition at line 198 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fd0

double BeamFitter::fd0
private

Definition at line 182 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fd0bs

double BeamFitter::fd0bs
private

Definition at line 183 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fdxdz

double BeamFitter::fdxdz
private

Definition at line 220 of file BeamFitter.h.

Referenced by BeamFitter(), and runFitterNoTxt().

◆ fdxdzErr

double BeamFitter::fdxdzErr
private

Definition at line 226 of file BeamFitter.h.

Referenced by BeamFitter(), and runFitterNoTxt().

◆ fdydz

double BeamFitter::fdydz
private

Definition at line 221 of file BeamFitter.h.

Referenced by BeamFitter(), and runFitterNoTxt().

◆ fdydzErr

double BeamFitter::fdydzErr
private

Definition at line 227 of file BeamFitter.h.

Referenced by BeamFitter(), and runFitterNoTxt().

◆ fendLumiOfFit

int BeamFitter::fendLumiOfFit
private

Definition at line 213 of file BeamFitter.h.

Referenced by BeamFitter(), dumpTxtFile(), getFitLSRange(), readEvent(), resetLSRange(), runFitterNoTxt(), and setFitLSRange().

◆ fendTimeOfFit

char BeamFitter::fendTimeOfFit[32]
private

Definition at line 215 of file BeamFitter.h.

Referenced by BeamFitter(), dumpTxtFile(), runFitterNoTxt(), runPVandTrkFitter(), and updateBTime().

◆ feta

double BeamFitter::feta
private

Definition at line 178 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ ffilename_changed

bool BeamFitter::ffilename_changed
private

Definition at line 163 of file BeamFitter.h.

Referenced by BeamFitter(), and dumpTxtFile().

◆ file_

TFile* BeamFitter::file_
private

Definition at line 174 of file BeamFitter.h.

Referenced by BeamFitter(), and ~BeamFitter().

◆ fitted_

bool BeamFitter::fitted_
private

Definition at line 162 of file BeamFitter.h.

Referenced by BeamFitter(), runFitterNoTxt(), and ~BeamFitter().

◆ flumi

int BeamFitter::flumi
private

Definition at line 202 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fnormchi2

double BeamFitter::fnormchi2
private

Definition at line 180 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fnPixelLayerMeas

int BeamFitter::fnPixelLayerMeas
private

Definition at line 188 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fnPXBLayerMeas

int BeamFitter::fnPXBLayerMeas
private

Definition at line 194 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fnPXFLayerMeas

int BeamFitter::fnPXFLayerMeas
private

Definition at line 195 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fnStripLayerMeas

int BeamFitter::fnStripLayerMeas
private

Definition at line 189 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fnTECLayerMeas

int BeamFitter::fnTECLayerMeas
private

Definition at line 193 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fnTIBLayerMeas

int BeamFitter::fnTIBLayerMeas
private

Definition at line 190 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fnTIDLayerMeas

int BeamFitter::fnTIDLayerMeas
private

Definition at line 191 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fnTOBLayerMeas

int BeamFitter::fnTOBLayerMeas
private

Definition at line 192 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fnTotLayerMeas

int BeamFitter::fnTotLayerMeas
private

Definition at line 187 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ ForDIPPV_

std::vector<float> BeamFitter::ForDIPPV_
private

Definition at line 166 of file BeamFitter.h.

Referenced by dumpTxtFile(), runPVandTrkFitter(), and SetPVInfo().

◆ fphi0

double BeamFitter::fphi0
private

Definition at line 181 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fpt

double BeamFitter::fpt
private

Definition at line 177 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fPVTree_

TTree* BeamFitter::fPVTree_
private

Definition at line 239 of file BeamFitter.h.

Referenced by BeamFitter(), and ~BeamFitter().

◆ fpvValid

bool BeamFitter::fpvValid
private

Definition at line 205 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fpvx

double BeamFitter::fpvx
private

Definition at line 206 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fpvy

double BeamFitter::fpvy
private

Definition at line 206 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fpvz

double BeamFitter::fpvz
private

Definition at line 206 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fquality

bool BeamFitter::fquality
private

Definition at line 203 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ freftime

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

Definition at line 207 of file BeamFitter.h.

Referenced by dumpTxtFile(), getRefTime(), readEvent(), resetRefTime(), runFitterNoTxt(), runPVandTrkFitter(), setRefTime(), and updateBTime().

◆ frun

int BeamFitter::frun
private

Definition at line 201 of file BeamFitter.h.

Referenced by BeamFitter(), dumpTxtFile(), getRunNumber(), readEvent(), and setRun().

◆ frunFit

int BeamFitter::frunFit
private

Definition at line 211 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fsigmad0

double BeamFitter::fsigmad0
private

Definition at line 184 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fsigmaZ

double BeamFitter::fsigmaZ
private

Definition at line 219 of file BeamFitter.h.

Referenced by BeamFitter(), and runFitterNoTxt().

◆ fsigmaz0

double BeamFitter::fsigmaz0
private

Definition at line 186 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fsigmaZErr

double BeamFitter::fsigmaZErr
private

Definition at line 225 of file BeamFitter.h.

Referenced by BeamFitter(), and runFitterNoTxt().

◆ ftheta

double BeamFitter::ftheta
private

Definition at line 176 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ ftotal_tracks

int BeamFitter::ftotal_tracks
private

Definition at line 159 of file BeamFitter.h.

Referenced by BeamFitter(), readEvent(), resetCutFlow(), resetTotTrk(), runBeamWidthFitter(), and runFitterNoTxt().

◆ ftree_

TTree* BeamFitter::ftree_
private

Definition at line 175 of file BeamFitter.h.

Referenced by BeamFitter(), readEvent(), and ~BeamFitter().

◆ ftreeFit_

TTree* BeamFitter::ftreeFit_
private

Definition at line 210 of file BeamFitter.h.

Referenced by BeamFitter(), runFitterNoTxt(), and ~BeamFitter().

◆ fvx

double BeamFitter::fvx
private

Definition at line 199 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fvy

double BeamFitter::fvy
private

Definition at line 200 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fwidthX

double BeamFitter::fwidthX
private

Definition at line 228 of file BeamFitter.h.

Referenced by BeamFitter(), and runFitterNoTxt().

◆ fwidthXErr

double BeamFitter::fwidthXErr
private

Definition at line 230 of file BeamFitter.h.

Referenced by BeamFitter(), and runFitterNoTxt().

◆ fwidthY

double BeamFitter::fwidthY
private

Definition at line 229 of file BeamFitter.h.

Referenced by BeamFitter(), and runFitterNoTxt().

◆ fwidthYErr

double BeamFitter::fwidthYErr
private

Definition at line 231 of file BeamFitter.h.

Referenced by BeamFitter(), and runFitterNoTxt().

◆ fx

double BeamFitter::fx
private

Definition at line 216 of file BeamFitter.h.

Referenced by BeamFitter(), and runFitterNoTxt().

◆ fxErr

double BeamFitter::fxErr
private

Definition at line 222 of file BeamFitter.h.

Referenced by BeamFitter(), and runFitterNoTxt().

◆ fy

double BeamFitter::fy
private

Definition at line 217 of file BeamFitter.h.

Referenced by BeamFitter(), and runFitterNoTxt().

◆ fyErr

double BeamFitter::fyErr
private

Definition at line 223 of file BeamFitter.h.

Referenced by BeamFitter(), and runFitterNoTxt().

◆ fz

double BeamFitter::fz
private

Definition at line 218 of file BeamFitter.h.

Referenced by BeamFitter(), and runFitterNoTxt().

◆ fz0

double BeamFitter::fz0
private

Definition at line 185 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ fzErr

double BeamFitter::fzErr
private

Definition at line 224 of file BeamFitter.h.

Referenced by BeamFitter(), and runFitterNoTxt().

◆ h1cutFlow

TH1F* BeamFitter::h1cutFlow
private

Definition at line 235 of file BeamFitter.h.

Referenced by BeamFitter(), getCutFlow(), readEvent(), resetCutFlow(), subtractFromCutFlow(), and ~BeamFitter().

◆ h1ntrks

TH1F* BeamFitter::h1ntrks
private

Definition at line 233 of file BeamFitter.h.

Referenced by BeamFitter(), readEvent(), and ~BeamFitter().

◆ h1vz_event

TH1F* BeamFitter::h1vz_event
private

Definition at line 234 of file BeamFitter.h.

Referenced by BeamFitter(), readEvent(), and ~BeamFitter().

◆ h1z

TH1F* BeamFitter::h1z
private

Definition at line 169 of file BeamFitter.h.

Referenced by runFitterNoTxt(), and ~BeamFitter().

◆ inputBeamWidth_

double BeamFitter::inputBeamWidth_
private

Definition at line 157 of file BeamFitter.h.

Referenced by BeamFitter(), runBeamWidthFitter(), runFitterNoTxt(), and runPVandTrkFitter().

◆ isMuon_

bool BeamFitter::isMuon_
private

Definition at line 161 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ min_Ntrks_

int BeamFitter::min_Ntrks_
private

Definition at line 160 of file BeamFitter.h.

Referenced by BeamFitter(), runBeamWidthFitter(), and runFitterNoTxt().

◆ MyPVFitter

PVFitter* BeamFitter::MyPVFitter
private

Definition at line 238 of file BeamFitter.h.

Referenced by BeamFitter(), dumpTxtFile(), getNPVs(), getNPVsperBX(), getPVvectorSize(), readEvent(), resetPVFitter(), resizePVvector(), runPVandTrkFitter(), and ~BeamFitter().

◆ outputDIPTxt_

std::string BeamFitter::outputDIPTxt_
private

Definition at line 145 of file BeamFitter.h.

Referenced by BeamFitter(), runFitter(), and runPVandTrkFitter().

◆ outputfilename_

std::string BeamFitter::outputfilename_
private

Definition at line 173 of file BeamFitter.h.

Referenced by BeamFitter().

◆ outputTxt_

std::string BeamFitter::outputTxt_
private

Definition at line 144 of file BeamFitter.h.

Referenced by BeamFitter(), dumpTxtFile(), runBeamWidthFitter(), runFitter(), and runPVandTrkFitter().

◆ quality_

std::vector<reco::TrackBase::TrackQuality> BeamFitter::quality_
private

Definition at line 155 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ saveBeamFit_

bool BeamFitter::saveBeamFit_
private

Definition at line 171 of file BeamFitter.h.

Referenced by BeamFitter(), runFitterNoTxt(), and ~BeamFitter().

◆ saveNtuple_

bool BeamFitter::saveNtuple_
private

Definition at line 170 of file BeamFitter.h.

Referenced by BeamFitter(), readEvent(), and ~BeamFitter().

◆ savePVVertices_

bool BeamFitter::savePVVertices_
private

Definition at line 172 of file BeamFitter.h.

Referenced by BeamFitter(), and ~BeamFitter().

◆ tracksToken_

edm::EDGetTokenT<reco::TrackCollection> BeamFitter::tracksToken_
private

Definition at line 138 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ trk_Algorithm_

std::vector<std::string> BeamFitter::trk_Algorithm_
private

Definition at line 153 of file BeamFitter.h.

Referenced by BeamFitter().

◆ trk_MaxEta_

double BeamFitter::trk_MaxEta_
private

Definition at line 148 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ trk_MaxIP_

double BeamFitter::trk_MaxIP_
private

Definition at line 149 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ trk_MaxNormChi2_

double BeamFitter::trk_MaxNormChi2_
private

Definition at line 152 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ trk_MaxZ_

double BeamFitter::trk_MaxZ_
private

Definition at line 147 of file BeamFitter.h.

Referenced by BeamFitter(), readEvent(), runBeamWidthFitter(), and runFitterNoTxt().

◆ trk_MinNPixLayers_

int BeamFitter::trk_MinNPixLayers_
private

Definition at line 151 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ trk_MinNTotLayers_

int BeamFitter::trk_MinNTotLayers_
private

Definition at line 150 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ trk_MinpT_

double BeamFitter::trk_MinpT_
private

Definition at line 146 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ trk_Quality_

std::vector<std::string> BeamFitter::trk_Quality_
private

Definition at line 154 of file BeamFitter.h.

Referenced by BeamFitter().

◆ vertexToken_

edm::EDGetTokenT<edm::View<reco::Vertex> > BeamFitter::vertexToken_
private

Definition at line 139 of file BeamFitter.h.

Referenced by BeamFitter(), and readEvent().

◆ writeDIPBadFit_

bool BeamFitter::writeDIPBadFit_
private

Definition at line 143 of file BeamFitter.h.

Referenced by BeamFitter(), runFitter(), and runPVandTrkFitter().

◆ writeDIPTxt_

bool BeamFitter::writeDIPTxt_
private

Definition at line 142 of file BeamFitter.h.

Referenced by BeamFitter(), runFitter(), and runPVandTrkFitter().

◆ writeTxt_

bool BeamFitter::writeTxt_
private

Definition at line 141 of file BeamFitter.h.

Referenced by BeamFitter(), dumpTxtFile(), runBeamWidthFitter(), runFitter(), and runPVandTrkFitter().

BeamFitter::beamSpotToken_
edm::EDGetTokenT< reco::BeamSpot > beamSpotToken_
Definition: BeamFitter.h:140
BeamFitter::fd0bs
double fd0bs
Definition: BeamFitter.h:183
PVFitter::resizepvStore
void resizepvStore(unsigned int rmSize)
Definition: PVFitter.h:112
BeamFitter::frunFit
int frunFit
Definition: BeamFitter.h:211
BeamFitter::fnTotLayerMeas
int fnTotLayerMeas
Definition: BeamFitter.h:187
beamspot::BeamSpotContainer::endTimeOfFit
char endTimeOfFit[32]
Definition: BeamSpotWrite2Txt.h:14
BSFitter::Fit_d0phi
reco::BeamSpot Fit_d0phi()
Definition: BSFitter.cc:410
BeamFitter::resetRefTime
void resetRefTime()
Definition: BeamFitter.h:61
PDWG_EXOHSCP_cff.tracks
tracks
Definition: PDWG_EXOHSCP_cff.py:28
BeamFitter::fwidthY
double fwidthY
Definition: BeamFitter.h:229
BeamFitter::trk_MaxIP_
double trk_MaxIP_
Definition: BeamFitter.h:149
mps_fire.i
i
Definition: mps_fire.py:428
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:95
reco::BeamSpot::setBeamWidthY
void setBeamWidthY(double v)
Definition: BeamSpot.h:105
HLT_FULL_cff.track
track
Definition: HLT_FULL_cff.py:11713
BeamFitter::quality_
std::vector< reco::TrackBase::TrackQuality > quality_
Definition: BeamFitter.h:155
reco::BeamSpot::emittanceX
double emittanceX() const
additional information
Definition: BeamSpot.h:129
BeamFitter::fz0
double fz0
Definition: BeamFitter.h:185
dqmMemoryStats.float
float
Definition: dqmMemoryStats.py:127
BeamFitter::updateBTime
void updateBTime()
Definition: BeamFitter.cc:29
reco::BeamSpot::BeamWidthXError
double BeamWidthXError() const
error on beam width X, assume error in X = Y
Definition: BeamSpot.h:99
reco::HitPattern::pixelEndcapLayersWithMeasurement
int pixelEndcapLayersWithMeasurement() const
Definition: HitPattern.cc:566
edm::Handle::product
T const * product() const
Definition: Handle.h:70
BeamFitter::fyErr
double fyErr
Definition: BeamFitter.h:223
BeamFitter::trk_MaxEta_
double trk_MaxEta_
Definition: BeamFitter.h:148
makeMuonMisalignmentScenario.matrix
list matrix
Definition: makeMuonMisalignmentScenario.py:141
cond::service::PoolDBOutputService::beginOfTime
cond::Time_t beginOfTime() const
Definition: PoolDBOutputService.cc:215
reco::HitPattern::stripTOBLayersWithMeasurement
int stripTOBLayersWithMeasurement() const
Definition: HitPattern.cc:597
filterCSVwithJSON.copy
copy
Definition: filterCSVwithJSON.py:36
BeamFitter::ftotal_tracks
int ftotal_tracks
Definition: BeamFitter.h:159
BeamFitter::vertexToken_
edm::EDGetTokenT< edm::View< reco::Vertex > > vertexToken_
Definition: BeamFitter.h:139
reco::BeamSpot::z0
double z0() const
z coordinate
Definition: BeamSpot.h:65
PVFitter::resetAll
void resetAll()
Definition: PVFitter.h:83
BeamFitter::fzErr
double fzErr
Definition: BeamFitter.h:224
BeamFitter::feta
double feta
Definition: BeamFitter.h:178
BeamFitter::fsigmaZ
double fsigmaZ
Definition: BeamFitter.h:219
gather_cfg.cout
cout
Definition: gather_cfg.py:144
BeamFitter::writeTxt_
bool writeTxt_
Definition: BeamFitter.h:141
HLT_FULL_cff.InputTag
InputTag
Definition: HLT_FULL_cff.py:89287
l1GtPatternGenerator_cfi.bx
bx
Definition: l1GtPatternGenerator_cfi.py:18
edm::ParameterSetDescription
Definition: ParameterSetDescription.h:52
BeamFitter::fasciiDIP
std::ofstream fasciiDIP
Definition: BeamFitter.h:134
cond::service::PoolDBOutputService::appendSinceTime
void appendSinceTime(const T *payloadObj, cond::Time_t sinceTime, const std::string &recordName)
Definition: PoolDBOutputService.h:141
BeamFitter::saveBeamFit_
bool saveBeamFit_
Definition: BeamFitter.h:171
BeamFitter::fpvx
double fpvx
Definition: BeamFitter.h:206
reco::BeamSpot::covariance
double covariance(int i, int j) const
(i,j)-th element of error matrix
Definition: BeamSpot.h:108
BeamFitter::fBSvector
std::vector< BSTrkParameters > fBSvector
Definition: BeamFitter.h:128
BeamFitter::fpt
double fpt
Definition: BeamFitter.h:177
BeamSpotObjects::SetBeamWidthY
void SetBeamWidthY(double val)
set average transverse beam width Y
Definition: BeamSpotObjects.h:50
BeamFitter::frun
int frun
Definition: BeamFitter.h:201
beamspot::BeamSpotContainer::reftime
std::time_t reftime[2]
Definition: BeamSpotWrite2Txt.h:17
BeamFitter::fxErr
double fxErr
Definition: BeamFitter.h:222
BeamFitter::h1vz_event
TH1F * h1vz_event
Definition: BeamFitter.h:234
BeamFitter::min_Ntrks_
int min_Ntrks_
Definition: BeamFitter.h:160
BeamFitter::falgo
bool falgo
Definition: BeamFitter.h:204
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
BeamFitter::fbeamWidthFit
reco::BeamSpot fbeamWidthFit
Definition: BeamFitter.h:130
BeamFitter::fdydz
double fdydz
Definition: BeamFitter.h:221
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
BeamFitter::trk_MinNTotLayers_
int trk_MinNTotLayers_
Definition: BeamFitter.h:150
MillePedeFileConverter_cfg.fileName
fileName
Definition: MillePedeFileConverter_cfg.py:32
reco::BeamSpot::dxdz
double dxdz() const
dxdz slope
Definition: BeamSpot.h:78
BSFitter::SetMaximumZ
void SetMaximumZ(double z)
Definition: BSFitter.h:55
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
reco::BeamSpot::dxdzError
double dxdzError() const
error on dxdz
Definition: BeamSpot.h:94
BeamFitter::fphi0
double fphi0
Definition: BeamFitter.h:181
BeamFitter::outputfilename_
std::string outputfilename_
Definition: BeamFitter.h:173
edm::Handle< reco::TrackCollection >
edm::Service::isAvailable
bool isAvailable() const
Definition: Service.h:40
BeamFitter::fbeginLumiOfFit
int fbeginLumiOfFit
Definition: BeamFitter.h:212
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
beamspot::BeamSpotContainer::beginLumiOfFit
int beginLumiOfFit
Definition: BeamSpotWrite2Txt.h:15
BeamFitter::fbeginTimeOfFit
char fbeginTimeOfFit[32]
Definition: BeamFitter.h:214
BeamFitter::flumi
int flumi
Definition: BeamFitter.h:202
reco::HitPattern::pixelLayersWithMeasurement
int pixelLayersWithMeasurement() const
Definition: HitPattern.cc:492
BeamSpotObjects::SetSigmaZ
void SetSigmaZ(double val)
set sigma Z, RMS bunch length
Definition: BeamSpotObjects.h:42
BeamFitter::ftree_
TTree * ftree_
Definition: BeamFitter.h:175
PVFitter::getNPVsperBX
const std::map< int, int > & getNPVsperBX()
Definition: PVFitter.h:122
BeamFitter::countPass
int countPass[9]
Definition: BeamFitter.h:236
reco::BeamSpot::y0Error
double y0Error() const
error on y
Definition: BeamSpot.h:88
reco::HitPattern::pixelBarrelLayersWithMeasurement
int pixelBarrelLayersWithMeasurement() const
Definition: HitPattern.cc:555
PVFitter::getpvStore
std::vector< BeamSpotFitPVData > getpvStore()
Definition: PVFitter.h:65
reco::BeamSpot::sigmaZ
double sigmaZ() const
sigma z
Definition: BeamSpot.h:76
BeamFitter::MyPVFitter
PVFitter * MyPVFitter
Definition: BeamFitter.h:238
BSFitter::SetInputBeamWidth
void SetInputBeamWidth(double val)
Definition: BSFitter.h:60
beamspot::BeamSpotContainer::beamspot
reco::BeamSpot beamspot
Definition: BeamSpotWrite2Txt.h:11
reco::HitPattern
Definition: HitPattern.h:147
BeamFitter::outputTxt_
std::string outputTxt_
Definition: BeamFitter.h:144
RandomServiceHelper.t1
t1
Definition: RandomServiceHelper.py:256
reco::HitPattern::stripLayersWithMeasurement
int stripLayersWithMeasurement() const
Definition: HitPattern.h:974
BeamFitter::file_
TFile * file_
Definition: BeamFitter.h:174
BSFitter::Fit
reco::BeamSpot Fit()
Definition: BSFitter.cc:102
beamspot::BeamSpotContainer::endLumiOfFit
int endLumiOfFit
Definition: BeamSpotWrite2Txt.h:16
BeamFitter::dumpBWTxtFile
void dumpBWTxtFile(std::string &)
Definition: BeamFitter.cc:656
beamspot::dumpBeamSpotTxt
void dumpBeamSpotTxt(std::ofstream &outFile, BeamSpotContainer const &bsContainer)
Definition: BeamSpotWrite2Txt.h:20
reco::HitPattern::trackerLayersWithMeasurement
int trackerLayersWithMeasurement() const
Definition: HitPattern.cc:513
BeamFitter::saveNtuple_
bool saveNtuple_
Definition: BeamFitter.h:170
BeamFitter::convergence_
double convergence_
Definition: BeamFitter.h:158
reco::BeamSpot::dydzError
double dydzError() const
error on dydz
Definition: BeamSpot.h:96
FrontierCondition_GT_autoExpress_cfi.t0
t0
Definition: FrontierCondition_GT_autoExpress_cfi.py:148
cond::service::PoolDBOutputService::createNewIOV
void createNewIOV(const T *firstPayloadObj, cond::Time_t firstSinceTime, cond::Time_t, const std::string &recordName)
Definition: PoolDBOutputService.h:116
BeamFitter::fbeamspot
reco::BeamSpot fbeamspot
Definition: BeamFitter.h:129
BeamFitter::ffilename_changed
bool ffilename_changed
Definition: BeamFitter.h:163
reco::BeamSpot::type
BeamType type() const
return beam type
Definition: BeamSpot.h:122
L1TEGammaOffline_cfi.PVCollection
PVCollection
Definition: L1TEGammaOffline_cfi.py:34
BeamFitter::fd0
double fd0
Definition: BeamFitter.h:182
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
PVFitter::readEvent
void readEvent(const edm::Event &iEvent)
Definition: PVFitter.cc:120
PVFitter::getBeamSpotMap
std::map< int, reco::BeamSpot > getBeamSpotMap()
Definition: PVFitter.h:95
reco::BeamSpot
Definition: BeamSpot.h:21
BeamFitter::fcharge
int fcharge
Definition: BeamFitter.h:179
BeamFitter::fquality
bool fquality
Definition: BeamFitter.h:203
BeamFitter::fvy
double fvy
Definition: BeamFitter.h:200
BeamFitter::h1z
TH1F * h1z
Definition: BeamFitter.h:169
BeamFitter::trk_Quality_
std::vector< std::string > trk_Quality_
Definition: BeamFitter.h:154
BeamFitter::savePVVertices_
bool savePVVertices_
Definition: BeamFitter.h:172
edm::ConsumesCollector::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: ConsumesCollector.h:55
reco::BeamSpot::Fake
Definition: BeamSpot.h:24
PVFitter::getWidthXerr
double getWidthXerr()
Definition: PVFitter.h:61
dqmdumpme.k
k
Definition: dqmdumpme.py:60
reco::BeamSpot::position
const Point & position() const
position
Definition: BeamSpot.h:59
BeamFitter::fnPXBLayerMeas
int fnPXBLayerMeas
Definition: BeamFitter.h:194
BSFitter::Fit_ited0phi
reco::BeamSpot Fit_ited0phi()
Definition: BSFitter.cc:360
BeamFitter::outputDIPTxt_
std::string outputDIPTxt_
Definition: BeamFitter.h:145
BeamSpotObjects::SetPosition
void SetPosition(double x, double y, double z)
set XYZ position
Definition: BeamSpotObjects.h:36
L1TdeCSCTF_cfi.outFile
outFile
Definition: L1TdeCSCTF_cfi.py:5
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
reco::BeamSpot::Unknown
Definition: BeamSpot.h:24
BeamFitter::fx
double fx
Definition: BeamFitter.h:216
reco::HitPattern::stripTIBLayersWithMeasurement
int stripTIBLayersWithMeasurement() const
Definition: HitPattern.cc:577
BeamFitter::appendRunTxt_
bool appendRunTxt_
Definition: BeamFitter.h:137
AlcaBeamMonitor_cfi.PVFitter
PVFitter
Definition: AlcaBeamMonitor_cfi.py:37
BeamFitter::fitted_
bool fitted_
Definition: BeamFitter.h:162
BeamFitter::fnStripLayerMeas
int fnStripLayerMeas
Definition: BeamFitter.h:189
edm::View
Definition: CaloClusterFwd.h:14
BeamFitter::writeDIPTxt_
bool writeDIPTxt_
Definition: BeamFitter.h:142
edm::ParameterSetDescription::addUntracked
ParameterDescriptionBase * addUntracked(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:100
BeamFitter::h1ntrks
TH1F * h1ntrks
Definition: BeamFitter.h:233
edm::ParameterSet
Definition: ParameterSet.h:47
BeamFitter::fdydzErr
double fdydzErr
Definition: BeamFitter.h:227
PVFitter::getNPVs
int getNPVs()
Definition: PVFitter.h:120
reco::BeamSpot::betaStar
double betaStar() const
Definition: BeamSpot.h:131
BeamFitter::fPVTree_
TTree * fPVTree_
Definition: BeamFitter.h:239
BSFitter::SetConvergence
void SetConvergence(double val)
Definition: BSFitter.h:56
BeamSpotObjects::SetCovariance
void SetCovariance(int i, int j, double val)
set i,j element of the full covariance matrix 7x7
Definition: BeamSpotObjects.h:56
BeamFitter::fpvy
double fpvy
Definition: BeamFitter.h:206
BeamFitter::trk_MaxNormChi2_
double trk_MaxNormChi2_
Definition: BeamFitter.h:152
mps_setup.append
append
Definition: mps_setup.py:85
BeamFitter::fsigmaZErr
double fsigmaZErr
Definition: BeamFitter.h:225
BeamFitter::fwidthX
double fwidthX
Definition: BeamFitter.h:228
BeamFitter::fz
double fz
Definition: BeamFitter.h:218
BeamFitter::fdxdz
double fdxdz
Definition: BeamFitter.h:220
BeamFitter::freftime
std::time_t freftime[2]
Definition: BeamFitter.h:207
edm::Service< cond::service::PoolDBOutputService >
createfilelist.int
int
Definition: createfilelist.py:10
MetAnalyzer.pv
def pv(vc)
Definition: MetAnalyzer.py:7
iEvent
int iEvent
Definition: GenABIO.cc:224
BeamFitter::algorithm_
std::vector< reco::TrackBase::TrackAlgorithm > algorithm_
Definition: BeamFitter.h:156
PVFitter::runBXFitter
bool runBXFitter()
Definition: PVFitter.cc:211
reco::BeamSpot::z0Error
double z0Error() const
error on z
Definition: BeamSpot.h:90
BSFitter::Setd0Cut_d0phi
void Setd0Cut_d0phi(double d0cut)
Definition: BSFitter.cc:576
BeamFitter::runFitterNoTxt
bool runFitterNoTxt()
Definition: BeamFitter.cc:546
BeamSpotObjects
Definition: BeamSpotObjects.h:20
BeamFitter::dumpTxtFile
void dumpTxtFile(std::string &, bool)
Definition: BeamFitter.cc:668
reco::BeamSpot::BeamWidthYError
double BeamWidthYError() const
error on beam width Y, assume error in X = Y
Definition: BeamSpot.h:101
beamspot::BeamSpotContainer::run
int run
Definition: BeamSpotWrite2Txt.h:12
reco::BeamSpot::CovarianceMatrix
math::Error< dimension >::type CovarianceMatrix
Definition: BeamSpot.h:29
reco::BeamSpot::BeamWidthX
double BeamWidthX() const
beam width X
Definition: BeamSpot.h:82
BSTrkParameters
Definition: BSTrkParameters.h:15
BSFitter::SetMinimumNTrks
void SetMinimumNTrks(int n)
Definition: BSFitter.h:57
BeamSpotObjects::Setdxdz
void Setdxdz(double val)
set dxdz slope, crossing angle
Definition: BeamSpotObjects.h:44
BeamFitter::fnTIDLayerMeas
int fnTIDLayerMeas
Definition: BeamFitter.h:191
BeamFitter::resetCutFlow
void resetCutFlow()
Definition: BeamFitter.h:104
BeamFitter::tracksToken_
edm::EDGetTokenT< reco::TrackCollection > tracksToken_
Definition: BeamFitter.h:138
BeamFitter::trk_Algorithm_
std::vector< std::string > trk_Algorithm_
Definition: BeamFitter.h:153
BeamFitter::fsigmad0
double fsigmad0
Definition: BeamFitter.h:184
BeamFitter::fnTIBLayerMeas
int fnTIBLayerMeas
Definition: BeamFitter.h:190
BeamFitter::fnormchi2
double fnormchi2
Definition: BeamFitter.h:180
reco::BeamSpot::dydz
double dydz() const
dydz slope
Definition: BeamSpot.h:80
cond::service::PoolDBOutputService::endOfTime
cond::Time_t endOfTime() const
Definition: PoolDBOutputService.cc:213
BeamFitter::fy
double fy
Definition: BeamFitter.h:217
reco::TrackBase::qualityByName
static TrackQuality qualityByName(const std::string &name)
Definition: TrackBase.cc:126
BeamFitter::h1cutFlow
TH1F * h1cutFlow
Definition: BeamFitter.h:235
BeamFitter::debug_
bool debug_
Definition: BeamFitter.h:136
reco::BeamSpot::x0
double x0() const
x coordinate
Definition: BeamSpot.h:61
BeamFitter::ftheta
double ftheta
Definition: BeamFitter.h:176
BeamFitter::fpvValid
bool fpvValid
Definition: BeamFitter.h:205
PVFitter::IsFitPerBunchCrossing
bool IsFitPerBunchCrossing()
Definition: PVFitter.h:96
BeamFitter::fnPixelLayerMeas
int fnPixelLayerMeas
Definition: BeamFitter.h:188
BeamFitter::isMuon_
bool isMuon_
Definition: BeamFitter.h:161
BeamFitter::fnTOBLayerMeas
int fnTOBLayerMeas
Definition: BeamFitter.h:192
reco::HitPattern::stripTIDLayersWithMeasurement
int stripTIDLayersWithMeasurement() const
Definition: HitPattern.cc:587
writedatasetfile.run
run
Definition: writedatasetfile.py:27
BeamFitter::fvx
double fvx
Definition: BeamFitter.h:199
BeamFitter::fcov
double fcov[7][7]
Definition: BeamFitter.h:198
BeamFitter::fsigmaz0
double fsigmaz0
Definition: BeamFitter.h:186
BeamFitter::trk_MaxZ_
double trk_MaxZ_
Definition: BeamFitter.h:147
BeamFitter::fendTimeOfFit
char fendTimeOfFit[32]
Definition: BeamFitter.h:215
reco::BeamSpot::x0Error
double x0Error() const
error on x
Definition: BeamSpot.h:86
reco::BeamSpot::sigmaZ0Error
double sigmaZ0Error() const
error on sigma z
Definition: BeamSpot.h:92
BeamFitter::trk_MinNPixLayers_
int trk_MinNPixLayers_
Definition: BeamFitter.h:151
BeamFitter::fnPXFLayerMeas
int fnPXFLayerMeas
Definition: BeamFitter.h:195
BeamSpotObjects::SetBeamWidthX
void SetBeamWidthX(double val)
set average transverse beam width X
Definition: BeamSpotObjects.h:48
BeamFitter::fnTECLayerMeas
int fnTECLayerMeas
Definition: BeamFitter.h:193
BSFitter::SetFitVariable
void SetFitVariable(std::string name)
Definition: BSFitter.h:41
BeamFitter::ForDIPPV_
std::vector< float > ForDIPPV_
Definition: BeamFitter.h:166
cond::service::PoolDBOutputService::isNewTagRequest
bool isNewTagRequest(const std::string &recordName)
Definition: PoolDBOutputService.cc:128
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
BeamFitter::fwidthYErr
double fwidthYErr
Definition: BeamFitter.h:231
AlignmentPI::index
index
Definition: AlignmentPayloadInspectorHelper.h:46
beamspot::BeamSpotContainer::beginTimeOfFit
char beginTimeOfFit[32]
Definition: BeamSpotWrite2Txt.h:13
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:29
beamspot::BeamSpotContainer
Definition: BeamSpotWrite2Txt.h:10
BeamFitter::fdxdzErr
double fdxdzErr
Definition: BeamFitter.h:226
BSFitter
Definition: BSFitter.h:30
BSFitter::SetFitType
void SetFitType(std::string type)
Definition: BSFitter.h:39
reco::TrackBase::algoByName
static TrackAlgorithm algoByName(const std::string &name)
Definition: TrackBase.cc:137
reco::BeamSpot::Point
math::XYZPoint Point
point in the space
Definition: BeamSpot.h:27
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
reco::HitPattern::stripTECLayersWithMeasurement
int stripTECLayersWithMeasurement() const
Definition: HitPattern.cc:607
PVFitter::runFitter
bool runFitter()
Definition: PVFitter.cc:330
reco::BeamSpot::setType
void setType(BeamType type)
set beam type
Definition: BeamSpot.h:124
dqmiolumiharvest.j
j
Definition: dqmiolumiharvest.py:66
BeamFitter::trk_MinpT_
double trk_MinpT_
Definition: BeamFitter.h:146
reco::BeamSpot::y0
double y0() const
y coordinate
Definition: BeamSpot.h:63
BeamFitter::fbspotPVMap
std::map< int, reco::BeamSpot > fbspotPVMap
Definition: BeamFitter.h:131
BeamSpotObjects::Setdydz
void Setdydz(double val)
set dydz slope, crossing angle in XZ
Definition: BeamSpotObjects.h:46
reco::BeamSpot::emittanceY
double emittanceY() const
Definition: BeamSpot.h:130
BeamFitter::fendLumiOfFit
int fendLumiOfFit
Definition: BeamFitter.h:213
BeamFitter::inputBeamWidth_
double inputBeamWidth_
Definition: BeamFitter.h:157
edm::InputTag
Definition: InputTag.h:15
reco::TrackCollection
std::vector< Track > TrackCollection
collection of Tracks
Definition: TrackFwd.h:14
cond::service::PoolDBOutputService::currentTime
cond::Time_t currentTime() const
Definition: PoolDBOutputService.cc:217
PVFitter::getBeamSpot
reco::BeamSpot getBeamSpot()
Definition: PVFitter.h:94
BeamFitter::fwidthXErr
double fwidthXErr
Definition: BeamFitter.h:230
reco::BeamSpot::BeamWidthY
double BeamWidthY() const
beam width Y
Definition: BeamSpot.h:84
edm::TimeValue_t
unsigned long long TimeValue_t
Definition: Timestamp.h:28
BeamFitter::fpvz
double fpvz
Definition: BeamFitter.h:206
BSFitter::GetVzHisto
TH1F * GetVzHisto()
Definition: BSFitter.h:87
BeamFitter::ftreeFit_
TTree * ftreeFit_
Definition: BeamFitter.h:210
BeamFitter::writeDIPBadFit_
bool writeDIPBadFit_
Definition: BeamFitter.h:143
findQualityFiles.size
size
Write out results.
Definition: findQualityFiles.py:443
PVFitter::setTree
void setTree(TTree *tree)
Definition: PVFitter.cc:206
Generated for CMSSW Reference Manual by doxygen 1.8.16