#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

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

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)

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	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 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

Definition at line 26 of file ReconstructerFP420.h.

Member Typedef Documentation

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

private

Definition at line 41 of file ReconstructerFP420.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 30 of file ReconstructerFP420.cc.

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

                                                                    :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_);
     
   }

cms::ReconstructerFP420::~ReconstructerFP420 ( )

override

Definition at line 57 of file ReconstructerFP420.cc.

References gather_cfg::cout, sFP420RecoMain_, and verbosity.

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

Member Function Documentation

void cms::ReconstructerFP420::beginJob ( void )

overridevirtual

Reimplemented from edm::EDProducer.

Definition at line 65 of file ReconstructerFP420.cc.

References gather_cfg::cout, and verbosity.

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

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

override

Definition at line 72 of file ReconstructerFP420.cc.

References gather_cfg::cout, edm::Event::getByLabel(), input, estimatePileup::inputRange, eostools::move(), edm::Event::put(), FP420RecoMain::run(), sFP420RecoMain_, trackerContainers, verbosity, and VtxFlag.

   {
     //  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::auto_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