#include <SiStripQuality.cc>

Inheritance diagram for SiStripQuality:

Classes
struct	BadComponent

class	BadComponentStrictWeakOrdering

Public Member Functions
void	add (const uint32_t &, const SiStripBadStrip::Range &)

void	add (const SiStripBadStrip *)

void	add (const SiStripDetCabling *)

void	add (const SiStripDetVOff *)

void	add (const RunInfo *)

void	addInvalidConnectionFromCabling ()

void	addNotConnectedConnectionFromCabling ()

SiStripDetCabling const *	cabling () const

bool	cleanUp (bool force=false)

void	clear ()

void	compact (unsigned int &, std::vector< unsigned int > &)

SiStripQuality	difference (const SiStripQuality &) const

void	fillBadComponents ()

short	getBadApvs (const uint32_t &detid) const

const std::vector< BadComponent > &	getBadComponentList () const

short	getBadFibers (const uint32_t &detid) const

bool	IsApvBad (const uint32_t &detid, const short &apvNb) const

bool	IsFiberBad (const uint32_t &detid, const short &fiberNb) const

bool	IsModuleBad (const uint32_t &detid) const

bool	IsModuleUsable (const uint32_t &detid) const

bool	IsStripBad (const uint32_t &detid, const short &strip) const

bool	IsStripBad (const Range &range, const short &strip) const

int	nBadStripsOnTheLeft (const Range &range, const short &strip) const

int	nBadStripsOnTheRight (const Range &range, const short &strip) const

void	ReduceGranularity (double)

void	setPrintDebugOutput (const bool printDebug)

void	setUseEmptyRunInfo (const bool useEmptyRunInfo)

	SiStripQuality ()=delete

	SiStripQuality (SiStripDetInfo)

	SiStripQuality (const SiStripQuality &)=default

	SiStripQuality (SiStripQuality &&)=default

	~SiStripQuality () override=default

Public Member Functions inherited from SiStripBadStrip
data	decode (const unsigned int &value) const

unsigned int	encode (const unsigned short &first, const unsigned short &NconsecutiveBadStrips, const unsigned short &flag=0)

ContainerIterator	getDataVectorBegin () const

ContainerIterator	getDataVectorEnd () const

void	getDetIds (std::vector< uint32_t > &DetIds_) const

const Range	getRange (const uint32_t detID) const

Range	getRangeByPos (unsigned short pos) const

RegistryIterator	getRegistryVectorBegin () const

RegistryIterator	getRegistryVectorEnd () const

void	printDebug (std::stringstream &ss, const TrackerTopology *trackerTopo) const

void	printSummary (std::stringstream &ss, const TrackerTopology *trackerTopo) const

bool	put (const uint32_t &detID, const InputVector &vect)

bool	put (const uint32_t &detID, Range input)

	SiStripBadStrip ()

	SiStripBadStrip (const SiStripBadStrip &orig)

virtual	~SiStripBadStrip ()

Private Member Functions
void	compact (std::vector< unsigned int > &, std::vector< unsigned int > &, unsigned short &)

void	printActiveFedsInfo (const std::vector< uint16_t > &activeFedsFromCabling, const std::vector< int > &activeFedsFromRunInfo, const std::vector< int > &differentFeds, const bool printDebug)

void	printDetInfo (const TrackerTopology *const tTopo, const uint32_t &detId, const uint32_t &apvPairNumber, std::stringstream &ss)
	Prints debug output for a given detId. More...

bool	put_replace (const uint32_t &DetId, Range input)

void	subtract (std::vector< unsigned int > &, const std::vector< unsigned int > &)

void	subtraction (std::vector< unsigned int > &, const unsigned int &)

void	turnOffFeds (const std::vector< int > &fedsList, const bool turnOffStrips, const bool printDebug)

Private Attributes
std::vector< BadComponent >	BadComponentVect

SiStripDetInfo	info_

bool	printDebug_

const SiStripDetCabling *	SiStripDetCabling_

bool	toCleanUp

bool	useEmptyRunInfo_

Additional Inherited Members
Public Types inherited from SiStripBadStrip
typedef std::vector< unsigned int >	Container

typedef std::vector< unsigned int >::const_iterator	ContainerIterator

typedef Container	InputVector

typedef std::pair < ContainerIterator, ContainerIterator >	Range

typedef std::vector< DetRegistry >	Registry

typedef Registry::const_iterator	RegistryIterator

Protected Attributes inherited from SiStripBadStrip
Registry	indexes

Container	v_badstrips

Detailed Description

Description: <one line="" class="" summary>="">

Usage: <usage>

Definition at line 32 of file SiStripQuality.h.

Constructor & Destructor Documentation

SiStripQuality::SiStripQuality ( )

delete

SiStripQuality::SiStripQuality ( SiStripDetInfo iInfo )

explicit

Definition at line 17 of file SiStripQuality.cc.

     : info_(std::move(iInfo)),
       toCleanUp(false),
       SiStripDetCabling_(nullptr),
       printDebug_(false),
       useEmptyRunInfo_(false) {}

SiStripQuality::SiStripQuality ( const SiStripQuality & )

default

SiStripQuality::SiStripQuality ( SiStripQuality && )

default

SiStripQuality::~SiStripQuality ( )

overridedefault

Member Function Documentation

void SiStripQuality::add	(	const uint32_t &	detid,
		const SiStripBadStrip::Range &	baserange
	)

Definition at line 243 of file SiStripQuality.cc.

References compact(), SiStripBadStrip::decode(), SiStripBadStrip::data::firstStrip, SiStripBadStrip::data::flag, SiStripDetInfo::getNumberOfApvsAndStripLength(), SiStripBadStrip::getRange(), info_, LogTrace, put_replace(), SiStripBadStrip::data::range, sistrip::SpyUtilities::range(), and createJobs::tmp.

Referenced by add(), addInvalidConnectionFromCabling(), addNotConnectedConnectionFromCabling(), ReduceGranularity(), counter.Counter::register(), SequenceTypes.Task::remove(), SequenceTypes.Task::replace(), and turnOffFeds().

                                                                                      {
   std::vector<unsigned int> vect, tmp;
 
   unsigned short Nstrips = info_.getNumberOfApvsAndStripLength(detid).first * 128;
 
   // Is this detid already in the collections owned by this class?
   SiStripBadStrip::Range range = getRange(detid);
 
   // Append bad strips
   tmp.clear();
   if (range.first == range.second) {
     LogTrace("SiStripQuality") << "new detid" << std::endl;
     // It's a new detid
     tmp.insert(tmp.end(), baserange.first, baserange.second);
     std::stable_sort(tmp.begin(), tmp.end());
     LogTrace("SiStripQuality") << "ordered" << std::endl;
   } else {
     LogTrace("SiStripQuality") << "already exists" << std::endl;
     // alredy existing detid
 
     // if full det is bad go to next detid
     SiStripBadStrip::data data_ = decode(*(range.first));
     if (range.second - range.first == 1 && data_.firstStrip == 0 && data_.range >= Nstrips) {
       LogTrace("SiStripQuality") << "full det is bad.. " << range.second - range.first << " "
                                  << decode(*(range.first)).firstStrip << " " << decode(*(range.first)).range << " "
                                  << decode(*(range.first)).flag << "\n"
                                  << std::endl;
       return;
     }
 
     tmp.insert(tmp.end(), baserange.first, baserange.second);
     tmp.insert(tmp.end(), range.first, range.second);
     std::stable_sort(tmp.begin(), tmp.end());
     LogTrace("SiStripQuality") << "ordered" << std::endl;
   }
   // Compact data
   compact(tmp, vect, Nstrips);
   SiStripBadStrip::Range newrange(vect.begin(), vect.end());
   if (!put_replace(detid, newrange))
     edm::LogError("SiStripQuality") << "[" << __PRETTY_FUNCTION__ << "] " << std::endl;
 }

void SiStripQuality::add ( const SiStripBadStrip * base )

Definition at line 226 of file SiStripQuality.cc.

References add(), SiStripBadStrip::getDataVectorBegin(), SiStripBadStrip::getRegistryVectorBegin(), SiStripBadStrip::getRegistryVectorEnd(), and LogTrace.

Referenced by counter.Counter::register(), SequenceTypes.Task::remove(), and SequenceTypes.Task::replace().

                                                     {
   SiStripBadStrip::RegistryIterator basebegin = base->getRegistryVectorBegin();
   SiStripBadStrip::RegistryIterator baseend = base->getRegistryVectorEnd();
 
   // the Registry already contains data
   // Loop on detids
   for (SiStripBadStrip::RegistryIterator basep = basebegin; basep != baseend; ++basep) {
     uint32_t detid = basep->detid;
     LogTrace("SiStripQuality") << "add detid " << detid << std::endl;
 
     SiStripBadStrip::Range baserange =
         SiStripBadStrip::Range(base->getDataVectorBegin() + basep->ibegin, base->getDataVectorBegin() + basep->iend);
 
     add(detid, baserange);
   }
 }

void SiStripQuality::add ( const SiStripDetCabling * cab )

Definition at line 158 of file SiStripQuality.cc.

References addInvalidConnectionFromCabling(), addNotConnectedConnectionFromCabling(), and SiStripDetCabling_.

Referenced by counter.Counter::register(), SequenceTypes.Task::remove(), and SequenceTypes.Task::replace().

                                                      {
   SiStripDetCabling_ = cab;
   addInvalidConnectionFromCabling();
   addNotConnectedConnectionFromCabling();
 }

void SiStripQuality::add ( const SiStripDetVOff * Voff )

Definition at line 63 of file SiStripQuality.cc.

References add(), SiStripBadStrip::encode(), SiStripDetVOff::getDetIds(), SiStripDetInfo::getNumberOfApvsAndStripLength(), info_, LogTrace, and sistrip::SpyUtilities::range().

Referenced by counter.Counter::register(), SequenceTypes.Task::remove(), and SequenceTypes.Task::replace().

                                                    {
   std::vector<unsigned int> vect;
   short firstStrip = 0;
   short range = 0;
 
   // Get vector of Voff dets
   std::vector<uint32_t> vdets;
   Voff->getDetIds(vdets);
   std::vector<uint32_t>::const_iterator iter = vdets.begin();
   std::vector<uint32_t>::const_iterator iterEnd = vdets.end();
 
   for (; iter != iterEnd; ++iter) {
     vect.clear();
     range = (short)(info_.getNumberOfApvsAndStripLength(*iter).first * 128.);
     LogTrace("SiStripQuality") << "[add Voff] add detid " << *iter << " first strip " << firstStrip << " range "
                                << range << std::endl;
     vect.push_back(encode(firstStrip, range));
     SiStripBadStrip::Range Range(vect.begin(), vect.end());
     add(*iter, Range);
   }
 }

void SiStripQuality::add ( const RunInfo * runInfo )

Used to get the cabling from RunInfo.
It compares the feds from cabling with those from runInfo and it turns off all the strips associated to feds that are off for the RunInfo but not for FedCabling.

Definition at line 85 of file SiStripQuality.cc.

References trackerTree::check(), SiStripDetCabling::fedCabling(), SiStripFedCabling::fedIds(), RunInfo::m_fed_in, FEDNumbering::MAXSiStripFEDID, FEDNumbering::MINSiStripFEDID, printDebug_, DBoxMetadataHelper::set_difference(), SiStripDetCabling_, contentValuesCheck::ss, turnOffFeds(), useEmptyRunInfo_, and x.

Referenced by counter.Counter::register(), SequenceTypes.Task::remove(), and SequenceTypes.Task::replace().

                                                {
   bool allFedsEmpty = runInfo->m_fed_in.empty();
   if (allFedsEmpty) {
     std::stringstream ss;
     ss << "WARNING: the full list of feds in RunInfo is empty. ";
     if (useEmptyRunInfo_) {
       ss << " SiStripQuality will still use it and all tracker will be off." << std::endl;
     } else {
       ss << " SiStripQuality will not use it." << std::endl;
     }
     edm::LogInfo("SiStripQuality") << ss.str();
   }
 
   if (!allFedsEmpty || useEmptyRunInfo_) {
     // Take the list of active feds from fedCabling
     auto ids = SiStripDetCabling_->fedCabling()->fedIds();
 
     std::vector<uint16_t> activeFedsFromCabling(ids.begin(), ids.end());
     // Take the list of active feds from RunInfo
     std::vector<int> activeFedsFromRunInfo;
     // Take only Tracker feds (remove all non Tracker)
     std::remove_copy_if(
         runInfo->m_fed_in.begin(), runInfo->m_fed_in.end(), std::back_inserter(activeFedsFromRunInfo), [&](int x) {
           return !((x >= int(FEDNumbering::MINSiStripFEDID)) && (x <= int(FEDNumbering::MAXSiStripFEDID)));
         });
 
     // Compare the two. If a fedId from RunInfo is not present in the fedCabling
     // we need to get all the corresponding fedChannels and then the single apv
     // pairs and use them to turn off the corresponding strips (apvNumber*256).
     // set_difference returns the set of elements that are in the first and not
     // in the second
     std::sort(activeFedsFromCabling.begin(), activeFedsFromCabling.end());
     std::sort(activeFedsFromRunInfo.begin(), activeFedsFromRunInfo.end());
     std::vector<int> differentFeds;
     // Take the feds active for cabling but not for runInfo
     std::set_difference(activeFedsFromCabling.begin(),
                         activeFedsFromCabling.end(),
                         activeFedsFromRunInfo.begin(),
                         activeFedsFromRunInfo.end(),
                         std::back_inserter(differentFeds));
 
     // IGNORE for time being.
     // printActiveFedsInfo(activeFedsFromCabling, activeFedsFromRunInfo,
     // differentFeds, printDebug_);
 
     // Feds in the differentFeds vector are now to be turned off as they are off
     // according to RunInfo but were not off in cabling and thus are still
     // active for the SiStripQuality. The "true" means that the strips are to be
     // set as bad.
     turnOffFeds(differentFeds, true, printDebug_);
 
     // Consistency check
     // -----------------
     std::vector<int> check;
     std::set_difference(activeFedsFromRunInfo.begin(),
                         activeFedsFromRunInfo.end(),
                         activeFedsFromCabling.begin(),
                         activeFedsFromCabling.end(),
                         std::back_inserter(check));
     // This must not happen
     if (!check.empty()) {
       // throw cms::Exception("LogicError")
       edm::LogInfo("SiStripQuality") << "The cabling should always include the active feds in runInfo and "
                                         "possibly have some more"
                                      << "there are instead " << check.size() << " feds only active in runInfo";
       // The "false" means that we are only printing the output, but not setting
       // the strips as bad. The second bool means that we always want the debug
       // output in this case.
       turnOffFeds(check, false, true);
     }
   }
 }

void SiStripQuality::addInvalidConnectionFromCabling ( )

Definition at line 182 of file SiStripQuality.cc.

References add(), SiStripDetCabling::addActiveDetectorsRawIds(), SiStripBadStrip::encode(), SiStripDetCabling::getConnections(), sistrip::invalid_, LogTrace, SiStripDetCabling::nApvPairs(), sistrip::SpyUtilities::range(), and SiStripDetCabling_.

Referenced by add().

                                                      {
   std::vector<uint32_t> connected_detids;
   SiStripDetCabling_->addActiveDetectorsRawIds(connected_detids);
   std::vector<uint32_t>::const_iterator itdet = connected_detids.begin();
   std::vector<uint32_t>::const_iterator itdetEnd = connected_detids.end();
   for (; itdet != itdetEnd; ++itdet) {
     // LogTrace("SiStripQuality") << "[addInvalidConnectionFromCabling] looking
     // at detid " <<*itdet << std::endl;
     const std::vector<const FedChannelConnection *> &fedconns = SiStripDetCabling_->getConnections(*itdet);
     std::vector<const FedChannelConnection *>::const_iterator itconns = fedconns.begin();
     std::vector<const FedChannelConnection *>::const_iterator itconnsEnd = fedconns.end();
 
     unsigned short nApvPairs = SiStripDetCabling_->nApvPairs(*itdet);
     short ngoodConn = 0, goodConn = 0;
     for (; itconns != itconnsEnd; ++itconns) {
       // LogTrace("SiStripQuality") << "[addInvalidConnectionFromCabling]
       // apvpair " << (*itconns)->apvPairNumber() << " napvpair " <<
       // (*itconns)->nApvPairs()<< " detid " << (*itconns)->detId() <<
       // std::endl;
       if ((*itconns == nullptr) || ((*itconns)->nApvPairs() == sistrip::invalid_))
         continue;
       ngoodConn++;
       goodConn = goodConn | (0x1 << (*itconns)->apvPairNumber());
     }
 
     if (ngoodConn != nApvPairs) {
       std::vector<unsigned int> vect;
       for (size_t idx = 0; idx < nApvPairs; ++idx) {
         if (!(goodConn & (0x1 << idx))) {
           short firstStrip = idx * 256;
           short range = 256;
           LogTrace("SiStripQuality") << "[addInvalidConnectionFromCabling] add detid " << *itdet << "firstStrip "
                                      << firstStrip << std::endl;
           vect.push_back(encode(firstStrip, range));
         }
       }
       if (!vect.empty()) {
         SiStripBadStrip::Range Range(vect.begin(), vect.end());
         add(*itdet, Range);
       }
     }
   }
 }

void SiStripQuality::addNotConnectedConnectionFromCabling ( )

Definition at line 164 of file SiStripQuality.cc.

References add(), SiStripBadStrip::encode(), SiStripDetInfo::getAllData(), info_, SiStripDetCabling::IsConnected(), LogTrace, sistrip::SpyUtilities::range(), and SiStripDetCabling_.

Referenced by add().

                                                           {
   auto allData = info_.getAllData();
   auto iter = allData.begin();
   auto iterEnd = allData.end();
   std::vector<unsigned int> vect;
   short firstStrip = 0;
   short range = 0;
   for (; iter != iterEnd; ++iter)
     if (!SiStripDetCabling_->IsConnected(iter->first)) {
       vect.clear();
       range = iter->second.nApvs * 128;
       LogTrace("SiStripQuality") << "[addNotConnectedConnectionFromCabling] add detid " << iter->first << std::endl;
       vect.push_back(encode(firstStrip, range));
       SiStripBadStrip::Range Range(vect.begin(), vect.end());
       add(iter->first, Range);
     }
 }

SiStripDetCabling const* SiStripQuality::cabling ( ) const

inline

Definition at line 121 of file SiStripQuality.h.

References SiStripDetCabling_.

121 { return SiStripDetCabling_; }

SiStripQuality::SiStripDetCabling_

const SiStripDetCabling * SiStripDetCabling_

Definition: SiStripQuality.h:151

bool SiStripQuality::cleanUp ( bool force = false )

Definition at line 459 of file SiStripQuality.cc.

References SiStripBadStrip::indexes, LogTrace, SiStripBadStrip::put(), sistrip::SpyUtilities::range(), toCleanUp, and SiStripBadStrip::v_badstrips.

Referenced by difference().

                                        {
   if (!toCleanUp && !force)
     return false;
 
   toCleanUp = false;
 
   std::vector<unsigned int> v_badstrips_tmp = v_badstrips;
   std::vector<DetRegistry> indexes_tmp = indexes;
 
   LogTrace("SiStripQuality") << "[SiStripQuality::cleanUp] before cleanUp v_badstrips.size()= " << v_badstrips.size()
                              << " indexes.size()=" << indexes.size() << std::endl;
 
   v_badstrips.clear();
   indexes.clear();
 
   SiStripBadStrip::RegistryIterator basebegin = indexes_tmp.begin();
   SiStripBadStrip::RegistryIterator baseend = indexes_tmp.end();
 
   for (SiStripBadStrip::RegistryIterator basep = basebegin; basep != baseend; ++basep) {
     if (basep->ibegin != basep->iend) {
       SiStripBadStrip::Range range(v_badstrips_tmp.begin() + basep->ibegin, v_badstrips_tmp.begin() + basep->iend);
       if (!put(basep->detid, range))
         edm::LogError("SiStripQuality") << "[" << __PRETTY_FUNCTION__ << "] " << std::endl;
     }
   }
 
   LogTrace("SiStripQuality") << "[SiStripQuality::cleanUp] after cleanUp v_badstrips.size()= " << v_badstrips.size()
                              << " indexes.size()=" << indexes.size() << std::endl;
   return true;
 }

void SiStripQuality::clear ( void )

inline

Definition at line 53 of file SiStripQuality.h.

References BadComponentVect, SiStripBadStrip::indexes, toCleanUp, and SiStripBadStrip::v_badstrips.

Referenced by SiStripBadAPVAlgorithmFromClusterOccupancy::extractBadAPVs(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::extractBadAPVSandStrips(), SiStripHotStripAlgorithmFromClusterOccupancy::extractBadStrips(), and BeautifulSoup.Tag::setString().

                {
     v_badstrips.clear();
     indexes.clear();
     BadComponentVect.clear();
     toCleanUp = false;
   }

void SiStripQuality::compact	(	unsigned int &	detid,
		std::vector< unsigned int > &	vect
	)

Definition at line 285 of file SiStripQuality.cc.

References SiStripDetInfo::getNumberOfApvsAndStripLength(), info_, and createJobs::tmp.

Referenced by add(), SiStripBadAPVAlgorithmFromClusterOccupancy::AnalyzeOccupancy(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::AnalyzeOccupancy(), SiStripHotStripAlgorithmFromClusterOccupancy::extractBadStrips(), SiStripHitEffFromCalibTree::makeSQLite(), and SiStripCondObjBuilderFromDb::storeQuality().

                                                                                {
   std::vector<unsigned int> tmp = vect;
   vect.clear();
   std::stable_sort(tmp.begin(), tmp.end());
   unsigned short Nstrips = info_.getNumberOfApvsAndStripLength(detid).first * 128;
   compact(tmp, vect, Nstrips);
 }

void SiStripQuality::compact	(	std::vector< unsigned int > &	tmp,
		std::vector< unsigned int > &	vect,
		unsigned short &	Nstrips
	)

private

Definition at line 381 of file SiStripQuality.cc.

References SiStripBadStrip::decode(), SiStripBadStrip::encode(), SiStripBadStrip::data::firstStrip, and SiStripBadStrip::data::range.

                                                                                                                  {
   SiStripBadStrip::data fs_0, fs_1;
   vect.clear();
 
   ContainerIterator it = tmp.begin();
   fs_0 = decode(*it);
 
   // Check if at the module end
   if (fs_0.firstStrip + fs_0.range >= Nstrips) {
     vect.push_back(encode(fs_0.firstStrip, Nstrips - fs_0.firstStrip));
     return;
   }
 
   ++it;
   for (; it != tmp.end(); ++it) {
     fs_1 = decode(*it);
 
     if (fs_0.firstStrip + fs_0.range >= fs_1.firstStrip + fs_1.range) {
       // fs_0 includes fs_1, go ahead
     } else if (fs_0.firstStrip + fs_0.range >= fs_1.firstStrip) {
       // contiguous or superimposed intervals
       // Check if at the module end
       if (fs_1.firstStrip + fs_1.range >= Nstrips) {
         vect.push_back(encode(fs_0.firstStrip, Nstrips - fs_0.firstStrip));
         return;
       } else {
         // create new fs_0
         fs_0.range = fs_1.firstStrip + fs_1.range - fs_0.firstStrip;
       }
     } else {
       // separated intervals
       vect.push_back(encode(fs_0.firstStrip, fs_0.range));
       fs_0 = fs_1;
     }
   }
   vect.push_back(encode(fs_0.firstStrip, fs_0.range));
 }

SiStripQuality SiStripQuality::difference ( const SiStripQuality & other ) const

Definition at line 24 of file SiStripQuality.cc.

References cleanUp(), SiStripBadStrip::decode(), fillBadComponents(), SiStripBadStrip::data::firstStrip, SiStripBadStrip::getDataVectorBegin(), SiStripDetInfo::getNumberOfApvsAndStripLength(), SiStripBadStrip::getRange(), SiStripBadStrip::getRegistryVectorBegin(), SiStripBadStrip::getRegistryVectorEnd(), info_, put_replace(), SiStripBadStrip::data::range, sistrip::SpyUtilities::range(), and subtract().

                                                                            {
   SiStripBadStrip::RegistryIterator rbegin = other.getRegistryVectorBegin();
   SiStripBadStrip::RegistryIterator rend = other.getRegistryVectorEnd();
   std::vector<unsigned int> ovect, vect;
   uint32_t detid;
   unsigned short Nstrips;
 
   SiStripQuality retValue(*this);
   for (SiStripBadStrip::RegistryIterator rp = rbegin; rp != rend; ++rp) {
     detid = rp->detid;
     Nstrips = info_.getNumberOfApvsAndStripLength(detid).first * 128;
 
     SiStripBadStrip::Range orange =
         SiStripBadStrip::Range(other.getDataVectorBegin() + rp->ibegin, other.getDataVectorBegin() + rp->iend);
 
     // Is this detid already in the collections owned by this class?
     SiStripBadStrip::Range range = retValue.getRange(detid);
     if (range.first != range.second) {  // yes, it is
 
       vect.clear();
       ovect.clear();
 
       // if other full det is bad, remove det from this
       SiStripBadStrip::data data_ = decode(*(orange.first));
       if (orange.second - orange.first != 1 || data_.firstStrip != 0 || data_.range < Nstrips) {
         ovect.insert(ovect.end(), orange.first, orange.second);
         vect.insert(vect.end(), range.first, range.second);
         retValue.subtract(vect, ovect);
       }
       SiStripBadStrip::Range newrange(vect.begin(), vect.end());
       if (!retValue.put_replace(detid, newrange))
         edm::LogError("SiStripQuality") << "[" << __PRETTY_FUNCTION__ << "] " << std::endl;
     }
   }
   retValue.cleanUp();
   retValue.fillBadComponents();
   return retValue;
 }

void SiStripQuality::fillBadComponents ( )

Definition at line 490 of file SiStripQuality.cc.

References SiStripQuality::BadComponent::BadApvs, BadComponentVect, SiStripQuality::BadComponent::BadFibers, SiStripQuality::BadComponent::BadModule, SiStripBadStrip::decode(), SiStripQuality::BadComponent::detid, SiStripBadStrip::data::firstStrip, SiStripDetInfo::getNumberOfApvsAndStripLength(), SiStripBadStrip::indexes, info_, SiStripBadStrip::data::range, sistrip::SpyUtilities::range(), mps_fire::result, and SiStripBadStrip::v_badstrips.

Referenced by difference(), SiStripBadAPVAlgorithmFromClusterOccupancy::extractBadAPVs(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::extractBadAPVSandStrips(), SiStripHotStripAlgorithmFromClusterOccupancy::extractBadStrips(), and SiStripHitEffFromCalibTree::makeSQLite().

                                        {
   BadComponentVect.clear();
 
   for (SiStripBadStrip::RegistryIterator basep = indexes.begin(); basep != indexes.end(); ++basep) {
     SiStripBadStrip::Range range(v_badstrips.begin() + basep->ibegin, v_badstrips.begin() + basep->iend);
 
     // Fill BadModules, BadFibers, BadApvs vectors
     unsigned short resultA = 0, resultF = 0;
     BadComponent result;
 
     SiStripBadStrip::data fs;
     unsigned short Nstrips = info_.getNumberOfApvsAndStripLength(basep->detid).first * 128;
 
     // BadModules
     fs = decode(*(range.first));
     if (basep->iend - basep->ibegin == 1 && fs.firstStrip == 0 && fs.range == Nstrips) {
       result.detid = basep->detid;
       result.BadModule = true;
       result.BadFibers = (1 << (Nstrips / 256)) - 1;
       result.BadApvs = (1 << (Nstrips / 128)) - 1;
 
       BadComponentVect.push_back(result);
 
     } else {
       // Bad Fibers and  Apvs
       for (SiStripBadStrip::ContainerIterator it = range.first; it != range.second; ++it) {
         fs = decode(*it);
 
         // BadApvs
         for (short apvNb = 0; apvNb < 6; ++apvNb) {
           if (fs.firstStrip <= apvNb * 128 && (apvNb + 1) * 128 <= fs.firstStrip + fs.range) {
             resultA = resultA | (1 << apvNb);
           }
         }
         // BadFibers
         for (short fiberNb = 0; fiberNb < 3; ++fiberNb) {
           if (fs.firstStrip <= fiberNb * 256 && (fiberNb + 1) * 256 <= fs.firstStrip + fs.range) {
             resultF = resultF | (1 << fiberNb);
           }
         }
       }
       if (resultA != 0) {
         result.detid = basep->detid;
         result.BadModule = false;
         result.BadFibers = resultF;
         result.BadApvs = resultA;
         BadComponentVect.push_back(result);
       }
     }
   }
 }

short SiStripQuality::getBadApvs ( const uint32_t & detid ) const

Definition at line 626 of file SiStripQuality.cc.

References BadComponentVect, cuda_std::lower_bound(), and AlCaHLTBitMon_ParallelJobs::p.

Referenced by SiStripRecHitConverterAlgorithm::fillBad128StripBlocks(), and MeasurementTrackerImpl::initializeStripStatus().

                                                             {
   std::vector<BadComponent>::const_iterator p = std::lower_bound(
       BadComponentVect.begin(), BadComponentVect.end(), detid, SiStripQuality::BadComponentStrictWeakOrdering());
   if (p != BadComponentVect.end() && p->detid == detid)
     return p->BadApvs;
   return 0;
 }

const std::vector<BadComponent>& SiStripQuality::getBadComponentList ( ) const

inline

Definition at line 114 of file SiStripQuality.h.

References BadComponentVect.

Referenced by SiStripHitEffFromCalibTree::algoAnalyze(), SiStripBadComponentInfo::fillBadComponentMaps(), SiStripPI::fillBCArrays(), and SiStripQualityStatistics::updateAndSave().

114 { return BadComponentVect; }

SiStripQuality::BadComponentVect

std::vector< BadComponent > BadComponentVect

Definition: SiStripQuality.h:149

short SiStripQuality::getBadFibers ( const uint32_t & detid ) const

Definition at line 634 of file SiStripQuality.cc.

References BadComponentVect, cuda_std::lower_bound(), and AlCaHLTBitMon_ParallelJobs::p.

Referenced by SiStripRecHitConverterAlgorithm::fillBad128StripBlocks(), and MeasurementTrackerImpl::initializeStripStatus().

                                                               {
   std::vector<BadComponent>::const_iterator p = std::lower_bound(
       BadComponentVect.begin(), BadComponentVect.end(), detid, SiStripQuality::BadComponentStrictWeakOrdering());
   if (p != BadComponentVect.end() && p->detid == detid)
     return p->BadFibers;
   return 0;
 }

bool SiStripQuality::IsApvBad	(	const uint32_t &	detid,
		const short &	apvNb
	)		const

Definition at line 574 of file SiStripQuality.cc.

References BadComponentVect, cuda_std::lower_bound(), and AlCaHLTBitMon_ParallelJobs::p.

Referenced by SiStripBadAPVAlgorithmFromClusterOccupancy::AnalyzeOccupancy(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::AnalyzeOccupancy(), SiStripHotStripAlgorithmFromClusterOccupancy::extractBadStrips(), SiStripAPVRestorer::inspect(), and SiStripClusterInfo::IsApvBad().

                                                                              {
   std::vector<BadComponent>::const_iterator p = std::lower_bound(
       BadComponentVect.begin(), BadComponentVect.end(), detid, SiStripQuality::BadComponentStrictWeakOrdering());
   if (p != BadComponentVect.end() && p->detid == detid)
     return ((p->BadApvs >> apvNb) & 0x1);
   return false;
 }

bool SiStripQuality::IsFiberBad	(	const uint32_t &	detid,
		const short &	fiberNb
	)		const

Definition at line 566 of file SiStripQuality.cc.

References BadComponentVect, cuda_std::lower_bound(), and AlCaHLTBitMon_ParallelJobs::p.

Referenced by SiStripClusterInfo::IsFiberBad().

                                                                                  {
   std::vector<BadComponent>::const_iterator p = std::lower_bound(
       BadComponentVect.begin(), BadComponentVect.end(), detid, SiStripQuality::BadComponentStrictWeakOrdering());
   if (p != BadComponentVect.end() && p->detid == detid)
     return ((p->BadFibers >> fiberNb) & 0x1);
   return false;
 }

bool SiStripQuality::IsModuleBad ( const uint32_t & detid ) const

Definition at line 558 of file SiStripQuality.cc.

References BadComponentVect, cuda_std::lower_bound(), and AlCaHLTBitMon_ParallelJobs::p.

Referenced by MCMisalignmentScaler::analyze(), TrackerGeometryCompare::fillTree(), LargeEvents< T >::filter(), SingleMultiplicity< T >::getEvent(), SiStripClusterizerConditions::isModuleBad(), SiStripClusterInfo::IsModuleBad(), and SiStripTrackerMapCreator::setTkMapFromAlarm().

                                                             {
   std::vector<BadComponent>::const_iterator p = std::lower_bound(
       BadComponentVect.begin(), BadComponentVect.end(), detid, SiStripQuality::BadComponentStrictWeakOrdering());
   if (p != BadComponentVect.end() && p->detid == detid)
     return p->BadModule;
   return false;
 }

bool SiStripQuality::IsModuleUsable ( const uint32_t & detid ) const

Definition at line 544 of file SiStripQuality.cc.

References BadComponentVect, SiStripDetCabling::IsConnected(), cuda_std::lower_bound(), AlCaHLTBitMon_ParallelJobs::p, and SiStripDetCabling_.

Referenced by MeasurementTrackerImpl::initializeStripStatus(), SiStripClusterizerConditions::isModuleUsable(), SiStripClusterInfo::IsModuleUsable(), SiStripQualityStatistics::updateAndSave(), and SiStripRecHitConverterAlgorithm::useModule().

                                                                {
   std::vector<BadComponent>::const_iterator p = std::lower_bound(
       BadComponentVect.begin(), BadComponentVect.end(), detid, SiStripQuality::BadComponentStrictWeakOrdering());
   if (p != BadComponentVect.end() && p->detid == detid)
     if (p->BadModule)
       return false;
 
   if (SiStripDetCabling_ != nullptr)
     if (!SiStripDetCabling_->IsConnected(detid))
       return false;
 
   return true;
 }

bool SiStripQuality::IsStripBad	(	const uint32_t &	detid,
		const short &	strip
	)		const

Definition at line 582 of file SiStripQuality.cc.

References SiStripBadStrip::getRange(), and sistrip::SpyUtilities::range().

Referenced by SiStripAPVRestorer::abnormalBaselineInspect(), SiStripMonitorCluster::analyze(), SiStripClusterizerConditions::Det::bad(), SiStripMonitorTrack::clusterInfos(), SiStripAPVRestorer::nullInspect(), SiStripClusterInfo::stripQualitiesBad(), TT6CMNSubtractor::subtract_(), and IteratedMedianCMNSubtractor::subtract_().

                                                                                {
   SiStripBadStrip::Range range = getRange(detid);
   return IsStripBad(range, strip);
 }

bool SiStripQuality::IsStripBad	(	const Range &	range,
		const short &	strip
	)		const

Definition at line 587 of file SiStripQuality.cc.

References SiStripBadStrip::decode(), SiStripBadStrip::data::firstStrip, SiStripBadStrip::data::range, and mps_fire::result.

                                                                             {
   bool result = false;
   SiStripBadStrip::data fs;
   for (SiStripBadStrip::ContainerIterator it = range.first; it != range.second; ++it) {
     fs = decode(*it);
     if ((fs.firstStrip <= strip) & (strip < fs.firstStrip + fs.range)) {
       result = true;
       break;
     }
   }
   return result;
 }

int SiStripQuality::nBadStripsOnTheLeft	(	const Range &	range,
		const short &	strip
	)		const

Definition at line 600 of file SiStripQuality.cc.

References SiStripBadStrip::decode(), SiStripBadStrip::data::firstStrip, SiStripBadStrip::data::range, and mps_fire::result.

                                                                                     {
   int result = 0;
   SiStripBadStrip::data fs;
   for (SiStripBadStrip::ContainerIterator it = range.first; it != range.second; ++it) {
     fs = decode(*it);
     if (fs.firstStrip <= strip && strip < fs.firstStrip + fs.range) {
       result = strip - fs.firstStrip + 1;
       break;
     }
   }
   return result;
 }

int SiStripQuality::nBadStripsOnTheRight	(	const Range &	range,
		const short &	strip
	)		const

Definition at line 613 of file SiStripQuality.cc.

References SiStripBadStrip::decode(), SiStripBadStrip::data::firstStrip, SiStripBadStrip::data::range, mps_fire::result, and digitizers_cfi::strip.

                                                                                      {
   int result = 0;
   SiStripBadStrip::data fs;
   for (SiStripBadStrip::ContainerIterator it = range.first; it != range.second; ++it) {
     fs = decode(*it);
     if (fs.firstStrip <= strip && strip < fs.firstStrip + fs.range) {
       result = fs.firstStrip + fs.range - strip;
       break;
     }
   }
   return result;
 }

void SiStripQuality::printActiveFedsInfo	(	const std::vector< uint16_t > &	activeFedsFromCabling,
		const std::vector< int > &	activeFedsFromRunInfo,
		const std::vector< int > &	differentFeds,
		const bool	printDebug
	)

private

Prints debug output for the active feds comparing the list in RunInfo and FedCabling

Definition at line 676 of file SiStripQuality.cc.

References filterCSVwithJSON::copy, and contentValuesCheck::ss.

                                                                 {
   std::ostringstream ss;
 
   if (printDebug) {
     ss << "activeFedsFromCabling:" << std::endl;
     std::copy(activeFedsFromCabling.begin(), activeFedsFromCabling.end(), std::ostream_iterator<uint16_t>(ss, " "));
     ss << std::endl;
     ss << "activeFedsFromRunInfo:" << std::endl;
     std::copy(activeFedsFromRunInfo.begin(), activeFedsFromRunInfo.end(), std::ostream_iterator<int>(ss, " "));
     ss << std::endl;
   }
   if (differentFeds.size() != 440) {
     ss << "differentFeds : " << std::endl;
     std::copy(differentFeds.begin(), differentFeds.end(), std::ostream_iterator<int>(ss, " "));
     ss << std::endl;
   } else {
     ss << "There are 440 feds (all) active for Cabling but off for RunInfo. "
           "Tracker was probably not in this run"
        << std::endl;
   }
   edm::LogInfo("SiStripQuality") << ss.str() << std::endl;
 }

void SiStripQuality::printDetInfo	(	const TrackerTopology *const	tTopo,
		const uint32_t &	detId,
		const uint32_t &	apvPairNumber,
		std::stringstream &	ss
	)

private

Prints debug output for a given detId.

Definition at line 642 of file SiStripQuality.cc.

References phase1PixelTopology::layer, TrackerTopology::layer(), AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), StripSubdetector::TEC, TrackerTopology::tecIsStereo(), StripSubdetector::TIB, TrackerTopology::tibIsStereo(), StripSubdetector::TID, TrackerTopology::tidIsStereo(), StripSubdetector::TOB, and TrackerTopology::tobIsStereo().

Referenced by turnOffFeds().

                                                        {
   std::string subDetName;
   DetId detid(detId);
   int layer = tTopo->layer(detid);
   int stereo = 0;
   switch (detid.subdetId()) {
     case StripSubdetector::TIB: {
       stereo = tTopo->tibIsStereo(detid);
       subDetName = "TIB";
       break;
     }
     case StripSubdetector::TOB: {
       stereo = tTopo->tobIsStereo(detid);
       subDetName = "TOB";
       break;
     }
     case StripSubdetector::TEC: {
       stereo = tTopo->tecIsStereo(detid);
       subDetName = "TEC";
       break;
     }
     case StripSubdetector::TID: {
       stereo = tTopo->tidIsStereo(detid);
       subDetName = "TID";
       break;
     }
   }
   ss << detId << " and apv = " << apvPairNumber << " of subDet = " << subDetName << ", layer = " << layer
      << " stereo = " << stereo << std::endl;
 }

bool SiStripQuality::put_replace	(	const uint32_t &	DetId,
		Range	input
	)

private

Definition at line 293 of file SiStripQuality.cc.

References SiStripBadStrip::DetRegistry::detid, SiStripBadStrip::DetRegistry::ibegin, SiStripBadStrip::DetRegistry::iend, SiStripBadStrip::indexes, LogTrace, cuda_std::lower_bound(), AlCaHLTBitMon_ParallelJobs::p, sd, toCleanUp, and SiStripBadStrip::v_badstrips.

Referenced by add(), and difference().

                                                                    {
   // put in SiStripQuality::v_badstrips of DetId
   Registry::iterator p = std::lower_bound(indexes.begin(), indexes.end(), DetId, SiStripBadStrip::StrictWeakOrdering());
 
   size_t sd = input.second - input.first;
   DetRegistry detregistry;
   detregistry.detid = DetId;
   detregistry.ibegin = v_badstrips.size();
   detregistry.iend = v_badstrips.size() + sd;
 
   v_badstrips.insert(v_badstrips.end(), input.first, input.second);
 
   if (p != indexes.end() && p->detid == DetId) {
     LogTrace("SiStripQuality") << "[SiStripQuality::put_replace]  Replacing "
                                   "SiStripQuality for already stored DetID "
                                << DetId << std::endl;
     toCleanUp = true;
     *p = detregistry;
   } else {
     indexes.insert(p, detregistry);
   }
 
   return true;
 }

void SiStripQuality::ReduceGranularity ( double threshold )

Definition at line 324 of file SiStripQuality.cc.

References add(), SiStripBadStrip::decode(), SiStripBadStrip::encode(), SiStripBadStrip::data::firstStrip, SiStripBadStrip::getDataVectorBegin(), SiStripBadStrip::getRegistryVectorBegin(), SiStripBadStrip::getRegistryVectorEnd(), mps_fire::i, LogTrace, SiStripPI::max, min(), and SiStripBadStrip::data::range.

                                                        {
   SiStripBadStrip::RegistryIterator rp = getRegistryVectorBegin();
   SiStripBadStrip::RegistryIterator rend = getRegistryVectorEnd();
   SiStripBadStrip::data data_;
   uint16_t BadStripPerApv[6], ipos;
   std::vector<unsigned int> vect;
 
   for (; rp != rend; ++rp) {
     uint32_t detid = rp->detid;
 
     BadStripPerApv[0] = 0;
     BadStripPerApv[1] = 0;
     BadStripPerApv[2] = 0;
     BadStripPerApv[3] = 0;
     BadStripPerApv[4] = 0;
     BadStripPerApv[5] = 0;
     ipos = 0;
 
     SiStripBadStrip::Range sqrange =
         SiStripBadStrip::Range(getDataVectorBegin() + rp->ibegin, getDataVectorBegin() + rp->iend);
 
     for (int it = 0; it < sqrange.second - sqrange.first; it++) {
       data_ = decode(*(sqrange.first + it));
       LogTrace("SiStripQuality") << "[SiStripQuality::ReduceGranularity] detid " << detid << " first strip "
                                  << data_.firstStrip << " lastStrip " << data_.firstStrip + data_.range - 1 << " range "
                                  << data_.range;
       ipos = data_.firstStrip / 128;
       while (ipos <= (data_.firstStrip + data_.range - 1) / 128) {
         BadStripPerApv[ipos] +=
             std::min(data_.firstStrip + data_.range, (ipos + 1) * 128) - std::max(data_.firstStrip * 1, ipos * 128);
         LogTrace("SiStripQuality") << "[SiStripQuality::ReduceGranularity] ipos " << ipos << " Counter "
                                    << BadStripPerApv[ipos] << " min "
                                    << std::min(data_.firstStrip + data_.range, (ipos + 1) * 128) << " max "
                                    << std::max(data_.firstStrip * 1, ipos * 128) << " added "
                                    << std::min(data_.firstStrip + data_.range, (ipos + 1) * 128) -
                                           std::max(data_.firstStrip * 1, ipos * 128);
         ipos++;
       }
     }
 
     LogTrace("SiStripQuality") << "[SiStripQuality::ReduceGranularity] Total for detid " << detid << " values "
                                << BadStripPerApv[0] << " " << BadStripPerApv[1] << " " << BadStripPerApv[2] << " "
                                << BadStripPerApv[3] << " " << BadStripPerApv[4] << " " << BadStripPerApv[5];
 
     vect.clear();
     for (size_t i = 0; i < 6; ++i) {
       if (BadStripPerApv[i] >= threshold * 128) {
         vect.push_back(encode(i * 128, 128));
       }
     }
     if (!vect.empty()) {
       SiStripBadStrip::Range Range(vect.begin(), vect.end());
       add(detid, Range);
     }
   }
 }

void SiStripQuality::setPrintDebugOutput ( const bool printDebug )

inline

Definition at line 118 of file SiStripQuality.h.

References SiStripBadStrip::printDebug(), and printDebug_.

118 { printDebug_ = printDebug; }

SiStripBadStrip::printDebug

void printDebug(std::stringstream &ss, const TrackerTopology *trackerTopo) const

Definition: SiStripBadStrip.cc:76

SiStripQuality::printDebug_

bool printDebug_

Definition: SiStripQuality.h:152

void SiStripQuality::setUseEmptyRunInfo ( const bool useEmptyRunInfo )

inline

Definition at line 119 of file SiStripQuality.h.

References useEmptyRunInfo_.

119 { useEmptyRunInfo_ = useEmptyRunInfo; }

SiStripQuality::useEmptyRunInfo_

bool useEmptyRunInfo_

Definition: SiStripQuality.h:153

void SiStripQuality::subtract	(	std::vector< unsigned int > &	A,
		const std::vector< unsigned int > &	B
	)

private

Definition at line 419 of file SiStripQuality.cc.

References subtraction().

Referenced by difference().

                                                                                           {
   ContainerIterator it = B.begin();
   ContainerIterator itend = B.end();
   for (; it != itend; ++it) {
     subtraction(A, *it);
   }
 }

void SiStripQuality::subtraction	(	std::vector< unsigned int > &	A,
		const unsigned int &	B
	)

private

Definition at line 427 of file SiStripQuality.cc.

References SiStripBadStrip::decode(), SiStripBadStrip::encode(), SiStripBadStrip::data::firstStrip, SiStripBadStrip::data::range, and createJobs::tmp.

Referenced by subtract().

                                                                                   {
   SiStripBadStrip::data fs_A, fs_B, fs_m, fs_M;
   std::vector<unsigned int> tmp;
 
   fs_B = decode(B);
   ContainerIterator jt = A.begin();
   ContainerIterator jtend = A.end();
   for (; jt != jtend; ++jt) {
     fs_A = decode(*jt);
     if (B < *jt) {
       fs_m = fs_B;
       fs_M = fs_A;
     } else {
       fs_m = fs_A;
       fs_M = fs_B;
     }
     // A) Verify the range to be subtracted crosses the new range
     if (fs_m.firstStrip + fs_m.range > fs_M.firstStrip) {
       if (*jt < B) {
         tmp.push_back(encode(fs_A.firstStrip, fs_B.firstStrip - fs_A.firstStrip));
       }
       if (fs_A.firstStrip + fs_A.range > fs_B.firstStrip + fs_B.range) {
         tmp.push_back(
             encode(fs_B.firstStrip + fs_B.range, fs_A.firstStrip + fs_A.range - (fs_B.firstStrip + fs_B.range)));
       }
     } else {
       tmp.push_back(*jt);
     }
   }
   A = tmp;
 }

void SiStripQuality::turnOffFeds	(	const std::vector< int > &	fedsList,
		const bool	turnOffStrips,
		const bool	printDebug
	)

private

Loop on all the fedIds, take the fedChannels and then the detId. Depending on the value of a bool turn off or not the strips in the list.

Definition at line 702 of file SiStripQuality.cc.

References add(), SiStripBadStrip::encode(), SiStripDetCabling::fedCabling(), SiStripFedCabling::fedConnections(), LogTrace, printDetInfo(), SiStripDetCabling_, contentValuesCheck::ss, and SiStripDetCabling::trackerTopology().

Referenced by add().

                                                                                                                 {
   std::stringstream ss;
   if (printDebug) {
     ss << "associated to detIds : " << std::endl;
   }
 
   std::vector<int>::const_iterator fedIdIt = fedsList.begin();
   for (; fedIdIt != fedsList.end(); ++fedIdIt) {
     std::vector<FedChannelConnection>::const_iterator fedChIt =
         SiStripDetCabling_->fedCabling()->fedConnections(*fedIdIt).begin();
     for (; fedChIt != SiStripDetCabling_->fedCabling()->fedConnections(*fedIdIt).end(); ++fedChIt) {
       uint32_t detId = fedChIt->detId();
       if (detId == 0 || detId == 0xFFFFFFFF)
         continue;
       uint16_t apvPairNumber = fedChIt->apvPairNumber();
 
       if (printDebug) {
         printDetInfo(SiStripDetCabling_->trackerTopology(), detId, apvPairNumber, ss);
       }
 
       if (turnOffStrips) {
         // apvPairNumber == i it means that the i*256 strips are to be set off
         std::vector<unsigned int> vect;
         vect.push_back(encode(apvPairNumber * 256, 256));
         SiStripBadStrip::Range Range(vect.begin(), vect.end());
         add(detId, Range);
         LogTrace("SiStripQuality") << "[addOffForRunInfo] adding apvPairNumber " << apvPairNumber << " for detId "
                                    << detId << " off according to RunInfo" << std::endl;
       }
     }
   }
   if (printDebug) {
     edm::LogInfo("SiStripQuality") << ss.str() << std::endl;
   }
 }