#include <OMTFPatternMaker.h>

Inheritance diagram for OMTFPatternMaker:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

virtual void	beginJob ()

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

virtual void	endJob ()

	OMTFPatternMaker (const edm::ParameterSet &cfg)

virtual	~OMTFPatternMaker ()

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

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
const SimTrack *	findSimMuon (const edm::Event &ev, const edm::EventSetup &es, const SimTrack *previous=0)

void	writeMergedGPs ()

Private Attributes
xercesc::DOMElement *	aTopElement

edm::InputTag	g4SimTrackSrc

edm::EDGetTokenT < CSCCorrelatedLCTDigiCollection >	inputTokenCSC

edm::EDGetTokenT < L1MuDTChambPhContainer >	inputTokenDTPh

edm::EDGetTokenT < L1MuDTChambThContainer >	inputTokenDTTh

edm::EDGetTokenT < RPCDigiCollection >	inputTokenRPC

edm::EDGetTokenT < edm::SimTrackContainer >	inputTokenSimHit

bool	makeConnectionsMaps

bool	makeGoldenPatterns

bool	mergeXMLFiles

OMTFinputMaker *	myInputMaker

OMTFProcessor *	myOMTF

OMTFConfiguration *	myOMTFConfig
	OMTF objects. More...

OMTFConfigMaker *	myOMTFConfigMaker

XMLConfigWriter *	myWriter

edm::ParameterSet	theConfig

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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 35 of file OMTFPatternMaker.h.

Constructor & Destructor Documentation

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

Definition at line 22 of file OMTFPatternMaker.cc.

References edm::ParameterSet::exists(), MainPageGenerator::fName, edm::ParameterSet::getParameter(), XMLConfigWriter::initialiseXMLDocument(), inputTokenCSC, inputTokenDTPh, inputTokenDTTh, inputTokenRPC, inputTokenSimHit, makeConnectionsMaps, makeGoldenPatterns, mergeXMLFiles, myInputMaker, myOMTFConfig, myWriter, AlCaHLTBitMon_QueryRunRegistry::string, and theConfig.

                                                             :
   theConfig(cfg),
   g4SimTrackSrc(cfg.getParameter<edm::InputTag>("g4SimTrackSrc")){
 
   inputTokenDTPh = consumes<L1MuDTChambPhContainer>(theConfig.getParameter<edm::InputTag>("srcDTPh"));
   inputTokenDTTh = consumes<L1MuDTChambThContainer>(theConfig.getParameter<edm::InputTag>("srcDTTh"));
   if(!theConfig.getParameter<bool>("dropCSCPrimitives"))
     inputTokenCSC = consumes<CSCCorrelatedLCTDigiCollection>(theConfig.getParameter<edm::InputTag>("srcCSC"));
   if(!theConfig.getParameter<bool>("dropRPCPrimitives"))  
     inputTokenRPC = consumes<RPCDigiCollection>(theConfig.getParameter<edm::InputTag>("srcRPC"));
   inputTokenSimHit = consumes<edm::SimTrackContainer>(theConfig.getParameter<edm::InputTag>("g4SimTrackSrc"));
 
     
   if(!theConfig.exists("omtf")){
     edm::LogError("OMTFPatternMaker")<<"omtf configuration not found in cfg.py";
   }
   
   myInputMaker = new OMTFinputMaker();
   
   myWriter = new XMLConfigWriter();
   std::string fName = "OMTF";
   myWriter->initialiseXMLDocument(fName);
 
   makeGoldenPatterns = theConfig.getParameter<bool>("makeGoldenPatterns");
   makeConnectionsMaps = theConfig.getParameter<bool>("makeConnectionsMaps");
   mergeXMLFiles = theConfig.getParameter<bool>("mergeXMLFiles");
 
   myOMTFConfig = 0;
 }

OMTFPatternMaker::~OMTFPatternMaker ( )

virtual

Definition at line 53 of file OMTFPatternMaker.cc.

References myOMTF, myOMTFConfig, and myOMTFConfigMaker.

                                    {
 
   delete myOMTFConfig;
   delete myOMTFConfigMaker;
   delete myOMTF;
 
 }

Member Function Documentation

void OMTFPatternMaker::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	evSetup
	)

virtual

Get the simulated muon parameters

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

Connections maps are made by hand. makeConnetionsMap method provides tables for checking their consistency.

Implements edm::EDAnalyzer.

Definition at line 231 of file OMTFPatternMaker.cc.

References OMTFinputMaker::buildInputForProcessor(), OMTFProcessor::fillCounts(), findSimMuon(), edm::Event::getByToken(), edm::ParameterSet::getParameter(), OMTFinputMaker::initialize(), inputTokenCSC, inputTokenDTPh, inputTokenDTTh, inputTokenRPC, makeConnectionsMaps, OMTFConfigMaker::makeConnetionsMap(), makeGoldenPatterns, mergeXMLFiles, myInputMaker, myOMTF, myOMTFConfigMaker, l1t::omtf_pos, edm::Handle< T >::product(), and theConfig.

                                                                                   {
 
   if(mergeXMLFiles) return;
 
   const SimTrack* aSimMuon = findSimMuon(iEvent,evSetup);
   if(!aSimMuon){
     edm::LogError("OMTFPatternMaker")<<"No SimMuon found in the event!";
     return;
   }
   
   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);
 
   //l1t::tftype mtfType = l1t::tftype::bmtf;
   l1t::tftype mtfType = l1t::tftype::omtf_pos;
   //l1t::tftype mtfType = l1t::tftype::emtf_pos;
  
   for(unsigned int iProcessor=0;iProcessor<6;++iProcessor){
 
     OMTFinput myInput = myInputMaker->buildInputForProcessor(dtPhDigis.product(),
                                  dtThDigis.product(),
                                  cscDigis.product(),
                                  rpcDigis.product(),                                       
                                  iProcessor,
                                  mtfType);
     
     if(makeConnectionsMaps) myOMTFConfigMaker->makeConnetionsMap(iProcessor,myInput);
   
     if(makeGoldenPatterns) myOMTF->fillCounts(iProcessor,myInput, aSimMuon);
     
   }
 }

void OMTFPatternMaker::beginJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 91 of file OMTFPatternMaker.cc.

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

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

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

virtual

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

For making the patterns use extended pdf width in phi /Ugly hack to modify configuration parameters at runtime.

Clear existing GoldenPatterns

Reimplemented from edm::EDAnalyzer.

Definition at line 62 of file OMTFPatternMaker.cc.

References OMTFProcessor::configure(), OMTFConfiguration::configure(), edm::EventSetup::get(), edm::eventsetup::EventSetupRecord::get(), edm::ParameterSet::getParameter(), OMTFProcessor::getPatterns(), myOMTF, myOMTFConfig, OMTFConfiguration::nPdfAddrBits, fakeOmtfParams_cff::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);
 
   const L1TMuonOverlapParams* omtfParams = omtfParamsHandle.product();
   if (!omtfParams) {
     edm::LogError("L1TMuonOverlapTrackProducer") << "Could not retrieve parameters from Event Setup" << std::endl;
   }
 
   myOMTFConfig->configure(omtfParams);
   myOMTF->configure(omtfParams);
 
   myOMTFConfig->nPdfAddrBits = 14;
 
   const std::map<Key,GoldenPattern*> & theGPs = myOMTF->getPatterns();
   for(auto itGP: theGPs) itGP.second->reset();
   
 }

void OMTFPatternMaker::endJob ( void )

virtual

Order important: printPhiMap updates global vector in OMTFConfiguration

Write GPs merged by 4 above iPt=71, and by 2 below// /////////////////////////////////////////////////// 4x merging

Reimplemented from edm::EDAnalyzer.

Definition at line 101 of file OMTFPatternMaker.cc.

References OMTFProcessor::averagePatterns(), gather_cfg::cout, XMLConfigWriter::finaliseXMLDocument(), MainPageGenerator::fName, edm::ParameterSet::getParameter(), OMTFProcessor::getPatterns(), XMLConfigWriter::initialiseXMLDocument(), OMTFConfiguration::instance(), makeConnectionsMaps, makeGoldenPatterns, mergeXMLFiles, myOMTF, myOMTFConfig, myOMTFConfigMaker, myWriter, OMTFConfiguration::nPdfAddrBits, OMTFConfigMaker::printConnections(), OMTFConfigMaker::printPhiMap(), RPCConst::ptFromIpt(), GoldenPattern::reset(), AlCaHLTBitMon_QueryRunRegistry::string, theConfig, XMLConfigWriter::writeConnectionsData(), XMLConfigWriter::writeGPData(), and writeMergedGPs().

                              {
 
   if(makeGoldenPatterns && !makeConnectionsMaps){
     std::string fName = "OMTF";
     myWriter->initialiseXMLDocument(fName);
     const std::map<Key,GoldenPattern*> & myGPmap = myOMTF->getPatterns();
     for(auto itGP: myGPmap){
       if(!itGP.second->hasCounts()) continue;
       itGP.second->normalise();
     }
 
     myOMTFConfig->nPdfAddrBits = 7;
     for(auto itGP: myGPmap){
       unsigned int iPt = theConfig.getParameter<int>("ptCode")+1;
       if(iPt>31) iPt = 200*2+1;
       else iPt = RPCConst::ptFromIpt(iPt)*2.0+1;//MicroGMT has 0.5 GeV step size, with lower bin edge  (uGMT_pt_code - 1)*step_size
       if(itGP.first.thePtCode==iPt && 
      itGP.first.theCharge==theConfig.getParameter<int>("charge")){ 
     std::cout<<*itGP.second<<std::endl;      
     myWriter->writeGPData(*itGP.second);
       }
     }
     fName = "GPs.xml";
     myWriter->finaliseXMLDocument(fName);        
   }
   
   if(makeConnectionsMaps && !makeGoldenPatterns){
     std::string fName = "Connections.xml";
     unsigned int iProcessor = 0;
     myOMTFConfigMaker->printPhiMap(std::cout);
     myOMTFConfigMaker->printConnections(std::cout,iProcessor,0);
     myOMTFConfigMaker->printConnections(std::cout,iProcessor,1);
     myOMTFConfigMaker->printConnections(std::cout,iProcessor,2);
     myOMTFConfigMaker->printConnections(std::cout,iProcessor,3);
     myOMTFConfigMaker->printConnections(std::cout,iProcessor,4);
     myOMTFConfigMaker->printConnections(std::cout,iProcessor,5);
     myWriter->writeConnectionsData(OMTFConfiguration::instance()->measurements4D);
     myWriter->finaliseXMLDocument(fName);
   }
 
   if(mergeXMLFiles){
 
     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){
       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);
   }  
 }

const SimTrack * OMTFPatternMaker::findSimMuon	(	const edm::Event &	ev,
		const edm::EventSetup &	es,
		const SimTrack *	previous = `0`
	)

private

Definition at line 282 of file OMTFPatternMaker.cc.

References HLT_FULL_cff::DeltaR, edm::Event::getByToken(), inputTokenSimHit, CoreSimTrack::momentum(), mps_fire::result, and CoreSimTrack::type().

Referenced by analyze().

                                                                                                                     {
 
   const SimTrack * result = 0;
   edm::Handle<edm::SimTrackContainer> simTks;
   ev.getByToken(inputTokenSimHit,simTks);
 
   for (std::vector<SimTrack>::const_iterator it=simTks->begin(); it< simTks->end(); it++) {
     const SimTrack & aTrack = *it;
     if ( !(aTrack.type() == 13 || aTrack.type() == -13) )continue;
     if(previous && ROOT::Math::VectorUtil::DeltaR(aTrack.momentum(),previous->momentum())<0.07) continue;
     if ( !result || aTrack.momentum().pt() > result->momentum().pt()) result = &aTrack;
   }
   return result;
 }

void OMTFPatternMaker::writeMergedGPs ( )

private

Write the opposite charge.

Definition at line 172 of file OMTFPatternMaker.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 = 9;
   Key aKey = Key(0, 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;
     while(myGPmap.find(aKey)==myGPmap.end() && aKey.thePtCode<=401) ++aKey.thePtCode;    
     if(aKey.thePtCode<=401 && myGPmap.find(aKey)!=myGPmap.end()) aGP2 =  myGPmap.find(aKey)->second;
 
     if(aKey.thePtCode>71){
       ++aKey.thePtCode;
       while(myGPmap.find(aKey)==myGPmap.end() && aKey.thePtCode<=401) ++aKey.thePtCode;    
       if(aKey.thePtCode<=401 && myGPmap.find(aKey)!=myGPmap.end()) aGP3 =  myGPmap.find(aKey)->second;
 
       ++aKey.thePtCode;
       while(myGPmap.find(aKey)==myGPmap.end() && aKey.thePtCode<=401) ++aKey.thePtCode;    
       if(aKey.thePtCode<=401 && myGPmap.find(aKey)!=myGPmap.end()) aGP4 =  myGPmap.find(aKey)->second;
     }
     ++aKey.thePtCode;
     while(myGPmap.find(aKey)==myGPmap.end() && aKey.thePtCode<=401) ++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);
     }
 }