CMS 3D CMS Logo
List of all members | Public Member Functions | Private Attributes
QjetsAdder Class Reference

#include <QjetsAdder.h>

Inheritance diagram for QjetsAdder:
edm::EDProducer edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

void produce (edm::Event &iEvent, const edm::EventSetup &iSetup)
 
 QjetsAdder (const edm::ParameterSet &iConfig)
 
virtual ~QjetsAdder ()
 
- Public Member Functions inherited from edm::EDProducer
 EDProducer ()
 
ModuleDescription const & moduleDescription () const
 
virtual ~EDProducer ()
 
- 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 ()
 
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, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)
 
virtual ~ProducerBase () noexcept(false)
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
 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 & itemsToGetFromEvent () 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
 
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
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)
 
- 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.

15  :
16  src_(iConfig.getParameter<edm::InputTag>("src")),
17  src_token_(consumes<edm::View<reco::Jet>>(src_)),
18  qjetsAlgo_( iConfig.getParameter<double>("zcut"),
19  iConfig.getParameter<double>("dcutfctr"),
20  iConfig.getParameter<double>("expmin"),
21  iConfig.getParameter<double>("expmax"),
22  iConfig.getParameter<double>("rigidity")),
23  ntrial_(iConfig.getParameter<int>("ntrial")),
24  cutoff_(iConfig.getParameter<double>("cutoff")),
25  jetRad_(iConfig.getParameter<double>("jetRad")),
26  mJetAlgo_(iConfig.getParameter<std::string>("jetAlgo")) ,
27  QjetsPreclustering_(iConfig.getParameter<int>("preclustering"))
28  {
29  produces<edm::ValueMap<float> >("QjetsVolatility");
30  }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
QjetsAdder::ntrial_
int ntrial_
Definition: QjetsAdder.h:40
QjetsAdder::jetRad_
double jetRad_
Definition: QjetsAdder.h:42
edm::EDConsumerBase::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: EDConsumerBase.h:101
QjetsAdder::qjetsAlgo_
QjetsPlugin qjetsAlgo_
Definition: QjetsAdder.h:39
QjetsAdder::QjetsPreclustering_
int QjetsPreclustering_
Definition: QjetsAdder.h:44
QjetsAdder::cutoff_
double cutoff_
Definition: QjetsAdder.h:41
QjetsAdder::src_token_
edm::EDGetTokenT< edm::View< reco::Jet > > src_token_
Definition: QjetsAdder.h:38
QjetsAdder::mJetAlgo_
std::string mJetAlgo_
Definition: QjetsAdder.h:43
QjetsAdder::src_
edm::InputTag src_
Definition: QjetsAdder.h:37
virtual QjetsAdder::~QjetsAdder ( )
inlinevirtual

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 
)

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, RecoTauDiscriminantConfiguration::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().

11  {
12  // read input collection
13  //edm::Handle<edm::View<pat::Jet> > jets;
14  edm::Handle<edm::View<reco::Jet> > jets;
15  iEvent.getByToken(src_token_, jets);
16 
17  // prepare room for output
18  std::vector<float> QjetsVolatility; QjetsVolatility.reserve(jets->size());
19 
20  for ( typename edm::View<reco::Jet>::const_iterator jetIt = jets->begin() ; jetIt != jets->end() ; ++jetIt ) {
21  reco::Jet newCand(*jetIt);
22 
23  if(newCand.pt()<cutoff_)
24  {
25  QjetsVolatility.push_back(-1);
26  continue;
27  }
28 
29  //refill and recluster
30  vector<fastjet::PseudoJet> allconstits;
31  //for (unsigned k=0; k < newCand.getPFConstituents().size(); k++){
32  for (unsigned k=0; k < newCand.getJetConstituents().size(); k++){
33  const edm::Ptr<reco::Candidate> thisParticle = newCand.getJetConstituents().at(k);
34  allconstits.push_back( fastjet::PseudoJet( thisParticle->px(), thisParticle->py(), thisParticle->pz(), thisParticle->energy() ) );
35  }
36 
37  fastjet::JetDefinition jetDef(fastjet::cambridge_algorithm, jetRad_);
38  if (mJetAlgo_== "AK") jetDef.set_jet_algorithm( fastjet::antikt_algorithm );
39  else if (mJetAlgo_ == "CA") jetDef.set_jet_algorithm( fastjet::cambridge_algorithm );
40  else throw cms::Exception("GroomedJetFiller") << " unknown jet algorithm " << std::endl;
41 
42  fastjet::ClusterSequence thisClustering_basic(allconstits, jetDef);
43  std::vector<fastjet::PseudoJet> out_jets_basic = sorted_by_pt(thisClustering_basic.inclusive_jets(cutoff_));
44  //std::cout << newCand.pt() << " " << out_jets_basic.size() <<std::endl;
45  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
46  QjetsVolatility.push_back(-1);
47  continue;
48  }
49 
50  // setup objects for qjets computation
51  fastjet::JetDefinition qjet_def(&qjetsAlgo_);
52 
53  std::vector<double> qjetmass;
54 
55  vector<fastjet::PseudoJet> constits;
56  unsigned int nqjetconstits = out_jets_basic.at(0).constituents().size(); // there should always be exactly one reclsutered jet => always "at(0)"
57  if (nqjetconstits < (unsigned int) QjetsPreclustering_) constits = out_jets_basic.at(0).constituents();
58  else constits = out_jets_basic.at(0).associated_cluster_sequence()->exclusive_subjets_up_to(out_jets_basic.at(0),QjetsPreclustering_);
59 
60  edm::Service<edm::RandomNumberGenerator> rng;
61  CLHEP::HepRandomEngine* engine = &rng->getEngine(iEvent.streamID());
62  qjetsAlgo_.SetRNEngine(engine);
63  // create probabilistic recusterings
64  for(unsigned int ii = 0 ; ii < (unsigned int) ntrial_ ; ii++){
65  //qjetsAlgo_.SetRandSeed(iEvent.id().event()*100 + (jetIt - jets->begin())*ntrial_ + ii );// set random seed for reprudcibility. We need a smarted scheme
66  fastjet::ClusterSequence qjet_seq(constits, qjet_def);
67  vector<fastjet::PseudoJet> inclusive_jets2 = sorted_by_pt(qjet_seq.inclusive_jets(cutoff_));
68  if (inclusive_jets2.size()>0){ // fill the massvalue only if the reclustering was successfull
69  qjetmass.push_back(inclusive_jets2[0].m());
70  }
71 
72  }//end loop on trials
73 
74  if (qjetmass.size()<1) {//protection against dummy case
75  QjetsVolatility.push_back(-1);
76  continue;
77  }
78 
79  double mean = std::accumulate( qjetmass.begin( ) , qjetmass.end( ) , 0 ) /qjetmass.size() ;
80  float totalsquared = 0.;
81  for (unsigned int i = 0; i < qjetmass.size(); i++){
82  totalsquared += (qjetmass[i] - mean)*(qjetmass[i] - mean) ;
83  }
84  float variance = sqrt( totalsquared/qjetmass.size() );
85 
86  QjetsVolatility.push_back(variance/mean);
87  }//end loop on jets
88 
89  auto outQJV = std::make_unique<edm::ValueMap<float>>();
90  edm::ValueMap<float>::Filler fillerQJV(*outQJV);
91  fillerQJV.insert(jets, QjetsVolatility.begin(), QjetsVolatility.end());
92  fillerQJV.fill();
93 
94  iEvent.put(std::move(outQJV),"QjetsVolatility");
95 }
reco::Candidate::pz
virtual double pz() const =0
z coordinate of momentum vector
edm::Event::put
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:122
QjetsPlugin::SetRNEngine
void SetRNEngine(CLHEP::HepRandomEngine *rnEngine)
Definition: QjetsPlugin.h:26
mps_fire.i
i
Definition: mps_fire.py:156
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:457
QjetsAdder::ntrial_
int ntrial_
Definition: QjetsAdder.h:40
reco::Jet
Base class for all types of Jets.
Definition: Jet.h:20
edm::RandomNumberGenerator::getEngine
virtual CLHEP::HepRandomEngine & getEngine(StreamID const &)=0
Use this engine in event methods.
QjetsAdder::jetRad_
double jetRad_
Definition: QjetsAdder.h:42
reco::Candidate::energy
virtual double energy() const =0
energy
reco::Candidate::py
virtual double py() const =0
y coordinate of momentum vector
QjetsAdder::qjetsAlgo_
QjetsPlugin qjetsAlgo_
Definition: QjetsAdder.h:39
QjetsAdder::QjetsPreclustering_
int QjetsPreclustering_
Definition: QjetsAdder.h:44
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
QjetsAdder::cutoff_
double cutoff_
Definition: QjetsAdder.h:41
fwrapper::jets
vector< PseudoJet > jets
Definition: fastjetfortran_madfks.cc:42
QjetsAdder::src_token_
edm::EDGetTokenT< edm::View< reco::Jet > > src_token_
Definition: QjetsAdder.h:38
cuy.ii
ii
Definition: cuy.py:588
gen::k
int k[5][pyjets_maxn]
Definition: Cascade2Hadronizer.cc:79
QjetsAdder::mJetAlgo_
std::string mJetAlgo_
Definition: QjetsAdder.h:43
edm::View::const_iterator
boost::indirect_iterator< typename seq_t::const_iterator > const_iterator
Definition: View.h:86
edm::Event::streamID
StreamID streamID() const
Definition: Event.h:81
reco::Candidate::px
virtual double px() const =0
x coordinate of momentum vector
eostools.move
def move(src, dest)
Definition: eostools.py:510

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.

Generated for CMSSW Reference Manual by doxygen 1.8.11