DTResolutionTest Class Reference

#include <DTResolutionTest.h>

Inheritance diagram for DTResolutionTest:

Public Member Functions
	DTResolutionTest (const edm::ParameterSet &ps)
	Constructor. More...
virtual	~DTResolutionTest ()
	Destructor. More...
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c)
	Analyze. More...
void	beginLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &context)
void	beginRun (const edm::Run &r, const edm::EventSetup &c)
	BeginRun. More...
void	bookHistos (const DTChamberId &ch)
	book the new ME More...
void	bookHistos (int wh)
	book the summary histograms More...
void	endJob ()
	Endjob. More...
void	endLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &c)
	DQM Client Diagnostic. More...
std::string	getMEName (const DTSuperLayerId &slID)
	Get the ME name. More...
std::string	getMEName2D (const DTSuperLayerId &slID)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Private Member Functions
int	slFromBin (int bin) const
int	stationFromBin (int bin) const

Private Attributes
std::map< std::pair< int, int > , int >	cmsMeanHistos
std::map< std::pair< int, int > , int >	cmsSigmaHistos
std::map< std::pair< int, int > , int >	cmsSlopeHistos
DQMStore *	dbe
std::map< std::pair< int, int > , bool >	MeanFilled
std::map< std::pair< int, int > , MonitorElement * >	MeanHistos
std::map< std::string, MonitorElement * >	MeanHistosSetRange
std::map< std::string, MonitorElement * >	MeanHistosSetRange2D
edm::ESHandle< DTGeometry >	muonGeom
int	nevents
unsigned int	nLumiSegs
edm::ParameterSet	parameters
int	percentual
int	prescaleFactor
int	run
std::map< std::pair< int, int > , bool >	SigmaFilled
std::map< std::pair< int, int > , MonitorElement * >	SigmaHistos
std::map< std::string, MonitorElement * >	SigmaHistosSetRange
std::map< std::string, MonitorElement * >	SigmaHistosSetRange2D
std::map< std::pair< int, int > , bool >	SlopeFilled
std::map< std::pair< int, int > , MonitorElement * >	SlopeHistos
std::map< std::string, MonitorElement * >	SlopeHistosSetRange
std::map< std::string, MonitorElement * >	SlopeHistosSetRange2D
std::map< int, MonitorElement * >	wheelMeanHistos
std::map< int, MonitorElement * >	wheelSigmaHistos
std::map< int, MonitorElement * >	wheelSlopeHistos

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)

Detailed Description

• DQM Test Client

  Date:

  2010/01/05 10:15:46

  Revision:

  1.13

  Author

  G. Mila - INFN Torino

Definition at line 41 of file DTResolutionTest.h.

Constructor & Destructor Documentation

DTResolutionTest::DTResolutionTest

(

const edm::ParameterSet &

ps

)

Constructor.

Definition at line 45 of file DTResolutionTest.cc.

References edm::ParameterSet::getUntrackedParameter(), cppFunctionSkipper::operator, Parameters::parameters, and dtDQMClient_cfg::prescaleFactor.

                                                           {

  edm::LogVerbatim ("resolution") << "[DTResolutionTest]: Constructor";
  parameters = ps;

  dbe = edm::Service<DQMStore>().operator->();
  if(ps.getUntrackedParameter<bool>("readFile", false))  
     dbe->open(ps.getUntrackedParameter<string>("inputFile", "residuals.root"));

  prescaleFactor = parameters.getUntrackedParameter<int>("diagnosticPrescale", 1);

  percentual = parameters.getUntrackedParameter<int>("BadSLpercentual", 10);

  //debug = parameters.getUntrackedParameter<bool>("debug", false);

}

DTResolutionTest::~DTResolutionTest ( )

virtual

Destructor.

Definition at line 63 of file DTResolutionTest.cc.

References nevents.

                                   {

  edm::LogVerbatim ("resolution") << "DTResolutionTest: analyzed " << nevents << " events";

}

Member Function Documentation

void DTResolutionTest::analyze ( const edm::Event &  e, const edm::EventSetup &  c  )

protectedvirtual

Analyze.

Implements edm::EDAnalyzer.

Definition at line 100 of file DTResolutionTest.cc.

References nevents.

                                                                             {

  nevents++;
  edm::LogVerbatim ("resolution") << "[DTResolutionTest]: "<<nevents<<" events";

}

void DTResolutionTest::beginLuminosityBlock ( edm::LuminosityBlock const &  lumiSeg, edm::EventSetup const &  context  )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 91 of file DTResolutionTest.cc.

References edm::LuminosityBlockBase::run(), and DTTTrigCorrFirst::run.

                                                                                                     {

  edm::LogVerbatim ("resolution") <<"[DTResolutionTest]: Begin of LS transition";

  // Get the run number
  run = lumiSeg.run();
}

void DTResolutionTest::beginRun ( const edm::Run &  r, const edm::EventSetup &  c  )

protectedvirtual

BeginRun.

Reimplemented from edm::EDAnalyzer.

Definition at line 70 of file DTResolutionTest.cc.

References bookHistos(), chambers, edm::EventSetup::get(), and nevents.

                                                                              {

  edm::LogVerbatim ("resolution") <<"[DTResolutionTest]: BeginRun";

  nevents = 0;
  // Get the geometry
  context.get<MuonGeometryRecord>().get(muonGeom);

  // book the histos
  for(int wheel=-2; wheel<3; wheel++){
    bookHistos(wheel);
  }
  vector<DTChamber*> chambers = muonGeom->chambers();
  for(vector<DTChamber*>::const_iterator chamber = chambers.begin();
      chamber != chambers.end(); ++chamber) {
    bookHistos((*chamber)->id());
  }

}

void DTResolutionTest::bookHistos ( const DTChamberId &  ch )

protected

book the new ME

Definition at line 504 of file DTResolutionTest.cc.

References Parameters::parameters, DTChamberId::sector(), and DTChamberId::wheel().

                                                        {

  stringstream wheel; wheel << ch.wheel();      
  stringstream sector; sector << ch.sector();   


  string MeanHistoName =  "MeanTest_" + parameters.getUntrackedParameter<string>("STEP", "STEP3") + "_W" + wheel.str() + "_Sec" + sector.str(); 
  string SigmaHistoName =  "SigmaTest_" + parameters.getUntrackedParameter<string>("STEP", "STEP3") + "_W" + wheel.str() + "_Sec" + sector.str(); 
  string SlopeHistoName =  "SlopeTest_" + parameters.getUntrackedParameter<string>("STEP", "STEP3") + "_W" + wheel.str() + "_Sec" + sector.str(); 

  dbe->setCurrentFolder("DT/Tests/DTResolution");

  // Book the histo for the mean value and set the axis labels
  MeanHistos[make_pair(ch.wheel(),ch.sector())] = dbe->book1D(MeanHistoName.c_str(),MeanHistoName.c_str(),11,0,11);
  (MeanHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(1,"MB1_SL1",1);
  (MeanHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(2,"MB1_SL2",1);
  (MeanHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(3,"MB1_SL3",1);
  (MeanHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(4,"MB2_SL1",1);
  (MeanHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(5,"MB2_SL2",1);
  (MeanHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(6,"MB2_SL3",1);
  (MeanHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(7,"MB3_SL1",1);
  (MeanHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(8,"MB3_SL2",1);
  (MeanHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(9,"MB3_SL3",1);
  (MeanHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(10,"MB4_SL1",1);
  (MeanHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(11,"MB4_SL3",1);


  // Book the histo for the sigma value and set the axis labels
  if(parameters.getUntrackedParameter<bool>("sigmaTest")){
    SigmaHistos[make_pair(ch.wheel(),ch.sector())] = dbe->book1D(SigmaHistoName.c_str(),SigmaHistoName.c_str(),11,0,11);
    (SigmaHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(1,"MB1_SL1",1);  
    (SigmaHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(2,"MB1_SL2",1);
    (SigmaHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(3,"MB1_SL3",1);
    (SigmaHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(4,"MB2_SL1",1);
    (SigmaHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(5,"MB2_SL2",1);
    (SigmaHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(6,"MB2_SL3",1);
    (SigmaHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(7,"MB3_SL1",1);
    (SigmaHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(8,"MB3_SL2",1);
    (SigmaHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(9,"MB3_SL3",1);
    (SigmaHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(10,"MB4_SL1",1);
    (SigmaHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(11,"MB4_SL3",1);
  }

  // Book the histo for the slope value and set the axis labels
  if(parameters.getUntrackedParameter<bool>("slopeTest")){
    SlopeHistos[make_pair(ch.wheel(),ch.sector())] = dbe->book1D(SlopeHistoName.c_str(),SlopeHistoName.c_str(),11,0,11);
    (SlopeHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(1,"MB1_SL1",1);  
    (SlopeHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(2,"MB1_SL2",1);
    (SlopeHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(3,"MB1_SL3",1);
    (SlopeHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(4,"MB2_SL1",1);
    (SlopeHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(5,"MB2_SL2",1);
    (SlopeHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(6,"MB2_SL3",1);
    (SlopeHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(7,"MB3_SL1",1);
    (SlopeHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(8,"MB3_SL2",1);
    (SlopeHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(9,"MB3_SL3",1);
    (SlopeHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(10,"MB4_SL1",1);
    (SlopeHistos[make_pair(ch.wheel(),ch.sector())])->setBinLabel(11,"MB4_SL3",1);
  }

  string HistoName = "W" + wheel.str() + "_Sec" + sector.str(); 

  if(parameters.getUntrackedParameter<bool>("meanWrongHisto")){
    string MeanHistoNameSetRange = "MeanWrong_" + parameters.getUntrackedParameter<string>("STEP", "STEP3") + "_W" + wheel.str() + "_Sec" + sector.str() + "_SetRange";
    MeanHistosSetRange[HistoName] = dbe->book1D(MeanHistoNameSetRange.c_str(),MeanHistoNameSetRange.c_str(),11,0.5,11.5);
    string MeanHistoNameSetRange2D = "MeanWrong_" + parameters.getUntrackedParameter<string>("STEP", "STEP3") + "_W" + wheel.str() + "_Sec" + sector.str() + "_SetRange" + "_2D";
    MeanHistosSetRange2D[HistoName] = dbe->book2D(MeanHistoNameSetRange2D.c_str(),MeanHistoNameSetRange2D.c_str(),11, 0.5, 11.5, 100, -0.05, 0.05);
  }
  
  if(parameters.getUntrackedParameter<bool>("sigmaTest")){
    string SigmaHistoNameSetRange =  "SigmaWrong_" + parameters.getUntrackedParameter<string>("STEP", "STEP3") + "_W" + wheel.str() + "_Sec" + sector.str()  + "_SetRange";
    SigmaHistosSetRange[HistoName] = dbe->book1D(SigmaHistoNameSetRange.c_str(),SigmaHistoNameSetRange.c_str(),11,0.5,11.5);
    string SigmaHistoNameSetRange2D =  "SigmaWrong_" + parameters.getUntrackedParameter<string>("STEP", "STEP3") + "_W" + wheel.str() + "_Sec" + sector.str()  + "_SetRange" + "_2D";
    SigmaHistosSetRange2D[HistoName] = dbe->book2D(SigmaHistoNameSetRange2D.c_str(),SigmaHistoNameSetRange2D.c_str(),11, 0.5, 11.5, 500, 0, 0.5);
  }

  if(parameters.getUntrackedParameter<bool>("slopeTest")){
    string SlopeHistoNameSetRange =  "SlopeWrong_" + parameters.getUntrackedParameter<string>("STEP", "STEP3") + "_W" + wheel.str() + "_Sec" + sector.str()  + "_SetRange";
    SlopeHistosSetRange[HistoName] = dbe->book1D(SlopeHistoNameSetRange.c_str(),SlopeHistoNameSetRange.c_str(),11,0.5,11.5);
    string SlopeHistoNameSetRange2D =  "SlopeWrong_" + parameters.getUntrackedParameter<string>("STEP", "STEP3") + "_W" + wheel.str() + "_Sec" + sector.str()  + "_SetRange" + "_2D";
    SlopeHistosSetRange2D[HistoName] = dbe->book2D(SlopeHistoNameSetRange2D.c_str(),SlopeHistoNameSetRange2D.c_str(),11, 0.5, 11.5, 200, -0.1, 0.1);
  }

}

void DTResolutionTest::bookHistos ( int  wh )

protected

book the summary histograms

Definition at line 589 of file DTResolutionTest.cc.

References Parameters::parameters.

                                        {
  
  dbe->setCurrentFolder("DT/CalibrationSummary");

  if(wheelMeanHistos.find(3) == wheelMeanHistos.end()){
    string histoName =  "MeanSummaryRes_testFailedByAtLeastBadSL_" + parameters.getUntrackedParameter<string>("STEP", "STEP3");
    wheelMeanHistos[3] = dbe->book2D(histoName.c_str(),histoName.c_str(),14,0,14,5,-2,3);
    wheelMeanHistos[3]->setBinLabel(1,"Sector1",1);
    wheelMeanHistos[3]->setBinLabel(1,"Sector1",1);
    wheelMeanHistos[3]->setBinLabel(2,"Sector2",1);
    wheelMeanHistos[3]->setBinLabel(3,"Sector3",1);
    wheelMeanHistos[3]->setBinLabel(4,"Sector4",1);
    wheelMeanHistos[3]->setBinLabel(5,"Sector5",1);
    wheelMeanHistos[3]->setBinLabel(6,"Sector6",1);
    wheelMeanHistos[3]->setBinLabel(7,"Sector7",1);
    wheelMeanHistos[3]->setBinLabel(8,"Sector8",1);
    wheelMeanHistos[3]->setBinLabel(9,"Sector9",1);
    wheelMeanHistos[3]->setBinLabel(10,"Sector10",1);
    wheelMeanHistos[3]->setBinLabel(11,"Sector11",1);
    wheelMeanHistos[3]->setBinLabel(12,"Sector12",1);
    wheelMeanHistos[3]->setBinLabel(13,"Sector13",1);
    wheelMeanHistos[3]->setBinLabel(14,"Sector14",1);
    wheelMeanHistos[3]->setBinLabel(1,"Wheel-2",2);
    wheelMeanHistos[3]->setBinLabel(2,"Wheel-1",2);
    wheelMeanHistos[3]->setBinLabel(3,"Wheel0",2);
    wheelMeanHistos[3]->setBinLabel(4,"Wheel+1",2);
    wheelMeanHistos[3]->setBinLabel(5,"Wheel+2",2);
  }

  if(parameters.getUntrackedParameter<bool>("sigmaTest")){
    if(wheelSigmaHistos.find(3) == wheelSigmaHistos.end()){
      string histoName =  "SigmaSummaryRes_testFailedByAtLeastBadSL_" + parameters.getUntrackedParameter<string>("STEP", "STEP3");
      wheelSigmaHistos[3] = dbe->book2D(histoName.c_str(),histoName.c_str(),14,0,14,5,-2,3);
      wheelSigmaHistos[3]->setBinLabel(1,"Sector1",1);
      wheelSigmaHistos[3]->setBinLabel(1,"Sector1",1);
      wheelSigmaHistos[3]->setBinLabel(2,"Sector2",1);
      wheelSigmaHistos[3]->setBinLabel(3,"Sector3",1);
      wheelSigmaHistos[3]->setBinLabel(4,"Sector4",1);
      wheelSigmaHistos[3]->setBinLabel(5,"Sector5",1);
      wheelSigmaHistos[3]->setBinLabel(6,"Sector6",1);
      wheelSigmaHistos[3]->setBinLabel(7,"Sector7",1);
      wheelSigmaHistos[3]->setBinLabel(8,"Sector8",1);
      wheelSigmaHistos[3]->setBinLabel(9,"Sector9",1);
      wheelSigmaHistos[3]->setBinLabel(10,"Sector10",1);
      wheelSigmaHistos[3]->setBinLabel(11,"Sector11",1);
      wheelSigmaHistos[3]->setBinLabel(12,"Sector12",1);
      wheelSigmaHistos[3]->setBinLabel(13,"Sector13",1);
      wheelSigmaHistos[3]->setBinLabel(14,"Sector14",1);
      wheelSigmaHistos[3]->setBinLabel(1,"Wheel-2",2);
      wheelSigmaHistos[3]->setBinLabel(2,"Wheel-1",2);
      wheelSigmaHistos[3]->setBinLabel(3,"Wheel0",2);
      wheelSigmaHistos[3]->setBinLabel(4,"Wheel+1",2);
      wheelSigmaHistos[3]->setBinLabel(5,"Wheel+2",2);
    }
  }

  if(parameters.getUntrackedParameter<bool>("slopeTest")){
    if(wheelSlopeHistos.find(3) == wheelSlopeHistos.end()){
      string histoName =  "SlopeSummaryRes_testFailedByAtLeastBadSL_" + parameters.getUntrackedParameter<string>("STEP", "STEP3");
      wheelSlopeHistos[3] = dbe->book2D(histoName.c_str(),histoName.c_str(),14,0,14,5,-2,3);
      wheelSlopeHistos[3]->setBinLabel(1,"Sector1",1);
      wheelSlopeHistos[3]->setBinLabel(1,"Sector1",1);
      wheelSlopeHistos[3]->setBinLabel(2,"Sector2",1);
      wheelSlopeHistos[3]->setBinLabel(3,"Sector3",1);
      wheelSlopeHistos[3]->setBinLabel(4,"Sector4",1);
      wheelSlopeHistos[3]->setBinLabel(5,"Sector5",1);
      wheelSlopeHistos[3]->setBinLabel(6,"Sector6",1);
      wheelSlopeHistos[3]->setBinLabel(7,"Sector7",1);
      wheelSlopeHistos[3]->setBinLabel(8,"Sector8",1);
      wheelSlopeHistos[3]->setBinLabel(9,"Sector9",1);
      wheelSlopeHistos[3]->setBinLabel(10,"Sector10",1);
      wheelSlopeHistos[3]->setBinLabel(11,"Sector11",1);
      wheelSlopeHistos[3]->setBinLabel(12,"Sector12",1);
      wheelSlopeHistos[3]->setBinLabel(13,"Sector13",1);
      wheelSlopeHistos[3]->setBinLabel(14,"Sector14",1);
      wheelSlopeHistos[3]->setBinLabel(1,"Wheel-2",2);
      wheelSlopeHistos[3]->setBinLabel(2,"Wheel-1",2);
      wheelSlopeHistos[3]->setBinLabel(3,"Wheel0",2);
      wheelSlopeHistos[3]->setBinLabel(4,"Wheel+1",2);
      wheelSlopeHistos[3]->setBinLabel(5,"Wheel+2",2);
    }
  }

  stringstream wheel; wheel <<wh;
  
  if(wheelMeanHistos.find(wh) == wheelMeanHistos.end()){
    string histoName =  "MeanSummaryRes_testFailed_" + parameters.getUntrackedParameter<string>("STEP", "STEP3") + "_W" + wheel.str();
    wheelMeanHistos[wh] = dbe->book2D(histoName.c_str(),histoName.c_str(),14,0,14,11,0,11);
    wheelMeanHistos[wh]->setBinLabel(1,"Sector1",1);
    wheelMeanHistos[wh]->setBinLabel(2,"Sector2",1);
    wheelMeanHistos[wh]->setBinLabel(3,"Sector3",1);
    wheelMeanHistos[wh]->setBinLabel(4,"Sector4",1);
    wheelMeanHistos[wh]->setBinLabel(5,"Sector5",1);
    wheelMeanHistos[wh]->setBinLabel(6,"Sector6",1);
    wheelMeanHistos[wh]->setBinLabel(7,"Sector7",1);
    wheelMeanHistos[wh]->setBinLabel(8,"Sector8",1);
    wheelMeanHistos[wh]->setBinLabel(9,"Sector9",1);
    wheelMeanHistos[wh]->setBinLabel(10,"Sector10",1);
    wheelMeanHistos[wh]->setBinLabel(11,"Sector11",1);
    wheelMeanHistos[wh]->setBinLabel(12,"Sector12",1);
    wheelMeanHistos[wh]->setBinLabel(13,"Sector13",1);
    wheelMeanHistos[wh]->setBinLabel(14,"Sector14",1);
    wheelMeanHistos[wh]->setBinLabel(1,"MB1_SL1",2);
    wheelMeanHistos[wh]->setBinLabel(2,"MB1_SL2",2);
    wheelMeanHistos[wh]->setBinLabel(3,"MB1_SL3",2);
    wheelMeanHistos[wh]->setBinLabel(4,"MB2_SL1",2);
    wheelMeanHistos[wh]->setBinLabel(5,"MB2_SL2",2);
    wheelMeanHistos[wh]->setBinLabel(6,"MB2_SL3",2);
    wheelMeanHistos[wh]->setBinLabel(7,"MB3_SL1",2);
    wheelMeanHistos[wh]->setBinLabel(8,"MB3_SL2",2);
    wheelMeanHistos[wh]->setBinLabel(9,"MB3_SL3",2);
    wheelMeanHistos[wh]->setBinLabel(10,"MB4_SL1",2);
    wheelMeanHistos[wh]->setBinLabel(11,"MB4_SL3",2);
  }

  if(parameters.getUntrackedParameter<bool>("sigmaTest")){
    if(wheelSigmaHistos.find(wh) == wheelSigmaHistos.end()){
      string histoName =  "SigmaSummaryRes_testFailed_" + parameters.getUntrackedParameter<string>("STEP", "STEP3") + "_W" + wheel.str();
      wheelSigmaHistos[wh] = dbe->book2D(histoName.c_str(),histoName.c_str(),14,0,14,11,0,11);
      wheelSigmaHistos[wh]->setBinLabel(1,"Sector1",1);
      wheelSigmaHistos[wh]->setBinLabel(2,"Sector2",1);
      wheelSigmaHistos[wh]->setBinLabel(3,"Sector3",1);
      wheelSigmaHistos[wh]->setBinLabel(4,"Sector4",1);
      wheelSigmaHistos[wh]->setBinLabel(5,"Sector5",1);
      wheelSigmaHistos[wh]->setBinLabel(6,"Sector6",1);
      wheelSigmaHistos[wh]->setBinLabel(7,"Sector7",1);
      wheelSigmaHistos[wh]->setBinLabel(8,"Sector8",1);
      wheelSigmaHistos[wh]->setBinLabel(9,"Sector9",1);
      wheelSigmaHistos[wh]->setBinLabel(10,"Sector10",1);
      wheelSigmaHistos[wh]->setBinLabel(11,"Sector11",1);
      wheelSigmaHistos[wh]->setBinLabel(12,"Sector12",1);
      wheelSigmaHistos[wh]->setBinLabel(13,"Sector13",1);
      wheelSigmaHistos[wh]->setBinLabel(14,"Sector14",1);
      wheelSigmaHistos[wh]->setBinLabel(1,"MB1_SL1",2);
      wheelSigmaHistos[wh]->setBinLabel(2,"MB1_SL2",2);
      wheelSigmaHistos[wh]->setBinLabel(3,"MB1_SL3",2);
      wheelSigmaHistos[wh]->setBinLabel(4,"MB2_SL1",2);
      wheelSigmaHistos[wh]->setBinLabel(5,"MB2_SL2",2);
      wheelSigmaHistos[wh]->setBinLabel(6,"MB2_SL3",2);
      wheelSigmaHistos[wh]->setBinLabel(7,"MB3_SL1",2);
      wheelSigmaHistos[wh]->setBinLabel(8,"MB3_SL2",2);
      wheelSigmaHistos[wh]->setBinLabel(9,"MB3_SL3",2);
      wheelSigmaHistos[wh]->setBinLabel(10,"MB4_SL1",2);
      wheelSigmaHistos[wh]->setBinLabel(11,"MB4_SL3",2);
    }  
  }

  if(parameters.getUntrackedParameter<bool>("slopeTest")){
    if(wheelSlopeHistos.find(wh) == wheelSlopeHistos.end()){
      string histoName =  "SlopeSummaryRes_testFailed_" + parameters.getUntrackedParameter<string>("STEP", "STEP3") + "_W" + wheel.str();
      wheelSlopeHistos[wh] = dbe->book2D(histoName.c_str(),histoName.c_str(),14,0,14,11,0,11);
      wheelSlopeHistos[wh]->setBinLabel(1,"Sector1",1);
      wheelSlopeHistos[wh]->setBinLabel(2,"Sector2",1);
      wheelSlopeHistos[wh]->setBinLabel(3,"Sector3",1);
      wheelSlopeHistos[wh]->setBinLabel(4,"Sector4",1);
      wheelSlopeHistos[wh]->setBinLabel(5,"Sector5",1);
      wheelSlopeHistos[wh]->setBinLabel(6,"Sector6",1);
      wheelSlopeHistos[wh]->setBinLabel(7,"Sector7",1);
      wheelSlopeHistos[wh]->setBinLabel(8,"Sector8",1);
      wheelSlopeHistos[wh]->setBinLabel(9,"Sector9",1);
      wheelSlopeHistos[wh]->setBinLabel(10,"Sector10",1);
      wheelSlopeHistos[wh]->setBinLabel(11,"Sector11",1);
      wheelSlopeHistos[wh]->setBinLabel(12,"Sector12",1);
      wheelSlopeHistos[wh]->setBinLabel(13,"Sector13",1);
      wheelSlopeHistos[wh]->setBinLabel(14,"Sector14",1);
      wheelSlopeHistos[wh]->setBinLabel(1,"MB1_SL1",2);
      wheelSlopeHistos[wh]->setBinLabel(2,"MB1_SL2",2);
      wheelSlopeHistos[wh]->setBinLabel(3,"MB1_SL3",2);
      wheelSlopeHistos[wh]->setBinLabel(4,"MB2_SL1",2);
      wheelSlopeHistos[wh]->setBinLabel(5,"MB2_SL2",2);
      wheelSlopeHistos[wh]->setBinLabel(6,"MB2_SL3",2);
      wheelSlopeHistos[wh]->setBinLabel(7,"MB3_SL1",2);
      wheelSlopeHistos[wh]->setBinLabel(8,"MB3_SL2",2);
      wheelSlopeHistos[wh]->setBinLabel(9,"MB3_SL3",2);
      wheelSlopeHistos[wh]->setBinLabel(10,"MB4_SL1",2);
      wheelSlopeHistos[wh]->setBinLabel(11,"MB4_SL3",2);
    }  
  }
  
}

void DTResolutionTest::endJob ( void  )

protectedvirtual

Endjob.

Reimplemented from edm::EDAnalyzer.

Definition at line 416 of file DTResolutionTest.cc.

References dumpDBToFile_GT_ttrig_cfg::outputFileName, and Parameters::parameters.

                             {

  edm::LogVerbatim ("resolution") << "[DTResolutionTest] endjob called!";
  //dbe->rmdir("DT/DTCalibValidation");
  //dbe->rmdir("DT/Tests/DTResolution");
  bool outputMEsInRootFile = parameters.getParameter<bool>("OutputMEsInRootFile");
  if(outputMEsInRootFile){
    std::string outputFileName = parameters.getParameter<std::string>("OutputFileName");
    dbe->save(outputFileName,"DT/CalibrationSummary");  
  } 

}

void DTResolutionTest::endLuminosityBlock ( edm::LuminosityBlock const &  lumiSeg, edm::EventSetup const &  c  )

protectedvirtual

DQM Client Diagnostic.

Reimplemented from edm::EDAnalyzer.

Definition at line 109 of file DTResolutionTest.cc.

References abs, QReport::getBadChannels(), MonitorElement::getQReport(), MonitorElement::getTH2F(), interpolateCardsSimple::histo, i, edm::LuminosityBlockBase::id(), j, edm::LuminosityBlockID::luminosityBlock(), python.multivaluedict::map(), timingPdfMaker::mean, Parameters::parameters, dtDQMClient_cfg::prescaleFactor, DTChamberId::sector(), slope, DTChamberId::station(), relativeConstraints::station, DTSuperLayerId::superLayer(), DTSuperLayerId::superlayer(), and DTChamberId::wheel().

                                                                                                   {
  
  // counts number of updats (online mode) or number of events (standalone mode)
  //nevents++;
  // if running in standalone perform diagnostic only after a reasonalbe amount of events
  //if ( parameters.getUntrackedParameter<bool>("runningStandalone", false) && 
  //   nevents%parameters.getUntrackedParameter<int>("diagnosticPrescale", 1000) != 0 ) return;
  //edm::LogVerbatim ("resolution") << "[DTResolutionTest]: "<<nevents<<" updates";

  edm::LogVerbatim ("resolution") <<"[DTResolutionTest]: End of LS transition, performing the DQM client operation";

  // counts number of lumiSegs 
  nLumiSegs = lumiSeg.id().luminosityBlock();

  // prescale factor
  if ( nLumiSegs%prescaleFactor != 0 ) return;

  for(map<int, MonitorElement*> ::const_iterator histo = wheelMeanHistos.begin();
      histo != wheelMeanHistos.end();
      histo++) {
    (*histo).second->Reset();
  }
  if(parameters.getUntrackedParameter<bool>("sigmaTest")){
    for(map<int, MonitorElement*> ::const_iterator histo = wheelSigmaHistos.begin();
    histo != wheelSigmaHistos.end();
    histo++) {
      (*histo).second->Reset();
    }
  }
  if(parameters.getUntrackedParameter<bool>("slopeTest")){
    for(map<int, MonitorElement*> ::const_iterator histo = wheelSlopeHistos.begin();
    histo != wheelSlopeHistos.end();
    histo++) {
      (*histo).second->Reset();
    }
  }

  cmsMeanHistos.clear();
  for(int i=-2; i<3; i++){
    for(int j=1; j<15; j++){
      MeanFilled[make_pair(i,j)]=false;
    }
  }
  if(parameters.getUntrackedParameter<bool>("sigmaTest")){
    cmsSigmaHistos.clear();
    for(int i=-2; i<3; i++){
      for(int j=1; j<15; j++){
    SigmaFilled[make_pair(i,j)]=false;
      }
    }
  }
  if(parameters.getUntrackedParameter<bool>("slopeTest")){
    cmsSlopeHistos.clear();
    for(int i=-2; i<3; i++){
      for(int j=1; j<15; j++){
    SlopeFilled[make_pair(i,j)]=false;
      }
    }
  }


  // summary histos initialization
  for(int wh=-2; wh<=3; wh++){
    if(wh!=3){
      for(int xBin=0; xBin<14; xBin++){
    for(int yBin=0; yBin<11; yBin++){
      wheelMeanHistos[wh]->setBinContent(xBin,yBin,0);
      if(parameters.getUntrackedParameter<bool>("sigmaTest"))
        wheelSigmaHistos[wh]->setBinContent(xBin,yBin,0);
      if(parameters.getUntrackedParameter<bool>("slopeTest"))
        wheelSlopeHistos[wh]->setBinContent(xBin,yBin,0);
      }
      }
    }
    else{
      for(int xBin=0; xBin<14; xBin++){
    for(int yBin=-2; yBin<3; yBin++){
      wheelMeanHistos[wh]->setBinContent(xBin,yBin,0);
      if(parameters.getUntrackedParameter<bool>("sigmaTest"))
        wheelSigmaHistos[wh]->setBinContent(xBin,yBin,0);
      if(parameters.getUntrackedParameter<bool>("slopeTest"))
        wheelSlopeHistos[wh]->setBinContent(xBin,yBin,0);
    }
      }
    }
  }


  edm::LogVerbatim ("resolution") <<"[DTResolutionTest]: "<<nLumiSegs<<" updates";

  vector<DTChamber*>::const_iterator ch_it = muonGeom->chambers().begin();
  vector<DTChamber*>::const_iterator ch_end = muonGeom->chambers().end();

  edm::LogVerbatim ("resolution") << "[DTResolutionTest]: Residual Distribution tests results";
  
  for (; ch_it != ch_end; ++ch_it) {

    DTChamberId chID = (*ch_it)->id();

    // Fill the test histos
    int entry=-1;
    if(chID.station() == 1) entry=0;
    if(chID.station() == 2) entry=3;
    if(chID.station() == 3) entry=6;
    if(chID.station() == 4) entry=9;

    vector<const DTSuperLayer*>::const_iterator sl_it = (*ch_it)->superLayers().begin(); 
    vector<const DTSuperLayer*>::const_iterator sl_end = (*ch_it)->superLayers().end();

    for(; sl_it != sl_end; ++sl_it) {

      DTSuperLayerId slID = (*sl_it)->id();

      edm::LogVerbatim ("resolution") << "[DTResolutionTest]: Superlayer: " << slID;

      stringstream wheel; wheel << slID.wheel();    
      stringstream station; station << slID.station();  
      stringstream sector; sector << slID.sector(); 
      stringstream superLayer; superLayer << slID.superlayer();
      
      string HistoName = "W" + wheel.str() + "_Sec" + sector.str(); 
      string supLayer = "W" + wheel.str() + "_St" + station.str() + "_Sec" + sector.str() + "_SL" + superLayer.str(); 

      MonitorElement * res_histo = dbe->get(getMEName(slID));
      if (res_histo) {
    // gaussian test
    string GaussianCriterionName = 
      parameters.getUntrackedParameter<string>("resDistributionTestName",
                           "ResidualsDistributionGaussianTest");
    const QReport * GaussianReport = res_histo->getQReport(GaussianCriterionName);
    if(GaussianReport){
      // FIXE ME: if the quality test fails this cout return a null pointer
      //edm::LogWarning ("resolution") << "-------- SuperLayer : "<<supLayer<<"  "<<GaussianReport->getMessage()<<" ------- "<<GaussianReport->getStatus();
    }
    int BinNumber = entry+slID.superLayer();
    if(BinNumber == 12) BinNumber=11;
    float mean = (*res_histo).getMean(1);
    float sigma = (*res_histo).getRMS(1);
    MeanHistos.find(make_pair(slID.wheel(),slID.sector()))->second->setBinContent(BinNumber, mean); 
    if(parameters.getUntrackedParameter<bool>("sigmaTest"))
      SigmaHistos.find(make_pair(slID.wheel(),slID.sector()))->second->setBinContent(BinNumber, sigma);
      }

      if(parameters.getUntrackedParameter<bool>("slopeTest")){
    MonitorElement * res_histo_2D = dbe->get(getMEName2D(slID));
    if (res_histo_2D) {
      TH2F * res_histo_2D_root = res_histo_2D->getTH2F();
      int BinNumber = entry+slID.superLayer();
      if(BinNumber == 12) BinNumber=11;
      TProfile* prof = res_histo_2D_root->ProfileX();
      prof->GetXaxis()->SetRangeUser(0,2);
      //prof->Fit("pol1","Q0");
      try {
        prof->Fit("pol1","Q0");
      } catch (...) {
        edm::LogError ("resolution") << "[DTResolutionTest]: Exception when fitting..."
                     << "SuperLayer : " << slID << "\n"
                     << "                    STEP : " << parameters.getUntrackedParameter<string>("STEP", "STEP3") << "\n"      
                     << "Filling slope histogram with standard value -99. for bin " << BinNumber;
        SlopeHistos.find(make_pair(slID.wheel(),slID.sector()))->second->setBinContent(BinNumber, -99.);
        continue;
      }
      TF1 *fitting = prof->GetFunction("pol1");
      double slope = fitting->GetParameter(1);
      SlopeHistos.find(make_pair(slID.wheel(),slID.sector()))->second->setBinContent(BinNumber, slope); 
    }
      }

    }
  }

  // Mean test 
  string MeanCriterionName = parameters.getUntrackedParameter<string>("meanTestName","ResidualsMeanInRange"); 
  for(map<pair<int,int>, MonitorElement*>::const_iterator hMean = MeanHistos.begin();
      hMean != MeanHistos.end();
      hMean++) {
    const QReport * theMeanQReport = (*hMean).second->getQReport(MeanCriterionName);
    stringstream wheel; wheel << (*hMean).first.first;
    stringstream sector; sector << (*hMean).first.second;
    // Report the channels failing the test on the mean
    if(theMeanQReport) { 
      vector<dqm::me_util::Channel> badChannels = theMeanQReport->getBadChannels();
      for (vector<dqm::me_util::Channel>::iterator channel = badChannels.begin(); 
       channel != badChannels.end(); channel++) {
    edm::LogError("resolution") << "Bad mean channel: wh: " << wheel.str()
                     << " st: " << stationFromBin((*channel).getBin())
                     << " sect: " <<sector.str()
                     << " sl: " << slFromBin((*channel).getBin())
                     << " mean (cm): " << (*channel).getContents();
    string HistoName = "W" + wheel.str() + "_Sec" + sector.str();
    if(parameters.getUntrackedParameter<bool>("meanWrongHisto")){
      MeanHistosSetRange.find(HistoName)->second->Fill((*channel).getBin());
      MeanHistosSetRange2D.find(HistoName)->second->Fill((*channel).getBin(),(*channel).getContents());
    }
    // fill the wheel summary histos if the SL has not passed the test
    if(abs((*channel).getContents())<parameters.getUntrackedParameter<double>("meanMaxLimit"))
      wheelMeanHistos[(*hMean).first.first]->Fill(((*hMean).first.second)-1,(*channel).getBin()-1,1);
    else
      wheelMeanHistos[(*hMean).first.first]->Fill(((*hMean).first.second)-1,(*channel).getBin()-1,2);
    // fill the cms summary histo if the percentual of SL which have not passed the test 
    // is more than a predefined treshold
    if(abs((*channel).getContents())>parameters.getUntrackedParameter<double>("meanMaxLimit")){
      cmsMeanHistos[make_pair((*hMean).first.first,(*hMean).first.second)]++;
      if(((*hMean).first.second<13 &&
          double(cmsMeanHistos[make_pair((*hMean).first.first,(*hMean).first.second)])/11>double(percentual)/100 &&
          MeanFilled[make_pair((*hMean).first.first,(*hMean).first.second)]==false) ||
         ((*hMean).first.first>=13 && 
        double(cmsMeanHistos[make_pair((*hMean).first.first,(*hMean).first.second)])/2>double(percentual)/100 &&
          MeanFilled[make_pair((*hMean).first.first,(*hMean).first.second)]==false)){
        MeanFilled[make_pair((*hMean).first.first,(*hMean).first.second)]=true;
        wheelMeanHistos[3]->Fill(((*hMean).first.second)-1,(*hMean).first.first);
      } 
    }
      }
      // FIXE ME: if the quality test fails this cout return a null pointer
      //edm::LogWarning ("resolution") << "-------- wheel: "<<wheel.str()<<" sector: "<<sector.str()<<"  "<<theMeanQReport->getMessage()<<" ------- "<<theMeanQReport->getStatus(); 
    }
  }
  
  // Sigma test
  if(parameters.getUntrackedParameter<bool>("sigmaTest")){
    string SigmaCriterionName = parameters.getUntrackedParameter<string>("sigmaTestName","ResidualsSigmaInRange"); 
    for(map<pair<int,int>, MonitorElement*>::const_iterator hSigma = SigmaHistos.begin();
    hSigma != SigmaHistos.end();
    hSigma++) {
      const QReport * theSigmaQReport = (*hSigma).second->getQReport(SigmaCriterionName);
      stringstream wheel; wheel << (*hSigma).first.first;
      stringstream sector; sector << (*hSigma).first.second;
      if(theSigmaQReport) {
    vector<dqm::me_util::Channel> badChannels = theSigmaQReport->getBadChannels();
    for (vector<dqm::me_util::Channel>::iterator channel = badChannels.begin(); 
         channel != badChannels.end(); channel++) {
      edm::LogError("resolution") << "Bad sigma: wh: " << wheel.str()
                      << " st: " << stationFromBin((*channel).getBin())
                      << " sect: " <<sector.str()
                      << " sl: " << slFromBin((*channel).getBin())
                      << " sigma (cm): " << (*channel).getContents();
      string HistoName = "W" + wheel.str() + "_Sec" + sector.str();
      SigmaHistosSetRange.find(HistoName)->second->Fill((*channel).getBin());
      SigmaHistosSetRange2D.find(HistoName)->second->Fill((*channel).getBin(),(*channel).getContents());
      // fill the wheel summary histos if the SL has not passed the test
      wheelSigmaHistos[(*hSigma).first.first]->Fill(((*hSigma).first.second)-1,(*channel).getBin()-1);
      // fill the cms summary histo if the percentual of SL which have not passed the test 
      // is more than a predefined treshold
      cmsSigmaHistos[make_pair((*hSigma).first.first,(*hSigma).first.second)]++;
      if(((*hSigma).first.second<13 &&
          double(cmsSigmaHistos[make_pair((*hSigma).first.first,(*hSigma).first.second)])/11>double(percentual)/100 &&
        SigmaFilled[make_pair((*hSigma).first.first,(*hSigma).first.second)]==false) ||
         ((*hSigma).first.first>=13 && 
          double(cmsSigmaHistos[make_pair((*hSigma).first.first,(*hSigma).first.second)])/2>double(percentual)/100 &&
          SigmaFilled[make_pair((*hSigma).first.first,(*hSigma).first.second)]==false)){
        SigmaFilled[make_pair((*hSigma).first.first,(*hSigma).first.second)]=true;
        wheelSigmaHistos[3]->Fill((*hSigma).first.second-1,(*hSigma).first.first);
      }
    }
    // FIXE ME: if the quality test fails this cout return a null pointer
    //edm::LogWarning ("resolution") << "-------- wheel: "<<wheel.str()<<" sector: "<<sector.str()<<"  "<<theSigmaQReport->getMessage()<<" ------- "<<theSigmaQReport->getStatus();
      }
    }
  }

  // Slope test
  if(parameters.getUntrackedParameter<bool>("slopeTest")){
    string SlopeCriterionName = parameters.getUntrackedParameter<string>("slopeTestName","ResidualsSlopeInRange"); 
    for(map<pair<int,int>, MonitorElement*>::const_iterator hSlope = SlopeHistos.begin();
    hSlope != SlopeHistos.end();
    hSlope++) {
      const QReport * theSlopeQReport = (*hSlope).second->getQReport(SlopeCriterionName);
      stringstream wheel; wheel << (*hSlope).first.first;
      stringstream sector; sector << (*hSlope).first.second;
      if(theSlopeQReport) {
    vector<dqm::me_util::Channel> badChannels = theSlopeQReport->getBadChannels();
    for (vector<dqm::me_util::Channel>::iterator channel = badChannels.begin(); 
         channel != badChannels.end(); channel++) {
      edm::LogError("resolution") << "Bad slope: wh: " << wheel.str()
                      << " st: " << stationFromBin((*channel).getBin())
                      << " sect: " <<sector.str()
                      << " sl: " << slFromBin((*channel).getBin())
                      << " slope: " << (*channel).getContents();
      string HistoName = "W" + wheel.str() + "_Sec" + sector.str();
      SlopeHistosSetRange.find(HistoName)->second->Fill((*channel).getBin());
      SlopeHistosSetRange2D.find(HistoName)->second->Fill((*channel).getBin(),(*channel).getContents());
      // fill the wheel summary histos if the SL has not passed the test
      wheelSlopeHistos[(*hSlope).first.first]->Fill(((*hSlope).first.second)-1,(*channel).getBin()-1);
      // fill the cms summary histo if the percentual of SL which have not passed the test 
      // is more than a predefined treshold
      cmsSlopeHistos[make_pair((*hSlope).first.first,(*hSlope).first.second)]++;
      if(((*hSlope).first.second<13 &&
          double(cmsSlopeHistos[make_pair((*hSlope).first.first,(*hSlope).first.second)])/11>double(percentual)/100 &&
          SlopeFilled[make_pair((*hSlope).first.first,(*hSlope).first.second)]==false) ||
         ((*hSlope).first.first>=13 && 
          double(cmsSlopeHistos[make_pair((*hSlope).first.first,(*hSlope).first.second)])/2>double(percentual)/100 &&
          SlopeFilled[make_pair((*hSlope).first.first,(*hSlope).first.second)]==false)){
        SlopeFilled[make_pair((*hSlope).first.first,(*hSlope).first.second)]=true;
        wheelSlopeHistos[3]->Fill((*hSlope).first.second-1,(*hSlope).first.first);
      }
    }
    // FIXE ME: if the quality test fails this cout return a null pointer
    //edm::LogWarning ("resolution") << "-------- wheel: "<<wheel.str()<<" sector: "<<sector.str()<<"  "<<theSlopeQReport->getMessage()<<" ------- "<<theSlopeQReport->getStatus();
      }
    }
  }

}

string DTResolutionTest::getMEName ( const DTSuperLayerId &  slID )

protected

Get the ME name.

Definition at line 431 of file DTResolutionTest.cc.

References Parameters::parameters, DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, DTSuperLayerId::superlayer(), and DTChamberId::wheel().

                                                              {
  
  stringstream wheel; wheel << slID.wheel();    
  stringstream station; station << slID.station();  
  stringstream sector; sector << slID.sector(); 
  stringstream superLayer; superLayer << slID.superlayer();
  
  string folderRoot = parameters.getUntrackedParameter<string>("folderRoot", "Collector/FU0/");
  string folderName;

  if(parameters.getUntrackedParameter<bool>("calibModule", false)){
    folderName = 
      folderRoot + "DT/DTCalibValidation/Wheel" +  wheel.str() +
      "/Station" + station.str() +
      "/Sector" + sector.str() + "/";
  }
  else{
    folderName = 
      folderRoot + "DT/DTResolutionAnalysisTask/Wheel" +  wheel.str() +
      "/Station" + station.str() +
      "/Sector" + sector.str() + "/";
  }

  string histoTag = parameters.getUntrackedParameter<string>("histoTag", "hResDist");

  string histoname = folderName + histoTag  
    + "_W" + wheel.str() 
    + "_St" + station.str() 
    + "_Sec" + sector.str() 
    + "_SL" + superLayer.str(); 
  
  return histoname;
  
}

string DTResolutionTest::getMEName2D ( const DTSuperLayerId &  slID )

protected

Definition at line 467 of file DTResolutionTest.cc.

References Parameters::parameters, DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, DTSuperLayerId::superlayer(), and DTChamberId::wheel().

                                                                {
  
  stringstream wheel; wheel << slID.wheel();    
  stringstream station; station << slID.station();  
  stringstream sector; sector << slID.sector(); 
  stringstream superLayer; superLayer << slID.superlayer();
  
  string folderRoot = parameters.getUntrackedParameter<string>("folderRoot", "Collector/FU0/");
  string folderName;

  if(parameters.getUntrackedParameter<bool>("calibModule", false)){
    folderName = 
      folderRoot + "DT/DTCalibValidation/Wheel" +  wheel.str() +
      "/Station" + station.str() +
      "/Sector" + sector.str() + "/";
  }
  else{
    folderName = 
      folderRoot + "DT/DTResolutionAnalysisTask/Wheel" +  wheel.str() +
      "/Station" + station.str() +
      "/Sector" + sector.str() + "/";
  }

  string histoTag2D = parameters.getUntrackedParameter<string>("histoTag2D", "hResDistVsDist");

  string histoname = folderName + histoTag2D  
    + "_W" + wheel.str() 
    + "_St" + station.str() 
    + "_Sec" + sector.str() 
    + "_SL" + superLayer.str(); 
  
  return histoname;
  
}

int DTResolutionTest::slFromBin ( int  bin ) const

private

Definition at line 775 of file DTResolutionTest.cc.

References run_regression::ret.

                                             {
  int ret = bin%3;
  if(ret == 0 || bin == 11) ret = 3;
  
  return ret;
}

int DTResolutionTest::stationFromBin ( int  bin ) const

private

Definition at line 770 of file DTResolutionTest.cc.

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

Member Data Documentation

std::map<std::pair<int,int>, int> DTResolutionTest::cmsMeanHistos

private

Definition at line 109 of file DTResolutionTest.h.

std::map<std::pair<int,int>, int> DTResolutionTest::cmsSigmaHistos

private

Definition at line 111 of file DTResolutionTest.h.

std::map<std::pair<int,int>, int> DTResolutionTest::cmsSlopeHistos

private

Definition at line 113 of file DTResolutionTest.h.

DQMStore* DTResolutionTest::dbe

private

Definition at line 87 of file DTResolutionTest.h.

std::map<std::pair<int,int>, bool> DTResolutionTest::MeanFilled

private

Definition at line 110 of file DTResolutionTest.h.

std::map< std::pair<int,int> , MonitorElement* > DTResolutionTest::MeanHistos

private

Definition at line 93 of file DTResolutionTest.h.

std::map< std::string , MonitorElement* > DTResolutionTest::MeanHistosSetRange

private

Definition at line 96 of file DTResolutionTest.h.

std::map< std::string , MonitorElement* > DTResolutionTest::MeanHistosSetRange2D

private

Definition at line 99 of file DTResolutionTest.h.

edm::ESHandle<DTGeometry> DTResolutionTest::muonGeom

private

Definition at line 90 of file DTResolutionTest.h.

int DTResolutionTest::nevents

private

Definition at line 81 of file DTResolutionTest.h.

unsigned int DTResolutionTest::nLumiSegs

private

Definition at line 82 of file DTResolutionTest.h.

edm::ParameterSet DTResolutionTest::parameters

private

Definition at line 89 of file DTResolutionTest.h.

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputTags(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::recursePSetProperties().

int DTResolutionTest::percentual

private

Definition at line 85 of file DTResolutionTest.h.

int DTResolutionTest::prescaleFactor

private

Definition at line 83 of file DTResolutionTest.h.

int DTResolutionTest::run

private

Definition at line 84 of file DTResolutionTest.h.

Referenced by generateEDF.LumiInfo::__str__(), Types.EventID::cppID(), Types.LuminosityBlockID::cppID(), generateEDF.LumiInfo::fixXingInfo(), and dqm_interfaces.DirWalkerFile::walk().

std::map<std::pair<int,int>, bool> DTResolutionTest::SigmaFilled

private

Definition at line 112 of file DTResolutionTest.h.

std::map< std::pair<int,int> , MonitorElement* > DTResolutionTest::SigmaHistos

private

Definition at line 94 of file DTResolutionTest.h.

std::map< std::string , MonitorElement* > DTResolutionTest::SigmaHistosSetRange

private

Definition at line 97 of file DTResolutionTest.h.

std::map< std::string , MonitorElement* > DTResolutionTest::SigmaHistosSetRange2D

private

Definition at line 100 of file DTResolutionTest.h.

std::map<std::pair<int,int>, bool> DTResolutionTest::SlopeFilled

private

Definition at line 114 of file DTResolutionTest.h.

std::map< std::pair<int,int> , MonitorElement* > DTResolutionTest::SlopeHistos

private

Definition at line 95 of file DTResolutionTest.h.

std::map< std::string , MonitorElement* > DTResolutionTest::SlopeHistosSetRange

private

Definition at line 98 of file DTResolutionTest.h.

std::map< std::string , MonitorElement* > DTResolutionTest::SlopeHistosSetRange2D

private

Definition at line 101 of file DTResolutionTest.h.

std::map< int, MonitorElement* > DTResolutionTest::wheelMeanHistos

private

Definition at line 104 of file DTResolutionTest.h.

std::map< int, MonitorElement* > DTResolutionTest::wheelSigmaHistos

private

Definition at line 105 of file DTResolutionTest.h.

std::map< int, MonitorElement* > DTResolutionTest::wheelSlopeHistos

private

Definition at line 106 of file DTResolutionTest.h.