CMS 3D CMS Logo

List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes | Static Private Attributes
TauSpinnerCMS Class Reference

#include <TauSpinnerCMS.h>

Inheritance diagram for TauSpinnerCMS:
edm::one::EDProducer< edm::one::SharedResources > edm::one::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

void beginJob () final
 
void endJob () final
 
virtual void initialize ()
 
void produce (edm::Event &, const edm::EventSetup &) final
 
void setRandomEngine (CLHEP::HepRandomEngine *v)
 
 TauSpinnerCMS (const edm::ParameterSet &)
 
 ~TauSpinnerCMS () override
 
- Public Member Functions inherited from edm::one::EDProducer< edm::one::SharedResources >
 EDProducer ()=default
 
SerialTaskQueueglobalLuminosityBlocksQueue () final
 
SerialTaskQueueglobalRunsQueue () final
 
bool hasAbilityToProduceInLumis () const final
 
bool hasAbilityToProduceInRuns () const final
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
- Public Member Functions inherited from edm::one::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~EDProducerBase () 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 () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () 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
 
EDConsumerBaseoperator= (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)
 

Static Public Member Functions

static double flat ()
 
- Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 

Private Member Functions

void GetLastSelf (const reco::GenParticle *Particle)
 
void GetRecoDaughters (const reco::GenParticle *Particle, std::vector< TauSpinner::SimpleParticle > &daughters, int parentpdgid)
 
bool isFirst (const reco::GenParticle *Particle)
 
int readParticlesfromReco (edm::Event &e, TauSpinner::SimpleParticle &X, TauSpinner::SimpleParticle &tau, TauSpinner::SimpleParticle &tau2, std::vector< TauSpinner::SimpleParticle > &tau_daughters, std::vector< TauSpinner::SimpleParticle > &tau2_daughters)
 

Private Attributes

double CMSEnergy_
 
edm::EDGetTokenT< reco::GenParticleCollectionGenParticleCollectionToken_
 
edm::InputTag gensrc_
 
edm::EDGetTokenT< edm::HepMCProducthepmcCollectionToken_
 
int Ipol_
 
bool isLHPDFConfigured_
 
bool isReco_
 
bool isTauolaConfigured_
 
std::string LHAPDFname_
 
int MotherPDGID_
 
int nonSM2_
 
int nonSMN_
 
double roundOff_
 

Static Private Attributes

static bool fInitialized =false
 
static CLHEP::HepRandomEngine * fRandomEngine = 0
 
static bool isTauSpinnerConfigure =false
 

Additional Inherited Members

- Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase 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
 
- 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 49 of file TauSpinnerCMS.h.

Constructor & Destructor Documentation

TauSpinnerCMS::TauSpinnerCMS ( const edm::ParameterSet pset)
explicit

Definition at line 23 of file TauSpinnerCMS.cc.

References GenParticleCollectionToken_, gensrc_, hepmcCollectionToken_, isReco_, and edm::SharedResourceNames::kTauola.

23  :
24  EDProducer()
25  ,isReco_(pset.getParameter<bool>("isReco"))
26  ,isTauolaConfigured_(pset.getParameter<bool>("isTauolaConfigured" ))
27  ,isLHPDFConfigured_(pset.getParameter<bool>("isLHPDFConfigured" ))
28  ,LHAPDFname_(pset.getUntrackedParameter("LHAPDFname",(string)("MSTW2008nnlo90cl.LHgrid")))
29  ,CMSEnergy_(pset.getParameter<double>("CMSEnergy"))//GeV
30  ,gensrc_(pset.getParameter<edm::InputTag>("gensrc"))
31  ,MotherPDGID_(pset.getUntrackedParameter("MotherPDGID",(int)(-1)))
32  ,Ipol_(pset.getUntrackedParameter("Ipol",(int)(0)))
33  ,nonSM2_(pset.getUntrackedParameter("nonSM2",(int)(0)))
34  ,nonSMN_(pset.getUntrackedParameter("nonSMN",(int)(0)))
35  ,roundOff_(pset.getUntrackedParameter("roundOff",(double)(0.01)))
36 {
37  //usesResource(edm::uniqueSharedResourceName());
39 
40  produces<bool>("TauSpinnerWTisValid").setBranchAlias("TauSpinnerWTisValid");
41  produces<double>("TauSpinnerWT").setBranchAlias("TauSpinnerWT");
42  produces<double>("TauSpinnerWTFlip").setBranchAlias("TauSpinnerWTFlip");
43  produces<double>("TauSpinnerWThplus").setBranchAlias("TauSpinnerWThplus");
44  produces<double>("TauSpinnerWThminus").setBranchAlias("TauSpinnerWThminus");
45 
46  if(isReco_){
47  GenParticleCollectionToken_=consumes<reco::GenParticleCollection>(gensrc_);
48  }
49  else{
50  hepmcCollectionToken_=consumes<HepMCProduct>(gensrc_);
51  }
52 }
T getParameter(std::string const &) const
T getUntrackedParameter(std::string const &, T const &) const
edm::EDGetTokenT< edm::HepMCProduct > hepmcCollectionToken_
Definition: TauSpinnerCMS.h:86
static const std::string kTauola
double roundOff_
Definition: TauSpinnerCMS.h:83
edm::EDGetTokenT< reco::GenParticleCollection > GenParticleCollectionToken_
Definition: TauSpinnerCMS.h:87
bool isTauolaConfigured_
Definition: TauSpinnerCMS.h:63
std::string LHAPDFname_
Definition: TauSpinnerCMS.h:65
edm::InputTag gensrc_
Definition: TauSpinnerCMS.h:67
double CMSEnergy_
Definition: TauSpinnerCMS.h:66
bool isLHPDFConfigured_
Definition: TauSpinnerCMS.h:64
TauSpinnerCMS::~TauSpinnerCMS ( )
inlineoverride

Definition at line 52 of file TauSpinnerCMS.h.

References beginJob(), endJob(), flat(), and produce().

52 {}; // no need to delete ROOT stuff

Member Function Documentation

void TauSpinnerCMS::beginJob ( void  )
finalvirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 54 of file TauSpinnerCMS.cc.

Referenced by ~TauSpinnerCMS().

54 {};
void TauSpinnerCMS::endJob ( void  )
finalvirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 55 of file TauSpinnerCMS.cc.

Referenced by ~TauSpinnerCMS().

55 {};
double TauSpinnerCMS::flat ( void  )
static

Definition at line 263 of file TauSpinnerCMS.cc.

References DEFINE_FWK_MODULE, Exception, and fRandomEngine.

Referenced by initialize(), produce(), and ~TauSpinnerCMS().

264 {
265  if ( !fRandomEngine ) {
266  throw cms::Exception("LogicError")
267  << "TauSpinnerCMS::flat: Attempt to generate random number when engine pointer is null\n"
268  << "This might mean that the code was modified to generate a random number outside the\n"
269  << "event and beginLuminosityBlock methods, which is not allowed.\n";
270  }
271  return fRandomEngine->flat();
272 }
static CLHEP::HepRandomEngine * fRandomEngine
Definition: TauSpinnerCMS.h:85
void TauSpinnerCMS::GetLastSelf ( const reco::GenParticle Particle)
private

Definition at line 231 of file TauSpinnerCMS.cc.

References reco::CompositeRefCandidateT< D >::daughter(), mps_fire::i, reco::CompositeRefCandidateT< D >::numberOfDaughters(), and reco::LeafCandidate::pdgId().

Referenced by readParticlesfromReco().

231  {
232  for (unsigned int i=0; i< Particle->numberOfDaughters(); i++){
233  const reco::GenParticle *dau=static_cast<const reco::GenParticle*>(Particle->daughter(i));
234  if(Particle->pdgId()==dau->pdgId()){
235  Particle=dau;
236  GetLastSelf(Particle);
237  }
238  }
239 }
int pdgId() const final
PDG identifier.
const Candidate * daughter(size_type) const override
return daughter at a given position, i = 0, ... numberOfDaughters() - 1 (read only mode) ...
void GetLastSelf(const reco::GenParticle *Particle)
size_t numberOfDaughters() const override
number of daughters
void TauSpinnerCMS::GetRecoDaughters ( const reco::GenParticle Particle,
std::vector< TauSpinner::SimpleParticle > &  daughters,
int  parentpdgid 
)
private

Definition at line 251 of file TauSpinnerCMS.cc.

References funct::abs(), reco::CompositeRefCandidateT< D >::daughter(), mps_fire::i, reco::CompositeRefCandidateT< D >::numberOfDaughters(), reco::LeafCandidate::p4(), and reco::LeafCandidate::pdgId().

Referenced by readParticlesfromReco().

251  {
252  if( Particle->numberOfDaughters()==0 || abs(Particle->pdgId())==111){
253  SimpleParticle tp(Particle->p4().Px(), Particle->p4().Py(), Particle->p4().Pz(), Particle->p4().E(), Particle->pdgId());
254  daughters.push_back(tp);
255  return;
256  }
257  for (unsigned int i=0; i< Particle->numberOfDaughters(); i++){
258  const reco::Candidate *dau=Particle->daughter(i);
259  GetRecoDaughters(static_cast<const reco::GenParticle*>(dau),daughters,Particle->pdgId());
260  }
261 }
int pdgId() const final
PDG identifier.
const Candidate * daughter(size_type) const override
return daughter at a given position, i = 0, ... numberOfDaughters() - 1 (read only mode) ...
void GetRecoDaughters(const reco::GenParticle *Particle, std::vector< TauSpinner::SimpleParticle > &daughters, int parentpdgid)
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
const LorentzVector & p4() const final
four-momentum Lorentz vector
Definition: LeafCandidate.h:99
size_t numberOfDaughters() const override
number of daughters
void TauSpinnerCMS::initialize ( )
virtual

Definition at line 57 of file TauSpinnerCMS.cc.

References CMSEnergy_, fInitialized, flat(), initialize(), Ipol_, isLHPDFConfigured_, isTauolaConfigured_, isTauSpinnerConfigure, LHAPDFname_, nonSM2_, and nonSMN_.

Referenced by produce(), and setRandomEngine().

57  {
58  // Now for Tauola and TauSpinner
59  Tauolapp::Tauola::setRandomGenerator(TauSpinnerCMS::flat);
62  }
63  if(!isLHPDFConfigured_){
64  LHAPDF::initPDFSetByName(LHAPDFname_);
65  }
68  bool Ipp = true; // for pp collisions
69  // Initialize TauSpinner
70  //Ipol - polarization of input sample
71  //nonSM2 - nonstandard model calculations
72  //nonSMN
73  TauSpinner::initialize_spinner(Ipp,Ipol_,nonSM2_,nonSMN_,CMSEnergy_);
74  }
75  fInitialized=true;
76 }
static AlgebraicMatrix initialize()
static double flat()
static bool isTauSpinnerConfigure
Definition: TauSpinnerCMS.h:69
bool isTauolaConfigured_
Definition: TauSpinnerCMS.h:63
std::string LHAPDFname_
Definition: TauSpinnerCMS.h:65
static bool fInitialized
Definition: TauSpinnerCMS.h:88
double CMSEnergy_
Definition: TauSpinnerCMS.h:66
bool isLHPDFConfigured_
Definition: TauSpinnerCMS.h:64
bool TauSpinnerCMS::isFirst ( const reco::GenParticle Particle)
private

Definition at line 241 of file TauSpinnerCMS.cc.

References mps_fire::i, reco::CompositeRefCandidateT< D >::mother(), reco::CompositeRefCandidateT< D >::numberOfMothers(), and reco::LeafCandidate::pdgId().

Referenced by readParticlesfromReco().

241  {
242  for (unsigned int i=0; i< Particle->numberOfMothers(); i++){
243  const reco::GenParticle *moth=static_cast<const reco::GenParticle*>(Particle->mother(i));
244  if(Particle->pdgId()==moth->pdgId()){
245  return false;
246  }
247  }
248  return true;
249 }
int pdgId() const final
PDG identifier.
size_t numberOfMothers() const override
number of mothers
const Candidate * mother(size_type=0) const override
return mother at a given position, i = 0, ... numberOfMothers() - 1 (read only mode) ...
void TauSpinnerCMS::produce ( edm::Event e,
const edm::EventSetup iSetup 
)
final

Definition at line 79 of file TauSpinnerCMS.cc.

References funct::abs(), gather_cfg::cout, fInitialized, flat(), edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), hepmcCollectionToken_, mps_fire::i, initialize(), isReco_, MotherPDGID_, eostools::move(), edm::Event::put(), readParticlesFromHepMC(), readParticlesfromReco(), reco::return(), roundOff_, trackingPlots::stat, edm::Event::streamID(), metsig::tau, jets_cff::tau2, and X.

Referenced by ~TauSpinnerCMS().

79  {
80  RandomEngineSentry<TauSpinnerCMS> randomEngineSentry(this, e.streamID());
81  if(!fInitialized) initialize();
82 
83  Tauolapp::Tauola::setRandomGenerator(TauSpinnerCMS::flat); // rest tauola++ random number incase other modules use tauola++
84  Tauolapp::jaki_.ktom = 1; // rest for when you run after tauola
85  double WT=1.0;
86  double WTFlip=1.0;
87  double polSM=-999; //range [-1,1]
88  SimpleParticle X, tau, tau2;
89  std::vector<SimpleParticle> tau_daughters, tau_daughters2;
90  int stat(0);
91  if(isReco_){
92  stat=readParticlesfromReco(e,X,tau,tau2,tau_daughters,tau_daughters2);
93  }
94  else{
97  //Get EVENT
98  HepMC::GenEvent *Evt = new HepMC::GenEvent(*(evt->GetEvent()));
99  stat=readParticlesFromHepMC(Evt,X,tau,tau2,tau_daughters,tau_daughters2);
100  }
101  if(MotherPDGID_<0 || abs(X.pdgid())==MotherPDGID_){
102  if(stat!=1){
103  // Determine the weight
104  if( abs(X.pdgid())==24 || abs(X.pdgid())==37 ){
105  TLorentzVector tau_1r(0,0,0,0);
106  TLorentzVector tau_1(tau.px(),tau.py(),tau.pz(),tau.e());
107  for(unsigned int i=0; i<tau_daughters.size();i++){
108  tau_1r+=TLorentzVector(tau_daughters.at(i).px(),tau_daughters.at(i).py(),tau_daughters.at(i).pz(),tau_daughters.at(i).e());
109  }
110  if(fabs(tau_1r.M()-tau_1.M())<roundOff_){
111  WT = TauSpinner::calculateWeightFromParticlesWorHpn(X, tau, tau2, tau_daughters); // note that tau2 is tau neutrino
112  polSM=getTauSpin();
113  WTFlip=(2.0-WT)/WT;
114  }
115  }
116  else if( X.pdgid()==25 || X.pdgid()==36 || X.pdgid()==22 || X.pdgid()==23 ){
117  TLorentzVector tau_1r(0,0,0,0), tau_2r(0,0,0,0);
118  TLorentzVector tau_1(tau.px(),tau.py(),tau.pz(),tau.e()), tau_2(tau2.px(),tau2.py(),tau2.pz(),tau2.e());
119  for(unsigned int i=0; i<tau_daughters.size();i++){
120  tau_1r+=TLorentzVector(tau_daughters.at(i).px(),tau_daughters.at(i).py(),tau_daughters.at(i).pz(),tau_daughters.at(i).e());
121  }
122  for(unsigned int i=0; i<tau_daughters2.size();i++){
123  tau_2r+=TLorentzVector(tau_daughters2.at(i).px(),tau_daughters2.at(i).py(),tau_daughters2.at(i).pz(),tau_daughters2.at(i).e());
124  }
125 
126  if(fabs(tau_1r.M()-tau_1.M())<roundOff_ && fabs(tau_2r.M()-tau_2.M())<roundOff_){
127  WT = TauSpinner::calculateWeightFromParticlesH(X, tau, tau2, tau_daughters,tau_daughters2);
128  //std::cout << "WT " << WT << std::endl;
129  polSM=getTauSpin();
130  if(X.pdgid()==25 || X.pdgid()==22 || X.pdgid()==23 ){
131  if(X.pdgid()==25) X.setPdgid(23);
132  if( X.pdgid()==22 || X.pdgid()==23 ) X.setPdgid(25);
133 
134  double WTother=TauSpinner::calculateWeightFromParticlesH(X, tau, tau2, tau_daughters,tau_daughters2);
135  WTFlip=WTother/WT;
136  }
137  }
138  }
139  else{
140  cout<<"TauSpinner: WARNING: Unexpected PDG for tau mother: "<<X.pdgid()<<endl;
141  }
142  }
143  }
144  bool isValid=true;
145  if(!(0<=WT && WT<10)){isValid=false; WT=1.0; WTFlip=1.0;}
146  std::unique_ptr<bool> TauSpinnerWeightisValid(new bool);
147  *TauSpinnerWeightisValid =isValid;
148  e.put(std::move(TauSpinnerWeightisValid),"TauSpinnerWTisValid");
149 
150  // regular weight
151  std::unique_ptr<double> TauSpinnerWeight(new double);
152  *TauSpinnerWeight =WT;
153  e.put(std::move(TauSpinnerWeight),"TauSpinnerWT");
154 
155  // flipped weight (ie Z->H or H->Z)
156  std::unique_ptr<double> TauSpinnerWeightFlip(new double);
157  *TauSpinnerWeightFlip =WTFlip;
158  e.put(std::move(TauSpinnerWeightFlip),"TauSpinnerWTFlip");
159 
160  // h+ polarization
161  double WThplus=WT;
162  if(polSM<0.0 && polSM!=-999 && isValid) WThplus=0;
163  std::unique_ptr<double> TauSpinnerWeighthplus(new double);
164  *TauSpinnerWeighthplus = WThplus;
165  e.put(std::move(TauSpinnerWeighthplus),"TauSpinnerWThplus");
166 
167  // h- polarization
168  double WThminus=WT;
169  if(polSM>0.0&& polSM!=-999 && isValid) WThminus=0;
170  std::unique_ptr<double> TauSpinnerWeighthminus(new double);
171  *TauSpinnerWeighthminus = WThminus;
172  e.put(std::move(TauSpinnerWeighthminus),"TauSpinnerWThminus");
173  return ;
174 }
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:136
edm::EDGetTokenT< edm::HepMCProduct > hepmcCollectionToken_
Definition: TauSpinnerCMS.h:86
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:519
virtual void initialize()
static double flat()
#define X(str)
Definition: MuonsGrabber.cc:48
double roundOff_
Definition: TauSpinnerCMS.h:83
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
int readParticlesfromReco(edm::Event &e, TauSpinner::SimpleParticle &X, TauSpinner::SimpleParticle &tau, TauSpinner::SimpleParticle &tau2, std::vector< TauSpinner::SimpleParticle > &tau_daughters, std::vector< TauSpinner::SimpleParticle > &tau2_daughters)
int readParticlesFromHepMC(const HepMC::GenEvent *event, TauSpinner::SimpleParticle &X, TauSpinner::SimpleParticle &tau, TauSpinner::SimpleParticle &tau2, std::vector< TauSpinner::SimpleParticle > &tau_daughters, std::vector< TauSpinner::SimpleParticle > &tau2_daughters)
const HepMC::GenEvent * GetEvent() const
Definition: HepMCProduct.h:38
static bool fInitialized
Definition: TauSpinnerCMS.h:88
return(e1-e2)*(e1-e2)+dp *dp
StreamID streamID() const
Definition: Event.h:95
def move(src, dest)
Definition: eostools.py:510
int TauSpinnerCMS::readParticlesfromReco ( edm::Event e,
TauSpinner::SimpleParticle &  X,
TauSpinner::SimpleParticle &  tau,
TauSpinner::SimpleParticle &  tau2,
std::vector< TauSpinner::SimpleParticle > &  tau_daughters,
std::vector< TauSpinner::SimpleParticle > &  tau2_daughters 
)
private

Definition at line 176 of file TauSpinnerCMS.cc.

References funct::abs(), reco::Candidate::daughter(), GenParticleCollectionToken_, GenHFHadronMatcher_cfi::genParticles, edm::Event::getByToken(), GetLastSelf(), GetRecoDaughters(), mps_fire::i, isFirst(), reco::LeafCandidate::p4(), BPhysicsValidation_cfi::pdgid, reco::Candidate::pdgId(), and reco::LeafCandidate::pdgId().

Referenced by produce().

177  {
179  e.getByToken(GenParticleCollectionToken_, genParticles);
180  for(reco::GenParticleCollection::const_iterator itr = genParticles->begin(); itr!= genParticles->end(); ++itr){
181  int pdgid=abs(itr->pdgId());
182  if(pdgid==24 || pdgid==37 || pdgid ==25 || pdgid==36 || pdgid==22 || pdgid==23 ){
183  const reco::GenParticle *hx=&(*itr);
184  if(!isFirst(hx)) continue;
185  GetLastSelf(hx);
186  const reco::GenParticle *recotau1=nullptr;
187  const reco::GenParticle *recotau2=nullptr;
188  unsigned int ntau(0),ntauornu(0);
189  for(unsigned int i=0; i<itr->numberOfDaughters(); i++){
190  const reco::Candidate *dau=itr->daughter(i);
191  if(abs(dau->pdgId())!=pdgid){
192  if(abs(dau->pdgId())==15 || abs(dau->pdgId())==16){
193  if(ntau==0 && abs(dau->pdgId())==15){
194  recotau1=static_cast<const reco::GenParticle*>(dau);
195  GetLastSelf(recotau1);
196  ntau++;
197  }
198  else if((ntau==1 && abs(dau->pdgId())==15) || abs(dau->pdgId())==16){
199  recotau2=static_cast<const reco::GenParticle*>(dau);
200  if(abs(dau->pdgId())==15){ntau++;GetLastSelf(recotau2);}
201  }
202  ntauornu++;
203  }
204  }
205  }
206  if((ntau==2 && ntauornu==2) || (ntau==1 && ntauornu==2)){
207  X.setPx(itr->p4().Px());
208  X.setPy(itr->p4().Py());
209  X.setPz(itr->p4().Pz());
210  X.setE (itr->p4().E());
211  X.setPdgid(itr->pdgId());
212  tau.setPx(recotau1->p4().Px());
213  tau.setPy(recotau1->p4().Py());
214  tau.setPz(recotau1->p4().Pz());
215  tau.setE (recotau1->p4().E());
216  tau.setPdgid(recotau1->pdgId());
217  GetRecoDaughters(recotau1,tau_daughters,recotau1->pdgId());
218  tau2.setPx(recotau2->p4().Px());
219  tau2.setPy(recotau2->p4().Py());
220  tau2.setPz(recotau2->p4().Pz());
221  tau2.setE (recotau2->p4().E());
222  tau2.setPdgid(recotau2->pdgId());
223  if(ntau==2)GetRecoDaughters(recotau2,tau2_daughters,recotau2->pdgId());
224  return 0;
225  }
226  }
227  }
228  return 1;
229 }
int pdgId() const final
PDG identifier.
void GetLastSelf(const reco::GenParticle *Particle)
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:519
virtual const Candidate * daughter(size_type i) const =0
return daughter at a given position, i = 0, ... numberOfDaughters() - 1 (read only mode) ...
bool isFirst(const reco::GenParticle *Particle)
#define X(str)
Definition: MuonsGrabber.cc:48
edm::EDGetTokenT< reco::GenParticleCollection > GenParticleCollectionToken_
Definition: TauSpinnerCMS.h:87
void GetRecoDaughters(const reco::GenParticle *Particle, std::vector< TauSpinner::SimpleParticle > &daughters, int parentpdgid)
virtual int pdgId() const =0
PDG identifier.
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
const LorentzVector & p4() const final
four-momentum Lorentz vector
Definition: LeafCandidate.h:99
void TauSpinnerCMS::setRandomEngine ( CLHEP::HepRandomEngine *  v)
inline

Definition at line 58 of file TauSpinnerCMS.h.

References fRandomEngine, initialize(), and findQualityFiles::v.

58 { fRandomEngine = v; }
static CLHEP::HepRandomEngine * fRandomEngine
Definition: TauSpinnerCMS.h:85

Member Data Documentation

double TauSpinnerCMS::CMSEnergy_
private

Definition at line 66 of file TauSpinnerCMS.h.

Referenced by initialize().

bool TauSpinnerCMS::fInitialized =false
staticprivate

Definition at line 88 of file TauSpinnerCMS.h.

Referenced by initialize(), and produce().

CLHEP::HepRandomEngine * TauSpinnerCMS::fRandomEngine = 0
staticprivate

Definition at line 85 of file TauSpinnerCMS.h.

Referenced by flat(), and setRandomEngine().

edm::EDGetTokenT<reco::GenParticleCollection> TauSpinnerCMS::GenParticleCollectionToken_
private

Definition at line 87 of file TauSpinnerCMS.h.

Referenced by readParticlesfromReco(), and TauSpinnerCMS().

edm::InputTag TauSpinnerCMS::gensrc_
private

Definition at line 67 of file TauSpinnerCMS.h.

Referenced by TauSpinnerCMS().

edm::EDGetTokenT<edm::HepMCProduct> TauSpinnerCMS::hepmcCollectionToken_
private

Definition at line 86 of file TauSpinnerCMS.h.

Referenced by produce(), and TauSpinnerCMS().

int TauSpinnerCMS::Ipol_
private

Definition at line 68 of file TauSpinnerCMS.h.

Referenced by initialize().

bool TauSpinnerCMS::isLHPDFConfigured_
private

Definition at line 64 of file TauSpinnerCMS.h.

Referenced by initialize().

bool TauSpinnerCMS::isReco_
private

Definition at line 62 of file TauSpinnerCMS.h.

Referenced by produce(), and TauSpinnerCMS().

bool TauSpinnerCMS::isTauolaConfigured_
private

Definition at line 63 of file TauSpinnerCMS.h.

Referenced by initialize().

bool TauSpinnerCMS::isTauSpinnerConfigure =false
staticprivate

Definition at line 69 of file TauSpinnerCMS.h.

Referenced by initialize().

std::string TauSpinnerCMS::LHAPDFname_
private

Definition at line 65 of file TauSpinnerCMS.h.

Referenced by initialize().

int TauSpinnerCMS::MotherPDGID_
private

Definition at line 68 of file TauSpinnerCMS.h.

Referenced by produce().

int TauSpinnerCMS::nonSM2_
private

Definition at line 68 of file TauSpinnerCMS.h.

Referenced by initialize().

int TauSpinnerCMS::nonSMN_
private

Definition at line 68 of file TauSpinnerCMS.h.

Referenced by initialize().

double TauSpinnerCMS::roundOff_
private

Definition at line 83 of file TauSpinnerCMS.h.

Referenced by produce().