SimHitShifter Class Reference

#include <simhitshifter/SimHitShifter/src/SimHitShifter.cc>

Inheritance diagram for SimHitShifter:

Public Member Functions
void	beginRun (const edm::Run &, const edm::EventSetup &) override
	SimHitShifter (const edm::ParameterSet &)
	~SimHitShifter () override
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
	~EDProducer () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
std::vector< bool > const &	recordProvenanceList () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...
	~ProducerBase () noexcept(false) 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 &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
ESProxyIndex const *	esGetTokenIndices (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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Public Attributes
std::map< int, float >	shiftinfo

Private Member Functions
void	beginJob () override
void	endJob () override
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
std::string	ShiftFileName

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex > >
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)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (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)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 106 of file SimHitShifter.cc.

Constructor & Destructor Documentation

SimHitShifter()

SimHitShifter::SimHitShifter ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 121 of file SimHitShifter.cc.

                                                           {
  std::cout << "in the constructor" << std::endl;
 
  ShiftFileName =
      iConfig.getUntrackedParameter<std::string>("ShiftFileName",
                                                 "/afs/cern.ch/user/c/carrillo/simhits/CMSSW_3_5_8_patch2/src/"
                                                 "simhitshifter/SimHitShifter/Merged_Muon_RawId_Shift.txt");
 
  //iSetup.get<MuonGeometryRecord>().get(rpcGeo);
 
  std::ifstream ifin(ShiftFileName.c_str());
 
  int rawId;
  float offset;
 
  std::cout << "In the constructor, The name of the file is " << ShiftFileName.c_str() << std::endl;
 
  if (!ifin)
    std::cout << "Problem reading the map rawId shift " << ShiftFileName.c_str() << std::endl;
  assert(ifin);
 
  while (ifin.good()) {
    ifin >> rawId >> offset;
    shiftinfo[rawId] = offset;
    std::cout << "rawId =" << rawId << " offset=" << offset << std::endl;
  }
 
  produces<edm::PSimHitContainer>("MuonCSCHits");
  produces<edm::PSimHitContainer>("MuonDTHits");
  produces<edm::PSimHitContainer>("MuonRPCHits");
}

References cms::cuda::assert(), gather_cfg::cout, edm::ParameterSet::getUntrackedParameter(), hltrates_dqm_sourceclient-live_cfg::offset, ShiftFileName, shiftinfo, and AlCaHLTBitMon_QueryRunRegistry::string.

~SimHitShifter()

SimHitShifter::~SimHitShifter ( )

override

Definition at line 153 of file SimHitShifter.cc.

{}

Member Function Documentation

beginJob()

void SimHitShifter::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 302 of file SimHitShifter.cc.

{}

beginRun()

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

overridevirtual

Reimplemented from edm::EDProducer.

Definition at line 299 of file SimHitShifter.cc.

{}

endJob()

void SimHitShifter::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 305 of file SimHitShifter.cc.

{}

produce()

void SimHitShifter::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 155 of file SimHitShifter.cc.

                                                                         {
  using namespace edm;
 
  //std::cout << " Getting the SimHits " <<std::endl;
  std::vector<edm::Handle<edm::PSimHitContainer> > theSimHitContainers;
  iEvent.getManyByType(theSimHitContainers);
  //std::cout << " The Number of sim Hits is  " << theSimHitContainers.size() <<std::endl;
 
  std::unique_ptr<edm::PSimHitContainer> pcsc(new edm::PSimHitContainer);
  std::unique_ptr<edm::PSimHitContainer> pdt(new edm::PSimHitContainer);
  std::unique_ptr<edm::PSimHitContainer> prpc(new edm::PSimHitContainer);
 
  std::vector<PSimHit> theSimHits;
 
  using std::dec;
  using std::oct;
 
  for (int i = 0; i < int(theSimHitContainers.size()); i++) {
    theSimHits.insert(theSimHits.end(), theSimHitContainers.at(i)->begin(), theSimHitContainers.at(i)->end());
  }
 
  for (std::vector<PSimHit>::const_iterator iHit = theSimHits.begin(); iHit != theSimHits.end(); iHit++) {
    DetId theDetUnitId((*iHit).detUnitId());
    DetId simdetid = DetId((*iHit).detUnitId());
 
    if (simdetid.det() != DetId::Muon)
      continue;
 
    float newtof = 0;
 
    if (simdetid.det() == DetId::Muon && simdetid.subdetId() == MuonSubdetId::RPC) {  //Only RPCs
      //std::cout<<"\t\t We have an RPC Sim Hit! in t="<<(*iHit).timeOfFlight()<<" DetId="<<(*iHit).detUnitId()<<std::endl;
      if (shiftinfo.find(simdetid.rawId()) == shiftinfo.end()) {
        std::cout << "RPC Warning the RawId = " << simdetid.det() << " | " << simdetid.rawId() << "is not in the map"
                  << std::endl;
        newtof = (*iHit).timeOfFlight();
      } else {
        newtof = (*iHit).timeOfFlight() + shiftinfo[simdetid.rawId()];
      }
 
      PSimHit hit((*iHit).entryPoint(),
                  (*iHit).exitPoint(),
                  (*iHit).pabs(),
                  newtof,
                  (*iHit).energyLoss(),
                  (*iHit).particleType(),
                  simdetid,
                  (*iHit).trackId(),
                  (*iHit).thetaAtEntry(),
                  (*iHit).phiAtEntry(),
                  (*iHit).processType());
      prpc->push_back(hit);
    } else if (simdetid.det() == DetId::Muon && simdetid.subdetId() == MuonSubdetId::DT) {  //Only DTs
      int RawId = simdetid.rawId();
      std::cout << "We found a DT simhit the RawId in Dec is";
      std::cout << dec << RawId << std::endl;
      std::cout << "and in oct" << std::endl;
      std::cout << oct << RawId << std::endl;
      std::cout << "once masked in oct " << std::endl;
      int compressedRawId = RawId / 8 / 8 / 8 / 8 / 8;
      std::cout << compressedRawId << std::endl;
      std::cout << "extendedRawId" << std::endl;
      int extendedRawId = compressedRawId * 8 * 8 * 8 * 8 * 8;
      std::cout << extendedRawId << std::endl;
      std::cout << "converted again in decimal" << std::endl;
      std::cout << dec << extendedRawId << std::endl;
 
      if (shiftinfo.find(extendedRawId) == shiftinfo.end()) {
        //std::cout<<"DT Warning the RawId = "<<extendedRawId<<"is not in the map"<<std::endl;
        newtof = (*iHit).timeOfFlight();
      } else {
        newtof = (*iHit).timeOfFlight() + shiftinfo[extendedRawId];
        std::cout << "RawId = " << extendedRawId << "is in the map " << (*iHit).timeOfFlight() << " " << newtof
                  << std::endl;
      }
 
      std::cout << "\t\t We have an DT Sim Hit! in t=" << (*iHit).timeOfFlight() << " DetId=" << (*iHit).detUnitId()
                << std::endl;
      PSimHit hit((*iHit).entryPoint(),
                  (*iHit).exitPoint(),
                  (*iHit).pabs(),
                  newtof,
                  (*iHit).energyLoss(),
                  (*iHit).particleType(),
                  simdetid,
                  (*iHit).trackId(),
                  (*iHit).thetaAtEntry(),
                  (*iHit).phiAtEntry(),
                  (*iHit).processType());
      pdt->push_back(hit);
    } else if (simdetid.det() == DetId::Muon && simdetid.subdetId() == MuonSubdetId::CSC) {  //Only CSCs
      //std::cout<<"\t\t We have an CSC Sim Hit! in t="<<(*iHit).timeOfFlight()<<" DetId="<<(*iHit).detUnitId()<<std::endl;
 
      CSCDetId TheCSCDetId = CSCDetId(simdetid);
      CSCDetId TheChamberDetId = TheCSCDetId.chamberId();
 
      if (shiftinfo.find(TheChamberDetId.rawId()) == shiftinfo.end()) {
        std::cout << "The RawId is not in the map,perhaps it is on the CSCs station 1 ring 4" << std::endl;
        if (TheChamberDetId.station() == 1 && TheChamberDetId.ring() == 4) {
          CSCDetId TheChamberDetIdNoring4 = CSCDetId(TheChamberDetId.endcap(),
                                                     TheChamberDetId.station(),
                                                     1  //1 instead of 4
                                                     ,
                                                     TheChamberDetId.chamber(),
                                                     TheChamberDetId.layer());
 
          if (shiftinfo.find(TheChamberDetIdNoring4.rawId()) == shiftinfo.end()) {
            std::cout << "CSC Warning the RawId = " << TheChamberDetIdNoring4 << " " << TheChamberDetIdNoring4.rawId()
                      << "is not in the map" << std::endl;
            newtof = (*iHit).timeOfFlight();
          } else {
            newtof = (*iHit).timeOfFlight() + shiftinfo[TheChamberDetIdNoring4.rawId()];
          }
        }
      } else {
        newtof = (*iHit).timeOfFlight() + shiftinfo[TheChamberDetId.rawId()];
      }
 
      PSimHit hit((*iHit).entryPoint(),
                  (*iHit).exitPoint(),
                  (*iHit).pabs(),
                  newtof,
                  (*iHit).energyLoss(),
                  (*iHit).particleType(),
                  simdetid,
                  (*iHit).trackId(),
                  (*iHit).thetaAtEntry(),
                  (*iHit).phiAtEntry(),
                  (*iHit).processType());
 
      std::cout << "CSC check newtof" << newtof << " " << (*iHit).timeOfFlight() << std::endl;
      if (newtof == (*iHit).timeOfFlight())
        std::cout << "Warning!!!" << std::endl;
      pcsc->push_back(hit);
    }
  }
 
  std::cout << "Putting collections in the event" << std::endl;
 
  iEvent.put(std::move(pcsc), "MuonCSCHits");
  iEvent.put(std::move(pdt), "MuonDTHits");
  iEvent.put(std::move(prpc), "MuonRPCHits");
}

References CSCDetId::chamber(), CSCDetId::chamberId(), gather_cfg::cout, MuonSubdetId::CSC, TauDecayModes::dec, MuonSubdetId::DT, CSCDetId::endcap(), mps_fire::i, iEvent, createfilelist::int, CSCDetId::layer(), eostools::move(), DetId::Muon, DetId::rawId(), CSCDetId::ring(), MuonSubdetId::RPC, shiftinfo, and CSCDetId::station().

Member Data Documentation

ShiftFileName

std::string SimHitShifter::ShiftFileName

private

Definition at line 115 of file SimHitShifter.cc.

Referenced by SimHitShifter().

shiftinfo

std::map<int, float> SimHitShifter::shiftinfo

Definition at line 112 of file SimHitShifter.cc.

Referenced by produce(), and SimHitShifter().