L1TDTTPG Class Reference

#include <L1TDTTPG.h>

Inheritance diagram for L1TDTTPG:

Public Member Functions
	L1TDTTPG (const edm::ParameterSet &ps)
	~L1TDTTPG () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final
virtual void	analyze (edm::Event const &, edm::EventSetup const &)
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
	DQMEDAnalyzer ()
	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete
void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override
void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) override
void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override
	~DQMEDAnalyzer () override=default
Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () const final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDProducerBase () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override
void	bookHistograms (DQMStore::IBooker &ibooker, edm::Run const &, edm::EventSetup const &) override
void	dqmBeginRun (edm::Run const &, edm::EventSetup const &) override
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 Member Functions
void	setMapPhLabel (MonitorElement *me)
void	setMapThLabel (MonitorElement *me)

Private Attributes
MonitorElement *	dttf_p_phi [3]
MonitorElement *	dttf_p_pt [3]
MonitorElement *	dttf_p_q [3]
MonitorElement *	dttf_p_qual [3]
MonitorElement *	dttpgphbestmap
MonitorElement *	dttpgphbestmapcorr
MonitorElement *	dttpgphbx [8]
MonitorElement *	dttpgphbxcomp
MonitorElement *	dttpgphmap
MonitorElement *	dttpgphmap2nd
MonitorElement *	dttpgphmapbx [3]
MonitorElement *	dttpgphmapcorr
MonitorElement *	dttpgphntrack
MonitorElement *	dttpgphquality [3]
MonitorElement *	dttpgphsector [3]
MonitorElement *	dttpgphstation [3]
MonitorElement *	dttpgphts2tag [3]
MonitorElement *	dttpgphwheel [3]
edm::InputTag	dttpgSource_
edm::EDGetTokenT< L1MuDTChambPhContainer >	dttpgSourcePhContainer_token_
edm::EDGetTokenT< L1MuDTChambThContainer >	dttpgSourceThContainer_token_
MonitorElement *	dttpgthbestmap
MonitorElement *	dttpgthbestmaph
MonitorElement *	dttpgthbx [3]
MonitorElement *	dttpgthmap
MonitorElement *	dttpgthmapbx [3]
MonitorElement *	dttpgthmaph
MonitorElement *	dttpgthntrack
MonitorElement *	dttpgthquality [3]
MonitorElement *	dttpgthsector [3]
MonitorElement *	dttpgthstation [3]
MonitorElement *	dttpgththeta [3]
MonitorElement *	dttpgthwheel [3]
std::ofstream	logFile_
bool	monitorDaemon_
int	nev_
std::string	outputFile_
std::string	trstring_
edm::EDGetTokenT< L1MuDTTrackContainer >	trToken_
bool	verbose_

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 40 of file L1TDTTPG.h.

Constructor & Destructor Documentation

L1TDTTPG::L1TDTTPG

(

const edm::ParameterSet &

ps

)

Definition at line 93 of file L1TDTTPG.cc.

References gather_cfg::cout, dttpgSource_, edm::ParameterSet::getUntrackedParameter(), edm::InputTag::label(), outputFile_, edm::InputTag::process(), trstring_, trToken_, and verbose_.

  : dttpgSourcePhContainer_token_( consumes<L1MuDTChambPhContainer>(ps.getParameter< InputTag >("dttpgSource") )),
    dttpgSourceThContainer_token_( consumes<L1MuDTChambThContainer>(ps.getParameter< InputTag >("dttpgSource") )),
    dttpgSource_( ps.getParameter< InputTag >("dttpgSource") )
{
  trstring_ = dttpgSource_.label()+":"+"DATA"+":"+dttpgSource_.process();
  trToken_ = consumes<L1MuDTTrackContainer>(trstring_);

  // verbosity switch
  verbose_ = ps.getUntrackedParameter<bool>("verbose", false);

  if(verbose_) cout << "L1TDTTPG: constructor...." << endl;

  outputFile_ = ps.getUntrackedParameter<string>("outputFile", "");
  if ( !outputFile_.empty() ) {
    cout << "L1T Monitoring histograms will be saved to " << outputFile_.c_str() << endl;
  }

  bool disable = ps.getUntrackedParameter<bool>("disableROOToutput", false);
  if(disable){
    outputFile_="";
  }
}

L1TDTTPG::~L1TDTTPG ( )

override

Definition at line 117 of file L1TDTTPG.cc.

{
}

Member Function Documentation

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

overrideprotected

Definition at line 288 of file L1TDTTPG.cc.

References L1MuDTChambPhDigi::code(), gather_cfg::cout, dttf_p_phi, dttf_p_pt, dttf_p_q, dttf_p_qual, dttpgphbestmap, dttpgphbestmapcorr, dttpgphbx, dttpgphbxcomp, dttpgphmap, dttpgphmap2nd, dttpgphmapbx, dttpgphmapcorr, dttpgphntrack, dttpgphquality, dttpgphsector, dttpgphstation, dttpgphts2tag, dttpgphwheel, dttpgSource_, dttpgSourcePhContainer_token_, dttpgSourceThContainer_token_, dttpgthbestmap, dttpgthbestmaph, dttpgthbx, dttpgthmap, dttpgthmapbx, dttpgthmaph, dttpgthntrack, dttpgthquality, dttpgthsector, dttpgthstation, dttpgththeta, dttpgthwheel, MonitorElement::Fill(), edm::Event::getByToken(), L1MuDTChambPhContainer::getContainer(), L1MuDTChambThContainer::getContainer(), L1MuDTTrackContainer::getContainer(), mps_fire::i, edm::HandleBase::isValid(), gen::k, edm::InputTag::label(), nev_, lumiQTWidget::t, trToken_, and verbose_.

{

  nev_++; 
  if(verbose_) cout << "L1TDTTPG: analyze...." << endl;

  edm::Handle<L1MuDTChambPhContainer > myL1MuDTChambPhContainer;  
  e.getByToken(dttpgSourcePhContainer_token_,myL1MuDTChambPhContainer);
  
  if (!myL1MuDTChambPhContainer.isValid()) {
    edm::LogInfo("DataNotFound") << "can't find L1MuDTChambPhContainer with label "
                 << dttpgSource_.label() ;
    return;
  }
  L1MuDTChambPhContainer::Phi_Container const *myPhContainer =  
    myL1MuDTChambPhContainer->getContainer();

  edm::Handle<L1MuDTChambThContainer > myL1MuDTChambThContainer;  
  e.getByToken(dttpgSourceThContainer_token_,myL1MuDTChambThContainer);
  
  if (!myL1MuDTChambThContainer.isValid()) {
    edm::LogInfo("DataNotFound") << "can't find L1MuDTChambThContainer with label "
                 << dttpgSource_.label() ;
    edm::LogInfo("DataNotFound") << "if this fails try to add DATA to the process name." ;

    return;
  }
  L1MuDTChambThContainer::The_Container const* myThContainer =  
    myL1MuDTChambThContainer->getContainer();

  int ndttpgphtrack = 0;
  int ndttpgthtrack = 0; 
  int NumberOfSegmentsPhi[3]={0,0,0};
  
  for( L1MuDTChambPhContainer::Phi_Container::const_iterator 
     DTPhDigiItr =  myPhContainer->begin() ;
       DTPhDigiItr != myPhContainer->end() ;
       ++DTPhDigiItr ) 
    {           
      int bx = DTPhDigiItr->bxNum() - DTPhDigiItr->Ts2Tag();
      if(bx == -1)
       NumberOfSegmentsPhi[0]++;
      if(bx == 0)
       NumberOfSegmentsPhi[1]++;
      if(bx == 1)
       NumberOfSegmentsPhi[2]++;   
    }
   /*Fill Histos for Segment counter for each bx separately */

  for(int k=0;k<3;k++){
     dttpgphbx[k+2]->Fill(NumberOfSegmentsPhi[k]);
   }
   int bxCounterDttfPhi=0; // = no. of BX's with non-zero data
   for (int k=0;k<3;k++){
     if (NumberOfSegmentsPhi[k]>0)
       bxCounterDttfPhi++;
   }

   /* the BX "code" */

   int bxCodePhi=0;
   if(bxCounterDttfPhi==0){
     bxCodePhi=0;
   }else if(bxCounterDttfPhi==1){
     for(int k=0;k<3;k++){
       if(NumberOfSegmentsPhi[k]>0)
        bxCodePhi=k+2;
     }
   }else if(bxCounterDttfPhi==2){
     for(int k=0;k<3;k++){
       if(NumberOfSegmentsPhi[k]==0)
        bxCodePhi=8-k;
     }
   }else if(bxCounterDttfPhi==3){
     bxCodePhi=10;
   }

   //The bx analyzer histo
   dttpgphbx[0]->Fill(bxCodePhi);


   const L1MuDTChambPhDigi* bestPhQualMap[5][12][4];
   memset(bestPhQualMap,0,240*sizeof(L1MuDTChambPhDigi*));

   for( L1MuDTChambPhContainer::Phi_Container::const_iterator 
     DTPhDigiItr =  myPhContainer->begin() ;
       DTPhDigiItr != myPhContainer->end() ;
       ++DTPhDigiItr ) 
    {           

      ndttpgphtrack++;

      int bxindex = DTPhDigiItr->bxNum() - DTPhDigiItr->Ts2Tag() + 1;

      dttpgphwheel[bxindex]->Fill(DTPhDigiItr->whNum());
      if (verbose_)
    {
      cout << "DTTPG phi wheel number " << DTPhDigiItr->whNum() << endl;
    }
      dttpgphstation[bxindex]->Fill(DTPhDigiItr->stNum());
      if (verbose_)
    {   
      cout << "DTTPG phi station number " << DTPhDigiItr->stNum() << endl;
    }
      dttpgphsector[bxindex]->Fill(DTPhDigiItr->scNum());
      if (verbose_)
    {
      cout << "DTTPG phi sector number " << DTPhDigiItr->scNum() << endl;
    }
      dttpgphquality[bxindex]->Fill(DTPhDigiItr->code());
      if (verbose_)
    {
      cout << "DTTPG phi quality " << DTPhDigiItr->code() << endl;
    }
      dttpgphts2tag[bxindex]->Fill(DTPhDigiItr->Ts2Tag());
      if (verbose_)
    {
      cout << "DTTPG phi ts2tag " << DTPhDigiItr->Ts2Tag() << endl;
    }
      int ypos = DTPhDigiItr->scNum();
      int xpos = DTPhDigiItr->stNum()+4*(DTPhDigiItr->whNum()+2);
      dttpgphmap->Fill(xpos,ypos);
      if (DTPhDigiItr->Ts2Tag())
    dttpgphmap2nd->Fill(xpos,ypos);
      dttpgphmapbx[bxindex]->Fill(xpos,ypos);
      if (DTPhDigiItr->code()>3)
    dttpgphmapcorr->Fill(xpos,ypos);

      if (bestPhQualMap[DTPhDigiItr->whNum()+2][ DTPhDigiItr->scNum()][DTPhDigiItr->stNum()-1]==nullptr ||
      bestPhQualMap[DTPhDigiItr->whNum()+2][ DTPhDigiItr->scNum()][DTPhDigiItr->stNum()-1]->code()<DTPhDigiItr->code())
    {
      bestPhQualMap[DTPhDigiItr->whNum()+2][ DTPhDigiItr->scNum()][DTPhDigiItr->stNum()-1]=&(*DTPhDigiItr);
    }

    }

   for (int iwh=0; iwh<5; iwh++){
     for (int isec=0; isec<12; isec++){
       for (int ist=0; ist<4; ist++){
     if (bestPhQualMap[iwh][isec][ist]){
       int xpos = iwh*4+ist+1;
       dttpgphbestmap->Fill(xpos,isec);
       if(bestPhQualMap[iwh][isec][ist]->code()>3)
         dttpgphbestmapcorr->Fill(xpos,isec);
     }
       }
     }
   }


   int bestThQualMap[5][12][3];
   memset(bestThQualMap,0,180*sizeof(int));
   //for( vector<L1MuDTChambThDigi>::const_iterator 
   for( L1MuDTChambThContainer::The_Container::const_iterator 
     DTThDigiItr =  myThContainer->begin() ;
       DTThDigiItr != myThContainer->end() ;
       ++DTThDigiItr ) 
    {                   
      ndttpgthtrack++;

      int bxindex = DTThDigiItr->bxNum() + 1;

      dttpgthwheel[bxindex]->Fill(DTThDigiItr->whNum());
      if (verbose_)
    {
      cout << "DTTPG theta wheel number " << DTThDigiItr->whNum() << endl;
    }
      dttpgthstation[bxindex]->Fill(DTThDigiItr->stNum());
      if (verbose_)
    {   
      cout << "DTTPG theta station number " << DTThDigiItr->stNum() << endl;
    }
      dttpgthsector[bxindex]->Fill(DTThDigiItr->scNum());
      if (verbose_)
    {
      cout << "DTTPG theta sector number " << DTThDigiItr->scNum() << endl;
    }
      dttpgthbx[bxindex]->Fill(DTThDigiItr->bxNum());
      if (verbose_)
    {
      cout << "DTTPG theta bx number " << DTThDigiItr->bxNum() << endl;
    }
      int thcode[7]= {0,0,0,0,0,0,0};
      for (int j = 0; j < 7; j++)
    {
      dttpgththeta[bxindex]->Fill(DTThDigiItr->position(j));
      if (verbose_)
        {
          cout << "DTTPG theta position " << DTThDigiItr->position(j) << endl;
        }
      thcode[j]=DTThDigiItr->code(j);
      dttpgthquality[bxindex]->Fill(thcode[j]);
      if (verbose_)
        {
          cout << "DTTPG theta quality " << DTThDigiItr->code(j) << endl;
        }
    }
      
      int ypos = DTThDigiItr->scNum();
      int xpos = DTThDigiItr->stNum()+4*(DTThDigiItr->whNum()+2);
      int bestqual=0;
      dttpgthmap->Fill(xpos,ypos);
      dttpgthmapbx[bxindex]->Fill(xpos,ypos);
      for (int pos = 0; pos < 7; pos++){
    if (thcode[pos]>bestqual)
      bestqual=thcode[pos];
    if(thcode[pos]==2)
      dttpgthmaph->Fill(xpos,ypos);
      }

      if (bestThQualMap[DTThDigiItr->whNum()+2][ DTThDigiItr->scNum()][DTThDigiItr->stNum()-1] < bestqual)
    {
      bestThQualMap[DTThDigiItr->whNum()+2][ DTThDigiItr->scNum()][DTThDigiItr->stNum()-1]=bestqual;
    }
    }

   for (int iwh=0; iwh<5; iwh++){
     for (int isec=0; isec<12; isec++){
       for (int ist=0; ist<3; ist++){
     if (bestThQualMap[iwh][isec][ist]){
       int xpos = iwh*4+ist+1;
       dttpgthbestmap->Fill(xpos,isec);
       if(bestThQualMap[iwh][isec][ist]==2)
         dttpgthbestmaph->Fill(xpos,isec);
     }
       }
     }
   }


  dttpgphntrack->Fill(ndttpgphtrack);
  if (verbose_)
    {
      cout << "DTTPG phi ntrack " << ndttpgphtrack << endl;
    }
  dttpgthntrack->Fill(ndttpgthtrack);
  if (verbose_) {
    cout << "DTTPG theta ntrack " << ndttpgthtrack << endl;
  }

  edm::Handle<L1MuDTTrackContainer > myL1MuDTTrackContainer;
  e.getByToken(trToken_,myL1MuDTTrackContainer);

  if (!myL1MuDTTrackContainer.isValid()) {
    edm::LogInfo("DataNotFound") << "can't find L1MuDTTrackContainer with label "
                               << dttpgSource_.label() ;
    return;
  }

  L1MuDTTrackContainer::TrackContainer const * t =  myL1MuDTTrackContainer->getContainer();



  int NumberOfSegmentsOut[3]={0,0,0};
  for ( L1MuDTTrackContainer::TrackContainer::const_iterator i 
      = t->begin(); i != t->end(); ++i ) {
    if(i->bx() ==-1)
       NumberOfSegmentsOut[0]++;
    if(i->bx() ==0)
       NumberOfSegmentsOut[1]++;
    if(i->bx() ==1)
       NumberOfSegmentsOut[2]++;
  }


   /*Fill Histos for Segment counter*/
   for(int k=0;k<3;k++){
     dttpgphbx[k+5]->Fill(NumberOfSegmentsOut[k]);
   }

   /*Bunch assigments*/

   int bxCounterDttfOut=0;
   for (int k=0;k<3;k++){
   if (NumberOfSegmentsOut[k]>0)
     bxCounterDttfOut++;
   }

   int bxCodeOut=0;
   if(bxCounterDttfOut==0){
     bxCodeOut=0;
   }else if(bxCounterDttfOut==1){
     for(int k=0;k<3;k++){
       if(NumberOfSegmentsOut[k]>0)
         bxCodeOut=k+2;
     }
   }else if(bxCounterDttfOut==2){
     for(int k=0;k<3;k++){
       if(NumberOfSegmentsOut[k]==0)
         bxCodeOut=8-k;
     }
   }else if(bxCounterDttfOut==3){
     bxCodeOut=10;
   }

   //The bx analyzer histo
   dttpgphbx[1]->Fill(bxCodeOut);

   /*End Dttf Output Bunch analysis*/

   // the 2-DIM histo with phi.input vs. output
   dttpgphbxcomp->Fill(bxCodePhi,bxCodeOut);


  for ( L1MuDTTrackContainer::TrackContainer::const_iterator i 
      = t->begin(); i != t->end(); ++i ) {
    if ( verbose_ ) {
      std::cout << "bx = " << i->bx() 
        << std::endl;
      std::cout << "quality (packed) = " << i->quality_packed() 
        << std::endl;
      std::cout << "pt      (packed) = " << i->pt_packed() 
        << std::endl;
      std::cout << "phi     (packed) = " << i->phi_packed() 
        << std::endl;
      std::cout << "charge  (packed) = " << i->charge_packed() 
        << std::endl;
    }


    int bxindex = i->bx() + 1;
    dttf_p_phi[bxindex]->Fill(i->phi_packed());
    dttf_p_qual[bxindex]->Fill(i->quality_packed());
    dttf_p_pt[bxindex]->Fill(i->pt_packed());
    dttf_p_q[bxindex]->Fill(i->charge_packed());
  }
    
}

void L1TDTTPG::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  , edm::EventSetup const &    )

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 128 of file L1TDTTPG.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), dttf_p_phi, dttf_p_pt, dttf_p_q, dttf_p_qual, dttpgphbestmap, dttpgphbestmapcorr, dttpgphbx, dttpgphbxcomp, dttpgphmap, dttpgphmap2nd, dttpgphmapbx, dttpgphmapcorr, dttpgphntrack, dttpgphquality, dttpgphsector, dttpgphstation, dttpgphts2tag, dttpgphwheel, dttpgthbestmap, dttpgthbestmaph, dttpgthbx, dttpgthmap, dttpgthmapbx, dttpgthmaph, dttpgthntrack, dttpgthquality, dttpgthsector, dttpgthstation, dttpgththeta, dttpgthwheel, mps_fire::i, nev_, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), setMapPhLabel(), and setMapThLabel().

{
  nev_ = 0;
  ibooker.setCurrentFolder("L1T/L1TDTTPG");


  //hist1[0]
  dttpgphbx[0] = ibooker.book1D("BxEncoding_PHI",
                 "Bunch encoding DTTF Phi",11,0,11);
  //hist1[1]
  dttpgphbx[1] = ibooker.book1D("BxEncoding_OUT",
                 "Bunch encoding DTTF Output",11,0,11);

      //hist1[2]
  dttpgphbx[2] = ibooker.book1D("NumberOfSegmentsPHI_BunchNeg1",
                 "Number of segments for bunch -1 Dttf Phi",
                 20,0,20);
      //hist1[3]
  dttpgphbx[3] = ibooker.book1D("NumberOfSegmentsPHI_Bunch0",
                 "Number of segments for bunch 0 Dttf Phi",
                 20,0,20);
      //hist1[4]
  dttpgphbx[4] = ibooker.book1D("NumberOfSegmentsPHI_Bunch1",
                 "Number of segments for bunch 1 Dttf Phi",
                 20,0,20);

      //hist1[5]
  dttpgphbx[5] = ibooker.book1D("NumberOfSegmentsOUT_BunchNeg1",
                 "Number of segments for bunch -1 Dttf Output",
                 20,0,20);
      //hist1[6] 
  dttpgphbx[6] = ibooker.book1D("NumberOfSegmentsOUT_Bunch0",
                 "Number of segments for bunch 0 Dttf Output",
                 20,0,20);
      //hist1[7]
  dttpgphbx[7] = ibooker.book1D("NumberOfSegmentsOUT_Bunch1",
                 "Number of segments for bunch 1 Dttf Output",
                 20,0,20);

  for(int i=0;i<2;i++){
    dttpgphbx[i]->setBinLabel(1,"None");
    dttpgphbx[i]->setBinLabel(3,"Only bx=-1");
    dttpgphbx[i]->setBinLabel(4,"Only bx= 0");
    dttpgphbx[i]->setBinLabel(5,"Only bx=+1");
    dttpgphbx[i]->setBinLabel(7,"Bx=-1,0");
    dttpgphbx[i]->setBinLabel(8,"Bx=-1,1");
    dttpgphbx[i]->setBinLabel(9,"Bx= 0,1");
    dttpgphbx[i]->setBinLabel(11,"All bx");
      }
   
  dttpgphbxcomp = ibooker.book2D("BxEncoding_PHI_OUT",
                  "Bunch encoding: DTTF Phi vs. Output",
                  11,0,11,11,0,11);
  dttpgphbxcomp->setAxisTitle("DTTF (output)",1);
  dttpgphbxcomp->setAxisTitle("PHI-TF",2);
  for(int i=1;i<=2;i++){
    dttpgphbxcomp->setBinLabel(1,"None",i);
    dttpgphbxcomp->setBinLabel(3,"Only bx=-1",i);
    dttpgphbxcomp->setBinLabel(4,"Only bx= 0",i);
    dttpgphbxcomp->setBinLabel(5,"Only bx=+1",i);
    dttpgphbxcomp->setBinLabel(7,"Bx=-1,0",i);
    dttpgphbxcomp->setBinLabel(8,"Bx=-1,1",i);
    dttpgphbxcomp->setBinLabel(9,"Bx= 0,1",i);
    dttpgphbxcomp->setBinLabel(11,"All bx",i);
  }

  dttpgphntrack = ibooker.book1D("DT_TPG_phi_ntrack", 
                  "DT TPG phi ntrack", 20, -0.5, 19.5 ) ;  
  dttpgthntrack = ibooker.book1D("DT_TPG_theta_ntrack", 
                  "DT TPG theta ntrack", 20, -0.5, 19.5 ) ;  

  for (int ibx=0 ; ibx<=2; ibx++) {
    
    ostringstream bxnum;
    bxnum << ibx-1;
    string bxn;
    if (ibx<2)
      bxn = bxnum.str();
    else
      bxn = "+" + bxnum.str();
    
    // Phi
    dttpgphwheel[ibx] = ibooker.book1D("DT_TPG_phi_wheel_number_"+bxn, 
                        "DT TPG phi wheel number "+bxn, 5, -2.5, 2.5 ) ;  
    dttpgphsector[ibx] = ibooker.book1D("DT_TPG_phi_sector_number_"+bxn, 
                     "DT TPG phi sector number "+bxn, 12, -0.5, 11.5 );  
    dttpgphstation[ibx] = ibooker.book1D("DT_TPG_phi_station_number_"+bxn, 
                      "DT TPG phi station number "+bxn, 5, 0.5, 4.5 ) ;
//  dttpgphphi[ibx] = dbe->book1D("DT_TPG_phi_"+bxn, 
//                    "DT TPG phi "+bxn, 100, -2100., 2100. ) ;  
//  dttpgphphiB[ibx] = dbe->book1D("DT_TPG_phiB_"+bxn, 
//                     "DT TPG phiB "+bxn, 100, -550., 550. ) ;  
    dttpgphquality[ibx] = ibooker.book1D("DT_TPG_phi_quality_"+bxn, 
                      "DT TPG phi quality "+bxn, 8, -0.5, 7.5 ) ;  
    dttpgphts2tag[ibx] = ibooker.book1D("DT_TPG_phi_Ts2Tag_"+bxn, 
                     "DT TPG phi Ts2Tag "+bxn, 2, -0.5, 1.5 ) ;  
//  dttpgphbxcnt[ibx] = dbe->book1D("DT_TPG_phi_BxCnt_"+bxn, 
//                  "DT TPG phi BxCnt "+bxn, 10, -0.5, 9.5 ) ;  
    dttpgphmapbx[ibx] = ibooker.book2D("DT_TPG_phi_map_bx"+bxn,
                      "Map of triggers per station (BX="+bxn+")",20,1,21,12,0,12);
    setMapPhLabel(dttpgphmapbx[ibx]);

    //Theta
    dttpgthbx[ibx] = ibooker.book1D("DT_TPG_theta_bx_"+bxn, 
                     "DT TPG theta bx "+bxn, 50, -24.5, 24.5 ) ;  
    dttpgthwheel[ibx] = ibooker.book1D("DT_TPG_theta_wheel_number_"+bxn, 
                    "DT TPG theta wheel number "+bxn, 5, -2.5, 2.5 ) ;  
    dttpgthsector[ibx] = ibooker.book1D("DT_TPG_theta_sector_number_"+bxn, 
                     "DT TPG theta sector number "+bxn, 12, -0.5, 11.5 ) ;  
    dttpgthstation[ibx] = ibooker.book1D("DT_TPG_theta_station_number_"+bxn, 
                      "DT TPG theta station number "+bxn, 5, -0.5, 4.5 ) ;  
    dttpgththeta[ibx] = ibooker.book1D("DT_TPG_theta_"+bxn, 
                    "DT TPG theta "+bxn, 20, -0.5, 19.5 ) ;  
    dttpgthquality[ibx] = ibooker.book1D("DT_TPG_theta_quality_"+bxn, 
                      "DT TPG theta quality "+bxn, 8, -0.5, 7.5 ) ;  
    dttpgthmapbx[ibx] = ibooker.book2D("DT_TPG_theta_map_bx_"+bxn,
                    "Map of triggers per station (BX="+bxn+")",15,1,16,12,0,12);
    setMapThLabel(dttpgthmapbx[ibx]);

    // Phi output
    dttf_p_phi[ibx] = ibooker.book1D("dttf_p_phi_"+bxn, "dttf phi output #phi "+bxn, 256, 
                      -0.5, 255.5);
    dttf_p_qual[ibx] = ibooker.book1D("dttf_p_qual_"+bxn, "dttf phi output qual "+bxn, 8, -0.5, 7.5);
    dttf_p_q[ibx] = ibooker.book1D("dttf_p_q_"+bxn, "dttf phi output q "+bxn, 2, -0.5, 1.5);
    dttf_p_pt[ibx] = ibooker.book1D("dttf_p_pt_"+bxn, "dttf phi output p_{t} "+bxn, 32, -0.5, 31.5);
      
  }

  dttpgphmap = ibooker.book2D("DT_TPG_phi_map",
                   "Map of triggers per station",20,1,21,12,0,12);
  dttpgphmapcorr = ibooker.book2D("DT_TPG_phi_map_corr",
                   "Map of correlated triggers per station",20,1,21,12,0,12);
  dttpgphmap2nd = ibooker.book2D("DT_TPG_phi_map_2nd",
                  "Map of second tracks per station",20,1,21,12,0,12);
  dttpgphbestmap = ibooker.book2D("DT_TPG_phi_best_map",
                   "Map of best triggers per station",20,1,21,12,0,12);
  dttpgphbestmapcorr = ibooker.book2D("DT_TPG_phi_best_map_corr",
                       "Map of correlated best triggers per station",20,1,21,12,0,12);
  setMapPhLabel(dttpgphmap);
  setMapPhLabel(dttpgphmapcorr);
  setMapPhLabel(dttpgphmap2nd);
  setMapPhLabel(dttpgphbestmap);
  setMapPhLabel(dttpgphbestmapcorr);
      


  dttpgthmap = ibooker.book2D("DT_TPG_theta_map",
                   "Map of triggers per station",15,1,16,12,0,12);
  dttpgthmaph = ibooker.book2D("DT_TPG_theta_map_h",
                "Map of H quality triggers per station",15,1,16,12,0,12);
  dttpgthbestmap = ibooker.book2D("DT_TPG_theta_best_map",
                   "Map of besttriggers per station",15,1,16,12,0,12);
  dttpgthbestmaph = ibooker.book2D("DT_TPG_theta_best_map_h",
                    "Map of H quality best triggers per station",15,1,16,12,0,12);
  setMapThLabel(dttpgthmap);
  setMapThLabel(dttpgthmaph);
  setMapThLabel(dttpgthbestmap);
  setMapThLabel(dttpgthbestmaph);
}

void L1TDTTPG::dqmBeginRun ( edm::Run const &  , edm::EventSetup const &    )

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 122 of file L1TDTTPG.cc.

                                                                 {
  //empty

}

void L1TDTTPG::setMapPhLabel ( MonitorElement *  me )

private

Definition at line 617 of file L1TDTTPG.cc.

References mps_fire::i, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), and makeMuonMisalignmentScenario::wheel.

Referenced by bookHistograms().

{

  me->setAxisTitle("DTTF Sector",2);
      for(int i=0;i<5;i++){
    ostringstream wheel;
    wheel << i-2;
    me->setBinLabel(1+i*4,"Wheel "+ wheel.str(),1);
      }
  
}

void L1TDTTPG::setMapThLabel ( MonitorElement *  me )

private

Definition at line 629 of file L1TDTTPG.cc.

References mps_fire::i, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), and makeMuonMisalignmentScenario::wheel.

Referenced by bookHistograms().

{

  me->setAxisTitle("DTTF Sector",2);
      for(int i=0;i<5;i++){
    ostringstream wheel;
    wheel << i-2;
    me->setBinLabel(1+i*3,"Wheel "+ wheel.str(),1);
      }
  
}

Member Data Documentation

MonitorElement* L1TDTTPG::dttf_p_phi[3]

private

Definition at line 98 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttf_p_pt[3]

private

Definition at line 99 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttf_p_q[3]

private

Definition at line 100 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttf_p_qual[3]

private

Definition at line 101 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphbestmap

private

Definition at line 81 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphbestmapcorr

private

Definition at line 82 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphbx[8]

private

Definition at line 66 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphbxcomp

private

Definition at line 67 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphmap

private

Definition at line 77 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphmap2nd

private

Definition at line 79 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphmapbx[3]

private

Definition at line 78 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphmapcorr

private

Definition at line 80 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphntrack

private

Definition at line 76 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphquality[3]

private

Definition at line 73 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphsector[3]

private

Definition at line 69 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphstation[3]

private

Definition at line 70 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphts2tag[3]

private

Definition at line 74 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphwheel[3]

private

Definition at line 68 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

edm::InputTag L1TDTTPG::dttpgSource_

private

Definition at line 110 of file L1TDTTPG.h.

Referenced by analyze(), and L1TDTTPG().

edm::EDGetTokenT<L1MuDTChambPhContainer> L1TDTTPG::dttpgSourcePhContainer_token_

private

Definition at line 108 of file L1TDTTPG.h.

Referenced by analyze().

edm::EDGetTokenT<L1MuDTChambThContainer> L1TDTTPG::dttpgSourceThContainer_token_

private

Definition at line 109 of file L1TDTTPG.h.

Referenced by analyze().

MonitorElement* L1TDTTPG::dttpgthbestmap

private

Definition at line 95 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgthbestmaph

private

Definition at line 96 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgthbx[3]

private

Definition at line 85 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgthmap

private

Definition at line 92 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgthmapbx[3]

private

Definition at line 93 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgthmaph

private

Definition at line 94 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgthntrack

private

Definition at line 91 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgthquality[3]

private

Definition at line 90 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgthsector[3]

private

Definition at line 87 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgthstation[3]

private

Definition at line 88 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgththeta[3]

private

Definition at line 89 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgthwheel[3]

private

Definition at line 86 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

std::ofstream L1TDTTPG::logFile_

private

Definition at line 107 of file L1TDTTPG.h.

bool L1TDTTPG::monitorDaemon_

private

Definition at line 106 of file L1TDTTPG.h.

int L1TDTTPG::nev_

private

Definition at line 103 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

std::string L1TDTTPG::outputFile_

private

Definition at line 104 of file L1TDTTPG.h.

Referenced by L1TDTTPG().

std::string L1TDTTPG::trstring_

private

Definition at line 112 of file L1TDTTPG.h.

Referenced by L1TDTTPG().

edm::EDGetTokenT<L1MuDTTrackContainer> L1TDTTPG::trToken_

private

Definition at line 113 of file L1TDTTPG.h.

Referenced by analyze(), and L1TDTTPG().

bool L1TDTTPG::verbose_

private

Definition at line 105 of file L1TDTTPG.h.

Referenced by analyze(), and L1TDTTPG().