#include <L1TDTTPG.h>

Inheritance diagram for L1TDTTPG:

Public Member Functions
	L1TDTTPG (const edm::ParameterSet &ps)

virtual	~L1TDTTPG ()

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	analyze (const edm::Event &e, const edm::EventSetup &c)

void	beginJob (void)

void	beginRun (edm::Run const &iRun, edm::EventSetup const &iSetup)

void	endJob (void)

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
DQMStore *	dbe

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::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 &)

Detailed Description

Definition at line 37 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, dbe, dttpgSource_, edm::ParameterSet::getUntrackedParameter(), edm::InputTag::label(), NULL, cppFunctionSkipper::operator, 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;
 
 
   dbe = NULL;
   if ( ps.getUntrackedParameter<bool>("DQMStore", false) ) 
     {
       dbe = Service<DQMStore>().operator->();
       dbe->setVerbose(0);
     }
   
   outputFile_ = ps.getUntrackedParameter<string>("outputFile", "");
   if ( outputFile_.size() != 0 ) {
     cout << "L1T Monitoring histograms will be saved to " << outputFile_.c_str() << endl;
   }
 
   bool disable = ps.getUntrackedParameter<bool>("disableROOToutput", false);
   if(disable){
     outputFile_="";
   }
 
 
   if ( dbe !=NULL ) {
     dbe->setCurrentFolder("L1T/L1TDTTPG");
   }
 
 
 }

L1TDTTPG::~L1TDTTPG ( )

virtual

Definition at line 132 of file L1TDTTPG.cc.

133 {

134 }

Member Function Documentation

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

protectedvirtual

Implements edm::EDAnalyzer.

Definition at line 323 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(), i, edm::HandleBase::isValid(), j, roll_playback::k, edm::InputTag::label(), nev_, edmStreamStallGrapher::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;
     }
 //       dttpgphphi[bxindex]->Fill(DTPhDigiItr->phi());
 //       if (verbose_)
 //  {
 //    cout << "DTTPG phi phi " << DTPhDigiItr->phi() << endl;
 //  }
 //       dttpgphphiB[bxindex]->Fill(DTPhDigiItr->phiB());
 //       if (verbose_)
 //  {
 //    cout << "DTTPG phi phiB " << DTPhDigiItr->phiB() << 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;
     }
 //       dttpgphbxcnt[bxindex]->Fill(DTPhDigiItr->BxCnt());
 //       if (verbose_)
 //  {
 //    cout << "DTTPG phi bxcnt " << DTPhDigiItr->BxCnt() << 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]==0 ||
       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::beginJob ( void )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 136 of file L1TDTTPG.cc.

References nev_.

 {
   nev_ = 0;
 }

void L1TDTTPG::beginRun	(	edm::Run const &	iRun,
		edm::EventSetup const &	iSetup
	)

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 152 of file L1TDTTPG.cc.

References dbe, 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, i, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), setMapPhLabel(), and setMapThLabel().

 {
   if ( dbe ) {
     dbe->setCurrentFolder("L1T/L1TDTTPG");
     dbe->rmdir("L1T/L1TDTTPG");
   }
 
 
   if ( dbe ) 
     {
       dbe->setCurrentFolder("L1T/L1TDTTPG");
 
 
       //hist1[0]
       dttpgphbx[0] = dbe->book1D("BxEncoding_PHI",
                  "Bunch encoding DTTF Phi",11,0,11);
       //hist1[1]
       dttpgphbx[1] = dbe->book1D("BxEncoding_OUT",
                  "Bunch encoding DTTF Output",11,0,11);
 
       //hist1[2]
       dttpgphbx[2] = dbe->book1D("NumberOfSegmentsPHI_BunchNeg1",
                  "Number of segments for bunch -1 Dttf Phi",
                  20,0,20);
       //hist1[3]
       dttpgphbx[3] = dbe->book1D("NumberOfSegmentsPHI_Bunch0",
                  "Number of segments for bunch 0 Dttf Phi",
                  20,0,20);
       //hist1[4]
       dttpgphbx[4] = dbe->book1D("NumberOfSegmentsPHI_Bunch1",
                  "Number of segments for bunch 1 Dttf Phi",
                  20,0,20);
 
       //hist1[5]
       dttpgphbx[5] = dbe->book1D("NumberOfSegmentsOUT_BunchNeg1",
                  "Number of segments for bunch -1 Dttf Output",
                  20,0,20);
       //hist1[6] 
       dttpgphbx[6] = dbe->book1D("NumberOfSegmentsOUT_Bunch0",
                  "Number of segments for bunch 0 Dttf Output",
                  20,0,20);
       //hist1[7]
       dttpgphbx[7] = dbe->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 = dbe->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 = dbe->book1D("DT_TPG_phi_ntrack", 
                   "DT TPG phi ntrack", 20, -0.5, 19.5 ) ;  
       dttpgthntrack = dbe->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] = dbe->book1D("DT_TPG_phi_wheel_number_"+bxn, 
                         "DT TPG phi wheel number "+bxn, 5, -2.5, 2.5 ) ;  
     dttpgphsector[ibx] = dbe->book1D("DT_TPG_phi_sector_number_"+bxn, 
                      "DT TPG phi sector number "+bxn, 12, -0.5, 11.5 );  
     dttpgphstation[ibx] = dbe->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] = dbe->book1D("DT_TPG_phi_quality_"+bxn, 
                       "DT TPG phi quality "+bxn, 8, -0.5, 7.5 ) ;  
     dttpgphts2tag[ibx] = dbe->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] = dbe->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] = dbe->book1D("DT_TPG_theta_bx_"+bxn, 
                      "DT TPG theta bx "+bxn, 50, -24.5, 24.5 ) ;  
     dttpgthwheel[ibx] = dbe->book1D("DT_TPG_theta_wheel_number_"+bxn, 
                     "DT TPG theta wheel number "+bxn, 5, -2.5, 2.5 ) ;  
     dttpgthsector[ibx] = dbe->book1D("DT_TPG_theta_sector_number_"+bxn, 
                      "DT TPG theta sector number "+bxn, 12, -0.5, 11.5 ) ;  
     dttpgthstation[ibx] = dbe->book1D("DT_TPG_theta_station_number_"+bxn, 
                       "DT TPG theta station number "+bxn, 5, -0.5, 4.5 ) ;  
     dttpgththeta[ibx] = dbe->book1D("DT_TPG_theta_"+bxn, 
                     "DT TPG theta "+bxn, 20, -0.5, 19.5 ) ;  
     dttpgthquality[ibx] = dbe->book1D("DT_TPG_theta_quality_"+bxn, 
                       "DT TPG theta quality "+bxn, 8, -0.5, 7.5 ) ;  
     dttpgthmapbx[ibx] = dbe->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] = dbe->book1D("dttf_p_phi_"+bxn, "dttf phi output #phi "+bxn, 256, 
                       -0.5, 255.5);
     dttf_p_qual[ibx] = dbe->book1D("dttf_p_qual_"+bxn, "dttf phi output qual "+bxn, 8, -0.5, 7.5);
     dttf_p_q[ibx] = dbe->book1D("dttf_p_q_"+bxn, "dttf phi output q "+bxn, 2, -0.5, 1.5);
     dttf_p_pt[ibx] = dbe->book1D("dttf_p_pt_"+bxn, "dttf phi output p_{t} "+bxn, 32, -0.5, 31.5);
       
       }
 
       dttpgphmap = dbe->book2D("DT_TPG_phi_map",
                    "Map of triggers per station",20,1,21,12,0,12);
       dttpgphmapcorr = dbe->book2D("DT_TPG_phi_map_corr",
                    "Map of correlated triggers per station",20,1,21,12,0,12);
       dttpgphmap2nd = dbe->book2D("DT_TPG_phi_map_2nd",
                   "Map of second tracks per station",20,1,21,12,0,12);
       dttpgphbestmap = dbe->book2D("DT_TPG_phi_best_map",
                    "Map of best triggers per station",20,1,21,12,0,12);
       dttpgphbestmapcorr = dbe->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 = dbe->book2D("DT_TPG_theta_map",
                    "Map of triggers per station",15,1,16,12,0,12);
       dttpgthmaph = dbe->book2D("DT_TPG_theta_map_h",
                 "Map of H quality triggers per station",15,1,16,12,0,12);
       dttpgthbestmap = dbe->book2D("DT_TPG_theta_best_map",
                    "Map of besttriggers per station",15,1,16,12,0,12);
       dttpgthbestmaph = dbe->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::endJob ( void )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 142 of file L1TDTTPG.cc.

References gather_cfg::cout, dbe, nev_, outputFile_, and verbose_.

 {
   if(verbose_) cout << "L1TDTTPG: end job...." << endl;
   LogInfo("EndJob") << "analyzed " << nev_ << " events"; 
 
   if ( outputFile_.size() != 0  && dbe ) dbe->save(outputFile_);
 
   return;
 }

void L1TDTTPG::setMapPhLabel ( MonitorElement * me )

private

Definition at line 668 of file L1TDTTPG.cc.

References i, MonitorElement::setAxisTitle(), and MonitorElement::setBinLabel().

Referenced by beginRun().

 {
 
   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 680 of file L1TDTTPG.cc.

References i, MonitorElement::setAxisTitle(), and MonitorElement::setBinLabel().

Referenced by beginRun().

 {
 
   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

DQMStore* L1TDTTPG::dbe

private

Definition at line 67 of file L1TDTTPG.h.

Referenced by beginRun(), endJob(), and L1TDTTPG().

MonitorElement* L1TDTTPG::dttf_p_phi[3]

private

Definition at line 101 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttf_p_pt[3]

private

Definition at line 102 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttf_p_q[3]

private

Definition at line 103 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttf_p_qual[3]

private

Definition at line 104 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgphbestmap

private

Definition at line 84 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgphbestmapcorr

private

Definition at line 85 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgphbx[8]

private

Definition at line 69 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgphbxcomp

private

Definition at line 70 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgphmap

private

Definition at line 80 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgphmap2nd

private

Definition at line 82 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgphmapbx[3]

private

Definition at line 81 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgphmapcorr

private

Definition at line 83 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgphntrack

private

Definition at line 79 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgphquality[3]

private

Definition at line 76 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgphsector[3]

private

Definition at line 72 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgphstation[3]

private

Definition at line 73 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgphts2tag[3]

private

Definition at line 77 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgphwheel[3]

private

Definition at line 71 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

edm::InputTag L1TDTTPG::dttpgSource_

private

Definition at line 113 of file L1TDTTPG.h.

Referenced by analyze(), and L1TDTTPG().

edm::EDGetTokenT<L1MuDTChambPhContainer> L1TDTTPG::dttpgSourcePhContainer_token_

private

Definition at line 111 of file L1TDTTPG.h.

Referenced by analyze().

edm::EDGetTokenT<L1MuDTChambThContainer> L1TDTTPG::dttpgSourceThContainer_token_

private

Definition at line 112 of file L1TDTTPG.h.

Referenced by analyze().

MonitorElement* L1TDTTPG::dttpgthbestmap

private

Definition at line 98 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgthbestmaph

private

Definition at line 99 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgthbx[3]

private

Definition at line 88 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgthmap

private

Definition at line 95 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgthmapbx[3]

private

Definition at line 96 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgthmaph

private

Definition at line 97 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgthntrack

private

Definition at line 94 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgthquality[3]

private

Definition at line 93 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgthsector[3]

private

Definition at line 90 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgthstation[3]

private

Definition at line 91 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgththeta[3]

private

Definition at line 92 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

MonitorElement* L1TDTTPG::dttpgthwheel[3]

private

Definition at line 89 of file L1TDTTPG.h.

Referenced by analyze(), and beginRun().

std::ofstream L1TDTTPG::logFile_

private

Definition at line 110 of file L1TDTTPG.h.

bool L1TDTTPG::monitorDaemon_

private

Definition at line 109 of file L1TDTTPG.h.

int L1TDTTPG::nev_

private

Definition at line 106 of file L1TDTTPG.h.

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

std::string L1TDTTPG::outputFile_

private

Definition at line 107 of file L1TDTTPG.h.

Referenced by endJob(), and L1TDTTPG().

std::string L1TDTTPG::trstring_

private

Definition at line 115 of file L1TDTTPG.h.

Referenced by L1TDTTPG().

edm::EDGetTokenT<L1MuDTTrackContainer> L1TDTTPG::trToken_

private

Definition at line 116 of file L1TDTTPG.h.

Referenced by analyze(), and L1TDTTPG().

bool L1TDTTPG::verbose_

private

Definition at line 108 of file L1TDTTPG.h.

Referenced by analyze(), endJob(), and L1TDTTPG().

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation