Inheritance diagram for CSCTFAnalyzer:

Public Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override

	CSCTFAnalyzer (const edm::ParameterSet &conf)

	~CSCTFAnalyzer (void) override

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
edm::EDGetTokenT< CSCCorrelatedLCTDigiCollection >	CSCCDC_Tok

edm::EDGetTokenT< CSCTriggerContainer< csctf::TrackStub > >	CSCTC_Tok

int	dtPhi [12][2]

TFile *	file

edm::EDGetTokenT< L1CSCStatusDigiCollection >	L1CSCS_Tok

edm::EDGetTokenT< L1CSCTrackCollection >	L1CST_Tok

edm::InputTag	lctProducer

edm::InputTag	mbProducer

edm::InputTag	statusProducer

edm::InputTag	trackProducer

TTree *	tree

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 18 of file CSCTFAnalyzer.cc.

Constructor & Destructor Documentation

◆ CSCTFAnalyzer()

CSCTFAnalyzer::CSCTFAnalyzer ( const edm::ParameterSet & conf )

explicit

Definition at line 49 of file CSCTFAnalyzer.cc.

References CSCCDC_Tok, CSCTC_Tok, dtPhi, file, edm::ParameterSet::getUntrackedParameter(), HLT_2022v11_cff::InputTag, edm::InputTag::instance(), L1CSCS_Tok, L1CST_Tok, edm::InputTag::label(), lctProducer, mbProducer, statusProducer, and trackProducer.

                                                         : edm::one::EDAnalyzer<>() {
   mbProducer = conf.getUntrackedParameter<edm::InputTag>("mbProducer", edm::InputTag("csctfunpacker"));
   lctProducer = conf.getUntrackedParameter<edm::InputTag>("lctProducer", edm::InputTag("csctfunpacker"));
   trackProducer = conf.getUntrackedParameter<edm::InputTag>("trackProducer", edm::InputTag("csctfunpacker"));
   statusProducer = conf.getUntrackedParameter<edm::InputTag>("statusProducer", edm::InputTag("csctfunpacker"));
   file = new TFile("ewq.root", "RECREATE");
   tree = new TTree("dy", "QWE");
   tree->Branch("dtPhi_1_plus", &dtPhi[0][0], "dtPhi_1_plus/I");
   tree->Branch("dtPhi_2_plus", &dtPhi[1][0], "dtPhi_2_plus/I");
   tree->Branch("dtPhi_3_plus", &dtPhi[2][0], "dtPhi_3_plus/I");
   tree->Branch("dtPhi_4_plus", &dtPhi[3][0], "dtPhi_4_plus/I");
   tree->Branch("dtPhi_5_plus", &dtPhi[4][0], "dtPhi_5_plus/I");
   tree->Branch("dtPhi_6_plus", &dtPhi[5][0], "dtPhi_6_plus/I");
   tree->Branch("dtPhi_7_plus", &dtPhi[6][0], "dtPhi_7_plus/I");
   tree->Branch("dtPhi_8_plus", &dtPhi[7][0], "dtPhi_8_plus/I");
   tree->Branch("dtPhi_9_plus", &dtPhi[8][0], "dtPhi_9_plus/I");
   tree->Branch("dtPhi_10_plus", &dtPhi[9][0], "dtPhi_10_plus/I");
   tree->Branch("dtPhi_11_plus", &dtPhi[10][0], "dtPhi_11_plus/I");
   tree->Branch("dtPhi_12_plus", &dtPhi[11][0], "dtPhi_12_plus/I");
   tree->Branch("dtPhi_1_minus", &dtPhi[0][1], "dtPhi_1_minus/I");
   tree->Branch("dtPhi_2_minus", &dtPhi[1][1], "dtPhi_2_minus/I");
   tree->Branch("dtPhi_3_minus", &dtPhi[2][1], "dtPhi_3_minus/I");
   tree->Branch("dtPhi_4_minus", &dtPhi[3][1], "dtPhi_4_minus/I");
   tree->Branch("dtPhi_5_minus", &dtPhi[4][1], "dtPhi_5_minus/I");
   tree->Branch("dtPhi_6_minus", &dtPhi[5][1], "dtPhi_6_minus/I");
   tree->Branch("dtPhi_7_minus", &dtPhi[6][1], "dtPhi_7_minus/I");
   tree->Branch("dtPhi_8_minus", &dtPhi[7][1], "dtPhi_8_minus/I");
   tree->Branch("dtPhi_9_minus", &dtPhi[8][1], "dtPhi_9_minus/I");
   tree->Branch("dtPhi_10_minus", &dtPhi[9][1], "dtPhi_10_minus/I");
   tree->Branch("dtPhi_11_minus", &dtPhi[10][1], "dtPhi_11_minus/I");
   tree->Branch("dtPhi_12_minus", &dtPhi[11][1], "dtPhi_12_minus/I");
 
   L1CSCS_Tok = consumes<L1CSCStatusDigiCollection>(edm::InputTag(statusProducer.label(), statusProducer.instance()));
   CSCTC_Tok =
       consumes<CSCTriggerContainer<csctf::TrackStub> >(edm::InputTag(mbProducer.label(), mbProducer.instance()));
   CSCCDC_Tok = consumes<CSCCorrelatedLCTDigiCollection>(edm::InputTag(lctProducer.label(), lctProducer.instance()));
   L1CST_Tok = consumes<L1CSCTrackCollection>(edm::InputTag(trackProducer.label(), trackProducer.instance()));
 }

◆ ~CSCTFAnalyzer()

CSCTFAnalyzer::~CSCTFAnalyzer ( void )

inlineoverride

Definition at line 34 of file CSCTFAnalyzer.cc.

References file.

                                 {
     file->cd();
     tree->Write();
     file->Close();
   }

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 88 of file CSCTFAnalyzer.cc.

References gather_cfg::cout, CSCCDC_Tok, CSCTC_Tok, dtPhi, MillePedeFileConverter_cfg::e, edm::HandleBase::isValid(), L1CSCS_Tok, L1CST_Tok, edm::InputTag::label(), lctProducer, mbProducer, edm::Handle< T >::product(), relativeConstraints::station, mps_update::status, statusProducer, trackProducer, tracks, and CSCTriggerNumbering::triggerSubSectorFromLabels().

                                                                      {
   /*    edm::Handle<FEDRawDataCollection> rawdata;
     e.getByToken("source","",rawdata);
 
     const FEDRawData& fedData = rawdata->FEDData(750);
     if( fedData.size()==0 ) return;
     unsigned short *data = (unsigned short *)fedData.data();
     unsigned int    size = fedData.size()/2;
     std::cout<<"New event:"<<std::endl;
     for(unsigned i=0; i<size/4; i++)
         std::cout<<std::hex<<" "<<std::setw(6)<<data[i*4+0]<<" "<<std::setw(6)<<data[i*4+1]<<" "<<std::setw(6)<<data[i*4+2]<<" "<<std::setw(6)<<data[i*4+3]<<std::dec<<std::endl;
     std::cout<<"End of event"<<std::endl;
     return;
 */
   if (statusProducer.label() != "null") {
     edm::Handle<L1CSCStatusDigiCollection> status;
     e.getByToken(L1CSCS_Tok, status);
     if (status.isValid()) {
       edm::LogInfo("CSCTFAnalyzer") << "  Unpacking Errors: " << status->first;
       for (std::vector<L1CSCSPStatusDigi>::const_iterator stat = status->second.begin(); stat != status->second.end();
            stat++) {
         //edm::LogInfo("CSCTFAnalyzer") << "   Status: SP in slot "<<stat->slot()<<"  FMM: "<<stat->FMM()<<" SE: 0x"<<std::hex<<stat->SEs()<<" VP: 0x"<<stat->VPs()<<std::dec;
       }
     } else
       edm::LogInfo("CSCTFAnalyzer") << "  No valid L1CSCStatusDigiCollection products found";
   }
 
   if (mbProducer.label() != "null") {
     bzero(dtPhi, sizeof(dtPhi));
     edm::Handle<CSCTriggerContainer<csctf::TrackStub> > dtStubs;
     e.getByToken(CSCTC_Tok, dtStubs);
     if (dtStubs.isValid()) {
       std::vector<csctf::TrackStub> vstubs = dtStubs->get();
       std::cout << "DT size=" << vstubs.end() - vstubs.begin() << std::endl;
       for (std::vector<csctf::TrackStub>::const_iterator stub = vstubs.begin(); stub != vstubs.end(); stub++) {
         int dtSector = (stub->sector() - 1) * 2 + stub->subsector() - 1;
         int dtEndcap = stub->endcap() - 1;
         if (dtSector >= 0 && dtSector < 12 && dtEndcap >= 0 && dtEndcap < 2) {
           dtPhi[dtSector][dtEndcap] = stub->phiPacked();
         } else {
           edm::LogInfo("CSCTFAnalyzer: DT digi are out of range: ")
               << " dtSector=" << dtSector << " dtEndcap=" << dtEndcap;
         }
         edm::LogInfo("CSCTFAnalyzer") << "   DT data: tbin=" << stub->BX() << " CSC sector=" << stub->sector()
                                       << " CSC subsector=" << stub->subsector() << " station=" << stub->station()
                                       << " endcap=" << stub->endcap() << " phi=" << stub->phiPacked()
                                       << " phiBend=" << stub->getBend() << " quality=" << stub->getQuality()
                                       << " mb_bxn=" << stub->cscid();
       }
     } else
       edm::LogInfo("CSCTFAnalyzer") << "  No valid CSCTriggerContainer<csctf::TrackStub> products found";
     tree->Fill();
   }
 
   if (lctProducer.label() != "null") {
     edm::Handle<CSCCorrelatedLCTDigiCollection> corrlcts;
     e.getByToken(CSCCDC_Tok, corrlcts);
     if (corrlcts.isValid()) {
       for (CSCCorrelatedLCTDigiCollection::DigiRangeIterator csc = corrlcts.product()->begin();
            csc != corrlcts.product()->end();
            csc++) {
         int lctId = 0;
         CSCCorrelatedLCTDigiCollection::Range range1 = corrlcts.product()->get((*csc).first);
         for (CSCCorrelatedLCTDigiCollection::const_iterator lct = range1.first; lct != range1.second; lct++, lctId++) {
           int station = (*csc).first.station() - 1;
           int cscId = (*csc).first.triggerCscId() - 1;
           int sector = (*csc).first.triggerSector() - 1;  // + ( (*csc).first.endcap()==1 ? 0 : 6 );
           int subSector = CSCTriggerNumbering::triggerSubSectorFromLabels((*csc).first);
           int tbin = lct->getBX();
           int fpga = (subSector ? subSector - 1 : station + 1);
           // If Det Id is within range
           if (sector < 0 || sector > 11 || station < 0 || station > 3 || cscId < 0 || cscId > 8 || lctId < 0 ||
               lctId > 1) {
             edm::LogInfo("CSCTFAnalyzer: CSC digi are out of range: ");
             continue;
           }
           edm::LogInfo("CSCTFAnalyzer") << "   Front data   endcap: " << (*csc).first.endcap()
                                         << "  station: " << (station + 1) << "  sector: " << (sector + 1)
                                         << "  subSector: " << subSector << "  tbin: " << tbin
                                         << "  cscId: " << (cscId + 1) << "  fpga: " << (fpga + 1) << " "
                                         << "LCT(vp=" << lct->isValid() << ",qual=" << lct->getQuality()
                                         << ",wg=" << lct->getKeyWG() << ",strip=" << lct->getStrip()
                                         << ",link=" << lct->getMPCLink() << ")";
         }
       }
     } else
       edm::LogInfo("CSCTFAnalyzer") << "  No valid CSCCorrelatedLCTDigiCollection products found";
   }
 
   if (trackProducer.label() != "null") {
     edm::Handle<L1CSCTrackCollection> tracks;
     e.getByToken(L1CST_Tok, tracks);
     if (tracks.isValid()) {
       int nTrk = 0;
       for (L1CSCTrackCollection::const_iterator trk = tracks->begin(); trk < tracks->end(); trk++, nTrk++) {
         int sector = 6 * (trk->first.endcap() - 1) + trk->first.sector() - 1;
         int tbin = trk->first.BX();
         edm::LogInfo("CSCTFAnalyzer") << "   Track sector: " << (sector + 1) << "  tbin: " << tbin << " "
                                       << "TRK(mode=" << ((trk->first.ptLUTAddress() >> 16) & 0xF)
                                       << ",eta=" << trk->first.eta_packed() << ",phi=" << trk->first.localPhi()
                                       << ") IDs:"
                                       << " me1D=" << trk->first.me1ID() << " t1=" << trk->first.me1Tbin()
                                       << " me2D=" << trk->first.me2ID() << " t2=" << trk->first.me2Tbin()
                                       << " me3D=" << trk->first.me3ID() << " t3=" << trk->first.me3Tbin()
                                       << " me4D=" << trk->first.me4ID() << " t4=" << trk->first.me4Tbin()
                                       << " mb1D=" << trk->first.mb1ID() << " tb=" << trk->first.mb1Tbin();
 
         for (CSCCorrelatedLCTDigiCollection::DigiRangeIterator csc = trk->second.begin(); csc != trk->second.end();
              csc++) {
           int lctId = 0;
           CSCCorrelatedLCTDigiCollection::Range range1 = trk->second.get((*csc).first);
           for (CSCCorrelatedLCTDigiCollection::const_iterator lct = range1.first; lct != range1.second;
                lct++, lctId++) {
             int station = (*csc).first.station() - 1;
             int cscId = (*csc).first.triggerCscId() - 1;
             int sector = (*csc).first.triggerSector() - 1;  // + ( (*csc).first.endcap()==1 ? 0 : 6 );
             int subSector = CSCTriggerNumbering::triggerSubSectorFromLabels((*csc).first);
             int tbin = lct->getBX();
             int fpga = (subSector ? subSector - 1 : station + 1);
             // If Det Id is within range
             if (sector < 0 || sector > 11 || station < 0 || station > 3 || cscId < 0 || cscId > 8 || lctId < 0 ||
                 lctId > 1) {
               edm::LogInfo("CSCTFAnalyzer: Digi are out of range: ");
               continue;
             }
             if (lct->getQuality() < 100)
               edm::LogInfo("CSCTFAnalyzer")
                   << "       Linked LCT: " << (*csc).first.endcap() << "  station: " << (station + 1)
                   << "  sector: " << (sector + 1) << "  subSector: " << subSector << "  tbin: " << tbin
                   << "  cscId: " << (cscId + 1) << "  fpga: " << (fpga + 1) << " LCT(vp=" << lct->isValid()
                   << ",qual=" << lct->getQuality() << ",wg=" << lct->getKeyWG() << ",strip=" << lct->getStrip() << ")";
             else
               edm::LogInfo("CSCTFAnalyzer")
                   << "       Linked MB stub: " << (*csc).first.endcap() << "  sector: " << (sector + 1)
                   << "  subSector: " << subSector << "  tbin: " << tbin << " MB(vp=" << lct->isValid()
                   << ",qual=" << (lct->getQuality() - 100) << ",cal=" << lct->getKeyWG() << ",flag=" << lct->getStrip()
                   << ",bc0=" << lct->getPattern() << ",phiBend=" << lct->getBend() << ",tbin=" << lct->getBX()
                   << ",id=" << lct->getMPCLink() << ",bx0=" << lct->getBX0() << ",se=" << lct->getSyncErr()
                   << ",bxn=" << lct->getCSCID() << ",phi=" << lct->getTrknmb() << ")";
           }
         }
       }
     } else
       edm::LogInfo("CSCTFAnalyzer") << "  No valid L1CSCTrackCollection products found";
   }
 }