#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

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 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::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 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 32 of file QjetsAdder.h.

References iEvent, and produce().

32 {}

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, fwrapper::jets, gen::k, funct::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 41 of file QjetsAdder.h.

double QjetsAdder::jetRad_

private

Definition at line 42 of file QjetsAdder.h.

std::string QjetsAdder::mJetAlgo_

private

Definition at line 43 of file QjetsAdder.h.

int QjetsAdder::ntrial_

private

Definition at line 40 of file QjetsAdder.h.

QjetsPlugin QjetsAdder::qjetsAlgo_

private

Definition at line 39 of file QjetsAdder.h.

int QjetsAdder::QjetsPreclustering_

private

Definition at line 44 of file QjetsAdder.h.

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

private

Definition at line 45 of file QjetsAdder.h.

edm::InputTag QjetsAdder::src_

private

Definition at line 37 of file QjetsAdder.h.

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

private

Definition at line 38 of file QjetsAdder.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation