EcalDccWeightBuilder Class Reference

#include <EcalDccWeightBuilder.h>

Inheritance diagram for EcalDccWeightBuilder:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &es) override
	EcalDccWeightBuilder (edm::ParameterSet const &ps)
	~EcalDccWeightBuilder () override
Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
bool	wantsInputProcessBlocks () const final
bool	wantsProcessBlocks () const final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Types
enum	mode_t { WEIGHTS_FROM_CONFIG, COMPUTE_WEIGHTS }

Private Member Functions
void	computeAllWeights (bool withIntercalib, const edm::EventSetup &es)
void	computeWeights (const EcalShapeBase &shape, int binOfMax, double timePhase, int iFirst0, int nWeights, int iSkip0, std::vector< double > &result)
void	dbId (const DetId &detId, int &fedId, int &smId, int &ruId, int &xtalId) const
double	decodeWeight (int W)
int	encodeWeight (double w)
double	intercalib (const DetId &detId)
template<class T >
void	sort (const std::vector< T > &a, std::vector< int > &s, bool decreasingOrder=false)
void	unbiasWeights (std::vector< double > &weights, std::vector< int32_t > *encodedWeigths)
void	writeWeightToAsciiFile ()
void	writeWeightToDB ()
void	writeWeightToRootFile ()

Private Attributes
std::string	asciiOutputFileName_
EcalIntercalibConstantMap &	calibMap_
std::string	dbPassword_
std::string	dbSid_
std::string	dbTag_
std::string	dbUser_
int	dbVersion_
int	dcc1stSample_
bool	dccWeightsWithIntercalib_
EBShape	ebShape_
const EcalElectronicsMapping *	ecalElectronicsMap_
EEShape	eeShape_
EcalIntercalibConstantMap	emptyCalibMap_
std::map< DetId, std::vector< int > >	encodedWeights_
edm::ESHandle< CaloGeometry >	geom_
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	geometryToken_
mode_t	imode_
std::vector< double >	inputWeights_
const edm::ESGetToken< EcalIntercalibConstants, EcalIntercalibConstantsRcd >	intercalibConstToken_
const edm::ESGetToken< EcalElectronicsMapping, EcalMappingRcd >	mappingToken_
std::string	mode_
int	nDccWeights_
std::string	rootOutputFileName_
int	sampleToSkip_
bool	sqlMode_
bool	writeToAsciiFile_
bool	writeToDB_
bool	writeToRootFile_

Static Private Attributes
static const int	ecalDccFedIdMax = 654
static const int	ecalDccFedIdMin = 601
static const int	nDccs = ecalDccFedIdMax - ecalDccFedIdMin + 1
static const double	weightScale_ = 1024.

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 27 of file EcalDccWeightBuilder.h.

Member Enumeration Documentation

mode_t

enum EcalDccWeightBuilder::mode_t

private

Enumerator
WEIGHTS_FROM_CONFIG
COMPUTE_WEIGHTS

Definition at line 29 of file EcalDccWeightBuilder.h.

{ WEIGHTS_FROM_CONFIG, COMPUTE_WEIGHTS };

Constructor & Destructor Documentation

EcalDccWeightBuilder()

EcalDccWeightBuilder::EcalDccWeightBuilder

(

edm::ParameterSet const &

ps

)

Constructs an EcalDccWeightBuilder

Parameters

ps	CMSSW mondule configuration

Definition at line 38 of file EcalDccWeightBuilder.cc.

References COMPUTE_WEIGHTS, Exception, imode_, inputWeights_, mode_, nDccWeights_, and WEIGHTS_FROM_CONFIG.

    : dcc1stSample_(ps.getParameter<int>("dcc1stSample")),
      sampleToSkip_(ps.getParameter<int>("sampleToSkip")),
      nDccWeights_(ps.getParameter<int>("nDccWeights")),
      inputWeights_(ps.getParameter<vector<double> >("inputWeights")),
      mode_(ps.getParameter<string>("mode")),
      dccWeightsWithIntercalib_(ps.getParameter<bool>("dccWeightsWithIntercalib")),
      writeToDB_(ps.getParameter<bool>("writeToDB")),
      writeToAsciiFile_(ps.getParameter<bool>("writeToAsciiFile")),
      writeToRootFile_(ps.getParameter<bool>("writeToRootFile")),
      asciiOutputFileName_(ps.getParameter<string>("asciiOutputFileName")),
      rootOutputFileName_(ps.getParameter<string>("rootOutputFileName")),
      dbSid_(ps.getParameter<string>("dbSid")),
      dbUser_(ps.getParameter<string>("dbUser")),
      dbPassword_(ps.getUntrackedParameter<string>("dbPassword", "")),
      dbTag_(ps.getParameter<string>("dbTag")),
      dbVersion_(ps.getParameter<int>("dbVersion")),
      sqlMode_(ps.getParameter<bool>("sqlMode")),
      geometryToken_(esConsumes()),
      mappingToken_(esConsumes()),
      intercalibConstToken_(esConsumes()),
      calibMap_(emptyCalibMap_),
      ebShape_(consumesCollector()),
      eeShape_(consumesCollector()) {
  if (mode_ == "weightsFromConfig") {
    imode_ = WEIGHTS_FROM_CONFIG;
    if (inputWeights_.size() != (unsigned)nDccWeights_) {
      throw cms::Exception("Config") << "Inconsistent configuration. 'nDccWeights' parameters indicated "
                                     << nDccWeights_ << " weights while parameter 'inputWeights_' contains "
                                     << inputWeights_.size() << " weight values!\n";
    }
  } else if (mode_ == "computeWeights") {
    imode_ = COMPUTE_WEIGHTS;
  } else {
    throw cms::Exception("Config") << "Invalid value ('" << mode_ << "') for parameter mode. "
                                   << "Valid values are: 'weightsFromConfig' and 'computeWeights'\n";
  }
}

~EcalDccWeightBuilder()

EcalDccWeightBuilder::~EcalDccWeightBuilder ( )

inlineoverride

Destructor

Definition at line 40 of file EcalDccWeightBuilder.h.

{};

Member Function Documentation

analyze()

void EcalDccWeightBuilder::analyze ( const edm::Event &  event, const edm::EventSetup &  es  )

overridevirtual

Analyze method called by the event loop.

Parameters

event	CMSSW event
es	event setup

Implements edm::one::EDAnalyzerBase.

Definition at line 77 of file EcalDccWeightBuilder.cc.

References calibMap_, computeAllWeights(), dccWeightsWithIntercalib_, ecalElectronicsMap_, geom_, geometryToken_, edm::EventSetup::getData(), edm::EventSetup::getHandle(), EcalCondObjectContainer< T >::getMap(), intercalibConstToken_, mappingToken_, writeToAsciiFile_, writeToDB_, writeToRootFile_, writeWeightToAsciiFile(), writeWeightToDB(), and writeWeightToRootFile().

                                                                                 {
  const auto mappingHandle = es.getHandle(mappingToken_);
  ecalElectronicsMap_ = mappingHandle.product();

  // Retrieval of intercalib constants
  if (dccWeightsWithIntercalib_) {
    const auto& intercalibConst = es.getData(intercalibConstToken_);
    calibMap_ = intercalibConst.getMap();
  }

  //gets geometry
  geom_ = es.getHandle(geometryToken_);

  //computes the weights:
  computeAllWeights(dccWeightsWithIntercalib_, es);

  //Writing out weights.
  if (writeToAsciiFile_)
    writeWeightToAsciiFile();
  if (writeToRootFile_)
    writeWeightToRootFile();
  if (writeToDB_)
    writeWeightToDB();
}

computeAllWeights()

void EcalDccWeightBuilder::computeAllWeights ( bool  withIntercalib, const edm::EventSetup &  es  )

private

Definition at line 102 of file EcalDccWeightBuilder.cc.

References cms::cuda::assert(), CaloSimParameters::binOfMaximum(), COMPUTE_WEIGHTS, computeWeights(), filterCSVwithJSON::copy, dcc1stSample_, MillePedeFileConverter_cfg::e, ebShape_, DetId::Ecal, EcalBarrel, EcalEndcap, eeShape_, encodedWeights_, cppFunctionSkipper::exception, Exception, geom_, CaloGeometry::getValidDetIds(), mps_fire::i, hcalRecHitTable_cff::ieta, imode_, inputWeights_, intercalib(), hcalRecHitTable_cff::iphi, SiStripPI::min, nDccWeights_, sampleToSkip_, EcalSimParameterMap::simParameters(), CaloSimParameters::timePhase(), unbiasWeights(), w(), and WEIGHTS_FROM_CONFIG.

Referenced by analyze().

                                                                                         {
  const int nw = nDccWeights_;
  int iSkip0_ = sampleToSkip_ >= 0 ? (sampleToSkip_ - dcc1stSample_) : -1;

  EcalSimParameterMap parameterMap;
  const vector<DetId>& ebDetIds = geom_->getValidDetIds(DetId::Ecal, EcalBarrel);
  const vector<DetId>& eeDetIds = geom_->getValidDetIds(DetId::Ecal, EcalEndcap);

  vector<DetId> detIds(ebDetIds.size() + eeDetIds.size());
  copy(ebDetIds.begin(), ebDetIds.end(), detIds.begin());
  copy(eeDetIds.begin(), eeDetIds.end(), detIds.begin() + ebDetIds.size());

  vector<double> baseWeights(nw);  //weight obtained from signal shape
  vector<double> w(nw);            //weight*intercalib
  vector<int> W(nw);               //weight in hw encoding (integrer)
  double prevPhase = numeric_limits<double>::min();

  if (imode_ == WEIGHTS_FROM_CONFIG) {
    assert(inputWeights_.size() == baseWeights.size());
    copy(inputWeights_.begin(), inputWeights_.end(), baseWeights.begin());
  }

  for (vector<DetId>::const_iterator it = detIds.begin(); it != detIds.end(); ++it) {
    double phase = parameterMap.simParameters(*it).timePhase();
    int binOfMax = parameterMap.simParameters(*it).binOfMaximum();

#if 0
    //for debugging...
    edm::LogVerbatim("EcalDccWeightBuilder") << __FILE__ << ":" << __LINE__ << ": ";
    if(it->subdetId()==EcalBarrel){
      edm::LogVerbatim("EcalDccWeightBuilder") << "ieta = " << setw(4) << ((EBDetId)(*it)).ieta()
           << " iphi = " << setw(4) << ((EBDetId)(*it)).iphi() << " ";
    } else if(it->subdetId()==EcalEndcap){
      edm::LogVerbatim("EcalDccWeightBuilder") << "ix = " << setw(3) << ((EEDetId)(*it)).ix()
           << " iy = " << setw(3) << ((EEDetId)(*it)).iy()
           << " iz = " << setw(1) << ((EEDetId)(*it)).iy() << " ";
    } else{
      throw cms::Exception("EcalDccWeightBuilder")
        << "Bug found in " << __FILE__ << ":" << __LINE__ << ": "
        << "Got a detId which is neither tagged as ECAL Barrel "
        << "not ECAL endcap while looping on ECAL cell detIds\n";
    }
    edm::LogVerbatim("EcalDccWeightBuilder") << " -> phase: "  << phase << "\n";
    edm::LogVerbatim("EcalDccWeightBuilder") << " -> binOfMax: " << binOfMax << "\n";
#endif

    try {
      EcalShapeBase* pShape;

      if (it->subdetId() == EcalBarrel) {
        pShape = &ebShape_;
      } else if (it->subdetId() == EcalEndcap) {
        pShape = &eeShape_;
      } else {
        throw cms::Exception("EcalDccWeightBuilder") << "Bug found in " << __FILE__ << ":" << __LINE__ << ": "
                                                     << "Got a detId which is neither tagged as ECAL Barrel "
                                                     << "not ECAL endcap while looping on ECAL cell detIds\n";
      }

      if (phase != prevPhase) {
        if (imode_ == COMPUTE_WEIGHTS) {
          if (it->subdetId() == EcalBarrel) {
            computeWeights(*pShape, binOfMax, phase, dcc1stSample_ - 1, nDccWeights_, iSkip0_, baseWeights);
          }
          prevPhase = phase;
        }
      }
      for (int i = 0; i < nw; ++i) {
        w[i] = baseWeights[i];
        if (withIntercalib)
          w[i] *= intercalib(*it);
      }
      unbiasWeights(w, &W);
      encodedWeights_[*it] = W;
    } catch (std::exception& e) {
      edm::LogVerbatim("EcalDccWeightBuilder") << __FILE__ << ":" << __LINE__ << ": ";
      if (it->subdetId() == EcalBarrel) {
        edm::LogVerbatim("EcalDccWeightBuilder") << "ieta = " << setw(4) << ((EBDetId)(*it)).ieta()
                                                 << " iphi = " << setw(4) << ((EBDetId)(*it)).iphi() << " ";
      } else if (it->subdetId() == EcalEndcap) {
        edm::LogVerbatim("EcalDccWeightBuilder")
            << "ix = " << setw(3) << ((EEDetId)(*it)).ix() << " iy = " << setw(3) << ((EEDetId)(*it)).iy()
            << " iz = " << setw(1) << ((EEDetId)(*it)).iy() << " ";
      } else {
        edm::LogVerbatim("EcalDccWeightBuilder") << "DetId " << (uint32_t)(*it);
      }
      edm::LogVerbatim("EcalDccWeightBuilder") << "phase: " << phase << "\n";
      throw;
    }
  }
}

computeWeights()

void EcalDccWeightBuilder::computeWeights ( const EcalShapeBase &  shape, int  binOfMax, double  timePhase, int  iFirst0, int  nWeights, int  iSkip0, std::vector< double > &  result  )

private

Weight computation

Parameters

shape	signal shape to use for weight computation
binOfMax	time bin which shall contain the pulse peak
timePhase	signal phase. First sample corresponds to t = ((binOfMax-1)*25+timePhase) ns in the time reference of EcalShapeBase.
iFirst	first sample the weights must be applied to
nWeights	number of weights
iSkip0	if greater than 0, the corresponding sample will not be used (weights forced to 0).
result	[out] vector filled with computed weights. The vector is resized to the number of weights

Definition at line 194 of file EcalDccWeightBuilder.cc.

References mps_fire::i, mps_fire::result, alignCSCRings::s, l1trig_cff::shape, ecalSimParameterMap_cff::timePhase, and tzero.

Referenced by computeAllWeights().

                                                                  {
  double sum2 = 0.;
  double sum = 0;
  result.resize(nWeights);

  int nActualWeights = 0;

  const double tzero = -(binOfMax - 1) * 25 + timePhase + shape.timeToRise();  //ns

  for (int i = 0; i < nWeights; ++i) {
    double t_ns = tzero + (iFirst + i) * 25;
    double s = shape(t_ns);
    if (i == iSkip) {
      continue;
    }
    result[i] = s;
    sum += s;
    sum2 += s * s;
    ++nActualWeights;
  }
  for (int i = 0; i < nWeights; ++i) {
    if (i == iSkip) {
      result[i] = 0;
    } else {
      result[i] = (result[i] - sum / nActualWeights) / (sum2 - sum * sum / nActualWeights);
    }
  }
}

dbId()

void EcalDccWeightBuilder::dbId ( const DetId &  detId, int &  fedId, int &  smId, int &  ruId, int &  xtalId  ) const

private

Definition at line 559 of file EcalDccWeightBuilder.cc.

References EcalElectronicsId::dccId(), hcalRecHitTable_cff::detId, EcalBarrel, ecalElectronicsMap_, l1tstage2_dqm_sourceclient-live_cfg::fedId, EcalElectronicsMapping::getElectronicsId(), EcalElectronicsId::stripId(), EcalElectronicsId::towerId(), and EcalElectronicsId::xtalId().

Referenced by writeWeightToAsciiFile(), writeWeightToDB(), and writeWeightToRootFile().

                                                                                                       {
  const EcalElectronicsId& elecId = ecalElectronicsMap_->getElectronicsId(detId);

  fedId = 600 + elecId.dccId();
  ruId = ecalElectronicsMap_->getElectronicsId(detId).towerId();

  if (detId.subdetId() == EcalBarrel) {
    smId = ((EBDetId)detId).ism();
  } else {
    smId = 10000 - fedId;  //no SM in EE. Use some unique value to satisfy
    //              current DB PK constraints.
  }
  const int stripLength = 5;  //1 strip = 5 crystals in a row
  xtalId = (elecId.stripId() - 1) * stripLength + elecId.xtalId();

#if 0
  edm::LogVerbatim("EcalDccWeightBuilder") << __FILE__ << ":" << __LINE__ << ": FED ID "
       <<  fedId << "\n";

  edm::LogVerbatim("EcalDccWeightBuilder") << __FILE__ << ":" << __LINE__ << ": SM logical ID "
       <<  smId << "\n";
  
  edm::LogVerbatim("EcalDccWeightBuilder") << __FILE__ << ":" << __LINE__ << ": RU ID (TT or SC): "
       <<  ruId << "\n";
  
  edm::LogVerbatim("EcalDccWeightBuilder") << __FILE__ << ":" << __LINE__ << ": strip:"
       <<  elecId.stripId() << "\n";
  
  edm::LogVerbatim("EcalDccWeightBuilder") << __FILE__ << ":" << __LINE__ << ": xtal in strip: "
       <<  elecId.xtalId() << "\n";

  edm::LogVerbatim("EcalDccWeightBuilder") << __FILE__ << ":" << __LINE__ << ": xtalId in RU: "
       <<  xtalId << "\n";
#endif
}

decodeWeight()

double EcalDccWeightBuilder::decodeWeight ( int  W )

private

Definition at line 231 of file EcalDccWeightBuilder.cc.

References weightScale_.

Referenced by unbiasWeights().

{ return ((double)W) / weightScale_; }

encodeWeight()

int EcalDccWeightBuilder::encodeWeight ( double  w )

private

Definition at line 229 of file EcalDccWeightBuilder.cc.

References w(), and weightScale_.

Referenced by unbiasWeights().

{ return lround(w * weightScale_); }

intercalib()

double EcalDccWeightBuilder::intercalib ( const DetId &  detId )

private

Retrieve intercalibration coefficent for channel detId.

Parameters

detId

ID of the channel the intercalibration coef. must be retrieved for.

Definition at line 308 of file EcalDccWeightBuilder.cc.

References calibMap_, hcalRecHitTable_cff::detId, EcalBarrel, EcalCondObjectContainer< T >::end(), and EcalCondObjectContainer< T >::find().

Referenced by computeAllWeights().

                                                          {
  // get current intercalibration coeff
  double coef;
  EcalIntercalibConstantMap::const_iterator itCalib = calibMap_.find(detId.rawId());
  if (itCalib != calibMap_.end()) {
    coef = (*itCalib);
  } else {
    coef = 1.;
    edm::LogVerbatim("EcalDccWeightBuilder") << (uint32_t)detId << " not found in EcalIntercalibConstantMap";
  }
#if 0
  edm::LogVerbatim("EcalDccWeightBuilder") << __FILE__ << ":" << __LINE__ << ": ";
  if(detId.subdetId()==EcalBarrel){
    edm::LogVerbatim("EcalDccWeightBuilder") <<  "ieta = " << ((EBDetId)detId).ieta()
     << " iphi = " << ((EBDetId)detId).iphi();
  } else{
    edm::LogVerbatim("EcalDccWeightBuilder") << "ix = " << ((EEDetId)detId).ix()
     << " iy = " << ((EEDetId)detId).iy()
     << " iz = " << ((EEDetId)detId).zside();
  }
  edm::LogVerbatim("EcalDccWeightBuilder") << " coef = " <<  coef << "\n";
#endif
  return coef;
}

sort()

template<class T >

void EcalDccWeightBuilder::sort ( const std::vector< T > &  a, std::vector< int > &  s, bool  decreasingOrder = false  )

private

Vector index sorting.

Parameters

	vector	of values
[out]	sorted	index of the vector
	decreasingOrder	switch to use decreasing order instead of default increasing order

Definition at line 234 of file EcalDccWeightBuilder.cc.

References a, mps_fire::i, dqmiolumiharvest::j, alignCSCRings::s, and edm::swap().

Referenced by unbiasWeights().

                                                                                              {
  //performs a bubble sort: adjacent elements are successively swapped 2 by 2
  //until the list is finally sorted.
  bool changed = false;
  s.resize(a.size());
  for (unsigned i = 0; i < a.size(); ++i)
    s[i] = i;
  if (a.empty())
    return;
  do {
    changed = false;
    for (unsigned i = 0; i < a.size() - 1; ++i) {
      const int j = s[i];
      const int nextj = s[i + 1];
      if ((decreasingOrder && (a[j] < a[nextj])) || (!decreasingOrder && (a[j] > a[nextj]))) {
        std::swap(s[i], s[i + 1]);
        changed = true;
      }
    }
  } while (changed);
}

unbiasWeights()

void EcalDccWeightBuilder::unbiasWeights ( std::vector< double > &  weights, std::vector< int32_t > *  encodedWeigths  )

private

Definition at line 256 of file EcalDccWeightBuilder.cc.

References decodeWeight(), encodeWeight(), mps_fire::i, sort(), and hltDeepSecondaryVertexTagInfosPFPuppi_cfi::weights.

Referenced by computeAllWeights().

                                                                                                   {
  const unsigned nw = weights.size();

  //computes integer weights, weights residuals and weight sum residual:
  vector<double> dw(nw);  //weight residuals due to interger encoding
  vector<int> W(nw);      //integer weights
  int wsum = 0;
  for (unsigned i = 0; i < nw; ++i) {
    W[i] = encodeWeight(weights[i]);
    dw[i] = decodeWeight(W[i]) - weights[i];
    wsum += W[i];
  }

  //sorts weight residuals in decreasing order:
  vector<int> iw(nw);
  sort(dw, iw, true);

  //compensates weight sum residual by adding or substracting 1 to weights
  //starting from:
  // 1) the weight with the minimal signed residual if the correction
  // is positive (wsum<0)
  // 2) the weight with the maximal signed residual if the correction
  // is negative (wsum>0)
  int wsumSign = wsum > 0 ? 1 : -1;
  int i = wsum > 0 ? 0 : (nw - 1);
  while (wsum != 0) {
    W[iw[i]] -= wsumSign;
    wsum -= wsumSign;
    i += wsumSign;
    if (i < 0 || i >= (int)nw) {  //recompute the residuals if a second iteration is
      // needed (in principle, it is not expected with usual input weights), :
      for (unsigned i = 0; i < nw; ++i) {
        dw[i] = decodeWeight(W[i]) - weights[i];
        sort(dw, iw, true);
      }
    }
    if (i < 0)
      i = nw - 1;
    if (i >= (int)nw)
      i = 0;
  }

  //copy result
  if (encodedWeights != nullptr)
    encodedWeights->resize(nw);
  for (unsigned i = 0; i < nw; ++i) {
    weights[i] = decodeWeight(W[i]);
    if (encodedWeights)
      (*encodedWeights)[i] = W[i];
  }
}

writeWeightToAsciiFile()

void EcalDccWeightBuilder::writeWeightToAsciiFile ( )

private

Definition at line 333 of file EcalDccWeightBuilder.cc.

References asciiOutputFileName_, dbId(), dbTag_, dbVersion_, dcc1stSample_, hcalRecHitTable_cff::detId, ecalDccFedIdMax, ecalDccFedIdMin, encodedWeights_, Exception, l1tstage2_dqm_sourceclient-live_cfg::fedId, geometryDiff::file, fileCollector2::fName, mps_fire::i, genParticles_cff::map, sqlMode_, and hltDeepSecondaryVertexTagInfosPFPuppi_cfi::weights.

Referenced by analyze().

                                                  {
  string fName = !asciiOutputFileName_.empty() ? asciiOutputFileName_.c_str() : "dccWeights.txt";
  ofstream file(fName.c_str());
  if (!file.good()) {
    throw cms::Exception("Output") << "Failed to open file '" << fName << "'for writing DCC weights\n";
  }

  const char* comment = sqlMode_ ? "-- " : "# ";

  file << comment << "List of weights for amplitude estimation to be used in DCC for\n"
       << comment << "zero suppresssion.\n\n";
  if (!sqlMode_) {
    file << comment << "Note: RU: trigger tower in EB, supercrystal in EE\n"
         << comment << "      xtl: crystal electronic channel id in RU, from 1 to 25\n\n"
         << comment << " DetId    SM  FED RU xtl weights[0..5]...\n";
  }

  if (sqlMode_) {
    file << "variable recid number;\n"
            "exec select COND2CONF_INFO_SQ.NextVal into :recid from DUAL;\n"
            "insert into weights_info (rec_id,tag,version) values (:recid,'"
         << dbTag_ << "'," << dbVersion_ << ");\n";
    file << "\n"
         << comment
         << "index of first sample used in the weighting sum\n"
            "begin\n"
            "  for fedid in "
         << ecalDccFedIdMin << ".." << ecalDccFedIdMax
         << " loop\n"
            "    insert into dcc_weightsample_dat (rec_id, logic_id, sample_id, \n"
            "    weight_number)\n"
            "    values(:recid,fedid,"
         << dcc1stSample_
         << ",1);\n"
            "  end loop;\n"
            "end;\n"
            "/\n";
  } else {
    file << "1st DCC sample: " << dcc1stSample_ << "\n";
  }

  file << "\n" << comment << "list of weights per crystal channel\n";

  for (map<DetId, std::vector<int32_t> >::const_iterator it = encodedWeights_.begin(); it != encodedWeights_.end();
       ++it) {
    const DetId& detId = it->first;

    int fedId;
    int smId;
    int ruId;
    int xtalId;

    //detId ->  fedId, smId, ruId, xtalId
    dbId(detId, fedId, smId, ruId, xtalId);

    char delim = sqlMode_ ? ',' : ' ';

    if (sqlMode_)
      file << "-- detId " << detId.rawId() << "\n"
           << "insert into dcc_weights_dat(rec_id,sm_id,fed_id,"
              "tt_id, cry_id,\n"
              "weight_0,weight_1,weight_2,weight_3,weight_4,weight_5) \n"
              "values ("
              ":recid";

    const vector<int>& weights = it->second;
    if (!sqlMode_)
      file << setw(10) << detId.rawId();
    file << delim << setw(2) << smId;
    file << delim << setw(3) << fedId;
    file << delim << setw(2) << ruId;
    file << delim << setw(2) << xtalId;

    for (unsigned i = 0; i < weights.size(); ++i) {
      file << delim << setw(5) << weights[i];
    }
    if (sqlMode_)
      file << ");";
    file << "\n";
  }
  if (!file.good()) {
    throw cms::Exception("Output") << "Error while writing DCC weights to '" << fName << "' file.";
  }
}

writeWeightToDB()

void EcalDccWeightBuilder::writeWeightToDB ( )

private

Definition at line 470 of file EcalDccWeightBuilder.cc.

References dbId(), dbPassword_, dbSid_, dbTag_, dbUser_, dbVersion_, dcc1stSample_, hcalRecHitTable_cff::detId, MillePedeFileConverter_cfg::e, ecalDccFedIdMax, ecalDccFedIdMin, test_db_connect::econn, encodedWeights_, beamvalidation::exit(), l1tstage2_dqm_sourceclient-live_cfg::fedId, MillePedeFileConverter_cfg::fileName, ODFEWeightsInfo::getId(), mps_fire::i, genParticles_cff::map, nDccs, IODConfig::setConfigTag(), ODWeightsDat::setCrystalId(), ODWeightsDat::setFedId(), ODWeightsDat::setId(), ODWeightsDat::setSMId(), ODWeightsDat::setTTId(), ODFEWeightsInfo::setVersion(), ODWeightsDat::setWeight0(), ODWeightsDat::setWeight1(), ODWeightsDat::setWeight2(), ODWeightsDat::setWeight3(), ODWeightsDat::setWeight4(), ODWeightsDat::setWeight5(), AlCaHLTBitMon_QueryRunRegistry::string, w(), and hltDeepSecondaryVertexTagInfosPFPuppi_cfi::weights.

Referenced by analyze().

                                           {
  edm::LogVerbatim("EcalDccWeightBuilder") << "going to write to the online DB " << dbSid_ << " user " << dbUser_;
  ;
  EcalCondDBInterface* econn;

  try {
    edm::LogVerbatim("EcalDccWeightBuilder") << "Making connection..." << flush;
    const string& filePrefix = string("file:");
    if (dbPassword_.find(filePrefix) == 0) {  //password must be read for a file
      string fileName = dbPassword_.substr(filePrefix.size());
      //substitute dbPassword_ value by the password read from the file
      PasswordReader pr;
      pr.readPassword(fileName, dbUser_, dbPassword_);
    }

    econn = new EcalCondDBInterface(dbSid_, dbUser_, dbPassword_);
    edm::LogVerbatim("EcalDccWeightBuilder") << "Done.";
  } catch (runtime_error& e) {
    edm::LogError("dbconnection") << e.what();
    exit(-1);
  }

  ODFEWeightsInfo weight_info;
  weight_info.setConfigTag(dbTag_);
  weight_info.setVersion(dbVersion_);
  edm::LogVerbatim("EcalDccWeightBuilder") << "Inserting in DB...";

  econn->insertConfigSet(&weight_info);

  int weight_id = weight_info.getId();
  edm::LogVerbatim("EcalDccWeightBuilder") << "WeightInfo inserted with ID " << weight_id;

  vector<ODWeightsDat> datadel;
  datadel.reserve(encodedWeights_.size());

  vector<ODWeightsSamplesDat> dcc1stSampleConfig(nDccs);
  for (int i = ecalDccFedIdMin; i <= ecalDccFedIdMax; ++i) {
    dcc1stSampleConfig[i].setId(weight_id);
    dcc1stSampleConfig[i].setFedId(601 + i);
    dcc1stSampleConfig[i].setSampleId(dcc1stSample_);
    dcc1stSampleConfig[i].setWeightNumber(-1);  //not used.
  }
  econn->insertConfigDataArraySet(dcc1stSampleConfig, &weight_info);

  for (map<DetId, std::vector<int32_t> >::const_iterator it = encodedWeights_.begin(); it != encodedWeights_.end();
       ++it) {
    const DetId& detId = it->first;
    const unsigned nWeights = 6;
    vector<int> weights(nWeights);

    for (unsigned i = 0; i < weights.size(); ++i) {
      //completing the weight vector with zeros in case it has
      //less than 6 elements:
      const vector<int>& w = it->second;
      weights[i] = i < w.size() ? w[i] : 0;
    }

    ODWeightsDat one_dat;
    one_dat.setId(weight_id);

    int fedId;
    int smId;
    int ruId;
    int xtalId;

    //detId ->  fedId, smId, ruId, xtalId
    dbId(detId, fedId, smId, ruId, xtalId);

    one_dat.setSMId(smId);
    one_dat.setFedId(fedId);
    one_dat.setTTId(ruId);
    one_dat.setCrystalId(xtalId);

    one_dat.setWeight0(weights[0]);
    one_dat.setWeight1(weights[1]);
    one_dat.setWeight2(weights[2]);
    one_dat.setWeight3(weights[3]);
    one_dat.setWeight4(weights[4]);
    one_dat.setWeight5(weights[5]);

    datadel.push_back(one_dat);
  }
  econn->insertConfigDataArraySet(datadel, &weight_info);
  edm::LogVerbatim("EcalDccWeightBuilder") << " .. done insertion in DB ";
  delete econn;
  edm::LogVerbatim("EcalDccWeightBuilder") << "closed DB connection ... done";
}

writeWeightToRootFile()

void EcalDccWeightBuilder::writeWeightToRootFile ( )

private

Definition at line 417 of file EcalDccWeightBuilder.cc.

References visDQMUpload::buf, filterCSVwithJSON::copy, dbId(), hcalRecHitTable_cff::detId, encodedWeights_, Exception, l1tstage2_dqm_sourceclient-live_cfg::fedId, geometryDiff::file, fileCollector2::fName, genParticles_cff::map, rootOutputFileName_, submitPVValidationJobs::t, and hltDeepSecondaryVertexTagInfosPFPuppi_cfi::weights.

Referenced by analyze().

                                                 {
  string fName = !rootOutputFileName_.empty() ? rootOutputFileName_.c_str() : "dccWeights.root";
  TFile file(fName.c_str(), "RECREATE");
  if (file.IsZombie()) {
    throw cms::Exception("Output") << "Failed to open file '" << fName << "'for writing DCC weights\n";
  }
  TTree t("dccWeights", "Weights for DCC ZS filter");
  const int nWeightMax = 20;  //normally n_weights = 6. A different might be used
  //                           used for test purposes.
  struct {
    Int_t detId;
    Int_t fedId;
    Int_t smId;
    Int_t ruId;
    Int_t xtalId;
    Int_t n_weights;
    Int_t weights[nWeightMax];
  } buf;
  t.Branch("weights",
           &buf,
           "rawDetId/I:"
           "feId/I:"
           "smSlotId/I:"
           "ruId/I:"
           "xtalInRuId/I:"
           "n_weights/I:"
           "weights[n_weights]/I");
  for (map<DetId, std::vector<int32_t> >::const_iterator it = encodedWeights_.begin(); it != encodedWeights_.end();
       ++it) {
    buf.detId = it->first.rawId();
    buf.n_weights = it->second.size();

    //detId ->  fedId, smId, ruId, xtalId
    dbId(buf.detId, buf.fedId, buf.smId, buf.ruId, buf.xtalId);

    if (buf.n_weights > nWeightMax) {
      throw cms::Exception("EcalDccWeight") << "Number of weights (" << buf.n_weights << ") for DetId " << buf.detId
                                            << " exceeded maximum limit (" << nWeightMax << ") of root output format. ";
    }
    copy(it->second.begin(), it->second.end(), buf.weights);
    t.Fill();
  }
  t.Write();
  file.Close();
}

Member Data Documentation

asciiOutputFileName_

std::string EcalDccWeightBuilder::asciiOutputFileName_

private

Definition at line 117 of file EcalDccWeightBuilder.h.

Referenced by writeWeightToAsciiFile().

calibMap_

EcalIntercalibConstantMap& EcalDccWeightBuilder::calibMap_

private

Definition at line 132 of file EcalDccWeightBuilder.h.

Referenced by analyze(), and intercalib().

dbPassword_

std::string EcalDccWeightBuilder::dbPassword_

private

Definition at line 121 of file EcalDccWeightBuilder.h.

Referenced by writeWeightToDB().

dbSid_

std::string EcalDccWeightBuilder::dbSid_

private

Definition at line 119 of file EcalDccWeightBuilder.h.

Referenced by writeWeightToDB().

dbTag_

std::string EcalDccWeightBuilder::dbTag_

private

Definition at line 122 of file EcalDccWeightBuilder.h.

Referenced by writeWeightToAsciiFile(), and writeWeightToDB().

dbUser_

std::string EcalDccWeightBuilder::dbUser_

private

Definition at line 120 of file EcalDccWeightBuilder.h.

Referenced by writeWeightToDB().

dbVersion_

int EcalDccWeightBuilder::dbVersion_

private

Definition at line 123 of file EcalDccWeightBuilder.h.

Referenced by writeWeightToAsciiFile(), and writeWeightToDB().

dcc1stSample_

int EcalDccWeightBuilder::dcc1stSample_

private

Definition at line 105 of file EcalDccWeightBuilder.h.

Referenced by computeAllWeights(), writeWeightToAsciiFile(), and writeWeightToDB().

dccWeightsWithIntercalib_

bool EcalDccWeightBuilder::dccWeightsWithIntercalib_

private

Definition at line 113 of file EcalDccWeightBuilder.h.

Referenced by analyze().

ebShape_

EBShape EcalDccWeightBuilder::ebShape_

private

Definition at line 136 of file EcalDccWeightBuilder.h.

Referenced by computeAllWeights().

ecalDccFedIdMax

const int EcalDccWeightBuilder::ecalDccFedIdMax = 654

staticprivate

Definition at line 143 of file EcalDccWeightBuilder.h.

Referenced by writeWeightToAsciiFile(), and writeWeightToDB().

ecalDccFedIdMin

const int EcalDccWeightBuilder::ecalDccFedIdMin = 601

staticprivate

Definition at line 142 of file EcalDccWeightBuilder.h.

Referenced by writeWeightToAsciiFile(), and writeWeightToDB().

ecalElectronicsMap_

const EcalElectronicsMapping* EcalDccWeightBuilder::ecalElectronicsMap_

private

Definition at line 140 of file EcalDccWeightBuilder.h.

Referenced by analyze(), and dbId().

eeShape_

EEShape EcalDccWeightBuilder::eeShape_

private

Definition at line 137 of file EcalDccWeightBuilder.h.

Referenced by computeAllWeights().

emptyCalibMap_

EcalIntercalibConstantMap EcalDccWeightBuilder::emptyCalibMap_

private

Definition at line 133 of file EcalDccWeightBuilder.h.

encodedWeights_

std::map<DetId, std::vector<int> > EcalDccWeightBuilder::encodedWeights_

private

Definition at line 134 of file EcalDccWeightBuilder.h.

Referenced by computeAllWeights(), writeWeightToAsciiFile(), writeWeightToDB(), and writeWeightToRootFile().

geom_

edm::ESHandle<CaloGeometry> EcalDccWeightBuilder::geom_

private

Definition at line 130 of file EcalDccWeightBuilder.h.

Referenced by analyze(), and computeAllWeights().

geometryToken_

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> EcalDccWeightBuilder::geometryToken_

private

Definition at line 126 of file EcalDccWeightBuilder.h.

Referenced by analyze().

imode_

mode_t EcalDccWeightBuilder::imode_

private

Definition at line 112 of file EcalDccWeightBuilder.h.

Referenced by computeAllWeights(), and EcalDccWeightBuilder().

inputWeights_

std::vector<double> EcalDccWeightBuilder::inputWeights_

private

weights used in weightFromConfig mode.

Definition at line 110 of file EcalDccWeightBuilder.h.

Referenced by computeAllWeights(), and EcalDccWeightBuilder().

intercalibConstToken_

const edm::ESGetToken<EcalIntercalibConstants, EcalIntercalibConstantsRcd> EcalDccWeightBuilder::intercalibConstToken_

private

Definition at line 128 of file EcalDccWeightBuilder.h.

Referenced by analyze().

mappingToken_

const edm::ESGetToken<EcalElectronicsMapping, EcalMappingRcd> EcalDccWeightBuilder::mappingToken_

private

Definition at line 127 of file EcalDccWeightBuilder.h.

Referenced by analyze().

mode_

std::string EcalDccWeightBuilder::mode_

private

Definition at line 111 of file EcalDccWeightBuilder.h.

Referenced by EcalDccWeightBuilder().

nDccs

const int EcalDccWeightBuilder::nDccs = ecalDccFedIdMax - ecalDccFedIdMin + 1

staticprivate

Definition at line 144 of file EcalDccWeightBuilder.h.

Referenced by writeWeightToDB().

nDccWeights_

int EcalDccWeightBuilder::nDccWeights_

private

Definition at line 107 of file EcalDccWeightBuilder.h.

Referenced by computeAllWeights(), and EcalDccWeightBuilder().

rootOutputFileName_

std::string EcalDccWeightBuilder::rootOutputFileName_

private

Definition at line 118 of file EcalDccWeightBuilder.h.

Referenced by writeWeightToRootFile().

sampleToSkip_

int EcalDccWeightBuilder::sampleToSkip_

private

Definition at line 106 of file EcalDccWeightBuilder.h.

Referenced by computeAllWeights().

sqlMode_

bool EcalDccWeightBuilder::sqlMode_

private

Definition at line 124 of file EcalDccWeightBuilder.h.

Referenced by writeWeightToAsciiFile().

weightScale_

const double EcalDccWeightBuilder::weightScale_ = 1024.

staticprivate

Definition at line 139 of file EcalDccWeightBuilder.h.

Referenced by decodeWeight(), and encodeWeight().

writeToAsciiFile_

bool EcalDccWeightBuilder::writeToAsciiFile_

private

Definition at line 115 of file EcalDccWeightBuilder.h.

Referenced by analyze().

writeToDB_

bool EcalDccWeightBuilder::writeToDB_

private

Definition at line 114 of file EcalDccWeightBuilder.h.

Referenced by analyze().

writeToRootFile_

bool EcalDccWeightBuilder::writeToRootFile_

private

Definition at line 116 of file EcalDccWeightBuilder.h.

Referenced by analyze().