#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

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)

	~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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

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
ProducesCollector	producesCollector ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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 106 of file SimHitShifter.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 121 of file SimHitShifter.cc.

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

                                                            {
   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");
 }

SimHitShifter::~SimHitShifter ( )

override

Definition at line 153 of file SimHitShifter.cc.

153 {}

Member Function Documentation

void SimHitShifter::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 302 of file SimHitShifter.cc.

302 {}

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

override

Definition at line 299 of file SimHitShifter.cc.

299 {}

void SimHitShifter::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 305 of file SimHitShifter.cc.

References DEFINE_FWK_MODULE.

305 {}

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

overrideprivate

Definition at line 155 of file SimHitShifter.cc.

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

                                                                          {
   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");
 }

Member Data Documentation

std::string SimHitShifter::ShiftFileName

private

Definition at line 115 of file SimHitShifter.cc.

Referenced by SimHitShifter().

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

Definition at line 112 of file SimHitShifter.cc.

Referenced by produce(), and SimHitShifter().

Public Member Functions

Public Attributes

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation