L1TDTTPG Class Reference

#include <L1TDTTPG.h>

Inheritance diagram for L1TDTTPG:

Public Member Functions
	L1TDTTPG (const edm::ParameterSet &ps)
virtual	~L1TDTTPG ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c)
void	beginJob (void)
void	endJob (void)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

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_
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]
ofstream	logFile_
bool	monitorDaemon_
int	nev_
std::string	outputFile_
bool	verbose_

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)

Detailed Description

Definition at line 32 of file L1TDTTPG.h.

Constructor & Destructor Documentation

L1TDTTPG::L1TDTTPG

(

const edm::ParameterSet &

ps

)

Definition at line 105 of file L1TDTTPG.cc.

References gather_cfg::cout, dbe, edm::ParameterSet::getUntrackedParameter(), NULL, cmsCodeRules.cppFunctionSkipper::operator, outputFile_, DQMStore::setCurrentFolder(), DQMStore::setVerbose(), and verbose_.

  : dttpgSource_( ps.getParameter< InputTag >("dttpgSource") )
{

  // 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 140 of file L1TDTTPG.cc.

{
}

Member Function Documentation

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

protectedvirtual

Implements edm::EDAnalyzer.

Definition at line 331 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_, dttpgthbestmap, dttpgthbestmaph, dttpgthbx, dttpgthmap, dttpgthmapbx, dttpgthmaph, dttpgthntrack, dttpgthquality, dttpgthsector, dttpgthstation, dttpgththeta, dttpgthwheel, MonitorElement::Fill(), edm::Event::getByLabel(), i, edm::HandleBase::isValid(), j, gen::k, edm::InputTag::label(), nev_, pos, edm::InputTag::process(), matplotRender::t, and verbose_.

{

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

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

  edm::Handle<L1MuDTChambThContainer > myL1MuDTChambThContainer;  
  e.getByLabel(dttpgSource_,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* 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;

  
    std::string trstring;
    trstring = dttpgSource_.label()+":"+"DATA"+":"+dttpgSource_.process();
    edm::InputTag trInputTag(trstring);
    e.getByLabel(trInputTag,myL1MuDTTrackContainer);

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

  L1MuDTTrackContainer::TrackContainer *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 144 of file L1TDTTPG.cc.

References DQMStore::book1D(), DQMStore::book2D(), 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, nev_, cmsCodeRules.cppFunctionSkipper::operator, DQMStore::rmdir(), MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), setMapPhLabel(), and setMapThLabel().

{

  nev_ = 0;

  // get hold of back-end interface
  DQMStore* dbe = 0;
  dbe = Service<DQMStore>().operator->();
  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 321 of file L1TDTTPG.cc.

References gather_cfg::cout, dbe, nev_, outputFile_, DQMStore::save(), 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 681 of file L1TDTTPG.cc.

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

Referenced by beginJob().

{

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

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

Referenced by beginJob().

{

  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 58 of file L1TDTTPG.h.

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

MonitorElement* L1TDTTPG::dttf_p_phi[3]

private

Definition at line 92 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttf_p_pt[3]

private

Definition at line 93 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttf_p_q[3]

private

Definition at line 94 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttf_p_qual[3]

private

Definition at line 95 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgphbestmap

private

Definition at line 75 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgphbestmapcorr

private

Definition at line 76 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgphbx[8]

private

Definition at line 60 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgphbxcomp

private

Definition at line 61 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgphmap

private

Definition at line 71 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgphmap2nd

private

Definition at line 73 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgphmapbx[3]

private

Definition at line 72 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgphmapcorr

private

Definition at line 74 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgphntrack

private

Definition at line 70 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgphquality[3]

private

Definition at line 67 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgphsector[3]

private

Definition at line 63 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgphstation[3]

private

Definition at line 64 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgphts2tag[3]

private

Definition at line 68 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgphwheel[3]

private

Definition at line 62 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

edm::InputTag L1TDTTPG::dttpgSource_

private

Definition at line 102 of file L1TDTTPG.h.

Referenced by analyze().

MonitorElement* L1TDTTPG::dttpgthbestmap

private

Definition at line 89 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgthbestmaph

private

Definition at line 90 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgthbx[3]

private

Definition at line 79 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgthmap

private

Definition at line 86 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgthmapbx[3]

private

Definition at line 87 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgthmaph

private

Definition at line 88 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgthntrack

private

Definition at line 85 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgthquality[3]

private

Definition at line 84 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgthsector[3]

private

Definition at line 81 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgthstation[3]

private

Definition at line 82 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgththeta[3]

private

Definition at line 83 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TDTTPG::dttpgthwheel[3]

private

Definition at line 80 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

ofstream L1TDTTPG::logFile_

private

Definition at line 101 of file L1TDTTPG.h.

bool L1TDTTPG::monitorDaemon_

private

Definition at line 100 of file L1TDTTPG.h.

int L1TDTTPG::nev_

private

Definition at line 97 of file L1TDTTPG.h.

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

std::string L1TDTTPG::outputFile_

private

Definition at line 98 of file L1TDTTPG.h.

Referenced by endJob(), and L1TDTTPG().

bool L1TDTTPG::verbose_

private

Definition at line 99 of file L1TDTTPG.h.

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