#include <MyLib/RPCPhiEff/src/RPCPhiEff.cc>

Inheritance diagram for RPCPhiEff:

Public Member Functions
	RPCPhiEff (const edm::ParameterSet &)

	~RPCPhiEff () 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

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

std::vector< ESResolverIndex > 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

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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	beginJob () override

void	endJob () override

std::string	fromCones (const edm::Event &iEvent)

std::string	fromRaw (const edm::Event &iEvent)

Private Attributes
RPCConst	m_const

edm::EDGetTokenT< edm::SimTrackContainer >	m_g4Token

std::ofstream	m_outfileC

std::ofstream	m_outfileR

edm::EDGetTokenT< std::vector< L1MuRegionalCand > >	m_rpcbToken

edm::EDGetTokenT< RPCDigiCollection >	m_rpcdigiToken

edm::EDGetTokenT< std::vector< L1MuRegionalCand > >	m_rpcfToken

RPCConst	rpcconst

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

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 40 of file RPCPhiEff.h.

Constructor & Destructor Documentation

◆ RPCPhiEff()

RPCPhiEff::RPCPhiEff ( const edm::ParameterSet & ps )

explicit

Definition at line 55 of file RPCPhiEff.cc.

References m_outfileR.

     : m_rpcbToken(consumes<std::vector<L1MuRegionalCand> >(ps.getParameter<edm::InputTag>("rpcb"))),
       m_rpcfToken(consumes<std::vector<L1MuRegionalCand> >(ps.getParameter<edm::InputTag>("rpcf"))),
       m_g4Token(consumes<edm::SimTrackContainer>(ps.getParameter<edm::InputTag>("g4"))),
       m_rpcdigiToken(consumes<RPCDigiCollection>(ps.getParameter<edm::InputTag>("rpcdigi"))) {
   //m_outfileC.open("phieffC.txt");
   m_outfileR.open("muons.txt");
 }

◆ ~RPCPhiEff()

RPCPhiEff::~RPCPhiEff ( )

override

Definition at line 64 of file RPCPhiEff.cc.

References m_outfileR.

                       {
   //m_outfileC.close();
   m_outfileR.close();
 }

Member Function Documentation

◆ analyze()

void RPCPhiEff::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 74 of file RPCPhiEff.cc.

References funct::abs(), reco::deltaR(), iEvent, m_const, m_g4Token, m_outfileR, m_rpcbToken, m_rpcfToken, muonTagProbeFilters_cff::matched, pi, l1tGTMenu_lepSeeds_cff::qual, TrackCandidateProducer_cfi::simTracks, RPCConst::towerNumFromEta(), and trackerHitRTTI::vector.

                                                                            {
   edm::Handle<std::vector<L1MuRegionalCand> > rpcBarrel;
   edm::Handle<std::vector<L1MuRegionalCand> > rpcForward;
 
   iEvent.getByToken(m_rpcbToken, rpcBarrel);
   iEvent.getByToken(m_rpcfToken, rpcForward);
 
   std::vector<edm::Handle<std::vector<L1MuRegionalCand> > > handleVec;
   handleVec.push_back(rpcBarrel);
   handleVec.push_back(rpcForward);
 
   // Fetch MCdata
   edm::Handle<edm::SimTrackContainer> simTracks;
   iEvent.getByToken(m_g4Token, simTracks);
 
   float etaGen = 0, phiGen = 0, ptGen = 0;
 
   int noOfRecMuons = 0;
   int noOfMatchedRecMuons = 0;
 
   bool firstRunForMuonMatchingCnt = true;
 
   // ask MC muons to be separated
   double deltarMin = -1;
   edm::SimTrackContainer::const_iterator simTrk = simTracks->begin();
   for (; simTrk != simTracks->end(); ++simTrk) {
     if (simTrk->type() != -13 && simTrk->type() != 13)
       continue;
 
     edm::SimTrackContainer::const_iterator simTrk2 = simTrk;
     ++simTrk2;
     for (; simTrk2 != simTracks->end(); ++simTrk2) {
       if (simTrk2->type() != -13 && simTrk2->type() != 13)
         continue;
       double drCand = reco::deltaR(simTrk2->momentum(), simTrk->momentum());
       if (drCand < deltarMin || deltarMin < 0)
         deltarMin = drCand;
     }
   }
 
   //std::cout << deltarMin << std::endl;
   if (deltarMin < 0.7 && deltarMin > 0)
     return;
 
   simTrk = simTracks->begin();
   for (; simTrk != simTracks->end(); simTrk++) {
     int type = simTrk->type();
     if (type == 13 || type == -13) {
       // Get the data
       const math::XYZTLorentzVectorD momentum = simTrk->momentum();
       etaGen = momentum.eta();
       ptGen = momentum.Pt();
       phiGen = momentum.phi();
 
       bool matched = false;
       int ptCodeRec = 0;
       int towerRec = 0;
       int phiRec = 0;
       //int muonsFound=0;
       int qual = 0;
       int ghost = 0;  // number of ghost for montecarlo muon
 
       // Iter rpc muon cands, perform delta R matching
       // todo perform matching also using eta...
       for (std::vector<edm::Handle<std::vector<L1MuRegionalCand> > >::iterator it = handleVec.begin();
            it != handleVec.end();
            ++it) {
         std::vector<L1MuRegionalCand>::const_iterator itRPC;
         for (itRPC = (*it)->begin(); itRPC != (*it)->end(); itRPC++) {
           int ptCode = itRPC->pt_packed();
           if (ptCode != 0) {
             if (firstRunForMuonMatchingCnt)
               ++noOfRecMuons;
             ptCodeRec = ptCode;
             phiRec = itRPC->phi_packed();
             qual = itRPC->quality();
             towerRec = itRPC->eta_packed();
             if (towerRec > 16) {
               towerRec = -((~towerRec & 63) + 1);
             }
             // match
             float pi = 3.14159265;
             //float phiPhys = 2*pi*float(phiRec)/144-pi;
             float phiScaled = phiGen;
             if (phiScaled < 0)
               phiScaled += 2 * pi;
             int phiHwGen = (phiScaled) / 2 / pi * 144;
 
             //std::cout << "phi " << phiGen << " " << phiHwGen << " " << phiRec
             //          << " eta " << etaGen << " " << m_const.towerNumFromEta(etaGen) << " "<< towerRec << std::endl;
             if ((std::abs(phiHwGen - phiRec) < 10) && (std::abs(m_const.towerNumFromEta(etaGen) - towerRec) < 1)) {
               if (matched) {  // we have matched m.c. earlier, this is ghost
                 ++ghost;
               }
               matched = true;
               ++noOfMatchedRecMuons;
               m_outfileR << etaGen << " " << phiGen << " " << ptGen << " " << towerRec << " " << phiRec << " "
                          << ptCodeRec << " " << qual << " " << ghost << std::endl;
             }
           }  // (ptCode != 0)
         }    // muon cands iter ends
       }      // barrell/fwd iter ends
       firstRunForMuonMatchingCnt = false;
       if (!matched) {
         m_outfileR << etaGen << " " << phiGen << " " << ptGen << " " << 0 << " " << 0 << " " << 0 << " " << 0 << " "
                    << 0 << std::endl;
       }
     }
   }
 
   edm::EventID id = iEvent.id();
   edm::EventNumber_t evNum = id.event();
   edm::RunNumber_t rnNum = id.run();
 
   if (noOfMatchedRecMuons != noOfRecMuons) {
     edm::LogInfo("RPCEffWarn") << " MuonCands " << noOfRecMuons << " matched " << noOfMatchedRecMuons << " in run "
                                << rnNum << " event " << evNum;
   }
 
   /*
     m_outfileC << etaGen << " " << phiGen << " " << ptGen << " "
     << phiRec << " " << towerRec << " " << muonsFound  << " "
     <<  fromCones(iEvent) <<std::endl;*/
   /*m_outfileR << etaGen << " " << phiGen << " " << ptGen << " "
     << phiRec << " " << towerRec << " " << muonsFound  << " "
     <<  fromRaw(iEvent) <<std::endl;*/
 
   /*
     m_outfileR << etaGen << " " << phiGen << " " << ptGen << " "
     << phiRec << " " << towerRec << " " << muonsFound  << " "
     << std::endl;
 
   */
 }

◆ beginJob()

void RPCPhiEff::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 270 of file RPCPhiEff.cc.

270 {}

◆ endJob()

void RPCPhiEff::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 273 of file RPCPhiEff.cc.

273 {}

◆ fromCones()

std::string RPCPhiEff::fromCones ( const edm::Event & iEvent )

private

Definition at line 209 of file RPCPhiEff.cc.

209 { return ""; }

◆ fromRaw()

std::string RPCPhiEff::fromRaw ( const edm::Event & iEvent )

private

Definition at line 212 of file RPCPhiEff.cc.

References gather_cfg::cout, iEvent, nano_mu_digi_cff::layer, m_rpcdigiToken, FastTimerService_cff::range, contentValuesCheck::ss, and relativeConstraints::station.

                                                    {
   std::stringstream ss;
 
   // Digi data.
   edm::Handle<RPCDigiCollection> rpcDigis;
   iEvent.getByToken(m_rpcdigiToken, rpcDigis);
 
   std::set<int> hwPlanes;
 
   RPCDigiCollection::DigiRangeIterator detUnitIt;
   for (detUnitIt = rpcDigis->begin(); detUnitIt != rpcDigis->end(); ++detUnitIt) {
     const RPCDigiCollection::Range& range = (*detUnitIt).second;
     bool hasBX0 = false;
     for (RPCDigiCollection::const_iterator digiIt = range.first; digiIt != range.second; ++digiIt) {
       if (digiIt->bx() == 0) {
         hasBX0 = true;
         break;
       }
     }
 
     if (!hasBX0)
       continue;
 
     const RPCDetId& id = (*detUnitIt).first;
     int station = id.station();
     int layer = id.layer();
     //int region = id.region();
 
     if (station == 3)
       hwPlanes.insert(5);
 
     else if (station == 4)
       hwPlanes.insert(6);
 
     else if (station == 1 && layer == 1)
       hwPlanes.insert(1);
 
     else if (station == 1 && layer == 2)
       hwPlanes.insert(2);
 
     else if (station == 2 && layer == 1)
       hwPlanes.insert(3);
 
     else if (station == 2 && layer == 2)
       hwPlanes.insert(4);
 
     else
       std::cout << "??????????????" << std::endl;
   }
 
   for (std::set<int>::iterator it = hwPlanes.begin(); it != hwPlanes.end(); ++it) {
     ss << " " << *it;
   }
 
   return ss.str();
 }

Member Data Documentation

◆ m_const

RPCConst RPCPhiEff::m_const

private

Definition at line 63 of file RPCPhiEff.h.

Referenced by analyze().

◆ m_g4Token

edm::EDGetTokenT<edm::SimTrackContainer> RPCPhiEff::m_g4Token

private

Definition at line 60 of file RPCPhiEff.h.

Referenced by analyze().

◆ m_outfileC

std::ofstream RPCPhiEff::m_outfileC

private

Definition at line 51 of file RPCPhiEff.h.

◆ m_outfileR

std::ofstream RPCPhiEff::m_outfileR

private

Definition at line 52 of file RPCPhiEff.h.

Referenced by analyze(), RPCPhiEff(), and ~RPCPhiEff().

◆ m_rpcbToken

edm::EDGetTokenT<std::vector<L1MuRegionalCand> > RPCPhiEff::m_rpcbToken

private

Definition at line 58 of file RPCPhiEff.h.

Referenced by analyze().

◆ m_rpcdigiToken

edm::EDGetTokenT<RPCDigiCollection> RPCPhiEff::m_rpcdigiToken

private

Definition at line 61 of file RPCPhiEff.h.

Referenced by fromRaw().

◆ m_rpcfToken

edm::EDGetTokenT<std::vector<L1MuRegionalCand> > RPCPhiEff::m_rpcfToken

private

Definition at line 59 of file RPCPhiEff.h.

Referenced by analyze().

◆ rpcconst

RPCConst RPCPhiEff::rpcconst

private

Definition at line 56 of file RPCPhiEff.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ RPCPhiEff()

◆ ~RPCPhiEff()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ endJob()

◆ fromCones()

◆ fromRaw()

Member Data Documentation

◆ m_const

◆ m_g4Token

◆ m_outfileC

◆ m_outfileR

◆ m_rpcbToken

◆ m_rpcdigiToken

◆ m_rpcfToken

◆ rpcconst