BeamSpotFakeConditions.cc

Go to the documentation of this file.

//

#include <iostream>
#include <fstream>
#include <sstream>
#include <cstring>
#include <string>
#include <vector>
#include <memory>

#include <TClass.h>

#include "FWCore/Utilities/interface/Exception.h"
#include "FWCore/Framework/interface/ESProducer.h"
#include "FWCore/Framework/interface/EventSetupRecordIntervalFinder.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/SourceFactory.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ParameterSet/interface/FileInPath.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "CondFormats/DataRecord/interface/BeamSpotObjectsRcd.h"
#include "CondFormats/BeamSpotObjects/interface/BeamSpotObjects.h"

class BeamSpotFakeConditions : public edm::ESProducer, public edm::EventSetupRecordIntervalFinder {
public:
  typedef std::unique_ptr<BeamSpotObjects> ReturnType;
  BeamSpotFakeConditions(const edm::ParameterSet &params);
  ~BeamSpotFakeConditions() override;
  ReturnType produce(const BeamSpotObjectsRcd &record);

private:
  void setIntervalFor(const edm::eventsetup::EventSetupRecordKey &key,
                      const edm::IOVSyncValue &syncValue,
                      edm::ValidityInterval &oValidity) override;
  edm::FileInPath inputFilename_;
  bool getDataFromFile_;
  double x, y, z, sigmaZ, dxdz, dydz, beamWidthX, beamWidthY, emittanceX, emittanceY, betastar;
  std::string tag;
  double cov[7][7];
  int type;
};

BeamSpotFakeConditions::BeamSpotFakeConditions(const edm::ParameterSet &params) {
  setWhatProduced(this);
  findingRecord<BeamSpotObjectsRcd>();
  getDataFromFile_ = params.getParameter<bool>("getDataFromFile");
  if (getDataFromFile_) {
    inputFilename_ = params.getParameter<edm::FileInPath>("InputFilename");
    std::ifstream fasciiFile(inputFilename_.fullPath().c_str());
    fasciiFile >> tag >> type;
    fasciiFile >> tag >> x;
    fasciiFile >> tag >> y;
    fasciiFile >> tag >> z;
    fasciiFile >> tag >> sigmaZ;
    fasciiFile >> tag >> dxdz;
    fasciiFile >> tag >> dydz;
    fasciiFile >> tag >> beamWidthX;
    fasciiFile >> tag >> beamWidthY;
    fasciiFile >> tag >> cov[0][0] >> cov[0][1] >> cov[0][2] >> cov[0][3] >> cov[0][4] >> cov[0][5] >> cov[0][6];
    fasciiFile >> tag >> cov[1][0] >> cov[1][1] >> cov[1][2] >> cov[1][3] >> cov[1][4] >> cov[1][5] >> cov[1][6];
    fasciiFile >> tag >> cov[2][0] >> cov[2][1] >> cov[2][2] >> cov[2][3] >> cov[2][4] >> cov[2][5] >> cov[2][6];
    fasciiFile >> tag >> cov[3][0] >> cov[3][1] >> cov[3][2] >> cov[3][3] >> cov[3][4] >> cov[3][5] >> cov[3][6];
    fasciiFile >> tag >> cov[4][0] >> cov[4][1] >> cov[4][2] >> cov[4][3] >> cov[4][4] >> cov[4][5] >> cov[4][6];
    fasciiFile >> tag >> cov[5][0] >> cov[5][1] >> cov[5][2] >> cov[5][3] >> cov[5][4] >> cov[5][5] >> cov[5][6];
    fasciiFile >> tag >> cov[6][0] >> cov[6][1] >> cov[6][2] >> cov[6][3] >> cov[6][4] >> cov[6][5] >> cov[6][6];
    fasciiFile >> tag >> emittanceX;
    fasciiFile >> tag >> emittanceY;
    fasciiFile >> tag >> betastar;

  }
  // input values by hand
  else {
    x = params.getParameter<double>("X0");
    y = params.getParameter<double>("Y0");
    z = params.getParameter<double>("Z0");
    dxdz = params.getParameter<double>("dxdz");
    dydz = params.getParameter<double>("dydz");
    sigmaZ = params.getParameter<double>("sigmaZ");
    beamWidthX = params.getParameter<double>("widthX");
    beamWidthY = params.getParameter<double>("widthY");
    emittanceX = params.getParameter<double>("emittanceX");
    emittanceY = params.getParameter<double>("emittanceY");
    betastar = params.getParameter<double>("betaStar");

    // first set all elements (esp. off-diagonal elements to zero)
    for (int i = 0; i < 7; i++) {
      for (int j = 0; j < 7; j++)
        cov[i][j] = 0.0;
    }

    // we ignore correlations when values are given by hand
    cov[0][0] = pow(params.getParameter<double>("errorX0"), 2);
    cov[1][1] = pow(params.getParameter<double>("errorY0"), 2);
    cov[2][2] = pow(params.getParameter<double>("errorZ0"), 2);
    cov[3][3] = pow(params.getParameter<double>("errorSigmaZ"), 2);
    cov[4][4] = pow(params.getParameter<double>("errordxdz"), 2);
    cov[5][5] = pow(params.getParameter<double>("errordydz"), 2);
    cov[6][6] = pow(params.getParameter<double>("errorWidth"), 2);
  }
}

BeamSpotFakeConditions::~BeamSpotFakeConditions() {}

BeamSpotFakeConditions::ReturnType BeamSpotFakeConditions::produce(const BeamSpotObjectsRcd &record) {
  ReturnType adummy = std::make_unique<BeamSpotObjects>();

  adummy->setPosition(x, y, z);
  adummy->setSigmaZ(sigmaZ);
  adummy->setdxdz(dxdz);
  adummy->setdydz(dydz);
  adummy->setBeamWidthX(beamWidthX);
  adummy->setBeamWidthY(beamWidthY);
  for (int i = 0; i < 7; i++) {
    for (int j = 0; j < 7; j++) {
      adummy->setCovariance(i, j, cov[i][j]);
    }
  }
  adummy->setEmittanceX(emittanceX);
  adummy->setEmittanceY(emittanceY);
  adummy->setBetaStar(betastar);

  return adummy;
}

void BeamSpotFakeConditions::setIntervalFor(const edm::eventsetup::EventSetupRecordKey &key,
                                            const edm::IOVSyncValue &syncValue,
                                            edm::ValidityInterval &oValidity) {
  oValidity = edm::ValidityInterval(edm::IOVSyncValue::beginOfTime(), edm::IOVSyncValue::endOfTime());
}

DEFINE_FWK_EVENTSETUP_SOURCE(BeamSpotFakeConditions);