HcalQIEDataCheck Class Reference

#include <HcalQIEDataCheck.h>

Inheritance diagram for HcalQIEDataCheck:

Public Member Functions
void	analyze (const edm::Event &ev, const edm::EventSetup &es)
	HcalQIEDataCheck (edm::ParameterSet const &ps)
	~HcalQIEDataCheck ()
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
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
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::vector< ModuleDescription const * > &modules, 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
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Private Attributes
bool	checkemapflag
std::string	dumprefs
std::string	dumpupdate
std::string	outfile
bool	validateflag

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
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 ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
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)

Detailed Description

Definition at line 38 of file HcalQIEDataCheck.h.

Constructor & Destructor Documentation

HcalQIEDataCheck::HcalQIEDataCheck

(

edm::ParameterSet const &

ps

)

Definition at line 3 of file HcalQIEDataCheck.cc.

References checkemapflag, dumprefs, dumpupdate, edm::ParameterSet::getUntrackedParameter(), outfile, AlCaHLTBitMon_QueryRunRegistry::string, and validateflag.

{
  outfile = ps.getUntrackedParameter<std::string>("outFile","null");
  dumprefs = ps.getUntrackedParameter<std::string>("dumpRefQIEsTo","null");
  dumpupdate = ps.getUntrackedParameter<std::string>("dumpUpdateQIEsTo","null");
  checkemapflag = ps.getUntrackedParameter<bool>("checkEmap",false);
  validateflag = ps.getUntrackedParameter<bool>("validateQIEs",false);
//  epsilon = ps.getUntrackedParameter<double>("deltaQIE",0);
}

HcalQIEDataCheck::~HcalQIEDataCheck

(

)

Definition at line 13 of file HcalQIEDataCheck.cc.

{}

Member Function Documentation

void HcalQIEDataCheck::analyze	(	const edm::Event &	ev,
		const edm::EventSetup &	es
	)

Definition at line 15 of file HcalQIEDataCheck.cc.

References HcalQIEData::addCoder(), HcalElectronicsMap::allPrecisionId(), checkemapflag, gather_cfg::cout, DEFINE_FWK_MODULE, HcalDbASCIIIO::dumpObject(), dumprefs, dumpupdate, Exception, spr::find(), plotBeamSpotDB::first, edm::EventSetup::get(), HcalCondObjectContainer< Item >::getAllChannels(), HcalQIEData::getCoder(), HcalQIECoder::offset(), outfile, edm::ESHandle< T >::product(), edm::second(), HcalQIECoder::slope(), HcalQIEData::sort(), HcalCondObjectContainerBase::topo(), and validateflag.

{
  using namespace edm::eventsetup;

  edm::ESHandle<HcalQIEData> newQIEs;
  es.get<HcalQIEDataRcd>().get("update",newQIEs);
  const HcalQIEData* myNewQIEs = newQIEs.product();

  edm::ESHandle<HcalQIEData> refQIEs;
  es.get<HcalQIEDataRcd>().get("reference",refQIEs);
  const HcalQIEData* myRefQIEs = refQIEs.product();

  edm::ESHandle<HcalElectronicsMap> refEMap;
  es.get<HcalElectronicsMapRcd>().get("reference",refEMap);
  const HcalElectronicsMap* myRefEMap = refEMap.product();

  if(dumpupdate.compare("null")!=0){
    std::ofstream outStream(dumpupdate.c_str());
    std::cout << "--- Dumping QIEs - update ---" << std::endl;
    HcalDbASCIIIO::dumpObject (outStream, (*myNewQIEs) );
  }
  if(dumprefs.compare("null")!=0){
    std::ofstream outStream2(dumprefs.c_str());
    std::cout << "--- Dumping QIEs - reference ---" << std::endl;
    HcalDbASCIIIO::dumpObject (outStream2, (*myRefQIEs) );
  }

    // first get the list of all channels from the update
    std::vector<DetId> listNewChan = myNewQIEs->getAllChannels();   

    HcalQIEData *resultQIEs = new HcalQIEData(myRefQIEs->topo());
    std::vector<DetId> listRefChan = myRefQIEs->getAllChannels();
    std::vector<DetId>::iterator cell;

    if(validateflag){
    for (std::vector<DetId>::iterator it = listRefChan.begin(); it != listRefChan.end(); it++)
      {
        DetId mydetid = *it;
        cell = std::find(listNewChan.begin(), listNewChan.end(), mydetid);
        if (cell == listNewChan.end()) // not present in new list
          {
        throw cms::Exception("DataDoesNotMatch") << "Value not found in reference" << std::endl;
          }
        else // present in new list
          {

            const HcalQIECoder* first = myNewQIEs->getCoder( mydetid );
            const HcalQIECoder* second = myRefQIEs->getCoder( mydetid );
        {
        bool failflag = false;
            if(first->offset(0,0) != second->offset(0,0)) failflag = true;
            if(first->offset(0,1) != second->offset(0,1)) failflag = true;
            if(first->offset(0,2) != second->offset(0,2)) failflag = true;
            if(first->offset(0,3) != second->offset(0,3)) failflag = true;
            if(first->offset(1,0) != second->offset(1,0)) failflag = true;
            if(first->offset(1,1) != second->offset(1,1)) failflag = true;
            if(first->offset(1,2) != second->offset(1,2)) failflag = true;
            if(first->offset(1,3) != second->offset(1,3)) failflag = true;
            if(first->offset(2,0) != second->offset(2,0)) failflag = true;
            if(first->offset(2,1) != second->offset(2,1)) failflag = true;
            if(first->offset(2,2) != second->offset(2,2)) failflag = true;
            if(first->offset(2,3) != second->offset(2,3)) failflag = true;
            if(first->offset(3,0) != second->offset(3,0)) failflag = true;
            if(first->offset(3,1) != second->offset(3,1)) failflag = true;
            if(first->offset(3,2) != second->offset(3,2)) failflag = true;
            if(first->offset(3,3) != second->offset(3,3)) failflag = true;
            if(first->slope(0,0) != second->slope(0,0)) failflag = true;
            if(first->slope(0,1) != second->slope(0,1)) failflag = true;
            if(first->slope(0,2) != second->slope(0,2)) failflag = true;
            if(first->slope(0,3) != second->slope(0,3)) failflag = true;
            if(first->slope(1,0) != second->slope(1,0)) failflag = true;
            if(first->slope(1,1) != second->slope(1,1)) failflag = true;
            if(first->slope(1,2) != second->slope(1,2)) failflag = true;
            if(first->slope(1,3) != second->slope(1,3)) failflag = true;
            if(first->slope(2,0) != second->slope(2,0)) failflag = true;
            if(first->slope(2,1) != second->slope(2,1)) failflag = true;
            if(first->slope(2,2) != second->slope(2,2)) failflag = true;
            if(first->slope(2,3) != second->slope(2,3)) failflag = true;
            if(first->slope(3,0) != second->slope(3,0)) failflag = true;
            if(first->slope(3,1) != second->slope(3,1)) failflag = true;
            if(first->slope(3,2) != second->slope(3,2)) failflag = true;
            if(first->slope(3,3) != second->slope(3,3)) failflag = true;
        if(failflag) throw cms::Exception("DataDoesNotMatch") << "Values are do not match";
            }
            listNewChan.erase(cell);  // fix 25.02.08
          }
      }
      std::cout << "These are identical" << std::endl;
    }



//  if(epsilon!=0){
    //implement compare qies -- different epsilon for slope and offset?
  //  }

   if(outfile.compare("null")!=0){
   for (std::vector<DetId>::iterator it = listRefChan.begin(); it != listRefChan.end(); it++)
      {
    DetId mydetid = *it;
    cell = std::find(listNewChan.begin(), listNewChan.end(), mydetid);
    if (cell == listNewChan.end()) // not present in new list
      {
        const HcalQIECoder* myCoder = myRefQIEs->getCoder( mydetid );
        std::cout << "o";
        resultQIEs->addCoder( *myCoder );
      }
    else // present in new list
      {
        const HcalQIECoder* myCoder = myNewQIEs->getCoder( mydetid );
        std::cout << "n";
        resultQIEs->addCoder( *myCoder );
        listNewChan.erase(cell);  // fix 25.02.08
      }
      }
    for (std::vector<DetId>::iterator it = listNewChan.begin(); it != listNewChan.end(); it++)  // fix 25.02.08
      {
    DetId mydetid = *it;
    const HcalQIECoder* myCoder = myNewQIEs->getCoder( mydetid );
    std::cout << "N";
    resultQIEs->addCoder( *myCoder );
      }

    std::ofstream outStream3(outfile.c_str());
    std::cout << "--- Dumping QIEs - the combined ones ---" << std::endl;
    resultQIEs->sort();
    HcalDbASCIIIO::dumpObject (outStream3, (*resultQIEs) );
    }

    std::cout << std::endl;
    if(checkemapflag){
    std::vector<DetId> listResult = resultQIEs->getAllChannels();
    // get the e-map list of channels
    std::vector<HcalGenericDetId> listEMap = myRefEMap->allPrecisionId();
    for (std::vector<HcalGenericDetId>::const_iterator it = listEMap.begin(); it != listEMap.end(); it++)
      {
      DetId mydetid = DetId(it->rawId());
    if (std::find(listResult.begin(), listResult.end(), mydetid ) == listResult.end()  )
      {
        std::cout << "Conditions not found for DetId = " << HcalGenericDetId(it->rawId()) << std::endl;
      }
      }
    }

}

Member Data Documentation

bool HcalQIEDataCheck::checkemapflag

private

Definition at line 51 of file HcalQIEDataCheck.h.

Referenced by analyze(), and HcalQIEDataCheck().

std::string HcalQIEDataCheck::dumprefs

private

Definition at line 49 of file HcalQIEDataCheck.h.

Referenced by analyze(), and HcalQIEDataCheck().

std::string HcalQIEDataCheck::dumpupdate

private

Definition at line 50 of file HcalQIEDataCheck.h.

Referenced by analyze(), and HcalQIEDataCheck().

std::string HcalQIEDataCheck::outfile

private

Definition at line 48 of file HcalQIEDataCheck.h.

Referenced by analyze(), and HcalQIEDataCheck().

bool HcalQIEDataCheck::validateflag

private

Definition at line 52 of file HcalQIEDataCheck.h.

Referenced by analyze(), and HcalQIEDataCheck().