#include <DTFineDelayCorr.h>

Inheritance diagram for DTFineDelayCorr:

Public Member Functions
	DTFineDelayCorr (const edm::ParameterSet &ps)
	Constructor. More...

virtual	~DTFineDelayCorr ()
	Destructor. More...

Public Member Functions inherited from DTLocalTriggerBaseTest
	DTLocalTriggerBaseTest ()
	Constructor. More...

virtual	~DTLocalTriggerBaseTest ()
	Destructor. More...

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Protected Member Functions
void	beginJob ()
	Begin Job. More...

void	beginRun (const edm::Run &run, const edm::EventSetup &evSU)
	Begin Run. More...

void	endJob ()
	End Job. More...

void	runClientDiagnostic ()
	DQM Client Diagnostic. More...

void	writeConsts (const std::string &outputFileName) const

Protected Member Functions inherited from DTLocalTriggerBaseTest
void	analyze (const edm::Event &e, const edm::EventSetup &c)
	Analyze. More...

void	beginJob ()
	BeginJob. More...

void	beginLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &context)
	Perform begin lumiblock operations. More...

void	beginRun (edm::Run const &run, edm::EventSetup const &context)
	BeginRun. More...

void	bookCmsHistos (std::string hTag, std::string folder="", bool isGlb=false)
	Book the new MEs (CMS summary) More...

void	bookSectorHistos (int wheel, int sector, std::string hTag, std::string folder="")
	Book the new MEs (for each sector) More...

void	bookWheelHistos (int wheel, std::string hTag, std::string folder="")
	Book the new MEs (for each wheel) More...

std::string	category ()
	Get message logger name. More...

void	endJob ()
	EndJob. More...

void	endLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &context)
	Perform client diagnostic in online. More...

void	endRun (edm::Run const &run, edm::EventSetup const &context)
	Perform client diagnostic in offline. More...

std::string	fullName (std::string htype)
	Create fullname from histo partial name. More...

template<class T >
T *	getHisto (MonitorElement *me)
	Convert ME to Histogram fo type T. More...

std::string	getMEName (std::string histoTag, std::string subfolder, const DTChamberId &chambid)
	Get the ME name (by chamber) More...

std::string	getMEName (std::string histoTag, std::string subfolder, int wh)
	Get the ME name (by wheel) More...

std::pair< float, float >	phiRange (const DTChamberId &id)
	Calculate phi range for histograms. More...

void	setConfig (const edm::ParameterSet &ps, std::string name)
	Set configuration variables. More...

std::string &	topFolder (bool isDCC)
	Get top folder name. More...

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Private Attributes
std::map< DTChamberId, std::vector< float > >	delayMap

edm::ESHandle< DTConfigManager >	dtConfig

bool	gaussMean

int	minEntries

int	nEvents

std::map< DTChamberId, std::pair< int, float > >	oldDelayMap

std::string	oldDelaysInputFile

std::string	outputFileName

bool	readOldFromDb

std::string	t0MeanHistoTag

std::string	trSource

edm::ESHandle< DTTPGParameters >	worstPhaseMap

bool	writeDB

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Attributes inherited from DTLocalTriggerBaseTest
std::string	baseFolderDCC

std::string	baseFolderDDU

std::map< std::string, MonitorElement * >	cmsME

DQMStore *	dbe

std::string	hwSource

std::vector< std::string >	hwSources

edm::ESHandle< DTGeometry >	muonGeom

int	nevents

unsigned int	nLumiSegs

edm::ParameterSet	parameters

int	prescaleFactor

int	run

bool	runOnline

std::map< int, std::map < std::string, MonitorElement * > >	secME

std::string	sourceFolder

std::string	testName

std::string	trigSource

std::vector< std::string >	trigSources

std::map< int, std::map < std::string, MonitorElement * > >	whME

Detailed Description

Class to calculate corrections to the fine delays loaded in DT MiniCrates

Author
M. Giunta, C. Battilana

Definition at line 26 of file DTFineDelayCorr.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 43 of file DTFineDelayCorr.cc.

                                                        {
 
   setConfig(ps,"DTFineDelayCorr");  // sets parameter values and name used in log file 
   baseFolderDCC = "DT/90-LocalTriggerSynch/";
   baseFolderDDU = "DT/90-LocalTriggerSynch/";
  
 }

DTFineDelayCorr::~DTFineDelayCorr ( )

virtual

Destructor.

Definition at line 52 of file DTFineDelayCorr.cc.

                                  {
 
 }

Member Function Documentation

void DTFineDelayCorr::beginJob ( void )

protectedvirtual

Begin Job.

Reimplemented from edm::EDAnalyzer.

Definition at line 57 of file DTFineDelayCorr.cc.

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

                               {
 
   // Tag for Hardware Source (DDU or DCC)
   hwSource = parameters.getParameter<string>("hwSource");
   // Tag for the t0Mean Histograms
   t0MeanHistoTag = parameters.getParameter<string>("t0MeanHistoTag");
   // Read old delays from file or from Db
   readOldFromDb = parameters.getParameter<bool>("readOldFromDb");
   // Input file name for old delays
   oldDelaysInputFile = parameters.getParameter<string>("oldDelaysInputFile"),
   // Write new delays to file or to Db
   writeDB = parameters.getParameter<bool>("writeDB");
   // Output File Name
   outputFileName = parameters.getParameter<string>("outputFile");
   // Choose to use Hist Mean or Gaussian Fit Mean
   gaussMean =  parameters.getParameter<bool>("gaussMean");
   // Require Minimum Number Of Entries in the t0Mean Histogram
   minEntries =  parameters.getUntrackedParameter<int>("minEntries",5);
 
 }

void DTFineDelayCorr::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	evSU
	)

protectedvirtual

Begin Run.

Reimplemented from edm::EDAnalyzer.

Definition at line 78 of file DTFineDelayCorr.cc.

References DTLocalTriggerBaseTest::beginRun(), and edm::EventSetup::get().

                                                                     {
 
   DTLocalTriggerBaseTest::beginRun(run,evSU);
   evSU.get< DTConfigManagerRcd >().get(dtConfig);
   evSU.get< DTTPGParametersRcd >().get(worstPhaseMap);
 
 }

void DTFineDelayCorr::endJob ( void )

protectedvirtual

End Job.

Reimplemented from edm::EDAnalyzer.

Definition at line 185 of file DTFineDelayCorr.cc.

References filterCSVwithJSON::copy, DTLocalTriggerBaseTest::endJob(), getDQMSummary::iter, python.multivaluedict::map(), GetRecoTauVFromDQM_MC_cff::outFile, and dumpDBToFile_GT_ttrig_cfg::outputFileName.

                             {
 
   DTLocalTriggerBaseTest::endJob();
 
    if (writeDB) {
      // to be added if needed
    }
    else { // write txt file
      // ** Open output file **
      ofstream outFile(outputFileName.c_str());
      for(map< DTChamberId,  vector<float> >::const_iterator iter = delayMap.begin();
      iter != delayMap.end() ; iter++) {
        // writing
        ostream_iterator<float> oit(outFile, " ");    
        copy((*iter).second.begin(), (*iter).second.end(), oit);
        outFile << endl;
      }   
      outFile.close();
    }
 }

void DTFineDelayCorr::runClientDiagnostic ( )

protectedvirtual

DQM Client Diagnostic.

Implements DTLocalTriggerBaseTest.

Definition at line 86 of file DTFineDelayCorr.cc.

References python.rootplot.argparse::category, gather_cfg::cout, DTChamberId, DTLayerId, DTTopology::firstChannel(), DTConfigPedestals::getOffset(), DTLayer::id(), geometryCSVtoXML::line, timingPdfMaker::mean, TrigPrimClient_cfi::minEntries, DTTimeUnits::ns, DetId::rawId(), DTChamberId::sector(), DTLayer::specificTopology(), DTChamberId::station(), relativeConstraints::station, and DTChamberId::wheel().

                                           {
   int coarseDelay = -999;
   float oldFineDelay = -999;
   if(!readOldFromDb) { // read old delays from txt file
     // **  Open and read old delays input file  ** 
     ifstream oldDelaysFile(oldDelaysInputFile.c_str());
     string line;
 
     while (getline(oldDelaysFile, line)) {
       if( line == "" || line[0] == '#' ) continue; 
       stringstream linestr;
       int wheelKey,sectorKey, stationKey;      
       linestr << line;
       
       linestr >> wheelKey
           >> sectorKey
           >> stationKey
           >> coarseDelay
           >> oldFineDelay;
       
       pair<int,float> oldDelays = make_pair(coarseDelay,oldFineDelay);
       DTChamberId oldDelayKey = DTChamberId(wheelKey,stationKey,sectorKey);
       oldDelayMap.insert(make_pair(oldDelayKey,oldDelays));
     }
   }
 
   //  ** Loop over the chambers ** 
   vector<const DTChamber*>::const_iterator chambIt  = muonGeom->chambers().begin();
   vector<const DTChamber*>::const_iterator chambEnd = muonGeom->chambers().end();
   for (; chambIt!=chambEnd; ++chambIt) { 
     DTChamberId chId = (*chambIt)->id();
     uint32_t indexCh = chId.rawId();
     int wheel = chId.wheel();
     int sector = chId.sector();
     int station = chId.station();
     
     // ** Compute corrected values and write them to file or database **
     vector<float> newDelays;
 
     // **  Retrieve Delays Loaded in MiniCrates ** 
     if(readOldFromDb) {    // read from db 
       const DTConfigPedestals *pedestals = dtConfig->getDTConfigPedestals();
       const DTLayer *layer = muonGeom->layer(DTLayerId(chId,1,1));
       float delay = pedestals->getOffset(DTWireId(layer->id(),layer->specificTopology().firstChannel())); 
       coarseDelay = int(delay/25.);
       oldFineDelay = delay - coarseDelay * 25.;
     }
     else {                 // read from map created from txt file
       coarseDelay = oldDelayMap[chId].first;
       oldFineDelay = oldDelayMap[chId].second;
     }
 
     // ** Retrieve t0Mean histograms **
     TH1F *t0H = getHisto<TH1F>(dbe->get(getMEName(t0MeanHistoTag,"", chId)));
     float newFineDelay = -999;   // initialize to dummy number
     cout <<"MG: " << getMEName(t0MeanHistoTag,"", chId) << " entries: " << t0H->GetEntries() << endl; 
     if (t0H->GetEntries() > minEntries) {
       Double_t mean;
       // ** Find Mean Value of the distribution ** 
       if(gaussMean) {
     TF1 *funct = t0H->GetFunction("gaus");
     mean = funct->GetParameter(1);
       }
       else {
     mean = t0H->GetMean();
       }
       
       // ** Retrieve Worst Phase values **
       int wpCoarseDelay;     
       float wpFineDelay;
       worstPhaseMap->get(chId, wpCoarseDelay,  wpFineDelay, DTTimeUnits::ns);
 //       cout << "wpFineDelay, oldFineDelay, mean: " << wpFineDelay << " " 
 //     << oldFineDelay << " " << mean << endl;
       float bpFineDelay = (wpFineDelay < 12.5)? (wpFineDelay + 12.5) : (wpFineDelay - 12.5);  // Best Phase: half BX far from the worst phase 
       // ** Calculate correction **
       float diffFineDelays  = oldFineDelay + (mean - bpFineDelay); // positive mean shift implies positive delay correction
       int bxDiff = (int) (diffFineDelays / 25);
       coarseDelay += bxDiff;
       newFineDelay = fmodf(diffFineDelays, 25);
 //       cout << "diffFineDelays, newFineDelay, bxDiff, coarseDelay: " << diffFineDelays 
 //     << " "<< newFineDelay << " " << bxDiff << " " << coarseDelay << endl;
     }
     else {
       LogProblem(category()) << "[" << testName << "Test]:  Not enough entries in hist for Chamber "  
                  << indexCh << endl;
     }
 
     newDelays.push_back(wheel);
     newDelays.push_back(sector);
     newDelays.push_back(station);
     newDelays.push_back(coarseDelay);
     newDelays.push_back(newFineDelay);    
     pair< DTChamberId,  vector<float> > chDelays;
     chDelays.first = chId; 
     chDelays.second = newDelays;
     delayMap.insert(chDelays);
    }
 }

void DTFineDelayCorr::writeConsts ( const std::string & outputFileName ) const

protected

Member Data Documentation

std::map< DTChamberId, std::vector<float> > DTFineDelayCorr::delayMap

private

Definition at line 70 of file DTFineDelayCorr.h.

edm::ESHandle< DTConfigManager > DTFineDelayCorr::dtConfig

private

Definition at line 63 of file DTFineDelayCorr.h.

bool DTFineDelayCorr::gaussMean

private

Definition at line 60 of file DTFineDelayCorr.h.

int DTFineDelayCorr::minEntries

private

Definition at line 61 of file DTFineDelayCorr.h.

int DTFineDelayCorr::nEvents

private

Definition at line 62 of file DTFineDelayCorr.h.

Referenced by looper.Looper::loop().

std::map< DTChamberId, std::pair<int,float> > DTFineDelayCorr::oldDelayMap

private

Definition at line 67 of file DTFineDelayCorr.h.

std::string DTFineDelayCorr::oldDelaysInputFile

private

Definition at line 56 of file DTFineDelayCorr.h.

std::string DTFineDelayCorr::outputFileName

private

Definition at line 55 of file DTFineDelayCorr.h.

bool DTFineDelayCorr::readOldFromDb

private

Definition at line 58 of file DTFineDelayCorr.h.

std::string DTFineDelayCorr::t0MeanHistoTag

private

Definition at line 54 of file DTFineDelayCorr.h.

std::string DTFineDelayCorr::trSource

private

Definition at line 57 of file DTFineDelayCorr.h.

edm::ESHandle< DTTPGParameters > DTFineDelayCorr::worstPhaseMap

private

Definition at line 64 of file DTFineDelayCorr.h.

bool DTFineDelayCorr::writeDB

private

Definition at line 59 of file DTFineDelayCorr.h.

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation