|
|
|
#include <SiStripBadModuleGenerator.h>
Public Member Functions | |
| void | getObj (SiStripBadStrip *&obj) |
| SiStripBadModuleGenerator (const edm::ParameterSet &, const edm::ActivityRegistry &) | |
| ~SiStripBadModuleGenerator () | |
Private Types | |
| typedef std::vector < edm::ParameterSet > | Parameters |
Private Member Functions | |
| void | createObject () |
| bool | isTECDetector (const uint32_t &therawid, uint32_t requested_side, uint32_t requested_wheel, uint32_t requested_petal_bkw_frw, uint32_t requested_petal, uint32_t requested_ring, uint32_t requested_ster, uint32_t requested_detid) const |
| bool | isTIBDetector (const uint32_t &therawid, uint32_t requested_layer, uint32_t requested_bkw_frw, uint32_t requested_int_ext, uint32_t requested_string, uint32_t requested_ster, uint32_t requested_detid) const |
| bool | isTIDDetector (const uint32_t &therawid, uint32_t requested_side, uint32_t requested_wheel, uint32_t requested_ring, uint32_t requested_ster, uint32_t requested_detid) const |
| bool | isTOBDetector (const uint32_t &therawid, uint32_t requested_layer, uint32_t requested_bkw_frw, uint32_t requested_rod, uint32_t requested_ster, uint32_t requested_detid) const |
| void | selectDetectors (const std::vector< uint32_t > &, std::vector< uint32_t > &) |
Private Attributes | |
| Parameters | BadComponentList_ |
| bool | printdebug_ |
Definition at line 10 of file SiStripBadModuleGenerator.h.
typedef std::vector< edm::ParameterSet > SiStripBadModuleGenerator::Parameters [private] |
Definition at line 56 of file SiStripBadModuleGenerator.h.
| SiStripBadModuleGenerator::SiStripBadModuleGenerator | ( | const edm::ParameterSet & | iConfig, |
| const edm::ActivityRegistry & | aReg | ||
| ) | [explicit] |
Definition at line 16 of file SiStripBadModuleGenerator.cc.
: SiStripCondObjBuilderBase<SiStripBadStrip>::SiStripCondObjBuilderBase(iConfig) { edm::LogInfo("SiStripBadModuleGenerator") << "[SiStripBadModuleGenerator::SiStripBadModuleGenerator]"; }
| SiStripBadModuleGenerator::~SiStripBadModuleGenerator | ( | ) |
Definition at line 23 of file SiStripBadModuleGenerator.cc.
{
edm::LogInfo("SiStripBadModuleGenerator") << "[SiStripBadModuleGenerator::~SiStripBadModuleGenerator]";
}
| void SiStripBadModuleGenerator::createObject | ( | ) | [private] |
Definition at line 28 of file SiStripBadModuleGenerator.cc.
References SiStripCondObjBuilderBase< SiStripBadStrip >::_pset, BadComponentList_, SiStripQuality::cleanUp(), SiStripBadStrip::encode(), edm::FileInPath::fullPath(), SiStripQuality::getBadComponentList(), edm::ParameterSet::getUntrackedParameter(), VarParsing::obj, SiStripCondObjBuilderBase< SiStripBadStrip >::obj_, printdebug_, SiStripBadStrip::put(), matplotRender::reader, selectDetectors(), and AlCaHarvesting_cff::SiStripQuality.
Referenced by getObj().
{
SiStripQuality* obj = new SiStripQuality();
edm::FileInPath fp_ = _pset.getUntrackedParameter<edm::FileInPath>("file",edm::FileInPath("CalibTracker/SiStripCommon/data/SiStripDetInfo.dat"));
printdebug_ = _pset.getUntrackedParameter<bool>("printDebug",false);
BadComponentList_ = _pset.getUntrackedParameter<Parameters>("BadComponentList");
SiStripDetInfoFileReader reader(fp_.fullPath());
const std::vector<uint32_t>& DetIds= reader.getAllDetIds();
std::vector<uint32_t> selDetIds;
selectDetectors(DetIds,selDetIds);
edm::LogInfo("SiStripQualityConfigurableFakeESSource")<<"[produce] number of selected dets to be removed " << selDetIds.size() <<std::endl;
std::stringstream ss;
std::vector<uint32_t>::const_iterator iter=selDetIds.begin();
std::vector<uint32_t>::const_iterator iterEnd=selDetIds.end();
for(;iter!=iterEnd;++iter){
SiStripQuality::InputVector theSiStripVector;
unsigned short firstBadStrip=0, NconsecutiveBadStrips=reader.getNumberOfApvsAndStripLength(*iter).first * 128;
unsigned int theBadStripRange;
theBadStripRange = obj->encode(firstBadStrip,NconsecutiveBadStrips);
if (printdebug_)
ss << "detid " << *iter << " \t"
<< " firstBadStrip " << firstBadStrip << "\t "
<< " NconsecutiveBadStrips " << NconsecutiveBadStrips << "\t "
<< " packed integer " << std::hex << theBadStripRange << std::dec
<< std::endl;
theSiStripVector.push_back(theBadStripRange);
SiStripBadStrip::Range range(theSiStripVector.begin(),theSiStripVector.end());
if ( ! obj->put(*iter,range) )
edm::LogError("SiStripQualityConfigurableFakeESSource")<<"[produce] detid already exists"<<std::endl;
}
if (printdebug_)
edm::LogInfo("SiStripQualityConfigurableFakeESSource") << ss.str();
obj->cleanUp();
//obj->fillBadComponents();
std::stringstream ss1;
if (printdebug_){
for (std::vector<SiStripQuality::BadComponent>::const_iterator iter=obj->getBadComponentList().begin();iter!=obj->getBadComponentList().end();++iter)
ss1 << "bad module " << iter->detid << " " << iter->BadModule << "\n";
edm::LogInfo("SiStripQualityConfigurableFakeESSource") << ss1.str();
}
obj_ = new SiStripBadStrip( *(dynamic_cast<SiStripBadStrip*> (obj)));
delete obj;
}
| void SiStripBadModuleGenerator::getObj | ( | SiStripBadStrip *& | obj | ) | [inline, virtual] |
Returns the CondObj
Reimplemented from SiStripCondObjBuilderBase< SiStripBadStrip >.
Definition at line 16 of file SiStripBadModuleGenerator.h.
References createObject(), and SiStripCondObjBuilderBase< SiStripBadStrip >::obj_.
{createObject(); obj=obj_;}
| bool SiStripBadModuleGenerator::isTECDetector | ( | const uint32_t & | therawid, |
| uint32_t | requested_side, | ||
| uint32_t | requested_wheel, | ||
| uint32_t | requested_petal_bkw_frw, | ||
| uint32_t | requested_petal, | ||
| uint32_t | requested_ring, | ||
| uint32_t | requested_ster, | ||
| uint32_t | requested_detid | ||
| ) | const [private] |
Definition at line 255 of file SiStripBadModuleGenerator.cc.
References TECDetId::isFrontPetal(), TECDetId::isStereo(), TECDetId::isZPlusSide(), TECDetId::petalNumber(), DetId::rawId(), TECDetId::ringNumber(), SiStripDetId::subDetector(), sistripsummary::TEC, and TECDetId::wheelNumber().
Referenced by selectDetectors().
{
TECDetId potentialDet = TECDetId(therawid); // build TECDetId, at this point is just DetId, but do not want to cast twice
if( potentialDet.subDetector() == SiStripDetId::TEC ){ // check if subdetector field is a TEC, both tested numbers are int
if( // check if TEC is from the ones requested
( (potentialDet.wheelNumber()==requested_wheel) || requested_wheel==0 ) // take everything if default value is 0
&&
( (potentialDet.isZPlusSide() && requested_side==2) || (!potentialDet.isZPlusSide() && requested_side==1) || requested_side==0)
&&
( (potentialDet.isStereo() && requested_ster==1) || (!potentialDet.isStereo() && requested_ster==2) || requested_ster==0 )
&&
( (potentialDet.isFrontPetal() && requested_petal_bkw_frw==2) || (!potentialDet.isFrontPetal() && requested_petal_bkw_frw==2) || requested_petal_bkw_frw==0 )
&&
( (potentialDet.petalNumber()==requested_petal) || requested_petal==0 )
&&
( (potentialDet.ringNumber()==requested_ring) || requested_ring==0 )
&&
( (potentialDet.rawId()==requested_detid) || requested_detid==0 )
)
return 1;
}
return 0;
}
| bool SiStripBadModuleGenerator::isTIBDetector | ( | const uint32_t & | therawid, |
| uint32_t | requested_layer, | ||
| uint32_t | requested_bkw_frw, | ||
| uint32_t | requested_int_ext, | ||
| uint32_t | requested_string, | ||
| uint32_t | requested_ster, | ||
| uint32_t | requested_detid | ||
| ) | const [private] |
Definition at line 178 of file SiStripBadModuleGenerator.cc.
References TIBDetId::isInternalString(), TIBDetId::isRPhi(), TIBDetId::isStereo(), TIBDetId::isZPlusSide(), TIBDetId::layerNumber(), DetId::rawId(), TIBDetId::stringNumber(), SiStripDetId::subDetector(), and sistripsummary::TIB.
Referenced by selectDetectors().
{
TIBDetId potentialDet = TIBDetId(therawid); // build TIBDetId, at this point is just DetId, but do not want to cast twice
if( potentialDet.subDetector() == SiStripDetId::TIB ){ // check if subdetector field is a TIB, both tested numbers are int
if( // check if TIB is from the ones requested
( (potentialDet.layerNumber()==requested_layer) || requested_layer==0 ) // take everything if default value is 0
&&
( (potentialDet.isZPlusSide() && requested_bkw_frw==2) || (!potentialDet.isZPlusSide() && requested_bkw_frw==1) || requested_bkw_frw==0)
&&
( (potentialDet.isInternalString() && requested_int_ext==1) || (!potentialDet.isInternalString() && requested_int_ext==2) || requested_int_ext==0 )
&&
( (potentialDet.isStereo() && requested_ster==1) || (potentialDet.isRPhi() && requested_ster==2) || requested_ster==0 )
&&
( (potentialDet.stringNumber()==requested_string) || requested_string==0 )
&&
( (potentialDet.rawId()==requested_detid) || requested_detid==0 )
)
return 1;
}
return 0;
}
| bool SiStripBadModuleGenerator::isTIDDetector | ( | const uint32_t & | therawid, |
| uint32_t | requested_side, | ||
| uint32_t | requested_wheel, | ||
| uint32_t | requested_ring, | ||
| uint32_t | requested_ster, | ||
| uint32_t | requested_detid | ||
| ) | const [private] |
Definition at line 230 of file SiStripBadModuleGenerator.cc.
References TIDDetId::diskNumber(), TIDDetId::isRPhi(), TIDDetId::isStereo(), TIDDetId::isZPlusSide(), DetId::rawId(), TIDDetId::ringNumber(), SiStripDetId::subDetector(), and sistripsummary::TID.
Referenced by selectDetectors().
{
TIDDetId potentialDet = TIDDetId(therawid); // build TIDDetId, at this point is just DetId, but do not want to cast twice
if( potentialDet.subDetector() == SiStripDetId::TID ){ // check if subdetector field is a TID, both tested numbers are int
if( // check if TID is from the ones requested
( (potentialDet.diskNumber()==requested_wheel) || requested_wheel==0 ) // take everything if default value is 0
&&
( (potentialDet.isZPlusSide() && requested_side==2) || (!potentialDet.isZPlusSide() && requested_side==1) || requested_side==0)
&&
( (potentialDet.isStereo() && requested_ster==1) || (potentialDet.isRPhi() && requested_ster==2) || requested_ster==0 )
&&
( (potentialDet.ringNumber()==requested_ring) || requested_ring==0 )
&&
( (potentialDet.rawId()==requested_detid) || requested_detid==0 )
)
return 1;
}
return 0;
}
| bool SiStripBadModuleGenerator::isTOBDetector | ( | const uint32_t & | therawid, |
| uint32_t | requested_layer, | ||
| uint32_t | requested_bkw_frw, | ||
| uint32_t | requested_rod, | ||
| uint32_t | requested_ster, | ||
| uint32_t | requested_detid | ||
| ) | const [private] |
Definition at line 205 of file SiStripBadModuleGenerator.cc.
References TOBDetId::isRPhi(), TOBDetId::isStereo(), TOBDetId::isZPlusSide(), TOBDetId::layerNumber(), DetId::rawId(), TOBDetId::rodNumber(), SiStripDetId::subDetector(), and sistripsummary::TOB.
Referenced by selectDetectors().
{
TOBDetId potentialDet = TOBDetId(therawid); // build TOBDetId, at this point is just DetId, but do not want to cast twice
if( potentialDet.subDetector() == SiStripDetId::TOB ){ // check if subdetector field is a TOB, both tested numbers are int
if( // check if TOB is from the ones requested
( (potentialDet.layerNumber()==requested_layer) || requested_layer==0 ) // take everything if default value is 0
&&
( (potentialDet.isZPlusSide() && requested_bkw_frw==2) || (!potentialDet.isZPlusSide() && requested_bkw_frw==1) || requested_bkw_frw==0)
&&
( (potentialDet.isStereo() && requested_ster==1) || (potentialDet.isRPhi() && requested_ster==2) || requested_ster==0 )
&&
( (potentialDet.rodNumber()==requested_rod) || requested_rod==0 )
&&
( (potentialDet.rawId()==requested_detid) || requested_detid==0 )
)
return 1;
}
return 0;
}
| void SiStripBadModuleGenerator::selectDetectors | ( | const std::vector< uint32_t > & | DetIds, |
| std::vector< uint32_t > & | list | ||
| ) | [private] |
Definition at line 87 of file SiStripBadModuleGenerator.cc.
References BadComponentList_, gather_cfg::cout, spr::find(), isTECDetector(), isTIBDetector(), isTIDDetector(), isTOBDetector(), printdebug_, DetId::rawId(), SiStripDetId::TEC, SiStripDetId::TIB, SiStripDetId::TID, SiStripDetId::TOB, DetId::Tracker, and pat::UNKNOWN.
Referenced by createObject().
{
SiStripSubStructure siStripSubStructure;
std::stringstream ss;
for(Parameters::iterator iBadComponent = BadComponentList_.begin(); iBadComponent != BadComponentList_.end(); ++iBadComponent ) {
if (printdebug_)
ss << "Bad SubDet " << iBadComponent->getParameter<std::string>("SubDet") << " \t";
SiStripDetId::SubDetector subDet=SiStripDetId::UNKNOWN;
if (iBadComponent->getParameter<std::string>("SubDet")=="TIB")
subDet=SiStripDetId::TIB;
else if (iBadComponent->getParameter<std::string>("SubDet")=="TID")
subDet=SiStripDetId::TID;
else if (iBadComponent->getParameter<std::string>("SubDet")=="TOB")
subDet=SiStripDetId::TOB;
else if (iBadComponent->getParameter<std::string>("SubDet")=="TEC")
subDet=SiStripDetId::TEC;
std::vector<uint32_t> genericBadDetIds( iBadComponent->getUntrackedParameter<std::vector<uint32_t> >("detidList", std::vector<uint32_t>()) );
bool anySubDet = true;
if( genericBadDetIds.empty() ) anySubDet = false;
std::cout << "genericBadDetIds.size() = " << genericBadDetIds.size() << std::endl;
uint32_t startDet=DetId(DetId::Tracker,subDet).rawId();
uint32_t stopDet=DetId(DetId::Tracker,subDet+1).rawId();
if( anySubDet ) {
startDet=DetId(DetId::Tracker,SiStripDetId::TIB).rawId();
stopDet=DetId(DetId::Tracker,SiStripDetId::TEC+1).rawId();
}
std::vector<uint32_t>::const_iterator iter=lower_bound(DetIds.begin(),DetIds.end(),startDet);
std::vector<uint32_t>::const_iterator iterEnd=lower_bound(DetIds.begin(),DetIds.end(),stopDet);
bool resp;
for ( ;iter!=iterEnd;++iter){
resp=false;
if (iBadComponent->getParameter<std::string>("SubDet")=="TIB")
resp=isTIBDetector(*iter,
iBadComponent->getParameter<uint32_t>("layer"),
iBadComponent->getParameter<uint32_t>("bkw_frw"),
iBadComponent->getParameter<uint32_t>("int_ext"),
iBadComponent->getParameter<uint32_t>("ster"),
iBadComponent->getParameter<uint32_t>("string_"),
iBadComponent->getParameter<uint32_t>("detid")
);
else if (iBadComponent->getParameter<std::string>("SubDet")=="TID")
resp=isTIDDetector(*iter,
iBadComponent->getParameter<uint32_t>("side"),
iBadComponent->getParameter<uint32_t>("wheel"),
iBadComponent->getParameter<uint32_t>("ring"),
iBadComponent->getParameter<uint32_t>("ster"),
iBadComponent->getParameter<uint32_t>("detid")
);
else if (iBadComponent->getParameter<std::string>("SubDet")=="TOB")
resp=isTOBDetector(*iter,
iBadComponent->getParameter<uint32_t>("layer"),
iBadComponent->getParameter<uint32_t>("bkw_frw"),
iBadComponent->getParameter<uint32_t>("rod"),
iBadComponent->getParameter<uint32_t>("ster"),
iBadComponent->getParameter<uint32_t>("detid")
);
else if (iBadComponent->getParameter<std::string>("SubDet")=="TEC")
resp=isTECDetector(*iter,
iBadComponent->getParameter<uint32_t>("side"),
iBadComponent->getParameter<uint32_t>("wheel"),
iBadComponent->getParameter<uint32_t>("petal_bkw_frw"),
iBadComponent->getParameter<uint32_t>("petal"),
iBadComponent->getParameter<uint32_t>("ring"),
iBadComponent->getParameter<uint32_t>("ster"),
iBadComponent->getParameter<uint32_t>("detid")
);
if( anySubDet ) {
std::cout << "AnySubDet" << *iter << std::endl;
if( std::find(genericBadDetIds.begin(), genericBadDetIds.end(), *iter) == genericBadDetIds.end() ) resp = false;
else resp = true;
}
if(resp)
list.push_back(*iter);
}
}
if (printdebug_)
edm::LogInfo("SiStripBadModuleGenerator") << ss.str();
}
Definition at line 57 of file SiStripBadModuleGenerator.h.
Referenced by createObject(), and selectDetectors().
bool SiStripBadModuleGenerator::printdebug_ [private] |
Definition at line 55 of file SiStripBadModuleGenerator.h.
Referenced by createObject(), and selectDetectors().
1.7.3