#include <ReconstructerFP420.h>

Inheritance diagram for cms::ReconstructerFP420:

Public Member Functions
void	beginJob () override

void	produce (edm::Event &e, const edm::EventSetup &c) override

	ReconstructerFP420 (const edm::ParameterSet &conf)

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

std::vector< ESProxyIndex > 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 &	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)

Private Types
typedef std::vector< std::string >	vstring

Private Attributes
edm::ParameterSet	conf_

std::string	m_genReadoutName

FP420RecoMain *	sFP420RecoMain_

vstring	trackerContainers

int	verbosity

int	VtxFlag

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	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

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 = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

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

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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

Definition at line 23 of file ReconstructerFP420.h.

Member Typedef Documentation

◆ vstring

typedef std::vector<std::string> cms::ReconstructerFP420::vstring

private

Definition at line 36 of file ReconstructerFP420.h.

Constructor & Destructor Documentation

◆ ReconstructerFP420()

cms::ReconstructerFP420::ReconstructerFP420 ( const edm::ParameterSet & conf )

explicit

Definition at line 29 of file ReconstructerFP420.cc.

                                                                     : conf_(conf) {
     edm::LogInfo("ReconstructerFP420 ") << "Enter the FP420 Reco constructer";
  
     verbosity = conf_.getUntrackedParameter<int>("VerbosityLevel");
     if (verbosity > 0) {
       std::cout << "Constructor of  ReconstructerFP420" << std::endl;
     }
  
     std::string alias(conf.getParameter<std::string>("@module_label"));
  
     produces<RecoCollectionFP420>().setBranchAlias(alias);
  
     trackerContainers.clear();
     trackerContainers = conf.getParameter<std::vector<std::string> >("ROUList");
     VtxFlag = conf.getParameter<int>("VtxFlagGenRec");
     m_genReadoutName = conf.getParameter<string>("genReadoutName");
  
     // Initialization:
     sFP420RecoMain_ = new FP420RecoMain(conf_);
   }

References SiStripOfflineCRack_cfg::alias, conf_, gather_cfg::cout, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), m_genReadoutName, sFP420RecoMain_, AlCaHLTBitMon_QueryRunRegistry::string, trackerContainers, verbosity, and VtxFlag.

◆ ~ReconstructerFP420()

cms::ReconstructerFP420::~ReconstructerFP420 ( )

override

Definition at line 51 of file ReconstructerFP420.cc.

                                           {
     if (verbosity > 0) {
       std::cout << "ReconstructerFP420:delete FP420RecoMain" << std::endl;
     }
     delete sFP420RecoMain_;
   }

References gather_cfg::cout, sFP420RecoMain_, and verbosity.

Member Function Documentation

◆ beginJob()

void cms::ReconstructerFP420::beginJob ( void )

overridevirtual

Reimplemented from edm::EDProducer.

Definition at line 59 of file ReconstructerFP420.cc.

                                     {
     if (verbosity > 0) {
       std::cout << "ReconstructerFP420:BeginJob method " << std::endl;
     }
   }

References gather_cfg::cout, and verbosity.

◆ produce()

void cms::ReconstructerFP420::produce	(	edm::Event &	e,
		const edm::EventSetup &	c
	)

overridevirtual

Implements edm::EDProducer.

Definition at line 65 of file ReconstructerFP420.cc.

                                                                                 {
     //  beginJob;
     // be lazy and include the appropriate namespaces
     using namespace edm;
     using namespace std;
  
     // Get input
  
     // Vtx info:
  
     // define GEN Vtx of Signal
     double vtxGenX = 0.;
     double vtxGenY = 0.;
     double vtxGenZ = 0.;
  
     /*
     if(VtxFlag == 1) {
       
       Handle<HepMCProduct> EvtHandle;
       try{
     iEvent.getByLabel(m_genReadoutName,EvtHandle);
       }catch(const Exception&){
     if(verbosity>0){
       std::cout << "no HepMCProduct found"<< std::endl;
     }
       }
       
       const HepMC::GenEvent* evt = EvtHandle->GetEvent() ;
       HepMC::GenParticle* proton1 = 0;
       HepMC::GenParticle* proton2 = 0;  
       double partmomcut=4000.;
       double pz1max = 0.;
       double pz2min = 0.;
       for ( HepMC::GenEvent::particle_const_iterator p = evt->particles_begin(); p != evt->particles_end(); ++p ) {
     double pz = (*p)->momentum().pz();
     //  if (((*p)->pdg_id() == ipdgproton)&&((*p)->status() == 1)&&(pz > partmomcut)){
     if( pz > partmomcut){
       if(pz > pz1max){
         proton1 = *p;pz1max=pz;
       }
     } 
     //  else if(( (*p)->pdg_id() == ipdgproton)&&((*p)->status() == 1)&&(pz < -1.*partmomcut)) {
     else if(pz < -1.*partmomcut) {
       if(pz < pz2min){
         proton2 = *p;pz2min=pz;
       }
     }
     
       }// for
       if(proton1 && !proton2){
     vtxGenX = (proton1)->production_vertex()->position().x();
     vtxGenY = (proton1)->production_vertex()->position().y();
     vtxGenZ = (proton1)->production_vertex()->position().z();
       }
       else if(proton2 && !proton1){
     vtxGenX = (proton2)->production_vertex()->position().x();
     vtxGenY = (proton2)->production_vertex()->position().y();
     vtxGenZ = (proton2)->production_vertex()->position().z();
       }
       else if(proton1 && proton2){
     if(abs((proton1)->momentum().pz()) >= abs((proton2)->momentum().pz()) ) {
       vtxGenX = (proton1)->production_vertex()->position().x();
       vtxGenY = (proton1)->production_vertex()->position().y();
       vtxGenZ = (proton1)->production_vertex()->position().z();
     }
     else {
       vtxGenX = (proton2)->production_vertex()->position().x();
       vtxGenY = (proton2)->production_vertex()->position().y();
       vtxGenZ = (proton2)->production_vertex()->position().z();
     }
       }
     }// if(VtxFlag == 1 
  
 */
  
     double VtxX = 0.;
     double VtxY = 0.;
     double VtxZ = 0.;
     if (VtxFlag == 1) {
       VtxX = vtxGenX;  // mm
       VtxY = vtxGenY;  // mm
       VtxZ = vtxGenZ;  // mm
     }
  
     // track collection:
     //A
     //   edm::Handle<ClusterCollectionFP420> icf_simhit;
     /*
     Handle<ClusterCollectionFP420> cf_simhit;
     std::vector<const ClusterCollectionFP420 *> cf_simhitvec;
     for(uint32_t i = 0; i< trackerContainers.size();i++){
       iEvent.getByLabel( trackerContainers[i], cf_simhit);
       cf_simhitvec.push_back(cf_simhit.product());   }
     std::unique_ptr<ClusterCollectionFP420 > input(new DigiCollectionFP420(cf_simhitvec));
     */
  
     //B
  
     Handle<TrackCollectionFP420> input;
     iEvent.getByLabel(trackerContainers[0], input);
  
     // Step C: create empty output collection
     auto toutput = std::make_unique<RecoCollectionFP420>();
  
     //    put zero to container info from the beginning (important! because not any detID is updated with coming of new event     !!!!!!
     // clean info of container from previous event
  
     std::vector<RecoFP420> collector;
     collector.clear();
     RecoCollectionFP420::Range inputRange;
     inputRange.first = collector.begin();
     inputRange.second = collector.end();
  
     unsigned int detID = 0;
     toutput->putclear(inputRange, detID);
  
     unsigned int StID = 1;
     toutput->putclear(inputRange, StID);
     StID = 2;
     toutput->putclear(inputRange, StID);
  
     //                                                                                                                      !!!!!!
     // if we want to keep Reco container/Collection for one event --->   uncomment the line below and vice versa
     toutput->clear();  //container_.clear() --> start from the beginning of the container
  
     //                                RUN now:                                                                                 !!!!!!
     //   startFP420RecoMain_.run(input, toutput);
     sFP420RecoMain_->run(input, toutput.get(), VtxX, VtxY, VtxZ);
     // std::cout <<"=======           ReconstructerFP420:                    end of produce     " << endl;
  
     // Step D: write output to file
     if (verbosity > 0) {
       std::cout << "ReconstructerFP420: iEvent.put(std::move(toutput)" << std::endl;
     }
     iEvent.put(std::move(toutput));
     if (verbosity > 0) {
       std::cout << "ReconstructerFP420: iEvent.put(std::move(toutput) DONE" << std::endl;
     }
   }  //produce