|
|
|
#include <L1TCSCTF.h>
Definition at line 48 of file L1TCSCTF.h.
| L1TCSCTF::L1TCSCTF | ( | const edm::ParameterSet & | ps | ) |
Definition at line 29 of file L1TCSCTF.cc.
References edm::ParameterSet::addUntrackedParameter(), dbe, Reference_intrackfit_cff::endcap, edm::ParameterSet::getUntrackedParameter(), NULL, cmsCodeRules::cppFunctionSkipper::operator, outputFile_, DQMStore::setCurrentFolder(), DQMStore::setVerbose(), srLUTs_, relativeConstraints::station, and verbose_.
: gmtProducer( ps.getParameter< InputTag >("gmtProducer") ), lctProducer( ps.getParameter< InputTag >("lctProducer") ), trackProducer( ps.getParameter< InputTag >("trackProducer") ), statusProducer( ps.getParameter< InputTag >("statusProducer") ), mbProducer( ps.getParameter< InputTag >("mbProducer") ) // KK_end { // verbosity switch verbose_ = ps.getUntrackedParameter<bool>("verbose", false); if(verbose_) edm::LogInfo("DataNotFound") << "L1TCSCTF: 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 ) { edm::LogInfo("DataNotFound") << "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/L1TCSCTF"); } // KK_start: instantiate standard on-fly SR LUTs from CSC TF emulator package bzero(srLUTs_,sizeof(srLUTs_)); int endcap=1, sector=1; // assume SR LUTs are all same for every sector in either of endcaps bool TMB07=true; // specific TMB firmware // Create a dumy pset for SR LUTs edm::ParameterSet srLUTset; srLUTset.addUntrackedParameter<bool>("ReadLUTs", false); srLUTset.addUntrackedParameter<bool>("Binary", false); srLUTset.addUntrackedParameter<std::string>("LUTPath", "./"); for(int station=1,fpga=0; station<=4 && fpga<5; station++) { if(station==1) for(int subSector=0; subSector<2 && fpga<5; subSector++) srLUTs_[fpga++] = new CSCSectorReceiverLUT(endcap, sector, subSector+1, station, srLUTset, TMB07); else srLUTs_[fpga++] = new CSCSectorReceiverLUT(endcap, sector, 0, station, srLUTset, TMB07); } // KK_end }
| L1TCSCTF::~L1TCSCTF | ( | ) | [virtual] |
| void L1TCSCTF::analyze | ( | const edm::Event & | e, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
Implements edm::EDAnalyzer.
Definition at line 276 of file L1TCSCTF.cc.
References BxInEvent_, csctfbx, csctfChamberOccupancies, csctfCoincL1ABXN, csctferrors, csctfHaloL1ABXN, csctfntrack, csctfoccupancies, csctfTrackEta, csctfTrackPhi, csctfTrackQ, cscTrackStubNumbers, DTstubsTimeTrackMenTimeArrival, Reference_intrackfit_cff::endcap, MonitorElement::Fill(), edm::EventSetup::get(), edm::Event::getByLabel(), L1MuGMTReadoutCollection::getRecords(), CSCSectorReceiverLUT::globalEtaME(), CSCSectorReceiverLUT::globalPhiME(), gmtProducer, haloDelEta23, i, edm::InputTag::instance(), isCSCcand_, edm::HandleBase::isValid(), L1ABXN, edm::InputTag::label(), lctProducer, CSCSectorReceiverLUT::localPhi(), M_PI, mbProducer, nev_, edm::Handle< T >::product(), CSCTriggerGeometry::setGeometry(), edm::shift, srLUTs_, relativeConstraints::station, ntuplemaker::status, statusProducer, trackProducer, testEve_cfg::tracks, CSCTriggerNumbering::triggerSubSectorFromLabels(), and verbose_.
{
int NumCSCTfTracksRep = 0;
nev_++;
if(verbose_) edm::LogInfo("DataNotFound") << "L1TCSCTF: analyze...." << endl;
edm::Handle<L1MuGMTReadoutCollection> pCollection;
if( gmtProducer.label() != "null" )
{ // GMT block
e.getByLabel(gmtProducer,pCollection);
if (!pCollection.isValid())
{
edm::LogInfo("DataNotFound") << "can't find L1MuGMTReadoutCollection with label "; // << csctfSource_.label() ;
return;
}
L1MuGMTReadoutCollection const* gmtrc = pCollection.product();
vector<L1MuGMTReadoutRecord> gmt_records = gmtrc->getRecords();
vector<L1MuGMTReadoutRecord>::const_iterator RRItr;
// GP
// Look if the readout window contains one (and only one CSC cands)
// to make it simpler I reject events with more than a CSC cand in the
// same readout window
// count non-empty candidates in this bx
int bxWindow = 0;
int nCands = 0;
for( RRItr = gmt_records.begin(); RRItr != gmt_records.end(); RRItr++ ) {
bxWindow++;
// get the csc candidates
vector<L1MuRegionalCand> INPCands = RRItr->getCSCCands();
vector<L1MuRegionalCand>::const_iterator INPItr;
BxInEvent_ = 0;
isCSCcand_ = false;
int nCandsBx = 0;
for( INPItr = INPCands.begin(); INPItr != INPCands.end(); ++INPItr ) {
if(!INPItr->empty()) {
nCandsBx++;
nCands++;
BxInEvent_ = RRItr->getBxInEvent();
if (verbose_)
edm::LogInfo("DataNotFound") << "cand " << nCandsBx << " -> assigned CSCTF bx: " << INPItr->bx() << endl;
}
}
if (verbose_)
if(nCandsBx)
edm::LogInfo("DataNotFound") << nCandsBx << " cands in bx: " << BxInEvent_ << endl;
}
if (nCands != 1) return;
else isCSCcand_ = true;
if (verbose_)
edm::LogInfo("DataNotFound") << "bxWindow: " << bxWindow << endl;
// GP_end
int ncsctftrack = 0;
if (verbose_)
{
edm::LogInfo("DataNotFound") << "\tCSCTFCand ntrack " << ncsctftrack << endl;
}
} // end of GMT block
L1ABXN = -999;
if( statusProducer.label() != "null" )
{
edm::Handle<L1CSCStatusDigiCollection> status;
e.getByLabel(statusProducer.label(),statusProducer.instance(),status);
bool integrity=status->first, se=false, sm=false, bx=false, af=false, fmm=false;
int nStat = 0;
for(std::vector<L1CSCSPStatusDigi>::const_iterator stat=status->second.begin(); stat!=status->second.end(); stat++)
{
se |= stat->SEs()&0xFFF;
sm |= stat->SMs()&0xFFF;
bx |= stat->BXs()&0xFFF;
af |= stat->AFs()&0xFFF;
fmm|= stat->FMM()!=8;
//GP_start
if(stat->VPs() != 0) {
L1ABXN += stat->BXN();
nStat++;
}
//GP_end
}
// GP: compute the average
if(nStat!=0) L1ABXN /= nStat;
if(integrity) csctferrors->Fill(0.5);
if(se) csctferrors->Fill(1.5);
if(sm) csctferrors->Fill(2.5);
if(bx) csctferrors->Fill(3.5);
if(af) csctferrors->Fill(4.5);
if(fmm) csctferrors->Fill(5.5);
}
if( lctProducer.label() != "null" )
{
edm::ESHandle<CSCGeometry> pDD;
c.get<MuonGeometryRecord>().get( pDD );
CSCTriggerGeometry::setGeometry(pDD);
edm::Handle<CSCCorrelatedLCTDigiCollection> corrlcts;
e.getByLabel(lctProducer.label(),lctProducer.instance(),corrlcts);
for(CSCCorrelatedLCTDigiCollection::DigiRangeIterator csc=corrlcts.product()->begin(); csc!=corrlcts.product()->end(); csc++)
{
CSCCorrelatedLCTDigiCollection::Range range1 = corrlcts.product()->get((*csc).first);
for(CSCCorrelatedLCTDigiCollection::const_iterator lct=range1.first; lct!=range1.second; lct++)
{
int endcap = (*csc).first.endcap()-1;
int station = (*csc).first.station()-1;
int sector = (*csc).first.triggerSector()-1;
int subSector = CSCTriggerNumbering::triggerSubSectorFromLabels((*csc).first);
int cscId = (*csc).first.triggerCscId()-1;
int fpga = ( subSector ? subSector-1 : station+1 );
//JAG
int endcapAssignment = 1;
int shift = 1;
float sectorArg = sector;
//float sectorArg = j;
if( endcap == 1 ){
endcapAssignment = -1;
shift = 2;
//sectorArg = sector - 6;
}
int signedStation = (station + shift)* endcapAssignment;
if( (station == 0) && (endcap == 0)) signedStation = subSector - 1;
if( (station == 0) && (endcap == 1)) signedStation = (-1)*subSector;
float chamberArg1 = cscId * 0.1 + sectorArg;
//float chamberArg1 = i*0.1 + sectorArg;
//std::cout << "First" << i << " " << sectorArg << " " << chamberArg1 << std::endl;
float chamberArg11 = chamberArg1;
if(sectorArg == 1) chamberArg1 = chamberArg11 - 0.1;
if(sectorArg == 2) chamberArg1 = chamberArg11 - 0.2;
if(sectorArg == 3) chamberArg1 = chamberArg11 - 0.3;
if(sectorArg == 4) chamberArg1 = chamberArg11 - 0.4;
if(sectorArg == 5) chamberArg1 = chamberArg11 - 0.5;
//std::cout << "cscId, station, sector, endcap, sectorArg, chamber Arg: " << cscId << ", " << station << ", " <<sector << ", " << endcap << ", " << chamberArg1 << ", " << signedStation << std::endl;
csctfChamberOccupancies->Fill(chamberArg1, signedStation);
int bunchX = ( (lct->getBX()) - 6 );
int timingSectorArg = 3*(sector) + (lct->getMPCLink());
if( endcap == 1) timingSectorArg = 3*(sector + 6) + (lct->getMPCLink());
//std::cout << "Sector, MPCLink, TSA, endcap: " << sector << ", " << lct->getMPCLink() << ", " << timingSectorArg << ", " << endcap << std::endl;
csctfbx->Fill(timingSectorArg, bunchX );
//std::cout << "LCT'S, encap: " << endcap << ", station: " << station << ", sector: " << sector << ", subSector: " << subSector << ", cscId: " << cscId << std:: endl;
//End JAG
// Check if Det Id is within pysical range:
if( endcap<0||endcap>1 || sector<0||sector>6 || station<0||station>3 || cscId<0||cscId>8 || fpga<0||fpga>4)
{
edm::LogError("L1CSCTF: CSC TP are out of range: ")<<" endcap: "<<(endcap+1)<<" station: "<<(station+1) <<" sector: "<<(sector+1)<<" subSector: "<<subSector<<" fpga: "<<fpga<<" cscId: "<<(cscId+1);
continue;
}
lclphidat lclPhi;
try {
lclPhi = srLUTs_[fpga]->localPhi(lct->getStrip(), lct->getPattern(), lct->getQuality(), lct->getBend());
} catch(...) {
bzero(&lclPhi,sizeof(lclPhi));
}
gblphidat gblPhi;
try {
gblPhi = srLUTs_[fpga]->globalPhiME(lclPhi.phi_local, lct->getKeyWG(), cscId+1);
} catch(...) {
bzero(&gblPhi,sizeof(gblPhi));
}
gbletadat gblEta;
try {
gblEta = srLUTs_[fpga]->globalEtaME(lclPhi.phi_bend_local, lclPhi.phi_local, lct->getKeyWG(), cscId+1);
} catch(...) {
bzero(&gblEta,sizeof(gblEta));
}
// SR LUT gives packed eta and phi values -> normilize them to 1 by scale them to 'max' and shift by 'min'
//float etaP = gblEta.global_eta/127*1.5 + 0.9;
//float phiP = (gblPhi.global_phi);// + ( sector )*4096 + station*4096*12) * 1./(4*4096*12);
//std::cout << "LCT Eta & Phi Coordinates: " << etaP << ", " << phiP << "." << std::endl;
//csctfoccupancies->Fill( gblEta.global_eta/127. * 1.5 + 0.9, (gblPhi.global_phi + ( sector + (endcap?0:6) )*4096 + station*4096*12) * 1./(4*4096*12) );
}//lct != range1.scond
}//csc!=corrlcts.product()->end()
}// lctProducer.label() != "null"
if( trackProducer.label() != "null" )
{
edm::Handle<L1CSCTrackCollection> tracks;
e.getByLabel(trackProducer.label(),trackProducer.instance(),tracks);
for(L1CSCTrackCollection::const_iterator trk=tracks->begin(); trk<tracks->end(); trk++)
{
//JAG_START
NumCSCTfTracksRep++;
long LUTAdd = trk->first.ptLUTAddress();
int trigMode = ( (LUTAdd)&0xf0000 ) >> 16;
float etaReal = (trk->first.eta_packed() );
if( trk->first.endcap() != 1)
{
float holder = etaReal;
etaReal = -1*holder;
etaReal -= 1;
}
float phiReal = ((trk->first.localPhi())*(62.5/24) + (trk->first.sector() - 1)*60)*M_PI/180 ;//+ 15;
if( phiReal > (2*M_PI) ) phiReal -= (2*M_PI);
csctfoccupancies->Fill( etaReal, phiReal);
csctfTrackPhi->Fill(phiReal);
csctfTrackEta->Fill( etaReal );
//std::cout << "Eta, phi, trigger mode, sector: " << etaReal << ", " << phiReal << ", " << trigMode << ", " << trk->first.sector() << "." << std::endl;
csctfTrackQ->Fill( trigMode );
if( trigMode == 15 )
{
//GP_start
csctfHaloL1ABXN->Fill(L1ABXN);
//GP_end
double haloVals[4][4];
for( int i = 0; i < 4; i++)
{
haloVals[i][0] = 0;
}
edm::Handle<CSCCorrelatedLCTDigiCollection> corrlcts;
e.getByLabel(lctProducer.label(),lctProducer.instance(),corrlcts);
for(CSCCorrelatedLCTDigiCollection::DigiRangeIterator csc=corrlcts.product()->begin(); csc!=corrlcts.product()->end(); csc++)
{
CSCCorrelatedLCTDigiCollection::Range range1 = corrlcts.product()->get((*csc).first);
for(CSCCorrelatedLCTDigiCollection::const_iterator lct=range1.first; lct!=range1.second; lct++)
{
int endcap = (*csc).first.endcap()-1;
int station = (*csc).first.station()-1;
int sector = (*csc).first.triggerSector()-1;
int cscId = (*csc).first.triggerCscId()-1;
int subSector = CSCTriggerNumbering::triggerSubSectorFromLabels((*csc).first);
int fpga = ( subSector ? subSector-1 : station+1 );
if( (station == 1) || (station == 2) )
{
int modEnd = 1;
if( endcap == 0 ) modEnd = -1;
//if( endcap == 1 ) modEnd = 1;
int indexHalo = modEnd + station;
if(haloVals[indexHalo][0] == 1.0) haloVals[indexHalo][3] = 1.0;
if(haloVals[indexHalo][0] == 0) haloVals[indexHalo][0] = 1.0;
haloVals[indexHalo][1] = sector*1.0;
lclphidat lclPhi;
lclPhi = srLUTs_[fpga]->localPhi(lct->getStrip(), lct->getPattern(), lct->getQuality(), lct->getBend());
gblphidat gblPhi;
gblPhi = srLUTs_[fpga]->globalPhiME(lclPhi.phi_local, lct->getKeyWG(), cscId+1);
gbletadat gblEta;
gblEta = srLUTs_[fpga]->globalEtaME(lclPhi.phi_bend_local, lclPhi.phi_local, lct->getKeyWG(), cscId+1);
haloVals[indexHalo][2] = gblEta.global_eta/127. * 1.5 + 0.9;
} //station1 or 2
} //lct first to second
} //corrlcts
if( (haloVals[0][0] == 1.) && (haloVals[1][0] == 1.) && (haloVals[0][3] != 1.) && (haloVals[1][3] != 1.) )
{
if( haloVals[0][1] == haloVals[1][1] ){
double delEta23 = haloVals[1][2] - haloVals[0][2];
haloDelEta23->Fill( delEta23 );
}
}
if( (haloVals[2][0] == 1.) && (haloVals[3][0] == 1.) && (haloVals[2][3] != 1.) && (haloVals[3][3] != 1.) )
{
if( haloVals[2][1] == haloVals[3][1] ){
double delEta23 = haloVals[3][2] - haloVals[2][2];
haloDelEta23->Fill( delEta23 );
}
}
} //halo trigger
int cscTrackStub = 0;
CSCCorrelatedLCTDigiCollection lctsOfTracks=trk->second;
for(CSCCorrelatedLCTDigiCollection::DigiRangeIterator trackStub=lctsOfTracks.begin(); trackStub!=lctsOfTracks.end(); trackStub++)
{
CSCCorrelatedLCTDigiCollection::Range range2 = lctsOfTracks.get((*trackStub).first);
for(CSCCorrelatedLCTDigiCollection::const_iterator lct=range2.first; lct!=range2.second; lct++)
{
cscTrackStub++;
}
}
//
//edm::Handle<L1CSCTrack> tfInf;
//e.getByLabel(trackProducer.label(), tfInf);
cscTrackStubNumbers->Fill(cscTrackStub);
//JAG_END
//GP_start: only coincidence
if( trigMode > 0 && trigMode !=15 && trigMode != 11 )
csctfCoincL1ABXN->Fill(L1ABXN);
//GP_end
}
}
csctfntrack->Fill(NumCSCTfTracksRep);
// KK_end ///////////////////////////////////
// GP_Start
if( mbProducer.label() != "null" ) {
// handle to needed collections
edm::Handle<CSCTriggerContainer<csctf::TrackStub> > dtStubs;
e.getByLabel(mbProducer.label(), mbProducer.instance(), dtStubs);
edm::Handle<L1CSCTrackCollection> tracks;
e.getByLabel(trackProducer.label(),trackProducer.instance(),tracks);
// loop on the DT stubs
std::vector<csctf::TrackStub> vstubs = dtStubs->get();
for(std::vector<csctf::TrackStub>::const_iterator stub=vstubs.begin();
stub!=vstubs.end(); stub++){
if (verbose_) {
edm::LogInfo("DataNotFound") << "\n mbEndcap: " << stub->endcap();
edm::LogInfo("DataNotFound") << "\n stub->getStrip()[FLAG]: " << stub->getStrip();
edm::LogInfo("DataNotFound") << "\n stub->getKeyWG()[CAL]: " << stub->getKeyWG();
edm::LogInfo("DataNotFound") << "\n stub->BX(): " << stub->BX();
edm::LogInfo("DataNotFound") << "\n stub->sector(): " << stub->sector();
edm::LogInfo("DataNotFound") << "\n stub->subsector(): " << stub->subsector();
edm::LogInfo("DataNotFound") << "\n stub->station(): " << stub->station();
edm::LogInfo("DataNotFound") << "\n stub->phiPacked(): " << stub->phiPacked();
edm::LogInfo("DataNotFound") << "\n stub->getBend(): " << stub->getBend();
edm::LogInfo("DataNotFound") << "\n stub->getQuality(): " << stub->getQuality();
edm::LogInfo("DataNotFound") << "\n stub->cscid(): " << stub->cscid() << endl;
}
// define the sector ID
int mbId = (stub->endcap()==2) ? 6 : 0;
mbId += stub->sector();
// *** do not fill if CalMB variable is set ***
// horrible! They used the same class to write up the LCT and MB info,
// but given the MB does not have strip and WG they replaced this two
// with the flag and cal bits... :S
if (stub->getKeyWG() == 0) {// !CAL as Janos adviced
// if FLAG =1, muon belong to previous BX
int bxDT = stub->BX()-stub->getStrip(); // correct by the FLAG
int subDT = stub->subsector();
// Fill the event only if CSC had or would have triggered
if (isCSCcand_){
// look for tracks in the event and compare the matching DT stubs
int trkBX = 0;
for(L1CSCTrackCollection::const_iterator trk=tracks->begin();
trk<tracks->end(); trk++){
trkBX = trk->first.BX();
int trkId = (trk->first.endcap()==2) ? 6 : 0;
trkId += trk->first.sector();
if (verbose_){
edm::LogInfo("DataNotFound") << "\n trk BX: " << trkBX
<< " Sector: " << trkId
<< " SubSector: " << trk->first.subsector()
<< " Endcap: " << trk->first.endcap();
edm::LogInfo("DataNotFound") << "\n DT BX: " << stub->BX()
<< " Sector: " << mbId
<< " SubSector: " << stub->subsector()
<< " Endcap: " << stub->endcap() << endl;
}
if (mbId == trkId) {
if (verbose_) {
edm::LogInfo("DataNotFound") << " --> MATCH" << endl;
edm::LogInfo("DataNotFound") << "Fill :" << trkBX+6-bxDT << " -- " << subDT << " -- cands" << endl;
}
// DT bx ranges from 3 to 9
// trk bx ranges from -3 to 3
DTstubsTimeTrackMenTimeArrival[mbId-1]->Fill(bxDT-trkBX-6,subDT);//subsec
}
}// loop on the tracks
}//if (isCSCcand_){
}//if (stub->getKeyWG() == 0) {
}
}
// GP_end
}
| void L1TCSCTF::beginJob | ( | void | ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 97 of file L1TCSCTF.cc.
References DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookProfile(), csctfbx, csctfChamberOccupancies, csctfCoincL1ABXN, csctferrors, csctfHaloL1ABXN, csctfntrack, csctfoccupancies, csctfTrackEta, csctfTrackPhi, csctfTrackQ, cscTrackStubNumbers, dbe, DTstubsTimeTrackMenTimeArrival, MonitorElement::getTH2F(), haloDelEta23, i, M_PI, nev_, cmsCodeRules::cppFunctionSkipper::operator, DQMStore::rmdir(), MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), and DQMStore::setCurrentFolder().
{
nev_ = 0;
// get hold of back-end interface
DQMStore* dbe = 0;
dbe = Service<DQMStore>().operator->();
if( dbe )
{
dbe->setCurrentFolder("L1T/L1TCSCTF");
dbe->rmdir("L1T/L1TCSCTF");
}
if( dbe )
{
dbe->setCurrentFolder("L1T/L1TCSCTF");
// KK_start: declaration of two monitoring histograms
// Error counting histogram:
// 1) checks TF data integrity (error rate - first bin),
// 2) monitors sychronization on input links (4 errors types: SE/SM/BX/AF; ORed for all time bins, links, and SPs),
// 3) reports FMM status (if in any SP FMM status != "Ready" - fill the last bin)
csctferrors = dbe->book1D("CSCTF_errors","CSCTF Errors",6,0,6);
csctferrors->setAxisTitle("Error type",1);
csctferrors->setAxisTitle("Number of Errors",2);
csctferrors->setBinLabel(1,"Corruptions",1);
csctferrors->setBinLabel(2,"Synch. Err.",1);
csctferrors->setBinLabel(3,"Synch. Mod.",1);
csctferrors->setBinLabel(4,"BX mismatch",1);
csctferrors->setBinLabel(5,"Time misalign.",1);
csctferrors->setBinLabel(6,"FMM != Ready",1);
// Occupancy histogram Eta x Y, where Y:
// 1) Phi_packed of input LCTs from 1st, 2nd, 3rd, and 4th stations
// 2) Phi_packed of output tracks
// (all 12 SPs - 360 degree coveradge)
//csctfoccupancies = dbe->book2D("CSCTF_occupancies","CSCTF Occupancies",100,0.8,2.5,1229,0,1.2);
//JAG MOD
csctfoccupancies = dbe->book2D("CSCTF_occupancies", "CSCTF Occupancies", 64,-32,31,32,0,6.2);
csctfoccupancies->setAxisTitle("#eta",1);
csctfoccupancies->setAxisTitle("#phi",2);
csctfoccupancies->setBinLabel(1, "-2.5", 1);
csctfoccupancies->setBinLabel(16,"-1.3", 1);
csctfoccupancies->setBinLabel(32,"-0.9", 1);
csctfoccupancies->setBinLabel(33,"0.9", 1);
csctfoccupancies->setBinLabel(45,"1.3", 1);
csctfoccupancies->setBinLabel(64,"2.5", 1);
// KK_end
//JAG
haloDelEta23 = dbe->book1D("CSCTF_Halo_Eta23","Delta station 2 to station 3 Eta for Halo Muons", 40, -0.20,0.30);
csctfTrackQ = dbe->book1D("CSCTF_Track_Q","CSC Track Mode", 16, -0.5, 15.5);
csctfTrackQ->setAxisTitle("Track Type", 1);
csctfTrackQ->setBinLabel(1,"No Track",1);
csctfTrackQ->setBinLabel(2,"Bad Phi Road",1);
csctfTrackQ->setBinLabel(3,"ME1-2-3",1);
csctfTrackQ->setBinLabel(4,"ME1-2-4",1);
csctfTrackQ->setBinLabel(5,"ME1-3-4",1);
csctfTrackQ->setBinLabel(6,"ME2-3-4",1);
csctfTrackQ->setBinLabel(7,"ME1-2",1);
csctfTrackQ->setBinLabel(8,"ME1-3",1);
csctfTrackQ->setBinLabel(9,"ME2-3",1);
csctfTrackQ->setBinLabel(10,"ME2-4",1);
csctfTrackQ->setBinLabel(11,"ME3-4",1);
csctfTrackQ->setBinLabel(12,"Singles",1);
csctfTrackQ->setBinLabel(13,"ME1-2,MB1",1);
csctfTrackQ->setBinLabel(14,"ME1-4",1);
csctfTrackQ->setBinLabel(15,"ME2,MB1",1);
csctfTrackQ->setBinLabel(16,"Halo Trigger",1);
csctfChamberOccupancies = dbe->book2D("CSCTF_Chamber_Occupancies","CSCTF Chamber Occupancies", 54, -0.05, 5.35, 10, -5.5, 4.5);
csctfChamberOccupancies->setAxisTitle("Sector (Endcap), (chambers 1-9 not labeled)",1);
csctfChamberOccupancies->setBinLabel(1,"ME-4",2);
csctfChamberOccupancies->setBinLabel(2,"ME-3",2);
csctfChamberOccupancies->setBinLabel(3,"ME-2",2);
csctfChamberOccupancies->setBinLabel(4,"ME-1b",2);
csctfChamberOccupancies->setBinLabel(5,"ME-1a",2);
csctfChamberOccupancies->setBinLabel(6,"ME+1a",2);
csctfChamberOccupancies->setBinLabel(7,"ME+1b",2);
csctfChamberOccupancies->setBinLabel(8,"ME+2",2);
csctfChamberOccupancies->setBinLabel(9,"ME+3",2);
csctfChamberOccupancies->setBinLabel(10,"ME+4",2);
csctfChamberOccupancies->setBinLabel(1, "1(+), 7(-)",1);
csctfChamberOccupancies->setBinLabel(10,"2(+), 8(-)",1);
csctfChamberOccupancies->setBinLabel(19,"3(+), 9(-)",1);
csctfChamberOccupancies->setBinLabel(28,"4(+), 10(-)",1);
csctfChamberOccupancies->setBinLabel(37,"5(+), 11(-)",1);
csctfChamberOccupancies->setBinLabel(46,"6(+), 12(-)",1);
csctfTrackPhi = dbe->book1D("CSCTF_Track_Phi", "CSCTF Track Phi",144,0,2*M_PI);
//csctfTrackPhi = dbe->book1D("CSCTF_Track_Phi", "CSCTF Track Phi",32,0,32);
csctfTrackPhi->setAxisTitle("Track #phi", 1);
csctfTrackEta = dbe->book1D("CSCTF_Track_Eta", "CSCTF Track Eta",64,-32,32);
//csctfTrackEta = dbe->book1D("CSCTF_Track_Eta", "CSCTF TrackEta",32,0,32);
csctfTrackEta->setAxisTitle("Track #eta", 1);
csctfTrackEta->setBinLabel(1, "-2.5", 1);
csctfTrackEta->setBinLabel(16,"-1.3", 1);
csctfTrackEta->setBinLabel(32,"-0.9", 1);
csctfTrackEta->setBinLabel(33,"0.9", 1);
csctfTrackEta->setBinLabel(45,"1.3", 1);
csctfTrackEta->setBinLabel(64,"2.5", 1);
csctfbx = dbe->bookProfile("CSCTF_bx","CSCTF bx", 36, 0.5, 36.5, 15, -5.5,9.5 ) ;
csctfbx->setAxisTitle("Sector, Endcap, MPC Link", 1);
csctfbx->setBinLabel(1,"1, +",1);
csctfbx->setBinLabel(4,"2, +",1);
csctfbx->setBinLabel(7,"3, +",1);
csctfbx->setBinLabel(10,"4, +",1);
csctfbx->setBinLabel(13,"5, +",1);
csctfbx->setBinLabel(16,"6, +",1);
csctfbx->setBinLabel(19,"1, -",1);
csctfbx->setBinLabel(22,"2, -",1);
csctfbx->setBinLabel(25,"3, -",1);
csctfbx->setBinLabel(28,"4, -",1);
csctfbx->setBinLabel(31,"5, -",1);
csctfbx->setBinLabel(34,"6, -",1);
cscTrackStubNumbers = dbe->book1D("CSCTF_TrackStubs", "Number of Stubs in CSC Tracks", 5, 0.5, 5.5);
csctfntrack = dbe->book1D("CSCTF_ntrack","Number of CSCTracks found per event", 5, 0.5, 5.5 ) ;
//JAG
}
//GP_start
char hname [200];
char htitle[200];
for(int i=0; i<12; i++) {
sprintf(hname ,"DTstubsTimeTrackMenTimeArrival_%d",i+1);
sprintf(htitle,"T_{track} - T_{DT stub} sector %d",i+1);
DTstubsTimeTrackMenTimeArrival[i] = dbe->book2D(hname,htitle, 7,-3,3, 2,1,3);
DTstubsTimeTrackMenTimeArrival[i]->getTH2F()->SetMinimum(0);
// axis makeup
DTstubsTimeTrackMenTimeArrival[i]->setAxisTitle("bx_{CSC track} - bx_{DT stub}",1);
DTstubsTimeTrackMenTimeArrival[i]->setAxisTitle("subsector" ,2);
DTstubsTimeTrackMenTimeArrival[i]->setBinLabel(1,"-3",1);
DTstubsTimeTrackMenTimeArrival[i]->setBinLabel(2,"-2",1);
DTstubsTimeTrackMenTimeArrival[i]->setBinLabel(3,"-1",1);
DTstubsTimeTrackMenTimeArrival[i]->setBinLabel(4, "0",1);
DTstubsTimeTrackMenTimeArrival[i]->setBinLabel(5,"+1",1);
DTstubsTimeTrackMenTimeArrival[i]->setBinLabel(6,"+2",1);
DTstubsTimeTrackMenTimeArrival[i]->setBinLabel(7,"+3",1);
DTstubsTimeTrackMenTimeArrival[i]->setBinLabel(1,"sub1",2);
DTstubsTimeTrackMenTimeArrival[i]->setBinLabel(2,"sub2",2);
}
csctfHaloL1ABXN = dbe->book1D("csctfHaloL1ABXN" ,"L1A BXN (HALO)", 3564,0,3564);
csctfHaloL1ABXN->setAxisTitle("L1A (BX)",1);
csctfCoincL1ABXN = dbe->book1D("csctfCoincL1ABXN","L1A BXN (COINCIDENCE)", 3564,0,3564);
csctfCoincL1ABXN->setAxisTitle("L1A (BX)",1);
//GP_end
}
| void L1TCSCTF::endJob | ( | void | ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 265 of file L1TCSCTF.cc.
References dbe, nev_, outputFile_, DQMStore::save(), and verbose_.
{
if(verbose_) edm::LogInfo("DataNotFound") << "L1TCSCTF: end job...." << endl;
LogInfo("EndJob") << "analyzed " << nev_ << " events";
if ( outputFile_.size() != 0 && dbe ) dbe->save(outputFile_);
return;
}
int L1TCSCTF::BxInEvent_ [private] |
Definition at line 100 of file L1TCSCTF.h.
Referenced by analyze().
MonitorElement* L1TCSCTF::csctfAFerror [private] |
Definition at line 93 of file L1TCSCTF.h.
MonitorElement* L1TCSCTF::csctfbx [private] |
Definition at line 78 of file L1TCSCTF.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TCSCTF::csctfChamberOccupancies [private] |
Definition at line 88 of file L1TCSCTF.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TCSCTF::csctfCoincL1ABXN [private] |
Definition at line 104 of file L1TCSCTF.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TCSCTF::csctferrors [private] |
Definition at line 81 of file L1TCSCTF.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TCSCTF::csctfHaloL1ABXN [private] |
Definition at line 103 of file L1TCSCTF.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TCSCTF::csctfntrack [private] |
Definition at line 77 of file L1TCSCTF.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TCSCTF::csctfoccupancies [private] |
Definition at line 82 of file L1TCSCTF.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TCSCTF::csctfTrackEta [private] |
Definition at line 90 of file L1TCSCTF.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TCSCTF::csctfTrackPhi [private] |
Definition at line 89 of file L1TCSCTF.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TCSCTF::csctfTrackQ [private] |
Definition at line 92 of file L1TCSCTF.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TCSCTF::cscTrackStubNumbers [private] |
Definition at line 91 of file L1TCSCTF.h.
Referenced by analyze(), and beginJob().
DQMStore* L1TCSCTF::dbe [private] |
Definition at line 70 of file L1TCSCTF.h.
Referenced by beginJob(), endJob(), and L1TCSCTF().
MonitorElement* L1TCSCTF::DTstubsTimeTrackMenTimeArrival[12] [private] |
Definition at line 99 of file L1TCSCTF.h.
Referenced by analyze(), and beginJob().
edm::InputTag L1TCSCTF::gmtProducer [private] |
Definition at line 113 of file L1TCSCTF.h.
Referenced by analyze().
MonitorElement* L1TCSCTF::haloDelEta23 [private] |
Definition at line 87 of file L1TCSCTF.h.
Referenced by analyze(), and beginJob().
bool L1TCSCTF::isCSCcand_ [private] |
Definition at line 101 of file L1TCSCTF.h.
Referenced by analyze().
int L1TCSCTF::L1ABXN [private] |
Definition at line 105 of file L1TCSCTF.h.
Referenced by analyze().
edm::InputTag L1TCSCTF::lctProducer [private] |
Definition at line 113 of file L1TCSCTF.h.
Referenced by analyze().
ofstream L1TCSCTF::logFile_ [private] |
Definition at line 112 of file L1TCSCTF.h.
edm::InputTag L1TCSCTF::mbProducer [private] |
Definition at line 113 of file L1TCSCTF.h.
Referenced by analyze().
bool L1TCSCTF::monitorDaemon_ [private] |
Definition at line 111 of file L1TCSCTF.h.
int L1TCSCTF::nev_ [private] |
Definition at line 108 of file L1TCSCTF.h.
Referenced by analyze(), beginJob(), and endJob().
std::string L1TCSCTF::outputFile_ [private] |
Definition at line 109 of file L1TCSCTF.h.
Referenced by endJob(), and L1TCSCTF().
CSCSectorReceiverLUT* L1TCSCTF::srLUTs_[5] [private] |
Definition at line 83 of file L1TCSCTF.h.
Referenced by analyze(), L1TCSCTF(), and ~L1TCSCTF().
edm::InputTag L1TCSCTF::statusProducer [private] |
Definition at line 113 of file L1TCSCTF.h.
Referenced by analyze().
edm::InputTag L1TCSCTF::trackProducer [private] |
Definition at line 113 of file L1TCSCTF.h.
Referenced by analyze().
bool L1TCSCTF::verbose_ [private] |
Definition at line 110 of file L1TCSCTF.h.
Referenced by analyze(), endJob(), and L1TCSCTF().
1.7.3