AlcaBeamSpotManager Class Reference

#include <AlcaBeamSpotManager.h>

Public Types
typedef std::map< edm::LuminosityBlockNumber_t, std::pair< edm::Timestamp, reco::BeamSpot > >::iterator	bsMap_iterator

Public Member Functions
	AlcaBeamSpotManager (void)
	AlcaBeamSpotManager (const edm::ParameterSet &, edm::ConsumesCollector &&)
void	createWeightedPayloads (void)
const std::map< edm::LuminosityBlockNumber_t, std::pair< edm::Timestamp, reco::BeamSpot > > &	getPayloads (void)
void	readLumi (const edm::LuminosityBlock &)
void	reset (void)
virtual	~AlcaBeamSpotManager (void)

Private Member Functions
std::pair< float, float >	delta (const float &x, const float &xError, const float &nextX, const float &nextXError)
float	deltaSig (const float &num, const float &den)
reco::BeamSpot	weight (const bsMap_iterator &begin, const bsMap_iterator &end)
void	weight (double &mean, double &meanError, const double &val, const double &valError)

Private Attributes
std::string	beamSpotLabel_
std::map< edm::LuminosityBlockNumber_t, std::pair< edm::Timestamp, reco::BeamSpot > >	beamSpotMap_
std::string	beamSpotModuleName_
std::string	beamSpotOutputBase_
edm::InputTag	beamSpotTag_
edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken_
double	sigmaZCut_

Detailed Description

No description available.

Author

L. Uplegger F. Yumiceva - Fermilab

Definition at line 18 of file AlcaBeamSpotManager.h.

Member Typedef Documentation

bsMap_iterator

typedef std::map<edm::LuminosityBlockNumber_t, std::pair<edm::Timestamp, reco::BeamSpot> >::iterator AlcaBeamSpotManager::bsMap_iterator

Definition at line 31 of file AlcaBeamSpotManager.h.

Constructor & Destructor Documentation

AlcaBeamSpotManager() [1/2]

AlcaBeamSpotManager::AlcaBeamSpotManager

(

void

)

Definition at line 19 of file AlcaBeamSpotManager.cc.

{}

AlcaBeamSpotManager() [2/2]

AlcaBeamSpotManager::AlcaBeamSpotManager	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &&	iC
	)

Definition at line 22 of file AlcaBeamSpotManager.cc.

References align::BeamSpot, beamSpotLabel_, beamSpotModuleName_, beamSpotOutputBase_, beamSpotTag_, beamSpotToken_, edm::InLumi, and reset().

    : beamSpotOutputBase_(iConfig.getParameter<ParameterSet>("AlcaBeamSpotHarvesterParameters")
                              .getUntrackedParameter<std::string>("BeamSpotOutputBase")),
      beamSpotModuleName_(iConfig.getParameter<ParameterSet>("AlcaBeamSpotHarvesterParameters")
                              .getUntrackedParameter<std::string>("BeamSpotModuleName")),
      beamSpotLabel_(iConfig.getParameter<ParameterSet>("AlcaBeamSpotHarvesterParameters")
                         .getUntrackedParameter<std::string>("BeamSpotLabel")),
      sigmaZCut_(iConfig.getParameter<ParameterSet>("AlcaBeamSpotHarvesterParameters")
                     .getUntrackedParameter<double>("SigmaZCut")) {
  edm::InputTag beamSpotTag_(beamSpotModuleName_, beamSpotLabel_);
  beamSpotToken_ = iC.consumes<reco::BeamSpot, edm::InLumi>(beamSpotTag_);
  LogInfo("AlcaBeamSpotManager") << "Output base: " << beamSpotOutputBase_ << std::endl;
  reset();
}

~AlcaBeamSpotManager()

AlcaBeamSpotManager::~AlcaBeamSpotManager ( void  )

virtual

Definition at line 38 of file AlcaBeamSpotManager.cc.

{}

Member Function Documentation

createWeightedPayloads()

void AlcaBeamSpotManager::createWeightedPayloads

(

void

)

Definition at line 61 of file AlcaBeamSpotManager.cc.

References beamSpotMap_, beamSpotOutputBase_, reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthXError(), reco::BeamSpot::BeamWidthY(), reco::BeamSpot::BeamWidthYError(), delta(), deltaSig(), l1tTrackerHTMiss_cfi::deltaZ, reco::BeamSpot::dxdz(), reco::BeamSpot::dxdzError(), reco::BeamSpot::dydz(), reco::BeamSpot::dydzError(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, align_cfg::iteration, lut2db_cfg::limit, reco::BeamSpot::sigmaZ(), reco::BeamSpot::sigmaZ0Error(), sigmaZCut_, ntuplemaker::tmprun, align::Tracker, funct::true, reco::Unknown, weight(), reco::BeamSpot::x0(), reco::BeamSpot::x0Error(), reco::BeamSpot::y0(), reco::BeamSpot::y0Error(), reco::BeamSpot::z0(), and reco::BeamSpot::z0Error().

Referenced by AlcaBeamSpotHarvester::endRun().

                                                     {
  vector<bsMap_iterator> listToErase;
  for (bsMap_iterator it = beamSpotMap_.begin(); it != beamSpotMap_.end(); it++) {
    if (it->second.second.type() != BeamSpot::Tracker || it->second.second.sigmaZ() < sigmaZCut_) {
      listToErase.push_back(it);
    }
  }
  for (vector<bsMap_iterator>::iterator it = listToErase.begin(); it != listToErase.end(); it++) {
    beamSpotMap_.erase(*it);
  }
  if (beamSpotMap_.size() <= 1) {
    return;
  }
  // Return only if lumibased since the collapsing alghorithm requires the next
  // and next to next lumi sections
  else if (beamSpotMap_.size() == 2 && beamSpotOutputBase_ == "lumibased") {
    return;
  }
  if (beamSpotOutputBase_ == "lumibased") {
    //    bsMap_iterator referenceBS = beamSpotMap_.begin();
    bsMap_iterator firstBS = beamSpotMap_.begin();
    //    bsMap_iterator lastBS      = beamSpotMap_.begin();
    bsMap_iterator currentBS = beamSpotMap_.begin();
    bsMap_iterator nextBS = ++beamSpotMap_.begin();
    bsMap_iterator nextNextBS = ++(++(beamSpotMap_.begin()));

    reco::BeamSpot currentBSObj;
    reco::BeamSpot nextBSObj;

    map<LuminosityBlockNumber_t, std::pair<edm::Timestamp, BeamSpot>> tmpBeamSpotMap_;
    bool docreate = true;
    bool endOfRun = false;  // Added
    bool docheck = true;    // Added
    bool foundShift = false;
    long countlumi = 0;   // Added
    string tmprun = "";   // Added
    long maxNlumis = 20;  // Added
    //    if weighted:
    //        maxNlumis = 999999999

    unsigned int iteration = 0;
    // while(nextNextBS!=beamSpotMap_.end()){
    while (nextBS != beamSpotMap_.end()) {
      LogInfo("AlcaBeamSpotManager") << "Iteration: " << iteration << " size: " << beamSpotMap_.size() << "\n"
                                     << "Lumi: " << currentBS->first << "\n"
                                     << currentBS->second.second << "\n"
                                     << nextBS->first << "\n"
                                     << nextBS->second.second << endl;
      if (nextNextBS != beamSpotMap_.end())
        LogInfo("AlcaBeamSpotManager") << nextNextBS->first << "\n" << nextNextBS->second.second << endl;

      if (docreate) {
        firstBS = currentBS;
        docreate = false;  // Added
      }
      // if(iteration >= beamSpotMap_.size()-3){
      if (iteration >= beamSpotMap_.size() - 2) {
        LogInfo("AlcaBeamSpotManager") << "Reached lumi " << currentBS->first
                                       << " now close payload because end of data has been reached.";
        docreate = true;
        endOfRun = true;
      }
      // check we run over the same run
      //      if (ibeam->first.Run() != inextbeam->first.Run()){
      //        LogInfo("AlcaBeamSpotManager")
      //          << "close payload because end of run.";
      //        docreate = true;
      //      }
      // check maximum lumi counts
      if (countlumi == maxNlumis - 1) {
        LogInfo("AlcaBeamSpotManager") << "close payload because maximum lumi "
                                          "sections accumulated within run ";
        docreate = true;
        countlumi = 0;
      }
      //      if weighted:
      //          docheck = False
      // check offsets
      if (docheck) {
        foundShift = false;
        LogInfo("AlcaBeamSpotManager") << "Checking checking!" << endl;
        float limit = 0;
        pair<float, float> adelta1;
        pair<float, float> adelta2;
        pair<float, float> adelta;
        pair<float, float> adelta1dxdz;
        pair<float, float> adelta2dxdz;
        pair<float, float> adelta1dydz;
        pair<float, float> adelta2dydz;
        pair<float, float> adelta1widthX;
        pair<float, float> adelta2widthX;
        pair<float, float> adelta1widthY;
        pair<float, float> adelta2widthY;
        pair<float, float> adelta1z0;
        pair<float, float> adelta1sigmaZ;

        // define minimum limit
        float min_limit = 0.0025;

        // limit for x and y
        limit = currentBS->second.second.BeamWidthX() / 2.;
        if (limit < min_limit) {
          limit = min_limit;
        }

        currentBSObj = currentBS->second.second;
        nextBSObj = nextBS->second.second;

        // check movements in X
        adelta1 = delta(currentBSObj.x0(), currentBSObj.x0Error(), nextBSObj.x0(), nextBSObj.x0Error());
        adelta2 = pair<float, float>(0., 1.e9);
        if (nextNextBS->second.second.type() != -1) {
          adelta2 = delta(nextBS->second.second.x0(),
                          nextBSObj.x0Error(),
                          nextNextBS->second.second.x0(),
                          nextNextBS->second.second.x0Error());
        }
        bool deltaX = (deltaSig(adelta1.first, adelta1.second) > 3.5 && adelta1.first >= limit) ? true : false;
        if (iteration < beamSpotMap_.size() - 2) {
          if (!deltaX && adelta1.first * adelta2.first > 0. && fabs(adelta1.first + adelta2.first) >= limit) {
            LogInfo("AlcaBeamSpotManager")
                << " positive, " << (adelta1.first + adelta2.first) << " limit=" << limit << endl;
            deltaX = true;
          } else if (deltaX && adelta1.first * adelta2.first < 0 && adelta2.first != 0 &&
                     fabs(adelta1.first / adelta2.first) > 0.33 && fabs(adelta1.first / adelta2.first) < 3) {
            LogInfo("AlcaBeamSpotManager") << " negative, " << adelta1.first / adelta2.first << endl;
            deltaX = false;
          }
        }

        // calculating all deltas
        adelta1dxdz = delta(currentBSObj.dxdz(), currentBSObj.dxdzError(), nextBSObj.dxdz(), nextBSObj.dxdzError());
        adelta2dxdz = pair<float, float>(0., 1.e9);
        adelta1dydz = delta(currentBSObj.dydz(), currentBSObj.dydzError(), nextBSObj.dydz(), nextBSObj.dydzError());
        adelta2dydz = pair<float, float>(0., 1.e9);
        adelta1widthX = delta(currentBSObj.BeamWidthX(),
                              currentBSObj.BeamWidthXError(),
                              nextBSObj.BeamWidthX(),
                              nextBSObj.BeamWidthXError());
        adelta2widthX = pair<float, float>(0., 1.e9);
        adelta1widthY = delta(currentBSObj.BeamWidthY(),
                              currentBSObj.BeamWidthYError(),
                              nextBSObj.BeamWidthY(),
                              nextBSObj.BeamWidthYError());
        adelta2widthY = pair<float, float>(0., 1.e9);
        adelta1z0 = delta(currentBSObj.z0(), currentBSObj.z0Error(), nextBSObj.z0(), nextBSObj.z0Error());
        adelta1sigmaZ =
            delta(currentBSObj.sigmaZ(), currentBSObj.sigmaZ0Error(), nextBSObj.sigmaZ(), nextBSObj.sigmaZ0Error());

        // check movements in Y
        adelta1 = delta(currentBSObj.y0(), currentBSObj.y0Error(), nextBSObj.y0(), nextBSObj.y0Error());
        adelta2 = pair<float, float>(0., 1.e9);
        if (nextNextBS->second.second.type() != BeamSpot::Unknown) {
          adelta2 = delta(
              nextBSObj.y0(), nextBSObj.y0Error(), nextNextBS->second.second.y0(), nextNextBS->second.second.y0Error());
          adelta2dxdz = delta(nextBSObj.dxdz(),
                              nextBSObj.dxdzError(),
                              nextNextBS->second.second.dxdz(),
                              nextNextBS->second.second.dxdzError());
          adelta2dydz = delta(nextBSObj.dydz(),
                              nextBSObj.dydzError(),
                              nextNextBS->second.second.dydz(),
                              nextNextBS->second.second.dydzError());
          adelta2widthX = delta(nextBSObj.BeamWidthX(),
                                nextBSObj.BeamWidthXError(),
                                nextNextBS->second.second.BeamWidthX(),
                                nextNextBS->second.second.BeamWidthXError());
          adelta2widthY = delta(nextBSObj.BeamWidthY(),
                                nextBSObj.BeamWidthYError(),
                                nextNextBS->second.second.BeamWidthY(),
                                nextNextBS->second.second.BeamWidthYError());
        }
        bool deltaY = (deltaSig(adelta1.first, adelta1.second) > 3.5 && adelta1.first >= limit) ? true : false;
        if (iteration < beamSpotMap_.size() - 2) {
          if (!deltaY && adelta1.first * adelta2.first > 0. && fabs(adelta1.first + adelta2.first) >= limit) {
            LogInfo("AlcaBeamSpotManager")
                << " positive, " << (adelta1.first + adelta2.first) << " limit=" << limit << endl;
            deltaY = true;
          } else if (deltaY && adelta1.first * adelta2.first < 0 && adelta2.first != 0 &&
                     fabs(adelta1.first / adelta2.first) > 0.33 && fabs(adelta1.first / adelta2.first) < 3) {
            LogInfo("AlcaBeamSpotManager") << " negative, " << adelta1.first / adelta2.first << endl;
            deltaY = false;
          }
        }

        limit = currentBSObj.sigmaZ() / 2.;
        bool deltaZ =
            (deltaSig(adelta1z0.first, adelta1z0.second) > 3.5 && fabs(adelta1z0.first) >= limit) ? true : false;
        adelta =
            delta(currentBSObj.sigmaZ(), currentBSObj.sigmaZ0Error(), nextBSObj.sigmaZ(), nextBSObj.sigmaZ0Error());
        bool deltasigmaZ = (deltaSig(adelta.first, adelta.second) > 5.0) ? true : false;
        bool deltadxdz = false;
        bool deltadydz = false;
        bool deltawidthX = false;
        bool deltawidthY = false;

        if (iteration < beamSpotMap_.size() - 2) {
          adelta = delta(currentBSObj.dxdz(), currentBSObj.dxdzError(), nextBSObj.dxdz(), nextBSObj.dxdzError());
          deltadxdz = (deltaSig(adelta.first, adelta.second) > 5.0) ? true : false;
          if (deltadxdz && (adelta1dxdz.first * adelta2dxdz.first) < 0 && adelta2dxdz.first != 0 &&
              fabs(adelta1dxdz.first / adelta2dxdz.first) > 0.33 && fabs(adelta1dxdz.first / adelta2dxdz.first) < 3) {
            deltadxdz = false;
          }

          adelta = delta(currentBSObj.dydz(), currentBSObj.dydzError(), nextBSObj.dydz(), nextBSObj.dydzError());
          deltadydz = (deltaSig(adelta.first, adelta.second) > 5.0) ? true : false;
          if (deltadydz && (adelta1dydz.first * adelta2dydz.first) < 0 && adelta2dydz.first != 0 &&
              fabs(adelta1dydz.first / adelta2dydz.first) > 0.33 && fabs(adelta1dydz.first / adelta2dydz.first) < 3) {
            deltadydz = false;
          }

          adelta = delta(currentBSObj.BeamWidthX(),
                         currentBSObj.BeamWidthXError(),
                         nextBSObj.BeamWidthX(),
                         nextBSObj.BeamWidthXError());
          deltawidthX = (deltaSig(adelta.first, adelta.second) > 5.0) ? true : false;
          if (deltawidthX && (adelta1widthX.first * adelta2widthX.first) < 0 && adelta2widthX.first != 0 &&
              fabs(adelta1widthX.first / adelta2widthX.first) > 0.33 &&
              fabs(adelta1widthX.first / adelta2widthX.first) < 3) {
            deltawidthX = false;
          }

          adelta = delta(currentBSObj.BeamWidthY(),
                         currentBSObj.BeamWidthYError(),
                         nextBSObj.BeamWidthY(),
                         nextBSObj.BeamWidthYError());
          deltawidthY = (deltaSig(adelta.first, adelta.second) > 5.0) ? true : false;
          if (deltawidthY && (adelta1widthY.first * adelta2widthY.first) < 0 && adelta2widthY.first != 0 &&
              fabs(adelta1widthY.first / adelta2widthY.first) > 0.33 &&
              fabs(adelta1widthY.first / adelta2widthY.first) < 3) {
            deltawidthY = false;
          }
        }
        if (deltaX || deltaY || deltaZ || deltasigmaZ || deltadxdz || deltadydz || deltawidthX || deltawidthY) {
          docreate = true;
          foundShift = true;
          LogInfo("AlcaBeamSpotManager") << "close payload because of movement in"
                                         << " X=" << deltaX << ", Y=" << deltaY << ", Z=" << deltaZ
                                         << ", sigmaZ=" << deltasigmaZ << ", dxdz=" << deltadxdz
                                         << ", dydz=" << deltadydz << ", widthX=" << deltawidthX
                                         << ", widthY=" << deltawidthY << endl;
        }
      }
      if (docreate) {
        std::pair<edm::Timestamp, reco::BeamSpot> thepair;
        if (foundShift) {
          thepair = std::make_pair(currentBS->second.first, weight(firstBS, nextBS));
          tmpBeamSpotMap_[firstBS->first] = thepair;
          if (endOfRun) {
            // if we're here, then we need to found a shift in the last LS
            // We already created a new IOV, now create one just for the last LS
            thepair = std::make_pair(nextBS->second.first, nextBSObj);
            tmpBeamSpotMap_[nextBS->first] = thepair;
          }
        } else if (!foundShift && !endOfRun) {  // maxLS reached
          thepair = std::make_pair(currentBS->second.first, weight(firstBS, nextBS));
          tmpBeamSpotMap_[firstBS->first] = thepair;
        } else {  // end of run with no shift detectred in last LS
          thepair = std::make_pair(currentBS->second.first, weight(firstBS, beamSpotMap_.end()));
          tmpBeamSpotMap_[firstBS->first] = thepair;
        }
        firstBS = nextBS;
        countlumi = 0;
      }
      // tmprun = currentBS->second.Run
      // increase the counter by one only if the IOV hasn't been closed
      if (!docreate)
        ++countlumi;

      currentBS = nextBS;
      nextBS = nextNextBS;
      nextNextBS++;
      ++iteration;
    }
    beamSpotMap_.clear();
    beamSpotMap_ = tmpBeamSpotMap_;
  } else if (beamSpotOutputBase_ == "runbased") {
    LuminosityBlockNumber_t firstLumi = beamSpotMap_.begin()->first;
    std::pair<edm::Timestamp, reco::BeamSpot> thepair(beamSpotMap_.begin()->second.first,
                                                      weight(beamSpotMap_.begin(), beamSpotMap_.end()));
    beamSpotMap_.clear();
    beamSpotMap_[firstLumi] = thepair;
  } else {
    LogInfo("AlcaBeamSpotManager") << "Unrecognized BeamSpotOutputBase parameter: " << beamSpotOutputBase_ << endl;
  }
}

delta()

pair< float, float > AlcaBeamSpotManager::delta ( const float &  x, const float &  xError, const float &  nextX, const float &  nextXError  )

private

Definition at line 404 of file AlcaBeamSpotManager.cc.

References funct::pow(), mathSSE::sqrt(), x, and run3scouting_cff::xError.

Referenced by createWeightedPayloads().

                                                                       {
  return pair<float, float>(x - nextX, sqrt(pow(xError, 2) + pow(nextXError, 2)));
}

deltaSig()

float AlcaBeamSpotManager::deltaSig ( const float &  num, const float &  den  )

private

Definition at line 412 of file AlcaBeamSpotManager.cc.

References ALCARECOEcalPhiSym_cff::float, and EgammaValidation_cff::num.

Referenced by createWeightedPayloads().

                                                                      {
  if (den != 0) {
    return fabs(num / den);
  } else {
    return float(LONG_MAX);
  }
}

getPayloads()

const std::map<edm::LuminosityBlockNumber_t, std::pair<edm::Timestamp, reco::BeamSpot> >& AlcaBeamSpotManager::getPayloads ( void  )

inline

Definition at line 27 of file AlcaBeamSpotManager.h.

References beamSpotMap_.

Referenced by AlcaBeamSpotHarvester::endRun().

                                                                                                   {
    return beamSpotMap_;
  }

readLumi()

void AlcaBeamSpotManager::readLumi

(

const edm::LuminosityBlock &

iLumi

)

Definition at line 43 of file AlcaBeamSpotManager.cc.

References beamSpotMap_, beamSpotToken_, Hcal_Conditions_forGlobalTag_cff::iLumi, edm::HandleBase::isValid(), and edm::Handle< T >::product().

Referenced by AlcaBeamSpotHarvester::endLuminosityBlock().

                                                               {
  Handle<BeamSpot> beamSpotHandle;
  iLumi.getByToken(beamSpotToken_, beamSpotHandle);

  if (beamSpotHandle.isValid()) {  // check the product
    std::pair<edm::Timestamp, reco::BeamSpot> time_bs(iLumi.beginTime(), *beamSpotHandle);
    beamSpotMap_[iLumi.luminosityBlock()] = time_bs;
    const BeamSpot *aBeamSpot = &beamSpotMap_[iLumi.luminosityBlock()].second;
    aBeamSpot = beamSpotHandle.product();
    LogInfo("AlcaBeamSpotManager") << "Lumi: " << iLumi.luminosityBlock() << std::endl;
    LogInfo("AlcaBeamSpotManager") << *aBeamSpot << std::endl;
  } else {
    LogInfo("AlcaBeamSpotManager") << "Lumi: " << iLumi.luminosityBlock() << std::endl;
    LogInfo("AlcaBeamSpotManager") << "   BS is not valid!" << std::endl;
  }
}

reset()

void AlcaBeamSpotManager::reset

(

void

)

Definition at line 41 of file AlcaBeamSpotManager.cc.

References beamSpotMap_.

Referenced by MatrixReader.MatrixReader::__init__(), data_sources.json_list::__next__(), AlcaBeamSpotManager(), AlcaBeamSpotHarvester::beginRun(), and MatrixReader.MatrixReader::showRaw().

{ beamSpotMap_.clear(); }

weight() [1/2]

BeamSpot AlcaBeamSpotManager::weight ( const bsMap_iterator &  begin, const bsMap_iterator &  end  )

private

Definition at line 349 of file AlcaBeamSpotManager.cc.

References SplitLinear::begin, dxdz, dydz, mps_fire::end, relativeConstraints::error, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, reco::BeamSpot::setBeamWidthY(), beamSpotPI::sigmaZ, globals_cff::sigmaZError, align::Tracker, reco::Unknown, ApeEstimator_cff::widthX, ApeEstimator_cff::widthY, x, run3scouting_cff::xError, y, run3scouting_cff::yError, z, and globals_cff::zError.

Referenced by createWeightedPayloads().

                                                                                           {
  double x = 0, xError = 0;
  double y = 0, yError = 0;
  double z = 0, zError = 0;
  double sigmaZ = 0, sigmaZError = 0;
  double dxdz = 0, dxdzError = 0;
  double dydz = 0, dydzError = 0;
  double widthX = 0, widthXError = 0;
  double widthY = 0, widthYError = 0;
  LogInfo("AlcaBeamSpotManager") << "Weighted BeamSpot will span lumi " << begin->first << " to " << end->first << endl;

  BeamSpot::BeamType type = BeamSpot::Unknown;
  for (bsMap_iterator it = begin; it != end; it++) {
    weight(x, xError, it->second.second.x0(), it->second.second.x0Error());
    weight(y, yError, it->second.second.y0(), it->second.second.y0Error());
    weight(z, zError, it->second.second.z0(), it->second.second.z0Error());
    weight(sigmaZ, sigmaZError, it->second.second.sigmaZ(), it->second.second.sigmaZ0Error());
    weight(dxdz, dxdzError, it->second.second.dxdz(), it->second.second.dxdzError());
    weight(dydz, dydzError, it->second.second.dydz(), it->second.second.dydzError());
    weight(widthX, widthXError, it->second.second.BeamWidthX(), it->second.second.BeamWidthXError());
    weight(widthY, widthYError, it->second.second.BeamWidthY(), it->second.second.BeamWidthYError());
    if (it->second.second.type() == BeamSpot::Tracker) {
      type = BeamSpot::Tracker;
    }
  }
  BeamSpot::Point bsPosition(x, y, z);
  BeamSpot::CovarianceMatrix error;
  error(0, 0) = xError * xError;
  error(1, 1) = yError * yError;
  error(2, 2) = zError * zError;
  error(3, 3) = sigmaZError * sigmaZError;
  error(4, 4) = dxdzError * dxdzError;
  error(5, 5) = dydzError * dydzError;
  error(6, 6) = widthXError * widthXError;
  BeamSpot weightedBeamSpot(bsPosition, sigmaZ, dxdz, dydz, widthX, error, type);
  weightedBeamSpot.setBeamWidthY(widthY);
  LogInfo("AlcaBeamSpotManager") << "Weighted BeamSpot will be:" << '\n' << weightedBeamSpot << endl;
  return weightedBeamSpot;
}

weight() [2/2]

void AlcaBeamSpotManager::weight ( double &  mean, double &  meanError, const double &  val, const double &  valError  )

private

Definition at line 390 of file AlcaBeamSpotManager.cc.

References MillePedeFileConverter_cfg::e, SiStripPI::mean, mathSSE::sqrt(), and heppy_batch::val.

                                                                                                           {
  double tmpError = 0;
  if (meanError < 1e-8) {
    tmpError = 1 / (valError * valError);
    mean = val * tmpError;
  } else {
    tmpError = 1 / (meanError * meanError) + 1 / (valError * valError);
    mean = mean / (meanError * meanError) + val / (valError * valError);
  }
  mean = mean / tmpError;
  meanError = sqrt(1 / tmpError);
}

Member Data Documentation

beamSpotLabel_

std::string AlcaBeamSpotManager::beamSpotLabel_

private

Definition at line 42 of file AlcaBeamSpotManager.h.

Referenced by AlcaBeamSpotManager().

beamSpotMap_

std::map<edm::LuminosityBlockNumber_t, std::pair<edm::Timestamp, reco::BeamSpot> > AlcaBeamSpotManager::beamSpotMap_

private

Definition at line 38 of file AlcaBeamSpotManager.h.

Referenced by createWeightedPayloads(), getPayloads(), readLumi(), and reset().

beamSpotModuleName_

std::string AlcaBeamSpotManager::beamSpotModuleName_

private

Definition at line 41 of file AlcaBeamSpotManager.h.

Referenced by AlcaBeamSpotManager().

beamSpotOutputBase_

std::string AlcaBeamSpotManager::beamSpotOutputBase_

private

Definition at line 40 of file AlcaBeamSpotManager.h.

Referenced by AlcaBeamSpotManager(), and createWeightedPayloads().

beamSpotTag_

edm::InputTag AlcaBeamSpotManager::beamSpotTag_

private

Definition at line 44 of file AlcaBeamSpotManager.h.

Referenced by AlcaBeamSpotManager().

beamSpotToken_

edm::EDGetTokenT<reco::BeamSpot> AlcaBeamSpotManager::beamSpotToken_

private

Definition at line 45 of file AlcaBeamSpotManager.h.

Referenced by AlcaBeamSpotManager(), and readLumi().

sigmaZCut_

double AlcaBeamSpotManager::sigmaZCut_

private

Definition at line 43 of file AlcaBeamSpotManager.h.

Referenced by createWeightedPayloads().