#include <L1TMuonOverlapTrackProducer.h>

Inheritance diagram for L1TMuonOverlapTrackProducer:

Public Member Functions
virtual void	beginJob ()

virtual void	beginRun (edm::Run const &run, edm::EventSetup const &iSetup)

virtual void	endJob ()

	L1TMuonOverlapTrackProducer (const edm::ParameterSet &)

virtual void	produce (edm::Event &, const edm::EventSetup &)

	~L1TMuonOverlapTrackProducer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

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

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

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

	EDConsumerBase ()

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
void	processCandidates (unsigned int iProcessor, int bx, std::auto_ptr< l1t::RegionalMuonCandBxCollection > &myCands, l1t::RegionalMuonCandBxCollection &myOTFCandidates, l1t::tftype mtfType)

void	writeMergedGPs ()

Private Attributes
xercesc::DOMElement *	aTopElement

bool	dumpDetailedResultToXML

bool	dumpGPToXML

bool	dumpResultToXML

edm::EDGetTokenT < CSCCorrelatedLCTDigiCollection >	inputTokenCSC

edm::EDGetTokenT < L1MuDTChambPhContainer >	inputTokenDTPh

edm::EDGetTokenT < L1MuDTChambThContainer >	inputTokenDTTh

edm::EDGetTokenT < RPCDigiCollection >	inputTokenRPC

OMTFinputMaker *	myInputMaker

OMTFinput *	myInputXML

OMTFProcessor *	myOMTF

OMTFConfiguration *	myOMTFConfig
	OMTF objects. More...

OMTFConfigMaker *	myOMTFConfigMaker

OMTFSorter *	mySorter

XMLConfigWriter *	myWriter

std::shared_ptr < L1TMuonOverlapParams >	omtfParams

edm::ParameterSet	theConfig

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDProducer
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)

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 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)

Detailed Description

Definition at line 40 of file L1TMuonOverlapTrackProducer.h.

Constructor & Destructor Documentation

L1TMuonOverlapTrackProducer::L1TMuonOverlapTrackProducer ( const edm::ParameterSet & cfg )

Definition at line 23 of file L1TMuonOverlapTrackProducer.cc.

References dumpDetailedResultToXML, dumpGPToXML, dumpResultToXML, edm::ParameterSet::exists(), MainPageGenerator::fName, edm::ParameterSet::getParameter(), XMLConfigWriter::initialiseXMLDocument(), inputTokenCSC, inputTokenDTPh, inputTokenDTTh, inputTokenRPC, myInputMaker, myOMTF, myOMTFConfig, mySorter, myWriter, AlCaHLTBitMon_QueryRunRegistry::string, and theConfig.

                                                                                   :theConfig(cfg){
 
   produces<l1t::RegionalMuonCandBxCollection >("OMTF");
 
   inputTokenDTPh = consumes<L1MuDTChambPhContainer>(theConfig.getParameter<edm::InputTag>("srcDTPh"));
   inputTokenDTTh = consumes<L1MuDTChambThContainer>(theConfig.getParameter<edm::InputTag>("srcDTTh"));
   inputTokenCSC = consumes<CSCCorrelatedLCTDigiCollection>(theConfig.getParameter<edm::InputTag>("srcCSC"));
   inputTokenRPC = consumes<RPCDigiCollection>(theConfig.getParameter<edm::InputTag>("srcRPC"));
 
   if(!theConfig.exists("omtf")){
     edm::LogError("L1TMuonOverlapTrackProducer")<<"omtf configuration not found in cfg.py";
   }
 
   myInputMaker = new OMTFinputMaker();
   mySorter = new OMTFSorter();
   myOMTF = 0;
   myWriter = 0;
 
   dumpResultToXML = theConfig.getParameter<bool>("dumpResultToXML");
   dumpDetailedResultToXML = theConfig.getParameter<bool>("dumpDetailedResultToXML");
   dumpGPToXML = theConfig.getParameter<bool>("dumpGPToXML");
   theConfig.getParameter<std::string>("XMLDumpFileName");
 
   if(dumpResultToXML || dumpGPToXML){
     myWriter = new XMLConfigWriter();
     std::string fName = "OMTF";
     myWriter->initialiseXMLDocument(fName);
   }
 
   myOMTFConfig = 0;
 }

L1TMuonOverlapTrackProducer::~L1TMuonOverlapTrackProducer ( )

Definition at line 56 of file L1TMuonOverlapTrackProducer.cc.

References myInputMaker, myOMTF, myOMTFConfig, mySorter, and myWriter.

                                                          {
 
   delete myOMTFConfig;
   delete myOMTF;
 
   delete myInputMaker;
   delete mySorter;
 
   if (myWriter) delete myWriter;
 
 }

Member Function Documentation

void L1TMuonOverlapTrackProducer::beginJob ( void )

virtual

Reimplemented from edm::EDProducer.

Definition at line 69 of file L1TMuonOverlapTrackProducer.cc.

References edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), myOMTF, myOMTFConfig, and theConfig.

                                           {
 
   if(theConfig.exists("omtf")){
     myOMTFConfig = new OMTFConfiguration(theConfig.getParameter<edm::ParameterSet>("omtf"));
     myOMTF = new OMTFProcessor(theConfig.getParameter<edm::ParameterSet>("omtf"));
   }
 }

void L1TMuonOverlapTrackProducer::beginRun	(	edm::Run const &	run,
		edm::EventSetup const &	iSetup
	)

virtual

If configuration is read from XML do not look at the DB.

Reimplemented from edm::EDProducer.

Definition at line 171 of file L1TMuonOverlapTrackProducer.cc.

References OMTFProcessor::configure(), OMTFConfiguration::configure(), edm::EventSetup::get(), edm::eventsetup::EventSetupRecord::get(), edm::ParameterSet::getParameter(), myOMTF, myOMTFConfig, omtfParams, and theConfig.

                                                                                         {
 
   if(theConfig.getParameter<edm::ParameterSet>("omtf").getParameter<bool>("configFromXML")) return;  
 
   const L1TMuonOverlapParamsRcd& omtfParamsRcd = iSetup.get<L1TMuonOverlapParamsRcd>();
   
   edm::ESHandle<L1TMuonOverlapParams> omtfParamsHandle;
   omtfParamsRcd.get(omtfParamsHandle);
 
   omtfParams = std::unique_ptr<L1TMuonOverlapParams>(new L1TMuonOverlapParams(*omtfParamsHandle.product()));
   if (!omtfParams) {
     edm::LogError("L1TMuonOverlapTrackProducer") << "Could not retrieve parameters from Event Setup" << std::endl;
   }
 
   myOMTFConfig->configure(omtfParams);
   myOMTF->configure(omtfParams);
   
 }

void L1TMuonOverlapTrackProducer::endJob ( void )

virtual

Write GPs merged by 4 above iPt19, and by 2 below// /////////////////////////////////////////////////// 4x merging

Reimplemented from edm::EDProducer.

Definition at line 78 of file L1TMuonOverlapTrackProducer.cc.

References OMTFProcessor::averagePatterns(), gather_cfg::cout, dumpGPToXML, dumpResultToXML, XMLConfigWriter::finaliseXMLDocument(), MainPageGenerator::fName, edm::ParameterSet::getParameter(), OMTFProcessor::getPatterns(), XMLConfigWriter::initialiseXMLDocument(), myOMTF, myWriter, GoldenPattern::reset(), AlCaHLTBitMon_QueryRunRegistry::string, theConfig, XMLConfigWriter::writeGPData(), and writeMergedGPs().

                                         {
 
   if(dumpResultToXML && !dumpGPToXML){
     std::string fName = theConfig.getParameter<std::string>("XMLDumpFileName");
     myWriter->finaliseXMLDocument(fName);
   }
 
   if(dumpGPToXML && !dumpResultToXML){
 
     GoldenPattern *dummy = new GoldenPattern(Key(0,0,0));
     dummy->reset();
 
     std::string fName = "OMTF";
     myWriter->initialiseXMLDocument(fName);
     const std::map<Key,GoldenPattern*> & myGPmap = myOMTF->getPatterns();
     for(auto itGP: myGPmap){
       if(itGP.second->key().thePtCode==6) std::cout<<*itGP.second<<std::endl;
       //myWriter->writeGPData(*itGP.second);
       if(itGP.second->key().thePtCode>5) myWriter->writeGPData(*itGP.second,*dummy, *dummy, *dummy);
     }
     fName = "GPs.xml";
     myWriter->finaliseXMLDocument(fName);
     fName = "OMTF";
     myWriter->initialiseXMLDocument(fName);
     myOMTF->averagePatterns(1);
     myOMTF->averagePatterns(-1);
 
     writeMergedGPs();
 
     fName = "GPs_4x.xml";
     myWriter->finaliseXMLDocument(fName);
   }
 }

void L1TMuonOverlapTrackProducer::processCandidates	(	unsigned int	iProcessor,
		int	bx,
		std::auto_ptr< l1t::RegionalMuonCandBxCollection > &	myCands,
		l1t::RegionalMuonCandBxCollection &	myOTFCandidates,
		l1t::tftype	mtfType
	)

private

Set local 0 at iProcessor x 15 deg

Definition at line 282 of file L1TMuonOverlapTrackProducer.cc.

References BXVector< T >::at(), OMTFConfiguration::globalPhiStart(), l1t::RegionalMuonCand::hwPhi(), l1t::RegionalMuonCand::hwPt(), OMTFConfiguration::nPhiBins, OMTFConfiguration::nPhiBits, funct::pow(), l1t::RegionalMuonCand::setHwPhi(), l1t::RegionalMuonCand::setTFIdentifiers(), and BXVector< T >::size().

Referenced by produce().

                                        {
 
   int procOffset = OMTFConfiguration::globalPhiStart(iProcessor);
   if(procOffset<0) procOffset+=(int)OMTFConfiguration::nPhiBins;
   procOffset-=(15+iProcessor*60)/360.0*OMTFConfiguration::nPhiBins;
   int lowScaleEnd = pow(2,OMTFConfiguration::nPhiBits-1);
 
     for(unsigned int iCand=0; iCand<myOTFCandidates.size(bx); ++iCand){
       // shift phi from processor to global coordinates
       l1t::RegionalMuonCand cand = myOTFCandidates.at(bx, iCand);
       int phiValue = (cand.hwPhi()+procOffset+lowScaleEnd);
       if(phiValue>=(int)OMTFConfiguration::nPhiBins) phiValue-=OMTFConfiguration::nPhiBins;
       phiValue/=10; //MicroGMT has 10x coarser scale than OMTF
 
       cand.setHwPhi(phiValue);
       cand.setTFIdentifiers(iProcessor,mtfType);
       // store candidate
       if(cand.hwPt()) myCands->push_back(bx, cand);
     }
 }

void L1TMuonOverlapTrackProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	evSetup
	)

virtual

Filter digis by dropping digis from selected (by cfg.py) subsystems

Loop over all processors, each covering 60 deg in phi

Input data with phi ranges shifted for each processor, so it fits 11 bits range

Results for each GP in each logic region of given processor

Shift phi scales, and put MicroGMT candidates into myCands collection

Write data to XML file

Dump only regions, where a candidate was found

Implements edm::EDProducer.

Definition at line 192 of file L1TMuonOverlapTrackProducer.cc.

References aTopElement, OMTFinputMaker::buildInputForProcessor(), dumpDetailedResultToXML, dumpResultToXML, edm::EventID::event(), edm::Event::getByToken(), edm::ParameterSet::getParameter(), edm::EventBase::id(), OMTFinputMaker::initialize(), inputTokenCSC, inputTokenDTPh, inputTokenDTTh, inputTokenRPC, myInputMaker, myOMTF, mySorter, myWriter, OMTFConfiguration::nTestRefHits, l1t::omtf_neg, l1t::omtf_pos, processCandidates(), OMTFProcessor::processInput(), edm::Handle< T >::product(), InternalObj::pt, edm::Event::put(), OMTFProcessor::shiftInput(), OMTFSorter::sortProcessor(), OMTFSorter::sortRefHitResults(), theConfig, XMLConfigWriter::writeCandidateData(), XMLConfigWriter::writeEventData(), XMLConfigWriter::writeEventHeader(), and XMLConfigWriter::writeResultsData().

                                                                                        {
 
   std::ostringstream myStr;
 
   myInputMaker->initialize(evSetup);
 
   edm::Handle<L1MuDTChambPhContainer> dtPhDigis;
   edm::Handle<L1MuDTChambThContainer> dtThDigis;
   edm::Handle<CSCCorrelatedLCTDigiCollection> cscDigis;
   edm::Handle<RPCDigiCollection> rpcDigis;
 
   if(!theConfig.getParameter<bool>("dropDTPrimitives")){
     iEvent.getByToken(inputTokenDTPh,dtPhDigis);
     iEvent.getByToken(inputTokenDTTh,dtThDigis);
   }
   if(!theConfig.getParameter<bool>("dropRPCPrimitives")) iEvent.getByToken(inputTokenRPC,rpcDigis);  
   if(!theConfig.getParameter<bool>("dropCSCPrimitives")) iEvent.getByToken(inputTokenCSC,cscDigis);
   
   std::auto_ptr<l1t::RegionalMuonCandBxCollection > myCands(new l1t::RegionalMuonCandBxCollection);
 
   if(dumpResultToXML) aTopElement = myWriter->writeEventHeader(iEvent.id().event());
 
   // NOTE: assuming all is for bx 0
   int bx = 0;
 
   for(unsigned int iProcessor=0;iProcessor<6;++iProcessor){
 
     //myStr<<" iProcessor: "<<iProcessor;
 
     const OMTFinput *myInputPos = myInputMaker->buildInputForProcessor(dtPhDigis.product(),
                                        dtThDigis.product(),
                                        cscDigis.product(),
                                        rpcDigis.product(),                                     
                                        iProcessor,
                                        l1t::tftype::omtf_pos);
 
     OMTFinput myShiftedInputPos =  myOMTF->shiftInput(iProcessor,*myInputPos);
 
     const OMTFinput *myInputNeg = myInputMaker->buildInputForProcessor(dtPhDigis.product(),
                                        dtThDigis.product(),
                                        cscDigis.product(),
                                        rpcDigis.product(),                                     
                                        iProcessor,
                                        l1t::tftype::omtf_neg);
 
     
     OMTFinput myShiftedInputNeg =  myOMTF->shiftInput(iProcessor,*myInputNeg);
 
     l1t::RegionalMuonCandBxCollection myOTFCandidatesPos, myOTFCandidatesNeg;
     //Retreive all candidates returned by sorter: upto 3 non empty ones with different phi or charge
     const std::vector<OMTFProcessor::resultsMap> & myResultsNeg = myOMTF->processInput(iProcessor,myShiftedInputNeg);
     mySorter->sortProcessor(myResultsNeg, myOTFCandidatesNeg, bx);
 
     const std::vector<OMTFProcessor::resultsMap> & myResultsPos = myOMTF->processInput(iProcessor,myShiftedInputPos);
     mySorter->sortProcessor(myResultsPos, myOTFCandidatesPos, bx);
 
     processCandidates(iProcessor, bx, myCands, myOTFCandidatesPos, l1t::tftype::omtf_pos);
     processCandidates(iProcessor, bx, myCands, myOTFCandidatesNeg, l1t::tftype::omtf_neg);
 
     if(dumpResultToXML){
       xercesc::DOMElement * aProcElement = myWriter->writeEventData(aTopElement,iProcessor,myShiftedInputPos);
       for(unsigned int iRefHit=0;iRefHit<OMTFConfiguration::nTestRefHits;++iRefHit){
     InternalObj myCand = mySorter->sortRefHitResults(myResultsPos[iRefHit],0);//charge=0 means ignore charge
     if(myCand.pt){
       myWriter->writeCandidateData(aProcElement,iRefHit,myCand);
       if(dumpDetailedResultToXML){
         for(auto & itKey: myResultsNeg[iRefHit]) myWriter->writeResultsData(aProcElement,
                                         iRefHit,
                                         itKey.first,itKey.second);
       }
     }
       }
     }
   }
 
   //dumpResultToXML = true;
   myStr<<" Number of candidates: "<<myCands->size(bx);
   edm::LogInfo("OMTFOMTFProducer")<<myStr.str();
 
   iEvent.put(myCands, "OMTF");
 }

void L1TMuonOverlapTrackProducer::writeMergedGPs ( )

private

Write the opposite charge.

Definition at line 116 of file L1TMuonOverlapTrackProducer.cc.

References OMTFProcessor::getPatterns(), GoldenPattern::key(), myOMTF, myWriter, GoldenPattern::reset(), Key::theCharge, Key::thePtCode, and XMLConfigWriter::writeGPData().

Referenced by endJob().

                                                 {
   
   const std::map<Key,GoldenPattern*> & myGPmap = myOMTF->getPatterns();
 
   GoldenPattern *dummy = new GoldenPattern(Key(0,0,0));
   dummy->reset();
 
   unsigned int iPtMin = 6;
   Key aKey = Key(1, iPtMin,-1);
   while(myGPmap.find(aKey)!=myGPmap.end()){
 
     GoldenPattern *aGP1 = myGPmap.find(aKey)->second;
     GoldenPattern *aGP2 = dummy;
     GoldenPattern *aGP3 = dummy;
     GoldenPattern *aGP4 = dummy;
 
     ++aKey.thePtCode;
     if(aKey.thePtCode<=31 && myGPmap.find(aKey)!=myGPmap.end()) aGP2 =  myGPmap.find(aKey)->second;
 
     if(aKey.thePtCode>19){
       ++aKey.thePtCode;
       if(aKey.thePtCode<=31 && myGPmap.find(aKey)!=myGPmap.end()) aGP3 =  myGPmap.find(aKey)->second;
 
       ++aKey.thePtCode;
       if(aKey.thePtCode<=31 && myGPmap.find(aKey)!=myGPmap.end()) aGP4 =  myGPmap.find(aKey)->second;
     }
     ++aKey.thePtCode;
     myWriter->writeGPData(*aGP1,*aGP2, *aGP3, *aGP4);
 
     Key aTmpKey = aGP1->key();
     aTmpKey.theCharge = 1;
     if(myGPmap.find(aTmpKey)!=myGPmap.end()) aGP1 =  myGPmap.find(aTmpKey)->second;
     else aGP1 = dummy;
 
     aTmpKey = aGP2->key();
     aTmpKey.theCharge = 1;
     if(myGPmap.find(aTmpKey)!=myGPmap.end()) aGP2 =  myGPmap.find(aTmpKey)->second;
     else aGP2 = dummy;
 
     aTmpKey = aGP3->key();
     aTmpKey.theCharge = 1;
     if(myGPmap.find(aTmpKey)!=myGPmap.end()) aGP3 =  myGPmap.find(aTmpKey)->second;
     else aGP3 = dummy;
 
     aTmpKey = aGP4->key();
     aTmpKey.theCharge = 1;
     if(myGPmap.find(aTmpKey)!=myGPmap.end()) aGP4 =  myGPmap.find(aTmpKey)->second;
     else aGP4 = dummy;
     
     myWriter->writeGPData(*aGP1,*aGP2, *aGP3, *aGP4);
     }
 }