cms::ReconstructerFP420 Class Reference

#include <ReconstructerFP420.h>

Inheritance diagram for cms::ReconstructerFP420:

Public Member Functions
void	produce (edm::Event &e, const edm::EventSetup &c) override
	ReconstructerFP420 (const edm::ParameterSet &conf)
	~ReconstructerFP420 () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

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::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

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 34 of file ReconstructerFP420.h.

Constructor & Destructor Documentation

ReconstructerFP420()

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

explicit

Definition at line 28 of file ReconstructerFP420.cc.

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

~ReconstructerFP420()

cms::ReconstructerFP420::~ReconstructerFP420 ( )

override

Definition at line 50 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

produce()

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

override

Definition at line 57 of file ReconstructerFP420.cc.

References gather_cfg::cout, iEvent, input, pileupCalc::inputRange, eostools::move(), 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::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

Member Data Documentation

conf_

edm::ParameterSet cms::ReconstructerFP420::conf_

private

Definition at line 35 of file ReconstructerFP420.h.

Referenced by ReconstructerFP420().

m_genReadoutName

std::string cms::ReconstructerFP420::m_genReadoutName

private

Definition at line 41 of file ReconstructerFP420.h.

Referenced by ReconstructerFP420().

sFP420RecoMain_

FP420RecoMain* cms::ReconstructerFP420::sFP420RecoMain_

private

Definition at line 38 of file ReconstructerFP420.h.

Referenced by produce(), ReconstructerFP420(), and ~ReconstructerFP420().

trackerContainers

vstring cms::ReconstructerFP420::trackerContainers

private

Definition at line 36 of file ReconstructerFP420.h.

Referenced by produce(), and ReconstructerFP420().

verbosity

int cms::ReconstructerFP420::verbosity

private

Definition at line 39 of file ReconstructerFP420.h.

Referenced by produce(), ReconstructerFP420(), and ~ReconstructerFP420().

VtxFlag

int cms::ReconstructerFP420::VtxFlag

private

Definition at line 40 of file ReconstructerFP420.h.

Referenced by produce(), and ReconstructerFP420().