|
|
|
#include <DQMOffline/Muon/src/CSCOfflineMonitor.h>
This program merely unpacks collections and fills a few simple histograms. The idea is to compare the histograms for one offline release and another and look for unexpected differences.
Michael Schmitt, Northwestern University, July 2007
Definition at line 60 of file CSCOfflineMonitor.h.
| CSCOfflineMonitor::CSCOfflineMonitor | ( | const edm::ParameterSet & | pset | ) |
Constructor.
Definition at line 17 of file CSCOfflineMonitor.cc.
References cscRecHitTag_, cscSegTag_, edm::ParameterSet::getParameter(), param, stripDigiTag_, and wireDigiTag_.
00017 { 00018 00019 param = pset; 00020 00021 stripDigiTag_ = pset.getParameter<edm::InputTag>("stripDigiTag"); 00022 wireDigiTag_ = pset.getParameter<edm::InputTag>("wireDigiTag"); 00023 cscRecHitTag_ = pset.getParameter<edm::InputTag>("cscRecHitTag"); 00024 cscSegTag_ = pset.getParameter<edm::InputTag>("cscSegTag"); 00025 00026 }
| CSCOfflineMonitor::~CSCOfflineMonitor | ( | ) | [virtual] |
| void CSCOfflineMonitor::analyze | ( | const edm::Event & | event, | |
| const edm::EventSetup & | eventSetup | |||
| ) | [virtual] |
Perform the real analysis.
Implements edm::EDAnalyzer.
Definition at line 426 of file CSCOfflineMonitor.cc.
References c, CSCDetId::chamber(), ChiSquaredProbability(), cscRecHitTag_, cscSegments_cfi::cscSegments, cscSegTag_, diff, e, CSCDetId::endcap(), MonitorElement::Fill(), fitX(), edm::EventSetup::get(), getSignal(), hORecHits, hOSegments, hOStrips, hOWires, hRHCodeBroad, hRHCodeNarrow, hRHLayer, hRHnrechits, hRHRatioQ, hRHSumQ, hRHTiming, hRHX, hRHY, hSChiSqProb, hSCodeBroad, hSCodeNarrow, hSGlobalPhi, hSGlobalTheta, hSnHits, hSnhitsAll, hSnSegments, hSResid, hSTheta, hStripAll, hStripCodeBroad, hStripCodeNarrow, hStripLayer, hStripNFired, hStripPed, hStripPedAll, hStripStrip, hWireAll, hWireCodeBroad, hWireCodeNarrow, hWireLayer, hWirenGroupsTotal, hWireTBinAll, hWireWire, i, index, it, j, prof2calltree::last, CSCDetId::layer(), PV3DBase< T, PVType, FrameType >::phi(), r, CSCDetId::ring(), s, edm::RangeMap< ID, C, P >::size(), funct::sqrt(), CSCDetId::station(), stripDigiTag_, theta, PV3DBase< T, PVType, FrameType >::theta(), dimuonsSequences_cff::threshold, GeomDet::toGlobal(), typeIndex(), wireDigiTag_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().
00426 { 00427 00428 // Variables for occupancy plots 00429 bool wireo[2][4][4][36]; 00430 bool stripo[2][4][4][36]; 00431 bool rechito[2][4][4][36]; 00432 bool segmento[2][4][4][36]; 00433 00434 for (int e = 0; e < 2; e++){ 00435 for (int s = 0; s < 4; s++){ 00436 for (int r = 0; r < 4; r++){ 00437 for (int c = 0; c < 36; c++){ 00438 wireo[e][s][r][c] = false; 00439 stripo[e][s][r][c] = false; 00440 rechito[e][s][r][c] = false; 00441 segmento[e][s][r][c] = false; 00442 } 00443 } 00444 } 00445 } 00446 00447 // ============================================== 00448 // 00449 // look at DIGIs 00450 // 00451 // ============================================== 00452 00453 00454 edm::Handle<CSCStripDigiCollection> strips; 00455 edm::Handle<CSCWireDigiCollection> wires; 00456 event.getByLabel( stripDigiTag_, strips); 00457 event.getByLabel( wireDigiTag_, wires); 00458 00459 // 00460 // WIRE GROUPS 00461 int nWireGroupsTotal = 0; 00462 for (CSCWireDigiCollection::DigiRangeIterator j=wires->begin(); j!=wires->end(); j++) { 00463 CSCDetId id = (CSCDetId)(*j).first; 00464 int index = typeIndex(id); 00465 int kEndcap = id.endcap(); 00466 int cEndcap = id.endcap(); 00467 if (kEndcap == 2) cEndcap = -1; 00468 int kRing = id.ring(); 00469 int kStation = id.station(); 00470 int kChamber = id.chamber(); 00471 int kLayer = id.layer(); 00472 std::vector<CSCWireDigi>::const_iterator digiItr = (*j).second.first; 00473 std::vector<CSCWireDigi>::const_iterator last = (*j).second.second; 00474 for( ; digiItr != last; ++digiItr) { 00475 int myWire = digiItr->getWireGroup(); 00476 int myTBin = digiItr->getTimeBin(); 00477 nWireGroupsTotal++; 00478 int kCodeBroad = cEndcap * ( 4*(kStation-1) + kRing) ; 00479 int kCodeNarrow = cEndcap * ( 100*(kRing-1) + kChamber) ; 00480 00481 // fill wire histos 00482 hWireTBinAll->Fill(myTBin); 00483 hWireAll->Fill(myWire); 00484 hWireCodeBroad->Fill(kCodeBroad); 00485 hWireCodeNarrow[kStation-1]->Fill(kCodeNarrow); 00486 hWireLayer[index]->Fill(kLayer); 00487 hWireWire[index]->Fill(myWire); 00488 00489 wireo[kEndcap-1][kStation-1][kRing-1][kChamber-1] = true; 00490 } 00491 } // end wire loop 00492 00493 00494 // 00495 // STRIPS 00496 // 00497 int nStripsFired = 0; 00498 for (CSCStripDigiCollection::DigiRangeIterator j=strips->begin(); j!=strips->end(); j++) { 00499 CSCDetId id = (CSCDetId)(*j).first; 00500 int index = typeIndex(id); 00501 int kEndcap = id.endcap(); 00502 int cEndcap = id.endcap(); 00503 if (kEndcap == 2) cEndcap = -1; 00504 int kRing = id.ring(); 00505 int kStation = id.station(); 00506 int kChamber = id.chamber(); 00507 int kLayer = id.layer(); 00508 std::vector<CSCStripDigi>::const_iterator digiItr = (*j).second.first; 00509 std::vector<CSCStripDigi>::const_iterator last = (*j).second.second; 00510 for( ; digiItr != last; ++digiItr) { 00511 int myStrip = digiItr->getStrip(); 00512 std::vector<int> myADCVals = digiItr->getADCCounts(); 00513 bool thisStripFired = false; 00514 float thisPedestal = 0.5*(float)(myADCVals[0]+myADCVals[1]); 00515 float threshold = 13.3 ; 00516 float diff = 0.; 00517 for (unsigned int iCount = 0; iCount < myADCVals.size(); iCount++) { 00518 diff = (float)myADCVals[iCount]-thisPedestal; 00519 if (diff > threshold) { thisStripFired = true; } 00520 } 00521 if (thisStripFired) { 00522 nStripsFired++; 00523 int kCodeBroad = cEndcap * ( 4*(kStation-1) + kRing) ; 00524 int kCodeNarrow = cEndcap * ( 100*(kRing-1) + kChamber) ; 00525 // fill strip histos 00526 hStripAll->Fill(myStrip); 00527 hStripCodeBroad->Fill(kCodeBroad); 00528 hStripCodeNarrow[kStation-1]->Fill(kCodeNarrow); 00529 hStripLayer[index]->Fill(kLayer); 00530 hStripStrip[index]->Fill(myStrip); 00531 00532 stripo[kEndcap-1][kStation-1][kRing-1][kChamber-1] = true; 00533 } 00534 } 00535 } // end strip loop 00536 00537 //======================================================= 00538 // 00539 // Look at the Pedestal Noise Distributions 00540 // 00541 //======================================================= 00542 00543 for (CSCStripDigiCollection::DigiRangeIterator j=strips->begin(); j!=strips->end(); j++) { 00544 CSCDetId id = (CSCDetId)(*j).first; 00545 int index = typeIndex(id); 00546 int kEndcap = id.endcap(); 00547 int cEndcap = id.endcap(); 00548 if (kEndcap == 2) cEndcap = -1; 00549 int kRing = id.ring(); 00550 int kStation = id.station(); 00551 //int kChamber = id.chamber(); 00552 //int kLayer = id.layer(); 00553 std::vector<CSCStripDigi>::const_iterator digiItr = (*j).second.first; 00554 std::vector<CSCStripDigi>::const_iterator last = (*j).second.second; 00555 for( ; digiItr != last; ++digiItr) { 00556 int myStrip = digiItr->getStrip(); 00557 std::vector<int> myADCVals = digiItr->getADCCounts(); 00558 float TotalADC = getSignal(*strips, id, myStrip); 00559 bool thisStripFired = false; 00560 float thisPedestal = 0.5*(float)(myADCVals[0]+myADCVals[1]); 00561 float thisSignal = (1./6)*(myADCVals[2]+myADCVals[3]+myADCVals[4]+myADCVals[5]+myADCVals[6]+myADCVals[7]); 00562 float threshold = 13.3; 00563 if(kStation == 1 && kRing == 4) 00564 { 00565 kRing = 1; 00566 if(myStrip <= 16) myStrip += 64; // no trapping for any bizarreness 00567 } 00568 //int globalStrip = cEndcap*( kStation*1000000 + kRing*100000 + kChamber*1000 + kLayer*100 + myStrip); 00569 if (TotalADC > threshold) { thisStripFired = true;} 00570 if (!thisStripFired){ 00571 float ADC = thisSignal - thisPedestal; 00572 hStripPedAll->Fill(ADC); 00573 //hPedvsStrip->Fill(globalStrip,ADC); 00574 hStripPed[index]->Fill(ADC); 00575 } 00576 } 00577 } 00578 00579 00580 00581 00582 // ============================================== 00583 // 00584 // look at RECHITs 00585 // 00586 // ============================================== 00587 00588 // Get the CSC Geometry : 00589 ESHandle<CSCGeometry> cscGeom; 00590 eventSetup.get<MuonGeometryRecord>().get(cscGeom); 00591 00592 // Get the RecHits collection : 00593 Handle<CSCRecHit2DCollection> recHits; 00594 event.getByLabel(cscRecHitTag_,recHits); 00595 int nRecHits = recHits->size(); 00596 00597 00598 // --------------------- 00599 // Loop over rechits 00600 // --------------------- 00601 int iHit = 0; 00602 00603 // Build iterator for rechits and loop : 00604 CSCRecHit2DCollection::const_iterator recIt; 00605 for (recIt = recHits->begin(); recIt != recHits->end(); recIt++) { 00606 iHit++; 00607 00608 // Find chamber with rechits in CSC 00609 CSCDetId idrec = (CSCDetId)(*recIt).cscDetId(); 00610 int index = typeIndex(idrec); 00611 int kEndcap = idrec.endcap(); 00612 int cEndcap = idrec.endcap(); 00613 if (kEndcap == 2) cEndcap = -1; 00614 int kRing = idrec.ring(); 00615 int kStation = idrec.station(); 00616 int kChamber = idrec.chamber(); 00617 int kLayer = idrec.layer(); 00618 00619 rechito[kEndcap-1][kStation-1][kRing-1][kChamber-1] = true; 00620 00621 // Store rechit as a Local Point: 00622 LocalPoint rhitlocal = (*recIt).localPosition(); 00623 float xreco = rhitlocal.x(); 00624 float yreco = rhitlocal.y(); 00625 00626 // Find the strip containing this hit 00627 //CSCRecHit2D::ChannelContainer hitstrips = (*recIt).channels(); 00628 //int nStrips = hitstrips.size(); 00629 //int centerid = nStrips/2 + 1; 00630 //int centerStrip = hitstrips[centerid - 1]; 00631 00632 // Find the charge associated with this hit 00633 00634 CSCRecHit2D::ADCContainer adcs = (*recIt).adcs(); 00635 int adcsize = adcs.size(); 00636 float rHSumQ = 0; 00637 float sumsides = 0; 00638 for (int i = 0; i < adcsize; i++){ 00639 if (i != 3 && i != 7 && i != 11){ 00640 rHSumQ = rHSumQ + adcs[i]; 00641 } 00642 if (adcsize == 12 && (i < 3 || i > 7) && i < 12){ 00643 sumsides = sumsides + adcs[i]; 00644 } 00645 } 00646 float rHratioQ = sumsides/rHSumQ; 00647 if (adcsize != 12) rHratioQ = -99; 00648 00649 // Get the signal timing of this hit 00650 float rHtime = (*recIt).tpeak()/50; 00651 //float rHtime = getTiming(*strips, idrec, centerStrip); 00652 00653 // Get pointer to the layer: 00654 //const CSCLayer* csclayer = cscGeom->layer( idrec ); 00655 00656 // Transform hit position from local chamber geometry to global CMS geom 00657 //GlobalPoint rhitglobal= csclayer->toGlobal(rhitlocal); 00658 //float grecx = rhitglobal.x(); 00659 //float grecy = rhitglobal.y(); 00660 00661 00662 // Simple occupancy variables 00663 int kCodeBroad = cEndcap * ( 4*(kStation-1) + kRing) ; 00664 int kCodeNarrow = cEndcap * ( 100*(kRing-1) + kChamber) ; 00665 00666 // Fill some histograms 00667 hRHCodeBroad->Fill(kCodeBroad); 00668 hRHCodeNarrow[kStation-1]->Fill(kCodeNarrow); 00669 //hRHGlobal[(kStation-1) + ((kEndcap - 1)*4)]->Fill(grecx,grecy); 00670 hRHLayer[index]->Fill(kLayer); 00671 hRHX[index]->Fill(xreco); 00672 hRHY[index]->Fill(yreco); 00673 hRHSumQ[index]->Fill(rHSumQ); 00674 hRHRatioQ[index]->Fill(rHratioQ); 00675 hRHTiming[index]->Fill(rHtime); 00676 00677 } //end rechit loop 00678 00679 // ============================================== 00680 // 00681 // look at SEGMENTs 00682 // 00683 // =============================================== 00684 00685 // get CSC segment collection 00686 Handle<CSCSegmentCollection> cscSegments; 00687 event.getByLabel(cscSegTag_, cscSegments); 00688 int nSegments = cscSegments->size(); 00689 00690 // ----------------------- 00691 // loop over segments 00692 // ----------------------- 00693 int iSegment = 0; 00694 for(CSCSegmentCollection::const_iterator it=cscSegments->begin(); it != cscSegments->end(); it++) { 00695 iSegment++; 00696 // 00697 CSCDetId id = (CSCDetId)(*it).cscDetId(); 00698 int index = typeIndex(id); 00699 int kEndcap = id.endcap(); 00700 int cEndcap = id.endcap(); 00701 if (kEndcap == 2) cEndcap = -1; 00702 int kRing = id.ring(); 00703 int kStation = id.station(); 00704 int kChamber = id.chamber(); 00705 segmento[kEndcap-1][kStation-1][kRing-1][kChamber-1] = true; 00706 00707 00708 // 00709 float chisq = (*it).chi2(); 00710 int nhits = (*it).nRecHits(); 00711 int nDOF = 2*nhits-4; 00712 double chisqProb = ChiSquaredProbability( (double)chisq, nDOF ); 00713 LocalPoint localPos = (*it).localPosition(); 00714 LocalVector segDir = (*it).localDirection(); 00715 double theta = segDir.theta(); 00716 00717 // 00718 // try to get the CSC recHits that contribute to this segment. 00719 std::vector<CSCRecHit2D> theseRecHits = (*it).specificRecHits(); 00720 int nRH = (*it).nRecHits(); 00721 int jRH = 0; 00722 HepMatrix sp(6,1); 00723 HepMatrix se(6,1); 00724 for ( vector<CSCRecHit2D>::const_iterator iRH = theseRecHits.begin(); iRH != theseRecHits.end(); iRH++) { 00725 jRH++; 00726 CSCDetId idRH = (CSCDetId)(*iRH).cscDetId(); 00727 int kRing = idRH.ring(); 00728 int kStation = idRH.station(); 00729 int kLayer = idRH.layer(); 00730 00731 // Find the strip containing this hit 00732 CSCRecHit2D::ChannelContainer hitstrips = (*iRH).channels(); 00733 int nStrips = hitstrips.size(); 00734 int centerid = nStrips/2 + 1; 00735 int centerStrip = hitstrips[centerid - 1]; 00736 00737 // If this segment has 6 hits, find the position of each hit on the strip in units of stripwidth and store values 00738 if (nRH == 6){ 00739 float stpos = (*iRH).positionWithinStrip(); 00740 se(kLayer,1) = (*iRH).errorWithinStrip(); 00741 // Take into account half-strip staggering of layers (ME1/1 has no staggering) 00742 if (kStation == 1 && (kRing == 1 || kRing == 4)) sp(kLayer,1) = stpos + centerStrip; 00743 else{ 00744 if (kLayer == 1 || kLayer == 3 || kLayer == 5) sp(kLayer,1) = stpos + centerStrip; 00745 if (kLayer == 2 || kLayer == 4 || kLayer == 6) sp(kLayer,1) = stpos - 0.5 + centerStrip; 00746 } 00747 } 00748 00749 } 00750 00751 float residual = -99; 00752 // Fit all points except layer 3, then compare expected value for layer 3 to reconstructed value 00753 if (nRH == 6){ 00754 float expected = fitX(sp,se); 00755 residual = expected - sp(3,1); 00756 } 00757 00758 // global transformation 00759 float globX = 0.; 00760 float globY = 0.; 00761 float globZ = 0.; 00762 float globpPhi = 0.; 00763 float globR = 0.; 00764 float globTheta = 0.; 00765 float globPhi = 0.; 00766 const CSCChamber* cscchamber = cscGeom->chamber(id); 00767 if (cscchamber) { 00768 GlobalPoint globalPosition = cscchamber->toGlobal(localPos); 00769 globX = globalPosition.x(); 00770 globY = globalPosition.y(); 00771 globZ = globalPosition.z(); 00772 globpPhi = globalPosition.phi(); 00773 globR = sqrt(globX*globX + globY*globY); 00774 GlobalVector globalDirection = cscchamber->toGlobal(segDir); 00775 globTheta = globalDirection.theta(); 00776 globPhi = globalDirection.phi(); 00777 } 00778 00779 // Simple occupancy variables 00780 int kCodeBroad = cEndcap * ( 4*(kStation-1) + kRing) ; 00781 int kCodeNarrow = cEndcap * ( 100*(kRing-1) + kChamber) ; 00782 00783 00784 // Fill histos 00785 hSCodeBroad->Fill(kCodeBroad); 00786 hSChiSqProb->Fill(chisqProb); 00787 hSGlobalTheta->Fill(globTheta); 00788 hSGlobalPhi->Fill(globPhi); 00789 hSnhitsAll->Fill(nhits); 00790 hSCodeNarrow[kStation-1]->Fill(kCodeNarrow); 00791 //hSGlobal[(kStation-1) + ((kEndcap - 1)*4)]->Fill(globX,globY); 00792 hSnHits[index]->Fill(nhits); 00793 hSTheta[index]->Fill(theta); 00794 hSResid[index]->Fill(residual); 00795 00796 } // end segment loop 00797 00798 // Fill # per even histos (how many stips/wires/rechits/segments per event) 00799 hWirenGroupsTotal->Fill(nWireGroupsTotal); 00800 hStripNFired->Fill(nStripsFired); 00801 hSnSegments->Fill(nSegments); 00802 hRHnrechits->Fill(nRecHits); 00803 00804 // Fill occupancy plots 00805 for (int e = 0; e < 2; e++){ 00806 for (int s = 0; s < 4; s++){ 00807 for (int r = 0; r < 4; r++){ 00808 for (int c = 0; c < 36; c++){ 00809 int type = 0; 00810 if ((s+1) == 1) type = (r+1); 00811 else type = (s+1)*2 + (r+1); 00812 if ((e+1) == 1) type = type + 10; 00813 if ((e+1) == 2) type = 11 - type; 00814 //int bin = hOWires->GetBin(chamber,type); 00815 //hOWires->AddBinContent(bin); 00816 if (wireo[e][s][r][c]) hOWires->Fill((c+1),type); 00817 if (stripo[e][s][r][c]) hOStrips->Fill((c+1),type); 00818 if (rechito[e][s][r][c]) hORecHits->Fill((c+1),type); 00819 if (segmento[e][s][r][c]) hOSegments->Fill((c+1),type); 00820 } 00821 } 00822 } 00823 } 00824 00825 00826 }
| void CSCOfflineMonitor::beginJob | ( | edm::EventSetup const & | iSetup | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 28 of file CSCOfflineMonitor.cc.
References DQMStore::book1D(), DQMStore::book2D(), dbe, hORecHits, hOSegments, hOStrips, hOWires, hRHCodeBroad, hRHCodeNarrow, hRHLayer, hRHnrechits, hRHRatioQ, hRHSumQ, hRHTiming, hRHX, hRHY, hSChiSqProb, hSCodeBroad, hSCodeNarrow, hSGlobalPhi, hSGlobalTheta, hSnHits, hSnhitsAll, hSnSegments, hSResid, hSTheta, hStripAll, hStripCodeBroad, hStripCodeNarrow, hStripLayer, hStripNFired, hStripPed, hStripPedAll, hStripStrip, hWireAll, hWireCodeBroad, hWireCodeNarrow, hWireLayer, hWirenGroupsTotal, hWireTBinAll, hWireWire, and DQMStore::setCurrentFolder().
00028 { 00029 dbe = Service<DQMStore>().operator->(); 00030 00031 // wire digis 00032 dbe->setCurrentFolder("CSC/CSCOfflineMonitor/Digis"); 00033 hWireAll = dbe->book1D("hWireAll","all wire group numbers",121,-0.5,120.5); 00034 hWireTBinAll = dbe->book1D("hWireTBinAll","time bins all wires",21,-0.5,20.5); 00035 hWirenGroupsTotal = dbe->book1D("hWirenGroupsTotal","total number of wire groups",101,-0.5,100.5); 00036 hWireCodeBroad = dbe->book1D("hWireCodeBroad","broad scope code for wires",33,-16.5,16.5); 00037 hWireCodeNarrow.push_back(dbe->book1D("hWireCodeNarrow1","narrow scope wire code station 1",801,-400.5,400.5)); 00038 hWireCodeNarrow.push_back(dbe->book1D("hWireCodeNarrow2","narrow scope wire code station 2",801,-400.5,400.5)); 00039 hWireCodeNarrow.push_back(dbe->book1D("hWireCodeNarrow3","narrow scope wire code station 3",801,-400.5,400.5)); 00040 hWireCodeNarrow.push_back(dbe->book1D("hWireCodeNarrow4","narrow scope wire code station 4",801,-400.5,400.5)); 00041 hWireWire.push_back(dbe->book1D("hWireWire_m42","wire number ME -4/2",113,-0.5,112.5)); 00042 hWireWire.push_back(dbe->book1D("hWireWire_m41","wire number ME -4/1",113,-0.5,112.5)); 00043 hWireWire.push_back(dbe->book1D("hWireWire_m32","wire number ME -3/2",113,-0.5,112.5)); 00044 hWireWire.push_back(dbe->book1D("hWireWire_m31","wire number ME -3/1",113,-0.5,112.5)); 00045 hWireWire.push_back(dbe->book1D("hWireWire_m22","wire number ME -2/2",113,-0.5,112.5)); 00046 hWireWire.push_back(dbe->book1D("hWireWire_m21","wire number ME -2/1",113,-0.5,112.5)); 00047 hWireWire.push_back(dbe->book1D("hWireWire_m11a","wire number ME -1/1a",113,-0.5,112.5)); 00048 hWireWire.push_back(dbe->book1D("hWireWire_m13","wire number ME -1/3",113,-0.5,112.5)); 00049 hWireWire.push_back(dbe->book1D("hWireWire_m12","wire number ME -1/2",113,-0.5,112.5)); 00050 hWireWire.push_back(dbe->book1D("hWireWire_m11b","wire number ME -1/1b",113,-0.5,112.5)); 00051 hWireWire.push_back(dbe->book1D("hWireWire_p11b","wire number ME +1/1b",113,-0.5,112.5)); 00052 hWireWire.push_back(dbe->book1D("hWireWire_p12","wire number ME +1/2",113,-0.5,112.5)); 00053 hWireWire.push_back(dbe->book1D("hWireWire_p13","wire number ME +1/3",113,-0.5,112.5)); 00054 hWireWire.push_back(dbe->book1D("hWireWire_p11a","wire number ME +1/1a",113,-0.5,112.5)); 00055 hWireWire.push_back(dbe->book1D("hWireWire_p21","wire number ME +2/1",113,-0.5,112.5)); 00056 hWireWire.push_back(dbe->book1D("hWireWire_p22","wire number ME +2/2",113,-0.5,112.5)); 00057 hWireWire.push_back(dbe->book1D("hWireWire_p31","wire number ME +3/1",113,-0.5,112.5)); 00058 hWireWire.push_back(dbe->book1D("hWireWire_p32","wire number ME +3/2",113,-0.5,112.5)); 00059 hWireWire.push_back(dbe->book1D("hWireWire_p41","wire number ME +4/1",113,-0.5,112.5)); 00060 hWireWire.push_back(dbe->book1D("hWireWire_p42","wire number ME +4/2",113,-0.5,112.5)); 00061 hWireLayer.push_back(dbe->book1D("hWireLayer_m42","layer wire ME -4/2",7,-0.5,6.5)); 00062 hWireLayer.push_back(dbe->book1D("hWireLayer_m41","layer wire ME -4/1",7,-0.5,6.5)); 00063 hWireLayer.push_back(dbe->book1D("hWireLayer_m32","layer wire ME -3/2",7,-0.5,6.5)); 00064 hWireLayer.push_back(dbe->book1D("hWireLayer_m31","layer wire ME -3/1",7,-0.5,6.5)); 00065 hWireLayer.push_back(dbe->book1D("hWireLayer_m22","layer wire ME -2/2",7,-0.5,6.5)); 00066 hWireLayer.push_back(dbe->book1D("hWireLayer_m21","layer wire ME -2/1",7,-0.5,6.5)); 00067 hWireLayer.push_back(dbe->book1D("hWireLayer_m11a","layer wire ME -1/1a",7,-0.5,6.5)); 00068 hWireLayer.push_back(dbe->book1D("hWireLayer_m13","layer wire ME -1/3",7,-0.5,6.5)); 00069 hWireLayer.push_back(dbe->book1D("hWireLayer_m12","layer wire ME -1/2",7,-0.5,6.5)); 00070 hWireLayer.push_back(dbe->book1D("hWireLayer_m11b","layer wire ME -1/1b",7,-0.5,6.5)); 00071 hWireLayer.push_back(dbe->book1D("hWireLayer_p11b","layer wire ME +1/1b",7,-0.5,6.5)); 00072 hWireLayer.push_back(dbe->book1D("hWireLayer_p12","layer wire ME +1/2",7,-0.5,6.5)); 00073 hWireLayer.push_back(dbe->book1D("hWireLayer_p13","layer wire ME +1/3",7,-0.5,6.5)); 00074 hWireLayer.push_back(dbe->book1D("hWireLayer_p11a","layer wire ME +1/1a",7,-0.5,6.5)); 00075 hWireLayer.push_back(dbe->book1D("hWireLayer_p21","layer wire ME +2/1",7,-0.5,6.5)); 00076 hWireLayer.push_back(dbe->book1D("hWireLayer_p22","layer wire ME +2/2",7,-0.5,6.5)); 00077 hWireLayer.push_back(dbe->book1D("hWireLayer_p31","layer wire ME +3/1",7,-0.5,6.5)); 00078 hWireLayer.push_back(dbe->book1D("hWireLayer_p32","layer wire ME +3/2",7,-0.5,6.5)); 00079 hWireLayer.push_back(dbe->book1D("hWireLayer_p41","layer wire ME +4/1",7,-0.5,6.5)); 00080 hWireLayer.push_back(dbe->book1D("hWireLayer_p42","layer wire ME +4/2",7,-0.5,6.5)); 00081 00082 // strip digis 00083 hStripAll = dbe->book1D("hStripAll","all strip numbers",81,-0.5,80.5); 00084 hStripNFired = dbe->book1D("hStripNFired","total number of fired strips",601,-0.5,600.5); 00085 hStripCodeBroad = dbe->book1D("hStripCodeBroad","broad scope code for strips",33,-16.5,16.5); 00086 hStripCodeNarrow.push_back(dbe->book1D("hStripCodeNarrow1","narrow scope strip code station 1",801,-400.5,400.5)); 00087 hStripCodeNarrow.push_back(dbe->book1D("hStripCodeNarrow2","narrow scope strip code station 2",801,-400.5,400.5)); 00088 hStripCodeNarrow.push_back(dbe->book1D("hStripCodeNarrow3","narrow scope strip code station 3",801,-400.5,400.5)); 00089 hStripCodeNarrow.push_back(dbe->book1D("hStripCodeNarrow4","narrow scope strip code station 4",801,-400.5,400.5)); 00090 hStripLayer.push_back(dbe->book1D("hStripLayer_m42","layer strip ME -4/2",7,-0.5,6.5)); 00091 hStripLayer.push_back(dbe->book1D("hStripLayer_m41","layer strip ME -4/1",7,-0.5,6.5)); 00092 hStripLayer.push_back(dbe->book1D("hStripLayer_m32","layer strip ME -3/2",7,-0.5,6.5)); 00093 hStripLayer.push_back(dbe->book1D("hStripLayer_m31","layer strip ME -3/1",7,-0.5,6.5)); 00094 hStripLayer.push_back(dbe->book1D("hStripLayer_m22","layer strip ME -2/2",7,-0.5,6.5)); 00095 hStripLayer.push_back(dbe->book1D("hStripLayer_m21","layer strip ME -2/1",7,-0.5,6.5)); 00096 hStripLayer.push_back(dbe->book1D("hStripLayer_m11a","layer strip ME -1/1a",7,-0.5,6.5)); 00097 hStripLayer.push_back(dbe->book1D("hStripLayer_m13","layer strip ME -1/3",7,-0.5,6.5)); 00098 hStripLayer.push_back(dbe->book1D("hStripLayer_m12","layer strip ME -1/2",7,-0.5,6.5)); 00099 hStripLayer.push_back(dbe->book1D("hStripLayer_m11b","layer strip ME -1/1b",7,-0.5,6.5)); 00100 hStripLayer.push_back(dbe->book1D("hStripLayer_p11b","layer strip ME +1/1b",7,-0.5,6.5)); 00101 hStripLayer.push_back(dbe->book1D("hStripLayer_p12","layer strip ME +1/2",7,-0.5,6.5)); 00102 hStripLayer.push_back(dbe->book1D("hStripLayer_p13","layer strip ME +1/3",7,-0.5,6.5)); 00103 hStripLayer.push_back(dbe->book1D("hStripLayer_p11a","layer strip ME +1/1a",7,-0.5,6.5)); 00104 hStripLayer.push_back(dbe->book1D("hStripLayer_p21","layer strip ME +2/1",7,-0.5,6.5)); 00105 hStripLayer.push_back(dbe->book1D("hStripLayer_p22","layer strip ME +2/2",7,-0.5,6.5)); 00106 hStripLayer.push_back(dbe->book1D("hStripLayer_p31","layer strip ME +3/1",7,-0.5,6.5)); 00107 hStripLayer.push_back(dbe->book1D("hStripLayer_p32","layer strip ME +3/2",7,-0.5,6.5)); 00108 hStripLayer.push_back(dbe->book1D("hStripLayer_p41","layer strip ME +4/1",7,-0.5,6.5)); 00109 hStripLayer.push_back(dbe->book1D("hStripLayer_p42","layer strip ME +4/2",7,-0.5,6.5)); 00110 hStripStrip.push_back(dbe->book1D("hStripStrip_m42","strip number ME -4/2",81,-0.5,80.5)); 00111 hStripStrip.push_back(dbe->book1D("hStripStrip_m41","strip number ME -4/1",81,-0.5,80.5)); 00112 hStripStrip.push_back(dbe->book1D("hStripStrip_m32","strip number ME -3/2",81,-0.5,80.5)); 00113 hStripStrip.push_back(dbe->book1D("hStripStrip_m31","strip number ME -3/1",81,-0.5,80.5)); 00114 hStripStrip.push_back(dbe->book1D("hStripStrip_m22","strip number ME -2/2",81,-0.5,80.5)); 00115 hStripStrip.push_back(dbe->book1D("hStripStrip_m21","strip number ME -2/1",81,-0.5,80.5)); 00116 hStripStrip.push_back(dbe->book1D("hStripStrip_m11a","strip number ME -1/1a",81,-0.5,80.5)); 00117 hStripStrip.push_back(dbe->book1D("hStripStrip_m13","strip number ME -1/3",81,-0.5,80.5)); 00118 hStripStrip.push_back(dbe->book1D("hStripStrip_m12","strip number ME -1/2",81,-0.5,80.5)); 00119 hStripStrip.push_back(dbe->book1D("hStripStrip_m11b","strip number ME -1/1b",81,-0.5,80.5)); 00120 hStripStrip.push_back(dbe->book1D("hStripStrip_p11b","strip number ME +1/1b",81,-0.5,80.5)); 00121 hStripStrip.push_back(dbe->book1D("hStripStrip_p12","strip number ME +1/2",81,-0.5,80.5)); 00122 hStripStrip.push_back(dbe->book1D("hStripStrip_p13","strip number ME +1/3",81,-0.5,80.5)); 00123 hStripStrip.push_back(dbe->book1D("hStripStrip_p11a","strip number ME +1/1a",81,-0.5,80.5)); 00124 hStripStrip.push_back(dbe->book1D("hStripStrip_p21","strip number ME +2/1",81,-0.5,80.5)); 00125 hStripStrip.push_back(dbe->book1D("hStripStrip_p22","strip number ME +2/2",81,-0.5,80.5)); 00126 hStripStrip.push_back(dbe->book1D("hStripStrip_p31","strip number ME +3/1",81,-0.5,80.5)); 00127 hStripStrip.push_back(dbe->book1D("hStripStrip_p32","strip number ME +3/2",81,-0.5,80.5)); 00128 hStripStrip.push_back(dbe->book1D("hStripStrip_p41","strip number ME +4/1",81,-0.5,80.5)); 00129 hStripStrip.push_back(dbe->book1D("hStripStrip_p42","strip number ME +4/2",81,-0.5,80.5)); 00130 00131 00132 //Pedestal Noise Plots 00133 dbe->setCurrentFolder("CSC/CSCOfflineMonitor/PedestalNoise"); 00134 00135 hStripPedAll = dbe->book1D("hStripPed","Pedestal Noise Distribution",50,-25.,25.); 00136 00137 hStripPed.push_back(dbe->book1D("hStripPedMEm42","Pedestal Noise Distribution Chamber ME -4/2 ",50,-25.,25.)); 00138 hStripPed.push_back(dbe->book1D("hStripPedMEm41","Pedestal Noise Distribution Chamber ME -4/1 ",50,-25.,25.)); 00139 hStripPed.push_back(dbe->book1D("hStripPedMEm32","Pedestal Noise Distribution Chamber ME -3/2 ",50,-25.,25.)); 00140 hStripPed.push_back(dbe->book1D("hStripPedMEm31","Pedestal Noise Distribution Chamber ME -3/1 ",50,-25.,25.)); 00141 hStripPed.push_back(dbe->book1D("hStripPedMEm22","Pedestal Noise Distribution Chamber ME -2/2 ",50,-25.,25.)); 00142 hStripPed.push_back(dbe->book1D("hStripPedMEm21","Pedestal Noise Distribution Chamber ME -2/1 ",50,-25.,25.)); 00143 hStripPed.push_back(dbe->book1D("hStripPedMEm11a","Pedestal Noise Distribution Chamber ME -1/1 ",50,-25.,25.)); 00144 hStripPed.push_back(dbe->book1D("hStripPedMEm13","Pedestal Noise Distribution Chamber ME -1/3 ",50,-25.,25.)); 00145 hStripPed.push_back(dbe->book1D("hStripPedMEm12","Pedestal Noise Distribution Chamber ME -1/2 ",50,-25.,25.)); 00146 hStripPed.push_back(dbe->book1D("hStripPedMEm11b","Pedestal Noise Distribution Chamber ME -1/1 ",50,-25.,25.)); 00147 hStripPed.push_back(dbe->book1D("hStripPedMEp11b","Pedestal Noise Distribution Chamber ME +1/1 ",50,-25.,25.)); 00148 hStripPed.push_back(dbe->book1D("hStripPedMEp12","Pedestal Noise Distribution Chamber ME +1/2 ",50,-25.,25.)); 00149 hStripPed.push_back(dbe->book1D("hStripPedMEp13","Pedestal Noise Distribution Chamber ME +1/3 ",50,-25.,25.)); 00150 hStripPed.push_back(dbe->book1D("hStripPedMEp11a","Pedestal Noise Distribution Chamber ME +1/1 ",50,-25.,25.)); 00151 hStripPed.push_back(dbe->book1D("hStripPedMEp21","Pedestal Noise Distribution Chamber ME +2/1 ",50,-25.,25.)); 00152 hStripPed.push_back(dbe->book1D("hStripPedMEp22","Pedestal Noise Distribution Chamber ME +2/2 ",50,-25.,25.)); 00153 hStripPed.push_back(dbe->book1D("hStripPedMEp31","Pedestal Noise Distribution Chamber ME +3/1 ",50,-25.,25.)); 00154 hStripPed.push_back(dbe->book1D("hStripPedMEp32","Pedestal Noise Distribution Chamber ME +3/2 ",50,-25.,25.)); 00155 hStripPed.push_back(dbe->book1D("hStripPedMEp41","Pedestal Noise Distribution Chamber ME +4/1 ",50,-25.,25.)); 00156 hStripPed.push_back(dbe->book1D("hStripPedMEp42","Pedestal Noise Distribution Chamber ME +4/2 ",50,-25.,25.)); 00157 00158 //hPedvsStrip = dbe->book2D("hPedvsStrip","Pedestal Noise Distribution",4000000,1000000.,5000000.,50,-25.,25.); 00159 00160 // recHits 00161 dbe->setCurrentFolder("CSC/CSCOfflineMonitor/recHits"); 00162 hRHCodeBroad = dbe->book1D("hRHCodeBroad","broad scope code for recHits",33,-16.5,16.5); 00163 hRHCodeNarrow.push_back(dbe->book1D("hRHCodeNarrow1","narrow scope recHit code station 1",801,-400.5,400.5)); 00164 hRHCodeNarrow.push_back(dbe->book1D("hRHCodeNarrow2","narrow scope recHit code station 2",801,-400.5,400.5)); 00165 hRHCodeNarrow.push_back(dbe->book1D("hRHCodeNarrow3","narrow scope recHit code station 3",801,-400.5,400.5)); 00166 hRHCodeNarrow.push_back(dbe->book1D("hRHCodeNarrow4","narrow scope recHit code station 4",801,-400.5,400.5)); 00167 00168 hRHLayer.push_back(dbe->book1D("hRHLayerm42","layer recHit ME -4/2",7,-0.5,6.5)); 00169 hRHLayer.push_back(dbe->book1D("hRHLayerm41","layer recHit ME -4/1",7,-0.5,6.5)); 00170 hRHLayer.push_back(dbe->book1D("hRHLayerm32","layer recHit ME -3/2",7,-0.5,6.5)); 00171 hRHLayer.push_back(dbe->book1D("hRHLayerm31","layer recHit ME -3/1",7,-0.5,6.5)); 00172 hRHLayer.push_back(dbe->book1D("hRHLayerm22","layer recHit ME -2/2",7,-0.5,6.5)); 00173 hRHLayer.push_back(dbe->book1D("hRHLayerm21","layer recHit ME -2/1",7,-0.5,6.5)); 00174 hRHLayer.push_back(dbe->book1D("hRHLayerm11a","layer recHit ME -1/1a",7,-0.5,6.5)); 00175 hRHLayer.push_back(dbe->book1D("hRHLayerm13","layer recHit ME -1/3",7,-0.5,6.5)); 00176 hRHLayer.push_back(dbe->book1D("hRHLayerm12","layer recHit ME -1/2",7,-0.5,6.5)); 00177 hRHLayer.push_back(dbe->book1D("hRHLayerm11b","layer recHit ME -1/1b",7,-0.5,6.5)); 00178 hRHLayer.push_back(dbe->book1D("hRHLayerp11b","layer recHit ME +1/1b",7,-0.5,6.5)); 00179 hRHLayer.push_back(dbe->book1D("hRHLayerp12","layer recHit ME +1/2",7,-0.5,6.5)); 00180 hRHLayer.push_back(dbe->book1D("hRHLayerp13","layer recHit ME +1/3",7,-0.5,6.5)); 00181 hRHLayer.push_back(dbe->book1D("hRHLayerp11a","layer recHit ME +1/1a",7,-0.5,6.5)); 00182 hRHLayer.push_back(dbe->book1D("hRHLayerp21","layer recHit ME +2/1",7,-0.5,6.5)); 00183 hRHLayer.push_back(dbe->book1D("hRHLayerp22","layer recHit ME +2/2",7,-0.5,6.5)); 00184 hRHLayer.push_back(dbe->book1D("hRHLayerp31","layer recHit ME +3/1",7,-0.5,6.5)); 00185 hRHLayer.push_back(dbe->book1D("hRHLayerp32","layer recHit ME +3/2",7,-0.5,6.5)); 00186 hRHLayer.push_back(dbe->book1D("hRHLayerp41","layer recHit ME +4/1",7,-0.5,6.5)); 00187 hRHLayer.push_back(dbe->book1D("hRHLayerp42","layer recHit ME +4/2",7,-0.5,6.5)); 00188 00189 00190 hRHX.push_back(dbe->book1D("hRHXm42","local X recHit ME -4/2",160,-80.,80.)); 00191 hRHX.push_back(dbe->book1D("hRHXm41","local X recHit ME -4/1",160,-80.,80.)); 00192 hRHX.push_back(dbe->book1D("hRHXm32","local X recHit ME -3/2",160,-80.,80.)); 00193 hRHX.push_back(dbe->book1D("hRHXm31","local X recHit ME -3/1",160,-80.,80.)); 00194 hRHX.push_back(dbe->book1D("hRHXm22","local X recHit ME -2/2",160,-80.,80.)); 00195 hRHX.push_back(dbe->book1D("hRHXm21","local X recHit ME -2/1",160,-80.,80.)); 00196 hRHX.push_back(dbe->book1D("hRHXm11a","local X recHit ME -1/1a",120,-60.,60.)); 00197 hRHX.push_back(dbe->book1D("hRHXm13","local X recHit ME -1/3",120,-60.,60.)); 00198 hRHX.push_back(dbe->book1D("hRHXm12","local X recHit ME -1/2",120,-60.,60.)); 00199 hRHX.push_back(dbe->book1D("hRHXm11b","local X recHit ME -1/1b",120,-60.,60.)); 00200 hRHX.push_back(dbe->book1D("hRHXp11b","local X recHit ME +1/1b",120,-60.,60.)); 00201 hRHX.push_back(dbe->book1D("hRHXp12","local X recHit ME +1/2",120,-60.,60.)); 00202 hRHX.push_back(dbe->book1D("hRHXp13","local X recHit ME +1/3",120,-60.,60.)); 00203 hRHX.push_back(dbe->book1D("hRHXp11a","local X recHit ME +1/1a",120,-60.,60.)); 00204 hRHX.push_back(dbe->book1D("hRHXp21","local X recHit ME +2/1",160,-80.,80.)); 00205 hRHX.push_back(dbe->book1D("hRHXp22","local X recHit ME +2/2",160,-80.,80.)); 00206 hRHX.push_back(dbe->book1D("hRHXp31","local X recHit ME +3/1",160,-80.,80.)); 00207 hRHX.push_back(dbe->book1D("hRHXp32","local X recHit ME +3/2",160,-80.,80.)); 00208 hRHX.push_back(dbe->book1D("hRHXp41","local X recHit ME +4/1",160,-80.,80.)); 00209 hRHX.push_back(dbe->book1D("hRHXp42","local X recHit ME +4/2",160,-80.,80.)); 00210 00211 00212 hRHY.push_back(dbe->book1D("hRHYm42","local Y recHit ME -4/2",60,-180.,180.)); 00213 hRHY.push_back(dbe->book1D("hRHYm41","local Y recHit ME -4/1",60,-180.,180.)); 00214 hRHY.push_back(dbe->book1D("hRHYm32","local Y recHit ME -3/2",60,-180.,180.)); 00215 hRHY.push_back(dbe->book1D("hRHYm31","local Y recHit ME -3/1",60,-180.,180.)); 00216 hRHY.push_back(dbe->book1D("hRHYm22","local Y recHit ME -2/2",60,-180.,180.)); 00217 hRHY.push_back(dbe->book1D("hRHYm21","local Y recHit ME -2/1",60,-180.,180.)); 00218 hRHY.push_back(dbe->book1D("hRHYm11a","local Y recHit ME -1/1a",50,-100.,100.)); 00219 hRHY.push_back(dbe->book1D("hRHYm13","local Y recHit ME -1/3",50,-100.,100.)); 00220 hRHY.push_back(dbe->book1D("hRHYm12","local Y recHit ME -1/2",50,-100.,100.)); 00221 hRHY.push_back(dbe->book1D("hRHYm11b","local Y recHit ME -1/1b",50,-100.,100.)); 00222 hRHY.push_back(dbe->book1D("hRHYp11b","local Y recHit ME +1/1b",50,-100.,100.)); 00223 hRHY.push_back(dbe->book1D("hRHYp12","local Y recHit ME +1/2",50,-100.,100.)); 00224 hRHY.push_back(dbe->book1D("hRHYp13","local Y recHit ME +1/3",50,-100.,100.)); 00225 hRHY.push_back(dbe->book1D("hRHYp11a","local Y recHit ME +1/1a",50,-100.,100.)); 00226 hRHY.push_back(dbe->book1D("hRHYp21","local Y recHit ME +2/1",60,-180.,180.)); 00227 hRHY.push_back(dbe->book1D("hRHYp22","local Y recHit ME +2/2",60,-180.,180.)); 00228 hRHY.push_back(dbe->book1D("hRHYp31","local Y recHit ME +3/1",60,-180.,180.)); 00229 hRHY.push_back(dbe->book1D("hRHYp32","local Y recHit ME +3/2",60,-180.,180.)); 00230 hRHY.push_back(dbe->book1D("hRHYp41","local Y recHit ME +4/1",60,-180.,180.)); 00231 hRHY.push_back(dbe->book1D("hRHYp42","local Y recHit ME +4/2",60,-180.,180.)); 00232 00233 00234 //hRHGlobal.push_back(dbe->book2D("hRHGlobalp1","recHit global X,Y station +1",400,-800.,800.,400,-800.,800.)); 00235 //hRHGlobal.push_back(dbe->book2D("hRHGlobalp2","recHit global X,Y station +2",400,-800.,800.,400,-800.,800.)); 00236 //hRHGlobal.push_back(dbe->book2D("hRHGlobalp3","recHit global X,Y station +3",400,-800.,800.,400,-800.,800.)); 00237 //hRHGlobal.push_back(dbe->book2D("hRHGlobalp4","recHit global X,Y station +4",400,-800.,800.,400,-800.,800.)); 00238 //hRHGlobal.push_back(dbe->book2D("hRHGlobalm1","recHit global X,Y station -1",400,-800.,800.,400,-800.,800.)); 00239 //hRHGlobal.push_back(dbe->book2D("hRHGlobalm2","recHit global X,Y station -2",400,-800.,800.,400,-800.,800.)); 00240 //hRHGlobal.push_back(dbe->book2D("hRHGlobalm3","recHit global X,Y station -3",400,-800.,800.,400,-800.,800.)); 00241 //hRHGlobal.push_back(dbe->book2D("hRHGlobalm4","recHit global X,Y station -4",400,-800.,800.,400,-800.,800.)); 00242 00243 00244 hRHSumQ.push_back(dbe->book1D("hRHSumQm42","Sum 3x3 recHit Charge (ME -4/2)",250,0,2000)); 00245 hRHSumQ.push_back(dbe->book1D("hRHSumQm41","Sum 3x3 recHit Charge (ME -4/1)",250,0,2000)); 00246 hRHSumQ.push_back(dbe->book1D("hRHSumQm32","Sum 3x3 recHit Charge (ME -3/2)",250,0,2000)); 00247 hRHSumQ.push_back(dbe->book1D("hRHSumQm31","Sum 3x3 recHit Charge (ME -3/1)",250,0,2000)); 00248 hRHSumQ.push_back(dbe->book1D("hRHSumQm22","Sum 3x3 recHit Charge (ME -2/2)",250,0,2000)); 00249 hRHSumQ.push_back(dbe->book1D("hRHSumQm21","Sum 3x3 recHit Charge (ME -2/1)",250,0,2000)); 00250 hRHSumQ.push_back(dbe->book1D("hRHSumQm11a","Sum 3x3 recHit Charge (ME -1/1a)",250,0,2000)); 00251 hRHSumQ.push_back(dbe->book1D("hRHSumQm13","Sum 3x3 recHit Charge (ME -1/3)",250,0,2000)); 00252 hRHSumQ.push_back(dbe->book1D("hRHSumQm12","Sum 3x3 recHit Charge (ME -1/2)",250,0,2000)); 00253 hRHSumQ.push_back(dbe->book1D("hRHSumQm11b","Sum 3x3 recHit Charge (ME -1/1b)",250,0,2000)); 00254 hRHSumQ.push_back(dbe->book1D("hRHSumQp11b","Sum 3x3 recHit Charge (ME +1/1b)",250,0,2000)); 00255 hRHSumQ.push_back(dbe->book1D("hRHSumQp12","Sum 3x3 recHit Charge (ME +1/2)",250,0,2000)); 00256 hRHSumQ.push_back(dbe->book1D("hRHSumQp13","Sum 3x3 recHit Charge (ME +1/3)",250,0,2000)); 00257 hRHSumQ.push_back(dbe->book1D("hRHSumQp11a","Sum 3x3 recHit Charge (ME +1/1a)",250,0,2000)); 00258 hRHSumQ.push_back(dbe->book1D("hRHSumQp21","Sum 3x3 recHit Charge (ME +2/1)",250,0,2000)); 00259 hRHSumQ.push_back(dbe->book1D("hRHSumQp22","Sum 3x3 recHit Charge (ME +2/2)",250,0,2000)); 00260 hRHSumQ.push_back(dbe->book1D("hRHSumQp31","Sum 3x3 recHit Charge (ME +3/1)",250,0,2000)); 00261 hRHSumQ.push_back(dbe->book1D("hRHSumQp32","Sum 3x3 recHit Charge (ME +3/2)",250,0,2000)); 00262 hRHSumQ.push_back(dbe->book1D("hRHSumQp41","Sum 3x3 recHit Charge (ME +4/1)",250,0,2000)); 00263 hRHSumQ.push_back(dbe->book1D("hRHSumQp42","Sum 3x3 recHit Charge (ME +4/2)",250,0,2000)); 00264 00265 00266 hRHRatioQ.push_back(dbe->book1D("hRHRatioQm42","Ratio (Ql+Qr)/Qt (ME -4/2)",120,-0.1,1.1)); 00267 hRHRatioQ.push_back(dbe->book1D("hRHRatioQm41","Ratio (Ql+Qr)/Qt (ME -4/1)",120,-0.1,1.1)); 00268 hRHRatioQ.push_back(dbe->book1D("hRHRatioQm32","Ratio (Ql+Qr)/Qt (ME -3/2)",120,-0.1,1.1)); 00269 hRHRatioQ.push_back(dbe->book1D("hRHRatioQm31","Ratio (Ql+Qr)/Qt (ME -3/1)",120,-0.1,1.1)); 00270 hRHRatioQ.push_back(dbe->book1D("hRHRatioQm22","Ratio (Ql+Qr)/Qt (ME -2/2)",120,-0.1,1.1)); 00271 hRHRatioQ.push_back(dbe->book1D("hRHRatioQm21","Ratio (Ql+Qr)/Qt (ME -2/1)",120,-0.1,1.1)); 00272 hRHRatioQ.push_back(dbe->book1D("hRHRatioQm11a","Ratio (Ql+Qr)/Qt (ME -1/1a)",120,-0.1,1.1)); 00273 hRHRatioQ.push_back(dbe->book1D("hRHRatioQm13","Ratio (Ql+Qr)/Qt (ME -1/3)",120,-0.1,1.1)); 00274 hRHRatioQ.push_back(dbe->book1D("hRHRatioQm12","Ratio (Ql+Qr)/Qt (ME -1/2)",120,-0.1,1.1)); 00275 hRHRatioQ.push_back(dbe->book1D("hRHRatioQm11b","Ratio (Ql+Qr)/Qt (ME -1/1b)",120,-0.1,1.1)); 00276 hRHRatioQ.push_back(dbe->book1D("hRHRatioQp11b","Ratio (Ql+Qr)/Qt (ME +1/1b)",120,-0.1,1.1)); 00277 hRHRatioQ.push_back(dbe->book1D("hRHRatioQp12","Ratio (Ql+Qr)/Qt (ME +1/2)",120,-0.1,1.1)); 00278 hRHRatioQ.push_back(dbe->book1D("hRHRatioQp13","Ratio (Ql+Qr)/Qt (ME +1/3)",120,-0.1,1.1)); 00279 hRHRatioQ.push_back(dbe->book1D("hRHRatioQp11a","Ratio (Ql+Qr)/Qt (ME +1/1a)",120,-0.1,1.1)); 00280 hRHRatioQ.push_back(dbe->book1D("hRHRatioQp21","Ratio (Ql+Qr)/Qt (ME +2/1)",120,-0.1,1.1)); 00281 hRHRatioQ.push_back(dbe->book1D("hRHRatioQp22","Ratio (Ql+Qr)/Qt (ME +2/2)",120,-0.1,1.1)); 00282 hRHRatioQ.push_back(dbe->book1D("hRHRatioQp31","Ratio (Ql+Qr)/Qt (ME +3/1)",120,-0.1,1.1)); 00283 hRHRatioQ.push_back(dbe->book1D("hRHRatioQp32","Ratio (Ql+Qr)/Qt (ME +3/2)",120,-0.1,1.1)); 00284 hRHRatioQ.push_back(dbe->book1D("hRHRatioQp41","Ratio (Ql+Qr)/Qt (ME +4/1)",120,-0.1,1.1)); 00285 hRHRatioQ.push_back(dbe->book1D("hRHRatioQp42","Ratio (Ql+Qr)/Qt (ME +4/2)",120,-0.1,1.1)); 00286 00287 00288 hRHTiming.push_back(dbe->book1D("hRHTimingm42","recHit Timing (ME -4/2)",100,0,10)); 00289 hRHTiming.push_back(dbe->book1D("hRHTimingm41","recHit Timing (ME -4/1)",100,0,10)); 00290 hRHTiming.push_back(dbe->book1D("hRHTimingm32","recHit Timing (ME -3/2)",100,0,10)); 00291 hRHTiming.push_back(dbe->book1D("hRHTimingm31","recHit Timing (ME -3/1)",100,0,10)); 00292 hRHTiming.push_back(dbe->book1D("hRHTimingm22","recHit Timing (ME -2/2)",100,0,10)); 00293 hRHTiming.push_back(dbe->book1D("hRHTimingm21","recHit Timing (ME -2/1)",100,0,10)); 00294 hRHTiming.push_back(dbe->book1D("hRHTimingm11a","recHit Timing (ME -1/1a)",100,0,10)); 00295 hRHTiming.push_back(dbe->book1D("hRHTimingm13","recHit Timing (ME -1/3)",100,0,10)); 00296 hRHTiming.push_back(dbe->book1D("hRHTimingm12","recHit Timing (ME -1/2)",100,0,10)); 00297 hRHTiming.push_back(dbe->book1D("hRHTimingm11b","recHit Timing (ME -1/1b)",100,0,10)); 00298 hRHTiming.push_back(dbe->book1D("hRHTimingp11b","recHit Timing (ME +1/1b)",100,0,10)); 00299 hRHTiming.push_back(dbe->book1D("hRHTimingp12","recHit Timing (ME +1/2)",100,0,10)); 00300 hRHTiming.push_back(dbe->book1D("hRHTimingp13","recHit Timing (ME +1/3)",100,0,10)); 00301 hRHTiming.push_back(dbe->book1D("hRHTimingp11a","recHit Timing (ME +1/1a)",100,0,10)); 00302 hRHTiming.push_back(dbe->book1D("hRHTimingp21","recHit Timing (ME +2/1)",100,0,10)); 00303 hRHTiming.push_back(dbe->book1D("hRHTimingp22","recHit Timing (ME +2/2)",100,0,10)); 00304 hRHTiming.push_back(dbe->book1D("hRHTimingp31","recHit Timing (ME +3/1)",100,0,10)); 00305 hRHTiming.push_back(dbe->book1D("hRHTimingp32","recHit Timing (ME +3/2)",100,0,10)); 00306 hRHTiming.push_back(dbe->book1D("hRHTimingp41","recHit Timing (ME +4/1)",100,0,10)); 00307 hRHTiming.push_back(dbe->book1D("hRHTimingp42","recHit Timing (ME +4/2)",100,0,10)); 00308 00309 hRHnrechits = dbe->book1D("hRHnrechits","recHits per Event (all chambers)",50,0,50); 00310 00311 // segments 00312 dbe->setCurrentFolder("CSC/CSCOfflineMonitor/Segments"); 00313 hSCodeBroad = dbe->book1D("hSCodeBroad","broad scope code for recHits",33,-16.5,16.5); 00314 hSCodeNarrow.push_back(dbe->book1D("hSCodeNarrow1","narrow scope Segment code station 1",801,-400.5,400.5)); 00315 hSCodeNarrow.push_back(dbe->book1D("hSCodeNarrow2","narrow scope Segment code station 2",801,-400.5,400.5)); 00316 hSCodeNarrow.push_back(dbe->book1D("hSCodeNarrow3","narrow scope Segment code station 3",801,-400.5,400.5)); 00317 hSCodeNarrow.push_back(dbe->book1D("hSCodeNarrow4","narrow scope Segment code station 4",801,-400.5,400.5)); 00318 00319 hSnHits.push_back(dbe->book1D("hSnHitsm42","N hits on Segments ME -4/2",7,-0.5,6.5)); 00320 hSnHits.push_back(dbe->book1D("hSnHitsm41","N hits on Segments ME -4/1",7,-0.5,6.5)); 00321 hSnHits.push_back(dbe->book1D("hSnHitsm32","N hits on Segments ME -3/2",7,-0.5,6.5)); 00322 hSnHits.push_back(dbe->book1D("hSnHitsm31","N hits on Segments ME -3/1",7,-0.5,6.5)); 00323 hSnHits.push_back(dbe->book1D("hSnHitsm22","N hits on Segments ME -2/2",7,-0.5,6.5)); 00324 hSnHits.push_back(dbe->book1D("hSnHitsm21","N hits on Segments ME -2/1",7,-0.5,6.5)); 00325 hSnHits.push_back(dbe->book1D("hSnHitsm11a","N hits on Segments ME -1/1a",7,-0.5,6.5)); 00326 hSnHits.push_back(dbe->book1D("hSnHitsm13","N hits on Segments ME -1/3",7,-0.5,6.5)); 00327 hSnHits.push_back(dbe->book1D("hSnHitsm12","N hits on Segments ME -1/2",7,-0.5,6.5)); 00328 hSnHits.push_back(dbe->book1D("hSnHitsm11b","N hits on Segments ME -1/1b",7,-0.5,6.5)); 00329 hSnHits.push_back(dbe->book1D("hSnHitsp11b","N hits on Segments ME +1/1b",7,-0.5,6.5)); 00330 hSnHits.push_back(dbe->book1D("hSnHitsp12","N hits on Segments ME +1/2",7,-0.5,6.5)); 00331 hSnHits.push_back(dbe->book1D("hSnHitsp13","N hits on Segments ME +1/3",7,-0.5,6.5)); 00332 hSnHits.push_back(dbe->book1D("hSnHitsp11a","N hits on Segments ME +1/1a",7,-0.5,6.5)); 00333 hSnHits.push_back(dbe->book1D("hSnHitsp21","N hits on Segments ME +2/1",7,-0.5,6.5)); 00334 hSnHits.push_back(dbe->book1D("hSnHitsp22","N hits on Segments ME +2/2",7,-0.5,6.5)); 00335 hSnHits.push_back(dbe->book1D("hSnHitsp31","N hits on Segments ME +3/1",7,-0.5,6.5)); 00336 hSnHits.push_back(dbe->book1D("hSnHitsp32","N hits on Segments ME +3/2",7,-0.5,6.5)); 00337 hSnHits.push_back(dbe->book1D("hSnHitsp41","N hits on Segments ME +4/1",7,-0.5,6.5)); 00338 hSnHits.push_back(dbe->book1D("hSnHitsp42","N hits on Segments ME +4/2",7,-0.5,6.5)); 00339 00340 hSTheta.push_back(dbe->book1D("hSThetam42","local theta segments in ME -4/2",128,-3.2,3.2)); 00341 hSTheta.push_back(dbe->book1D("hSThetam41","local theta segments in ME -4/1",128,-3.2,3.2)); 00342 hSTheta.push_back(dbe->book1D("hSThetam32","local theta segments in ME -3/2",128,-3.2,3.2)); 00343 hSTheta.push_back(dbe->book1D("hSThetam31","local theta segments in ME -3/1",128,-3.2,3.2)); 00344 hSTheta.push_back(dbe->book1D("hSThetam22","local theta segments in ME -2/2",128,-3.2,3.2)); 00345 hSTheta.push_back(dbe->book1D("hSThetam21","local theta segments in ME -2/1",128,-3.2,3.2)); 00346 hSTheta.push_back(dbe->book1D("hSThetam11a","local theta segments ME -1/1a",128,-3.2,3.2)); 00347 hSTheta.push_back(dbe->book1D("hSThetam13","local theta segments ME -1/3",128,-3.2,3.2)); 00348 hSTheta.push_back(dbe->book1D("hSThetam12","local theta segments ME -1/2",128,-3.2,3.2)); 00349 hSTheta.push_back(dbe->book1D("hSThetam11b","local theta segments ME -1/1b",128,-3.2,3.2)); 00350 hSTheta.push_back(dbe->book1D("hSThetap11b","local theta segments ME +1/1b",128,-3.2,3.2)); 00351 hSTheta.push_back(dbe->book1D("hSThetap12","local theta segments ME +1/2",128,-3.2,3.2)); 00352 hSTheta.push_back(dbe->book1D("hSThetap13","local theta segments ME +1/3",128,-3.2,3.2)); 00353 hSTheta.push_back(dbe->book1D("hSThetap11a","local theta segments ME +1/1a",128,-3.2,3.2)); 00354 hSTheta.push_back(dbe->book1D("hSThetap21","local theta segments in ME +2/1",128,-3.2,3.2)); 00355 hSTheta.push_back(dbe->book1D("hSThetap22","local theta segments in ME +2/2",128,-3.2,3.2)); 00356 hSTheta.push_back(dbe->book1D("hSThetap31","local theta segments in ME +3/1",128,-3.2,3.2)); 00357 hSTheta.push_back(dbe->book1D("hSThetap32","local theta segments in ME +3/2",128,-3.2,3.2)); 00358 hSTheta.push_back(dbe->book1D("hSThetap41","local theta segments in ME +4/1",128,-3.2,3.2)); 00359 hSTheta.push_back(dbe->book1D("hSThetap42","local theta segments in ME +4/2",128,-3.2,3.2)); 00360 00361 //hSGlobal.push_back(dbe->book2D("hSGlobalp1","segment global X,Y station +1",400,-800.,800.,400,-800.,800.)); 00362 //hSGlobal.push_back(dbe->book2D("hSGlobalp2","segment global X,Y station +2",400,-800.,800.,400,-800.,800.)); 00363 //hSGlobal.push_back(dbe->book2D("hSGlobalp3","segment global X,Y station +3",400,-800.,800.,400,-800.,800.)); 00364 //hSGlobal.push_back(dbe->book2D("hSGlobalp4","segment global X,Y station +4",400,-800.,800.,400,-800.,800.)); 00365 //hSGlobal.push_back(dbe->book2D("hSGlobalm1","segment global X,Y station -1",400,-800.,800.,400,-800.,800.)); 00366 //hSGlobal.push_back(dbe->book2D("hSGlobalm2","segment global X,Y station -2",400,-800.,800.,400,-800.,800.)); 00367 //hSGlobal.push_back(dbe->book2D("hSGlobalm3","segment global X,Y station -3",400,-800.,800.,400,-800.,800.)); 00368 //hSGlobal.push_back(dbe->book2D("hSGlobalm4","segment global X,Y station -4",400,-800.,800.,400,-800.,800.)); 00369 00370 hSnhitsAll = dbe->book1D("hSnhits","N hits on Segments",7,-0.5,6.5); 00371 hSChiSqProb = dbe->book1D("hSChiSqProb","segments chi-squared probability",100,0.,1.); 00372 hSGlobalTheta = dbe->book1D("hSGlobalTheta","segment global theta",64,0,1.6); 00373 hSGlobalPhi = dbe->book1D("hSGlobalPhi", "segment global phi", 128,-3.2,3.2); 00374 hSnSegments = dbe->book1D("hSnSegments","number of segments per event",11,-0.5,10.5); 00375 00376 hSResid.push_back(dbe->book1D("hSResidm42","Fitted Position on Strip - Reconstructed for Layer 3 (ME -4/2)",100,-0.5,0.5)); 00377 hSResid.push_back(dbe->book1D("hSResidm41","Fitted Position on Strip - Reconstructed for Layer 3 (ME -4/1)",100,-0.5,0.5)); 00378 hSResid.push_back(dbe->book1D("hSResidm32","Fitted Position on Strip - Reconstructed for Layer 3 (ME -3/2)",100,-0.5,0.5)); 00379 hSResid.push_back(dbe->book1D("hSResidm31","Fitted Position on Strip - Reconstructed for Layer 3 (ME -3/1)",100,-0.5,0.5)); 00380 hSResid.push_back(dbe->book1D("hSResidm22","Fitted Position on Strip - Reconstructed for Layer 3 (ME -2/2)",100,-0.5,0.5)); 00381 hSResid.push_back(dbe->book1D("hSResidm21","Fitted Position on Strip - Reconstructed for Layer 3 (ME -2/1)",100,-0.5,0.5)); 00382 hSResid.push_back(dbe->book1D("hSResidm11a","Fitted Position on Strip - Reconstructed for Layer 3 (ME -1/1a)",100,-0.5,0.5)); 00383 hSResid.push_back(dbe->book1D("hSResidm13","Fitted Position on Strip - Reconstructed for Layer 3 (ME -1/3)",100,-0.5,0.5)); 00384 hSResid.push_back(dbe->book1D("hSResidm12","Fitted Position on Strip - Reconstructed for Layer 3 (ME -1/2)",100,-0.5,0.5)); 00385 hSResid.push_back(dbe->book1D("hSResidm11b","Fitted Position on Strip - Reconstructed for Layer 3 (ME -1/1b)",100,-0.5,0.5)); 00386 hSResid.push_back(dbe->book1D("hSResidp11b","Fitted Position on Strip - Reconstructed for Layer 3 (ME +1/1b)",100,-0.5,0.5)); 00387 hSResid.push_back(dbe->book1D("hSResidp12","Fitted Position on Strip - Reconstructed for Layer 3 (ME +1/2)",100,-0.5,0.5)); 00388 hSResid.push_back(dbe->book1D("hSResidp13","Fitted Position on Strip - Reconstructed for Layer 3 (ME +1/3)",100,-0.5,0.5)); 00389 hSResid.push_back(dbe->book1D("hSResidp11a","Fitted Position on Strip - Reconstructed for Layer 3 (ME +1/1a)",100,-0.5,0.5)); 00390 hSResid.push_back(dbe->book1D("hSResidp21","Fitted Position on Strip - Reconstructed for Layer 3 (ME +2/1)",100,-0.5,0.5)); 00391 hSResid.push_back(dbe->book1D("hSResidp22","Fitted Position on Strip - Reconstructed for Layer 3 (ME +2/2)",100,-0.5,0.5)); 00392 hSResid.push_back(dbe->book1D("hSResidp31","Fitted Position on Strip - Reconstructed for Layer 3 (ME +3/1)",100,-0.5,0.5)); 00393 hSResid.push_back(dbe->book1D("hSResidp32","Fitted Position on Strip - Reconstructed for Layer 3 (ME +3/2)",100,-0.5,0.5)); 00394 hSResid.push_back(dbe->book1D("hSResidp41","Fitted Position on Strip - Reconstructed for Layer 3 (ME +4/1)",100,-0.5,0.5)); 00395 hSResid.push_back(dbe->book1D("hSResidp42","Fitted Position on Strip - Reconstructed for Layer 3 (ME +4/2)",100,-0.5,0.5)); 00396 00397 //occupancy plots 00398 dbe->setCurrentFolder("CSC/CSCOfflineMonitor/LocalRecoOccupancy"); 00399 hOWires = dbe->book2D("hOWires","Wire Digi Occupancy",36,0.5,36.5,20,0.5,20.5); 00400 hOStrips = dbe->book2D("hOStrips","Strip Digi Occupancy",36,0.5,36.5,20,0.5,20.5); 00401 hORecHits = dbe->book2D("hORecHits","RecHit Occupancy",36,0.5,36.5,20,0.5,20.5); 00402 hOSegments = dbe->book2D("hOSegments","Segment Occupancy",36,0.5,36.5,20,0.5,20.5); 00403 00404 00405 00406 }
Reimplemented from edm::EDAnalyzer.
Definition at line 415 of file CSCOfflineMonitor.cc.
References dbe, edm::ParameterSet::getParameter(), iggi_31X_cfg::outputFileName, param, and DQMStore::save().
00415 { 00416 bool saveHistos = param.getParameter<bool>("saveHistos"); 00417 string outputFileName = param.getParameter<string>("outputFileName"); 00418 if(saveHistos){ 00419 dbe->save(outputFileName); 00420 } 00421 }
| float CSCOfflineMonitor::fitX | ( | HepMatrix | sp, | |
| HepMatrix | ep | |||
| ) | [private] |
Definition at line 833 of file CSCOfflineMonitor.cc.
References HLT_VtxMuL3::errors, i, and slope.
Referenced by analyze().
00833 { 00834 00835 float S = 0; 00836 float Sx = 0; 00837 float Sy = 0; 00838 float Sxx = 0; 00839 float Sxy = 0; 00840 float sigma2 = 0; 00841 00842 for (int i=1;i<7;i++){ 00843 if (i != 3){ 00844 sigma2 = errors(i,1)*errors(i,1); 00845 S = S + (1/sigma2); 00846 Sy = Sy + (points(i,1)/sigma2); 00847 Sx = Sx + ((i)/sigma2); 00848 Sxx = Sxx + (i*i)/sigma2; 00849 Sxy = Sxy + (((i)*points(i,1))/sigma2); 00850 } 00851 } 00852 00853 float delta = S*Sxx - Sx*Sx; 00854 float intercept = (Sxx*Sy - Sx*Sxy)/delta; 00855 float slope = (S*Sxy - Sx*Sy)/delta; 00856 00857 float chi = 0; 00858 float chi2 = 0; 00859 00860 // calculate chi2 (not currently used) 00861 for (int i=1;i<7;i++){ 00862 chi = (points(i,1) - intercept - slope*i)/(errors(i,1)); 00863 chi2 = chi2 + chi*chi; 00864 } 00865 00866 return (intercept + slope*3); 00867 00868 }
| float CSCOfflineMonitor::getSignal | ( | const CSCStripDigiCollection & | stripdigis, | |
| CSCDetId | idRH, | |||
| int | centerStrip | |||
| ) | [private] |
Definition at line 924 of file CSCOfflineMonitor.cc.
References if(), and prof2calltree::last.
Referenced by analyze().
00925 { 00926 00927 float SigADC[5]; 00928 float TotalADC = 0; 00929 SigADC[0] = 0; 00930 SigADC[1] = 0; 00931 SigADC[2] = 0; 00932 SigADC[3] = 0; 00933 SigADC[4] = 0; 00934 00935 00936 // Loop over strip digis 00937 CSCStripDigiCollection::DigiRangeIterator sIt; 00938 00939 for (sIt = stripdigis.begin(); sIt != stripdigis.end(); sIt++){ 00940 CSCDetId id = (CSCDetId)(*sIt).first; 00941 if (id == idCS){ 00942 00943 // First, find the Signal-Pedestal for center strip 00944 vector<CSCStripDigi>::const_iterator digiItr = (*sIt).second.first; 00945 vector<CSCStripDigi>::const_iterator last = (*sIt).second.second; 00946 for ( ; digiItr != last; ++digiItr ) { 00947 int thisStrip = digiItr->getStrip(); 00948 if (thisStrip == (centerStrip)){ 00949 std::vector<int> myADCVals = digiItr->getADCCounts(); 00950 float thisPedestal = 0.5*(float)(myADCVals[0]+myADCVals[1]); 00951 float thisSignal = (myADCVals[2]+myADCVals[3]+myADCVals[4]+myADCVals[5]+myADCVals[6]+myADCVals[7]); 00952 SigADC[0] = thisSignal - 6*thisPedestal; 00953 } 00954 // Now,find the Signal-Pedestal for neighbouring 4 strips 00955 if (thisStrip == (centerStrip+1)){ 00956 std::vector<int> myADCVals = digiItr->getADCCounts(); 00957 float thisPedestal = 0.5*(float)(myADCVals[0]+myADCVals[1]); 00958 float thisSignal = (myADCVals[2]+myADCVals[3]+myADCVals[4]+myADCVals[5]+myADCVals[6]+myADCVals[7]); 00959 SigADC[1] = thisSignal - 6*thisPedestal; 00960 } 00961 if (thisStrip == (centerStrip+2)){ 00962 std::vector<int> myADCVals = digiItr->getADCCounts(); 00963 float thisPedestal = 0.5*(float)(myADCVals[0]+myADCVals[1]); 00964 float thisSignal = (myADCVals[2]+myADCVals[3]+myADCVals[4]+myADCVals[5]+myADCVals[6]+myADCVals[7]); 00965 SigADC[2] = thisSignal - 6*thisPedestal; 00966 } 00967 if (thisStrip == (centerStrip-1)){ 00968 std::vector<int> myADCVals = digiItr->getADCCounts(); 00969 float thisPedestal = 0.5*(float)(myADCVals[0]+myADCVals[1]); 00970 float thisSignal = (myADCVals[2]+myADCVals[3]+myADCVals[4]+myADCVals[5]+myADCVals[6]+myADCVals[7]); 00971 SigADC[3] = thisSignal - 6*thisPedestal; 00972 } 00973 if (thisStrip == (centerStrip-2)){ 00974 std::vector<int> myADCVals = digiItr->getADCCounts(); 00975 float thisPedestal = 0.5*(float)(myADCVals[0]+myADCVals[1]); 00976 float thisSignal = (myADCVals[2]+myADCVals[3]+myADCVals[4]+myADCVals[5]+myADCVals[6]+myADCVals[7]); 00977 SigADC[4] = thisSignal - 6*thisPedestal; 00978 } 00979 } 00980 TotalADC = 0.2*(SigADC[0]+SigADC[1]+SigADC[2]+SigADC[3]+SigADC[4]); 00981 } 00982 } 00983 return TotalADC; 00984 }
| float CSCOfflineMonitor::getTiming | ( | const CSCStripDigiCollection & | stripdigis, | |
| CSCDetId | idRH, | |||
| int | centerStrip | |||
| ) | [private] |
Definition at line 876 of file CSCOfflineMonitor.cc.
References diff, i, if(), CSCDetId::ring(), and CSCDetId::station().
00876 { 00877 00878 float ADC[8]; 00879 for (int i = 0; i < 8; i++){ 00880 ADC[i] = 0; 00881 } 00882 float timing = -1; 00883 00884 if (idRH.station() == 1 && idRH.ring() == 4){ 00885 while(centerStrip> 16) centerStrip -= 16; 00886 } 00887 00888 // Loop over strip digis responsible for this recHit and sum charge 00889 CSCStripDigiCollection::DigiRangeIterator gTstripIter; 00890 for (gTstripIter = stripdigis.begin(); gTstripIter != stripdigis.end(); gTstripIter++){ 00891 CSCDetId id = (CSCDetId)(*gTstripIter).first; 00892 if (id == idRH){ 00893 vector<CSCStripDigi>::const_iterator STiter = (*gTstripIter).second.first; 00894 vector<CSCStripDigi>::const_iterator lastS = (*gTstripIter).second.second; 00895 for ( ; STiter != lastS; ++STiter ) { 00896 int thisStrip = STiter->getStrip(); 00897 if (thisStrip == (centerStrip)){ 00898 float diff = 0; 00899 vector<int> myADCVals = STiter->getADCCounts(); 00900 float thisPedestal = 0.5*(float)(myADCVals[0]+myADCVals[1]); 00901 for (unsigned int iCount = 0; iCount < myADCVals.size(); iCount++) { 00902 diff = (float)myADCVals[iCount]-thisPedestal; 00903 ADC[iCount] = diff; 00904 } 00905 } 00906 } 00907 } 00908 } 00909 00910 timing = (ADC[2]*2 + ADC[3]*3 + ADC[4]*4 + ADC[5]*5 + ADC[6]*6)/(ADC[2] + ADC[3] + ADC[4] + ADC[5] + ADC[6]); 00911 00912 return timing; 00913 00914 }
Definition at line 986 of file CSCOfflineMonitor.cc.
References GeomDetEnumerators::endcap, and index.
Referenced by analyze().
00986 { 00987 00988 // linearlized index bases on endcap, station, and ring 00989 int index = 0; 00990 if (id.station() == 1) index = id.ring(); 00991 else index = id.station()*2 + id.ring(); 00992 if (id.endcap() == 1) index = index + 9; 00993 if (id.endcap() == 2) index = 10 - index; 00994 return index; 00995 00996 }
Definition at line 83 of file CSCOfflineMonitor.h.
Referenced by analyze(), and CSCOfflineMonitor().
edm::InputTag CSCOfflineMonitor::cscSegTag_ [private] |
Definition at line 84 of file CSCOfflineMonitor.h.
Referenced by analyze(), and CSCOfflineMonitor().
DQMStore* CSCOfflineMonitor::dbe [private] |
MonitorElement* CSCOfflineMonitor::hORecHits [private] |
MonitorElement* CSCOfflineMonitor::hOSegments [private] |
MonitorElement* CSCOfflineMonitor::hOStrips [private] |
MonitorElement* CSCOfflineMonitor::hOWires [private] |
MonitorElement* CSCOfflineMonitor::hRHCodeBroad [private] |
std::vector<MonitorElement*> CSCOfflineMonitor::hRHCodeNarrow [private] |
std::vector<MonitorElement*> CSCOfflineMonitor::hRHLayer [private] |
MonitorElement* CSCOfflineMonitor::hRHnrechits [private] |
std::vector<MonitorElement*> CSCOfflineMonitor::hRHRatioQ [private] |
std::vector<MonitorElement*> CSCOfflineMonitor::hRHResid [private] |
Definition at line 126 of file CSCOfflineMonitor.h.
std::vector<MonitorElement*> CSCOfflineMonitor::hRHSumQ [private] |
std::vector<MonitorElement*> CSCOfflineMonitor::hRHTiming [private] |
std::vector<MonitorElement*> CSCOfflineMonitor::hRHX [private] |
std::vector<MonitorElement*> CSCOfflineMonitor::hRHY [private] |
MonitorElement* CSCOfflineMonitor::hSChiSqProb [private] |
MonitorElement* CSCOfflineMonitor::hSCodeBroad [private] |
std::vector<MonitorElement*> CSCOfflineMonitor::hSCodeNarrow [private] |
MonitorElement* CSCOfflineMonitor::hSGlobalPhi [private] |
MonitorElement* CSCOfflineMonitor::hSGlobalTheta [private] |
std::vector<MonitorElement*> CSCOfflineMonitor::hSnHits [private] |
MonitorElement* CSCOfflineMonitor::hSnhitsAll [private] |
MonitorElement* CSCOfflineMonitor::hSnSegments [private] |
std::vector<MonitorElement*> CSCOfflineMonitor::hSResid [private] |
std::vector<MonitorElement*> CSCOfflineMonitor::hSTheta [private] |
MonitorElement* CSCOfflineMonitor::hStripAll [private] |
MonitorElement* CSCOfflineMonitor::hStripCodeBroad [private] |
std::vector<MonitorElement*> CSCOfflineMonitor::hStripCodeNarrow [private] |
std::vector<MonitorElement*> CSCOfflineMonitor::hStripLayer [private] |
MonitorElement* CSCOfflineMonitor::hStripNFired [private] |
std::vector<MonitorElement*> CSCOfflineMonitor::hStripPed [private] |
MonitorElement* CSCOfflineMonitor::hStripPedAll [private] |
std::vector<MonitorElement*> CSCOfflineMonitor::hStripStrip [private] |
MonitorElement* CSCOfflineMonitor::hWireAll [private] |
MonitorElement* CSCOfflineMonitor::hWireCodeBroad [private] |
std::vector<MonitorElement*> CSCOfflineMonitor::hWireCodeNarrow [private] |
std::vector<MonitorElement*> CSCOfflineMonitor::hWireLayer [private] |
MonitorElement* CSCOfflineMonitor::hWireTBinAll [private] |
std::vector<MonitorElement*> CSCOfflineMonitor::hWireWire [private] |
edm::ParameterSet CSCOfflineMonitor::param [private] |
Definition at line 81 of file CSCOfflineMonitor.h.
Referenced by analyze(), and CSCOfflineMonitor().
edm::InputTag CSCOfflineMonitor::wireDigiTag_ [private] |
Definition at line 82 of file CSCOfflineMonitor.h.
Referenced by analyze(), and CSCOfflineMonitor().
1.5.4