CMS 3D CMS Logo

Public Member Functions | Protected Member Functions | Private Member Functions | Private Attributes

L1TDTTPG Class Reference

#include <L1TDTTPG.h>

Inheritance diagram for L1TDTTPG:
edm::EDAnalyzer

List of all members.

Public Member Functions

 L1TDTTPG (const edm::ParameterSet &ps)
virtual ~L1TDTTPG ()

Protected Member Functions

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

Private Member Functions

void setMapPhLabel (MonitorElement *me)
void setMapThLabel (MonitorElement *me)

Private Attributes

DQMStoredbe
MonitorElementdttf_p_phi [3]
MonitorElementdttf_p_pt [3]
MonitorElementdttf_p_q [3]
MonitorElementdttf_p_qual [3]
MonitorElementdttpgphbestmap
MonitorElementdttpgphbestmapcorr
MonitorElementdttpgphbx [8]
MonitorElementdttpgphbxcomp
MonitorElementdttpgphmap
MonitorElementdttpgphmap2nd
MonitorElementdttpgphmapbx [3]
MonitorElementdttpgphmapcorr
MonitorElementdttpgphntrack
MonitorElementdttpgphquality [3]
MonitorElementdttpgphsector [3]
MonitorElementdttpgphstation [3]
MonitorElementdttpgphts2tag [3]
MonitorElementdttpgphwheel [3]
edm::InputTag dttpgSource_
MonitorElementdttpgthbestmap
MonitorElementdttpgthbestmaph
MonitorElementdttpgthbx [3]
MonitorElementdttpgthmap
MonitorElementdttpgthmapbx [3]
MonitorElementdttpgthmaph
MonitorElementdttpgthntrack
MonitorElementdttpgthquality [3]
MonitorElementdttpgthsector [3]
MonitorElementdttpgthstation [3]
MonitorElementdttpgththeta [3]
MonitorElementdttpgthwheel [3]
ofstream logFile_
bool monitorDaemon_
int nev_
std::string outputFile_
bool verbose_

Detailed Description

Definition at line 32 of file L1TDTTPG.h.


Constructor & Destructor Documentation

L1TDTTPG::L1TDTTPG ( const edm::ParameterSet ps)

Definition at line 111 of file L1TDTTPG.cc.

References gather_cfg::cout, dbe, edm::ParameterSet::getUntrackedParameter(), NULL, 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 146 of file L1TDTTPG.cc.

{
}

Member Function Documentation

void L1TDTTPG::analyze ( const edm::Event e,
const edm::EventSetup c 
) [protected, virtual]

Implements edm::EDAnalyzer.

Definition at line 337 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(), lumiQTWidget::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  ) [protected, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 150 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_, 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  ) [protected, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 327 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 687 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 699 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().

Definition at line 92 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 93 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 94 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 95 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 75 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 76 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 60 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 61 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 71 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 73 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 72 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 74 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 70 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 67 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 63 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 64 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 68 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 62 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 102 of file L1TDTTPG.h.

Referenced by analyze().

Definition at line 89 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 90 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 79 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 86 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 87 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 88 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 85 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 84 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 81 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 82 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

Definition at line 83 of file L1TDTTPG.h.

Referenced by analyze(), and beginJob().

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().