|
|
|
#include <L1TDTTPG.h>
Definition at line 32 of file L1TDTTPG.h.
| L1TDTTPG::L1TDTTPG | ( | const edm::ParameterSet & | ps | ) |
Definition at line 111 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 146 of file L1TDTTPG.cc.
{
}
| 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(), 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 | ) | [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_, 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 | ) | [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);
}
}
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().
1.7.3