Go to the documentation of this file.00001
00002
00003
00004
00005
00013
00014
00015
00016
00017
00018
00019
00020
00021 #include <memory>
00022
00023
00024 #include <vector>
00025 #include <string>
00026
00027 #include "FWCore/Framework/interface/Frameworkfwd.h"
00028 #include "FWCore/Framework/interface/EDAnalyzer.h"
00029
00030 #include "FWCore/Framework/interface/Event.h"
00031 #include "FWCore/Framework/interface/MakerMacros.h"
00032
00033 #include "FWCore/ParameterSet/interface/ParameterSet.h"
00034
00035 #include "FWCore/MessageLogger/interface/MessageLogger.h"
00036
00037 #include "FWCore/ServiceRegistry/interface/Service.h"
00038 #include "CommonTools/UtilAlgos/interface/TFileService.h"
00039
00040 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
00041 #include "Geometry/TrackerGeometryBuilder/interface/PixelGeomDetUnit.h"
00042 #include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
00043
00044 #include "CalibTracker/SiStripCommon/interface/SiStripDetInfoFileReader.h"
00045 #include "CalibTracker/SiPixelESProducers/interface/SiPixelDetInfoFileReader.h"
00046 #include "CommonTools/UtilAlgos/interface/DetIdSelector.h"
00047 #include "FWCore/ParameterSet/interface/FileInPath.h"
00048
00049
00050 #include "DQM/SiStripCommon/interface/TkHistoMap.h"
00051 #include "CommonTools/TrackerMap/interface/TrackerMap.h"
00052
00053
00054
00055
00056
00057
00058 class DetIdSelectorTest : public edm::EDAnalyzer {
00059 public:
00060 explicit DetIdSelectorTest(const edm::ParameterSet&);
00061 ~DetIdSelectorTest();
00062
00063
00064 private:
00065 virtual void beginJob() ;
00066 virtual void analyze(const edm::Event&, const edm::EventSetup&) override;
00067 virtual void endJob() ;
00068
00069
00070
00071
00072 std::vector<DetIdSelector> detidsels_;
00073 TkHistoMap *tkhisto_;
00074 TrackerMap tkmap_;
00075
00076
00077 };
00078
00079
00080
00081
00082
00083
00084
00085
00086
00087
00088
00089
00090 DetIdSelectorTest::DetIdSelectorTest(const edm::ParameterSet& iConfig):
00091 detidsels_(), tkhisto_(new TkHistoMap("SelectorTest","SelectorTest",-1)),tkmap_()
00092 {
00093
00094
00095 std::vector<edm::ParameterSet> selconfigs = iConfig.getParameter<std::vector<edm::ParameterSet> >("selections");
00096
00097 for(std::vector<edm::ParameterSet>::const_iterator selconfig=selconfigs.begin();selconfig!=selconfigs.end();++selconfig) {
00098 DetIdSelector selection(*selconfig);
00099 detidsels_.push_back(selection);
00100 }
00101
00102 tkmap_.setPalette(1);
00103 tkmap_.addPixel(true);
00104
00105 }
00106
00107
00108 DetIdSelectorTest::~DetIdSelectorTest()
00109 {
00110
00111
00112
00113
00114 }
00115
00116
00117
00118
00119
00120
00121
00122 void
00123 DetIdSelectorTest::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
00124 {
00125 using namespace edm;
00126
00127 {
00128 SiStripDetInfoFileReader * reader=edm::Service<SiStripDetInfoFileReader>().operator->();
00129
00130
00131
00132
00133 const std::vector<uint32_t>& detids = reader->getAllDetIds();
00134
00135 for(std::vector<uint32_t>::const_iterator detid=detids.begin();detid!=detids.end();++detid) {
00136
00137 LogDebug("DetID") << *detid;
00138 int index=0;
00139 for(std::vector<DetIdSelector>::const_iterator detidsel=detidsels_.begin();detidsel!=detidsels_.end();++detidsel) {
00140 if(detidsel->isSelected(*detid)) {
00141 LogDebug("selected") << "Selected by selection " << index;
00142 unsigned int det = *detid;
00143 tkhisto_->add(det,index);
00144 tkmap_.fill_current_val(det,index);
00145 }
00146 ++index;
00147 }
00148
00149 }
00150 }
00151
00152 {
00153 edm::FileInPath fp("CalibTracker/SiPixelESProducers/data/PixelSkimmedGeometry.txt");
00154
00155 SiPixelDetInfoFileReader pxlreader(fp.fullPath());
00156 const std::vector<uint32_t>& detids = pxlreader.getAllDetIds();
00157
00158 for(std::vector<uint32_t>::const_iterator detid=detids.begin();detid!=detids.end();++detid) {
00159
00160 LogDebug("DetID") << *detid;
00161 int index=0;
00162 for(std::vector<DetIdSelector>::const_iterator detidsel=detidsels_.begin();detidsel!=detidsels_.end();++detidsel) {
00163 if(detidsel->isSelected(*detid)) {
00164 LogDebug("selected") << "Selected by selection " << index;
00165 unsigned int det = *detid;
00166
00167 tkmap_.fill_current_val(det,index);
00168 }
00169 ++index;
00170 }
00171
00172 }
00173 }
00174
00175
00176
00177
00178
00179
00180
00181
00182
00183
00184
00185
00186
00187
00188
00189
00190
00191
00192
00193
00194
00195
00196
00197
00198
00199 }
00200
00201
00202 void
00203 DetIdSelectorTest::beginJob()
00204 {
00205
00206 }
00207
00208
00209 void
00210 DetIdSelectorTest::endJob() {
00211
00212
00213
00214 std::string mapname = "SelectorTest.png";
00215 tkmap_.save(true,0,0,mapname,5700,2400);
00216
00217 std::string rootmapname = "TKMap_Selectortest.root";
00218 tkhisto_->save(rootmapname);
00219 }
00220
00221
00222
00223 DEFINE_FWK_MODULE(DetIdSelectorTest);