#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.

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_2022v12_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_.

                                                                                    : fPVTree_(nullptr) {
   debug_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<bool>("Debug");
   tracksToken_ = iColl.consumes<reco::TrackCollection>(
       iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<edm::InputTag>("TrackCollection"));
   vertexToken_ = iColl.consumes<edm::View<reco::Vertex> >(
       iConfig.getUntrackedParameter<edm::InputTag>("primaryVertex", edm::InputTag("offlinePrimaryVertices")));
   beamSpotToken_ = iColl.consumes<reco::BeamSpot>(
       iConfig.getUntrackedParameter<edm::InputTag>("beamSpot", edm::InputTag("offlineBeamSpot")));
   writeTxt_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<bool>("WriteAscii");
   outputTxt_ =
       iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<std::string>("AsciiFileName");
   appendRunTxt_ =
       iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<bool>("AppendRunToFileName");
   writeDIPTxt_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<bool>("WriteDIPAscii");
   // Specify whether we want to write the DIP file even if the fit is failed.
   writeDIPBadFit_ =
       iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<bool>("WriteDIPOnBadFit", true);
   outputDIPTxt_ =
       iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<std::string>("DIPFileName");
   saveNtuple_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<bool>("SaveNtuple");
   saveBeamFit_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<bool>("SaveFitResults");
   savePVVertices_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<bool>("SavePVVertices");
   isMuon_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<bool>("IsMuonCollection");
 
   trk_MinpT_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<double>("MinimumPt");
   trk_MaxEta_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<double>("MaximumEta");
   trk_MaxIP_ =
       iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<double>("MaximumImpactParameter");
   trk_MaxZ_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<double>("MaximumZ");
   trk_MinNTotLayers_ =
       iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<int>("MinimumTotalLayers");
   trk_MinNPixLayers_ =
       iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<int>("MinimumPixelLayers");
   trk_MaxNormChi2_ =
       iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<double>("MaximumNormChi2");
   trk_Algorithm_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter")
                        .getUntrackedParameter<std::vector<std::string> >("TrackAlgorithm");
   trk_Quality_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter")
                      .getUntrackedParameter<std::vector<std::string> >("TrackQuality");
   min_Ntrks_ = iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<int>("MinimumInputTracks");
   convergence_ =
       iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<double>("FractionOfFittedTrks");
   inputBeamWidth_ =
       iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<double>("InputBeamWidth", -1.);
 
   for (unsigned int j = 0; j < trk_Algorithm_.size(); j++)
     algorithm_.push_back(reco::TrackBase::algoByName(trk_Algorithm_[j]));
   for (unsigned int j = 0; j < trk_Quality_.size(); j++)
     quality_.push_back(reco::TrackBase::qualityByName(trk_Quality_[j]));
 
   if (saveNtuple_ || saveBeamFit_ || savePVVertices_) {
     outputfilename_ =
         iConfig.getParameter<edm::ParameterSet>("BeamFitter").getUntrackedParameter<std::string>("OutputFileName");
     file_ = TFile::Open(outputfilename_.c_str(), "RECREATE");
   }
   if (saveNtuple_) {
     ftree_ = new TTree("mytree", "mytree");
     ftree_->AutoSave();
 
     ftree_->Branch("pt", &fpt, "fpt/D");
     ftree_->Branch("d0", &fd0, "fd0/D");
     ftree_->Branch("d0bs", &fd0bs, "fd0bs/D");
     ftree_->Branch("sigmad0", &fsigmad0, "fsigmad0/D");
     ftree_->Branch("phi0", &fphi0, "fphi0/D");
     ftree_->Branch("z0", &fz0, "fz0/D");
     ftree_->Branch("sigmaz0", &fsigmaz0, "fsigmaz0/D");
     ftree_->Branch("theta", &ftheta, "ftheta/D");
     ftree_->Branch("eta", &feta, "feta/D");
     ftree_->Branch("charge", &fcharge, "fcharge/I");
     ftree_->Branch("normchi2", &fnormchi2, "fnormchi2/D");
     ftree_->Branch("nTotLayerMeas", &fnTotLayerMeas, "fnTotLayerMeas/i");
     ftree_->Branch("nStripLayerMeas", &fnStripLayerMeas, "fnStripLayerMeas/i");
     ftree_->Branch("nPixelLayerMeas", &fnPixelLayerMeas, "fnPixelLayerMeas/i");
     ftree_->Branch("nTIBLayerMeas", &fnTIBLayerMeas, "fnTIBLayerMeas/i");
     ftree_->Branch("nTOBLayerMeas", &fnTOBLayerMeas, "fnTOBLayerMeas/i");
     ftree_->Branch("nTIDLayerMeas", &fnTIDLayerMeas, "fnTIDLayerMeas/i");
     ftree_->Branch("nTECLayerMeas", &fnTECLayerMeas, "fnTECLayerMeas/i");
     ftree_->Branch("nPXBLayerMeas", &fnPXBLayerMeas, "fnPXBLayerMeas/i");
     ftree_->Branch("nPXFLayerMeas", &fnPXFLayerMeas, "fnPXFLayerMeas/i");
     ftree_->Branch("cov", &fcov, "fcov[7][7]/D");
     ftree_->Branch("vx", &fvx, "fvx/D");
     ftree_->Branch("vy", &fvy, "fvy/D");
     ftree_->Branch("quality", &fquality, "fquality/O");
     ftree_->Branch("algo", &falgo, "falgo/O");
     ftree_->Branch("run", &frun, "frun/i");
     ftree_->Branch("lumi", &flumi, "flumi/i");
     ftree_->Branch("pvValid", &fpvValid, "fpvValid/O");
     ftree_->Branch("pvx", &fpvx, "fpvx/D");
     ftree_->Branch("pvy", &fpvy, "fpvy/D");
     ftree_->Branch("pvz", &fpvz, "fpvz/D");
   }
   if (saveBeamFit_) {
     ftreeFit_ = new TTree("fitResults", "fitResults");
     ftreeFit_->AutoSave();
     ftreeFit_->Branch("run", &frunFit, "frunFit/i");
     ftreeFit_->Branch("beginLumi", &fbeginLumiOfFit, "fbeginLumiOfFit/i");
     ftreeFit_->Branch("endLumi", &fendLumiOfFit, "fendLumiOfFit/i");
     ftreeFit_->Branch("beginTime", fbeginTimeOfFit, "fbeginTimeOfFit/C");
     ftreeFit_->Branch("endTime", fendTimeOfFit, "fendTimeOfFit/C");
     ftreeFit_->Branch("x", &fx, "fx/D");
     ftreeFit_->Branch("y", &fy, "fy/D");
     ftreeFit_->Branch("z", &fz, "fz/D");
     ftreeFit_->Branch("sigmaZ", &fsigmaZ, "fsigmaZ/D");
     ftreeFit_->Branch("dxdz", &fdxdz, "fdxdz/D");
     ftreeFit_->Branch("dydz", &fdydz, "fdydz/D");
     ftreeFit_->Branch("xErr", &fxErr, "fxErr/D");
     ftreeFit_->Branch("yErr", &fyErr, "fyErr/D");
     ftreeFit_->Branch("zErr", &fzErr, "fzErr/D");
     ftreeFit_->Branch("sigmaZErr", &fsigmaZErr, "fsigmaZErr/D");
     ftreeFit_->Branch("dxdzErr", &fdxdzErr, "fdxdzErr/D");
     ftreeFit_->Branch("dydzErr", &fdydzErr, "fdydzErr/D");
     ftreeFit_->Branch("widthX", &fwidthX, "fwidthX/D");
     ftreeFit_->Branch("widthY", &fwidthY, "fwidthY/D");
     ftreeFit_->Branch("widthXErr", &fwidthXErr, "fwidthXErr/D");
     ftreeFit_->Branch("widthYErr", &fwidthYErr, "fwidthYErr/D");
   }
 
   fBSvector.clear();
   ftotal_tracks = 0;
   fnTotLayerMeas = fnPixelLayerMeas = fnStripLayerMeas = fnTIBLayerMeas = 0;
   fnTIDLayerMeas = fnTOBLayerMeas = fnTECLayerMeas = fnPXBLayerMeas = fnPXFLayerMeas = 0;
   frun = flumi = -1;
   frunFit = fbeginLumiOfFit = fendLumiOfFit = -1;
   fquality = falgo = true;
   fpvValid = true;
   fpvx = fpvy = fpvz = 0;
   fitted_ = false;
   resetRefTime();
 
   //debug histograms
   h1ntrks = new TH1F("h1ntrks", "number of tracks per event", 50, 0, 50);
   h1vz_event = new TH1F("h1vz_event", "track Vz", 50, -30, 30);
   h1cutFlow = new TH1F("h1cutFlow", "Cut flow table of track selection", 9, 0, 9);
   h1cutFlow->GetXaxis()->SetBinLabel(1, "No cut");
   h1cutFlow->GetXaxis()->SetBinLabel(2, "Traker hits");
   h1cutFlow->GetXaxis()->SetBinLabel(3, "Pixel hits");
   h1cutFlow->GetXaxis()->SetBinLabel(4, "norm. #chi^{2}");
   h1cutFlow->GetXaxis()->SetBinLabel(5, "algo");
   h1cutFlow->GetXaxis()->SetBinLabel(6, "quality");
   h1cutFlow->GetXaxis()->SetBinLabel(7, "d_{0}");
   h1cutFlow->GetXaxis()->SetBinLabel(8, "z_{0}");
   h1cutFlow->GetXaxis()->SetBinLabel(9, "p_{T}");
   resetCutFlow();
 
   // Primary vertex fitter
   MyPVFitter = new PVFitter(iConfig, iColl);
   MyPVFitter->resetAll();
   if (savePVVertices_) {
     fPVTree_ = new TTree("PrimaryVertices", "PrimaryVertices");
     MyPVFitter->setTree(fPVTree_);
   }
 
   // check filename
   ffilename_changed = false;
   if (writeDIPTxt_)
     fasciiDIP.open(outputDIPTxt_.c_str());
 }

◆ ~BeamFitter()

BeamFitter::~BeamFitter ( )

virtual

Definition at line 197 of file BeamFitter.cc.

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

                         {
   if (saveNtuple_) {
     file_->cd();
     if (fitted_ && h1z)
       h1z->Write();
     h1ntrks->Write();
     h1vz_event->Write();
     if (h1cutFlow)
       h1cutFlow->Write();
     ftree_->Write();
   }
   if (saveBeamFit_) {
     file_->cd();
     ftreeFit_->Write();
   }
   if (savePVVertices_) {
     file_->cd();
     fPVTree_->Write();
   }
 
   if (saveNtuple_ || saveBeamFit_ || savePVVertices_) {
     file_->Close();
     delete file_;
   }
   delete MyPVFitter;
 }

Member Function Documentation

◆ dumpBWTxtFile()

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

Definition at line 656 of file BeamFitter.cc.

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

Referenced by runBeamWidthFitter().

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

◆ dumpTxtFile()

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

Definition at line 668 of file BeamFitter.cc.

References mps_setup::append, appendRunTxt_, beamspot::BeamSpotContainer::beamspot, reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), beamspot::BeamSpotContainer::beginLumiOfFit, beamspot::BeamSpotContainer::beginTimeOfFit, reco::BeamSpot::betaStar(), simKBmtfDigis_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().

                                                              {
   std::ofstream outFile;
 
   std::string tmpname = outputTxt_;
   char index[15];
   if (appendRunTxt_ && writeTxt_ && !ffilename_changed) {
     sprintf(index, "%s%i", "_Run", frun);
     tmpname.insert(outputTxt_.length() - 4, index);
     fileName = tmpname;
     ffilename_changed = true;
   }
 
   if (!append)
     outFile.open(fileName.c_str());
   else
     outFile.open(fileName.c_str(), std::ios::app);
 
   if (MyPVFitter->IsFitPerBunchCrossing()) {
     for (std::map<int, reco::BeamSpot>::const_iterator abspot = fbspotPVMap.begin(); abspot != fbspotPVMap.end();
          ++abspot) {
       reco::BeamSpot beamspottmp = abspot->second;
       int bx = abspot->first;
 
       outFile << "Runnumber " << frun << " bx " << bx << std::endl;
       outFile << "BeginTimeOfFit " << fbeginTimeOfFit << " " << freftime[0] << std::endl;
       outFile << "EndTimeOfFit " << fendTimeOfFit << " " << freftime[1] << std::endl;
       outFile << "LumiRange " << fbeginLumiOfFit << " - " << fendLumiOfFit << std::endl;
       outFile << "Type " << beamspottmp.type() << std::endl;
       outFile << "X0 " << beamspottmp.x0() << std::endl;
       outFile << "Y0 " << beamspottmp.y0() << std::endl;
       outFile << "Z0 " << beamspottmp.z0() << std::endl;
       outFile << "sigmaZ0 " << beamspottmp.sigmaZ() << std::endl;
       outFile << "dxdz " << beamspottmp.dxdz() << std::endl;
       outFile << "dydz " << beamspottmp.dydz() << std::endl;
       outFile << "BeamWidthX " << beamspottmp.BeamWidthX() << std::endl;
       outFile << "BeamWidthY " << beamspottmp.BeamWidthY() << std::endl;
       for (int i = 0; i < 6; ++i) {
         outFile << "Cov(" << i << ",j) ";
         for (int j = 0; j < 7; ++j) {
           outFile << beamspottmp.covariance(i, j) << " ";
         }
         outFile << std::endl;
       }
       outFile << "Cov(6,j) 0 0 0 0 0 0 " << beamspottmp.covariance(6, 6) << std::endl;
       //}
       outFile << "EmittanceX " << beamspottmp.emittanceX() << std::endl;
       outFile << "EmittanceY " << beamspottmp.emittanceY() << std::endl;
       outFile << "BetaStar " << beamspottmp.betaStar() << std::endl;
     }
   }  //if bx results needed
   else {
     beamspot::BeamSpotContainer currentBS;
 
     currentBS.beamspot = fbeamspot;
     currentBS.run = frun;
     std::copy(fbeginTimeOfFit, fbeginTimeOfFit + 32, currentBS.beginTimeOfFit);
     std::copy(fendTimeOfFit, fendTimeOfFit + 32, currentBS.endTimeOfFit);
     currentBS.beginLumiOfFit = fbeginLumiOfFit;
     currentBS.endLumiOfFit = fendLumiOfFit;
     std::copy(freftime, freftime + 2, currentBS.reftime);
 
     beamspot::dumpBeamSpotTxt(outFile, currentBS);
 
     //write here Pv info for DIP only: This added only if append is false, which happen for DIP only :)
     if (!append) {
       outFile << "events " << (int)ForDIPPV_[0] << std::endl;
       outFile << "meanPV " << ForDIPPV_[1] << std::endl;
       outFile << "meanErrPV " << ForDIPPV_[2] << std::endl;
       outFile << "rmsPV " << ForDIPPV_[3] << std::endl;
       outFile << "rmsErrPV " << ForDIPPV_[4] << std::endl;
       outFile << "maxPV " << (int)ForDIPPV_[5] << std::endl;
       outFile << "nPV " << (int)ForDIPPV_[6] << std::endl;
     }  //writeDIPPVInfo_
   }    //else end  here
 
   outFile.close();
 }

◆ fillDescription()

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

static

Definition at line 224 of file BeamFitter.cc.

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

Referenced by BeamSpotAnalyzer::fillDescriptions(), and AlcaBeamMonitor::fillDescriptions().

                                                                   {
   edm::ParameterSetDescription beamFitter;
 
   beamFitter.addUntracked<bool>("Debug");
   beamFitter.addUntracked<edm::InputTag>("TrackCollection");
   iDesc.addUntracked<edm::InputTag>("primaryVertex", edm::InputTag("offlinePrimaryVertices"));
   iDesc.addUntracked<edm::InputTag>("beamSpot", edm::InputTag("offlineBeamSpot"));
   beamFitter.addUntracked<bool>("WriteAscii");
   beamFitter.addUntracked<std::string>("AsciiFileName");
   beamFitter.addUntracked<bool>("AppendRunToFileName");
   beamFitter.addUntracked<bool>("WriteDIPAscii");
   // Specify whether we want to write the DIP file even if the fit is failed.
   beamFitter.addUntracked<bool>("WriteDIPOnBadFit", true);
   beamFitter.addUntracked<std::string>("DIPFileName");
   beamFitter.addUntracked<bool>("SaveNtuple");
   beamFitter.addUntracked<bool>("SaveFitResults");
   beamFitter.addUntracked<bool>("SavePVVertices");
   beamFitter.addUntracked<bool>("IsMuonCollection");
 
   beamFitter.addUntracked<double>("MinimumPt");
   beamFitter.addUntracked<double>("MaximumEta");
   beamFitter.addUntracked<double>("MaximumImpactParameter");
   beamFitter.addUntracked<double>("MaximumZ");
   beamFitter.addUntracked<int>("MinimumTotalLayers");
   beamFitter.addUntracked<int>("MinimumPixelLayers");
   beamFitter.addUntracked<double>("MaximumNormChi2");
   beamFitter.addUntracked<std::vector<std::string> >("TrackAlgorithm");
   beamFitter.addUntracked<std::vector<std::string> >("TrackQuality");
   beamFitter.addUntracked<int>("MinimumInputTracks");
   beamFitter.addUntracked<double>("FractionOfFittedTrks");
   beamFitter.addUntracked<double>("InputBeamWidth", -1.);
 
   beamFitter.addUntracked<std::string>("OutputFileName", "");
 
   iDesc.add<edm::ParameterSetDescription>("BeamFitter", beamFitter);
 }

◆ getBeamSpot()

reco::BeamSpot BeamFitter::getBeamSpot ( )

inline

Definition at line 93 of file BeamFitter.h.

References fbeamspot.

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

93 { return fbeamspot; }

BeamFitter::fbeamspot

reco::BeamSpot fbeamspot

Definition: BeamFitter.h:129

◆ getBeamSpotMap()

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

inline

Definition at line 94 of file BeamFitter.h.

References fbspotPVMap.

Referenced by BeamMonitorBx::FitAndFill().

94 { return fbspotPVMap; }

BeamFitter::fbspotPVMap

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

Definition: BeamFitter.h:131

◆ getBeamWidth()

reco::BeamSpot BeamFitter::getBeamWidth ( )

inline

Definition at line 56 of file BeamFitter.h.

References fbeamWidthFit.

Referenced by BeamSpotAnalyzer::endJob().

56 { return fbeamWidthFit; }

BeamFitter::fbeamWidthFit

reco::BeamSpot fbeamWidthFit

Definition: BeamFitter.h:130

◆ getBSvector()

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

inline

Definition at line 95 of file BeamFitter.h.

References fBSvector.

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

95 { return fBSvector; }

BeamFitter::fBSvector

std::vector< BSTrkParameters > fBSvector

Definition: BeamFitter.h:128

◆ getCutFlow()

TH1F* BeamFitter::getCutFlow ( )

inline

Definition at line 96 of file BeamFitter.h.

References h1cutFlow.

96 { return h1cutFlow; }

BeamFitter::h1cutFlow

TH1F * h1cutFlow

Definition: BeamFitter.h:235

◆ getFitLSRange()

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

inline

Definition at line 114 of file BeamFitter.h.

References fbeginLumiOfFit, and fendLumiOfFit.

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

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

BeamFitter::fendLumiOfFit

int fendLumiOfFit

Definition: BeamFitter.h:213

BeamFitter::fbeginLumiOfFit

int fbeginLumiOfFit

Definition: BeamFitter.h:212

◆ getNPVs()

int BeamFitter::getNPVs ( )

inline

Definition at line 122 of file BeamFitter.h.

References PVFitter::getNPVs(), and MyPVFitter.

Referenced by BeamSpotAnalyzer::endLuminosityBlock().

122 { return MyPVFitter->getNPVs(); }

BeamFitter::MyPVFitter

PVFitter * MyPVFitter

Definition: BeamFitter.h:238

PVFitter::getNPVs

int getNPVs()

Definition: PVFitter.h:120

◆ getNPVsperBX()

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

inline

Definition at line 123 of file BeamFitter.h.

References PVFitter::getNPVsperBX(), and MyPVFitter.

Referenced by BeamMonitorBx::FitAndFill().

123 { return MyPVFitter->getNPVsperBX(); }

PVFitter::getNPVsperBX

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

Definition: PVFitter.h:122

BeamFitter::MyPVFitter

PVFitter * MyPVFitter

Definition: BeamFitter.h:238

◆ getNTracks()

int BeamFitter::getNTracks ( )

inline

Definition at line 121 of file BeamFitter.h.

References fBSvector.

Referenced by BeamSpotAnalyzer::endLuminosityBlock().

121 { return fBSvector.size(); }

BeamFitter::fBSvector

std::vector< BSTrkParameters > fBSvector

Definition: BeamFitter.h:128

◆ getPVvectorSize()

std::size_t BeamFitter::getPVvectorSize ( )

inline

Definition at line 75 of file BeamFitter.h.

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

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

findQualityFiles.size

size

Write out results.

Definition: findQualityFiles.py:443

BeamFitter::MyPVFitter

PVFitter * MyPVFitter

Definition: BeamFitter.h:238

PVFitter::getpvStore

std::vector< BeamSpotFitPVData > getpvStore()

Definition: PVFitter.h:65

◆ getRefTime()

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

inline

Definition at line 70 of file BeamFitter.h.

References freftime.

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

BeamFitter::freftime

std::time_t freftime[2]

Definition: BeamFitter.h:207

◆ getRunNumber()

int BeamFitter::getRunNumber ( )

inline

Definition at line 112 of file BeamFitter.h.

References frun.

112 { return frun; }

BeamFitter::frun

int frun

Definition: BeamFitter.h:201

◆ readEvent()

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

Definition at line 261 of file BeamFitter.cc.

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_2022v12_cff::track, reco::HitPattern::trackerLayersWithMeasurement(), 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().

                                                  {
   frun = iEvent.id().run();
   const edm::TimeValue_t ftimestamp = iEvent.time().value();
   const std::time_t ftmptime = ftimestamp >> 32;
 
   if (fbeginLumiOfFit == -1)
     freftime[0] = freftime[1] = ftmptime;
   if (freftime[0] == 0 || ftmptime < freftime[0])
     freftime[0] = ftmptime;
   if (freftime[1] == 0 || ftmptime > freftime[1])
     freftime[1] = ftmptime;
   // Update the human-readable string versions of the time
   updateBTime();
 
   flumi = iEvent.luminosityBlock();
   frunFit = frun;
   if (fbeginLumiOfFit == -1 || fbeginLumiOfFit > flumi)
     fbeginLumiOfFit = flumi;
   if (fendLumiOfFit == -1 || fendLumiOfFit < flumi)
     fendLumiOfFit = flumi;
 
   edm::Handle<reco::TrackCollection> TrackCollection;
   iEvent.getByToken(tracksToken_, TrackCollection);
 
   //------ Primary Vertices
   edm::Handle<edm::View<reco::Vertex> > PVCollection;
   bool hasPVs = false;
   edm::View<reco::Vertex> pv;
   if (iEvent.getByToken(vertexToken_, PVCollection)) {
     pv = *PVCollection;
     hasPVs = true;
   }
   //------
 
   //------ Beam Spot in current event
   edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
   const reco::BeamSpot *refBS = nullptr;
   if (iEvent.getByToken(beamSpotToken_, recoBeamSpotHandle))
     refBS = recoBeamSpotHandle.product();
   //-------
 
   const reco::TrackCollection *tracks = TrackCollection.product();
 
   double eventZ = 0;
   double averageZ = 0;
 
   for (reco::TrackCollection::const_iterator track = tracks->begin(); track != tracks->end(); ++track) {
     if (!isMuon_) {
       const reco::HitPattern &trkHP = track->hitPattern();
 
       fnPixelLayerMeas = trkHP.pixelLayersWithMeasurement();
       fnStripLayerMeas = trkHP.stripLayersWithMeasurement();
       fnTotLayerMeas = trkHP.trackerLayersWithMeasurement();
       fnPXBLayerMeas = trkHP.pixelBarrelLayersWithMeasurement();
       fnPXFLayerMeas = trkHP.pixelEndcapLayersWithMeasurement();
       fnTIBLayerMeas = trkHP.stripTIBLayersWithMeasurement();
       fnTIDLayerMeas = trkHP.stripTIDLayersWithMeasurement();
       fnTOBLayerMeas = trkHP.stripTOBLayersWithMeasurement();
       fnTECLayerMeas = trkHP.stripTECLayersWithMeasurement();
     } else {
       fnTotLayerMeas = track->numberOfValidHits();
     }
 
     fpt = track->pt();
     feta = track->eta();
     fphi0 = track->phi();
     fcharge = track->charge();
     fnormchi2 = track->normalizedChi2();
     fd0 = track->d0();
     if (refBS)
       fd0bs = -1 * track->dxy(refBS->position());
     else
       fd0bs = 0.;
 
     fsigmad0 = track->d0Error();
     fz0 = track->dz();
     fsigmaz0 = track->dzError();
     ftheta = track->theta();
     fvx = track->vx();
     fvy = track->vy();
 
     for (int i = 0; i < 5; ++i) {
       for (int j = 0; j < 5; ++j) {
         fcov[i][j] = track->covariance(i, j);
       }
     }
 
     fquality = true;
     falgo = true;
 
     if (!isMuon_) {
       if (!quality_.empty()) {
         fquality = false;
         for (unsigned int i = 0; i < quality_.size(); ++i) {
           if (debug_)
             edm::LogInfo("BeamFitter") << "quality_[" << i << "] = " << track->qualityName(quality_[i]) << std::endl;
           if (track->quality(quality_[i])) {
             fquality = true;
             break;
           }
         }
       }
 
       // Track algorithm
 
       if (!algorithm_.empty()) {
         if (std::find(algorithm_.begin(), algorithm_.end(), track->algo()) == algorithm_.end())
           falgo = false;
       }
     }
 
     // check if we have a valid PV
     fpvValid = false;
 
     if (hasPVs) {
       for (size_t ipv = 0; ipv != pv.size(); ++ipv) {
         if (!pv[ipv].isFake())
           fpvValid = true;
 
         if (ipv == 0 && !pv[0].isFake()) {
           fpvx = pv[0].x();
           fpvy = pv[0].y();
           fpvz = pv[0].z();
         }  // fix this later
       }
     }
 
     if (saveNtuple_)
       ftree_->Fill();
     ftotal_tracks++;
 
     countPass[0] = ftotal_tracks;
     // Track selection
     if (fnTotLayerMeas >= trk_MinNTotLayers_) {
       countPass[1] += 1;
       if (fnPixelLayerMeas >= trk_MinNPixLayers_) {
         countPass[2] += 1;
         if (fnormchi2 < trk_MaxNormChi2_) {
           countPass[3] += 1;
           if (falgo) {
             countPass[4] += 1;
             if (fquality) {
               countPass[5] += 1;
               if (std::abs(fd0) < trk_MaxIP_) {
                 countPass[6] += 1;
                 if (std::abs(fz0) < trk_MaxZ_) {
                   countPass[7] += 1;
                   if (fpt > trk_MinpT_) {
                     countPass[8] += 1;
                     if (std::abs(feta) < trk_MaxEta_
                         //&& fpvValid
                     ) {
                       if (debug_) {
                         edm::LogInfo("BeamFitter") << "Selected track quality = " << track->qualityMask()
                                                    << "; track algorithm = " << track->algoName()
                                                    << "= TrackAlgorithm: " << track->algo() << std::endl;
                       }
                       BSTrkParameters BSTrk(fz0, fsigmaz0, fd0, fsigmad0, fphi0, fpt, 0., 0.);
                       BSTrk.setVx(fvx);
                       BSTrk.setVy(fvy);
                       fBSvector.push_back(BSTrk);
                       averageZ += fz0;
                     }
                   }
                 }
               }
             }
           }
         }
       }
     }  // track selection
 
   }  // tracks
 
   averageZ = averageZ / (float)(fBSvector.size());
 
   for (std::vector<BSTrkParameters>::const_iterator iparam = fBSvector.begin(); iparam != fBSvector.end(); ++iparam) {
     eventZ += fabs(iparam->z0() - averageZ);
   }
 
   h1ntrks->Fill(fBSvector.size());
   h1vz_event->Fill(eventZ / (float)(fBSvector.size()));
   for (unsigned int i = 0; i < sizeof(countPass) / sizeof(countPass[0]); i++)
     h1cutFlow->SetBinContent(i + 1, countPass[i]);
 
   MyPVFitter->readEvent(iEvent);
 }

◆ resetCutFlow()

void BeamFitter::resetCutFlow ( )

inline

Definition at line 104 of file BeamFitter.h.

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

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

                       {
     h1cutFlow->Reset();
     ftotal_tracks = 0;
     for (unsigned int i = 0; i < sizeof(countPass) / sizeof(countPass[0]); i++)
       countPass[i] = 0;
   }

◆ resetLSRange()

void BeamFitter::resetLSRange ( )

inline

Definition at line 60 of file BeamFitter.h.

References fbeginLumiOfFit, and fendLumiOfFit.

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

60 { fbeginLumiOfFit = fendLumiOfFit = -1; }

BeamFitter::fendLumiOfFit

int fendLumiOfFit

Definition: BeamFitter.h:213

BeamFitter::fbeginLumiOfFit

int fbeginLumiOfFit

Definition: BeamFitter.h:212

◆ resetPVFitter()

void BeamFitter::resetPVFitter ( )

inline

Definition at line 72 of file BeamFitter.h.

References MyPVFitter, and PVFitter::resetAll().

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

72 { MyPVFitter->resetAll(); }

BeamFitter::MyPVFitter

PVFitter * MyPVFitter

Definition: BeamFitter.h:238

PVFitter::resetAll

void resetAll()

Definition: PVFitter.h:83

◆ resetRefTime()

void BeamFitter::resetRefTime ( )

inline

Definition at line 61 of file BeamFitter.h.

References freftime.

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

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

BeamFitter::freftime

std::time_t freftime[2]

Definition: BeamFitter.h:207

◆ resetTotTrk()

void BeamFitter::resetTotTrk ( )

inline

Definition at line 59 of file BeamFitter.h.

References ftotal_tracks.

59 { ftotal_tracks = 0; }

BeamFitter::ftotal_tracks

int ftotal_tracks

Definition: BeamFitter.h:159

◆ resetTrkVector()

void BeamFitter::resetTrkVector ( )

inline

Definition at line 58 of file BeamFitter.h.

References fBSvector.

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

58 { fBSvector.clear(); }

BeamFitter::fBSvector

std::vector< BSTrkParameters > fBSvector

Definition: BeamFitter.h:128

◆ resizeBSvector()

void BeamFitter::resizeBSvector ( unsigned int nsize )

inline

Definition at line 77 of file BeamFitter.h.

References fBSvector.

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

BeamFitter::fBSvector

std::vector< BSTrkParameters > fBSvector

Definition: BeamFitter.h:128

◆ resizePVvector()

void BeamFitter::resizePVvector ( unsigned int npvsize )

inline

Definition at line 80 of file BeamFitter.h.

References MyPVFitter, and PVFitter::resizepvStore().

80 { MyPVFitter->resizepvStore(npvsize); }

BeamFitter::MyPVFitter

PVFitter * MyPVFitter

Definition: BeamFitter.h:238

PVFitter::resizepvStore

void resizepvStore(unsigned int rmSize)

Definition: PVFitter.h:112

◆ runAllFitter()

void BeamFitter::runAllFitter ( )

Definition at line 783 of file BeamFitter.cc.

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().

                               {
   if (!fBSvector.empty()) {
     BSFitter *myalgo = new BSFitter(fBSvector);
     fbeamspot = myalgo->Fit_d0phi();
 
     // iterative
     if (debug_)
       std::cout << " d0-phi Iterative:" << std::endl;
     BSFitter *myitealgo = new BSFitter(fBSvector);
     myitealgo->Setd0Cut_d0phi(4.0);
     reco::BeamSpot beam_ite = myitealgo->Fit_ited0phi();
     if (debug_) {
       std::cout << beam_ite << std::endl;
       std::cout << "\n Now run tests of the different fits\n";
     }
     // from here are just tests
     std::string fit_type = "chi2";
     myalgo->SetFitVariable(std::string("z"));
     myalgo->SetFitType(std::string("chi2"));
     reco::BeamSpot beam_fit_z_chi2 = myalgo->Fit();
     if (debug_) {
       std::cout << " z Chi2 Fit ONLY:" << std::endl;
       std::cout << beam_fit_z_chi2 << std::endl;
     }
 
     fit_type = "combined";
     myalgo->SetFitVariable("z");
     myalgo->SetFitType("combined");
     reco::BeamSpot beam_fit_z_lh = myalgo->Fit();
     if (debug_) {
       std::cout << " z Log-Likelihood Fit ONLY:" << std::endl;
       std::cout << beam_fit_z_lh << std::endl;
     }
 
     myalgo->SetFitVariable("d");
     myalgo->SetFitType("d0phi");
     reco::BeamSpot beam_fit_dphi = myalgo->Fit();
     if (debug_) {
       std::cout << " d0-phi0 Fit ONLY:" << std::endl;
       std::cout << beam_fit_dphi << std::endl;
     }
     /*
     myalgo->SetFitVariable(std::string("d*z"));
     myalgo->SetFitType(std::string("likelihood"));
     reco::BeamSpot beam_fit_dz_lh = myalgo->Fit();
     if (debug_){
       std::cout << " Log-Likelihood Fit:" << std::endl;
       std::cout << beam_fit_dz_lh << std::endl;
     }
 
     myalgo->SetFitVariable(std::string("d*z"));
     myalgo->SetFitType(std::string("resolution"));
     reco::BeamSpot beam_fit_dresz_lh = myalgo->Fit();
     if(debug_){
       std::cout << " IP Resolution Fit" << std::endl;
       std::cout << beam_fit_dresz_lh << std::endl;
 
       std::cout << "c0 = " << myalgo->GetResPar0() << " +- " << myalgo->GetResPar0Err() << std::endl;
       std::cout << "c1 = " << myalgo->GetResPar1() << " +- " << myalgo->GetResPar1Err() << std::endl;
     }
     */
   } else if (debug_)
     std::cout << "No good track selected! No beam fit!" << std::endl;
 }

◆ runBeamWidthFitter()

bool BeamFitter::runBeamWidthFitter ( )

Definition at line 620 of file BeamFitter.cc.

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().

                                     {
   bool widthfit_ok = false;
   // default fit to extract beam spot info
   if (fBSvector.size() > 1) {
     edm::LogInfo("BeamFitter") << "Calculating beam spot positions('d0-phi0' method) and width using llh Fit"
                                << std::endl
                                << "We will use " << fBSvector.size() << " good tracks out of " << ftotal_tracks
                                << std::endl;
 
     BSFitter *myalgo = new BSFitter(fBSvector);
     myalgo->SetMaximumZ(trk_MaxZ_);
     myalgo->SetConvergence(convergence_);
     myalgo->SetMinimumNTrks(min_Ntrks_);
     if (inputBeamWidth_ > 0)
       myalgo->SetInputBeamWidth(inputBeamWidth_);
 
     myalgo->SetFitVariable(std::string("d*z"));
     myalgo->SetFitType(std::string("likelihood"));
     fbeamWidthFit = myalgo->Fit();
 
     //Add to .txt file
     if (writeTxt_)
       dumpBWTxtFile(outputTxt_);
 
     delete myalgo;
 
     // not fake
     if (fbeamspot.type() != 0)
       widthfit_ok = true;
   } else {
     fbeamspot.setType(reco::BeamSpot::Fake);
     edm::LogWarning("BeamFitter") << "Not enough good tracks selected! No beam fit!" << std::endl;
   }
   return widthfit_ok;
 }

◆ runFitter()

bool BeamFitter::runFitter ( )

Definition at line 609 of file BeamFitter.cc.

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

                            {
   bool fit_ok = runFitterNoTxt();
 
   if (writeTxt_)
     dumpTxtFile(outputTxt_, true);  // all reaults
   if (writeDIPTxt_ && (fit_ok || writeDIPBadFit_)) {
     dumpTxtFile(outputDIPTxt_, false);  // for DQM/DIP
   }
   return fit_ok;
 }

◆ runFitterNoTxt()

bool BeamFitter::runFitterNoTxt ( )

Definition at line 546 of file BeamFitter.cc.

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().

                                 {
   edm::LogInfo("BeamFitter") << " [BeamFitterDebugTime] freftime[0] = " << freftime[0]
                              << "; address =  " << &freftime[0] << " = " << fbeginTimeOfFit << std::endl;
   edm::LogInfo("BeamFitter") << " [BeamFitterDebugTime] freftime[1] = " << freftime[1]
                              << "; address =  " << &freftime[1] << " = " << fendTimeOfFit << std::endl;
 
   if (fbeginLumiOfFit == -1 || fendLumiOfFit == -1) {
     edm::LogWarning("BeamFitter") << "No event read! No Fitting!" << std::endl;
     return false;
   }
 
   bool fit_ok = false;
   // default fit to extract beam spot info
   if (fBSvector.size() > 1) {
     edm::LogInfo("BeamFitter") << "Calculating beam spot..." << std::endl
                                << "We will use " << fBSvector.size() << " good tracks out of " << ftotal_tracks
                                << std::endl;
 
     BSFitter *myalgo = new BSFitter(fBSvector);
     myalgo->SetMaximumZ(trk_MaxZ_);
     myalgo->SetConvergence(convergence_);
     myalgo->SetMinimumNTrks(min_Ntrks_);
     if (inputBeamWidth_ > 0)
       myalgo->SetInputBeamWidth(inputBeamWidth_);
 
     fbeamspot = myalgo->Fit();
 
     // retrieve histogram for Vz
     h1z = (TH1F *)myalgo->GetVzHisto();
 
     delete myalgo;
     if (fbeamspot.type() > 0) {  // save all results except for Fake and Unknown (all 0.)
       fit_ok = true;
       if (saveBeamFit_) {
         fx = fbeamspot.x0();
         fy = fbeamspot.y0();
         fz = fbeamspot.z0();
         fsigmaZ = fbeamspot.sigmaZ();
         fdxdz = fbeamspot.dxdz();
         fdydz = fbeamspot.dydz();
         fwidthX = fbeamspot.BeamWidthX();
         fwidthY = fbeamspot.BeamWidthY();
         fxErr = fbeamspot.x0Error();
         fyErr = fbeamspot.y0Error();
         fzErr = fbeamspot.z0Error();
         fsigmaZErr = fbeamspot.sigmaZ0Error();
         fdxdzErr = fbeamspot.dxdzError();
         fdydzErr = fbeamspot.dydzError();
         fwidthXErr = fbeamspot.BeamWidthXError();
         fwidthYErr = fbeamspot.BeamWidthYError();
         ftreeFit_->Fill();
       }
     }
   } else {  // tracks <= 1
     reco::BeamSpot tmpbs;
     fbeamspot = tmpbs;
     fbeamspot.setType(reco::BeamSpot::Fake);
     edm::LogInfo("BeamFitter") << "Not enough good tracks selected! No beam fit!" << std::endl;
   }
   fitted_ = true;
   return fit_ok;
 }

◆ runPVandTrkFitter()

bool BeamFitter::runPVandTrkFitter ( )

Definition at line 449 of file BeamFitter.cc.

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

                                    {
   // run both PV and track fitters
   bool fit_ok = false;
   bool pv_fit_ok = false;
   reco::BeamSpot bspotPV;
   reco::BeamSpot bspotTrk;
 
   // First run PV fitter
   if (MyPVFitter->IsFitPerBunchCrossing()) {
     edm::LogInfo("BeamFitter") << " [BeamFitterBxDebugTime] freftime[0] = " << freftime[0]
                                << "; address =  " << &freftime[0] << " = " << fbeginTimeOfFit << std::endl;
     edm::LogInfo("BeamFitter") << " [BeamFitterBxDebugTime] freftime[1] = " << freftime[1]
                                << "; address =  " << &freftime[1] << " = " << fendTimeOfFit << std::endl;
 
     if (MyPVFitter->runBXFitter()) {
       fbspotPVMap = MyPVFitter->getBeamSpotMap();
       pv_fit_ok = true;
     }
     if (writeTxt_)
       dumpTxtFile(outputTxt_, true);  // all reaults
     if (writeDIPTxt_ && (pv_fit_ok || writeDIPBadFit_)) {
       dumpTxtFile(outputDIPTxt_, false);  // for DQM/DIP
     }
     return pv_fit_ok;
   }
 
   if (MyPVFitter->runFitter()) {
     bspotPV = MyPVFitter->getBeamSpot();
 
     // take beam width from PV fit and pass it to track fitter
     // assume circular transverse beam width
     inputBeamWidth_ = sqrt(pow(bspotPV.BeamWidthX(), 2) + pow(bspotPV.BeamWidthY(), 2)) / sqrt(2);
     pv_fit_ok = true;
 
   } else {
     // problems with PV fit
     bspotPV.setType(reco::BeamSpot::Unknown);
     bspotTrk.setType(reco::BeamSpot::Unknown);  //propagate error to trk beam spot
   }
 
   if (runFitterNoTxt()) {
     bspotTrk = fbeamspot;
     fit_ok = true;
   } else {
     // beamfit failed, flagged as empty beam spot
     bspotTrk.setType(reco::BeamSpot::Fake);
     fit_ok = false;
   }
 
   // combined results into one single beam spot
 
   // to pass errors I have to create another beam spot object
   reco::BeamSpot::CovarianceMatrix matrix;
   for (int j = 0; j < 7; ++j) {
     for (int k = j; k < 7; ++k) {
       matrix(j, k) = bspotTrk.covariance(j, k);
     }
   }
   // change beam width error to one from PV
   if (pv_fit_ok && fit_ok) {
     matrix(6, 6) = MyPVFitter->getWidthXerr() * MyPVFitter->getWidthXerr();
 
     // get Z and sigmaZ from PV fit
     matrix(2, 2) = bspotPV.covariance(2, 2);
     matrix(3, 3) = bspotPV.covariance(3, 3);
     reco::BeamSpot tmpbs(reco::BeamSpot::Point(bspotTrk.x0(), bspotTrk.y0(), bspotPV.z0()),
                          bspotPV.sigmaZ(),
                          bspotTrk.dxdz(),
                          bspotTrk.dydz(),
                          bspotPV.BeamWidthX(),
                          matrix,
                          bspotPV.type());
     tmpbs.setBeamWidthY(bspotPV.BeamWidthY());
     // overwrite beam spot result
     fbeamspot = tmpbs;
   }
   if (pv_fit_ok && fit_ok) {
     fbeamspot.setType(bspotPV.type());
   } else if (!pv_fit_ok && fit_ok) {
     fbeamspot.setType(reco::BeamSpot::Unknown);
   } else if (pv_fit_ok && !fit_ok) {
     fbeamspot.setType(reco::BeamSpot::Unknown);
   } else if (!pv_fit_ok && !fit_ok) {
     fbeamspot.setType(reco::BeamSpot::Fake);
   }
 
   if (writeTxt_)
     dumpTxtFile(outputTxt_, true);  // all reaults
   if (writeDIPTxt_ && ((fit_ok && pv_fit_ok) || writeDIPBadFit_)) {
     dumpTxtFile(outputDIPTxt_, false);  // for DQM/DIP
     for (size_t i = 0; i < 7; i++)
       ForDIPPV_.push_back(0.0);
   }
 
   return fit_ok && pv_fit_ok;
 }

◆ setFitLSRange()

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

inline

Definition at line 115 of file BeamFitter.h.

References fbeginLumiOfFit, and fendLumiOfFit.

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

                                        {
     fbeginLumiOfFit = ls0;
     fendLumiOfFit = ls1;
   }

◆ SetPVInfo()

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

inline

Definition at line 83 of file BeamFitter.h.

References ForDIPPV_.

                                               {
     ForDIPPV_.clear();
     ForDIPPV_.assign(v1_.begin(), v1_.end());
   }

◆ setRefTime()

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

inline

Definition at line 62 of file BeamFitter.h.

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

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

                                         {
     freftime[0] = t0;
     freftime[1] = t1;
     // Make sure the string representation of the time
     // is up-to-date
     updateBTime();
   }

◆ setRun()

void BeamFitter::setRun ( int run )

inline

Definition at line 119 of file BeamFitter.h.

References frun, and writedatasetfile::run.

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

119 { frun = run; }

BeamFitter::frun

int frun

Definition: BeamFitter.h:201

writedatasetfile.run

run

Definition: writedatasetfile.py:27

◆ subtractFromCutFlow()

void BeamFitter::subtractFromCutFlow ( const TH1F * toSubtract )

inline

Definition at line 97 of file BeamFitter.h.

References countPass, h1cutFlow, and mps_fire::i.

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

◆ 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.

References fbeginTimeOfFit, fendTimeOfFit, and freftime.

Referenced by readEvent(), and setRefTime().

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

◆ write2DB()

void BeamFitter::write2DB ( )

Definition at line 746 of file BeamFitter.cc.

References cond::service::PoolDBOutputService::appendOneIOV(), reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), cond::service::PoolDBOutputService::beginOfTime(), gather_cfg::cout, reco::BeamSpot::covariance(), cond::service::PoolDBOutputService::createOneIOV(), cond::service::PoolDBOutputService::currentTime(), reco::BeamSpot::dxdz(), reco::BeamSpot::dydz(), 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().

                           {
   BeamSpotObjects pBSObjects;
 
   pBSObjects.setPosition(fbeamspot.position().X(), fbeamspot.position().Y(), fbeamspot.position().Z());
   //std::cout << " wrote: x= " << fbeamspot.position().X() << " y= "<< fbeamspot.position().Y() << " z= " << fbeamspot.position().Z() << std::endl;
   pBSObjects.setSigmaZ(fbeamspot.sigmaZ());
   pBSObjects.setdxdz(fbeamspot.dxdz());
   pBSObjects.setdydz(fbeamspot.dydz());
   //if (inputBeamWidth_ > 0 ) {
   //  std::cout << " beam width value forced to be " << inputBeamWidth_ << std::endl;
   //  pBSObjects->SetBeamWidthX(inputBeamWidth_);
   //  pBSObjects->SetBeamWidthY(inputBeamWidth_);
   //} else {
   // need to fix this
   //std::cout << " using default value, 15e-4, for beam width!!!"<<std::endl;
   pBSObjects.setBeamWidthX(fbeamspot.BeamWidthX());
   pBSObjects.setBeamWidthY(fbeamspot.BeamWidthY());
   //}
 
   for (int i = 0; i < 7; ++i) {
     for (int j = 0; j < 7; ++j) {
       pBSObjects.setCovariance(i, j, fbeamspot.covariance(i, j));
     }
   }
   edm::Service<cond::service::PoolDBOutputService> poolDbService;
   if (poolDbService.isAvailable()) {
     std::cout << "poolDBService available" << std::endl;
     if (poolDbService->isNewTagRequest("BeamSpotObjectsRcd")) {
       std::cout << "new tag requested" << std::endl;
       poolDbService->createOneIOV<BeamSpotObjects>(pBSObjects, poolDbService->beginOfTime(), "BeamSpotObjectsRcd");
     } else {
       std::cout << "no new tag requested" << std::endl;
       poolDbService->appendOneIOV<BeamSpotObjects>(pBSObjects, poolDbService->currentTime(), "BeamSpotObjectsRcd");
     }
   }
 }

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().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ BeamFitter() [1/2]

◆ BeamFitter() [2/2]

◆ ~BeamFitter()

Member Function Documentation

◆ dumpBWTxtFile()

◆ dumpTxtFile()

◆ fillDescription()

◆ getBeamSpot()

◆ getBeamSpotMap()

◆ getBeamWidth()

◆ getBSvector()

◆ getCutFlow()

◆ getFitLSRange()

◆ getNPVs()

◆ getNPVsperBX()

◆ getNTracks()

◆ getPVvectorSize()

◆ getRefTime()

◆ getRunNumber()

◆ readEvent()

◆ resetCutFlow()

◆ resetLSRange()

◆ resetPVFitter()

◆ resetRefTime()

◆ resetTotTrk()

◆ resetTrkVector()

◆ resizeBSvector()

◆ resizePVvector()

◆ runAllFitter()

◆ runBeamWidthFitter()

◆ runFitter()

◆ runFitterNoTxt()

◆ runPVandTrkFitter()

◆ setFitLSRange()

◆ SetPVInfo()

◆ setRefTime()

◆ setRun()

◆ subtractFromCutFlow()

◆ updateBTime()

◆ write2DB()

Member Data Documentation

◆ algorithm_

◆ appendRunTxt_

◆ beamSpotToken_

◆ convergence_

◆ countPass

◆ debug_

◆ falgo

◆ fasciiDIP

◆ fasciiFile

◆ fbeamspot

◆ fbeamWidthFit

◆ fbeginLumiOfFit

◆ fbeginTimeOfFit

◆ fbspotPVMap

◆ fBSvector

◆ fcharge

◆ fcov

◆ fd0

◆ fd0bs

◆ fdxdz

◆ fdxdzErr

◆ fdydz

◆ fdydzErr

◆ fendLumiOfFit

◆ fendTimeOfFit

◆ feta

◆ ffilename_changed

◆ file_

◆ fitted_

◆ flumi

◆ fnormchi2

◆ fnPixelLayerMeas

◆ fnPXBLayerMeas