41 edm::LogInfo (
"PixelHLTDQM") <<
"SiPixelMuonHLT::SiPixelMuonHLT: Got DQM BackEnd interface"<<std::endl;
48 if (theDMBE !=
NULL) theDMBE->setCurrentFolder (monitorName_);
57 edm::LogInfo (
"PixelHLTDQM") <<
"SiPixelMuonHLT::~SiPixelMuonHLT: Destructor"<<std::endl;
63 edm::LogInfo (
"PixelHLTDQM") <<
" SiPixelMuonHLT::beginJob - Initialisation ... " << std::endl;
71 edm::LogInfo (
"PixelHLTDQM") <<
" SiPixelMuonHLT::endJob - Saving Root File " << std::endl;
88 reco::RecoChargedCandidateCollection::const_iterator cand;
92 bool GotClusters =
true;
93 bool GotRecHits =
true;
94 bool GotL3Muons =
true;
96 iEvent.
getByLabel(
"hltSiPixelClusters", clusters);
101 iEvent.
getByLabel(
"hltSiPixelRecHits", rechits);
115 int NBarrel[4] = {0,0,0,0};
116 int NEndcap[5] = {0,0,0,0,0};
117 for (
size_t i = 0;
i < clusters->size(); ++
i){
121 uint detID = clusters->
id(
i);
149 clustDisk = clustDisk +2;
155 ++NEndcap[clustDisk];
160 for (
int lay = 1; lay < 4; ++lay)
163 for (
int disk = 1; disk < 5; ++disk)
167 bool doRecHits =
false;
169 if (GotRecHits && doRecHits){
172 for (
size_t i = 0;
i < rechits->size(); ++
i){
174 uint detID = rechits->
id(
i);
182 edm::LogInfo(
"PixelHLTDQM") <<
" (From SiPixelRecHit) Hit Eta: " << hitgp.
eta() <<
" Hit Phi: " << hitgp.
phi() << std::endl;
190 int NBarrel[4] = {0,0,0,0};
191 int NEndcap[5] = {0,0,0,0,0};
192 for (cand = l3mucands->begin (); cand != l3mucands->end (); ++cand){
194 for (
size_t hit = 0;
hit < l3tk->recHitsSize ();
hit++){
195 if (l3tk->recHit (
hit)->isValid () ==
true && l3tk->recHit (
hit)->geographicalId ().det () ==
DetId::Tracker){
196 int detID = l3tk->recHit(
hit)->geographicalId().rawId();
201 if((*pixhit).isValid() ==
true){
203 if (!(*pixhit).cluster().isAvailable())
209 if(l3tk->recHit(
hit)->geographicalId().subdetId() == 1){
225 if(l3tk->recHit(
hit)->geographicalId().subdetId() == 2){
228 clustDisk = clustDisk +2;
234 ++NEndcap[clustDisk];
242 for (
int lay = 1; lay < 4; ++lay)
245 for (
int disk = 1; disk < 5; ++disk)
260 float PhiMax = 3.142;
267 std::string histoname;
271 std::string layerLabel[4] = {
"All_Layers",
"Layer1",
"Layer2",
"Layer3"};
272 for (
unsigned int i = 0;
i < 4;
i++)
277 histoname =
"EtaPhiAllBarrelMap_" + layerLabel[
i];
278 title =
"#eta-#phi Pixel Barrel Cluster Occupancy Map for " + layerLabel[
i];
281 histoname =
"ZPhiAllBarrelMap_" + layerLabel[
i];
282 title =
"Z-#phi Pixel Barrel Cluster Occupancy Map for " + layerLabel[
i];
285 histoname =
"EtaAllBarrelMap_" + layerLabel[
i];
286 title =
"#eta Pixel Barrel Cluster Occupancy Map for " + layerLabel[
i];
289 histoname =
"ZAllBarrelMap_" + layerLabel[
i];
290 title =
"Z Pixel Barrel Cluster Occupancy Map for " + layerLabel[
i];
293 histoname =
"PhiAllBarrelMap_" + layerLabel[
i];
294 title =
"#phi Pixel Barrel Cluster Occupancy Map for " + layerLabel[
i];
297 histoname =
"NAllBarrelMap_" + layerLabel[
i];
298 title =
"#phi Pixel Barrel Cluster Occupancy Map for " + layerLabel[
i];
303 histoname =
"EtaPhiOnTrackBarrelMap_" + layerLabel[
i];
304 title =
"#eta-#phi On Track Pixel Barrel Cluster Occupancy Map for " + layerLabel[
i];
307 histoname =
"ZPhiOnTrackBarrelMap_" + layerLabel[
i];
308 title =
"Z-#phi On Track Pixel Barrel Cluster Occupancy Map for " + layerLabel[
i];
311 histoname =
"EtaOnTrackBarrelMap_" + layerLabel[
i];
312 title =
"#eta On Track Pixel Barrel Cluster Occupancy Map for " + layerLabel[
i];
315 histoname =
"ZOnTrackBarrelMap_" + layerLabel[
i];
316 title =
"Z On Track Pixel Barrel Cluster Occupancy Map for " + layerLabel[
i];
319 histoname =
"PhiOnTrackBarrelMap_" + layerLabel[
i];
320 title =
"#phi On Track Pixel Barrel Cluster Occupancy Map for " + layerLabel[
i];
323 histoname =
"NOnTrackBarrelMap_" + layerLabel[
i];
324 title =
"N_{Clusters} On Track Pixel Barrel Cluster Occupancy Map for " + layerLabel[
i];
329 std::string diskLabel[5] = {
"All_Disks",
"InnerPosZ",
"OuterPosZ",
"InnerNegZ",
"OuterNegZ"};
330 for (
int i = 0;
i < 5; ++
i)
335 histoname =
"XYAllEndcapMap_" + diskLabel[
i];
336 title =
"X-Y Pixel Endcap Cluster Occupancy Map for " + diskLabel[
i];
339 histoname =
"PhiAllEndcapMap_" + diskLabel[
i];
340 title =
"#phi Pixel Endcap Cluster Occupancy Map for " + diskLabel[
i];
343 histoname =
"NAllEndcapMap_" + diskLabel[
i];
344 title =
"#phi Pixel Endcap Cluster Occupancy Map for " + diskLabel[
i];
349 histoname =
"XYOnTrackEndcapMap_" + diskLabel[
i];
350 title =
"X-Y Pixel Endcap On Track Cluster Occupancy Map for " + diskLabel[
i];
353 histoname =
"PhiOnTrackEndcapMap_" + diskLabel[
i];
354 title =
"#phi Pixel Endcap On Track Cluster Occupancy Map for " + diskLabel[
i];
357 histoname =
"NOnTrackEndcapMap_" + diskLabel[
i];
358 title =
"#phi Pixel Endcap On Track Cluster Occupancy Map for " + diskLabel[
i];
SiPixelMuonHLT(const edm::ParameterSet &conf)
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
edm::InputTag rechitsCollectionTag_
T getUntrackedParameter(std::string const &, T const &) const
std::map< int, MonitorElement * > MEContainerAllBarrelZ
edm::InputTag l3MuonCollectionTag_
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
std::map< int, MonitorElement * > MEContainerAllEndcapPhi
std::map< int, MonitorElement * > MEContainerAllBarrelN
#define DEFINE_FWK_MODULE(type)
bool verbose_
every n events
void save(const std::string &filename, const std::string &path="", const std::string &pattern="", const std::string &rewrite="", SaveReferenceTag ref=SaveWithReference, int minStatus=dqm::qstatus::STATUS_OK, const std::string &fileupdate="RECREATE")
std::map< int, MonitorElement * > MEContainerOnTrackBarrelPhi
Geom::Phi< T > phi() const
std::map< int, MonitorElement * > MEContainerAllBarrelZPhi
std::map< int, MonitorElement * > MEContainerOnTrackEndcapN
edm::ParameterSet parameters_
const double EtaMax[kNumberCalorimeter]
std::map< int, MonitorElement * > MEContainerAllBarrelEta
std::map< int, MonitorElement * > MEContainerOnTrackBarrelN
std::map< int, MonitorElement * > MEContainerAllBarrelPhi
virtual LocalPoint localPosition() const
virtual void analyze(const edm::Event &, const edm::EventSetup &)
unsigned int layer() const
layer id
std::map< int, MonitorElement * > MEContainerAllEndcapXY
std::map< int, MonitorElement * > MEContainerOnTrackBarrelEta
Measurement2DPoint MeasurementPoint
Measurement points are two-dimensional by default.
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
bool hasPositionAndError() const
edm::InputTag clusterCollectionTag_
std::map< int, MonitorElement * > MEContainerAllEndcapN
DetId geographicalId() const
The label of this GeomDet.
std::map< int, MonitorElement * > MEContainerAllBarrelEtaPhi
int subdetId() const
get the contents of the subdetector field (not cast into any detector's numbering enum) ...
virtual const GeomDet * idToDet(DetId) const
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
unsigned int disk() const
disk id
EventAuxiliary const & eventAuxiliary() const
virtual const PixelTopology & specificTopology() const
Returns a reference to the pixel proxy topology.
T const * product() const
std::map< int, MonitorElement * > MEContainerOnTrackBarrelZPhi
Pixel cluster – collection of neighboring pixels above threshold.
std::map< int, MonitorElement * > MEContainerOnTrackBarrelZ
virtual void Histo_init()
std::map< int, MonitorElement * > MEContainerOnTrackBarrelEtaPhi
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
T const * get() const
Returns C++ pointer to the item.
virtual const BoundPlane & surface() const
The nominal surface of the GeomDet.
EventNumber_t event() const
void setCurrentFolder(const std::string &fullpath)
std::map< int, MonitorElement * > MEContainerOnTrackEndcapXY
std::map< int, MonitorElement * > MEContainerOnTrackEndcapPhi