#include <QjetsAdder.h>

Inheritance diagram for QjetsAdder:

Public Member Functions
void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override

	QjetsAdder (const edm::ParameterSet &iConfig)

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

Private Attributes
double	cutoff_

double	jetRad_

std::string	mJetAlgo_

int	ntrial_

QjetsPlugin	qjetsAlgo_

int	QjetsPreclustering_

edm::Service< edm::RandomNumberGenerator >	rng_

edm::InputTag	src_

edm::EDGetTokenT< edm::View< reco::Jet > >	src_token_

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

Definition at line 13 of file QjetsAdder.h.

Constructor & Destructor Documentation

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

inlineexplicit

Definition at line 15 of file QjetsAdder.h.

       : src_(iConfig.getParameter<edm::InputTag>("src")),
         src_token_(consumes<edm::View<reco::Jet>>(src_)),
         qjetsAlgo_(iConfig.getParameter<double>("zcut"),
                    iConfig.getParameter<double>("dcutfctr"),
                    iConfig.getParameter<double>("expmin"),
                    iConfig.getParameter<double>("expmax"),
                    iConfig.getParameter<double>("rigidity")),
         ntrial_(iConfig.getParameter<int>("ntrial")),
         cutoff_(iConfig.getParameter<double>("cutoff")),
         jetRad_(iConfig.getParameter<double>("jetRad")),
         mJetAlgo_(iConfig.getParameter<std::string>("jetAlgo")),
         QjetsPreclustering_(iConfig.getParameter<int>("preclustering")) {
     produces<edm::ValueMap<float>>("QjetsVolatility");
   }

QjetsAdder::~QjetsAdder ( )

inlineoverride

Definition at line 31 of file QjetsAdder.h.

References iEvent, and produce().

31 {}

Member Function Documentation

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

override

Definition at line 11 of file QjetsAdder.cc.

References DEFINE_FWK_MODULE, reco::Candidate::energy(), Exception, edm::Event::getByToken(), edm::RandomNumberGenerator::getEngine(), reco::Jet::getJetConstituents(), mps_fire::i, cuy::ii, createfilelist::int, singleTopDQM_cfi::jets, dqmdumpme::k, visualization-live-secondInstance_cfg::m, SiStripPI::mean, eostools::move(), reco::LeafCandidate::pt(), edm::Event::put(), reco::Candidate::px(), reco::Candidate::py(), reco::Candidate::pz(), mathSSE::sqrt(), and edm::Event::streamID().

Referenced by ~QjetsAdder().

                                                                       {
   // read input collection
   //edm::Handle<edm::View<pat::Jet> > jets;
   edm::Handle<edm::View<reco::Jet>> jets;
   iEvent.getByToken(src_token_, jets);
 
   // prepare room for output
   std::vector<float> QjetsVolatility;
   QjetsVolatility.reserve(jets->size());
 
   for (typename edm::View<reco::Jet>::const_iterator jetIt = jets->begin(); jetIt != jets->end(); ++jetIt) {
     reco::Jet newCand(*jetIt);
 
     if (newCand.pt() < cutoff_) {
       QjetsVolatility.push_back(-1);
       continue;
     }
 
     //refill and recluster
     vector<fastjet::PseudoJet> allconstits;
     //for (unsigned k=0; k < newCand.getPFConstituents().size(); k++){
     for (unsigned k = 0; k < newCand.getJetConstituents().size(); k++) {
       const edm::Ptr<reco::Candidate> thisParticle = newCand.getJetConstituents().at(k);
       allconstits.push_back(
           fastjet::PseudoJet(thisParticle->px(), thisParticle->py(), thisParticle->pz(), thisParticle->energy()));
     }
 
     fastjet::JetDefinition jetDef(fastjet::cambridge_algorithm, jetRad_);
     if (mJetAlgo_ == "AK")
       jetDef.set_jet_algorithm(fastjet::antikt_algorithm);
     else if (mJetAlgo_ == "CA")
       jetDef.set_jet_algorithm(fastjet::cambridge_algorithm);
     else
       throw cms::Exception("GroomedJetFiller") << " unknown jet algorithm " << std::endl;
 
     fastjet::ClusterSequence thisClustering_basic(allconstits, jetDef);
     std::vector<fastjet::PseudoJet> out_jets_basic = sorted_by_pt(thisClustering_basic.inclusive_jets(cutoff_));
     //std::cout << newCand.pt() << " " << out_jets_basic.size() <<std::endl;
     if (out_jets_basic.size() !=
         1) {  // jet reclustering did not return exactly 1 jet, most likely due to the higher cutoff or large cone size. Use a recognizeable default value for this jet
       QjetsVolatility.push_back(-1);
       continue;
     }
 
     // setup objects for qjets computation
     fastjet::JetDefinition qjet_def(&qjetsAlgo_);
 
     std::vector<double> qjetmass;
 
     vector<fastjet::PseudoJet> constits;
     unsigned int nqjetconstits = out_jets_basic.at(0)
                                      .constituents()
                                      .size();  // there should always be exactly one reclsutered jet => always "at(0)"
     if (nqjetconstits < (unsigned int)QjetsPreclustering_)
       constits = out_jets_basic.at(0).constituents();
     else
       constits = out_jets_basic.at(0).associated_cluster_sequence()->exclusive_subjets_up_to(out_jets_basic.at(0),
                                                                                              QjetsPreclustering_);
 
     edm::Service<edm::RandomNumberGenerator> rng;
     CLHEP::HepRandomEngine* engine = &rng->getEngine(iEvent.streamID());
     qjetsAlgo_.SetRNEngine(engine);
     // create probabilistic recusterings
     for (unsigned int ii = 0; ii < (unsigned int)ntrial_; ii++) {
       //qjetsAlgo_.SetRandSeed(iEvent.id().event()*100 + (jetIt - jets->begin())*ntrial_ + ii );// set random seed for reprudcibility. We need a smarted scheme
       fastjet::ClusterSequence qjet_seq(constits, qjet_def);
       vector<fastjet::PseudoJet> inclusive_jets2 = sorted_by_pt(qjet_seq.inclusive_jets(cutoff_));
       if (!inclusive_jets2.empty()) {  // fill the massvalue only if the reclustering was successfull
         qjetmass.push_back(inclusive_jets2[0].m());
       }
 
     }  //end loop on trials
 
     if (qjetmass.empty()) {  //protection against dummy case
       QjetsVolatility.push_back(-1);
       continue;
     }
 
     double mean = std::accumulate(qjetmass.begin(), qjetmass.end(), 0) / qjetmass.size();
     float totalsquared = 0.;
     for (unsigned int i = 0; i < qjetmass.size(); i++) {
       totalsquared += (qjetmass[i] - mean) * (qjetmass[i] - mean);
     }
     float variance = sqrt(totalsquared / qjetmass.size());
 
     QjetsVolatility.push_back(variance / mean);
   }  //end loop on jets
 
   auto outQJV = std::make_unique<edm::ValueMap<float>>();
   edm::ValueMap<float>::Filler fillerQJV(*outQJV);
   fillerQJV.insert(jets, QjetsVolatility.begin(), QjetsVolatility.end());
   fillerQJV.fill();
 
   iEvent.put(std::move(outQJV), "QjetsVolatility");
 }

Member Data Documentation

double QjetsAdder::cutoff_

private

Definition at line 40 of file QjetsAdder.h.

double QjetsAdder::jetRad_

private

Definition at line 41 of file QjetsAdder.h.

std::string QjetsAdder::mJetAlgo_

private

Definition at line 42 of file QjetsAdder.h.

int QjetsAdder::ntrial_

private

Definition at line 39 of file QjetsAdder.h.

QjetsPlugin QjetsAdder::qjetsAlgo_

private

Definition at line 38 of file QjetsAdder.h.

int QjetsAdder::QjetsPreclustering_

private

Definition at line 43 of file QjetsAdder.h.

edm::Service<edm::RandomNumberGenerator> QjetsAdder::rng_

private

Definition at line 44 of file QjetsAdder.h.

edm::InputTag QjetsAdder::src_

private

Definition at line 36 of file QjetsAdder.h.

edm::EDGetTokenT<edm::View<reco::Jet> > QjetsAdder::src_token_

private

Definition at line 37 of file QjetsAdder.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation