CMS 3D CMS Logo

DTResolutionAnalysisTest Class Reference

* DQM Test Client More...

#include <DQM/DTMonitorClient/src/DTResolutionAnalysisTest.h>

Inheritance diagram for DTResolutionAnalysisTest:

edm::EDAnalyzer

List of all members.

Public Member Functions

void analyze (const edm::Event &e, const edm::EventSetup &c)
 Analyze.
void beginJob (const edm::EventSetup &c)
 BeginJob.
void beginLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &context)
void bookHistos (int wh, int sect)
void bookHistos (int wh)
 book the summary histograms
 DTResolutionAnalysisTest (const edm::ParameterSet &ps)
 Constructor.
void endLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &c)
 DQM Client Diagnostic.
std::string getMEName (const DTSuperLayerId &slID)
 Get the ME name.
virtual ~DTResolutionAnalysisTest ()
 Destructor.

Private Member Functions

int slFromBin (int bin) const
int stationFromBin (int bin) const

Private Attributes

DQMStore * dbe
std::string folderRoot
std::map< std::pair< int, int >,
MonitorElement * > 		MeanHistos
edm::ESHandle< DTGeometry > muonGeom
int nevents
unsigned int nLumiSegs
int percentual
double permittedMeanRange
double permittedSigmaRange
int prescaleFactor
int run
std::map< std::pair< int, int >,
MonitorElement * > 		SigmaHistos
std::map< int, MonitorElement * > wheelMeanHistos
std::map< int, MonitorElement * > wheelSigmaHistos

Detailed Description

* DQM Test Client

Date
2008/12/10 10:28:23
Revision
1.5
Author:
G. Mila - INFN Torino

Definition at line 31 of file DTResolutionAnalysisTest.h.

Constructor & Destructor Documentation

DTResolutionAnalysisTest::DTResolutionAnalysisTest ( const edm::ParameterSet &  ps  ) 

Constructor.

Definition at line 38 of file DTResolutionAnalysisTest.cc.

References dbe, folderRoot, edm::ParameterSet::getUntrackedParameter(), permittedMeanRange, permittedSigmaRange, and prescaleFactor.

00038                                                                            {
00039 
00040   edm::LogVerbatim ("DTDQM|DTMonitorClient|DTResolutionAnalysisTest") << "[DTResolutionAnalysisTest]: Constructor";
00041 
00042   dbe = edm::Service<DQMStore>().operator->();
00043 
00044   prescaleFactor = ps.getUntrackedParameter<int>("diagnosticPrescale", 1);
00045   folderRoot = ps.getUntrackedParameter<string>("folderRoot", "Collector/FU0/");
00046 
00047   // permitted test range
00048   permittedMeanRange = ps.getUntrackedParameter<double>("permittedMeanRange",0.005); 
00049   permittedSigmaRange = ps.getUntrackedParameter<double>("permittedSigmaRange",0.01); 
00050   
00051  }

DTResolutionAnalysisTest::~DTResolutionAnalysisTest (  )  [virtual]

Destructor.

Definition at line 54 of file DTResolutionAnalysisTest.cc.

References nevents.

00054                                                    {
00055 
00056   edm::LogVerbatim ("DTDQM|DTMonitorClient|DTResolutionAnalysisTest") << "DTResolutionAnalysisTest: analyzed " << nevents << " events";
00057 
00058 }

Member Function Documentation

void DTResolutionAnalysisTest::analyze ( const edm::Event &  e,
const edm::EventSetup &  c 
) [virtual]

Analyze.

Implements edm::EDAnalyzer.

Definition at line 103 of file DTResolutionAnalysisTest.cc.

References nevents.

00103                                                                                      {
00104 
00105   nevents++;
00106 
00107 }

void DTResolutionAnalysisTest::beginJob ( const edm::EventSetup &  c  )  [virtual]

BeginJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 61 of file DTResolutionAnalysisTest.cc.

References bookHistos(), edm::EventSetup::get(), muonGeom, nevents, and muonGeometry::wheel.

00061                                                                    {
00062 
00063   edm::LogVerbatim ("DTDQM|DTMonitorClient|DTResolutionAnalysisTest") <<"[DTResolutionAnalysisTest]: BeginJob"; 
00064 
00065   nevents = 0;
00066 
00067   // Get the geometry
00068   context.get<MuonGeometryRecord>().get(muonGeom);
00069 
00070   // loop over all the CMS wheels, sectors & book the summary histos
00071   for (int wheel=-2; wheel<=2; wheel++){
00072     bookHistos(wheel);
00073     for (int sector=1; sector<=12; sector++){
00074        bookHistos(wheel, sector);
00075     }
00076   }
00077 
00078 }

void DTResolutionAnalysisTest::beginLuminosityBlock ( edm::LuminosityBlock const &  lumiSeg,
edm::EventSetup const &  context 
) [virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 81 of file DTResolutionAnalysisTest.cc.

References histo, run, edm::LuminosityBlock::run(), wheelMeanHistos, and wheelSigmaHistos.

00081                                                                                                              {
00082 
00083   edm::LogVerbatim ("DTDQM|DTMonitorClient|DTResolutionAnalysisTest") <<"[DTResolutionAnalysisTest]: Begin of LS transition";
00084 
00085   // Get the run number
00086   run = lumiSeg.run();
00087 
00088   // Reset the summary histo
00089   for(map<int, MonitorElement*> ::const_iterator histo = wheelMeanHistos.begin();
00090       histo != wheelMeanHistos.end();
00091       histo++) {
00092     (*histo).second->Reset();
00093   }
00094   for(map<int, MonitorElement*> ::const_iterator histo = wheelSigmaHistos.begin();
00095       histo != wheelSigmaHistos.end();
00096       histo++) {
00097     (*histo).second->Reset();
00098   }
00099 
00100 }

void DTResolutionAnalysisTest::bookHistos ( int  wh,
int  sect 
)

Definition at line 362 of file DTResolutionAnalysisTest.cc.

References DQMStore::book1D(), dbe, MeanHistos, DQMStore::setCurrentFolder(), SigmaHistos, and muonGeometry::wheel.

00362                                                           {
00363 
00364   stringstream wheel; wheel << wh;              
00365   stringstream sector; sector << sect;  
00366 
00367 
00368   string MeanHistoName =  "MeanTest_W" + wheel.str() + "_Sec" + sector.str(); 
00369   string SigmaHistoName =  "SigmaTest_W" + wheel.str() + "_Sec" + sector.str(); 
00370  
00371   string folder = "DT/02-Segments/Wheel" + wheel.str() + "/Sector" + sector.str();
00372   dbe->setCurrentFolder(folder);
00373 
00374   if(sect!=4 && sect!=10)
00375     MeanHistos[make_pair(wh,sect)] = dbe->book1D(MeanHistoName.c_str(),"Mean (from gaussian fit) of the residuals distribution",11,1,12);
00376   else
00377     MeanHistos[make_pair(wh,sect)] = dbe->book1D(MeanHistoName.c_str(),"Mean (from gaussian fit) of the residuals distribution",13,1,14);
00378   (MeanHistos[make_pair(wh,sect)])->setBinLabel(1,"MB1_SL1",1);
00379   (MeanHistos[make_pair(wh,sect)])->setBinLabel(2,"MB1_SL2",1);
00380   (MeanHistos[make_pair(wh,sect)])->setBinLabel(3,"MB1_SL3",1);
00381   (MeanHistos[make_pair(wh,sect)])->setBinLabel(4,"MB2_SL1",1);
00382   (MeanHistos[make_pair(wh,sect)])->setBinLabel(5,"MB2_SL2",1);
00383   (MeanHistos[make_pair(wh,sect)])->setBinLabel(6,"MB2_SL3",1);
00384   (MeanHistos[make_pair(wh,sect)])->setBinLabel(7,"MB3_SL1",1);
00385   (MeanHistos[make_pair(wh,sect)])->setBinLabel(8,"MB3_SL2",1);
00386   (MeanHistos[make_pair(wh,sect)])->setBinLabel(9,"MB3_SL3",1);
00387   (MeanHistos[make_pair(wh,sect)])->setBinLabel(10,"MB4_SL1",1);
00388   (MeanHistos[make_pair(wh,sect)])->setBinLabel(11,"MB4_SL3",1);
00389   if(sect==4){
00390     (MeanHistos[make_pair(wh,sect)])->setBinLabel(12,"MB4S4_SL1",1);
00391     (MeanHistos[make_pair(wh,sect)])->setBinLabel(13,"MB4S4_SL3",1);
00392   }
00393   if(sect==10){
00394     (MeanHistos[make_pair(wh,sect)])->setBinLabel(12,"MB4S10_SL1",1);
00395     (MeanHistos[make_pair(wh,sect)])->setBinLabel(13,"MB4S10_SL3",1);
00396   }
00397 
00398   if(sect!=4 && sect!=10)
00399     SigmaHistos[make_pair(wh,sect)] = dbe->book1D(SigmaHistoName.c_str(),"Sigma (from gaussian fit) of the residuals distribution",11,1,12);
00400   else
00401     SigmaHistos[make_pair(wh,sect)] = dbe->book1D(SigmaHistoName.c_str(),"Sigma (from gaussian fit) of the residuals distribution",13,1,14);
00402   (SigmaHistos[make_pair(wh,sect)])->setBinLabel(1,"MB1_SL1",1);  
00403   (SigmaHistos[make_pair(wh,sect)])->setBinLabel(2,"MB1_SL2",1);
00404   (SigmaHistos[make_pair(wh,sect)])->setBinLabel(3,"MB1_SL3",1);
00405   (SigmaHistos[make_pair(wh,sect)])->setBinLabel(4,"MB2_SL1",1);
00406   (SigmaHistos[make_pair(wh,sect)])->setBinLabel(5,"MB2_SL2",1);
00407   (SigmaHistos[make_pair(wh,sect)])->setBinLabel(6,"MB2_SL3",1);
00408   (SigmaHistos[make_pair(wh,sect)])->setBinLabel(7,"MB3_SL1",1);
00409   (SigmaHistos[make_pair(wh,sect)])->setBinLabel(8,"MB3_SL2",1);
00410   (SigmaHistos[make_pair(wh,sect)])->setBinLabel(9,"MB3_SL3",1);
00411   (SigmaHistos[make_pair(wh,sect)])->setBinLabel(10,"MB4_SL1",1);
00412   (SigmaHistos[make_pair(wh,sect)])->setBinLabel(11,"MB4_SL3",1);
00413   if(sect==4){
00414     (SigmaHistos[make_pair(wh,sect)])->setBinLabel(12,"MB4S13_SL1",1);
00415     (SigmaHistos[make_pair(wh,sect)])->setBinLabel(13,"MB4S13_SL3",1);
00416   }
00417   if(sect==10){
00418     (SigmaHistos[make_pair(wh,sect)])->setBinLabel(12,"MB4S14_SL1",1);
00419     (SigmaHistos[make_pair(wh,sect)])->setBinLabel(13,"MB4S14_SL3",1);
00420   }
00421 
00422 
00423 }

void DTResolutionAnalysisTest::bookHistos ( int  wh  ) 

book the summary histograms

Definition at line 253 of file DTResolutionAnalysisTest.cc.

References DQMStore::book2D(), dbe, permittedMeanRange, permittedSigmaRange, DQMStore::setCurrentFolder(), muonGeometry::wheel, wheelMeanHistos, and wheelSigmaHistos.

Referenced by beginJob().

00253                                                 { 
00254 
00255   stringstream wheel; wheel <<wh;
00256 
00257   dbe->setCurrentFolder("DT/02-Segments/00-MeanRes");
00258   string histoName =  "MeanSummaryRes_W" + wheel.str();
00259   stringstream meanRange; meanRange << (permittedMeanRange*10000);
00260   string histoTitle = "# of SL with |mean of res.| > " + meanRange.str() + "#mum (Wheel " + wheel.str() + ")";
00261   wheelMeanHistos[wh] = dbe->book2D(histoName.c_str(),histoTitle.c_str(),12,1,13,11,1,12);
00262   wheelMeanHistos[wh]->setAxisTitle("Sector",1);
00263   wheelMeanHistos[wh]->setBinLabel(1,"1",1);
00264   wheelMeanHistos[wh]->setBinLabel(2,"2",1);
00265   wheelMeanHistos[wh]->setBinLabel(3,"3",1);
00266   wheelMeanHistos[wh]->setBinLabel(4,"4",1);
00267   wheelMeanHistos[wh]->setBinLabel(5,"5",1);
00268   wheelMeanHistos[wh]->setBinLabel(6,"6",1);
00269   wheelMeanHistos[wh]->setBinLabel(7,"7",1);
00270   wheelMeanHistos[wh]->setBinLabel(8,"8",1);
00271   wheelMeanHistos[wh]->setBinLabel(9,"9",1);
00272   wheelMeanHistos[wh]->setBinLabel(10,"10",1);
00273   wheelMeanHistos[wh]->setBinLabel(11,"11",1);
00274   wheelMeanHistos[wh]->setBinLabel(12,"12",1);
00275   wheelMeanHistos[wh]->setBinLabel(1,"MB1_SL1",2);
00276   wheelMeanHistos[wh]->setBinLabel(2,"MB1_SL2",2);
00277   wheelMeanHistos[wh]->setBinLabel(3,"MB1_SL3",2);
00278   wheelMeanHistos[wh]->setBinLabel(4,"MB2_SL1",2);
00279   wheelMeanHistos[wh]->setBinLabel(5,"MB2_SL2",2);
00280   wheelMeanHistos[wh]->setBinLabel(6,"MB2_SL3",2);
00281   wheelMeanHistos[wh]->setBinLabel(7,"MB3_SL1",2);
00282   wheelMeanHistos[wh]->setBinLabel(8,"MB3_SL2",2);
00283   wheelMeanHistos[wh]->setBinLabel(9,"MB3_SL3",2);
00284   wheelMeanHistos[wh]->setBinLabel(10,"MB4_SL1",2);
00285   wheelMeanHistos[wh]->setBinLabel(11,"MB4_SL3",2); 
00286   
00287   if(wheelMeanHistos.find(3) == wheelMeanHistos.end()){
00288     string histoName =  "MeanSummaryRes";
00289     histoTitle = "# of SL with |mean of res.| > " + meanRange.str() + "#mum";
00290     wheelMeanHistos[3] = dbe->book2D(histoName.c_str(),histoTitle.c_str(),12,1,13,5,-2,3);
00291     wheelMeanHistos[3]->setAxisTitle("Sector",1);
00292     wheelMeanHistos[3]->setBinLabel(1,"1",1);
00293     wheelMeanHistos[3]->setBinLabel(2,"2",1);
00294     wheelMeanHistos[3]->setBinLabel(3,"3",1);
00295     wheelMeanHistos[3]->setBinLabel(4,"4",1);
00296     wheelMeanHistos[3]->setBinLabel(5,"5",1);
00297     wheelMeanHistos[3]->setBinLabel(6,"6",1);
00298     wheelMeanHistos[3]->setBinLabel(7,"7",1);
00299     wheelMeanHistos[3]->setBinLabel(8,"8",1);
00300     wheelMeanHistos[3]->setBinLabel(9,"9",1);
00301     wheelMeanHistos[3]->setBinLabel(10,"10",1);
00302     wheelMeanHistos[3]->setBinLabel(11,"11",1);
00303     wheelMeanHistos[3]->setBinLabel(12,"12",1);
00304     wheelMeanHistos[3]->setAxisTitle("Wheel",2);
00305   }
00306 
00307 
00308   dbe->setCurrentFolder("DT/02-Segments/01-SigmaRes");
00309   histoName =  "SigmaSummaryRes_W" + wheel.str();
00310   stringstream sigmaRange; sigmaRange << (permittedSigmaRange*10000);
00311   histoTitle = "# of SL with #sigma res. > " + sigmaRange.str() + "#mum (Wheel " + wheel.str() + ")";
00312   wheelSigmaHistos[wh] = dbe->book2D(histoName.c_str(),histoTitle.c_str(),12,1,13,11,1,12);
00313   wheelSigmaHistos[wh]->setAxisTitle("Sector",1);
00314   wheelSigmaHistos[wh]->setBinLabel(1,"1",1);
00315   wheelSigmaHistos[wh]->setBinLabel(2,"2",1);
00316   wheelSigmaHistos[wh]->setBinLabel(3,"3",1);
00317   wheelSigmaHistos[wh]->setBinLabel(4,"4",1);
00318   wheelSigmaHistos[wh]->setBinLabel(5,"5",1);
00319   wheelSigmaHistos[wh]->setBinLabel(6,"6",1);
00320   wheelSigmaHistos[wh]->setBinLabel(7,"7",1);
00321   wheelSigmaHistos[wh]->setBinLabel(8,"8",1);
00322   wheelSigmaHistos[wh]->setBinLabel(9,"9",1);
00323   wheelSigmaHistos[wh]->setBinLabel(10,"10",1);
00324   wheelSigmaHistos[wh]->setBinLabel(11,"11",1);
00325   wheelSigmaHistos[wh]->setBinLabel(12,"12",1);
00326   wheelSigmaHistos[wh]->setBinLabel(1,"MB1_SL1",2);
00327   wheelSigmaHistos[wh]->setBinLabel(2,"MB1_SL2",2);
00328   wheelSigmaHistos[wh]->setBinLabel(3,"MB1_SL3",2);
00329   wheelSigmaHistos[wh]->setBinLabel(4,"MB2_SL1",2);
00330   wheelSigmaHistos[wh]->setBinLabel(5,"MB2_SL2",2);
00331   wheelSigmaHistos[wh]->setBinLabel(6,"MB2_SL3",2);
00332   wheelSigmaHistos[wh]->setBinLabel(7,"MB3_SL1",2);
00333   wheelSigmaHistos[wh]->setBinLabel(8,"MB3_SL2",2);
00334   wheelSigmaHistos[wh]->setBinLabel(9,"MB3_SL3",2);
00335   wheelSigmaHistos[wh]->setBinLabel(10,"MB4_SL1",2);
00336   wheelSigmaHistos[wh]->setBinLabel(11,"MB4_SL3",2);
00337 
00338   if(wheelSigmaHistos.find(3) == wheelSigmaHistos.end()){
00339     string histoName =  "SigmaSummaryRes";
00340     histoTitle = "# of SL with #sigma res. > " + sigmaRange.str() + "#mum";
00341     wheelSigmaHistos[3] = dbe->book2D(histoName.c_str(),histoTitle.c_str(),14,1,15,5,-2,3);
00342     wheelSigmaHistos[3]->setAxisTitle("Sector",1);
00343     wheelSigmaHistos[3]->setBinLabel(1,"1",1);
00344     wheelSigmaHistos[3]->setBinLabel(2,"2",1);
00345     wheelSigmaHistos[3]->setBinLabel(3,"3",1);
00346     wheelSigmaHistos[3]->setBinLabel(4,"4",1);
00347     wheelSigmaHistos[3]->setBinLabel(5,"5",1);
00348     wheelSigmaHistos[3]->setBinLabel(6,"6",1);
00349     wheelSigmaHistos[3]->setBinLabel(7,"7",1);
00350     wheelSigmaHistos[3]->setBinLabel(8,"8",1);
00351     wheelSigmaHistos[3]->setBinLabel(9,"9",1);
00352     wheelSigmaHistos[3]->setBinLabel(10,"10",1);
00353     wheelSigmaHistos[3]->setBinLabel(11,"11",1);
00354     wheelSigmaHistos[3]->setBinLabel(12,"12",1);
00355     wheelSigmaHistos[3]->setBinLabel(13,"13",1);
00356     wheelSigmaHistos[3]->setBinLabel(14,"14",1);
00357     wheelSigmaHistos[3]->setAxisTitle("Wheel",2);
00358   }  
00359 }

void DTResolutionAnalysisTest::endLuminosityBlock ( edm::LuminosityBlock const &  lumiSeg,
edm::EventSetup const &  c 
) [virtual]

DQM Client Diagnostic.

Reimplemented from edm::EDAnalyzer.

Definition at line 111 of file DTResolutionAnalysisTest.cc.

References dbe, DQMStore::get(), MonitorElement::getMean(), getMEName(), MonitorElement::getRMS(), MonitorElement::getTH1F(), edm::LuminosityBlock::id(), edm::LuminosityBlockID::luminosityBlock(), mean(), MeanHistos, muonGeom, nLumiSegs, permittedMeanRange, permittedSigmaRange, prescaleFactor, DTChamberId::sector(), SigmaHistos, DTChamberId::station(), DTSuperLayerId::superLayer(), muonGeometry::wheel, DTChamberId::wheel(), wheelMeanHistos, and wheelSigmaHistos.

00111                                                                                                            {
00112   
00113 
00114   edm::LogVerbatim ("DTDQM|DTMonitorClient|DTResolutionAnalysisTest") <<"[DTResolutionAnalysisTest]: End of LS transition, performing the DQM client operation";
00115 
00116   // counts number of lumiSegs 
00117   nLumiSegs = lumiSeg.id().luminosityBlock();
00118   
00119   // prescale factor
00120   if ( nLumiSegs%prescaleFactor != 0 ) return;
00121 
00122    vector<DTChamber*>::const_iterator ch_it = muonGeom->chambers().begin();
00123   vector<DTChamber*>::const_iterator ch_end = muonGeom->chambers().end();
00124 
00125   for (; ch_it != ch_end; ++ch_it) {
00126 
00127     DTChamberId chID = (*ch_it)->id();
00128 
00129     // Fill the test histos
00130     int entry=-1;
00131     if(chID.station() == 1) entry=0;
00132     if(chID.station() == 2) entry=3;
00133     if(chID.station() == 3) entry=6;
00134     if(chID.station() == 4) entry=9;
00135 
00136     vector<const DTSuperLayer*>::const_iterator sl_it = (*ch_it)->superLayers().begin(); 
00137     vector<const DTSuperLayer*>::const_iterator sl_end = (*ch_it)->superLayers().end();
00138 
00139     for(; sl_it != sl_end; ++sl_it) {
00140 
00141       DTSuperLayerId slID = (*sl_it)->id();
00142       MonitorElement * res_histo = dbe->get(getMEName(slID));
00143 
00144       if(res_histo){
00145 
00146         // Gaussian Fit
00147         float statMean = res_histo->getMean(1);
00148         float statSigma = res_histo->getRMS(1);
00149         Double_t mean = -1;
00150         Double_t sigma = -1;
00151         TH1F * histo_root = res_histo->getTH1F();
00152         if(histo_root->GetEntries()>20){
00153           TF1 *gfit = new TF1("Gaussian","gaus",(statMean-(2*statSigma)),(statMean+(2*statSigma)));
00154           try {
00155             histo_root->Fit(gfit, "Q0");
00156           } catch (...) {
00157             edm::LogWarning ("DTDQM|DTMonitorModule|DTResolutionAnalysisTask")
00158               << "[DTResolutionAnalysisTask]: Exception when fitting SL : " << slID;
00159             continue;
00160           }
00161           if(gfit){
00162             mean = gfit->GetParameter(1); 
00163             sigma = gfit->GetParameter(2);
00164           }
00165           delete gfit;
00166         }
00167         else{
00168           edm::LogVerbatim ("DTDQM|DTMonitorModule|DTResolutionAnalysisTask")
00169             << "[DTResolutionAnalysisTask] Fit of " << slID
00170             << " not performed because # entries < 20 ";
00171         }
00172         
00173         int BinNumber = entry+slID.superLayer();
00174         if(BinNumber == 12) BinNumber=11;
00175         if((slID.sector()==13 || slID.sector()==14)  && slID.superLayer()==1) BinNumber=12;
00176         if((slID.sector()==13 || slID.sector()==14) && slID.superLayer()==3) BinNumber=13;
00177 
00178         // Fill the summary histos
00179         if(slID.sector()<13){
00180           MeanHistos[make_pair(slID.wheel(),slID.sector())]->setBinContent(BinNumber, mean);    
00181           SigmaHistos[make_pair(slID.wheel(),slID.sector())]->setBinContent(BinNumber, sigma);
00182         }
00183         if(slID.sector()==13){
00184           MeanHistos[make_pair(slID.wheel(),4)]->setBinContent(BinNumber, mean);        
00185           SigmaHistos[make_pair(slID.wheel(),4)]->setBinContent(BinNumber, sigma);
00186         }
00187         if(slID.sector()==14){
00188           MeanHistos[make_pair(slID.wheel(),10)]->setBinContent(BinNumber, mean);       
00189           SigmaHistos[make_pair(slID.wheel(),10)]->setBinContent(BinNumber, sigma);
00190         }
00191       }
00192 
00193     } // loop on SLs
00194   } // Loop on Stations
00195   
00196 
00197   for(int wheel=-2; wheel<=2; wheel++){
00198     for(int sector=1; sector<=12; sector++){
00199 
00200       int lastBin=-1;
00201       if(sector!=4 && sector!=10) lastBin=11;
00202       else lastBin=13;
00203 
00204       for (int bin=1; bin<=lastBin; bin++){
00205 
00206         // Mean test
00207         double mean = MeanHistos.find(make_pair(wheel,sector))->second->getBinContent(bin);
00208         if(mean<(-permittedMeanRange) || mean>permittedMeanRange){
00209 //        edm::LogError("DTDQM|DTMonitorClient|DTResolutionAnalysisTest") << "Bad mean channel: wh: " << wheel
00210 //                                                                        << " st: " << stationFromBin(bin)
00211 //                                                                        << " sect: " <<sector
00212 //                                                                        << " sl: " << slFromBin(bin)
00213 //                                                                        << " mean (cm): " << mean;
00214           // fill the wheel summary histos
00215           if(bin<12){
00216             wheelMeanHistos[wheel]->Fill(sector,bin);
00217             wheelMeanHistos[3]->Fill(sector,wheel);
00218           }
00219           else{
00220             wheelMeanHistos[wheel]->Fill(sector,11);
00221             wheelMeanHistos[3]->Fill(sector,wheel);
00222           }
00223         }
00224         
00225         // Sigma test
00226         double sigma = SigmaHistos.find(make_pair(wheel,sector))->second->getBinContent(bin);
00227         if(sigma>permittedSigmaRange){
00228 //        edm::LogError("DTDQM|DTMonitorClient|DTResolutionAnalysisTest") << "Bad sigma: wh: " << wheel
00229 //                                                                        << " st: " << stationFromBin(bin)
00230 //                                                                        << " sect: " <<sector
00231 //                                                                        << " sl: " << slFromBin(bin)
00232 //                                                                        << " sigma (cm): " << sigma;
00233           // fill the wheel summary histos
00234           if(bin<12){
00235             wheelSigmaHistos[wheel]->Fill(sector,bin);
00236             wheelSigmaHistos[3]->Fill(sector,wheel);
00237           }
00238           else{
00239             wheelSigmaHistos[wheel]->Fill(sector,11);
00240             wheelSigmaHistos[3]->Fill(sector,wheel);
00241           }
00242         }
00243 
00244       } // loop over bins
00245 
00246     } // loop over sectors
00247   } // loop over wheels
00248 
00249 }

string DTResolutionAnalysisTest::getMEName ( const DTSuperLayerId &  slID  ) 

Get the ME name.

Definition at line 426 of file DTResolutionAnalysisTest.cc.

References folderRoot, DTChamberId::sector(), DTChamberId::station(), DTSuperLayerId::superlayer(), muonGeometry::wheel, and DTChamberId::wheel().

Referenced by endLuminosityBlock().

00426                                                                       {
00427   
00428   stringstream wheel; wheel << slID.wheel();    
00429   stringstream station; station << slID.station();      
00430   stringstream sector; sector << slID.sector(); 
00431   stringstream superLayer; superLayer << slID.superlayer();
00432   
00433   string folderName = 
00434     folderRoot + "DT/02-Segments/Wheel" +  wheel.str() +
00435     "/Sector" + sector.str() +
00436     "/Station" + station.str() + "/";
00437   
00438   string histoname = folderName + "hResDist" 
00439     + "_W" + wheel.str() 
00440     + "_St" + station.str() 
00441     + "_Sec" + sector.str() 
00442     + "_SL" + superLayer.str(); 
00443   
00444   return histoname;
00445   
00446 }

int DTResolutionAnalysisTest::slFromBin ( int  bin  )  const [private]

Definition at line 455 of file DTResolutionAnalysisTest.cc.

00455                                                      {
00456   int ret = bin%3;
00457   if(ret == 0 || bin == 11) ret = 3;
00458   
00459   return ret;
00460 }

int DTResolutionAnalysisTest::stationFromBin ( int  bin  )  const [private]

Definition at line 450 of file DTResolutionAnalysisTest.cc.

00450                                                           {
00451   return (int) (bin /3.1)+1;
00452 }

Member Data Documentation

DQMStore* DTResolutionAnalysisTest::dbe [private]

Definition at line 70 of file DTResolutionAnalysisTest.h.

Referenced by bookHistos(), DTResolutionAnalysisTest(), and endLuminosityBlock().

std::string DTResolutionAnalysisTest::folderRoot [private]

Definition at line 68 of file DTResolutionAnalysisTest.h.

Referenced by DTResolutionAnalysisTest(), and getMEName().

std::map< std::pair<int,int> , MonitorElement* > DTResolutionAnalysisTest::MeanHistos [private]

Definition at line 79 of file DTResolutionAnalysisTest.h.

Referenced by bookHistos(), and endLuminosityBlock().

edm::ESHandle<DTGeometry> DTResolutionAnalysisTest::muonGeom [private]

Definition at line 76 of file DTResolutionAnalysisTest.h.

Referenced by beginJob(), and endLuminosityBlock().

int DTResolutionAnalysisTest::nevents [private]

Definition at line 63 of file DTResolutionAnalysisTest.h.

Referenced by analyze(), beginJob(), and ~DTResolutionAnalysisTest().

unsigned int DTResolutionAnalysisTest::nLumiSegs [private]

Definition at line 64 of file DTResolutionAnalysisTest.h.

Referenced by endLuminosityBlock().

int DTResolutionAnalysisTest::percentual [private]

Definition at line 67 of file DTResolutionAnalysisTest.h.

double DTResolutionAnalysisTest::permittedMeanRange [private]

Definition at line 73 of file DTResolutionAnalysisTest.h.

Referenced by bookHistos(), DTResolutionAnalysisTest(), and endLuminosityBlock().

double DTResolutionAnalysisTest::permittedSigmaRange [private]

Definition at line 74 of file DTResolutionAnalysisTest.h.

Referenced by bookHistos(), DTResolutionAnalysisTest(), and endLuminosityBlock().

int DTResolutionAnalysisTest::prescaleFactor [private]

Definition at line 65 of file DTResolutionAnalysisTest.h.

Referenced by DTResolutionAnalysisTest(), and endLuminosityBlock().

int DTResolutionAnalysisTest::run [private]

Definition at line 66 of file DTResolutionAnalysisTest.h.

Referenced by beginLuminosityBlock().

std::map< std::pair<int,int> , MonitorElement* > DTResolutionAnalysisTest::SigmaHistos [private]

Definition at line 80 of file DTResolutionAnalysisTest.h.

Referenced by bookHistos(), and endLuminosityBlock().

std::map< int, MonitorElement* > DTResolutionAnalysisTest::wheelMeanHistos [private]

Definition at line 82 of file DTResolutionAnalysisTest.h.

Referenced by beginLuminosityBlock(), bookHistos(), and endLuminosityBlock().

std::map< int, MonitorElement* > DTResolutionAnalysisTest::wheelSigmaHistos [private]

Definition at line 83 of file DTResolutionAnalysisTest.h.

Referenced by beginLuminosityBlock(), bookHistos(), and endLuminosityBlock().

The documentation for this class was generated from the following files:
Generated on Tue Jun 9 18:19:01 2009 for CMSSW by  doxygen 1.5.4