CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Member Functions | Private Attributes
SiStripQualityHotStripIdentifierRoot Class Reference

#include <SiStripQualityHotStripIdentifierRoot.h>

Inheritance diagram for SiStripQualityHotStripIdentifierRoot:
ConditionDBWriter< SiStripBadStrip > edm::EDAnalyzer edm::EDConsumerBase

Public Member Functions

 SiStripQualityHotStripIdentifierRoot (const edm::ParameterSet &)
 
 ~SiStripQualityHotStripIdentifierRoot ()
 
- Public Member Functions inherited from ConditionDBWriter< SiStripBadStrip >
 ConditionDBWriter (const edm::ParameterSet &iConfig)
 
virtual ~ConditionDBWriter ()
 
- Public Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzer ()
 
ModuleDescription const & moduleDescription () const
 
std::string workerType () const
 
virtual ~EDAnalyzer ()
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
 EDConsumerBase ()
 
ProductHolderIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
std::vector
< ProductHolderIndexAndSkipBit >
const & 		itemsToGetFromEvent () const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
bool registeredToConsume (ProductHolderIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
 
virtual ~EDConsumerBase ()
 

Private Member Functions

void algoAnalyze (const edm::Event &, const edm::EventSetup &)
 
void algoBeginLuminosityBlock (const edm::LuminosityBlock &, const edm::EventSetup &)
 
void algoBeginRun (const edm::Run &, const edm::EventSetup &)
 
void algoEndJob ()
 
void bookHistos ()
 
SiStripBadStrip * getNewObject ()
 

Private Attributes

const TrackerGeometry * _tracker
 
uint32_t calibrationthreshold
 
SiStrip::QualityHistosMap ClusterPositionHistoMap
 
const edm::ParameterSet conf_
 
std::string dataLabel_
 
std::string dirpath
 
DQMStore * dqmStore_
 
TFile * file0
 
std::string filename
 
edm::FileInPath fp_
 
unsigned long long m_cacheID_
 
unsigned short MaxClusterWidth_
 
double MeanNumberOfCluster
 
unsigned short MinClusterWidth_
 
SiStripDetInfoFileReader * reader
 
edm::ESHandle< SiStripQuality > SiStripQuality_
 
SiStripHotStripAlgorithmFromClusterOccupancy * theIdentifier
 
SiStripBadAPVAlgorithmFromClusterOccupancy * theIdentifier2
 
SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy * theIdentifier3
 
edm::ESHandle< TrackerGeometry > theTrackerGeom
 
double TotNumberOfEvents
 
const TrackerTopology * tTopo
 
bool UseInputDB_
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &)
 
- Protected Member Functions inherited from ConditionDBWriter< SiStripBadStrip >
void setDoStore (const bool doStore)
 When set to false the payload will not be written to the db. More...
 
void storeOnDbNow ()
 
cond::Time_t timeOfLastIOV ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Definition at line 48 of file SiStripQualityHotStripIdentifierRoot.h.

Constructor & Destructor Documentation

SiStripQualityHotStripIdentifierRoot::SiStripQualityHotStripIdentifierRoot ( const edm::ParameterSet &  iConfig)
explicit

Definition at line 22 of file SiStripQualityHotStripIdentifierRoot.cc.

References dqmStore_, filename, fp_, edm::FileInPath::fullPath(), edm::ParameterSet::getUntrackedParameter(), cppFunctionSkipper::operator, reader, and OfflineSiStripConfigDb_cff::SiStripDetInfoFileReader.

22  : ConditionDBWriter<SiStripBadStrip>(iConfig),
23  m_cacheID_(0),
24  dataLabel_(iConfig.getUntrackedParameter<std::string>("dataLabel","")),
25  UseInputDB_(iConfig.getUntrackedParameter<bool>("UseInputDB",false)),
26  conf_(iConfig),
27  fp_(iConfig.getUntrackedParameter<edm::FileInPath>("file",edm::FileInPath("CalibTracker/SiStripCommon/data/SiStripDetInfo.dat"))),
28  _tracker(nullptr),
29  tTopo(nullptr),
30  filename(iConfig.getUntrackedParameter<std::string>("rootFilename","CondDB_TKCC_20X_v3_hlt_50822.root")),
31  dirpath(iConfig.getUntrackedParameter<std::string>("rootDirPath","")),
32  TotNumberOfEvents(0),
33  MeanNumberOfCluster(0),
34  calibrationthreshold(iConfig.getUntrackedParameter<uint32_t>("CalibrationThreshold",10000))
35 {
36  reader = new SiStripDetInfoFileReader(fp_.fullPath());
37 
38  dqmStore_ = edm::Service<DQMStore>().operator->();
39  dqmStore_->setVerbose(iConfig.getUntrackedParameter<uint32_t>("verbosity",0));
40 
41  if(filename!=""){
42  edm::LogInfo("SiStripQualityHotStripIdentifierRoot") << " before opening file " << filename.c_str();
43  dqmStore_->open(filename.c_str(), false);
44  edm::LogInfo("SiStripQualityHotStripIdentifierRoot") << " after opening file ";
45  }
46 }
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
DQMStore::setVerbose
void setVerbose(unsigned level)
Definition: DQMStore.cc:691
OfflineSiStripConfigDb_cff.SiStripDetInfoFileReader
tuple SiStripDetInfoFileReader
Definition: OfflineSiStripConfigDb_cff.py:20
DQMStore::open
bool open(const std::string &filename, bool overwrite=false, const std::string &path="", const std::string &prepend="", OpenRunDirs stripdirs=KeepRunDirs, bool fileMustExist=true)
Definition: DQMStore.cc:2993
edm::FileInPath::fullPath
std::string fullPath() const
Definition: FileInPath.cc:184
SiStripQualityHotStripIdentifierRoot::~SiStripQualityHotStripIdentifierRoot ( )

Definition at line 49 of file SiStripQualityHotStripIdentifierRoot.cc.

49  {
50 }

Member Function Documentation

void SiStripQualityHotStripIdentifierRoot::algoAnalyze ( const edm::Event &  ,
const edm::EventSetup &   
)
inlineprivatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 66 of file SiStripQualityHotStripIdentifierRoot.h.

66 {};
void SiStripQualityHotStripIdentifierRoot::algoBeginLuminosityBlock ( const edm::LuminosityBlock &  ,
const edm::EventSetup &   
)
inlineprivatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 60 of file SiStripQualityHotStripIdentifierRoot.h.

60 { }
void SiStripQualityHotStripIdentifierRoot::algoBeginRun ( const edm::Run &  iRun,
const edm::EventSetup &  iSetup 
)
privatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 201 of file SiStripQualityHotStripIdentifierRoot.cc.

References _tracker, dataLabel_, edm::EventSetup::get(), m_cacheID_, edm::ESHandle< class >::product(), SiStripQuality_, theTrackerGeom, tTopo, and UseInputDB_.

201  {
202  //Retrieve tracker topology from geometry
203  edm::ESHandle<TrackerTopology> tTopoHandle;
204  iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
205  tTopo = tTopoHandle.product();
206 
207  iSetup.get<TrackerDigiGeometryRecord> ().get (theTrackerGeom);
208  _tracker=&(* theTrackerGeom);
209 
210  if(UseInputDB_){
211  unsigned long long cacheID = iSetup.get<SiStripQualityRcd>().cacheIdentifier();
212 
213  if (m_cacheID_ == cacheID)
214  return;
215 
216  m_cacheID_ = cacheID;
217 
218  iSetup.get<SiStripQualityRcd>().get(dataLabel_,SiStripQuality_);
219  }
220 }
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:56
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
void SiStripQualityHotStripIdentifierRoot::algoEndJob ( )
privatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 222 of file SiStripQualityHotStripIdentifierRoot.cc.

References ClusterPositionHistoMap.

222  {
223  //Clear map
224  ClusterPositionHistoMap.clear();
225 }
void SiStripQualityHotStripIdentifierRoot::bookHistos ( )
private

Definition at line 228 of file SiStripQualityHotStripIdentifierRoot.cc.

References EnergyCorrector::c, ClusterPositionHistoMap, conf_, dirpath, dqmStore_, DQMStore::getAllContents(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), LogDebug, MeanNumberOfCluster, NULL, HLT_25ns10e33_v2_cff::parameters, AlCaHLTBitMon_QueryRunRegistry::string, indexGen::title, and TotNumberOfEvents.

Referenced by getNewObject().

228  {
229  edm::LogInfo("SiStripQualityHotStripIdentifierRoot") <<" [SiStripQualityHotStripIdentifierRoot::bookHistos] " << dirpath << std::endl;
230  std::vector<MonitorElement*> MEs = dqmStore_->getAllContents(dirpath);
231  //"DQMData/Run 50908/SiStrip/MechanicalView");
232  // "/DQMData/Run 50908/SiStrip/Run summary/MechanicalView/TID/side_2/wheel_3/ring_2/mono_modules/module_402676874");
233  edm::LogInfo("SiStripQualityHotStripIdentifierRoot") <<" [SiStripQualityHotStripIdentifierRoot::bookHistos] vector size " << MEs.size()<< std::endl;
234 
235  std::vector<MonitorElement*>::const_iterator iter=MEs.begin();
236  std::vector<MonitorElement*>::const_iterator iterEnd=MEs.end();
237 
238  edm::ParameterSet parameters=conf_.getParameter<edm::ParameterSet>("AlgoParameters");
239  bool gotNentries=true;
240  if( parameters.getUntrackedParameter<uint32_t>("NumberOfEvents",0)==0 && parameters.getUntrackedParameter<double>("OccupancyThreshold",0)!=0)
241  gotNentries=false;
242  edm::LogInfo("SiStripQualityHotStripIdentifierRoot")<< "[SiStripQualityHotStripIdentifierRoot::bookHistos] gotNentries flag " << gotNentries
243  << " number of " << parameters.getUntrackedParameter<uint32_t>("NumberOfEvents",0)
244  << " occup " << parameters.getUntrackedParameter<double>("OccupancyThreshold",0)
245  << " OccupancyHisto" << parameters.getUntrackedParameter<std::string>("OccupancyHisto") << std::endl;
246 
247  // Check Number of Events
248  for (; iter!=iterEnd;++iter) {
249  std::string me_name = (*iter)->getName();
250 
251  if (!gotNentries && strstr(me_name.c_str(),"TotalNumberOfCluster__T")!=NULL && strstr(me_name.c_str(),"Profile")==NULL ){
252 
253  TotNumberOfEvents = ((TH1F*)(*iter)->getTH1F())->GetEntries();
254  MeanNumberOfCluster = ((TH1F*)(*iter)->getTH1F())->GetMean();
255  edm::LogInfo("SiStripQualityHotStripIdentifierRoot")<< "Total Number of Events: " << TotNumberOfEvents << std::endl;
256 
257  gotNentries=true;
258  edm::LogInfo("SiStripQualityHotStripIdentifierRoot")<< "[SiStripQualityHotStripIdentifierRoot::bookHistos] gotNentries flag " << gotNentries << std::endl;
259  break;
260  }
261  }
262  if (!gotNentries) {
263  edm::LogWarning("SiStripQualityHotStripIdentifierRoot") <<" [SiStripQualityHotStripIdentifierRoot::bookHistos] :: Histogram with to check # of evemnts missing" <<std::endl;
264  }
265  for (; iter!=iterEnd;++iter) {
266  std::string me_name = (*iter)->getName();
267 
268  if (strstr(me_name.c_str(),(parameters.getUntrackedParameter<std::string>("OccupancyHisto")).c_str())==NULL)
269  continue;
270 
271  unsigned int detid=0;
272  char title[128];
273  sprintf(title,"%s",me_name.c_str());
274  char *ptr=strtok(title,"__");
275  int c=0;
276  while (ptr!=NULL){
277  if (c==2){
278  detid=atol(ptr);
279  break;
280  }
281  ptr=strtok(NULL,"_");
282  c++;
283  }
284  LogDebug("SiStripQualityHotStripIdentifierRoot") <<" [SiStripQualityHotStripIdentifierRoot::bookHistos] detid " << detid<< std::endl;
285 
286  ClusterPositionHistoMap[detid]=boost::shared_ptr<TH1F>(new TH1F(*(*iter)->getTH1F()));
287 
288  }
289 
290 }
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:501
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
NULL
#define NULL
Definition: scimark2.h:8
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
DQMStore::getAllContents
std::vector< MonitorElement * > getAllContents(const std::string &path, uint32_t runNumber=0, uint32_t lumi=0) const
Definition: DQMStore.cc:1988
HLT_25ns10e33_v2_cff.parameters
tuple parameters
Definition: HLT_25ns10e33_v2_cff.py:3864
SiStripBadStrip * SiStripQualityHotStripIdentifierRoot::getNewObject ( )
privatevirtual

Implements ConditionDBWriter< SiStripBadStrip >.

Definition at line 52 of file SiStripQualityHotStripIdentifierRoot.cc.

References _tracker, a, HLT_25ns10e33_v2_cff::AlgoName, bookHistos(), calibrationthreshold, ClusterPositionHistoMap, conf_, SiStripBadAPVAlgorithmFromClusterOccupancy::extractBadAPVs(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::extractBadAPVSandStrips(), SiStripHotStripAlgorithmFromClusterOccupancy::extractBadStrips(), SiStripBadStrip::getDataVectorBegin(), edm::ParameterSet::getParameter(), SiStripBadStrip::getRegistryVectorBegin(), SiStripBadStrip::getRegistryVectorEnd(), edm::ParameterSet::getUntrackedParameter(), MeanNumberOfCluster, getGTfromDQMFile::obj, HLT_25ns10e33_v2_cff::parameters, SiStripBadStrip::put(), SiStripBadAPVAlgorithmFromClusterOccupancy::setAbsoluteLowThreshold(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setAbsoluteLowThreshold(), SiStripBadAPVAlgorithmFromClusterOccupancy::setAbsoluteOccupancyThreshold(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setAbsoluteOccupancyThreshold(), ConditionDBWriter< SiStripBadStrip >::setDoStore(), SiStripBadAPVAlgorithmFromClusterOccupancy::setHighOccupancyThreshold(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setHighOccupancyThreshold(), SiStripBadAPVAlgorithmFromClusterOccupancy::setLowOccupancyThreshold(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setLowOccupancyThreshold(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setMinNumEntries(), SiStripHotStripAlgorithmFromClusterOccupancy::setMinNumEntries(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setMinNumEntriesPerStrip(), SiStripHotStripAlgorithmFromClusterOccupancy::setMinNumEntriesPerStrip(), SiStripBadAPVAlgorithmFromClusterOccupancy::setMinNumOfEvents(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setMinNumOfEvents(), SiStripBadAPVAlgorithmFromClusterOccupancy::setNumberIterations(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setNumberIterations(), SiStripBadAPVAlgorithmFromClusterOccupancy::setNumberOfEvents(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setNumberOfEvents(), SiStripHotStripAlgorithmFromClusterOccupancy::setNumberOfEvents(), SiStripHotStripAlgorithmFromClusterOccupancy::setOccupancyThreshold(), SiStripBadAPVAlgorithmFromClusterOccupancy::setOutputFileName(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setOutputFileName(), SiStripHotStripAlgorithmFromClusterOccupancy::setOutputFileName(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setProbabilityThreshold(), SiStripHotStripAlgorithmFromClusterOccupancy::setProbabilityThreshold(), SiStripBadAPVAlgorithmFromClusterOccupancy::setTrackerGeometry(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setTrackerGeometry(), SiStripHotStripAlgorithmFromClusterOccupancy::setTrackerGeometry(), AlCaHarvesting_cff::SiStripQuality, SiStripQuality_, contentValuesCheck::ss, AlCaHLTBitMon_QueryRunRegistry::string, theIdentifier, theIdentifier2, theIdentifier3, TotNumberOfEvents, and tTopo.

52  {
53 
54  edm::LogInfo("SiStripQualityHotStripIdentifierRoot") <<"SiStripQualityHotStripIdentifierRoot::getNewObject called"<<std::endl;
55 
56  SiStripBadStrip* obj=new SiStripBadStrip();
57 
58  edm::ParameterSet parameters=conf_.getParameter<edm::ParameterSet>("AlgoParameters");
59  std::string AlgoName = parameters.getParameter<std::string>("AlgoName");
60 
61  bookHistos();
62 
63  if (TotNumberOfEvents>=calibrationthreshold)
64  {
65  edm::LogInfo("SiStripQualityHotStripIdentifierRoot") <<" [SiStripQualityHotStripIdentifierRoot::getNewObject] Total number of events is " << TotNumberOfEvents << " which is bigger than "<<calibrationthreshold<<". Calibration is preparing to be launched." <<std::endl;
66  if (MeanNumberOfCluster > 0.01)
67  {
68  edm::LogInfo("SiStripQualityHotStripIdentifierRoot") <<"Mean Number of Cluster is "<<MeanNumberOfCluster<<". Calibration is launched." <<std::endl;
69 
70  if (AlgoName=="SiStripHotStripAlgorithmFromClusterOccupancy")
71  {
72 
73  edm::LogInfo("SiStripQualityHotStripIdentifierRoot") <<" [SiStripQualityHotStripIdentifierRoot::getNewObject] call to SiStripHotStripAlgorithmFromClusterOccupancy"<<std::endl;
74 
75  theIdentifier= new SiStripHotStripAlgorithmFromClusterOccupancy(conf_, tTopo);
76  theIdentifier->setProbabilityThreshold(parameters.getUntrackedParameter<double>("ProbabilityThreshold",1.E-7));
77  theIdentifier->setMinNumEntries(parameters.getUntrackedParameter<uint32_t>("MinNumEntries",100));
78  theIdentifier->setMinNumEntriesPerStrip(parameters.getUntrackedParameter<uint32_t>("MinNumEntriesPerStrip",5));
79  theIdentifier->setOccupancyThreshold(parameters.getUntrackedParameter<double>("OccupancyThreshold",1.E-5));
80  theIdentifier->setNumberOfEvents(TotNumberOfEvents);
81  theIdentifier->setOutputFileName(conf_.getUntrackedParameter<std::string>("OccupancyRootFile","Occupancy.root"),conf_.getUntrackedParameter<bool>("WriteOccupancyRootFile",false));
82  theIdentifier->setTrackerGeometry(_tracker);
83 
84  SiStripQuality* qobj = new SiStripQuality();
85  theIdentifier->extractBadStrips(qobj,ClusterPositionHistoMap,SiStripQuality_);//here I insert SiStripQuality as input and get qobj as output
86 
87  //----------
88 
89  edm::LogInfo("SiStripQualityHotStripIdentifierRoot") <<" [SiStripQualityHotStripIdentifierRoot::getNewObject] copy SiStripObject in SiStripBadStrip"<<std::endl;
90 
91  std::stringstream ss;
92 
93  SiStripBadStrip::RegistryIterator rIter=qobj->getRegistryVectorBegin();
94  SiStripBadStrip::RegistryIterator rIterEnd=qobj->getRegistryVectorEnd();
95  for(;rIter!=rIterEnd;++rIter){
96  SiStripBadStrip::Range range(qobj->getDataVectorBegin()+rIter->ibegin,qobj->getDataVectorBegin()+rIter->iend);
97  if ( ! obj->put(rIter->detid,range) )
98  edm::LogError("SiStripQualityHotStripIdentifierRoot")<<"[SiStripQualityHotStripIdentifierRoot::getNewObject] detid already exists"<<std::endl;
99  }
100  edm::LogInfo("SiStripQualityHotStripIdentifierRoot") <<" [SiStripQualityHotStripIdentifierRoot::getNewObject] " << ss.str() << std::endl;
101  }
102  else if (AlgoName=="SiStripBadAPVAlgorithmFromClusterOccupancy")
103  {
104 
105  edm::LogInfo("SiStripQualityHotStripIdentifierRoot") <<" [SiStripQualityHotStripIdentifierRoot::getNewObject] call to SiStripBadAPVAlgorithmFromClusterOccupancy"<<std::endl;
106 
107  theIdentifier2 = new SiStripBadAPVAlgorithmFromClusterOccupancy(conf_, tTopo);
108  theIdentifier2->setLowOccupancyThreshold(parameters.getUntrackedParameter<double>("LowOccupancyThreshold",5));
109  theIdentifier2->setHighOccupancyThreshold(parameters.getUntrackedParameter<double>("HighOccupancyThreshold",10));
110  theIdentifier2->setAbsoluteLowThreshold(parameters.getUntrackedParameter<double>("AbsoluteLowThreshold",0));
111  theIdentifier2->setNumberIterations(parameters.getUntrackedParameter<uint32_t>("NumberIterations",2));
112  theIdentifier2->setAbsoluteOccupancyThreshold(parameters.getUntrackedParameter<double>("OccupancyThreshold",1.E-5));
113  theIdentifier2->setNumberOfEvents(TotNumberOfEvents);
114  theIdentifier2->setMinNumOfEvents();
115  theIdentifier2->setOutputFileName(conf_.getUntrackedParameter<std::string>("OccupancyRootFile","Occupancy.root"),conf_.getUntrackedParameter<bool>("WriteOccupancyRootFile",false));
116  theIdentifier2->setTrackerGeometry(_tracker);
117 
118  SiStripQuality* qobj = new SiStripQuality();
119  theIdentifier2->extractBadAPVs(qobj,ClusterPositionHistoMap,SiStripQuality_);
120 
121  //----------
122 
123  edm::LogInfo("SiStripQualityHotStripIdentifierRoot") <<" [SiStripQualityHotStripIdentifierRoot::getNewObject] copy SiStripObject in SiStripBadStrip"<<std::endl;
124 
125  std::stringstream ss;
126 
127  SiStripBadStrip::RegistryIterator rIter=qobj->getRegistryVectorBegin();
128  SiStripBadStrip::RegistryIterator rIterEnd=qobj->getRegistryVectorEnd();
129  for(;rIter!=rIterEnd;++rIter){
130  SiStripBadStrip::Range range(qobj->getDataVectorBegin()+rIter->ibegin,qobj->getDataVectorBegin()+rIter->iend);
131  if ( ! obj->put(rIter->detid,range) )
132  edm::LogError("SiStripQualityHotStripIdentifierRoot")<<"[SiStripQualityHotStripIdentifierRoot::getNewObject] detid already exists"<<std::endl;
133  }
134  edm::LogInfo("SiStripQualityHotStripIdentifierRoot") <<" [SiStripQualityHotStripIdentifierRoot::getNewObject] " << ss.str() << std::endl;
135 
136 
137  }
138  else if (AlgoName=="SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy")
139  {
140 
141  edm::LogInfo("SiStripQualityHotStripIdentifierRoot") <<" [SiStripQualityHotStripIdentifierRoot::getNewObject] call to SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy"<<std::endl;
142 
143  theIdentifier3= new SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy(conf_, tTopo);
144  theIdentifier3->setProbabilityThreshold(parameters.getUntrackedParameter<double>("ProbabilityThreshold",1.E-7));
145  theIdentifier3->setMinNumEntries(parameters.getUntrackedParameter<uint32_t>("MinNumEntries",100));
146  theIdentifier3->setMinNumEntriesPerStrip(parameters.getUntrackedParameter<uint32_t>("MinNumEntriesPerStrip",5));
147  theIdentifier3->setNumberOfEvents(TotNumberOfEvents);
148  theIdentifier3->setOutputFileName(conf_.getUntrackedParameter<std::string>("OccupancyRootFile","Occupancy.root"),conf_.getUntrackedParameter<bool>("WriteOccupancyRootFile",false),conf_.getUntrackedParameter<std::string>("DQMHistoOutputFile","DQMHistos.root"),conf_.getUntrackedParameter<bool>("WriteDQMHistoOutputFile",false));
149  theIdentifier3->setTrackerGeometry(_tracker);
150  theIdentifier3->setLowOccupancyThreshold(parameters.getUntrackedParameter<double>("LowOccupancyThreshold",5));
151  theIdentifier3->setHighOccupancyThreshold(parameters.getUntrackedParameter<double>("HighOccupancyThreshold",10));
152  theIdentifier3->setAbsoluteLowThreshold(parameters.getUntrackedParameter<double>("AbsoluteLowThreshold",0));
153  theIdentifier3->setNumberIterations(parameters.getUntrackedParameter<uint32_t>("NumberIterations",2));
154  theIdentifier3->setAbsoluteOccupancyThreshold(parameters.getUntrackedParameter<double>("OccupancyThreshold",1.E-5));
155  theIdentifier3->setMinNumOfEvents();
156 
157  SiStripQuality* qobj = new SiStripQuality();
158  theIdentifier3->extractBadAPVSandStrips(qobj,ClusterPositionHistoMap,SiStripQuality_);//here I insert SiStripQuality as input and get qobj as output
159 
160  //----------
161 
162  edm::LogInfo("SiStripQualityHotStripIdentifierRoot") <<" [SiStripQualityHotStripIdentifierRoot::getNewObject] copy SiStripObject in SiStripBadStrip"<<std::endl;
163 
164  std::stringstream ss;
165 
166  SiStripBadStrip::RegistryIterator rIter=qobj->getRegistryVectorBegin();
167  SiStripBadStrip::RegistryIterator rIterEnd=qobj->getRegistryVectorEnd();
168  for(;rIter!=rIterEnd;++rIter){
169  SiStripBadStrip::Range range(qobj->getDataVectorBegin()+rIter->ibegin,qobj->getDataVectorBegin()+rIter->iend);
170  if ( ! obj->put(rIter->detid,range) )
171  edm::LogError("SiStripQualityHotStripIdentifierRoot")<<"[SiStripQualityHotStripIdentifierRoot::getNewObject] detid already exists"<<std::endl;
172  }
173  edm::LogInfo("SiStripQualityHotStripIdentifierRoot") <<" [SiStripQualityHotStripIdentifierRoot::getNewObject] " << ss.str() << std::endl;
174  }
175  else
176  {
177  edm::LogError("SiStripQualityHotStripIdentifierRoot") <<" [SiStripQualityHotStripIdentifierRoot::getNewObject] call for a unknow HotStrip identification algoritm"<<std::endl;
178 
179  std::vector<uint32_t> a;
180  SiStripBadStrip::Range range(a.begin(),a.end());
181  if ( ! obj->put(0xFFFFFFFF,range) )
182  edm::LogError("SiStripQualityHotStripIdentifierRoot")<<"[SiStripQualityHotStripIdentifierRoot::getNewObject] detid already exists"<<std::endl;
183  }
184  }
185 
186  else
187  {
188  edm::LogWarning("SiStripQualityHotStripIdentifierRoot") <<" [SiStripQualityHotStripIdentifierRoot::getNewObject] :: Empty Events :: " <<"Mean Number of Cluster is " << MeanNumberOfCluster << " Calibration NOT launched" << std::endl;
189  setDoStore(false); // Don't put anything in the sqlite-file!
190  }
191  }
192  else
193  {
194  edm::LogWarning("SiStripQualityHotStripIdentifierRoot") <<" [SiStripQualityHotStripIdentifierRoot::getNewObject] :: Not Enough Events :: " <<"Total number of events is " << TotNumberOfEvents << " Calibration NOT launched" << std::endl;
195  setDoStore(false); // Don't put anything in the sqlite-file!
196  }
197 
198  return obj;
199 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
SiStripHotStripAlgorithmFromClusterOccupancy::setOutputFileName
void setOutputFileName(std::string OutputFileName, bool WriteOutputFile)
Definition: SiStripHotStripAlgorithmFromClusterOccupancy.h:67
SiStripHotStripAlgorithmFromClusterOccupancy::extractBadStrips
void extractBadStrips(SiStripQuality *, HistoMap &, edm::ESHandle< SiStripQuality > &)
Definition: SiStripHotStripAlgorithmFromClusterOccupancy.cc:26
SiStripQualityHotStripIdentifierRoot::theIdentifier3
SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy * theIdentifier3
Definition: SiStripQualityHotStripIdentifierRoot.h:97
SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setOutputFileName
void setOutputFileName(std::string OutputFileName, bool WriteOutputFile, std::string DQMOutfileName, bool WriteDQMHistograms)
Definition: SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy.h:63
HLT_25ns10e33_v2_cff.AlgoName
tuple AlgoName
Definition: HLT_25ns10e33_v2_cff.py:2729
SiStripBadStrip::RegistryIterator
Registry::const_iterator RegistryIterator
Definition: SiStripBadStrip.h:56
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
AlCaHarvesting_cff.SiStripQuality
tuple SiStripQuality
Definition: AlCaHarvesting_cff.py:90
SiStripBadStrip::getRegistryVectorEnd
RegistryIterator getRegistryVectorEnd() const
Definition: SiStripBadStrip.h:77
SiStripQualityHotStripIdentifierRoot::theIdentifier
SiStripHotStripAlgorithmFromClusterOccupancy * theIdentifier
Definition: SiStripQualityHotStripIdentifierRoot.h:95
HLT_25ns10e33_v2_cff.parameters
tuple parameters
Definition: HLT_25ns10e33_v2_cff.py:3864
SiStripBadStrip::getDataVectorBegin
ContainerIterator getDataVectorBegin() const
Definition: SiStripBadStrip.h:74
SiStripBadAPVAlgorithmFromClusterOccupancy::setOutputFileName
void setOutputFileName(std::string OutputFileName, bool WriteOutputFile)
Definition: SiStripBadAPVAlgorithmFromClusterOccupancy.h:57
SiStripBadStrip::getRegistryVectorBegin
RegistryIterator getRegistryVectorBegin() const
Definition: SiStripBadStrip.h:76
ConditionDBWriter< SiStripBadStrip >::setDoStore
void setDoStore(const bool doStore)
When set to false the payload will not be written to the db.
Definition: ConditionDBWriter.h:379
a
double a
Definition: hdecay.h:121
SiStripBadStrip::Range
std::pair< ContainerIterator, ContainerIterator > Range
Definition: SiStripBadStrip.h:54
SiStripBadAPVAlgorithmFromClusterOccupancy::extractBadAPVs
void extractBadAPVs(SiStripQuality *, HistoMap &, edm::ESHandle< SiStripQuality > &)
Definition: SiStripBadAPVAlgorithmFromClusterOccupancy.cc:30
SiStripQualityHotStripIdentifierRoot::theIdentifier2
SiStripBadAPVAlgorithmFromClusterOccupancy * theIdentifier2
Definition: SiStripQualityHotStripIdentifierRoot.h:96
SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::extractBadAPVSandStrips
void extractBadAPVSandStrips(SiStripQuality *, HistoMap &, edm::ESHandle< SiStripQuality > &)
Definition: SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy.cc:32
SiStripBadStrip::put
bool put(const uint32_t &detID, const InputVector &vect)
Definition: SiStripBadStrip.h:66

Member Data Documentation

const TrackerGeometry* SiStripQualityHotStripIdentifierRoot::_tracker
private

Definition at line 82 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by algoBeginRun(), and getNewObject().

uint32_t SiStripQualityHotStripIdentifierRoot::calibrationthreshold
private

Definition at line 92 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by getNewObject().

SiStrip::QualityHistosMap SiStripQualityHotStripIdentifierRoot::ClusterPositionHistoMap
private

Definition at line 94 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by algoEndJob(), bookHistos(), and getNewObject().

const edm::ParameterSet SiStripQualityHotStripIdentifierRoot::conf_
private

Definition at line 77 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by bookHistos(), and getNewObject().

std::string SiStripQualityHotStripIdentifierRoot::dataLabel_
private

Definition at line 74 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by algoBeginRun().

std::string SiStripQualityHotStripIdentifierRoot::dirpath
private

Definition at line 88 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by bookHistos().

DQMStore* SiStripQualityHotStripIdentifierRoot::dqmStore_
private

Definition at line 85 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by bookHistos(), and SiStripQualityHotStripIdentifierRoot().

TFile* SiStripQualityHotStripIdentifierRoot::file0
private

Definition at line 87 of file SiStripQualityHotStripIdentifierRoot.h.

std::string SiStripQualityHotStripIdentifierRoot::filename
private

Definition at line 88 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by cuy.ValElement::__init__(), python.rootplot.rootmath.Target::__repr__(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), utils.unpickler::run(), and SiStripQualityHotStripIdentifierRoot().

edm::FileInPath SiStripQualityHotStripIdentifierRoot::fp_
private

Definition at line 78 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by SiStripQualityHotStripIdentifierRoot().

unsigned long long SiStripQualityHotStripIdentifierRoot::m_cacheID_
private

Definition at line 73 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by algoBeginRun().

unsigned short SiStripQualityHotStripIdentifierRoot::MaxClusterWidth_
private

Definition at line 89 of file SiStripQualityHotStripIdentifierRoot.h.

double SiStripQualityHotStripIdentifierRoot::MeanNumberOfCluster
private

Definition at line 91 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by bookHistos(), and getNewObject().

unsigned short SiStripQualityHotStripIdentifierRoot::MinClusterWidth_
private

Definition at line 89 of file SiStripQualityHotStripIdentifierRoot.h.

SiStripDetInfoFileReader* SiStripQualityHotStripIdentifierRoot::reader
private

Definition at line 79 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by SiStripQualityHotStripIdentifierRoot().

edm::ESHandle<SiStripQuality> SiStripQualityHotStripIdentifierRoot::SiStripQuality_
private

Definition at line 75 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by algoBeginRun(), and getNewObject().

SiStripHotStripAlgorithmFromClusterOccupancy* SiStripQualityHotStripIdentifierRoot::theIdentifier
private

Definition at line 95 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by getNewObject().

SiStripBadAPVAlgorithmFromClusterOccupancy* SiStripQualityHotStripIdentifierRoot::theIdentifier2
private

Definition at line 96 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by getNewObject().

SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy* SiStripQualityHotStripIdentifierRoot::theIdentifier3
private

Definition at line 97 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by getNewObject().

edm::ESHandle<TrackerGeometry> SiStripQualityHotStripIdentifierRoot::theTrackerGeom
private

Definition at line 81 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by algoBeginRun().

double SiStripQualityHotStripIdentifierRoot::TotNumberOfEvents
private

Definition at line 90 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by bookHistos(), and getNewObject().

const TrackerTopology* SiStripQualityHotStripIdentifierRoot::tTopo
private

Definition at line 83 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by algoBeginRun(), and getNewObject().

bool SiStripQualityHotStripIdentifierRoot::UseInputDB_
private

Definition at line 76 of file SiStripQualityHotStripIdentifierRoot.h.

Referenced by algoBeginRun().

Generated for CMSSW Reference Manual by doxygen 1.8.5