CMS 3D CMS Logo

List of all members | Public Member Functions | Private Attributes
popcon::EcalPulseShapesHandler Class Reference

#include <EcalPulseShapesHandler.h>

Inheritance diagram for popcon::EcalPulseShapesHandler:
popcon::PopConSourceHandler< EcalPulseShapes >

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 EcalPulseShapesmypulseshapes
 

Additional Inherited Members

- Public Types inherited from popcon::PopConSourceHandler< EcalPulseShapes >
typedef std::map< Time_t, std::shared_ptr< EcalPulseShapes > > Container
 
typedef std::unique_ptr< EcalPulseShapesRef
 
typedef PopConSourceHandler< EcalPulseShapesself
 
typedef cond::Time_t Time_t
 
typedef EcalPulseShapes value_type
 
- Protected Member Functions inherited from popcon::PopConSourceHandler< EcalPulseShapes >
cond::persistency::SessiondbSession () 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.

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

◆ ~EcalPulseShapesHandler()

popcon::EcalPulseShapesHandler::~EcalPulseShapesHandler ( )
override

Definition at line 16 of file EcalPulseShapesHandler.cc.

16 {}

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.

18  {
19  // true means all is standard and OK
20  bool result = true;
21  for (int s = 0; s < EcalPulseShape::TEMPLATESAMPLES; ++s) {
22  if (item->pdfval[s] > 1 || item->pdfval[s] < 0)
23  result = false;
24  }
25  return result;
26 }
static const int TEMPLATESAMPLES

◆ 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.

28  {
29  for (int s = 0; s < EcalPulseShape::TEMPLATESAMPLES; ++s) {
30  item->pdfval[s] = isbarrel ? m_EBPulseShapeTemplate[s] : m_EEPulseShapeTemplate[s];
31  }
32 }
std::vector< double > m_EBPulseShapeTemplate
static const int TEMPLATESAMPLES
std::vector< double > m_EEPulseShapeTemplate

◆ 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, DetId::rawId(), alignCSCRings::s, AlCaHLTBitMon_QueryRunRegistry::string, EcalPulseShape::TEMPLATESAMPLES, EBDetId::validDetId(), and EEDetId::validDetId().

34  {
35  std::cout << "------- Ecal - > getNewObjects\n";
36 
37  // create the object pukse shapes
38  EcalPulseShapes* pulseshapes = new EcalPulseShapes();
39 
40  // read the templates from a text file
41  std::ifstream inputfile;
42  inputfile.open(m_filename.c_str());
43  float templatevals[EcalPulseShape::TEMPLATESAMPLES];
44  unsigned int rawId;
45  int isbarrel;
47 
48  // keep track of bad crystals
49  int nEBbad(0), nEEbad(0);
50  int nEBgood(0), nEEgood(0);
51  std::vector<EBDetId> ebgood;
52  std::vector<EEDetId> eegood;
53 
54  // fill with the measured shapes only for data
55  if (m_firstRun > 1) {
56  while (std::getline(inputfile, line)) {
57  std::istringstream linereader(line);
58  linereader >> isbarrel >> rawId;
59  // std::cout << "Inserting template for crystal with rawId = " << rawId << " (isbarrel = " << isbarrel << ") " << std::endl;
60  for (int s = 0; s < EcalPulseShape::TEMPLATESAMPLES; ++s) {
61  linereader >> templatevals[s];
62  // std::cout << templatevals[s] << "\t";
63  }
64  // std::cout << std::endl;
65 
66  if (!linereader) {
67  std::cout << "Wrong format of the text file. Exit." << std::endl;
68  return;
69  }
71  for (int s = 0; s < EcalPulseShape::TEMPLATESAMPLES; ++s)
72  item.pdfval[s] = templatevals[s];
73 
74  if (isbarrel) {
75  EBDetId ebdetid(rawId);
76  if (!checkPulseShape(&item))
77  nEBbad++;
78  else {
79  nEBgood++;
80  ebgood.push_back(ebdetid);
81  pulseshapes->insert(std::make_pair(ebdetid.rawId(), item));
82  }
83  } else {
84  EEDetId eedetid(rawId);
85  if (!checkPulseShape(&item))
86  nEEbad++;
87  else {
88  nEEgood++;
89  eegood.push_back(eedetid);
90  pulseshapes->insert(std::make_pair(eedetid.rawId(), item));
91  }
92  }
93  }
94  }
95 
96  // now fill the bad crystals and simulation with the simulation values (from TB)
97  std::cout << "Filled the DB with the good measured ECAL templates. Now filling the others with the TB values"
98  << std::endl;
99  for (int iEta = -EBDetId::MAX_IETA; iEta <= EBDetId::MAX_IETA; ++iEta) {
100  if (iEta == 0)
101  continue;
102  for (int iPhi = EBDetId::MIN_IPHI; iPhi <= EBDetId::MAX_IPHI; ++iPhi) {
103  if (EBDetId::validDetId(iEta, iPhi)) {
104  EBDetId ebdetid(iEta, iPhi, EBDetId::ETAPHIMODE);
105 
106  std::vector<EBDetId>::iterator it = find(ebgood.begin(), ebgood.end(), ebdetid);
107  if (it == ebgood.end()) {
109  fillSimPulseShape(&item, true);
110  pulseshapes->insert(std::make_pair(ebdetid.rawId(), item));
111  }
112  }
113  }
114  }
115 
116  for (int iZ = -1; iZ < 2; iZ += 2) {
117  for (int iX = EEDetId::IX_MIN; iX <= EEDetId::IX_MAX; ++iX) {
118  for (int iY = EEDetId::IY_MIN; iY <= EEDetId::IY_MAX; ++iY) {
119  if (EEDetId::validDetId(iX, iY, iZ)) {
120  EEDetId eedetid(iX, iY, iZ);
121 
122  std::vector<EEDetId>::iterator it = find(eegood.begin(), eegood.end(), eedetid);
123  if (it == eegood.end()) {
125  fillSimPulseShape(&item, false);
126  pulseshapes->insert(std::make_pair(eedetid.rawId(), item));
127  }
128  }
129  }
130  }
131  }
132 
133  std::cout << "Inserted the pulse shapes into the new item object" << std::endl;
134 
135  unsigned int irun = m_firstRun;
136  Time_t snc = (Time_t)irun;
137 
138  m_to_transfer.push_back(std::make_pair((EcalPulseShapes*)pulseshapes, snc));
139 
140  std::cout << "Ecal - > end of getNewObjects -----------" << std::endl;
141  std::cout << "N. bad shapes for EB = " << nEBbad << std::endl;
142  std::cout << "N. bad shapes for EE = " << nEEbad << std::endl;
143  std::cout << "Written the object" << std::endl;
144 }
EcalPulseShapesMap EcalPulseShapes
static const int MIN_IPHI
Definition: EBDetId.h:135
static const int TEMPLATESAMPLES
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
static const int IX_MIN
Definition: EEDetId.h:290
static const int IY_MIN
Definition: EEDetId.h:294
static bool validDetId(int i, int j)
check if a valid index combination
Definition: EBDetId.h:118
bool checkPulseShape(EcalPulseShapes::Item *item)
void fillSimPulseShape(EcalPulseShapes::Item *item, bool isbarrel)
std::vector< std::pair< EcalPulseShapes *, Time_t > > m_to_transfer
static const int ETAPHIMODE
Definition: EBDetId.h:158
void insert(std::pair< uint32_t, Item > const &a)
static const int IX_MAX
Definition: EEDetId.h:298
static const int MAX_IPHI
Definition: EBDetId.h:137
static bool validDetId(int crystal_ix, int crystal_iy, int iz)
Definition: EEDetId.h:248
static const int MAX_IETA
Definition: EBDetId.h:136
static const int IY_MAX
Definition: EEDetId.h:302

◆ id()

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

Implements popcon::PopConSourceHandler< EcalPulseShapes >.

Definition at line 47 of file EcalPulseShapesHandler.h.

References m_name.

47 { 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.