popcon::EcalPulseShapesHandler Class Reference

#include <EcalPulseShapesHandler.h>

Inheritance diagram for popcon::EcalPulseShapesHandler:

Public Member Functions
bool	checkPulseShape (EcalPulseShapes::Item *item)
	EcalPulseShapesHandler (edm::ParameterSet const &)
void	fillSimPulseShape (EcalPulseShapes::Item *item, bool isbarrel)
void	getNewObjects () override
std::string	id () const override
	~EcalPulseShapesHandler () override
Public Member Functions inherited from popcon::PopConSourceHandler< EcalPulseShapes >
void	initialize (const cond::persistency::Session &dbSession, cond::TagInfo_t const &tagInfo, cond::LogDBEntry_t const &logDBEntry)
Ref	lastPayload () const
cond::LogDBEntry_t const &	logDBEntry () const
std::pair< Container const *, std::string const >	operator() (const cond::persistency::Session &session, cond::TagInfo_t const &tagInfo, cond::LogDBEntry_t const &logDBEntry) const
	PopConSourceHandler ()
Container const &	returnData ()
cond::TagInfo_t const &	tagInfo () const
std::string const &	userTextLog () const
virtual	~PopConSourceHandler ()

Private Attributes
std::vector< double >	m_EBPulseShapeTemplate
std::vector< double >	m_EEPulseShapeTemplate
std::string	m_filename
unsigned int	m_firstRun
std::string	m_gentag
unsigned int	m_lastRun
std::string	m_name
const EcalPulseShapes *	mypulseshapes

Additional Inherited Members
Public Types inherited from popcon::PopConSourceHandler< EcalPulseShapes >
typedef std::map< Time_t, std::shared_ptr< EcalPulseShapes > >	Container
typedef std::unique_ptr< EcalPulseShapes >	Ref
typedef PopConSourceHandler< EcalPulseShapes >	self
typedef cond::Time_t	Time_t
typedef EcalPulseShapes	value_type
Protected Member Functions inherited from popcon::PopConSourceHandler< EcalPulseShapes >
cond::persistency::Session &	dbSession () const
Protected Attributes inherited from popcon::PopConSourceHandler< EcalPulseShapes >
Container	m_iovs
std::vector< std::pair< EcalPulseShapes *, Time_t > >	m_to_transfer
std::string	m_userTextLog

Detailed Description

Definition at line 40 of file EcalPulseShapesHandler.h.

Constructor & Destructor Documentation

EcalPulseShapesHandler()

popcon::EcalPulseShapesHandler::EcalPulseShapesHandler

(

edm::ParameterSet const &

ps

)

Definition at line 7 of file EcalPulseShapesHandler.cc.

References gather_cfg::cout, edm::ParameterSet::getParameter(), m_EBPulseShapeTemplate, m_EEPulseShapeTemplate, m_filename, m_firstRun, and AlCaHLTBitMon_QueryRunRegistry::string.

    : m_name(ps.getUntrackedParameter<std::string>("name", "EcalPulseShapesHandler")) {
  std::cout << "EcalPulseShape Source handler constructor\n" << std::endl;
  m_firstRun = static_cast<unsigned int>(atoi(ps.getParameter<std::string>("firstRun").c_str()));
  m_filename = ps.getParameter<std::string>("inputFileName");
  m_EBPulseShapeTemplate = ps.getParameter<std::vector<double> >("EBPulseShapeTemplate");
  m_EEPulseShapeTemplate = ps.getParameter<std::vector<double> >("EEPulseShapeTemplate");
}

~EcalPulseShapesHandler()

popcon::EcalPulseShapesHandler::~EcalPulseShapesHandler ( )

override

Definition at line 16 of file EcalPulseShapesHandler.cc.

{}

Member Function Documentation

checkPulseShape()

bool popcon::EcalPulseShapesHandler::checkPulseShape

(

EcalPulseShapes::Item *

item

)

Definition at line 18 of file EcalPulseShapesHandler.cc.

References B2GTnPMonitor_cfi::item, mps_fire::result, alignCSCRings::s, and EcalPulseShape::TEMPLATESAMPLES.

                                                                            {
  // true means all is standard and OK
  bool result = true;
  for (int s = 0; s < EcalPulseShape::TEMPLATESAMPLES; ++s) {
    if (item->pdfval[s] > 1 || item->pdfval[s] < 0)
      result = false;
  }
  return result;
}

fillSimPulseShape()

void popcon::EcalPulseShapesHandler::fillSimPulseShape	(	EcalPulseShapes::Item *	item,
		bool	isbarrel
	)

Definition at line 28 of file EcalPulseShapesHandler.cc.

References B2GTnPMonitor_cfi::item, alignCSCRings::s, and EcalPulseShape::TEMPLATESAMPLES.

                                                                                             {
  for (int s = 0; s < EcalPulseShape::TEMPLATESAMPLES; ++s) {
    item->pdfval[s] = isbarrel ? m_EBPulseShapeTemplate[s] : m_EEPulseShapeTemplate[s];
  }
}

getNewObjects()

void popcon::EcalPulseShapesHandler::getNewObjects ( )

overridevirtual

Implements popcon::PopConSourceHandler< EcalPulseShapes >.

Definition at line 34 of file EcalPulseShapesHandler.cc.

References gather_cfg::cout, EBDetId::ETAPHIMODE, spr::find(), l1tTowerCalibrationProducer_cfi::iEta, caHitNtupletGeneratorKernels::if(), EcalCondObjectContainer< T >::insert(), B2GTnPMonitor_cfi::item, EEDetId::IX_MAX, EEDetId::IX_MIN, EEDetId::IY_MAX, EEDetId::IY_MIN, mps_splice::line, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, nano_mu_digi_cff::rawId, DetId::rawId(), alignCSCRings::s, AlCaHLTBitMon_QueryRunRegistry::string, EcalPulseShape::TEMPLATESAMPLES, EBDetId::validDetId(), and EEDetId::validDetId().

                                                 {
  std::cout << "------- Ecal - > getNewObjects\n";

  // create the object pukse shapes
  EcalPulseShapes* pulseshapes = new EcalPulseShapes();

  // read the templates from a text file
  std::ifstream inputfile;
  inputfile.open(m_filename.c_str());
  float templatevals[EcalPulseShape::TEMPLATESAMPLES];
  unsigned int rawId;
  int isbarrel;
  std::string line;

  // keep track of bad crystals
  int nEBbad(0), nEEbad(0);
  std::vector<EBDetId> ebgood;
  std::vector<EEDetId> eegood;

  // fill with the measured shapes only for data
  if (m_firstRun > 1) {
    while (std::getline(inputfile, line)) {
      std::istringstream linereader(line);
      linereader >> isbarrel >> rawId;
      // std::cout << "Inserting template for crystal with rawId = " << rawId << " (isbarrel = " << isbarrel << ") " << std::endl;
      for (int s = 0; s < EcalPulseShape::TEMPLATESAMPLES; ++s) {
        linereader >> templatevals[s];
        // std::cout << templatevals[s] << "\t";
      }
      // std::cout << std::endl;

      if (!linereader) {
        std::cout << "Wrong format of the text file. Exit." << std::endl;
        return;
      }
      EcalPulseShapes::Item item;
      for (int s = 0; s < EcalPulseShape::TEMPLATESAMPLES; ++s)
        item.pdfval[s] = templatevals[s];

      if (isbarrel) {
        EBDetId ebdetid(rawId);
        if (!checkPulseShape(&item))
          nEBbad++;
        else {
          ebgood.push_back(ebdetid);
          pulseshapes->insert(std::make_pair(ebdetid.rawId(), item));
        }
      } else {
        EEDetId eedetid(rawId);
        if (!checkPulseShape(&item))
          nEEbad++;
        else {
          eegood.push_back(eedetid);
          pulseshapes->insert(std::make_pair(eedetid.rawId(), item));
        }
      }
    }
  }

  // now fill the bad crystals and simulation with the simulation values (from TB)
  std::cout << "Filled the DB with the good measured ECAL templates. Now filling the others with the TB values"
            << std::endl;
  for (int iEta = -EBDetId::MAX_IETA; iEta <= EBDetId::MAX_IETA; ++iEta) {
    if (iEta == 0)
      continue;
    for (int iPhi = EBDetId::MIN_IPHI; iPhi <= EBDetId::MAX_IPHI; ++iPhi) {
      if (EBDetId::validDetId(iEta, iPhi)) {
        EBDetId ebdetid(iEta, iPhi, EBDetId::ETAPHIMODE);

        std::vector<EBDetId>::iterator it = find(ebgood.begin(), ebgood.end(), ebdetid);
        if (it == ebgood.end()) {
          EcalPulseShapes::Item item;
          fillSimPulseShape(&item, true);
          pulseshapes->insert(std::make_pair(ebdetid.rawId(), item));
        }
      }
    }
  }

  for (int iZ = -1; iZ < 2; iZ += 2) {
    for (int iX = EEDetId::IX_MIN; iX <= EEDetId::IX_MAX; ++iX) {
      for (int iY = EEDetId::IY_MIN; iY <= EEDetId::IY_MAX; ++iY) {
        if (EEDetId::validDetId(iX, iY, iZ)) {
          EEDetId eedetid(iX, iY, iZ);

          std::vector<EEDetId>::iterator it = find(eegood.begin(), eegood.end(), eedetid);
          if (it == eegood.end()) {
            EcalPulseShapes::Item item;
            fillSimPulseShape(&item, false);
            pulseshapes->insert(std::make_pair(eedetid.rawId(), item));
          }
        }
      }
    }
  }

  std::cout << "Inserted the pulse shapes into the new item object" << std::endl;

  unsigned int irun = m_firstRun;
  Time_t snc = (Time_t)irun;

  m_to_transfer.push_back(std::make_pair((EcalPulseShapes*)pulseshapes, snc));

  std::cout << "Ecal - > end of getNewObjects -----------" << std::endl;
  std::cout << "N. bad shapes for EB = " << nEBbad << std::endl;
  std::cout << "N. bad shapes for EE = " << nEEbad << std::endl;
  std::cout << "Written the object" << std::endl;
}

id()

std::string popcon::EcalPulseShapesHandler::id ( ) const

inlineoverridevirtual

Implements popcon::PopConSourceHandler< EcalPulseShapes >.

Definition at line 47 of file EcalPulseShapesHandler.h.

References m_name.

{ return m_name; }

Member Data Documentation

m_EBPulseShapeTemplate

std::vector<double> popcon::EcalPulseShapesHandler::m_EBPulseShapeTemplate

private

Definition at line 58 of file EcalPulseShapesHandler.h.

Referenced by EcalPulseShapesHandler().

m_EEPulseShapeTemplate

std::vector<double> popcon::EcalPulseShapesHandler::m_EEPulseShapeTemplate

private

Definition at line 58 of file EcalPulseShapesHandler.h.

Referenced by EcalPulseShapesHandler().

m_filename

std::string popcon::EcalPulseShapesHandler::m_filename

private

Definition at line 56 of file EcalPulseShapesHandler.h.

Referenced by EcalPulseShapesHandler().

m_firstRun

unsigned int popcon::EcalPulseShapesHandler::m_firstRun

private

Definition at line 52 of file EcalPulseShapesHandler.h.

Referenced by EcalPulseShapesHandler().

m_gentag

std::string popcon::EcalPulseShapesHandler::m_gentag

private

Definition at line 55 of file EcalPulseShapesHandler.h.

m_lastRun

unsigned int popcon::EcalPulseShapesHandler::m_lastRun

private

Definition at line 53 of file EcalPulseShapesHandler.h.

m_name

std::string popcon::EcalPulseShapesHandler::m_name

private

Definition at line 57 of file EcalPulseShapesHandler.h.

Referenced by id().

mypulseshapes

const EcalPulseShapes* popcon::EcalPulseShapesHandler::mypulseshapes

private

Definition at line 50 of file EcalPulseShapesHandler.h.