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

Constructor & Destructor Documentation

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

explicit

Definition at line 125 of file SimHitShifter.cc.

References gather_cfg::cout, edm::ParameterSet::getUntrackedParameter(), PFRecoTauDiscriminationByIsolation_cfi::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 155 of file SimHitShifter.cc.

{
}

Member Function Documentation

void SimHitShifter::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 280 of file SimHitShifter.cc.

{

}

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

override

Definition at line 273 of file SimHitShifter.cc.

{

}

void SimHitShifter::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 287 of file SimHitShifter.cc.

References DEFINE_FWK_MODULE.

                      {
}

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

overrideprivate

Definition at line 159 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::oct;
   using std::dec;
   
   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 118 of file SimHitShifter.cc.

Referenced by SimHitShifter().

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

Definition at line 114 of file SimHitShifter.cc.

Referenced by produce(), and SimHitShifter().