|
|
|
#include <DataMixingHcalDigiWorker.h>
Definition at line 38 of file DataMixingHcalDigiWorker.h.
typedef std::multimap<DetId, CaloSamples> edm::DataMixingHcalDigiWorker::HBHEDigiMap [private] |
Definition at line 74 of file DataMixingHcalDigiWorker.h.
typedef std::multimap<DetId, CaloSamples> edm::DataMixingHcalDigiWorker::HFDigiMap [private] |
Definition at line 75 of file DataMixingHcalDigiWorker.h.
typedef std::multimap<DetId, CaloSamples> edm::DataMixingHcalDigiWorker::HODigiMap [private] |
Definition at line 76 of file DataMixingHcalDigiWorker.h.
typedef std::multimap<DetId, CaloSamples> edm::DataMixingHcalDigiWorker::ZDCDigiMap [private] |
Definition at line 77 of file DataMixingHcalDigiWorker.h.
| DataMixingHcalDigiWorker::DataMixingHcalDigiWorker | ( | ) |
Definition at line 34 of file DataMixingHcalDigiWorker.cc.
{ }
| DataMixingHcalDigiWorker::DataMixingHcalDigiWorker | ( | const edm::ParameterSet & | ps | ) | [explicit] |
standard constructor
Definition at line 37 of file DataMixingHcalDigiWorker.cc.
References DoZDC_, edm::ParameterSet::getParameter(), HBHEDigiCollectionDM_, HBHEdigiCollectionSig_, HBHEPileInputTag_, HFDigiCollectionDM_, HFdigiCollectionSig_, HFPileInputTag_, HODigiCollectionDM_, HOdigiCollectionSig_, HOPileInputTag_, edm::InputTag::label(), AlCaHLTBitMon_QueryRunRegistry::string, ZDCDigiCollectionDM_, ZDCdigiCollectionSig_, and ZDCPileInputTag_.
:
label_(ps.getParameter<std::string>("Label"))
{
// get the subdetector names
// this->getSubdetectorNames(); //something like this may be useful to check what we are supposed to do...
// declare the products to produce
// Hcal
HBHEdigiCollectionSig_ = ps.getParameter<edm::InputTag>("HBHEdigiCollectionSig");
HOdigiCollectionSig_ = ps.getParameter<edm::InputTag>("HOdigiCollectionSig");
HFdigiCollectionSig_ = ps.getParameter<edm::InputTag>("HFdigiCollectionSig");
ZDCdigiCollectionSig_ = ps.getParameter<edm::InputTag>("ZDCdigiCollectionSig");
HBHEPileInputTag_ = ps.getParameter<edm::InputTag>("HBHEPileInputTag");
HOPileInputTag_ = ps.getParameter<edm::InputTag>("HOPileInputTag");
HFPileInputTag_ = ps.getParameter<edm::InputTag>("HFPileInputTag");
ZDCPileInputTag_ = ps.getParameter<edm::InputTag>("ZDCPileInputTag");
DoZDC_ = false;
if(ZDCPileInputTag_.label() != "") DoZDC_ = true;
HBHEDigiCollectionDM_ = ps.getParameter<std::string>("HBHEDigiCollectionDM");
HODigiCollectionDM_ = ps.getParameter<std::string>("HODigiCollectionDM");
HFDigiCollectionDM_ = ps.getParameter<std::string>("HFDigiCollectionDM");
ZDCDigiCollectionDM_ = ps.getParameter<std::string>("ZDCDigiCollectionDM");
}
| DataMixingHcalDigiWorker::~DataMixingHcalDigiWorker | ( | ) | [virtual] |
| void DataMixingHcalDigiWorker::addHcalPileups | ( | const int | bcr, |
| const edm::EventPrincipal * | ep, | ||
| unsigned int | EventId, | ||
| const edm::EventSetup & | ES | ||
| ) |
Definition at line 263 of file DataMixingHcalDigiWorker.cc.
References HcalCoderDb::adc2fC(), edm::SortedCollection< T, SORT >::begin(), DoZDC_, edm::SortedCollection< T, SORT >::end(), edm::EventSetup::get(), HBHEDigiStorage_, HBHEPileInputTag_, HFDigiStorage_, HFPileInputTag_, HODigiStorage_, HOPileInputTag_, edm::EventPrincipal::id(), LogDebug, edm::SortedCollection< T, SORT >::size(), ZDCDigiStorage_, and ZDCPileInputTag_.
Referenced by edm::DataMixingModule::pileWorker().
{
LogDebug("DataMixingHcalDigiWorker") <<"\n===============> adding pileups from event "<<ep->id()<<" for bunchcrossing "<<bcr;
// get conditions
edm::ESHandle<HcalDbService> conditions;
ES.get<HcalDbRecord>().get(conditions);
// fill in maps of hits; same code as addSignals, except now applied to the pileup events
// HBHE first
boost::shared_ptr<Wrapper<HBHEDigiCollection> const> HBHEDigisPTR =
getProductByTag<HBHEDigiCollection>(*ep, HBHEPileInputTag_ );
if(HBHEDigisPTR ) {
const HBHEDigiCollection* HBHEDigis = const_cast< HBHEDigiCollection * >(HBHEDigisPTR->product());
LogInfo("DataMixingHcalDigiWorker") << "total # HBHE digis: " << HBHEDigis->size();
// loop over digis, adding these to the existing maps
for(HBHEDigiCollection::const_iterator it = HBHEDigis->begin();
it != HBHEDigis->end(); ++it) {
// calibration, for future reference: (same block for all Hcal types)
HcalDetId cell = it->id();
// const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); // this one is generic
HcalCoderDb coder (*channelCoder, *shape);
CaloSamples tool;
coder.adc2fC((*it),tool);
HBHEDigiStorage_.insert(HBHEDigiMap::value_type( (it->id()), tool ));
#ifdef DEBUG
LogDebug("DataMixingHcalDigiWorker") << "processed HBHEDigi with rawId: "
<< it->id() << "\n"
<< " digi energy: " << it->energy();
#endif
}
}
// HO Next
boost::shared_ptr<Wrapper<HODigiCollection> const> HODigisPTR =
getProductByTag<HODigiCollection>(*ep, HOPileInputTag_ );
if(HODigisPTR ) {
const HODigiCollection* HODigis = const_cast< HODigiCollection * >(HODigisPTR->product());
LogDebug("DataMixingHcalDigiWorker") << "total # HO digis: " << HODigis->size();
// loop over digis, adding these to the existing maps
for(HODigiCollection::const_iterator it = HODigis->begin();
it != HODigis->end(); ++it) {
// calibration, for future reference: (same block for all Hcal types)
HcalDetId cell = it->id();
// const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); // this one is generic
HcalCoderDb coder (*channelCoder, *shape);
CaloSamples tool;
coder.adc2fC((*it),tool);
HODigiStorage_.insert(HODigiMap::value_type( (it->id()), tool ));
#ifdef DEBUG
LogDebug("DataMixingHcalDigiWorker") << "processed HODigi with rawId: "
<< it->id() << "\n"
<< " digi energy: " << it->energy();
#endif
}
}
// HF Next
boost::shared_ptr<Wrapper<HFDigiCollection> const> HFDigisPTR =
getProductByTag<HFDigiCollection>(*ep, HFPileInputTag_ );
if(HFDigisPTR ) {
const HFDigiCollection* HFDigis = const_cast< HFDigiCollection * >(HFDigisPTR->product());
LogDebug("DataMixingHcalDigiWorker") << "total # HF digis: " << HFDigis->size();
// loop over digis, adding these to the existing maps
for(HFDigiCollection::const_iterator it = HFDigis->begin();
it != HFDigis->end(); ++it) {
// calibration, for future reference: (same block for all Hcal types)
HcalDetId cell = it->id();
// const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); // this one is generic
HcalCoderDb coder (*channelCoder, *shape);
CaloSamples tool;
coder.adc2fC((*it),tool);
HFDigiStorage_.insert(HFDigiMap::value_type( (it->id()), tool ));
#ifdef DEBUG
LogDebug("DataMixingHcalDigiWorker") << "processed HFDigi with rawId: "
<< it->id() << "\n"
<< " digi energy: " << it->energy();
#endif
}
}
// ZDC Next
if(DoZDC_) {
boost::shared_ptr<Wrapper<ZDCDigiCollection> const> ZDCDigisPTR =
getProductByTag<ZDCDigiCollection>(*ep, ZDCPileInputTag_ );
if(ZDCDigisPTR ) {
const ZDCDigiCollection* ZDCDigis = const_cast< ZDCDigiCollection * >(ZDCDigisPTR->product());
LogDebug("DataMixingHcalDigiWorker") << "total # ZDC digis: " << ZDCDigis->size();
// loop over digis, adding these to the existing maps
for(ZDCDigiCollection::const_iterator it = ZDCDigis->begin();
it != ZDCDigis->end(); ++it) {
// calibration, for future reference: (same block for all Hcal types)
HcalDetId cell = it->id();
// const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); // this one is generic
HcalCoderDb coder (*channelCoder, *shape);
CaloSamples tool;
coder.adc2fC((*it),tool);
ZDCDigiStorage_.insert(ZDCDigiMap::value_type( (it->id()), tool ));
#ifdef DEBUG
LogDebug("DataMixingHcalDigiWorker") << "processed ZDCDigi with rawId: "
<< it->id() << "\n"
<< " digi energy: " << it->energy();
#endif
}
}
}
}
| void DataMixingHcalDigiWorker::addHcalSignals | ( | const edm::Event & | e, |
| const edm::EventSetup & | ES | ||
| ) |
Definition at line 74 of file DataMixingHcalDigiWorker.cc.
References HcalCoderDb::adc2fC(), edm::SortedCollection< T, SORT >::begin(), DoZDC_, edm::SortedCollection< T, SORT >::end(), edm::EventSetup::get(), edm::Event::getByLabel(), HBHEdigiCollectionSig_, HBHEDigiStorage_, HFdigiCollectionSig_, HFDigiStorage_, HOdigiCollectionSig_, HODigiStorage_, edm::EventBase::id(), LogDebug, edm::Handle< T >::product(), edm::SortedCollection< T, SORT >::size(), ZDCdigiCollectionSig_, and ZDCDigiStorage_.
Referenced by edm::DataMixingModule::addSignals().
{
// Calibration stuff will look like this:
// get conditions
edm::ESHandle<HcalDbService> conditions;
ES.get<HcalDbRecord>().get(conditions);
// fill in maps of hits
LogInfo("DataMixingHcalDigiWorker")<<"===============> adding MC signals for "<<e.id();
// HBHE first
Handle< HBHEDigiCollection > pHBHEDigis;
const HBHEDigiCollection* HBHEDigis = 0;
if( e.getByLabel( HBHEdigiCollectionSig_, pHBHEDigis) ) {
HBHEDigis = pHBHEDigis.product(); // get a ptr to the product
LogDebug("DataMixingHcalDigiWorker") << "total # HBHE digis: " << HBHEDigis->size();
}
// else { std::cout << "NO HBHE Digis " << HBHEdigiCollectionSig_.label() << std::endl;}
if (HBHEDigis)
{
// loop over digis, storing them in a map so we can add pileup later
for(HBHEDigiCollection::const_iterator it = HBHEDigis->begin();
it != HBHEDigis->end(); ++it) {
// calibration, for future reference: (same block for all Hcal types)
HcalDetId cell = it->id();
// const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); // this one is generic
HcalCoderDb coder (*channelCoder, *shape);
CaloSamples tool;
coder.adc2fC((*it),tool);
// put sample values back into digi?
// RecHit MyHit = reco_.reconstruct(*it,coder,calibrations));
//... can now fish calibrated information from RecHit
HBHEDigiStorage_.insert(HBHEDigiMap::value_type( ( it->id() ), tool ));
#ifdef DEBUG
LogDebug("DataMixingHcalDigiWorker") << "processed HBHEDigi with rawId: "
<< it->id() << "\n"
<< " digi energy: " << it->energy();
#endif
}
}
// HO next
Handle< HODigiCollection > pHODigis;
const HODigiCollection* HODigis = 0;
if( e.getByLabel( HOdigiCollectionSig_, pHODigis) ){
HODigis = pHODigis.product(); // get a ptr to the product
#ifdef DEBUG
LogDebug("DataMixingHcalDigiWorker") << "total # HO digis: " << HODigis->size();
#endif
}
if (HODigis)
{
// loop over digis, storing them in a map so we can add pileup later
for(HODigiCollection::const_iterator it = HODigis->begin();
it != HODigis->end(); ++it) {
// calibration, for future reference: (same block for all Hcal types)
HcalDetId cell = it->id();
// const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); // this one is generic
HcalCoderDb coder (*channelCoder, *shape);
CaloSamples tool;
coder.adc2fC((*it),tool);
HODigiStorage_.insert(HODigiMap::value_type( ( it->id() ), tool ));
#ifdef DEBUG
LogDebug("DataMixingHcalDigiWorker") << "processed HODigi with rawId: "
<< it->id() << "\n"
<< " digi energy: " << it->energy();
#endif
}
}
// HF next
Handle< HFDigiCollection > pHFDigis;
const HFDigiCollection* HFDigis = 0;
if( e.getByLabel( HFdigiCollectionSig_, pHFDigis) ) {
HFDigis = pHFDigis.product(); // get a ptr to the product
#ifdef DEBUG
LogDebug("DataMixingHcalDigiWorker") << "total # HF digis: " << HFDigis->size();
#endif
}
if (HFDigis)
{
// loop over digis, storing them in a map so we can add pileup later
for(HFDigiCollection::const_iterator it = HFDigis->begin();
it != HFDigis->end(); ++it) {
// calibration, for future reference: (same block for all Hcal types)
HcalDetId cell = it->id();
// const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); // this one is generic
HcalCoderDb coder (*channelCoder, *shape);
CaloSamples tool;
coder.adc2fC((*it),tool);
HFDigiStorage_.insert(HFDigiMap::value_type( ( it->id() ), tool ));
#ifdef DEBUG
LogDebug("DataMixingHcalDigiWorker") << "processed HFDigi with rawId: "
<< it->id() << "\n"
<< " digi energy: " << it->energy();
#endif
}
}
// ZDC next
if(DoZDC_){
Handle< ZDCDigiCollection > pZDCDigis;
const ZDCDigiCollection* ZDCDigis = 0;
if( e.getByLabel( ZDCdigiCollectionSig_, pZDCDigis) ) {
ZDCDigis = pZDCDigis.product(); // get a ptr to the product
#ifdef DEBUG
LogDebug("DataMixingHcalDigiWorker") << "total # ZDC digis: " << ZDCDigis->size();
#endif
}
if (ZDCDigis)
{
// loop over digis, storing them in a map so we can add pileup later
for(ZDCDigiCollection::const_iterator it = ZDCDigis->begin();
it != ZDCDigis->end(); ++it) {
// calibration, for future reference: (same block for all Hcal types)
HcalDetId cell = it->id();
// const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); // this one is generic
HcalCoderDb coder (*channelCoder, *shape);
CaloSamples tool;
coder.adc2fC((*it),tool);
ZDCDigiStorage_.insert(ZDCDigiMap::value_type( ( it->id() ), tool ));
#ifdef DEBUG
LogDebug("DataMixingHcalDigiWorker") << "processed ZDCDigi with rawId: "
<< it->id() << "\n"
<< " digi energy: " << it->energy();
#endif
}
}
}
} // end of addHCalSignals
| void DataMixingHcalDigiWorker::putHcal | ( | edm::Event & | e, |
| const edm::EventSetup & | ES | ||
| ) |
Definition at line 421 of file DataMixingHcalDigiWorker.cc.
References HcalCoderDb::fC2adc(), edm::EventSetup::get(), HBHEDigiCollectionDM_, HBHEDigiStorage_, HFDigiCollectionDM_, HFDigiStorage_, HODigiCollectionDM_, HODigiStorage_, CaloSamples::id(), edm::max(), edm::Event::put(), findQualityFiles::size, CaloSamples::size(), ZDCDigiCollectionDM_, and ZDCDigiStorage_.
Referenced by edm::DataMixingModule::put().
{
// collection of digis to put in the event
std::auto_ptr< HBHEDigiCollection > HBHEdigis( new HBHEDigiCollection );
std::auto_ptr< HODigiCollection > HOdigis( new HODigiCollection );
std::auto_ptr< HFDigiCollection > HFdigis( new HFDigiCollection );
std::auto_ptr< ZDCDigiCollection > ZDCdigis( new ZDCDigiCollection );
// get conditions
edm::ESHandle<HcalDbService> conditions;
ES.get<HcalDbRecord>().get(conditions);
// loop over the maps we have, re-making individual hits or digis if necessary.
DetId formerID = 0;
DetId currentID;
CaloSamples HB_old;
double fC_new;
double fC_old;
double fC_sum;
// HB first...
HBHEDigiMap::const_iterator iHBchk;
for(HBHEDigiMap::const_iterator iHB = HBHEDigiStorage_.begin();
iHB != HBHEDigiStorage_.end(); ++iHB) {
currentID = iHB->first;
if (currentID == formerID) { // we have to add these digis together
//loop over digi samples in each CaloSample
unsigned int sizenew = (iHB->second).size();
unsigned int sizeold = HB_old.size();
bool usenew = false;
if(sizenew > sizeold) usenew = true;
unsigned int max_samp = std::max(sizenew, sizeold);
CaloSamples HB_bigger(currentID,max_samp);
// HB_old.setSize(max_samp); --> can't do this...
// samples from different events can be of different lengths - sum all
// that overlap.
for(unsigned int isamp = 0; isamp<max_samp; isamp++) {
if(isamp < sizenew) {
fC_new = (iHB->second)[isamp]; // should return nominal_fC();
}
else { fC_new = 0;}
if(isamp < sizeold) {
fC_old = HB_old[isamp];
}
else { fC_old = 0;}
// add values
fC_sum = fC_new + fC_old;
//uint fCS = int(fC_sum);
//const HcalQIESample fC(fCS);
//HB_old.setSample(isamp, fC);
if(usenew) {HB_bigger[isamp] = fC_sum; }
else { HB_old[isamp] = fC_sum; } // overwrite old sample, adding new info
}
if(usenew) HB_old = HB_bigger; // save new, larger sized sample in "old" slot
}
else {
if(formerID>0) {
// make new digi
HBHEdigis->push_back(HBHEDataFrame(formerID));
// set up information to convert back
HcalDetId cell = HB_old.id();
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); // this one is generic
HcalCoderDb coder (*channelCoder, *shape);
unsigned int sizeold = HB_old.size();
for(unsigned int isamp = 0; isamp<sizeold; isamp++) {
coder.fC2adc(HB_old,(HBHEdigis->back()), 0 ); // as per simulation, capid=0???
}
}
//save pointers for next iteration
formerID = currentID;
HB_old = iHB->second;
//OldUpAdd = HB_old.id();
}
iHBchk = iHB;
if((++iHBchk) == HBHEDigiStorage_.end()) { //make sure not to lose the last one
// make new digi
HBHEdigis->push_back(HBHEDataFrame(currentID));
// set up information to convert back
HcalDetId cell = (iHB->second).id();
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); // this one is generic
HcalCoderDb coder (*channelCoder, *shape);
unsigned int sizenew = (iHB->second).size();
for(unsigned int isamp = 0; isamp<sizenew; isamp++) {
coder.fC2adc(HB_old,(HBHEdigis->back()), 0 ); // as per simulation, capid=0???
}
}
}
// HO next...
// loop over the maps we have, re-making individual hits or digis if necessary.
formerID = 0;
CaloSamples HO_old;
HODigiMap::const_iterator iHOchk;
for(HODigiMap::const_iterator iHO = HODigiStorage_.begin();
iHO != HODigiStorage_.end(); ++iHO) {
currentID = iHO->first;
if (currentID == formerID) { // we have to add these digis together
//loop over digi samples in each CaloSample
unsigned int sizenew = (iHO->second).size();
unsigned int sizeold = HO_old.size();
unsigned int max_samp = std::max(sizenew, sizeold);
CaloSamples HO_bigger(currentID,max_samp);
bool usenew = false;
if(sizenew > sizeold) usenew = true;
// samples from different events can be of different lengths - sum all
// that overlap.
for(unsigned int isamp = 0; isamp<max_samp; isamp++) {
if(isamp < sizenew) {
fC_new = (iHO->second)[isamp];
}
else { fC_new = 0;}
if(isamp < sizeold) {
fC_old = HO_old[isamp];
}
else { fC_old = 0;}
// add values
fC_sum = fC_new + fC_old;
if(usenew) {HO_bigger[isamp] = fC_sum; }
else { HO_old[isamp] = fC_sum; } // overwrite old sample, adding new info
}
if(usenew) HO_old = HO_bigger; // save new, larger sized sample in "old" slot
}
else {
if(formerID>0) {
// make new digi
HOdigis->push_back(HODataFrame(formerID));
// set up information to convert back
HcalDetId cell = HO_old.id();
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); // this one is generic
HcalCoderDb coder (*channelCoder, *shape);
unsigned int sizeold = HO_old.size();
for(unsigned int isamp = 0; isamp<sizeold; isamp++) {
coder.fC2adc(HO_old,(HOdigis->back()), 0 ); // as per simulation, capid=0???
}
}
//save pointers for next iteration
formerID = currentID;
HO_old = iHO->second;
}
iHOchk = iHO;
if((++iHOchk) == HODigiStorage_.end()) { //make sure not to lose the last one
// make new digi
HOdigis->push_back(HODataFrame(currentID));
// set up information to convert back
HcalDetId cell = (iHO->second).id();
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); // this one is generic
HcalCoderDb coder (*channelCoder, *shape);
unsigned int sizeold = (iHO->second).size();
for(unsigned int isamp = 0; isamp<sizeold; isamp++) {
coder.fC2adc(HO_old,(HOdigis->back()), 0 ); // as per simulation, capid=0???
}
}
}
// HF next...
// loop over the maps we have, re-making individual hits or digis if necessary.
formerID = 0;
CaloSamples HF_old;
HFDigiMap::const_iterator iHFchk;
for(HFDigiMap::const_iterator iHF = HFDigiStorage_.begin();
iHF != HFDigiStorage_.end(); ++iHF) {
currentID = iHF->first;
if (currentID == formerID) { // we have to add these digis together
//loop over digi samples in each CaloSample
unsigned int sizenew = (iHF->second).size();
unsigned int sizeold = HF_old.size();
unsigned int max_samp = std::max(sizenew, sizeold);
CaloSamples HF_bigger(currentID,max_samp);
bool usenew = false;
if(sizenew > sizeold) usenew = true;
// samples from different events can be of different lengths - sum all
// that overlap.
for(unsigned int isamp = 0; isamp<max_samp; isamp++) {
if(isamp < sizenew) {
fC_new = (iHF->second)[isamp];
}
else { fC_new = 0;}
if(isamp < sizeold) {
fC_old = HF_old[isamp];
}
else { fC_old = 0;}
// add values
fC_sum = fC_new + fC_old;
if(usenew) {HF_bigger[isamp] = fC_sum; }
else { HF_old[isamp] = fC_sum; } // overwrite old sample, adding new info
}
if(usenew) HF_old = HF_bigger; // save new, larger sized sample in "old" slot
}
else {
if(formerID>0) {
// make new digi
HFdigis->push_back(HFDataFrame(formerID));
// set up information to convert back
HcalDetId cell = HF_old.id();
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); // this one is generic
HcalCoderDb coder (*channelCoder, *shape);
unsigned int sizeold = HF_old.size();
for(unsigned int isamp = 0; isamp<sizeold; isamp++) {
coder.fC2adc(HF_old,(HFdigis->back()), 0 ); // as per simulation, capid=0???
}
}
//save pointers for next iteration
formerID = currentID;
HF_old = iHF->second;
}
iHFchk = iHF;
if((++iHFchk) == HFDigiStorage_.end()) { //make sure not to lose the last one
// make new digi
HFdigis->push_back(HFDataFrame(currentID));
// set up information to convert back
HcalDetId cell = (iHF->second).id();
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); // this one is generic
HcalCoderDb coder (*channelCoder, *shape);
unsigned int sizeold = (iHF->second).size();
for(unsigned int isamp = 0; isamp<sizeold; isamp++) {
coder.fC2adc(HF_old,(HFdigis->back()), 0 ); // as per simulation, capid=0???
}
}
}
// ZDC next...
// loop over the maps we have, re-making individual hits or digis if necessary.
formerID = 0;
CaloSamples ZDC_old;
ZDCDigiMap::const_iterator iZDCchk;
for(ZDCDigiMap::const_iterator iZDC = ZDCDigiStorage_.begin();
iZDC != ZDCDigiStorage_.end(); ++iZDC) {
currentID = iZDC->first;
if (currentID == formerID) { // we have to add these digis together
//loop over digi samples in each CaloSample
unsigned int sizenew = (iZDC->second).size();
unsigned int sizeold = ZDC_old.size();
unsigned int max_samp = std::max(sizenew, sizeold);
CaloSamples ZDC_bigger(currentID,max_samp);
bool usenew = false;
if(sizenew > sizeold) usenew = true;
// samples from different events can be of different lengths - sum all
// that overlap.
for(unsigned int isamp = 0; isamp<max_samp; isamp++) {
if(isamp < sizenew) {
fC_new = (iZDC->second)[isamp];
}
else { fC_new = 0;}
if(isamp < sizeold) {
fC_old = ZDC_old[isamp];
}
else { fC_old = 0;}
// add values
fC_sum = fC_new + fC_old;
if(usenew) {ZDC_bigger[isamp] = fC_sum; }
else { ZDC_old[isamp] = fC_sum; } // overwrite old sample, adding new info
}
if(usenew) ZDC_old = ZDC_bigger; // save new, larger sized sample in "old" slot
}
else {
if(formerID>0) {
// make new digi
ZDCdigis->push_back(ZDCDataFrame(formerID));
// set up information to convert back
HcalDetId cell = ZDC_old.id();
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); // this one is generic
HcalCoderDb coder (*channelCoder, *shape);
unsigned int sizeold = ZDC_old.size();
for(unsigned int isamp = 0; isamp<sizeold; isamp++) {
coder.fC2adc(ZDC_old,(ZDCdigis->back()), 0 ); // as per simulation, capid=0???
}
}
//save pointers for next iteration
formerID = currentID;
ZDC_old = iZDC->second;
}
iZDCchk = iZDC;
if((++iZDCchk) == ZDCDigiStorage_.end()) { //make sure not to lose the last one
// make new digi
ZDCdigis->push_back(ZDCDataFrame(currentID));
// set up information to convert back
HcalDetId cell = (iZDC->second).id();
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); // this one is generic
HcalCoderDb coder (*channelCoder, *shape);
unsigned int sizeold = (iZDC->second).size();
for(unsigned int isamp = 0; isamp<sizeold; isamp++) {
coder.fC2adc(ZDC_old,(ZDCdigis->back()), 0 ); // as per simulation, capid=0???
}
}
}
//done merging
// put the collection of recunstructed hits in the event
LogInfo("DataMixingHcalDigiWorker") << "total # HBHE Merged digis: " << HBHEdigis->size() ;
LogInfo("DataMixingHcalDigiWorker") << "total # HO Merged digis: " << HOdigis->size() ;
LogInfo("DataMixingHcalDigiWorker") << "total # HF Merged digis: " << HFdigis->size() ;
LogInfo("DataMixingHcalDigiWorker") << "total # ZDC Merged digis: " << ZDCdigis->size() ;
e.put( HBHEdigis, HBHEDigiCollectionDM_ );
e.put( HOdigis, HODigiCollectionDM_ );
e.put( HFdigis, HFDigiCollectionDM_ );
e.put( ZDCdigis, ZDCDigiCollectionDM_ );
// clear local storage after this event
HBHEDigiStorage_.clear();
HODigiStorage_.clear();
HFDigiStorage_.clear();
ZDCDigiStorage_.clear();
}
bool edm::DataMixingHcalDigiWorker::DoZDC_ [private] |
Definition at line 89 of file DataMixingHcalDigiWorker.h.
Referenced by addHcalPileups(), addHcalSignals(), and DataMixingHcalDigiWorker().
std::string edm::DataMixingHcalDigiWorker::HBHEDigiCollectionDM_ [private] |
Definition at line 69 of file DataMixingHcalDigiWorker.h.
Referenced by DataMixingHcalDigiWorker(), and putHcal().
Definition at line 59 of file DataMixingHcalDigiWorker.h.
Referenced by addHcalSignals(), and DataMixingHcalDigiWorker().
Definition at line 84 of file DataMixingHcalDigiWorker.h.
Referenced by addHcalPileups(), addHcalSignals(), and putHcal().
Definition at line 64 of file DataMixingHcalDigiWorker.h.
Referenced by addHcalPileups(), and DataMixingHcalDigiWorker().
std::string edm::DataMixingHcalDigiWorker::HFDigiCollectionDM_ [private] |
Definition at line 71 of file DataMixingHcalDigiWorker.h.
Referenced by DataMixingHcalDigiWorker(), and putHcal().
Definition at line 61 of file DataMixingHcalDigiWorker.h.
Referenced by addHcalSignals(), and DataMixingHcalDigiWorker().
Definition at line 85 of file DataMixingHcalDigiWorker.h.
Referenced by addHcalPileups(), addHcalSignals(), and putHcal().
Definition at line 66 of file DataMixingHcalDigiWorker.h.
Referenced by addHcalPileups(), and DataMixingHcalDigiWorker().
std::string edm::DataMixingHcalDigiWorker::HODigiCollectionDM_ [private] |
Definition at line 70 of file DataMixingHcalDigiWorker.h.
Referenced by DataMixingHcalDigiWorker(), and putHcal().
Definition at line 60 of file DataMixingHcalDigiWorker.h.
Referenced by addHcalSignals(), and DataMixingHcalDigiWorker().
Definition at line 86 of file DataMixingHcalDigiWorker.h.
Referenced by addHcalPileups(), addHcalSignals(), and putHcal().
Definition at line 65 of file DataMixingHcalDigiWorker.h.
Referenced by addHcalPileups(), and DataMixingHcalDigiWorker().
std::string edm::DataMixingHcalDigiWorker::label_ [private] |
Definition at line 93 of file DataMixingHcalDigiWorker.h.
std::string edm::DataMixingHcalDigiWorker::ZDCDigiCollectionDM_ [private] |
Definition at line 72 of file DataMixingHcalDigiWorker.h.
Referenced by DataMixingHcalDigiWorker(), and putHcal().
Definition at line 62 of file DataMixingHcalDigiWorker.h.
Referenced by addHcalSignals(), and DataMixingHcalDigiWorker().
Definition at line 87 of file DataMixingHcalDigiWorker.h.
Referenced by addHcalPileups(), addHcalSignals(), and putHcal().
Definition at line 67 of file DataMixingHcalDigiWorker.h.
Referenced by addHcalPileups(), and DataMixingHcalDigiWorker().
1.7.3