DTFineDelayCorr Class Reference

#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.

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.