RPCPhiEff Class Reference

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

Inheritance diagram for RPCPhiEff:

Public Member Functions
	RPCPhiEff (const edm::ParameterSet &)
	~RPCPhiEff ()
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 (const std::string &iProcessName, std::vector< const char * > &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
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob (const edm::EventSetup &)
virtual void	endJob ()
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::EDAnalyzer
typedef EDAnalyzer	ModuleType
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

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

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

Definition at line 40 of file RPCPhiEff.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 57 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

(

)

Definition at line 75 of file RPCPhiEff.cc.

References m_outfileR.

{

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

}

Member Function Documentation

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

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 90 of file RPCPhiEff.cc.

References funct::abs(), reco::deltaR(), edm::EventID::event(), edm::Event::getByToken(), edm::EventBase::id(), m_const, m_g4Token, m_outfileR, m_rpcbToken, m_rpcfToken, pi, and RPCConst::towerNumFromEta().

{


  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 noOfMuons = 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();
      noOfMuons++;

      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;

  */


}

void RPCPhiEff::beginJob ( const edm::EventSetup &  )

privatevirtual

Definition at line 328 of file RPCPhiEff.cc.

{
}

void RPCPhiEff::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 333 of file RPCPhiEff.cc.

{
}

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

private

Definition at line 246 of file RPCPhiEff.cc.

                                                     {

  return "";
}

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

private

Definition at line 252 of file RPCPhiEff.cc.

References gather_cfg::cout, edm::Event::getByToken(), m_rpcdigiToken, 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

RPCConst RPCPhiEff::m_const

private

Definition at line 65 of file RPCPhiEff.h.

Referenced by analyze().

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

private

Definition at line 62 of file RPCPhiEff.h.

Referenced by analyze().

std::ofstream RPCPhiEff::m_outfileC

private

Definition at line 52 of file RPCPhiEff.h.

std::ofstream RPCPhiEff::m_outfileR

private

Definition at line 53 of file RPCPhiEff.h.

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

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

private

Definition at line 60 of file RPCPhiEff.h.

Referenced by analyze().

edm::EDGetTokenT<RPCDigiCollection> RPCPhiEff::m_rpcdigiToken

private

Definition at line 63 of file RPCPhiEff.h.

Referenced by fromRaw().

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

private

Definition at line 61 of file RPCPhiEff.h.

Referenced by analyze().

RPCConst RPCPhiEff::rpcconst

private

Definition at line 58 of file RPCPhiEff.h.