DTLocalTriggerSynchTask Class Reference

#include <DTLocalTriggerSynchTask.h>

Inheritance diagram for DTLocalTriggerSynchTask:

Public Member Functions
	DTLocalTriggerSynchTask (const edm::ParameterSet &ps)
	Constructor. More...
virtual	~DTLocalTriggerSynchTask ()
	Destructor. More...
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Protected Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &context)
	Analyze. More...
std::string &	baseDir ()
void	beginJob ()
void	beginRun (const edm::Run &run, const edm::EventSetup &context)
	Begin Run. More...
void	bookHistos (const DTChamberId &dtCh)
	Book the histograms. More...
void	endJob (void)
	EndJob. More...
Protected Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
CurrentProcessingContext const *	currentContext () const
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
float	angleRange
std::string	baseDirectory
float	bxTime
DQMStore *	dbe
MonitorElement *	dcc_IDDataErrorPlot
int	fineDelay
float	minHitsPhi
edm::ESHandle< DTGeometry >	muonGeom
int	nBXHigh
int	nBXLow
int	nevents
edm::ParameterSet	parameters
int	phCodeBestDCC [6][5][13]
int	phCodeBestDDU [6][5][13]
int	phCodeBXDCC [6][5][13][3]
bool	rangeInBX
int	segHitBest [6][5][13]
int	thCodeBestDDU [6][5][13]
std::map< uint32_t, std::map < std::string, MonitorElement * > >	triggerHistos
DTTTrigBaseSync *	tTrigSync

Friends
class	DTMonitorModule

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

Definition at line 44 of file DTLocalTriggerSynchTask.h.

Constructor & Destructor Documentation

DTLocalTriggerSynchTask::DTLocalTriggerSynchTask

(

const edm::ParameterSet &

ps

)

Constructor.

Definition at line 50 of file DTLocalTriggerSynchTask.cc.

References parameters.

                                                                          : nevents(0) {
  
  edm::LogVerbatim ("DTLocalTriggerSynchTask")  << "[DTLocalTriggerSynchTask]: Constructor" << endl;
  parameters = ps;

}

DTLocalTriggerSynchTask::~DTLocalTriggerSynchTask ( )

virtual

Destructor.

Definition at line 58 of file DTLocalTriggerSynchTask.cc.

References nevents.

                                                  {

  edm::LogVerbatim ("DTLocalTriggerSynchTask")  << "[DTLocalTriggerSynchTask]: analyzed " << nevents << " events" << endl;

}

Member Function Documentation

void DTLocalTriggerSynchTask::analyze ( const edm::Event &  event, const edm::EventSetup &  context  )

protectedvirtual

Analyze.

Implements edm::EDAnalyzer.

Definition at line 113 of file DTLocalTriggerSynchTask.cc.

References angleRange, bxTime, DTRecSegment2D::degreesOfFreedom(), dir, edm::ParameterSet::getParameter(), i, DTRecSegment2D::ist0Valid(), j, gen::k, minHitsPhi, nevents, DTTTrigBaseSync::offset(), parameters, phCodeBestDCC, phCodeBestDDU, DTRecSegment4D::phiSegment(), Geom::pi(), rangeInBX, DetId::rawId(), DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, cond::rpcobgas::time, triggerHistos, tTrigSync, and DTChamberId::wheel().

                                                                                        {

  nevents++;

  InputTag inputTagDCC  = parameters.getParameter<edm::InputTag>("DCCInputTag");
  InputTag inputTagDDU  = parameters.getParameter<edm::InputTag>("DDUInputTag");
  InputTag inputTagSEG  = parameters.getParameter<edm::InputTag>("SEGInputTag");

  for (int i=0;i<5;++i){
    for (int j=0;j<6;++j){
      for (int k=0;k<13;++k){
    phCodeBestDCC[j][i][k] = -1;
    phCodeBestDDU[j][i][k] = -1;
      }
    }
  }

  // Get best DCC triggers
  edm::Handle<L1MuDTChambPhContainer> l1DTTPGPh;
  event.getByLabel(inputTagDCC,l1DTTPGPh);
  vector<L1MuDTChambPhDigi>*  phTrigs = l1DTTPGPh->getContainer();
  
  vector<L1MuDTChambPhDigi>::const_iterator iph  = phTrigs->begin();
  vector<L1MuDTChambPhDigi>::const_iterator iphe = phTrigs->end();
  for(; iph !=iphe ; ++iph) {
    
    int phwheel = iph->whNum();
    int phsec   = iph->scNum() + 1; // DTTF[0-11] -> DT[1-12] Sector Numbering
    int phst    = iph->stNum();
    int phcode  = iph->code();

    if(phcode>phCodeBestDCC[phwheel+3][phst][phsec] && phcode<7) {
      phCodeBestDCC[phwheel+3][phst][phsec]=phcode; 
    }
    
  }

  // Get best DDU triggers
  Handle<DTLocalTriggerCollection> trigsDDU;
  event.getByLabel(inputTagDDU,trigsDDU);
  DTLocalTriggerCollection::DigiRangeIterator detUnitIt;

  for (detUnitIt=trigsDDU->begin();detUnitIt!=trigsDDU->end();++detUnitIt){
      
    const DTChamberId& id = (*detUnitIt).first;
    const DTLocalTriggerCollection::Range& range = (*detUnitIt).second;

    int wh = id.wheel();
    int sec = id.sector();
    int st = id.station();

    for (DTLocalTriggerCollection::const_iterator trigIt = range.first; trigIt!=range.second;++trigIt){
    
      int quality = trigIt->quality();
      
      if(quality>-1 && quality<7 &&
     quality>phCodeBestDDU[wh+3][st][sec]) {
    phCodeBestDDU[wh+3][st][sec]=quality;
      }
    }
  }
  
  //Get best segments (highest number of phi hits)
  vector<const DTRecSegment4D*> bestSegments4D;
  Handle<DTRecSegment4DCollection> segments4D;
  event.getByLabel(inputTagSEG, segments4D);
  DTRecSegment4DCollection::const_iterator track;
  DTRecSegment4DCollection::id_iterator chambIdIt;

  for (chambIdIt = segments4D->id_begin(); chambIdIt != segments4D->id_end(); ++chambIdIt){
    
    DTRecSegment4DCollection::range  range = segments4D->get(*chambIdIt);
    const DTRecSegment4D* best=0;
    int hitsBest = 0;
    int hits = 0;

    for ( track = range.first; track != range.second; ++track){
      if( (*track).hasPhi() ) {
    hits = (*track).phiSegment()->degreesOfFreedom()+2;
    if ( hits>hitsBest ){
      best = &(*track);
      hitsBest = hits;
    }
      }
    }
    if (best) {
      bestSegments4D.push_back(best);
    }
  }
    
  
  // Filling histos
  vector<const DTRecSegment4D*>::const_iterator bestSegIt  = bestSegments4D.begin();
  vector<const DTRecSegment4D*>::const_iterator bestSegEnd = bestSegments4D.end();
  for (; bestSegIt!=bestSegEnd; ++bestSegIt ){

    float dir = atan((*bestSegIt)->localDirection().x()/ (*bestSegIt)->localDirection().z())*180/Geom::pi(); // CB cerca un modo migliore x farlo
    const DTRecSegment2D* seg2D = (*bestSegIt)->phiSegment();
    int nHitsPhi = seg2D->degreesOfFreedom()+2; 
    DTChamberId chambId = (*bestSegIt)->chamberId();
    map<string, MonitorElement*> &innerME = triggerHistos[chambId.rawId()];
    
    if (fabs(dir)<angleRange && 
    nHitsPhi>=minHitsPhi && 
    seg2D->ist0Valid()){
      
      float t0seg = (*bestSegIt)->phiSegment()->t0();
      float tTrig = (tTrigSync->offset(DTWireId(chambId,1,1,2)) + tTrigSync->offset(DTWireId(chambId,3,1,2)) )/2;
      float time = tTrig+t0seg;
      float htime = rangeInBX ? time-int(time/bxTime)*bxTime : time-int(tTrig/bxTime)*bxTime;
      
      int wheel   = chambId.wheel();
      int sector  = chambId.sector();
      int station = chambId.station();
      int scsector = sector>12 ? sector==13 ? 4 : 10 : sector;

      int qualDCC = phCodeBestDCC[wheel+3][station][scsector];
      int qualDDU = phCodeBestDDU[wheel+3][station][scsector];

      if (fabs(t0seg)>0.01) {
    innerME.find("SEG_TrackCrossingTime")->second->Fill(htime);
    if ( qualDCC>=0 ) innerME.find("DCC_TrackCrossingTimeAll")->second->Fill(htime);      
    if ( qualDCC==6 ) innerME.find("DCC_TrackCrossingTimeHH")->second->Fill(htime);
    if ( qualDDU>=0 ) innerME.find("DDU_TrackCrossingTimeAll")->second->Fill(htime);      
    if ( qualDDU==6 ) innerME.find("DDU_TrackCrossingTimeHH")->second->Fill(htime);
      }

    }
  }

}

std::string& DTLocalTriggerSynchTask::baseDir ( )

inlineprotected

Definition at line 73 of file DTLocalTriggerSynchTask.h.

References baseDirectory.

Referenced by beginJob(), bookHistos(), and endJob().

{ return baseDirectory; }

void DTLocalTriggerSynchTask::beginJob ( void  )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 65 of file DTLocalTriggerSynchTask.cc.

References angleRange, baseDir(), baseDirectory, bxTime, dbe, edm::ParameterSet::getParameter(), minHitsPhi, nBXHigh, nBXLow, cppFunctionSkipper::operator, parameters, and rangeInBX.

                                      {

  edm::LogVerbatim ("DTLocalTriggerSynchTask") << "[DTLocalTriggerSynchTask]: BeginJob" << endl;

  bxTime        = parameters.getParameter<double>("bxTimeInterval");   // CB move this to static const or DB
  rangeInBX     = parameters.getParameter<bool>("rangeWithinBX");
  nBXLow        = parameters.getParameter<int>("nBXLow");
  nBXHigh       = parameters.getParameter<int>("nBXHigh");
  angleRange    = parameters.getParameter<double>("angleRange");
  minHitsPhi    = parameters.getParameter<int>("minHitsPhi");
  baseDirectory = parameters.getParameter<string>("baseDir");

  dbe = edm::Service<DQMStore>().operator->();
  dbe->setCurrentFolder(baseDir());
  dbe->bookFloat("BXTimeSpacing")->Fill(bxTime);
  
}

void DTLocalTriggerSynchTask::beginRun ( const edm::Run &  run, const edm::EventSetup &  context  )

protectedvirtual

Begin Run.

Reimplemented from edm::EDAnalyzer.

Definition at line 83 of file DTLocalTriggerSynchTask.cc.

References bookHistos(), edm::EventSetup::get(), reco::get(), edm::ParameterSet::getParameter(), muonGeom, DTTTrigBaseSync::offset(), parameters, DTTTrigBaseSync::setES(), AlCaHLTBitMon_QueryRunRegistry::string, triggerHistos, and tTrigSync.

                                                                                {

  edm::LogVerbatim ("DTLocalTriggerSynchTask") <<"[DTLocalTriggerSynchTask]: Begin Run"<<endl;

  context.get<MuonGeometryRecord>().get(muonGeom);
  tTrigSync = DTTTrigSyncFactory::get()->create(parameters.getParameter<std::string>("tTrigMode"),
                        parameters.getParameter<edm::ParameterSet>("tTrigModeConfig"));
  tTrigSync->setES(context);


  std::vector<DTChamber*>::const_iterator chambIt  = muonGeom->chambers().begin();
  std::vector<DTChamber*>::const_iterator chambEnd = muonGeom->chambers().end();

  for (; chambIt!=chambEnd; ++chambIt) { 
    bookHistos((*chambIt)->id());
    triggerHistos[(*chambIt)->id().rawId()]["tTrig_SL1"]->Fill(tTrigSync->offset(DTWireId((*chambIt)->id(),1,1,2)));
    triggerHistos[(*chambIt)->id().rawId()]["tTrig_SL3"]->Fill(tTrigSync->offset(DTWireId((*chambIt)->id(),3,1,2)));
  } 
  
}

void DTLocalTriggerSynchTask::bookHistos ( const DTChamberId &  dtCh )

protected

Book the histograms.

Definition at line 245 of file DTLocalTriggerSynchTask.cc.

References baseDir(), DQMStore::book1D(), DQMStore::bookFloat(), bxTime, dbe, max(), min, pileupCalc::nbins, nBXHigh, nBXLow, rangeInBX, DetId::rawId(), DTChamberId::sector(), DQMStore::setCurrentFolder(), DTChamberId::station(), relativeConstraints::station, triggerHistos, and DTChamberId::wheel().

Referenced by beginRun().

                                                                  {
  
  stringstream wheel; wheel << dtChId.wheel();  
  stringstream station; station << dtChId.station();    
  stringstream sector; sector << dtChId.sector();
  uint32_t chRawId = dtChId.rawId();
  
  dbe->setCurrentFolder(baseDir() + "/Wheel" + wheel.str() + "/Sector" + sector.str() + "/Station" + station.str() );
  

  string histoTag[5] = { "SEG_TrackCrossingTime", "DCC_TrackCrossingTimeAll", "DCC_TrackCrossingTimeHH", "DDU_TrackCrossingTimeAll", "DDU_TrackCrossingTimeHH" };
 
  float min = rangeInBX ?      0 : nBXLow*bxTime;
  float max = rangeInBX ? bxTime : nBXHigh*bxTime;
  int nbins = static_cast<int>(ceil( rangeInBX ? bxTime : (nBXHigh-nBXLow)*bxTime));

  for (int iHisto=0;iHisto<5;++iHisto) {
    string histoName = histoTag[iHisto] + (rangeInBX ? "InBX" : "") + "_W" + wheel.str() + "_Sec" + sector.str() + "_St" + station.str();
    edm::LogVerbatim ("DTLocalTriggerSynchTask") << "[DTLocalTriggerSynchTask]: booking " 
                         << baseDir() + "/Wheel" << wheel.str()
                        << "/Sector" << sector.str()
                        << "/Station"<< station.str() 
                        << "/" << histoName << endl;

    triggerHistos[chRawId][histoTag[iHisto]] = dbe->book1D(histoName.c_str(),"Track time distribution",nbins,min,max);
  }

  string floatTag[2] = { "tTrig_SL1", "tTrig_SL3" };

  for (int iFloat=0;iFloat<2;++iFloat) { 
    string floatName = floatTag[iFloat] + "_W" + wheel.str() + "_Sec" + sector.str() + "_St" + station.str();
    triggerHistos[chRawId][floatTag[iFloat]] = dbe->bookFloat(floatName);
  }

}

void DTLocalTriggerSynchTask::endJob ( void  )

protectedvirtual

EndJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 105 of file DTLocalTriggerSynchTask.cc.

References baseDir(), dbe, nevents, and DQMStore::rmdir().

                                     {

  edm::LogVerbatim ("DTLocalTriggerSynchTask")  << "[DTLocalTriggerSynchTask]: analyzed " << nevents << " events" << endl;
  dbe->rmdir(baseDir());

}

Friends And Related Function Documentation

friend class DTMonitorModule

friend

Definition at line 46 of file DTLocalTriggerSynchTask.h.

Member Data Documentation

float DTLocalTriggerSynchTask::angleRange

private

Definition at line 89 of file DTLocalTriggerSynchTask.h.

Referenced by analyze(), and beginJob().

std::string DTLocalTriggerSynchTask::baseDirectory

private

Definition at line 94 of file DTLocalTriggerSynchTask.h.

Referenced by baseDir(), and beginJob().

float DTLocalTriggerSynchTask::bxTime

private

Definition at line 85 of file DTLocalTriggerSynchTask.h.

Referenced by analyze(), beginJob(), and bookHistos().

DQMStore* DTLocalTriggerSynchTask::dbe

private

Definition at line 96 of file DTLocalTriggerSynchTask.h.

Referenced by beginJob(), bookHistos(), and endJob().

MonitorElement* DTLocalTriggerSynchTask::dcc_IDDataErrorPlot

private

Definition at line 100 of file DTLocalTriggerSynchTask.h.

int DTLocalTriggerSynchTask::fineDelay

private

Definition at line 91 of file DTLocalTriggerSynchTask.h.

float DTLocalTriggerSynchTask::minHitsPhi

private

Definition at line 90 of file DTLocalTriggerSynchTask.h.

Referenced by analyze(), and beginJob().

edm::ESHandle<DTGeometry> DTLocalTriggerSynchTask::muonGeom

private

Definition at line 98 of file DTLocalTriggerSynchTask.h.

Referenced by beginRun().

int DTLocalTriggerSynchTask::nBXHigh

private

Definition at line 88 of file DTLocalTriggerSynchTask.h.

Referenced by beginJob(), and bookHistos().

int DTLocalTriggerSynchTask::nBXLow

private

Definition at line 87 of file DTLocalTriggerSynchTask.h.

Referenced by beginJob(), and bookHistos().

int DTLocalTriggerSynchTask::nevents

private

Definition at line 77 of file DTLocalTriggerSynchTask.h.

Referenced by analyze(), endJob(), and ~DTLocalTriggerSynchTask().

edm::ParameterSet DTLocalTriggerSynchTask::parameters

private

Definition at line 97 of file DTLocalTriggerSynchTask.h.

Referenced by analyze(), beginJob(), beginRun(), DTLocalTriggerSynchTask(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputTags(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::recursePSetProperties().

int DTLocalTriggerSynchTask::phCodeBestDCC[6][5][13]

private

Definition at line 79 of file DTLocalTriggerSynchTask.h.

Referenced by analyze().

int DTLocalTriggerSynchTask::phCodeBestDDU[6][5][13]

private

Definition at line 81 of file DTLocalTriggerSynchTask.h.

Referenced by analyze().

int DTLocalTriggerSynchTask::phCodeBXDCC[6][5][13][3]

private

Definition at line 80 of file DTLocalTriggerSynchTask.h.

bool DTLocalTriggerSynchTask::rangeInBX

private

Definition at line 86 of file DTLocalTriggerSynchTask.h.

Referenced by analyze(), beginJob(), and bookHistos().

int DTLocalTriggerSynchTask::segHitBest[6][5][13]

private

Definition at line 83 of file DTLocalTriggerSynchTask.h.

int DTLocalTriggerSynchTask::thCodeBestDDU[6][5][13]

private

Definition at line 82 of file DTLocalTriggerSynchTask.h.

std::map<uint32_t, std::map<std::string, MonitorElement*> > DTLocalTriggerSynchTask::triggerHistos

private

Definition at line 99 of file DTLocalTriggerSynchTask.h.

Referenced by analyze(), beginRun(), and bookHistos().

DTTTrigBaseSync* DTLocalTriggerSynchTask::tTrigSync

private

Definition at line 92 of file DTLocalTriggerSynchTask.h.

Referenced by analyze(), and beginRun().