---

CSCTFValidator.cc

Go to the documentation of this file.

#include <iostream>
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "EventFilter/CSCTFRawToDigi/interface/CSCTFValidator.h"
#include "DataFormats/CSCDigi/interface/CSCCorrelatedLCTDigi.h"
#include "DataFormats/CSCDigi/interface/CSCCorrelatedLCTDigiCollection.h"
#include "DataFormats/FEDRawData/interface/FEDRawData.h"
#include "DataFormats/FEDRawData/interface/FEDNumbering.h"
#include "DataFormats/FEDRawData/interface/FEDRawDataCollection.h"
#include "DataFormats/MuonDetId/interface/CSCDetId.h"
#include "EventFilter/CSCTFRawToDigi/src/CSCTFTBEventData.h"
#include "CondFormats/CSCObjects/interface/CSCTriggerMappingFromFile.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "DataFormats/L1CSCTrackFinder/interface/L1CSCTrackCollection.h"


CSCTFValidator::CSCTFValidator(edm::ParameterSet const& conf) {

  // If your module takes parameters, here is where you would define
  // their names and types, and access them to initialize internal
  // variables. Example as follows:
  //
  testBeam = conf.getUntrackedParameter<bool>("TestBeam",true);
  std::string mapPath = "/" + conf.getUntrackedParameter<std::string>("MappingFile","");

  if(testBeam)
    {
      TBFEDid = conf.getUntrackedParameter<int>("TBFedId");
      TBendcap = conf.getUntrackedParameter<int>("TBEndcap");
      TBsector = conf.getUntrackedParameter<int>("TBSector");
    }
  else
    {
      TBFEDid = 0;
      TBsector = 0;
      TBendcap = 0;
    }

  TFMapping = new CSCTriggerMappingFromFile(getenv("CMSSW_BASE") + mapPath);

  eventNumber = 0;
}

void CSCTFValidator::analyze(edm::Event const& e, edm::EventSetup const& iSetup) {

  // These declarations create handles to the types of records that you want
  // to retrieve from event "e".
  //
  edm::Handle<CSCCorrelatedLCTDigiCollection> corrlcts;
  edm::Handle<L1CSCTrackCollection> tracks;
  edm::Handle<FEDRawDataCollection> rawdata;

  // Pass the handle to the method "getByType", which is used to retrieve
  // one and only one instance of the type in question out of event "e". If
  // zero or more than one instance exists in the event an exception is thrown.
  //

  e.getByLabel("csctfunpacker","MuonCSCTFCorrelatedLCTDigi",corrlcts);
  e.getByLabel("csctfunpacker","MuonL1CSCTrackCollection",tracks);
  //e.getByLabel("DaqSource",rawdata);

  // read digi collections and print digis
  // compare to raw data
  CSCTFTBEventData *tbdata = NULL;

  for(int fedid = FEDNumbering::getCSCTFFEDIds().first;
      fedid <= ((testBeam) ? (FEDNumbering::getCSCTFFEDIds().first) : (FEDNumbering::getCSCTFFEDIds().second));
      ++fedid)
    {
/*      const FEDRawData& fedData = rawdata->FEDData(fedid);
      if(fedData.size())
        {
      if(testBeam){
            CSCTFTBEventData *tbdata = new CSCTFTBEventData(reinterpret_cast<unsigned short*>(fedData.data()));

                eventNumber = tbdata->eventHeader().getLvl1num();

            CSCTFTBFrontBlock aFB;
            CSCTFTBSPBlock aSPB;
            CSCTFTBSPData aSPD;

            for(int BX = 1; BX <= tbdata->eventHeader().numBX(); ++BX)
              {
                if(testBeam) aFB = tbdata->frontDatum(BX);
                else edm::LogInfo("CSCTFValidator|analyze") << "Not implemented yet\n";
                for(int FPGA = 1; FPGA <= 5; ++FPGA)
                  {
                  for(int MPClink = 1; MPClink <= 3; ++MPClink)
                    {
                          int subsector = 0;
                          int station = 0;

                          if(FPGA == 1) subsector = 1;
                          if(FPGA == 2) subsector = 1;
                          station = (((FPGA - 1)==0) ? 1 : FPGA - 1);

                          int cscid = aFB.frontData(FPGA,MPClink).CSCIDPacked();
                          if(cscid)
                            {
                              try
                                {
                                  CSCDetId id = TFMapping->detId(TBendcap,station,TBsector,subsector,cscid,3);
                                  int match = 0;

                                  CSCCorrelatedLCTDigiCollection::Range range = corrlcts->get(id);
                                  for(CSCCorrelatedLCTDigiCollection::const_iterator j = range.first;
                                      j != range.second; j++)
                                    if(j->getTrknmb() == aFB.frontDigiData(FPGA,MPClink).getTrknmb() &&
                                       j->isValid() == aFB.frontDigiData(FPGA,MPClink).isValid() &&
                                       j->getQuality() == aFB.frontDigiData(FPGA,MPClink).getQuality() &&
                                       j->getKeyWG() == aFB.frontDigiData(FPGA,MPClink).getKeyWG() &&
                                       j->getStrip() == aFB.frontDigiData(FPGA,MPClink).getStrip() &&
                                       j->getCLCTPattern() == aFB.frontDigiData(FPGA,MPClink).getCLCTPattern() &&
                                       j->getStripType() == aFB.frontDigiData(FPGA,MPClink).getStripType() &&
                                       j->getBend() == aFB.frontDigiData(FPGA,MPClink).getBend() &&
                                       j->getBX() == aFB.frontDigiData(FPGA,MPClink).getBX()) ++match;


                                  if(match < 1)
                                    {
                                      edm::LogError("CSCTFBValidator|analyze") << "NO DIGI IN DET\n";
                                      assert(1==0);
                                    }
                                }
                              catch (cms::Exception &e)
                                {
                                  edm::LogInfo("CSCTFValidator|analyze") << e.what() << "\ndigi not processed.";
                                }
                            }
                  }
              }
            }
          }
      delete tbdata;

        } else { // Not a Test Beam
*/              for(unsigned int FPGA=0; FPGA<5; FPGA++){
                        int station = ( FPGA ? FPGA : 1 );
                        int subsector = ( FPGA>1 ? 0 : FPGA+1 );

                        for(int endcap=1; endcap<=2; endcap++)
                                for(int sector=1; sector<=6; sector++)
                                        for(int cscid=1; cscid<10; cscid++){

                                                try {
                                                        CSCDetId id = TFMapping->detId(endcap,station,sector,subsector,cscid,0);
                                                        CSCCorrelatedLCTDigiCollection::Range range = corrlcts->get(id);
                                                        for(CSCCorrelatedLCTDigiCollection::const_iterator j=range.first; j!=range.second; j++)
                                                                std::cout<<"Trk: "<<j->getTrknmb()<<" valid: "<<j->isValid()<<" quality: "<<j->getQuality()
                                                                        <<" WG: "<<j->getKeyWG()<<" strip: "<<j->getStrip()<<" CLCT: "<<j->getCLCTPattern()
                                                                        <<" type: "<<j->getStripType()<<" bend: "<<j->getBend()<<" BX: "<<j->getBX()<<std::endl;

                                                } catch (cms::Exception &e) {}
                                        }

                        for(L1CSCTrackCollection::const_iterator trk=tracks->begin(); trk<tracks->end(); trk++)
                                std::cout<<"Trk: "<<trk->first.endcap()<<" endcap: "<<trk->first.endcap()<<" sector: "<<trk->first.sector()<<" BX: "<<trk->first.BX()<<std::endl;
                }
/*      }

*/

        }

  LogDebug ("CSCTFValidator::analyze")  << "end of event number " << eventNumber;
}

---