Reference Manual

Go to the documentation of this file.
00001 // -*- C++ -*-
00002 // Package:    SiStripMonitorCluster
00003 // Class:      SiStripMonitorCluster
00006 // Original Author:  Dorian Kcira
00007 //         Created:  Wed Feb  1 16:42:34 CET 2006
00008 // $Id: SiStripMonitorCluster.cc,v 1.84 2012/11/20 14:45:33 eulisse Exp $
00009 #include <vector>
00010 #include <numeric>
00011 #include <fstream>
00012 #include <math.h>
00013 #include "TNamed.h"
00014 #include "FWCore/MessageLogger/interface/MessageLogger.h"
00015 #include "FWCore/ServiceRegistry/interface/Service.h"
00016 
00017 #include "CalibTracker/Records/interface/SiStripDetCablingRcd.h"
00018 #include "CalibFormats/SiStripObjects/interface/SiStripDetCabling.h"
00019 #include "CondFormats/DataRecord/interface/SiStripCondDataRecords.h"
00020 #include "CondFormats/SiStripObjects/interface/SiStripNoises.h"
00021 #include "CalibFormats/SiStripObjects/interface/SiStripGain.h"
00022 #include "CalibFormats/SiStripObjects/interface/SiStripQuality.h"
00023 #include "DataFormats/SiStripCluster/interface/SiStripClusterCollection.h"
00024 #include "DataFormats/Common/interface/DetSetVector.h"
00025 #include "DataFormats/Common/interface/DetSetVectorNew.h"
00026 #include "DataFormats/SiStripCluster/interface/SiStripCluster.h"
00027 #include "DQM/SiStripCommon/interface/SiStripFolderOrganizer.h"
00028 #include "DQM/SiStripCommon/interface/SiStripHistoId.h"
00029 #include "DQM/SiStripMonitorCluster/interface/SiStripMonitorCluster.h"
00030 #include "DQMServices/Core/interface/DQMStore.h"
00031 #include "DQMServices/Core/interface/MonitorElement.h"
00032 #include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
00033 #include "DataFormats/SiStripDetId/interface/SiStripSubStructure.h"
00034 #include "DataFormats/SiStripDetId/interface/TECDetId.h"
00035 #include "DataFormats/SiStripDetId/interface/TIBDetId.h"
00036 #include "DataFormats/SiStripDetId/interface/TIDDetId.h"
00037 #include "DataFormats/SiStripDetId/interface/TOBDetId.h"
00038 #include "CalibTracker/SiStripCommon/interface/SiStripDCSStatus.h"
00039 #include "DataFormats/SiPixelCluster/interface/SiPixelCluster.h"
00040 
00041 #include "DPGAnalysis/SiStripTools/interface/APVCyclePhaseCollection.h"
00042 #include "DPGAnalysis/SiStripTools/interface/EventWithHistory.h"
00043 
00044 #include "CommonTools/TriggerUtils/interface/GenericTriggerEventFlag.h"
00045 
00046 #include "TMath.h"
00047 #include <iostream>
00048 
00049 //--------------------------------------------------------------------------------------------
00050 SiStripMonitorCluster::SiStripMonitorCluster(const edm::ParameterSet& iConfig)
00051   : dqmStore_(edm::Service<DQMStore>().operator->()), conf_(iConfig), show_mechanical_structure_view(true), show_readout_view(false), show_control_view(false), select_all_detectors(false), reset_each_run(false), m_cacheID_(0)
00052                                             //  , genTriggerEventFlag_(new GenericTriggerEventFlag(iConfig))
00053 {
00054 
00055   // initialize
00056   passBPTXfilter_ = true;
00057 
00058   // initialize GenericTriggerEventFlag by specific configuration
00059   // in this way, one can set specific selections for different MEs
00060   genTriggerEventFlagBPTXfilter_     = new GenericTriggerEventFlag(iConfig.getParameter<edm::ParameterSet>("BPTXfilter")     );
00061   genTriggerEventFlagPixelDCSfilter_ = new GenericTriggerEventFlag(iConfig.getParameter<edm::ParameterSet>("PixelDCSfilter") );
00062   genTriggerEventFlagStripDCSfilter_ = new GenericTriggerEventFlag(iConfig.getParameter<edm::ParameterSet>("StripDCSfilter") );
00063 
00064   firstEvent = -1;
00065   eventNb = 0;
00066 
00067   // Detector Partitions
00068   SubDetPhasePartMap["TIB"]        = "TI";
00069   SubDetPhasePartMap["TID__side__1"] = "TI";
00070   SubDetPhasePartMap["TID__side__2"] = "TI";
00071   SubDetPhasePartMap["TOB"]        = "TO";
00072   SubDetPhasePartMap["TEC__side__1"] = "TM";
00073   SubDetPhasePartMap["TEC__side__2"] = "TP";
00074 
00075   //get on/off option for every cluster from cfi
00076   edm::ParameterSet ParametersnClusters =  conf_.getParameter<edm::ParameterSet>("TH1nClusters");
00077   layerswitchncluson = ParametersnClusters.getParameter<bool>("layerswitchon");
00078   moduleswitchncluson = ParametersnClusters.getParameter<bool>("moduleswitchon");
00079   
00080   edm::ParameterSet ParametersClusterCharge =  conf_.getParameter<edm::ParameterSet>("TH1ClusterCharge");
00081   layerswitchcluschargeon = ParametersClusterCharge.getParameter<bool>("layerswitchon");
00082   moduleswitchcluschargeon = ParametersClusterCharge.getParameter<bool>("moduleswitchon");
00083   
00084   edm::ParameterSet ParametersClusterStoN =  conf_.getParameter<edm::ParameterSet>("TH1ClusterStoN");
00085   layerswitchclusstonon = ParametersClusterStoN.getParameter<bool>("layerswitchon");
00086   moduleswitchclusstonon = ParametersClusterStoN.getParameter<bool>("moduleswitchon");
00087 
00088   edm::ParameterSet ParametersClusterStoNVsPos =  conf_.getParameter<edm::ParameterSet>("TH1ClusterStoNVsPos");
00089   layerswitchclusstonVsposon = ParametersClusterStoNVsPos.getParameter<bool>("layerswitchon");
00090   moduleswitchclusstonVsposon = ParametersClusterStoNVsPos.getParameter<bool>("moduleswitchon");
00091   
00092   edm::ParameterSet ParametersClusterPos =  conf_.getParameter<edm::ParameterSet>("TH1ClusterPos");
00093   layerswitchclusposon = ParametersClusterPos.getParameter<bool>("layerswitchon");
00094   moduleswitchclusposon = ParametersClusterPos.getParameter<bool>("moduleswitchon");
00095 
00096   edm::ParameterSet ParametersClusterDigiPos =  conf_.getParameter<edm::ParameterSet>("TH1ClusterDigiPos");
00097   layerswitchclusdigiposon = ParametersClusterDigiPos.getParameter<bool>("layerswitchon");
00098   moduleswitchclusdigiposon = ParametersClusterDigiPos.getParameter<bool>("moduleswitchon");
00099   
00100   edm::ParameterSet ParametersClusterNoise =  conf_.getParameter<edm::ParameterSet>("TH1ClusterNoise");
00101   layerswitchclusnoiseon = ParametersClusterNoise.getParameter<bool>("layerswitchon");
00102   moduleswitchclusnoiseon = ParametersClusterNoise.getParameter<bool>("moduleswitchon");
00103   
00104   edm::ParameterSet ParametersClusterWidth =  conf_.getParameter<edm::ParameterSet>("TH1ClusterWidth");
00105   layerswitchcluswidthon = ParametersClusterWidth.getParameter<bool>("layerswitchon");
00106   moduleswitchcluswidthon = ParametersClusterWidth.getParameter<bool>("moduleswitchon");
00107   
00108   edm::ParameterSet ParametersModuleLocalOccupancy =  conf_.getParameter<edm::ParameterSet>("TH1ModuleLocalOccupancy");
00109   layerswitchlocaloccupancy = ParametersModuleLocalOccupancy.getParameter<bool>("layerswitchon");
00110   moduleswitchlocaloccupancy = ParametersModuleLocalOccupancy.getParameter<bool>("moduleswitchon");
00111 
00112   edm::ParameterSet ParametersNrOfClusterizedStrips =  conf_.getParameter<edm::ParameterSet>("TH1NrOfClusterizedStrips");
00113   layerswitchnrclusterizedstrip = ParametersNrOfClusterizedStrips.getParameter<bool>("layerswitchon");
00114   moduleswitchnrclusterizedstrip = ParametersNrOfClusterizedStrips.getParameter<bool>("moduleswitchon");
00115 
00116   edm::ParameterSet ParametersClusterProf = conf_.getParameter<edm::ParameterSet>("TProfNumberOfCluster");
00117   layerswitchnumclusterprofon = ParametersClusterProf.getParameter<bool>("layerswitchon");
00118 
00119   edm::ParameterSet ParametersClusterWidthProf = conf_.getParameter<edm::ParameterSet>("TProfClusterWidth");
00120   layerswitchclusterwidthprofon = ParametersClusterWidthProf.getParameter<bool>("layerswitchon");
00121 
00122   edm::ParameterSet ParametersTotClusterProf = conf_.getParameter<edm::ParameterSet>("TProfTotalNumberOfClusters");
00123   subdetswitchtotclusprofon = ParametersTotClusterProf.getParameter<bool>("subdetswitchon");
00124 
00125   edm::ParameterSet ParametersTotClusterTH1 = conf_.getParameter<edm::ParameterSet>("TH1TotalNumberOfClusters");
00126   subdetswitchtotclusth1on = ParametersTotClusterTH1.getParameter<bool>("subdetswitchon");
00127 
00128   edm::ParameterSet ParametersClusterApvProf = conf_.getParameter<edm::ParameterSet>("TProfClustersApvCycle");
00129   subdetswitchapvcycleprofon = ParametersClusterApvProf.getParameter<bool>("subdetswitchon");
00130 
00131   edm::ParameterSet ParametersClustersApvTH2 = conf_.getParameter<edm::ParameterSet>("TH2ClustersApvCycle");
00132   subdetswitchapvcycleth2on = ParametersClustersApvTH2.getParameter<bool>("subdetswitchon");
00133 
00134   edm::ParameterSet ParametersApvCycleDBxProf2 = conf_.getParameter<edm::ParameterSet>("TProf2ApvCycleVsDBx");
00135   subdetswitchapvcycledbxprof2on = ParametersApvCycleDBxProf2.getParameter<bool>("subdetswitchon");
00136 
00137   edm::ParameterSet ParametersDBxCycleProf = conf_.getParameter<edm::ParameterSet>("TProfClustersVsDBxCycle");
00138   subdetswitchdbxcycleprofon = ParametersDBxCycleProf.getParameter<bool>("subdetswitchon");
00139 
00140   edm::ParameterSet ParametersCStripVsCPix = conf_.getParameter<edm::ParameterSet>("TH2CStripVsCpixel");
00141   globalswitchcstripvscpix = ParametersCStripVsCPix.getParameter<bool>("globalswitchon");
00142 
00143   edm::ParameterSet ParametersMultiplicityRegionsTH1 = conf_.getParameter<edm::ParameterSet>("TH1MultiplicityRegions");
00144   globalswitchMultiRegions =  ParametersMultiplicityRegionsTH1.getParameter<bool>("globalswitchon");
00145 
00146   edm::ParameterSet ParametersApvCycleVsDBxGlobalTH2 = conf_.getParameter<edm::ParameterSet>("TH2ApvCycleVsDBxGlobal");
00147   globalswitchapvcycledbxth2on = ParametersApvCycleVsDBxGlobalTH2.getParameter<bool>("globalswitchon");
00148 
00149   edm::ParameterSet ParametersNoiseStrip2ApvCycle = conf_.getParameter<edm::ParameterSet>("TH1StripNoise2ApvCycle");
00150   globalswitchstripnoise2apvcycle = ParametersNoiseStrip2ApvCycle.getParameter<bool>("globalswitchon");
00151 
00152   edm::ParameterSet ParametersNoiseStrip3ApvCycle = conf_.getParameter<edm::ParameterSet>("TH1StripNoise3ApvCycle");
00153   globalswitchstripnoise3apvcycle = ParametersNoiseStrip3ApvCycle.getParameter<bool>("globalswitchon");
00154 
00155   edm::ParameterSet ParametersMainDiagonalPosition = conf_.getParameter<edm::ParameterSet>("TH1MainDiagonalPosition");
00156   globalswitchmaindiagonalposition= ParametersMainDiagonalPosition.getParameter<bool>("globalswitchon");
00157 
00158   edm::ParameterSet ClusterMultiplicityRegions = conf_.getParameter<edm::ParameterSet>("MultiplicityRegions");
00159   k0 = ClusterMultiplicityRegions.getParameter<double>("k0");
00160   q0 = ClusterMultiplicityRegions.getParameter<double>("q0");
00161   dk0 = ClusterMultiplicityRegions.getParameter<double>("dk0");
00162   maxClus = ClusterMultiplicityRegions.getParameter<double>("MaxClus");
00163   minPix = ClusterMultiplicityRegions.getParameter<double>("MinPix");
00164 
00165   clustertkhistomapon = conf_.getParameter<bool>("TkHistoMap_On");
00166   createTrendMEs = conf_.getParameter<bool>("CreateTrendMEs");
00167   Mod_On_ = conf_.getParameter<bool>("Mod_On");
00168   ClusterHisto_ = conf_.getParameter<bool>("ClusterHisto");
00169 
00170   topFolderName_ = conf_.getParameter<std::string>("TopFolderName");
00171 
00172 
00173   // Poducer name of input StripClusterCollection
00174   clusterProducerStrip_ = conf_.getParameter<edm::InputTag>("ClusterProducerStrip");
00175   clusterProducerPix_ = conf_.getParameter<edm::InputTag>("ClusterProducerPix");
00176   // SiStrip Quality Label
00177   qualityLabel_  = conf_.getParameter<std::string>("StripQualityLabel");
00178   // cluster quality conditions 
00179   edm::ParameterSet cluster_condition = conf_.getParameter<edm::ParameterSet>("ClusterConditions");
00180   applyClusterQuality_ = cluster_condition.getParameter<bool>("On");
00181   sToNLowerLimit_      = cluster_condition.getParameter<double>("minStoN");
00182   sToNUpperLimit_      = cluster_condition.getParameter<double>("maxStoN");
00183   widthLowerLimit_     = cluster_condition.getParameter<double>("minWidth"); 
00184   widthUpperLimit_     = cluster_condition.getParameter<double>("maxWidth"); 
00185 
00186   // Event History Producer
00187   historyProducer_ = conf_.getParameter<edm::InputTag>("HistoryProducer");
00188   // Apv Phase Producer
00189   apvPhaseProducer_ = conf_.getParameter<edm::InputTag>("ApvPhaseProducer");
00190 
00191   // Create DCS Status
00192   bool checkDCS    = conf_.getParameter<bool>("UseDCSFiltering");
00193   if (checkDCS) dcsStatus_ = new SiStripDCSStatus();
00194   else dcsStatus_ = 0; 
00195 
00196 } 
00197 
00198 SiStripMonitorCluster::~SiStripMonitorCluster() { 
00199   if (dcsStatus_)           delete dcsStatus_;
00200   if (genTriggerEventFlagBPTXfilter_    ) delete genTriggerEventFlagBPTXfilter_;
00201   if (genTriggerEventFlagPixelDCSfilter_) delete genTriggerEventFlagPixelDCSfilter_;
00202   if (genTriggerEventFlagStripDCSfilter_) delete genTriggerEventFlagStripDCSfilter_;
00203 }
00204 
00205 //--------------------------------------------------------------------------------------------
00206 void SiStripMonitorCluster::beginRun(const edm::Run& run, const edm::EventSetup& es){
00207 
00208   // Initialize the GenericTriggerEventFlag
00209   if ( genTriggerEventFlagBPTXfilter_->on() )
00210     genTriggerEventFlagBPTXfilter_->initRun( run, es );
00211   if ( genTriggerEventFlagPixelDCSfilter_->on() )
00212     genTriggerEventFlagPixelDCSfilter_->initRun( run, es );
00213   if ( genTriggerEventFlagStripDCSfilter_->on() )
00214     genTriggerEventFlagStripDCSfilter_->initRun( run, es );
00215 
00216   if (show_mechanical_structure_view) {
00217     unsigned long long cacheID = es.get<SiStripDetCablingRcd>().cacheIdentifier();
00218     if (m_cacheID_ != cacheID) {
00219       m_cacheID_ = cacheID;       
00220      edm::LogInfo("SiStripMonitorCluster") <<"SiStripMonitorCluster::beginRun: " 
00221                                             << " Creating MEs for new Cabling ";     
00222 
00223       createMEs(es);
00224     } 
00225   } else if (reset_each_run) {
00226     edm::LogInfo("SiStripMonitorCluster") <<"SiStripMonitorCluster::beginRun: " 
00227                                           << " Resetting MEs ";        
00228     for (std::map<uint32_t, ModMEs >::const_iterator idet = ModuleMEsMap.begin() ; idet!=ModuleMEsMap.end() ; idet++) {
00229       ResetModuleMEs(idet->first);
00230     }
00231   }
00232 }
00233 
00234 //--------------------------------------------------------------------------------------------
00235 void SiStripMonitorCluster::createMEs(const edm::EventSetup& es){
00236 
00237   if ( show_mechanical_structure_view ){
00238     // take from eventSetup the SiStripDetCabling object - here will use SiStripDetControl later on
00239     es.get<SiStripDetCablingRcd>().get(SiStripDetCabling_);
00240     
00241     // get list of active detectors from SiStripDetCabling 
00242     std::vector<uint32_t> activeDets;
00243     SiStripDetCabling_->addActiveDetectorsRawIds(activeDets);
00244     
00245     SiStripSubStructure substructure;
00246 
00247     SiStripFolderOrganizer folder_organizer;
00248     folder_organizer.setSiStripFolderName(topFolderName_);
00249     folder_organizer.setSiStripFolder();
00250 
00251 
00252     // Create TkHistoMap for Cluster
00253     if (clustertkhistomapon) {
00254       if (topFolderName_ == "SiStrip") tkmapcluster = new TkHistoMap("SiStrip/TkHistoMap","TkHMap_NumberOfCluster",0.,1);
00255       else tkmapcluster = new TkHistoMap(topFolderName_+"/TkHistoMap","TkHMap_NumberOfCluster",0.,0);
00256     }    
00257 
00258     // loop over detectors and book MEs
00259     edm::LogInfo("SiStripTkDQM|SiStripMonitorCluster")<<"nr. of activeDets:  "<<activeDets.size();
00260     for(std::vector<uint32_t>::iterator detid_iterator = activeDets.begin(); detid_iterator!=activeDets.end(); detid_iterator++){
00261       uint32_t detid = (*detid_iterator);
00262       // remove any eventual zero elements - there should be none, but just in case
00263       if(detid == 0) {
00264         activeDets.erase(detid_iterator);
00265         continue;
00266       }
00267       
00268       if (Mod_On_) {
00269         ModMEs mod_single;
00270         // set appropriate folder using SiStripFolderOrganizer
00271         folder_organizer.setDetectorFolder(detid); // pass the detid to this method
00272         if (reset_each_run) ResetModuleMEs(detid);
00273         createModuleMEs(mod_single, detid);
00274         // append to ModuleMEsMap
00275         ModuleMEsMap.insert( std::make_pair(detid, mod_single));
00276       }
00277       
00278       // Create Layer Level MEs if they are not created already
00279       std::pair<std::string,int32_t> det_layer_pair = folder_organizer.GetSubDetAndLayer(detid);
00280       SiStripHistoId hidmanager;
00281       std::string label = hidmanager.getSubdetid(detid,false);
00282       
00283       std::map<std::string, LayerMEs>::iterator iLayerME  = LayerMEsMap.find(label);
00284       if(iLayerME==LayerMEsMap.end()) {
00285         
00286         // get detids for the layer
00287         int32_t lnumber = det_layer_pair.second;
00288         std::vector<uint32_t> layerDetIds;        
00289         if (det_layer_pair.first == "TIB") {
00290           substructure.getTIBDetectors(activeDets,layerDetIds,lnumber,0,0,0);
00291         } else if (det_layer_pair.first == "TOB") {
00292           substructure.getTOBDetectors(activeDets,layerDetIds,lnumber,0,0);
00293         } else if (det_layer_pair.first == "TID" && lnumber > 0) {
00294           substructure.getTIDDetectors(activeDets,layerDetIds,2,abs(lnumber),0,0);
00295         } else if (det_layer_pair.first == "TID" && lnumber < 0) {
00296           substructure.getTIDDetectors(activeDets,layerDetIds,1,abs(lnumber),0,0);
00297         } else if (det_layer_pair.first == "TEC" && lnumber > 0) {
00298           substructure.getTECDetectors(activeDets,layerDetIds,2,abs(lnumber),0,0,0,0);
00299         } else if (det_layer_pair.first == "TEC" && lnumber < 0) {
00300           substructure.getTECDetectors(activeDets,layerDetIds,1,abs(lnumber),0,0,0,0);
00301         }
00302         LayerDetMap[label] = layerDetIds;
00303 
00304         // book Layer MEs 
00305         folder_organizer.setLayerFolder(detid,det_layer_pair.second);
00306         createLayerMEs(label, layerDetIds.size());
00307       }
00308       // book sub-detector plots
00309       std::pair<std::string,std::string> sdet_pair = folder_organizer.getSubDetFolderAndTag(detid);
00310       if (SubDetMEsMap.find(sdet_pair.second) == SubDetMEsMap.end()){
00311         dqmStore_->setCurrentFolder(sdet_pair.first);
00312         createSubDetMEs(sdet_pair.second);        
00313       }
00314     }//end of loop over detectors
00315 
00316     // Create Global Histogram
00317     if (globalswitchapvcycledbxth2on) {
00318       dqmStore_->setCurrentFolder(topFolderName_+"/MechanicalView/");
00319       edm::ParameterSet GlobalTH2Parameters =  conf_.getParameter<edm::ParameterSet>("TH2ApvCycleVsDBxGlobal");
00320       std::string HistoName = "DeltaBx_vs_ApvCycle";
00321       GlobalApvCycleDBxTH2 = dqmStore_->book2D(HistoName,HistoName,
00322                                                GlobalTH2Parameters.getParameter<int32_t>("Nbinsx"),
00323                                                GlobalTH2Parameters.getParameter<double>("xmin"),
00324                                                GlobalTH2Parameters.getParameter<double>("xmax"),
00325                                                GlobalTH2Parameters.getParameter<int32_t>("Nbinsy"),
00326                                                GlobalTH2Parameters.getParameter<double>("ymin"),
00327                                                GlobalTH2Parameters.getParameter<double>("ymax"));
00328       GlobalApvCycleDBxTH2->setAxisTitle("APV Cycle (Corrected Absolute Bx % 70)",1);
00329       GlobalApvCycleDBxTH2->setAxisTitle("Delta Bunch Crossing Cycle",2);
00330     }
00331 
00332     if (globalswitchcstripvscpix) {
00333       dqmStore_->setCurrentFolder(topFolderName_+"/MechanicalView/");
00334       edm::ParameterSet GlobalTH2Parameters =  conf_.getParameter<edm::ParameterSet>("TH2CStripVsCpixel");
00335       std::string HistoName = "StripClusVsPixClus";
00336       GlobalCStripVsCpix = dqmStore_->book2D(HistoName,HistoName,
00337                                                GlobalTH2Parameters.getParameter<int32_t>("Nbinsx"),
00338                                                GlobalTH2Parameters.getParameter<double>("xmin"),
00339                                                GlobalTH2Parameters.getParameter<double>("xmax"),
00340                                                GlobalTH2Parameters.getParameter<int32_t>("Nbinsy"),
00341                                                GlobalTH2Parameters.getParameter<double>("ymin"),
00342                                                GlobalTH2Parameters.getParameter<double>("ymax"));
00343       GlobalCStripVsCpix->setAxisTitle("Strip Clusters",1);
00344       GlobalCStripVsCpix->setAxisTitle("Pix Clusters",2);
00345     }
00346     
00347     if (globalswitchMultiRegions){
00348       dqmStore_->setCurrentFolder(topFolderName_+"/MechanicalView/");
00349       edm::ParameterSet GlobalTH2Parameters =  conf_.getParameter<edm::ParameterSet>("TH1MultiplicityRegions");
00350       std::string HistoName = "ClusterMultiplicityRegions";
00351       PixVsStripMultiplicityRegions = dqmStore_->book1D(HistoName,HistoName,
00352                                                GlobalTH2Parameters.getParameter<int32_t>("Nbinx"),
00353                                                GlobalTH2Parameters.getParameter<double>("xmin"),
00354                                                GlobalTH2Parameters.getParameter<double>("xmax"));
00355       PixVsStripMultiplicityRegions->setAxisTitle("");
00356       PixVsStripMultiplicityRegions->setBinLabel(1,"Main Diagonal");
00357       PixVsStripMultiplicityRegions->setBinLabel(2,"Strip Noise");
00358       PixVsStripMultiplicityRegions->setBinLabel(3,"High Strip Noise");
00359       PixVsStripMultiplicityRegions->setBinLabel(4,"Beam Background");
00360       PixVsStripMultiplicityRegions->setBinLabel(5,"No Strip Clusters");
00361     } 
00362 
00363     if (globalswitchmaindiagonalposition){
00364       dqmStore_->setCurrentFolder(topFolderName_+"/MechanicalView/");
00365       edm::ParameterSet GlobalTH1Parameters =  conf_.getParameter<edm::ParameterSet>("TH1MainDiagonalPosition");
00366       std::string HistoName = "MainDiagonal Position";
00367       GlobalMainDiagonalPosition = dqmStore_->book1D(HistoName,HistoName,
00368                                              GlobalTH1Parameters.getParameter<int32_t>("Nbinsx"),
00369                                              GlobalTH1Parameters.getParameter<double>("xmin"),
00370                                              GlobalTH1Parameters.getParameter<double>("xmax"));
00371       GlobalMainDiagonalPosition->setAxisTitle("atan(NPix/(k*NStrip))");
00372     }
00373 
00374     // TO BE ADDED !!!
00375     /*
00376     if ( globalswitchapvcycledbxth2on or globalswitchcstripvscpix or globalswitchMultiRegions or ClusterHisto_ ) {
00377       dqmStore_->setCurrentFolder(topFolderName_+"/MechanicalView/");
00378       std::string HistoName = "BPTX rate";
00379       BPTXrateTrend = dqmStore_->bookProfile(HistoName,HistoName, LSBin, LSMin, LSMax, 0, 10000.,"");
00380       BPTXrateTrend->getTH1()->SetBit(TH1::kCanRebin);
00381       BPTXrateTrend->setAxisTitle("#Lumi section",1);
00382       BPTXrateTrend->setAxisTitle("Number of BPTX events per LS",2);
00383     }
00384     */
00385 
00386     if (globalswitchstripnoise2apvcycle){
00387       dqmStore_->setCurrentFolder(topFolderName_+"/MechanicalView/");
00388       edm::ParameterSet GlobalTH1Parameters =  conf_.getParameter<edm::ParameterSet>("TH1StripNoise2ApvCycle");
00389       std::string HistoName = "StripNoise_ApvCycle";
00390       StripNoise2Cycle = dqmStore_->book1D(HistoName,HistoName,
00391                                              GlobalTH1Parameters.getParameter<int32_t>("Nbinsx"),
00392                                              GlobalTH1Parameters.getParameter<double>("xmin"),
00393                                              GlobalTH1Parameters.getParameter<double>("xmax"));
00394       StripNoise2Cycle->setAxisTitle("APV Cycle");
00395     }
00396 
00397     if (globalswitchstripnoise3apvcycle){
00398       dqmStore_->setCurrentFolder(topFolderName_+"/MechanicalView/");
00399       edm::ParameterSet GlobalTH1Parameters =  conf_.getParameter<edm::ParameterSet>("TH1StripNoise3ApvCycle");
00400       std::string HistoName = "HighStripNoise_ApvCycle";
00401       StripNoise3Cycle = dqmStore_->book1D(HistoName,HistoName,
00402                                              GlobalTH1Parameters.getParameter<int32_t>("Nbinsx"),
00403                                              GlobalTH1Parameters.getParameter<double>("xmin"),
00404                                              GlobalTH1Parameters.getParameter<double>("xmax"));
00405       StripNoise3Cycle->setAxisTitle("APV Cycle");
00406     }
00407 
00408     if (ClusterHisto_){
00409       dqmStore_->setCurrentFolder(topFolderName_+"/MechanicalView/");
00410       edm::ParameterSet PixelCluster =  conf_.getParameter<edm::ParameterSet>("TH1NClusPx");
00411       std::string HistoName = "NumberOfClustersInPixel";
00412       NumberOfPixelClus = dqmStore_->book1D(HistoName, HistoName, 
00413                                             PixelCluster.getParameter<int32_t>("Nbinsx"),
00414                                             PixelCluster.getParameter<double>("xmin"),
00415                                             PixelCluster.getParameter<double>("xmax"));
00416       NumberOfPixelClus->setAxisTitle("# of Clusters in Pixel", 1);
00417       NumberOfPixelClus->setAxisTitle("Number of Events", 2);
00418       //
00419       edm::ParameterSet StripCluster =  conf_.getParameter<edm::ParameterSet>("TH1NClusStrip");
00420       HistoName = "NumberOfClustersInStrip";
00421       NumberOfStripClus = dqmStore_->book1D(HistoName, HistoName, 
00422                                             StripCluster.getParameter<int32_t>("Nbinsx"),
00423                                             StripCluster.getParameter<double>("xmin"),
00424                                             StripCluster.getParameter<double>("xmax"));
00425       NumberOfStripClus->setAxisTitle("# of Clusters in Strip", 1);
00426       NumberOfStripClus->setAxisTitle("Number of Events", 2);
00427     }
00428 
00429  
00430   }//end of if
00431 }//end of method
00432 
00433 //--------------------------------------------------------------------------------------------
00434 void SiStripMonitorCluster::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
00435 {
00436 
00437   // Filter out events if Trigger Filtering is requested
00438   passBPTXfilter_     = ( iEvent.isRealData() and genTriggerEventFlagBPTXfilter_->on()     ) ? genTriggerEventFlagBPTXfilter_->accept( iEvent, iSetup)     : true;
00439   passPixelDCSfilter_ = ( iEvent.isRealData() and genTriggerEventFlagPixelDCSfilter_->on() ) ? genTriggerEventFlagPixelDCSfilter_->accept( iEvent, iSetup) : true;
00440   passStripDCSfilter_ = ( iEvent.isRealData() and genTriggerEventFlagStripDCSfilter_->on() ) ? genTriggerEventFlagStripDCSfilter_->accept( iEvent, iSetup) : true;
00441   //  std::cout << "passBPTXfilter_ ? " << passBPTXfilter_ << std::endl;
00442 
00443   // Filter out events if DCS Event if requested
00444   if (dcsStatus_ && !dcsStatus_->getStatus(iEvent,iSetup)) return;
00445 
00446   runNb   = iEvent.id().run();
00447   eventNb++;
00448   float iOrbitSec      = iEvent.orbitNumber()/11223.0;
00449 
00450   int NPixClusters=0, NStripClusters=0, MultiplicityRegion=0;
00451   bool isPixValid=false;
00452 
00453   edm::ESHandle<SiStripNoises> noiseHandle;
00454   iSetup.get<SiStripNoisesRcd>().get(noiseHandle);
00455 
00456   edm::ESHandle<SiStripGain> gainHandle;
00457   iSetup.get<SiStripGainRcd>().get(gainHandle);
00458 
00459   edm::ESHandle<SiStripQuality> qualityHandle;
00460   iSetup.get<SiStripQualityRcd>().get(qualityLabel_,qualityHandle);
00461 
00462   iSetup.get<SiStripDetCablingRcd>().get(SiStripDetCabling_);
00463 
00464   // get collection of DetSetVector of clusters from Event
00465   edm::Handle< edmNew::DetSetVector<SiStripCluster> > cluster_detsetvektor;
00466   iEvent.getByLabel(clusterProducerStrip_, cluster_detsetvektor);
00467 
00468   //get pixel clusters
00469   edm::Handle< edmNew::DetSetVector<SiPixelCluster> > cluster_detsetvektor_pix;
00470   iEvent.getByLabel(clusterProducerPix_, cluster_detsetvektor_pix);
00471 
00472   if (!cluster_detsetvektor.isValid()) return;
00473   
00474   const edmNew::DetSetVector<SiStripCluster> * StrC= cluster_detsetvektor.product();
00475   NStripClusters= StrC->data().size(); 
00476   
00477   if (cluster_detsetvektor_pix.isValid()){
00478     const edmNew::DetSetVector<SiPixelCluster> * PixC= cluster_detsetvektor_pix.product();
00479     NPixClusters= PixC->data().size();
00480     isPixValid=true;
00481     MultiplicityRegion=FindRegion(NStripClusters,NPixClusters);  
00482 
00483     if ( passBPTXfilter_ and passPixelDCSfilter_ and passStripDCSfilter_ ) {
00484       if (globalswitchcstripvscpix) GlobalCStripVsCpix->Fill(NStripClusters,NPixClusters);
00485       if (globalswitchmaindiagonalposition && NStripClusters > 0) GlobalMainDiagonalPosition->Fill(atan(NPixClusters/(k0*NStripClusters)));
00486       if (globalswitchMultiRegions) PixVsStripMultiplicityRegions->Fill(MultiplicityRegion);
00487     }
00488    
00489     if (ClusterHisto_){
00490       if ( passBPTXfilter_ and passPixelDCSfilter_ )
00491         NumberOfPixelClus->Fill(NPixClusters);
00492       if ( passBPTXfilter_ and passStripDCSfilter_ )
00493         NumberOfStripClus->Fill(NStripClusters);
00494     }
00495   }
00496   // initialise # of clusters to zero
00497   for (std::map<std::string, SubDetMEs>::iterator iSubdet  = SubDetMEsMap.begin();
00498        iSubdet != SubDetMEsMap.end(); iSubdet++) {
00499     iSubdet->second.totNClusters = 0;
00500   }
00501 
00502   SiStripFolderOrganizer folder_organizer;
00503   bool found_layer_me = false;
00504   for (std::map<std::string, std::vector< uint32_t > >::const_iterator iterLayer = LayerDetMap.begin();
00505        iterLayer != LayerDetMap.end(); iterLayer++) {
00506     
00507     std::string layer_label = iterLayer->first;
00508     
00509     int ncluster_layer = 0;
00510     std::map<std::string, LayerMEs>::iterator iLayerME = LayerMEsMap.find(layer_label);
00511     
00512     //get Layer MEs 
00513     LayerMEs layer_single;
00514     if(iLayerME != LayerMEsMap.end()) {
00515        layer_single = iLayerME->second; 
00516        found_layer_me = true;
00517      } 
00518 
00519     bool found_module_me = false;
00520     uint16_t iDet = 0;
00521     std::string subdet_label = ""; 
00522     // loop over all modules in the layer
00523     for (std::vector< uint32_t >::const_iterator iterDets = iterLayer->second.begin() ; 
00524          iterDets != iterLayer->second.end() ; iterDets++) {
00525       iDet++;
00526       // detid and type of ME
00527       uint32_t detid = (*iterDets);
00528 
00529       // Get SubDet label once
00530       if (subdet_label.size() == 0) subdet_label = folder_organizer.getSubDetFolderAndTag(detid).second;
00531 
00532       // DetId and corresponding set of MEs
00533       ModMEs mod_single;
00534       if (Mod_On_) {
00535         std::map<uint32_t, ModMEs >::iterator imodME = ModuleMEsMap.find(detid);
00536         if (imodME != ModuleMEsMap.end()) {
00537           mod_single = imodME->second;
00538           found_module_me = true;
00539         } 
00540       } else found_module_me = false;
00541 
00542       edmNew::DetSetVector<SiStripCluster>::const_iterator isearch = cluster_detsetvektor->find(detid); // search  clusters of detid
00543     
00544       if(isearch==cluster_detsetvektor->end()){
00545         if(found_module_me && moduleswitchncluson && (mod_single.NumberOfClusters)){
00546           (mod_single.NumberOfClusters)->Fill(0.); // no clusters for this detector module,fill histogram with 0
00547         }
00548         if(clustertkhistomapon) tkmapcluster->fill(detid,0.);
00549         if (found_layer_me && layerswitchnumclusterprofon) layer_single.LayerNumberOfClusterProfile->Fill(iDet, 0.0);
00550         continue; // no clusters for this detid => jump to next step of loop
00551       }
00552       
00553       //cluster_detset is a structure, cluster_detset.data is a std::vector<SiStripCluster>, cluster_detset.id is uint32_t
00554       edmNew::DetSet<SiStripCluster> cluster_detset = (*cluster_detsetvektor)[detid]; // the statement above makes sure there exists an element with 'detid'
00555       
00556       // Filling TkHistoMap with number of clusters for each module 
00557       if(clustertkhistomapon) {
00558         tkmapcluster->fill(detid,static_cast<float>(cluster_detset.size()));
00559       }
00560 
00561       if(moduleswitchncluson && found_module_me && (mod_single.NumberOfClusters != NULL)){ // nr. of clusters per module
00562         (mod_single.NumberOfClusters)->Fill(static_cast<float>(cluster_detset.size()));
00563       }
00564 
00565       if (found_layer_me && layerswitchnumclusterprofon) 
00566         layer_single.LayerNumberOfClusterProfile->Fill(iDet, static_cast<float>(cluster_detset.size()));
00567       ncluster_layer +=  cluster_detset.size();
00568       
00569       short total_clusterized_strips = 0;
00570       
00571       SiStripNoises::Range detNoiseRange = noiseHandle->getRange(detid);
00572       SiStripApvGain::Range detGainRange = gainHandle->getRange(detid); 
00573       SiStripQuality::Range qualityRange = qualityHandle->getRange(detid);
00574       
00575       for(edmNew::DetSet<SiStripCluster>::const_iterator clusterIter = cluster_detset.begin(); clusterIter!= cluster_detset.end(); clusterIter++){
00576 
00577         const std::vector<uint8_t>& ampls = clusterIter->amplitudes();
00578         // cluster position
00579         float cluster_position = clusterIter->barycenter();
00580         // start defined as nr. of first strip beloning to the cluster
00581         short cluster_start    = clusterIter->firstStrip();
00582         // width defined as nr. of strips that belong to cluster
00583         short cluster_width    = ampls.size(); 
00584         // add nr of strips of this cluster to total nr. of clusterized strips
00585         total_clusterized_strips = total_clusterized_strips + cluster_width; 
00586         
00587         // cluster signal and noise from the amplitudes
00588         float cluster_signal = 0.0;
00589         float cluster_noise  = 0.0;
00590         int nrnonzeroamplitudes = 0;
00591         float noise2 = 0.0;
00592         float noise  = 0.0;
00593         for(uint iamp=0; iamp<ampls.size(); iamp++){
00594           if(ampls[iamp]>0){ // nonzero amplitude
00595             cluster_signal += ampls[iamp];
00596             if(!qualityHandle->IsStripBad(qualityRange, clusterIter->firstStrip()+iamp)){
00597               noise = noiseHandle->getNoise(clusterIter->firstStrip()+iamp,detNoiseRange)/gainHandle->getStripGain(clusterIter->firstStrip()+iamp, detGainRange);
00598             }
00599             noise2 += noise*noise;
00600             nrnonzeroamplitudes++;
00601           }
00602         } // End loop over cluster amplitude
00603         
00604         if (nrnonzeroamplitudes > 0) cluster_noise = sqrt(noise2/nrnonzeroamplitudes);
00605         
00606         if( applyClusterQuality_ &&
00607             (cluster_signal/cluster_noise < sToNLowerLimit_ ||
00608              cluster_signal/cluster_noise > sToNUpperLimit_ ||
00609              cluster_width < widthLowerLimit_ ||
00610              cluster_width > widthUpperLimit_) ) continue;  
00611         
00612         ClusterProperties cluster_properties;
00613         cluster_properties.charge    = cluster_signal;
00614         cluster_properties.position  = cluster_position;
00615         cluster_properties.start     = cluster_start;
00616         cluster_properties.width     = cluster_width;
00617         cluster_properties.noise     = cluster_noise;
00618         
00619         // Fill Module Level MEs
00620         if (found_module_me) fillModuleMEs(mod_single, cluster_properties);
00621         
00622         // Fill Layer Level MEs
00623         if (found_layer_me) {
00624           fillLayerMEs(layer_single, cluster_properties, iOrbitSec);
00625           if (layerswitchclusterwidthprofon) 
00626             layer_single.LayerClusterWidthProfile->Fill(iDet, cluster_width);
00627         }
00628       } // end loop over clusters
00629       
00630       short total_nr_strips = SiStripDetCabling_->nApvPairs(detid) * 2 * 128; // get correct # of avp pairs
00631       float local_occupancy = static_cast<float>(total_clusterized_strips)/static_cast<float>(total_nr_strips);
00632       if (found_module_me) {
00633         if(moduleswitchnrclusterizedstrip && mod_single.NrOfClusterizedStrips ){ // nr of clusterized strips
00634           mod_single.NrOfClusterizedStrips->Fill(static_cast<float>(total_clusterized_strips));
00635         }
00636         
00637         if(moduleswitchlocaloccupancy && mod_single.ModuleLocalOccupancy ){ // Occupancy
00638           mod_single.ModuleLocalOccupancy->Fill(local_occupancy);
00639         }
00640       }
00641       if (layerswitchlocaloccupancy && found_layer_me && layer_single.LayerLocalOccupancy) {
00642         fillME(layer_single.LayerLocalOccupancy,local_occupancy);
00643         if (createTrendMEs) fillME(layer_single.LayerLocalOccupancyTrend,iOrbitSec,local_occupancy);
00644       }
00645     }
00646     std::map<std::string, SubDetMEs>::iterator iSubdet  = SubDetMEsMap.find(subdet_label);
00647     if(iSubdet != SubDetMEsMap.end()) iSubdet->second.totNClusters += ncluster_layer; 
00648   }
00649   
00650   //  EventHistory 
00651   edm::Handle<EventWithHistory> event_history;
00652   iEvent.getByLabel(historyProducer_,event_history);
00653   
00654   // Phase of APV
00655   edm::Handle<APVCyclePhaseCollection> apv_phase_collection;
00656   iEvent.getByLabel(apvPhaseProducer_,apv_phase_collection);
00657 
00658   if (event_history.isValid() 
00659         && !event_history.failedToGet()
00660         && apv_phase_collection.isValid() 
00661         && !apv_phase_collection.failedToGet()) {
00662 
00663 
00664     long long dbx        = event_history->deltaBX();
00665     long long tbx        = event_history->absoluteBX();    
00666 
00667     bool global_histo_filled = false;
00668     bool MultiplicityRegion_Vs_APVcycle_filled=false;
00669 
00670     for (std::map<std::string, SubDetMEs>::iterator it = SubDetMEsMap.begin();
00671        it != SubDetMEsMap.end(); it++) {
00672       std::string sdet = it->first;
00673       //std::string sdet = sdet_tag.substr(0,sdet_tag.find_first_of("_"));
00674       SubDetMEs sdetmes = it->second;
00675 
00676       int the_phase = APVCyclePhaseCollection::invalid;
00677       long long tbx_corr = tbx;
00678 
00679       if (SubDetPhasePartMap.find(sdet) != SubDetPhasePartMap.end()) the_phase = apv_phase_collection->getPhase(SubDetPhasePartMap[sdet]);
00680       if(the_phase==APVCyclePhaseCollection::nopartition ||
00681          the_phase==APVCyclePhaseCollection::multiphase ||
00682          the_phase==APVCyclePhaseCollection::invalid) {
00683         the_phase=30;
00684         //std::cout << " subdet " << it->first << " not valid" << " MR " << MultiplicityRegion <<std::endl;
00685       }
00686       tbx_corr  -= the_phase;
00687       long long dbxincycle = event_history->deltaBXinCycle(the_phase);
00688       if (globalswitchapvcycledbxth2on && !global_histo_filled) { 
00689         GlobalApvCycleDBxTH2->Fill(tbx_corr%70,dbx);
00690         global_histo_filled = true;
00691       }
00692 
00693       if (isPixValid && !MultiplicityRegion_Vs_APVcycle_filled){        
00694         if (globalswitchstripnoise2apvcycle && MultiplicityRegion==2) {StripNoise2Cycle->Fill(tbx_corr%70);}
00695         if (globalswitchstripnoise3apvcycle && MultiplicityRegion==3) {StripNoise3Cycle->Fill(tbx_corr%70);}
00696         MultiplicityRegion_Vs_APVcycle_filled=true;
00697       }
00698 
00699       if (subdetswitchtotclusth1on)
00700           sdetmes.SubDetTotClusterTH1->Fill(sdetmes.totNClusters);
00701       if (subdetswitchtotclusprofon)
00702           sdetmes.SubDetTotClusterProf->Fill(iOrbitSec,sdetmes.totNClusters);
00703       if (subdetswitchapvcycleprofon)
00704         sdetmes.SubDetClusterApvProf->Fill(tbx_corr%70,sdetmes.totNClusters);
00705       if (subdetswitchapvcycleth2on)
00706         sdetmes.SubDetClusterApvTH2->Fill(tbx_corr%70,sdetmes.totNClusters);
00707       if (subdetswitchdbxcycleprofon)
00708         sdetmes.SubDetClusterDBxCycleProf->Fill(dbxincycle,sdetmes.totNClusters);
00709       if (subdetswitchapvcycledbxprof2on)
00710         sdetmes.SubDetApvDBxProf2->Fill(tbx_corr%70,dbx,sdetmes.totNClusters);
00711     }
00712   }
00713 }
00714 //
00715 // -- EndJob
00716 //
00717 void SiStripMonitorCluster::endJob(void){
00718   bool outputMEsInRootFile = conf_.getParameter<bool>("OutputMEsInRootFile");
00719   std::string outputFileName = conf_.getParameter<std::string>("OutputFileName");
00720  
00721   // save histos in a file
00722   if(outputMEsInRootFile) dqmStore_->save(outputFileName);
00723 }
00724 //
00725 // -- Reset MEs
00726 //------------------------------------------------------------------------------
00727 void SiStripMonitorCluster::ResetModuleMEs(uint32_t idet){
00728   std::map<uint32_t, ModMEs >::iterator pos = ModuleMEsMap.find(idet);
00729   ModMEs mod_me = pos->second;
00730 
00731   if (moduleswitchncluson)            mod_me.NumberOfClusters->Reset();
00732   if (moduleswitchclusposon)          mod_me.ClusterPosition->Reset();
00733   if (moduleswitchclusdigiposon)      mod_me.ClusterDigiPosition->Reset();
00734   if (moduleswitchclusstonVsposon)    mod_me.ClusterSignalOverNoiseVsPos->Reset();
00735   if (moduleswitchcluswidthon)        mod_me.ClusterWidth->Reset();
00736   if (moduleswitchcluschargeon)       mod_me.ClusterCharge->Reset();
00737   if (moduleswitchclusnoiseon)        mod_me.ClusterNoise->Reset();
00738   if (moduleswitchclusstonon)         mod_me.ClusterSignalOverNoise->Reset();
00739   if (moduleswitchlocaloccupancy)     mod_me.ModuleLocalOccupancy->Reset();
00740   if (moduleswitchnrclusterizedstrip) mod_me.NrOfClusterizedStrips->Reset(); 
00741 }
00742 //
00743 // -- Create Module Level MEs
00744 //
00745 void SiStripMonitorCluster::createModuleMEs(ModMEs& mod_single, uint32_t detid) {
00746 
00747   // use SistripHistoId for producing histogram id (and title)
00748   SiStripHistoId hidmanager;
00749   std::string hid;
00750   
00751   //nr. of clusters per module
00752   if(moduleswitchncluson) {
00753     hid = hidmanager.createHistoId("NumberOfClusters","det",detid);
00754     mod_single.NumberOfClusters = bookME1D("TH1nClusters", hid.c_str());
00755     dqmStore_->tag(mod_single.NumberOfClusters, detid);
00756     mod_single.NumberOfClusters->setAxisTitle("number of clusters in one detector module");
00757     mod_single.NumberOfClusters->getTH1()->StatOverflows(kTRUE);  // over/underflows in Mean calculation
00758   }
00759   
00760   //ClusterPosition
00761   if(moduleswitchclusposon) {
00762     short total_nr_strips = SiStripDetCabling_->nApvPairs(detid) * 2 * 128; // get correct # of avp pairs    
00763     hid = hidmanager.createHistoId("ClusterPosition","det",detid);
00764     mod_single.ClusterPosition = dqmStore_->book1D(hid, hid, total_nr_strips, 0.5, total_nr_strips+0.5);
00765     dqmStore_->tag(mod_single.ClusterPosition, detid);
00766     mod_single.ClusterPosition->setAxisTitle("cluster position [strip number +0.5]");
00767   }
00768 
00769   //ClusterDigiPosition
00770   if(moduleswitchclusdigiposon) {
00771     short total_nr_strips = SiStripDetCabling_->nApvPairs(detid) * 2 * 128; // get correct # of avp pairs    
00772     hid = hidmanager.createHistoId("ClusterDigiPosition","det",detid);
00773     mod_single.ClusterDigiPosition = dqmStore_->book1D(hid, hid, total_nr_strips, 0.5, total_nr_strips+0.5);
00774     dqmStore_->tag(mod_single.ClusterDigiPosition, detid);
00775     mod_single.ClusterDigiPosition->setAxisTitle("digi in cluster position [strip number +0.5]");
00776   }
00777   
00778   //ClusterWidth
00779   if(moduleswitchcluswidthon) {
00780     hid = hidmanager.createHistoId("ClusterWidth","det",detid);
00781     mod_single.ClusterWidth = bookME1D("TH1ClusterWidth", hid.c_str());
00782     dqmStore_->tag(mod_single.ClusterWidth, detid);
00783     mod_single.ClusterWidth->setAxisTitle("cluster width [nr strips]");
00784   }
00785   
00786   //ClusterCharge
00787   if(moduleswitchcluschargeon) {
00788     hid = hidmanager.createHistoId("ClusterCharge","det",detid);
00789     mod_single.ClusterCharge = bookME1D("TH1ClusterCharge", hid.c_str());
00790     dqmStore_->tag(mod_single.ClusterCharge, detid);
00791     mod_single.ClusterCharge->setAxisTitle("cluster charge [ADC]");
00792   }
00793   
00794   //ClusterNoise
00795   if(moduleswitchclusnoiseon) {
00796     hid = hidmanager.createHistoId("ClusterNoise","det",detid);
00797     mod_single.ClusterNoise = bookME1D("TH1ClusterNoise", hid.c_str());
00798     dqmStore_->tag(mod_single.ClusterNoise, detid);
00799     mod_single.ClusterNoise->setAxisTitle("cluster noise");
00800   }
00801   
00802   //ClusterSignalOverNoise
00803   if(moduleswitchclusstonon) {
00804     hid = hidmanager.createHistoId("ClusterSignalOverNoise","det",detid);
00805     mod_single.ClusterSignalOverNoise = bookME1D("TH1ClusterStoN", hid.c_str());
00806     dqmStore_->tag(mod_single.ClusterSignalOverNoise, detid);
00807     mod_single.ClusterSignalOverNoise->setAxisTitle("ratio of signal to noise for each cluster");
00808   }
00809 
00810   //ClusterSignalOverNoiseVsPos
00811   if(moduleswitchclusstonVsposon) {
00812     hid = hidmanager.createHistoId("ClusterSignalOverNoiseVsPos","det",detid);
00813     Parameters =  conf_.getParameter<edm::ParameterSet>("TH1ClusterStoNVsPos");
00814     mod_single.ClusterSignalOverNoiseVsPos= dqmStore_->bookProfile(hid.c_str(),hid.c_str(),
00815                                                                    Parameters.getParameter<int32_t>("Nbinx"),
00816                                                                    Parameters.getParameter<double>("xmin"),
00817                                                                    Parameters.getParameter<double>("xmax"),
00818                                                                    Parameters.getParameter<int32_t>("Nbiny"),
00819                                                                    Parameters.getParameter<double>("ymin"),
00820                                                                    Parameters.getParameter<double>("ymax")
00821                                                                    );
00822     dqmStore_->tag(mod_single.ClusterSignalOverNoiseVsPos, detid);
00823     mod_single.ClusterSignalOverNoiseVsPos->setAxisTitle("pos");
00824   }
00825   
00826   //ModuleLocalOccupancy
00827   if (moduleswitchlocaloccupancy) {
00828     hid = hidmanager.createHistoId("ClusterLocalOccupancy","det",detid);
00829     mod_single.ModuleLocalOccupancy = bookME1D("TH1ModuleLocalOccupancy", hid.c_str());
00830     dqmStore_->tag(mod_single.ModuleLocalOccupancy, detid);
00831     mod_single.ModuleLocalOccupancy->setAxisTitle("module local occupancy [% of clusterized strips]");
00832   }
00833   
00834   //NrOfClusterizedStrips
00835   if (moduleswitchnrclusterizedstrip) {
00836     hid = hidmanager.createHistoId("NrOfClusterizedStrips","det",detid);
00837     mod_single.NrOfClusterizedStrips = bookME1D("TH1NrOfClusterizedStrips", hid.c_str());
00838     dqmStore_->tag(mod_single.NrOfClusterizedStrips, detid);
00839     mod_single.NrOfClusterizedStrips->setAxisTitle("number of clusterized strips");
00840   }
00841 }  
00842 //
00843 // -- Create Module Level MEs
00844 //  
00845 void SiStripMonitorCluster::createLayerMEs(std::string label, int ndets) {
00846 
00847   SiStripHistoId hidmanager;
00848   
00849   LayerMEs layerMEs; 
00850   layerMEs.LayerClusterStoN = 0;
00851   layerMEs.LayerClusterStoNTrend = 0;
00852   layerMEs.LayerClusterCharge = 0;
00853   layerMEs.LayerClusterChargeTrend = 0;
00854   layerMEs.LayerClusterNoise = 0;
00855   layerMEs.LayerClusterNoiseTrend = 0;
00856   layerMEs.LayerClusterWidth = 0;
00857   layerMEs.LayerClusterWidthTrend = 0;
00858   layerMEs.LayerLocalOccupancy = 0;
00859   layerMEs.LayerLocalOccupancyTrend = 0;
00860   layerMEs.LayerNumberOfClusterProfile = 0;
00861   layerMEs.LayerClusterWidthProfile = 0;
00862   
00863   //Cluster Width
00864   if(layerswitchcluswidthon) {
00865     layerMEs.LayerClusterWidth=bookME1D("TH1ClusterWidth", hidmanager.createHistoLayer("Summary_ClusterWidth","layer",label,"").c_str()); 
00866     if (createTrendMEs) layerMEs.LayerClusterWidthTrend=bookMETrend("TH1ClusterWidth", hidmanager.createHistoLayer("Trend_ClusterWidth","layer",label,"").c_str()); 
00867   }
00868   
00869   //Cluster Noise
00870   if(layerswitchclusnoiseon) {
00871     layerMEs.LayerClusterNoise=bookME1D("TH1ClusterNoise", hidmanager.createHistoLayer("Summary_ClusterNoise","layer",label,"").c_str()); 
00872     if (createTrendMEs) layerMEs.LayerClusterNoiseTrend=bookMETrend("TH1ClusterNoise", hidmanager.createHistoLayer("Trend_ClusterNoise","layer",label,"").c_str()); 
00873   }
00874   
00875   //Cluster Charge
00876   if(layerswitchcluschargeon) {
00877     layerMEs.LayerClusterCharge=bookME1D("TH1ClusterCharge", hidmanager.createHistoLayer("Summary_ClusterCharge","layer",label,"").c_str());
00878     if (createTrendMEs) layerMEs.LayerClusterChargeTrend=bookMETrend("TH1ClusterCharge", hidmanager.createHistoLayer("Trend_ClusterCharge","layer",label,"").c_str());
00879   }
00880   
00881   //Cluster StoN
00882   if(layerswitchclusstonon) {
00883     layerMEs.LayerClusterStoN=bookME1D("TH1ClusterStoN", hidmanager.createHistoLayer("Summary_ClusterSignalOverNoise","layer",label,"").c_str());
00884     if (createTrendMEs) layerMEs.LayerClusterStoNTrend=bookMETrend("TH1ClusterStoN", hidmanager.createHistoLayer("Trend_ClusterSignalOverNoise","layer",label,"").c_str());
00885   }
00886   
00887   //Cluster Occupancy
00888   if(layerswitchlocaloccupancy) {
00889     layerMEs.LayerLocalOccupancy=bookME1D("TH1ModuleLocalOccupancy", hidmanager.createHistoLayer("Summary_ClusterLocalOccupancy","layer",label,"").c_str());  
00890     if (createTrendMEs) layerMEs.LayerLocalOccupancyTrend=bookMETrend("TH1ModuleLocalOccupancy", hidmanager.createHistoLayer("Trend_ClusterLocalOccupancy","layer",label,"").c_str());  
00891     
00892   }
00893   
00894   // # of Cluster Profile 
00895   if(layerswitchnumclusterprofon) {
00896     std::string hid = hidmanager.createHistoLayer("NumberOfClusterProfile","layer",label,"");
00897     layerMEs.LayerNumberOfClusterProfile = dqmStore_->bookProfile(hid, hid, ndets, 0.5, ndets+0.5,21, -0.5, 20.5);      
00898   }
00899   
00900   // Cluster Width Profile 
00901   if(layerswitchclusterwidthprofon) {
00902     std::string hid = hidmanager.createHistoLayer("ClusterWidthProfile","layer",label,"");      
00903     layerMEs.LayerClusterWidthProfile = dqmStore_->bookProfile(hid, hid, ndets, 0.5, ndets+0.5, 20, -0.5, 19.5);      
00904   }
00905 
00906   LayerMEsMap[label]=layerMEs;
00907 }
00908 //
00909 // -- Create SubDetector MEs
00910 //
00911 void SiStripMonitorCluster::createSubDetMEs(std::string label) {
00912 
00913   SubDetMEs subdetMEs;
00914   subdetMEs.totNClusters              = 0;
00915   subdetMEs.SubDetTotClusterTH1       = 0;
00916   subdetMEs.SubDetTotClusterProf      = 0;
00917   subdetMEs.SubDetClusterApvProf      = 0;
00918   subdetMEs.SubDetClusterApvTH2       = 0;
00919   subdetMEs.SubDetClusterDBxCycleProf = 0;
00920   subdetMEs.SubDetApvDBxProf2         = 0;
00921   
00922   std::string HistoName;
00923   // Total Number of Cluster - 1D 
00924   if (subdetswitchtotclusth1on){
00925     HistoName = "TotalNumberOfCluster__" + label;
00926     subdetMEs.SubDetTotClusterTH1 = bookME1D("TH1TotalNumberOfClusters",HistoName.c_str());
00927     subdetMEs.SubDetTotClusterTH1->setAxisTitle("Total number of clusters in subdetector");
00928     subdetMEs.SubDetTotClusterTH1->getTH1()->StatOverflows(kTRUE);  // over/underflows in Mean calculation
00929   }
00930   // Total Number of Cluster vs Time - Profile
00931   if (subdetswitchtotclusprofon){
00932     edm::ParameterSet Parameters =  conf_.getParameter<edm::ParameterSet>("TProfTotalNumberOfClusters");
00933     HistoName = "TotalNumberOfClusterProfile__" + label;
00934     subdetMEs.SubDetTotClusterProf = dqmStore_->bookProfile(HistoName,HistoName,
00935                                                             Parameters.getParameter<int32_t>("Nbins"),
00936                                                             Parameters.getParameter<double>("xmin"),
00937                                                             Parameters.getParameter<double>("xmax"),
00938                                                             100, //that parameter should not be there !?
00939                                                             Parameters.getParameter<double>("ymin"),
00940                                                             Parameters.getParameter<double>("ymax"),
00941                                                             "" );
00942     subdetMEs.SubDetTotClusterProf->setAxisTitle("Event Time (Seconds)",1);
00943     if (subdetMEs.SubDetTotClusterProf->kind() == MonitorElement::DQM_KIND_TPROFILE) subdetMEs.SubDetTotClusterProf->getTH1()->SetBit(TH1::kCanRebin);
00944   }
00945   // Total Number of Cluster vs APV cycle - Profile
00946   if(subdetswitchapvcycleprofon){
00947     edm::ParameterSet Parameters =  conf_.getParameter<edm::ParameterSet>("TProfClustersApvCycle");
00948     HistoName = "Cluster_vs_ApvCycle__" + label;
00949     subdetMEs.SubDetClusterApvProf=dqmStore_->bookProfile(HistoName,HistoName,
00950                                                           Parameters.getParameter<int32_t>("Nbins"),
00951                                                           Parameters.getParameter<double>("xmin"),
00952                                                           Parameters.getParameter<double>("xmax"),
00953                                                           200, //that parameter should not be there !?
00954                                                           Parameters.getParameter<double>("ymin"),
00955                                                           Parameters.getParameter<double>("ymax"),
00956                                                           "" );
00957     subdetMEs.SubDetClusterApvProf->setAxisTitle("Apv Cycle (Corrected Absolute Bx % 70)",1);
00958   }
00959   
00960   // Total Number of Clusters vs ApvCycle - 2D 
00961   if(subdetswitchapvcycleth2on){
00962     edm::ParameterSet Parameters =  conf_.getParameter<edm::ParameterSet>("TH2ClustersApvCycle");
00963     HistoName = "Cluster_vs_ApvCycle_2D__" + label;
00964     // Adjusting the scale for 2D histogram
00965     double h2ymax = 9999.0;     
00966     double yfact = Parameters.getParameter<double>("yfactor");
00967     if(label.find("TIB") != std::string::npos) h2ymax = (6984.*256.)*yfact;
00968     else if (label.find("TID") != std::string::npos) h2ymax = (2208.*256.)*yfact;
00969     else if (label.find("TOB") != std::string::npos) h2ymax = (12906.*256.)*yfact;
00970     else if (label.find("TEC") != std::string::npos) h2ymax = (7552.*2.*256.)*yfact;
00971     
00972     subdetMEs.SubDetClusterApvTH2=dqmStore_->book2D(HistoName,HistoName,
00973                                                     Parameters.getParameter<int32_t>("Nbinsx"),
00974                                                     Parameters.getParameter<double>("xmin"),
00975                                                     Parameters.getParameter<double>("xmax"),
00976                                                     Parameters.getParameter<int32_t>("Nbinsy"),
00977                                                     Parameters.getParameter<double>("ymin"),
00978                                                     h2ymax);
00979     subdetMEs.SubDetClusterApvTH2->setAxisTitle("Apv Cycle (Corrected Absolute Bx % 70))",1);
00980     subdetMEs.SubDetClusterApvTH2->setAxisTitle("Total # of Clusters",2);
00981    
00982   }
00983   // Total Number of Cluster vs DeltaBxCycle - Profile
00984   if(subdetswitchdbxcycleprofon){
00985     edm::ParameterSet Parameters =  conf_.getParameter<edm::ParameterSet>("TProfClustersVsDBxCycle");
00986     HistoName = "Cluster_vs_DeltaBxCycle__" + label;
00987     subdetMEs.SubDetClusterDBxCycleProf = dqmStore_->bookProfile(HistoName,HistoName,
00988                                                                  Parameters.getParameter<int32_t>("Nbins"),
00989                                                                  Parameters.getParameter<double>("xmin"),
00990                                                                  Parameters.getParameter<double>("xmax"),
00991                                                                  200, //that parameter should not be there !?
00992                                                                  Parameters.getParameter<double>("ymin"),
00993                                                                  Parameters.getParameter<double>("ymax"),
00994                                                                  "" );
00995     subdetMEs.SubDetClusterDBxCycleProf->setAxisTitle("Delta Bunch Crossing Cycle",1);
00996   }
00997   // DeltaBx vs ApvCycle - 2DProfile
00998   if(subdetswitchapvcycledbxprof2on){
00999     edm::ParameterSet Parameters =  conf_.getParameter<edm::ParameterSet>("TProf2ApvCycleVsDBx");
01000     HistoName = "DeltaBx_vs_ApvCycle__" + label;
01001     subdetMEs.SubDetApvDBxProf2 = dqmStore_->bookProfile2D(HistoName,HistoName,
01002                                                            Parameters.getParameter<int32_t>("Nbinsx"),
01003                                                            Parameters.getParameter<double>("xmin"),
01004                                                            Parameters.getParameter<double>("xmax"),
01005                                                            Parameters.getParameter<int32_t>("Nbinsy"),
01006                                                            Parameters.getParameter<double>("ymin"),
01007                                                            Parameters.getParameter<double>("ymax"),
01008                                                            Parameters.getParameter<double>("zmin"),
01009                                                            Parameters.getParameter<double>("zmax"),
01010                                                            "" );
01011     subdetMEs.SubDetApvDBxProf2->setAxisTitle("APV Cycle (Corrected Absolute Bx % 70)",1);
01012     subdetMEs.SubDetApvDBxProf2->setAxisTitle("Delta Bunch Crossing Cycle",2);
01013   }
01014   SubDetMEsMap[label]=subdetMEs;
01015 }
01016 
01017 //
01018 // -- Fill Module Level Histograms
01019 //
01020 void SiStripMonitorCluster::fillModuleMEs(ModMEs& mod_mes, ClusterProperties& cluster) {
01021   
01022   if(moduleswitchclusposon && (mod_mes.ClusterPosition)) // position of cluster
01023     (mod_mes.ClusterPosition)->Fill(cluster.position);
01024   
01025   // position of digis in cluster
01026   if(moduleswitchclusdigiposon && (mod_mes.ClusterDigiPosition)) {
01027     for(int ipos=cluster.start+1; ipos<=cluster.start+cluster.width; ipos++){
01028       (mod_mes.ClusterDigiPosition)->Fill(ipos);
01029     }
01030   }
01031 
01032   if(moduleswitchcluswidthon && (mod_mes.ClusterWidth)) // width of cluster
01033     (mod_mes.ClusterWidth)->Fill(static_cast<float>(cluster.width));
01034  
01035   if(moduleswitchclusstonon && (mod_mes.ClusterSignalOverNoise)) {// SignalToNoise
01036     if (cluster.noise > 0) 
01037       (mod_mes.ClusterSignalOverNoise)->Fill(cluster.charge/cluster.noise);
01038   }
01039 
01040   if(moduleswitchclusstonVsposon && (mod_mes.ClusterSignalOverNoiseVsPos)) {// SignalToNoise
01041     if (cluster.noise > 0) 
01042       (mod_mes.ClusterSignalOverNoiseVsPos)->Fill(cluster.position,cluster.charge/cluster.noise);
01043   }
01044 
01045   if(moduleswitchclusnoiseon && (mod_mes.ClusterNoise))  // Noise
01046     (mod_mes.ClusterNoise)->Fill(cluster.noise);
01047 
01048   if(moduleswitchcluschargeon && (mod_mes.ClusterCharge)) // charge of cluster
01049     (mod_mes.ClusterCharge)->Fill(cluster.charge);
01050   
01051 } 
01052 //
01053 // -- Fill Layer Level MEs
01054 //
01055 void SiStripMonitorCluster::fillLayerMEs(LayerMEs& layerMEs, ClusterProperties& cluster, float timeinorbit) { 
01056   if(layerswitchclusstonon) {
01057     fillME(layerMEs.LayerClusterStoN  ,cluster.charge/cluster.noise);
01058     if (createTrendMEs) fillME(layerMEs.LayerClusterStoNTrend,timeinorbit,cluster.charge/cluster.noise);
01059   }
01060   
01061   if(layerswitchcluschargeon) {
01062     fillME(layerMEs.LayerClusterCharge,cluster.charge);
01063     if (createTrendMEs) fillME(layerMEs.LayerClusterChargeTrend,timeinorbit,cluster.charge);
01064   }
01065   
01066   if(layerswitchclusnoiseon) {
01067     fillME(layerMEs.LayerClusterNoise ,cluster.noise);
01068     if (createTrendMEs) fillME(layerMEs.LayerClusterNoiseTrend,timeinorbit,cluster.noise);
01069   }
01070   
01071   if(layerswitchcluswidthon) {
01072     fillME(layerMEs.LayerClusterWidth ,cluster.width);
01073     if (createTrendMEs) fillME(layerMEs.LayerClusterWidthTrend,timeinorbit,cluster.width);
01074   }
01075   
01076 }
01077 //------------------------------------------------------------------------------------------
01078 MonitorElement* SiStripMonitorCluster::bookMETrend(const char* ParameterSetLabel, const char* HistoName)
01079 {
01080   Parameters =  conf_.getParameter<edm::ParameterSet>(ParameterSetLabel);
01081   edm::ParameterSet ParametersTrend =  conf_.getParameter<edm::ParameterSet>("Trending");
01082   MonitorElement* me = dqmStore_->bookProfile(HistoName,HistoName,
01083                                               ParametersTrend.getParameter<int32_t>("Nbins"),
01084                                               //                                              0,
01085                                               ParametersTrend.getParameter<double>("xmin"),
01086                                               ParametersTrend.getParameter<double>("xmax"),
01087                                               //                                              ParametersTrend.getParameter<int32_t>("Nbins"),
01088                                               100, //that parameter should not be there !?
01089                                               ParametersTrend.getParameter<double>("ymin"),
01090                                               ParametersTrend.getParameter<double>("ymax"),
01091                                               "" );
01092   if(!me) return me;
01093   me->setAxisTitle("Event Time in Seconds",1);
01094   if (me->kind() == MonitorElement::DQM_KIND_TPROFILE) me->getTH1()->SetBit(TH1::kCanRebin);
01095   return me;
01096 }
01097 
01098 //------------------------------------------------------------------------------------------
01099 MonitorElement* SiStripMonitorCluster::bookME1D(const char* ParameterSetLabel, const char* HistoName)
01100 {
01101   Parameters =  conf_.getParameter<edm::ParameterSet>(ParameterSetLabel);
01102   return dqmStore_->book1D(HistoName,HistoName,
01103                            Parameters.getParameter<int32_t>("Nbinx"),
01104                            Parameters.getParameter<double>("xmin"),
01105                            Parameters.getParameter<double>("xmax")
01106                            );
01107 }
01108 
01109 int SiStripMonitorCluster::FindRegion(int nstrip,int npix){
01110   
01111   double kplus= k0*(1+dk0/100);
01112   double kminus=k0*(1-dk0/100);
01113   int region=0;
01114   
01115   if (nstrip!=0 && npix >= (nstrip*kminus-q0) && npix <=(nstrip*kplus+q0)) region=1; 
01116   else if (nstrip!=0 && npix < (nstrip*kminus-q0) &&  nstrip <= maxClus) region=2;
01117   else if (nstrip!=0 && npix < (nstrip*kminus-q0) &&  nstrip > maxClus) region=3;
01118   else if (nstrip!=0 && npix > (nstrip*kplus+q0)) region=4;
01119   else if (npix > minPix && nstrip==0) region=5;
01120   return region;
01121 
01122 }
01123 
01124 
01125
/data/refman/pasoursint/CMSSW_6_1_2_SLHC2/src/DQM/SiStripMonitorCluster/src/SiStripMonitorCluster.cc
