TauValidation Class Reference

#include <TauValidation.h>

Inheritance diagram for TauValidation:

Public Types
enum	{   undetermined, electron, muon, pi,   rho, a1, K, Kstar,   pi1pi0, pinpi0, tripi, tripinpi0,   stable }
enum	{   other, B, D, gamma,   Z, W, HSM, H0,   A0, Hpm }
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...>	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T...>	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
virtual void	analyze (edm::Event const &, edm::EventSetup const &) override
virtual void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &) override
virtual void	dqmBeginRun (const edm::Run &r, const edm::EventSetup &c) override
	TauValidation (const edm::ParameterSet &)
virtual	~TauValidation ()
Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final
virtual void	beginStream (edm::StreamID id) final
	DQMEDAnalyzer (void)
virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
uint32_t	streamId () const
Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default
Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDAnalyzerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
void	countParticles (const reco::GenParticle *p, int &allCount, int &eCount, int &muCount, int &pi0Count, int &piCount, int &rhoCount, int &a1Count, int &KCount, int &KstarCount)
void	findFSRandBrem (const reco::GenParticle *p, bool doBrem, std::vector< const reco::GenParticle * > &ListofFSR, std::vector< const reco::GenParticle * > &ListofBrem)
int	findMother (const reco::GenParticle *)
void	FindPhotosFSR (const reco::GenParticle *p, std::vector< const reco::GenParticle * > &ListofFSR, double &BosonScale)
void	findTauList (const reco::GenParticle *tau, std::vector< const reco::GenParticle * > &TauList)
const reco::GenParticle *	GetMother (const reco::GenParticle *tau)
const std::vector< const reco::GenParticle * >	GetMothers (const reco::GenParticle *boson)
bool	isLastTauinChain (const reco::GenParticle *tau)
double	leadingPionMomentum (const reco::GenParticle *, double weight)
TLorentzVector	leadingPionP4 (const reco::GenParticle *)
TLorentzVector	motherP4 (const reco::GenParticle *)
void	photons (const reco::GenParticle *, double weight)
void	spinEffectsWHpm (const reco::GenParticle *, int, int, std::vector< const reco::GenParticle * > &part, double weight)
void	spinEffectsZH (const reco::GenParticle *boson, double weight)
int	tauDecayChannel (const reco::GenParticle *tau, int jak_id, unsigned int TauBitMask, double weight)
int	tauMother (const reco::GenParticle *, double weight)
int	tauProngs (const reco::GenParticle *, double weight)
double	visibleTauEnergy (const reco::GenParticle *)
double	Zstoa (double zs)

Private Attributes
MonitorElement *	DecayLength
edm::ESHandle < HepPDT::ParticleDataTable >	fPDGTable
	PDT table. More...
edm::InputTag	genparticleCollection_
edm::EDGetTokenT < reco::GenParticleCollection >	genparticleCollectionToken_
MonitorElement *	LifeTime
MonitorElement *	MODEID
std::vector< std::vector < MonitorElement * > >	MODEInvMass
unsigned int	NMODEID
MonitorElement *	nPrimeTaus
MonitorElement *	nTaus
MonitorElement *	TauBremPhotonsN
MonitorElement *	TauBremPhotonsPt
MonitorElement *	TauBremPhotonsPtSum
MonitorElement *	TauDecayChannels
MonitorElement *	TauEta
MonitorElement *	TauFSRPhotonsN
MonitorElement *	TauFSRPhotonsPt
MonitorElement *	TauFSRPhotonsPtSum
MonitorElement *	TauMothers
MonitorElement *	TauPhi
MonitorElement *	TauProngs
MonitorElement *	TauPt
MonitorElement *	TauSpinEffectsH_eX
MonitorElement *	TauSpinEffectsH_muX
MonitorElement *	TauSpinEffectsH_MVis
MonitorElement *	TauSpinEffectsH_pipiAcollinearity
MonitorElement *	TauSpinEffectsH_pipiAcollinearityzoom
MonitorElement *	TauSpinEffectsH_pipiAcoplanarity
MonitorElement *	TauSpinEffectsH_rhorhoAcoplanarityminus
MonitorElement *	TauSpinEffectsH_rhorhoAcoplanarityplus
MonitorElement *	TauSpinEffectsH_X
MonitorElement *	TauSpinEffectsH_Xb
MonitorElement *	TauSpinEffectsH_Xf
MonitorElement *	TauSpinEffectsH_Zs
MonitorElement *	TauSpinEffectsHpm_eX
MonitorElement *	TauSpinEffectsHpm_muX
MonitorElement *	TauSpinEffectsHpm_UpsilonA1
MonitorElement *	TauSpinEffectsHpm_UpsilonRho
MonitorElement *	TauSpinEffectsHpm_X
MonitorElement *	TauSpinEffectsW_eX
MonitorElement *	TauSpinEffectsW_muX
MonitorElement *	TauSpinEffectsW_UpsilonA1
MonitorElement *	TauSpinEffectsW_UpsilonRho
MonitorElement *	TauSpinEffectsW_X
MonitorElement *	TauSpinEffectsZ_eX
MonitorElement *	TauSpinEffectsZ_muX
MonitorElement *	TauSpinEffectsZ_MVis
MonitorElement *	TauSpinEffectsZ_X
MonitorElement *	TauSpinEffectsZ_X100to120
MonitorElement *	TauSpinEffectsZ_X120UP
MonitorElement *	TauSpinEffectsZ_X50to75
MonitorElement *	TauSpinEffectsZ_X75to88
MonitorElement *	TauSpinEffectsZ_X88to100
MonitorElement *	TauSpinEffectsZ_Xb
MonitorElement *	TauSpinEffectsZ_Xf
MonitorElement *	TauSpinEffectsZ_Zs
int	zsbins
double	zsmax
double	zsmin

Additional Inherited Members
Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
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 26 of file TauValidation.h.

Member Enumeration Documentation

anonymous enum

Enumerator
undetermined
electron
muon
pi
rho
a1
K
Kstar
pi1pi0
pinpi0
tripi
tripinpi0
stable

Definition at line 30 of file TauValidation.h.

              {undetermined,
               electron,
               muon,
               pi,
               rho,
           a1,
               K,
           Kstar,
           pi1pi0,
               pinpi0,
               tripi,
               tripinpi0,
           stable};

anonymous enum

Enumerator
other
B
D
gamma
Z
W
HSM
H0
A0
Hpm

Definition at line 44 of file TauValidation.h.

          {other,
           B,
           D,
           gamma,
           Z,
           W,
           HSM,
           H0,
           A0,
           Hpm};

Constructor & Destructor Documentation

TauValidation::TauValidation ( const edm::ParameterSet &  iPSet )

explicit

Definition at line 18 of file TauValidation.cc.

References genparticleCollection_, and genparticleCollectionToken_.

                                                        : 
  //  wmanager_(iPSet,consumesCollector())
  genparticleCollection_(iPSet.getParameter<edm::InputTag>("genparticleCollection"))
  ,NMODEID(TauDecay::NMODEID-1)// fortran to C++ index
  ,zsbins(20)
  ,zsmin(-0.5)
  ,zsmax(0.5)
{    
  genparticleCollectionToken_=consumes<reco::GenParticleCollection>(genparticleCollection_);
}

TauValidation::~TauValidation ( )

virtual

Definition at line 29 of file TauValidation.cc.

{}

Member Function Documentation

void TauValidation::analyze ( edm::Event const &  iEvent, edm::EventSetup const &  iSetup  )

overridevirtual

Gathering the reco::GenParticleCollection information

Implements edm::stream::EDAnalyzerBase.

Definition at line 156 of file TauValidation.cc.

References funct::abs(), TauDecay_GenParticle::AnalyzeTau(), beta, EnergyCorrector::c, DecayLength, MonitorElement::Fill(), genparticleCollectionToken_, genParticleCandidates2GenParticles_cfi::genParticles, TauDecay_GenParticle::Get_TauDecayProducts(), edm::Event::getByToken(), PdtPdgMini::Higgs0, i, isLastTauinChain(), TauDecay::isTauFinalStateParticle(), PdtPdgMini::K_plus, LifeTime, TauDecay::MODE_3PI, TauDecay::MODE_KPIK, TauDecay::MODE_KPIPI, TauDecay::MODE_PI2PI0, MODEID, MODEInvMass, NMODEID, nPrimeTaus, TauDecay::nProng(), nTaus, PdtPdgMini::nu_e, PdtPdgMini::nu_mu, PdtPdgMini::nu_tau, photons(), PdtPdgMini::pi_plus, spinEffectsWHpm(), spinEffectsZH(), tauDecayChannel(), TauEta, tauMother(), TauPhi, TauProngs, TauPt, histoStyle::weight, and PdtPdgMini::Z0.

                                                                             { 
  edm::Handle<reco::GenParticleCollection> genParticles;
  iEvent.getByToken(genparticleCollectionToken_, genParticles );
  
  double weight = 1.0;//=   wmanager_.weight(iEvent);
  // find taus
  for (reco::GenParticleCollection::const_iterator iter = genParticles->begin(); iter != genParticles->end(); ++iter) {
    if(abs(iter->pdgId())==PdtPdgMini::Z0 || abs(iter->pdgId())==PdtPdgMini::Higgs0){
      spinEffectsZH(&(*iter),weight);
    }
    if(abs(iter->pdgId())==15){
      if(isLastTauinChain(&(*iter))){
    nTaus->Fill(0.5,weight);
    int mother  = tauMother(&(*iter),weight);
    if(mother>-1){ // exclude B, D and other non-signal decay modes
      nPrimeTaus->Fill(0.5,weight);
      TauPt->Fill(iter->pt(),weight);
      TauEta->Fill(iter->eta(),weight);
      TauPhi->Fill(iter->phi(),weight);
      photons(&(*iter),weight);
      // Adding MODEID and Mass information
      TauDecay_GenParticle TD;
      unsigned int jak_id, TauBitMask;
      if(TD.AnalyzeTau(&(*iter),jak_id,TauBitMask,false,false)){
        MODEID->Fill(jak_id,weight);
        TauProngs->Fill(TD.nProng(TauBitMask),weight);
        tauDecayChannel(&(*iter),jak_id,TauBitMask,weight);
        if(jak_id<=NMODEID){
          int tcharge=iter->pdgId()/abs(iter->pdgId());
          std::vector<const reco::GenParticle*> part=TD.Get_TauDecayProducts();
          spinEffectsWHpm(&(*iter),mother,jak_id,part,weight);
          TLorentzVector LVQ(0,0,0,0);
          TLorentzVector LVS12(0,0,0,0);
          TLorentzVector LVS13(0,0,0,0);
          TLorentzVector LVS23(0,0,0,0);
          bool haspart1=false;
          TVector3 PV,SV;
          bool hasDL(false);
          for(unsigned int i=0;i<part.size();i++){
        if(abs(part.at(i)->pdgId())!=PdtPdgMini::nu_tau && TD.isTauFinalStateParticle(part.at(i)->pdgId()) && !hasDL){
          hasDL=true;
          TLorentzVector tlv(iter->px(),iter->py(),iter->pz(),iter->energy());
          PV=TVector3(iter->vx(),iter->vy(),iter->vz());
          SV=TVector3(part.at(i)->vx(),part.at(i)->vy(),part.at(i)->vz());
          TVector3 DL=SV-PV;
          DecayLength->Fill(DL.Dot(tlv.Vect())/tlv.P(),weight);
          double c(2.99792458E8),Ltau(DL.Mag()/100)/*cm->m*/,beta(iter->p()/iter->mass());
          LifeTime->Fill( Ltau/(c*beta),weight);
        }
        
        if(TD.isTauFinalStateParticle(part.at(i)->pdgId()) && abs(part.at(i)->pdgId())!=PdtPdgMini::nu_e &&
           abs(part.at(i)->pdgId())!=PdtPdgMini::nu_mu && abs(part.at(i)->pdgId())!=PdtPdgMini::nu_tau ){
          TLorentzVector LV(part.at(i)->px(),part.at(i)->py(),part.at(i)->pz(),part.at(i)->energy());
          LVQ+=LV;
          if(jak_id==TauDecay::MODE_3PI || jak_id==TauDecay::MODE_PI2PI0 ||
             jak_id==TauDecay::MODE_KPIK ||
             jak_id==TauDecay::MODE_KPIPI
             ){
            if((tcharge==part.at(i)->pdgId()/abs(part.at(i)->pdgId()) && TD.nProng(TauBitMask)==3) || ((jak_id==TauDecay::MODE_3PI || jak_id==TauDecay::MODE_PI2PI0) && TD.nProng(TauBitMask)==1 && abs(part.at(i)->pdgId())==PdtPdgMini::pi_plus) ){
              LVS13+=LV;
              LVS23+=LV;
            }
            else{
              LVS12+=LV;
              if(!haspart1 && ((jak_id==TauDecay::MODE_3PI || jak_id==TauDecay::MODE_PI2PI0)  || ((jak_id!=TauDecay::MODE_3PI || jak_id==TauDecay::MODE_PI2PI0) && abs(part.at(i)->pdgId())==PdtPdgMini::K_plus) )){
            LVS13+=LV;
            haspart1=true;
              }
              else{
            LVS23+=LV;
              }
            }
          }
        }
          }
          part.clear();
          MODEInvMass.at(jak_id).at(0)->Fill(LVQ.M(),weight);
          if(jak_id==TauDecay::MODE_3PI || jak_id==TauDecay::MODE_PI2PI0 ||
         jak_id==TauDecay::MODE_KPIK ||
         jak_id==TauDecay::MODE_KPIPI
         ){
        MODEInvMass.at(jak_id).at(1)->Fill(LVS13.M(),weight);
        MODEInvMass.at(jak_id).at(2)->Fill(LVS23.M(),weight);
        MODEInvMass.at(jak_id).at(3)->Fill(LVS12.M(),weight);
          }
        }
      }
      else{
        MODEID->Fill(jak_id,weight);  
      }
    }
      }
    }
  }
}//analyze

void TauValidation::bookHistograms ( DQMStore::IBooker &  i, edm::Run const &  , edm::EventSetup const &    )

overridevirtual

Setting the DQM top directories

Implements DQMEDAnalyzer.

Definition at line 35 of file TauValidation.cc.

References A0, a1, DQMHelper::book1dHisto(), D, DecayLength, TauDecay::DecayMode(), electron, gamma, H0, Hpm, HSM, j, K, Kstar, LifeTime, TauDecay::MODE_3PI, TauDecay::MODE_KPIK, TauDecay::MODE_KPIPI, TauDecay::MODE_PI2PI0, MODEID, MODEInvMass, muon, NMODEID, nPrimeTaus, nTaus, other, Pi, pi, pi1pi0, pinpi0, rho, MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), stable, AlCaHLTBitMon_QueryRunRegistry::string, TauBremPhotonsN, TauBremPhotonsPt, TauBremPhotonsPtSum, TauDecayChannels, TauEta, TauFSRPhotonsN, TauFSRPhotonsPt, TauFSRPhotonsPtSum, TauMothers, TauPhi, TauProngs, TauPt, TauSpinEffectsH_eX, TauSpinEffectsH_muX, TauSpinEffectsH_MVis, TauSpinEffectsH_pipiAcollinearity, TauSpinEffectsH_pipiAcollinearityzoom, TauSpinEffectsH_pipiAcoplanarity, TauSpinEffectsH_rhorhoAcoplanarityminus, TauSpinEffectsH_rhorhoAcoplanarityplus, TauSpinEffectsH_X, TauSpinEffectsH_Xb, TauSpinEffectsH_Xf, TauSpinEffectsH_Zs, TauSpinEffectsHpm_eX, TauSpinEffectsHpm_muX, TauSpinEffectsHpm_UpsilonA1, TauSpinEffectsHpm_UpsilonRho, TauSpinEffectsHpm_X, TauSpinEffectsW_eX, TauSpinEffectsW_muX, TauSpinEffectsW_UpsilonA1, TauSpinEffectsW_UpsilonRho, TauSpinEffectsW_X, TauSpinEffectsZ_eX, TauSpinEffectsZ_muX, TauSpinEffectsZ_MVis, TauSpinEffectsZ_X, TauSpinEffectsZ_X100to120, TauSpinEffectsZ_X120UP, TauSpinEffectsZ_X50to75, TauSpinEffectsZ_X75to88, TauSpinEffectsZ_X88to100, TauSpinEffectsZ_Xb, TauSpinEffectsZ_Xf, TauSpinEffectsZ_Zs, tmp, cond::to_string(), tripi, tripinpi0, undetermined, W, Z, zsbins, zsmax, and zsmin.

                                                                                           {
    DQMHelper dqm(&i); i.setCurrentFolder("Generator/Tau");
    // Number of analyzed events
    nTaus = dqm.book1dHisto("nTaus", "n analyzed Taus", 1, 0., 1.,"bin","Number of #tau's found");
    nPrimeTaus = dqm.book1dHisto("nPrimeTaus", "n analyzed prime Taus", 1, 0., 1.,"bin","Number of #tau's from Gauge Bosons");

    //Kinematics
    TauPt            = dqm.book1dHisto("TauPt","Tau pT", 100 ,0,100,"P_{T}^{#tau}","Number of #tau's from Gauge Bosons"); 
    TauEta           = dqm.book1dHisto("TauEta","Tau eta", 100 ,-2.5,2.5,"#eta^{#tau}","Number of #tau's from Gauge Bosons");
    TauPhi           = dqm.book1dHisto("TauPhi","Tau phi", 100 ,-3.14,3.14,"#phi^{#tau}","Number of #tau's from Gauge Bosons");
    TauProngs        = dqm.book1dHisto("TauProngs","Tau n prongs", 7 ,0,7,"N_{prongs}","Number of #tau's from Gauge Bosons");
    TauDecayChannels = dqm.book1dHisto("TauDecayChannels","Tau decay channels", 13 ,0,13,"Tau POG Decay Mode","Number of #tau's from Gauge Bosons");
    TauDecayChannels->setBinLabel(1+undetermined,"?");
    TauDecayChannels->setBinLabel(1+electron,"e");
    TauDecayChannels->setBinLabel(1+muon,"mu");
    TauDecayChannels->setBinLabel(1+pi,"#pi^{#pm}");
    TauDecayChannels->setBinLabel(1+rho,"#rho^{#pm}");
    TauDecayChannels->setBinLabel(1+a1,"a_{1}^{#pm}");
    TauDecayChannels->setBinLabel(1+pi1pi0,"#pi^{#pm}#pi^{0}");
    TauDecayChannels->setBinLabel(1+pinpi0,"#pi^{#pm}n#pi^{0}");
    TauDecayChannels->setBinLabel(1+tripi,"3#pi^{#pm}");
    TauDecayChannels->setBinLabel(1+tripinpi0,"3#pi^{#pm}n#pi^{0}");
    TauDecayChannels->setBinLabel(1+K,"K");
    TauDecayChannels->setBinLabel(1+Kstar,"K^{*}");
    TauDecayChannels->setBinLabel(1+stable,"Stable");
    
    TauMothers        = dqm.book1dHisto("TauMothers","Tau mother particles", 10 ,0,10,"Mother of #tau","Number of #tau's");

    TauMothers->setBinLabel(1+other,"?");
    TauMothers->setBinLabel(1+B,"B Decays");
    TauMothers->setBinLabel(1+D,"D Decays");
    TauMothers->setBinLabel(1+gamma,"#gamma");
    TauMothers->setBinLabel(1+Z,"Z");
    TauMothers->setBinLabel(1+W,"W");
    TauMothers->setBinLabel(1+HSM,"H_{SM}/h^{0}");
    TauMothers->setBinLabel(1+H0,"H^{0}");
    TauMothers->setBinLabel(1+A0,"A^{0}");
    TauMothers->setBinLabel(1+Hpm,"H^{#pm}");

    DecayLength = dqm.book1dHisto("DecayLength","#tau Decay Length", 100 ,-20,20,"L_{#tau} (cm)","Number of #tau's from Gauge Bosons");
    LifeTime =  dqm.book1dHisto("LifeTime","#tau LifeTime ", 500 ,0,10000E-15,"#tau_{#tau} (s)","Number of #tau's from Gauge Bosons");
    
    TauSpinEffectsW_X   = dqm.book1dHisto("TauSpinEffectsWX","X for pion", 50 ,0,1,"X","Number of #tau#rightarrow#pi#nu from W^{#pm} Bosons");
    TauSpinEffectsHpm_X = dqm.book1dHisto("TauSpinEffectsHpmX","X for pion", 50 ,0,1,"X","Number of #tau#rightarrow#pi#nu from H^{#pm} Bosons");
    
    TauSpinEffectsW_eX   = dqm.book1dHisto("TauSpinEffectsWeX","X for e", 50 ,0,1,"X","Number of #tau#rightarrowe#nu#nu from W^{#pm} Bosons");
    TauSpinEffectsHpm_eX = dqm.book1dHisto("TauSpinEffectsHpmeX","X for e", 50 ,0,1,"X","Number of #tau#rightarrowe#nu#nu from H^{#pm} Bosons");

    TauSpinEffectsW_muX   = dqm.book1dHisto("TauSpinEffectsWmuX","X for mu", 50 ,0,1,"X","Number of #tau#rightarrow#mu#nu#nu from W^{#pm} Bosons");
    TauSpinEffectsHpm_muX = dqm.book1dHisto("TauSpinEffectsHpmmuX","X for mue", 50 ,0,1,"X","Number of #tau#rightarrow#mu#nu#nu from H^{#pm} Bosons");

    TauSpinEffectsW_UpsilonRho   = dqm.book1dHisto("TauSpinEffectsWUpsilonRho","#Upsilon for #rho", 50 ,-1,1,"#Upsilon","Number of #tau#rightarrow#rho#nu from Gauge Bosons");
    TauSpinEffectsHpm_UpsilonRho = dqm.book1dHisto("TauSpinEffectsHpmUpsilonRho","#Upsilon for #rho", 50 ,-1,1,"#Upsilon","Number of #tau#rightarrow#rho#nu from Gauge Bosons");
    
    TauSpinEffectsW_UpsilonA1   = dqm.book1dHisto("TauSpinEffectsWUpsilonA1","#Upsilon for a1", 50 ,-1,1,"#Upsilon","Number of #tau#rightarrow#pi#pi#pi#nu from Gauge Bosons");
    TauSpinEffectsHpm_UpsilonA1 = dqm.book1dHisto("TauSpinEffectsHpmUpsilonA1","#Upsilon for a1", 50 ,-1,1,"#Upsilon","Number of #tau#rightarrow#pi#pi#pi#nu from Gauge Bosons");

    TauSpinEffectsH_pipiAcoplanarity = dqm.book1dHisto("TauSpinEffectsH_pipiAcoplanarity","H Acoplanarity for #pi^{-}#pi^{+}", 50 ,0,2*TMath::Pi(),"Acoplanarity","Number of H#rightarrow#tau^{-}(#rightarrow#pi^{-}#nu)#tau^{+}(#rightarrow#pi^{+}#nu) Events"); 

    TauSpinEffectsH_pipiAcollinearity = dqm.book1dHisto("TauSpinEffectsH_pipiAcollinearity","H Acollinearity for #pi^{-}#pi^{+}", 50 ,0,TMath::Pi(),"Acollinearity","Number of H#rightarrow#tau^{-}(#rightarrow#pi^{-}#nu)#tau^{+}(#rightarrow#pi^{+}#nu) Events");
    TauSpinEffectsH_pipiAcollinearityzoom = dqm.book1dHisto("TauSpinEffectsH_pipiAcollinearityzoom","H Acollinearity for #pi^{-}#pi^{+}", 50 ,3,TMath::Pi(),"Acollinearity","Number of H#rightarrow#tau^{-}(#rightarrow#pi^{-}#nu)#tau^{+}(#rightarrow#pi^{+}#nu) Events");
    
    TauSpinEffectsZ_MVis   = dqm.book1dHisto("TauSpinEffectsZMVis","Mass of pi+ pi-", 25 ,0,1.1,"M_{#pi^{+}#pi^{-}} (GeV)","Number of Z#rightarrow#tau^{-}(#rightarrow#pi^{-}#nu)#tau^{+}(#rightarrow#pi^{+}#nu) Events");
    TauSpinEffectsH_MVis   = dqm.book1dHisto("TauSpinEffectsHMVis","Mass of pi+ pi-", 25 ,0,1.1,"M_{#pi^{+}#pi^{-}} (GeV)","Number of H#rightarrow#tau^{-}(#rightarrow#pi^{-}#nu)#tau^{+}(#rightarrow#pi^{+}#nu) Events");
    
    TauSpinEffectsZ_Zs   = dqm.book1dHisto("TauSpinEffectsZZs","Z_{s}", zsbins ,zsmin,zsmax,"Z_{s}","Number of Z#rightarrow#tau^{-}(#rightarrow#pi^{-}#nu)#tau^{+}(#rightarrow#pi^{+}#nu Events");
    TauSpinEffectsH_Zs   = dqm.book1dHisto("TauSpinEffectsHZs","Z_{s}", zsbins ,zsmin,zsmax,"Z_{s}","Number of H#rightarrow#tau^{-}(#rightarrow#pi^{-}#nu)#tau^{+}(#rightarrow#pi^{+}#nu Events");
    
    TauSpinEffectsZ_X= dqm.book1dHisto("TauSpinEffectsZX","X for pion of #tau^{-}", 25 ,0,1.0,"X","Number of #tau#rightarrow#pi#nu from Z Bosons");
    TauSpinEffectsZ_X50to75= dqm.book1dHisto("TauSpinEffectsZX50to75","X for pion of #tau^{-} (50GeV-75GeV)", 10 ,0,1.0,"X","Number of #tau#rightarrow#pi#nu from Z(50GeV<M<75GeV) Bosons");
    TauSpinEffectsZ_X75to88= dqm.book1dHisto("TauSpinEffectsZX75to88","X for pion of #tau^{-} (75GeV-88GeV)", 10 ,0,1.0,"X","Number of #tau#rightarrow#pi#nu from Z(75GeV<M<88GeV) Bosons");
    TauSpinEffectsZ_X88to100= dqm.book1dHisto("TauSpinEffectsZX88to100","X for pion of #tau^{-} (88GeV-100GeV)", 10 ,0,1.0,"X","Number of #tau#rightarrow#pi#nu from Z(88GeV<M<100GeV) Bosons");
    TauSpinEffectsZ_X100to120= dqm.book1dHisto("TauSpinEffectsZX100to120","X for pion of #tau^{-} (100GeV-120GeV)", 10 ,0,1.0,"X","Number of #tau#rightarrow#pi#nu from Z(100GeV<M<120GeV) Bosons");
    TauSpinEffectsZ_X120UP= dqm.book1dHisto("TauSpinEffectsZX120UP","X for pion of #tau^{-} (>120GeV)", 10 ,0,1.0,"X","Number of #tau#rightarrow#pi#nu from Z(120GeV<MGeV) Bosons");


    TauSpinEffectsH_X= dqm.book1dHisto("TauSpinEffectsH_X","X for pion of #tau^{-}", 25 ,0,1.0,"X","Number of #tau#rightarrow#pi#nu from H Bosons");
    
    TauSpinEffectsZ_Xf   = dqm.book1dHisto("TauSpinEffectsZXf","X for pion of forward emitted #tau^{-}", 25 ,0,1.0,"X_{f}","Number of #tau#rightarrow#pi#nu from Z Bosons");
    TauSpinEffectsH_Xf   = dqm.book1dHisto("TauSpinEffectsHXf","X for pion of forward emitted #tau^{-}", 25 ,0,1.0,"X_{f}","Number of #tau#rightarrow#pi#nu from H Bosons");
    
    TauSpinEffectsZ_Xb   = dqm.book1dHisto("TauSpinEffectsZXb","X for pion of backward emitted #tau^{-}", 25 ,0,1.0,"X_{b}","Number of #tau#rightarrow#pi#nu from Z Bosons");
    TauSpinEffectsH_Xb   = dqm.book1dHisto("TauSpinEffectsHXb","X for pion of backward emitted #tau^{-}", 25 ,0,1.0,"X_{b}","Number of #tau#rightarrow#pi#nu from H Bosons");

    TauSpinEffectsZ_eX   = dqm.book1dHisto("TauSpinEffectsZeX","X for e", 50 ,0,1,"X","Number of #tau#rightarrowe#nu#nu from Gauge Bosons");
    TauSpinEffectsH_eX = dqm.book1dHisto("TauSpinEffectsHeX","X for e", 50 ,0,1,"X","Number of #tau#rightarrowe#nu#nu from Gauge Bosons");

    TauSpinEffectsZ_muX   = dqm.book1dHisto("TauSpinEffectsZmuX","X for mu", 50 ,0,1,"X","Number of #tau#rightarrow#mu#nu#nu from Gauge Bosons");
    TauSpinEffectsH_muX = dqm.book1dHisto("TauSpinEffectsHmuX","X for mu", 50 ,0,1,"X","Number of #tau#rightarrow#mu#nu#nu from Gauge Bosons");
    
    TauSpinEffectsH_rhorhoAcoplanarityminus = dqm.book1dHisto("TauSpinEffectsH_rhorhoAcoplanarityminus","#phi^{*-} (acoplanarity) for Higgs #rightarrow #rho-#rho (y_{1}*y_{2}<0)", 32 ,0,2*TMath::Pi(),"#phi^{*-} (Acoplanarity)","Number of H#rightarrow#tau^{-}(#rightarrow#rho^{-}#nu)#tau^{+}(#rightarrow#rho^{+}#nu) Events");
    TauSpinEffectsH_rhorhoAcoplanarityplus = dqm.book1dHisto("TauSpinEffectsH_rhorhoAcoplanarityplus","#phi^{*+} (acoplanarity) for Higgs #rightarrow #rho-#rho (y_{1}*y_{2}>0)", 32 ,0,2*TMath::Pi(),"#phi^{*+} (Acoplanarity)","Number of H#rightarrow#tau^{-}(#rightarrow#rho^{-}#nu)#tau^{+}(#rightarrow#rho^{+}#nu) Events");
    
    TauFSRPhotonsN=dqm.book1dHisto("TauFSRPhotonsN","FSR Photons radiating from/with tau (Gauge Boson)", 5 ,-0.5,4.5,"N^{FSR Photons radiating from/with #tau}","Number of #tau's from Gauge Bosons");
    TauFSRPhotonsPt=dqm.book1dHisto("TauFSRPhotonsPt","Pt of FSR Photons radiating from/with tau (Gauge Boson)", 100 ,0,100,"P_{t}^{FSR Photons radiating from/with #tau [per #tau]} (GeV)","Number of #tau's from Gauge Bosons");
    TauFSRPhotonsPtSum=dqm.book1dHisto("TauFSRPhotonsPtSum","Pt of FSR Photons radiating from/with tau (Gauge Boson)", 100 ,0,100,"P_{t}^{FSR Photons radiating from/with #tau [per #tau]} (GeV)","Number of #tau's from Gauge Bosons");

    TauBremPhotonsN=dqm.book1dHisto("TauBremPhotonsN","Brem. Photons radiating in tau decay", 5 ,-0.5,4.5,"N FSR Photons radiating from/with tau","Number of #tau's from Gauge Bosons");
    TauBremPhotonsPt=dqm.book1dHisto("TauBremPhotonsPt","Sum Brem Pt ", 100 ,0,100,"P_{t}^{Brem. Photons radiating in tau decay} (GeV)","Number of #tau's from Gauge Bosons");    
    TauBremPhotonsPtSum =dqm.book1dHisto("TauBremPhotonsPtSum","Sum of Brem Pt ", 100 ,0,100,"Sum P_{t}^{Brem. Photons radiating in tau decay} (GeV)","Number of #tau's from Gauge Bosons");
    
    MODEID =dqm.book1dHisto("JAKID","JAK ID",NMODEID+1,-0.5,NMODEID+0.5);
    for(unsigned int j=0; j<NMODEID+1;j++){
      MODEInvMass.push_back(std::vector<MonitorElement *>());
      std::string tmp="JAKID";
      tmp+=std::to_string(j);
      MODEInvMass.at(j).push_back(dqm.book1dHisto("M"+tmp,"M_{"+TauDecay::DecayMode(j)+"} (GeV)", 80 ,0,2.0,"M_{"+TauDecay::DecayMode(j)+"} (GeV)","Number of #tau's from Gauge Bosons"));
      MODEID->setBinLabel(1+j,TauDecay::DecayMode(j));
      if(j==TauDecay::MODE_3PI || j==TauDecay::MODE_PI2PI0 ||
     j==TauDecay::MODE_KPIK ||
     j==TauDecay::MODE_KPIPI ){
    MODEInvMass.at(j).push_back(dqm.book1dHisto("M13"+tmp,"M_{13,"+TauDecay::DecayMode(j)+"} (GeV)", 80 ,0,2.0,"M_{13,"+TauDecay::DecayMode(j)+"} (GeV)","Number of #tau's from Gauge Bosons"));
    MODEInvMass.at(j).push_back(dqm.book1dHisto("M23"+tmp,"M_{23,"+TauDecay::DecayMode(j)+"} (GeV)", 80 ,0,2.0,"M_{23,"+TauDecay::DecayMode(j)+"} (GeV)","Number of #tau's from Gauge Bosons"));
    MODEInvMass.at(j).push_back(dqm.book1dHisto("M12"+tmp,"M_{12,"+TauDecay::DecayMode(j)+"} (GeV)", 80 ,0,2.0,"M_{12,"+TauDecay::DecayMode(j)+"} (GeV)","Number of #tau's from Gauge Bosons"));
      }
    }
  return;
}

void TauValidation::countParticles ( const reco::GenParticle *  p, int &  allCount, int &  eCount, int &  muCount, int &  pi0Count, int &  piCount, int &  rhoCount, int &  a1Count, int &  KCount, int &  KstarCount  )

private

Definition at line 388 of file TauValidation.cc.

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

Referenced by tauDecayChannel().

                                                                                                     {
  for(unsigned int i=0; i<p->numberOfDaughters(); i++){
    const reco::GenParticle *dau=static_cast<const reco::GenParticle*>(p->daughter(i));
    int pid = dau->pdgId();
    allCount++;
    if(abs(pid) == 11)    eCount++;
    if(abs(pid) == 13)    muCount++;
    if(abs(pid) == 111)   pi0Count++;
    if(abs(pid) == 211)   piCount++;
    if(abs(pid) == 213)   rhoCount++;
    if(abs(pid) == 20213) a1Count++;
    if(abs(pid) == 321)   KCount++;
    if(abs(pid) == 323)   KstarCount++;
    countParticles(dau,allCount,eCount,muCount,pi0Count,piCount,rhoCount,a1Count,KCount,KstarCount);
  }
}

void TauValidation::dqmBeginRun ( const edm::Run &  r, const edm::EventSetup &  c  )

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 31 of file TauValidation.cc.

References fPDGTable, and edm::EventSetup::getData().

                                                                       {
  c.getData( fPDGTable );
}

void TauValidation::findFSRandBrem ( const reco::GenParticle *  p, bool  doBrem, std::vector< const reco::GenParticle * > &  ListofFSR, std::vector< const reco::GenParticle * > &  ListofBrem  )

private

Definition at line 295 of file TauValidation.cc.

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

Referenced by photons().

                                                                {
  // note this code split the FSR and Brem based one if the tau decays into a tau+photon or not with the Fortran Tauola Interface, this is not 100% correct because photos puts the tau with the regular tau decay products. 
  if(abs(p->pdgId())==15){
    if(isLastTauinChain(p)){ doBrem=true;}
    else{ doBrem=false;}
  }
  int photo_ID=22;
  for(unsigned int i = 0; i <p->numberOfDaughters(); i++){
    const reco::GenParticle *dau=static_cast<const reco::GenParticle*>(p->daughter(i));
    if(abs((dau)->pdgId()) == abs(photo_ID) && !doBrem){ListofFSR.push_back(dau);}
    if(abs((dau)->pdgId()) == abs(photo_ID) && doBrem){ListofBrem.push_back(dau);}
    if(abs((dau)->pdgId()) != 111 && abs((dau)->pdgId()) != 221){ // remove pi0 and eta decays
      findFSRandBrem(dau,doBrem,ListofFSR,ListofBrem);
    }
  }
}

int TauValidation::findMother ( const reco::GenParticle *  tau )

private

Definition at line 274 of file TauValidation.cc.

References GetMother(), and reco::LeafCandidate::pdgId().

Referenced by spinEffectsZH(), and tauMother().

                                                       {
  return TauValidation::GetMother(tau)->pdgId();
}

void TauValidation::FindPhotosFSR ( const reco::GenParticle *  p, std::vector< const reco::GenParticle * > &  ListofFSR, double &  BosonScale  )

private

Definition at line 313 of file TauValidation.cc.

References funct::abs(), reco::CompositeRefCandidateT< D >::daughter(), GetMother(), i, visualization-live-secondInstance_cfg::m, reco::CompositeRefCandidateT< D >::numberOfDaughters(), and reco::LeafCandidate::pdgId().

Referenced by photons().

                                                                                                                         {
  BosonScale=0.0;
  const reco::GenParticle* m=GetMother(p);
  int mother_pid=m->pdgId();
  if(m->pdgId()!=p->pdgId()){
    for(unsigned int i=0; i <m->numberOfDaughters(); i++){
      const reco::GenParticle *dau=static_cast<const reco::GenParticle*>(m->daughter(i));
      if(abs(dau->pdgId()) == 22) {
    ListofFSR.push_back(dau);
      }
    }
  }
  if(abs(mother_pid) == 24) BosonScale=1.0; // W
  if(abs(mother_pid) == 23) BosonScale=2.0; // Z;
  if(abs(mother_pid) == 22) BosonScale=2.0; // gamma;
  if(abs(mother_pid) == 25) BosonScale=2.0; // HSM;
  if(abs(mother_pid) == 35) BosonScale=2.0; // H0;
  if(abs(mother_pid) == 36) BosonScale=2.0; // A0;
  if(abs(mother_pid) == 37) BosonScale=1.0; //Hpm;
}

void TauValidation::findTauList ( const reco::GenParticle *  tau, std::vector< const reco::GenParticle * > &  TauList  )

private

Definition at line 285 of file TauValidation.cc.

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

Referenced by photons().

                                                                                                    {
  TauList.insert(TauList.begin(),tau);
  for(unsigned int i=0;i<tau->numberOfMothers();i++) {
    const reco::GenParticle *mother=static_cast<const reco::GenParticle*>(tau->mother(i));
    if(mother->pdgId() == tau->pdgId()){
      findTauList(mother,TauList);
    }
  }
}

const reco::GenParticle * TauValidation::GetMother ( const reco::GenParticle *  tau )

private

Definition at line 255 of file TauValidation.cc.

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

Referenced by findMother(), FindPhotosFSR(), and motherP4().

                                                                         {
  for (unsigned int i=0;i<tau->numberOfMothers();i++) {
    const reco::GenParticle *mother=static_cast<const reco::GenParticle*>(tau->mother(i));
    if(mother->pdgId() == tau->pdgId()) return GetMother(mother);
    return mother;
  }
  return tau;
}

const std::vector< const reco::GenParticle * > TauValidation::GetMothers ( const reco::GenParticle *  boson )

private

Definition at line 264 of file TauValidation.cc.

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

Referenced by spinEffectsZH().

                                                                                             {
  std::vector<const reco::GenParticle*> mothers;
  for (unsigned int i=0;i<boson->numberOfMothers();i++) {
    const reco::GenParticle *mother=static_cast<const reco::GenParticle*>(boson->mother(i));
    if(mother->pdgId() == boson->pdgId()) return GetMothers(mother);
    mothers.push_back(mother);
  }
  return mothers;
}

bool TauValidation::isLastTauinChain ( const reco::GenParticle *  tau )

private

Definition at line 278 of file TauValidation.cc.

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

Referenced by analyze(), and findFSRandBrem().

                                                              {
  for(unsigned int i = 0; i <tau->numberOfDaughters(); i++){
    if(tau->daughter(i)->pdgId() == tau->pdgId()) return false;
  }
  return true;
}

double TauValidation::leadingPionMomentum ( const reco::GenParticle *  tau, double  weight  )

private

Definition at line 678 of file TauValidation.cc.

References leadingPionP4().

                                                                                  {
  return leadingPionP4(tau).P();
}

TLorentzVector TauValidation::leadingPionP4 ( const reco::GenParticle *  tau )

private

Definition at line 682 of file TauValidation.cc.

References funct::abs(), reco::CompositeRefCandidateT< D >::daughter(), reco::LeafCandidate::energy(), i, reco::CompositeRefCandidateT< D >::numberOfDaughters(), reco::LeafCandidate::p(), p4, reco::LeafCandidate::pdgId(), sysUtil::pid, reco::LeafCandidate::px(), reco::LeafCandidate::py(), and reco::LeafCandidate::pz().

Referenced by leadingPionMomentum(), spinEffectsWHpm(), and spinEffectsZH().

                                                                     {
  TLorentzVector p4(0,0,0,0);
  for(unsigned int i = 0; i <tau->numberOfDaughters(); i++){
    const reco::GenParticle *dau=static_cast<const reco::GenParticle*>(tau->daughter(i));
    int pid = dau->pdgId();
    if(abs(pid) == 15) return leadingPionP4(dau);
    if(!(abs(pid)==211 || abs(pid)==13 || abs(pid)==11)) continue;
    if(dau->p() > p4.P()) p4 = TLorentzVector(dau->px(),dau->py(),dau->pz(),dau->energy());
  }
  return p4;
}

TLorentzVector TauValidation::motherP4 ( const reco::GenParticle *  tau )

private

Definition at line 694 of file TauValidation.cc.

References reco::LeafCandidate::energy(), GetMother(), visualization-live-secondInstance_cfg::m, reco::LeafCandidate::px(), reco::LeafCandidate::py(), and reco::LeafCandidate::pz().

Referenced by spinEffectsWHpm().

                                                                {
  const reco::GenParticle* m=GetMother(tau);
  return TLorentzVector(m->px(),m->py(),m->pz(),m->energy());
}

void TauValidation::photons ( const reco::GenParticle *  tau, double  weight  )

private

Definition at line 712 of file TauValidation.cc.

References MonitorElement::Fill(), findFSRandBrem(), FindPhotosFSR(), findTauList(), i, TauBremPhotonsN, TauBremPhotonsPt, TauBremPhotonsPtSum, TauFSRPhotonsN, TauFSRPhotonsPt, TauFSRPhotonsPtSum, and histoStyle::weight.

Referenced by analyze().

                                                                    {
  // Find First tau in chain
  std::vector<const reco::GenParticle*> TauList;
  findTauList(tau,TauList);

  // Get List of Gauge Boson to tau(s) FSR and Brem
  bool passedW=false;
  std::vector<const reco::GenParticle*> ListofFSR;  ListofFSR.clear();
  std::vector<const reco::GenParticle*> ListofBrem; ListofBrem.clear();
  std::vector<const reco::GenParticle*> FSR_photos; FSR_photos.clear();
  double BosonScale(1);
  if(TauList.size()>0){
    TauValidation::findFSRandBrem(TauList.at(0),passedW,ListofFSR,ListofBrem);
    TauValidation::FindPhotosFSR(TauList.at(0),FSR_photos,BosonScale);

    // Add the Tau Brem. information
    TauBremPhotonsN->Fill(ListofBrem.size(),weight);
    double photonPtSum=0;
    for(unsigned int i=0;i<ListofBrem.size();i++){
      photonPtSum+=ListofBrem.at(i)->pt();
      TauBremPhotonsPt->Fill(ListofBrem.at(i)->pt(),weight);
    }
    TauBremPhotonsPtSum->Fill(photonPtSum,weight);
        
    // Now add the Gauge Boson FSR information
    if(BosonScale!=0){
      TauFSRPhotonsN->Fill(ListofFSR.size(),weight);
      photonPtSum=0;
      for(unsigned int i=0;i<ListofFSR.size();i++){
    photonPtSum+=ListofFSR.at(i)->pt();
    TauFSRPhotonsPt->Fill(ListofFSR.at(i)->pt(),weight);
      }
      double FSR_photosSum(0);
      for(unsigned int i=0;i<FSR_photos.size();i++){
    FSR_photosSum+=FSR_photos.at(i)->pt();
    TauFSRPhotonsPt->Fill(FSR_photos.at(i)->pt()/BosonScale,weight*BosonScale);
      }
      TauFSRPhotonsPtSum->Fill(photonPtSum+FSR_photosSum/BosonScale,weight);
    }
  }
}

void TauValidation::spinEffectsWHpm ( const reco::GenParticle *  tau, int  mother, int  decay, std::vector< const reco::GenParticle * > &  part, double  weight  )

private

Definition at line 408 of file TauValidation.cc.

References a1, funct::abs(), relval_parameters_module::energy, MonitorElement::Fill(), gamma, i, leadingPionP4(), TauDecay::MODE_3PI, TauDecay::MODE_ELECTRON, TauDecay::MODE_MUON, TauDecay::MODE_PI2PI0, TauDecay::MODE_PION, TauDecay::MODE_PIPI0, motherP4(), pi, PdtPdgMini::pi0, PdtPdgMini::pi_minus, PdtPdgMini::pi_plus, rho, TauSpinEffectsHpm_eX, TauSpinEffectsHpm_muX, TauSpinEffectsHpm_UpsilonA1, TauSpinEffectsHpm_UpsilonRho, TauSpinEffectsHpm_X, TauSpinEffectsW_eX, TauSpinEffectsW_muX, TauSpinEffectsW_UpsilonA1, TauSpinEffectsW_UpsilonRho, TauSpinEffectsW_X, and histoStyle::weight.

Referenced by analyze().

                                                                                                                                          {
  if(decay == TauDecay::MODE_PION || decay == TauDecay::MODE_MUON || decay == TauDecay::MODE_ELECTRON){  // polarization only for 1-prong hadronic taus with no neutral pions
    TLorentzVector momP4 = motherP4(tau);
    TLorentzVector pionP4 = leadingPionP4(tau);
    pionP4.Boost(-1*momP4.BoostVector());
    double energy = pionP4.E()/(momP4.M()/2);
    if(decay == TauDecay::MODE_PION){
      if(abs(mother) == 24) TauSpinEffectsW_X->Fill(energy,weight); 
      if(abs(mother) == 37) TauSpinEffectsHpm_X->Fill(energy,weight);
    }
    if(decay == TauDecay::MODE_MUON){
      if(abs(mother) == 24) TauSpinEffectsW_muX->Fill(energy,weight);
      if(abs(mother) == 37) TauSpinEffectsHpm_muX->Fill(energy,weight);
    }
    if(decay == TauDecay::MODE_ELECTRON){
      if(abs(mother) == 24) TauSpinEffectsW_eX->Fill(energy,weight);
      if(abs(mother) == 37) TauSpinEffectsHpm_eX->Fill(energy,weight);
    }
  }
  else if(decay==TauDecay::MODE_PIPI0){
    TLorentzVector rho(0,0,0,0),pi(0,0,0,0);
    for(unsigned int i=0;i<part.size();i++){
      TLorentzVector LV(part.at(i)->px(),part.at(i)->py(),part.at(i)->pz(),part.at(i)->energy());
      if(abs(part.at(i)->pdgId())==PdtPdgMini::pi_plus){pi+=LV; rho+=LV;}
      if(abs(part.at(i)->pdgId())==PdtPdgMini::pi0){rho+=LV;}
    }
    if(abs(mother) == 24) TauSpinEffectsW_UpsilonRho->Fill(2*pi.P()/rho.P()-1,weight);
    if(abs(mother) == 37) TauSpinEffectsHpm_UpsilonRho->Fill(2*pi.P()/rho.P()-1,weight);
  }
  else if(decay==TauDecay::MODE_3PI || decay==TauDecay::MODE_PI2PI0){ // only for pi2pi0 for now
    TLorentzVector a1(0,0,0,0),pi_p(0,0,0,0),pi_m(0,0,0,0); 
    int nplus(0),nminus(0);
    for(unsigned int i=0;i<part.size();i++){
      TLorentzVector LV(part.at(i)->px(),part.at(i)->py(),part.at(i)->pz(),part.at(i)->energy());
      if(part.at(i)->pdgId()==PdtPdgMini::pi_plus){ pi_p+=LV; a1+=LV; nplus++;}
      if(part.at(i)->pdgId()==PdtPdgMini::pi_minus){pi_m+=LV; a1+=LV; nminus++;}
    }
    double gamma=0;
    if(nplus+nminus==3 && nplus==1)  gamma=2*pi_p.P()/a1.P()-1;
    if(nplus+nminus==3 && nminus==1) gamma=2*pi_m.P()/a1.P()-1;
    else{
      pi_p+=pi_m; gamma=2*pi_p.P()/a1.P()-1;
    }
    if(abs(mother) == 24) TauSpinEffectsW_UpsilonA1->Fill(gamma,weight);
    if(abs(mother) == 37) TauSpinEffectsHpm_UpsilonA1->Fill(gamma,weight);
  }
}

void TauValidation::spinEffectsZH ( const reco::GenParticle *  boson, double  weight  )

private

Definition at line 456 of file TauValidation.cc.

References funct::abs(), TauDecay_GenParticle::AnalyzeTau(), PdtPdgMini::anti_d, PdtPdgMini::anti_u, RecoTauCleanerPlugins::charge, PdtPdgMini::d, reco::CompositeRefCandidateT< D >::daughter(), reco::LeafCandidate::energy(), MonitorElement::Fill(), findMother(), fPDGTable, TauDecay_GenParticle::Get_TauDecayProducts(), GetMothers(), PdtPdgMini::Higgs0, i, leadingPionP4(), visualization-live-secondInstance_cfg::m, TauDecay::MODE_ELECTRON, TauDecay::MODE_MUON, TauDecay::MODE_PION, TauDecay::MODE_PIPI0, gen::n, reco::CompositeRefCandidateT< D >::numberOfDaughters(), reco::LeafCandidate::pdgId(), Pi, sysUtil::pid, reco::LeafCandidate::px(), reco::LeafCandidate::py(), reco::LeafCandidate::pz(), lumiQueryAPI::q, TauSpinEffectsH_eX, TauSpinEffectsH_muX, TauSpinEffectsH_MVis, TauSpinEffectsH_pipiAcollinearity, TauSpinEffectsH_pipiAcollinearityzoom, TauSpinEffectsH_pipiAcoplanarity, TauSpinEffectsH_rhorhoAcoplanarityminus, TauSpinEffectsH_rhorhoAcoplanarityplus, TauSpinEffectsH_X, TauSpinEffectsH_Xb, TauSpinEffectsH_Xf, TauSpinEffectsH_Zs, TauSpinEffectsZ_eX, TauSpinEffectsZ_muX, TauSpinEffectsZ_MVis, TauSpinEffectsZ_X, TauSpinEffectsZ_X100to120, TauSpinEffectsZ_X120UP, TauSpinEffectsZ_X50to75, TauSpinEffectsZ_X75to88, TauSpinEffectsZ_X88to100, TauSpinEffectsZ_Xb, TauSpinEffectsZ_Xf, TauSpinEffectsZ_Zs, PdtPdgMini::u, histoStyle::weight, Z, PdtPdgMini::Z0, zsbins, zsmax, zsmin, and Zstoa().

Referenced by analyze().

                                                                            {
  int ntau(0);
  for(unsigned int i = 0; i <boson->numberOfDaughters(); i++){
    const reco::GenParticle *dau=static_cast<const reco::GenParticle*>(boson->daughter(i));
    if(ntau==1 && dau->pdgId() == 15)return;
    if(boson->pdgId()!= 15 && abs(dau->pdgId()) == 15)ntau++;
  }
  if(ntau!=2) return; 
  if(abs(boson->pdgId())==PdtPdgMini::Z0 || abs(boson->pdgId())==PdtPdgMini::Higgs0){
    TLorentzVector tautau(0,0,0,0);
    TLorentzVector pipi(0,0,0,0);
    TLorentzVector taum(0,0,0,0);
    TLorentzVector taup(0,0,0,0);
    TLorentzVector rho_plus,rho_minus,pi_rhominus,pi0_rhominus,pi_rhoplus,pi0_rhoplus,pi_plus,pi_minus;
    bool hasrho_minus(false),hasrho_plus(false),haspi_minus(false),haspi_plus(false);
    int nSinglePionDecays(0),nSingleMuonDecays(0),nSingleElectronDecays(0);
    double x1(0),x2(0); 
    TLorentzVector Zboson(boson->px(),boson->py(),boson->pz(),boson->energy());
    for(unsigned int i = 0; i <boson->numberOfDaughters(); i++){
      const reco::GenParticle *dau=static_cast<const reco::GenParticle*>(boson->daughter(i));
      int pid = dau->pdgId();
      if(abs(findMother(dau)) != 15 && abs(pid) == 15){
    TauDecay_GenParticle TD;
    unsigned int jak_id, TauBitMask;
    if(TD.AnalyzeTau(dau,jak_id,TauBitMask,false,false)){
      std::vector<const reco::GenParticle*> part=TD.Get_TauDecayProducts();
      if(jak_id==TauDecay::MODE_PION || jak_id==TauDecay::MODE_MUON || jak_id==TauDecay::MODE_ELECTRON){
        if(jak_id==TauDecay::MODE_PION)     nSinglePionDecays++;
        if(jak_id==TauDecay::MODE_MUON)     nSingleMuonDecays++;
        if(jak_id==TauDecay::MODE_ELECTRON) nSingleElectronDecays++;
        TLorentzVector LVtau(dau->px(),dau->py(),dau->pz(),dau->energy());
        tautau += LVtau;
        TLorentzVector LVpi=leadingPionP4(dau);
        pipi+=LVpi;
        const HepPDT::ParticleData*  pd = fPDGTable->particle(dau->pdgId ());
        int charge = (int) pd->charge();
        LVtau.Boost(-1*Zboson.BoostVector());
        LVpi.Boost(-1*Zboson.BoostVector());
     


        if(jak_id==TauDecay::MODE_MUON){
          if(abs(boson->pdgId())==PdtPdgMini::Z0)     TauSpinEffectsZ_muX->Fill(LVpi.P()/LVtau.E(),weight);
          if(abs(boson->pdgId())==PdtPdgMini::Higgs0) TauSpinEffectsH_muX->Fill(LVpi.P()/LVtau.E(),weight);

        }
        if(jak_id==TauDecay::MODE_ELECTRON){
          if(abs(boson->pdgId())==PdtPdgMini::Z0)     TauSpinEffectsZ_eX->Fill(LVpi.P()/LVtau.E(),weight);
          if(abs(boson->pdgId())==PdtPdgMini::Higgs0) TauSpinEffectsH_eX->Fill(LVpi.P()/LVtau.E(),weight);
        }


        if(jak_id==TauDecay::MODE_PION){
          if(abs(boson->pdgId())==PdtPdgMini::Z0){
        TauSpinEffectsZ_X->Fill(LVpi.P()/LVtau.E(),weight);
        if(50.0<Zboson.M() && Zboson.M()<75.0) TauSpinEffectsZ_X50to75->Fill(LVpi.P()/LVtau.E(),weight);
        if(75.0<Zboson.M() && Zboson.M()<88.0) TauSpinEffectsZ_X75to88->Fill(LVpi.P()/LVtau.E(),weight);
        if(88.0<Zboson.M() && Zboson.M()<100.0) TauSpinEffectsZ_X88to100->Fill(LVpi.P()/LVtau.E(),weight);
        if(100.0<Zboson.M() && Zboson.M()<120.0) TauSpinEffectsZ_X100to120->Fill(LVpi.P()/LVtau.E(),weight);
        if(120.0<Zboson.M()) TauSpinEffectsZ_X120UP->Fill(LVpi.P()/LVtau.E(),weight);
          }
          if(abs(boson->pdgId())==PdtPdgMini::Higgs0) TauSpinEffectsH_X->Fill(LVpi.P()/LVtau.E(),weight);
        }
        if(charge<0){x1=LVpi.P()/LVtau.E(); taum=LVtau;}
        else{ x2=LVpi.P()/LVtau.E();}
      }
      TLorentzVector LVtau(dau->px(),dau->py(),dau->pz(),dau->energy());
      if(pid == 15)taum=LVtau;
      if(pid ==-15)taup=LVtau;
      if(jak_id==TauDecay::MODE_PIPI0){
        for(unsigned int i=0; i<part.size();i++){
          int pid_d = part.at(i)->pdgId();
          if(abs(pid_d)==211 || abs(pid_d)==111){
        TLorentzVector LV(part.at(i)->px(),part.at(i)->py(),part.at(i)->pz(),part.at(i)->energy());
        if(pid==15){
          hasrho_minus=true;
          if(pid_d==-211 ){ pi_rhominus=LV;}
          if(abs(pid_d)==111 ){ pi0_rhominus=LV;}
        }
        if(pid==-15){
          hasrho_plus=true;
          if(pid_d==211 ){pi_rhoplus=LV;}
          if(abs(pid_d)==111 ){pi0_rhoplus=LV;} 
        }
          }
        }
      }
      if(jak_id==TauDecay::MODE_PION){
        for(unsigned int i=0; i<part.size();i++){
          int pid_d = part.at(i)->pdgId();
          if(abs(pid_d)==211 ){
        TLorentzVector LV(part.at(i)->px(),part.at(i)->py(),part.at(i)->pz(),part.at(i)->energy());
        if(pid==15){
          haspi_minus=true;
          if(pid_d==-211 ){ pi_minus=LV;}
        }
        if(pid==-15){
          haspi_plus=true;
          if(pid_d==211 ){pi_plus=LV;}
        }
          }
        }
      }
    }
      }
    }
    if(hasrho_minus && hasrho_plus){
      //compute rhorho
      rho_minus=pi_rhominus;
      rho_minus+=pi0_rhominus;
      rho_plus=pi_rhoplus;
      rho_plus+=pi0_rhoplus;
      TLorentzVector rhorho=rho_minus;rhorho+=rho_plus;

      // boost to rhorho cm
      TLorentzVector pi_rhoplusb=pi_rhoplus;     pi_rhoplusb.Boost(-1*rhorho.BoostVector());
      TLorentzVector pi0_rhoplusb=pi0_rhoplus;   pi0_rhoplusb.Boost(-1*rhorho.BoostVector());
      TLorentzVector pi_rhominusb=pi_rhominus;   pi_rhominusb.Boost(-1*rhorho.BoostVector());
      TLorentzVector pi0_rhominusb=pi0_rhominus; pi0_rhominusb.Boost(-1*rhorho.BoostVector());
  
      // compute n+/-
      TVector3 n_plus=pi_rhoplusb.Vect().Cross(pi0_rhoplusb.Vect());
      TVector3 n_minus=pi_rhominusb.Vect().Cross(pi0_rhominusb.Vect());
      
      // compute the acoplanarity
      double Acoplanarity=acos(n_plus.Dot(n_minus)/(n_plus.Mag()*n_minus.Mag()));
      if(pi_rhominusb.Vect().Dot(n_plus)>0){Acoplanarity*=-1;Acoplanarity+=2*TMath::Pi();}
  
      // now boost to tau frame
      pi_rhoplus.Boost(-1*taup.BoostVector());
      pi0_rhoplus.Boost(-1*taup.BoostVector());
      pi_rhominus.Boost(-1*taum.BoostVector());
      pi0_rhominus.Boost(-1*taum.BoostVector());
      
      // compute y1 and y2
      double y1=(pi_rhoplus.E()-pi0_rhoplus.E())/(pi_rhoplus.E()+pi0_rhoplus.E());
      double y2=(pi_rhominus.E()-pi0_rhominus.E())/(pi_rhominus.E()+pi0_rhominus.E());
      
      // fill histograms
      if(abs(boson->pdgId())==PdtPdgMini::Higgs0 && y1*y2<0) TauSpinEffectsH_rhorhoAcoplanarityminus->Fill(Acoplanarity,weight);
      if(abs(boson->pdgId())==PdtPdgMini::Higgs0 && y1*y2>0) TauSpinEffectsH_rhorhoAcoplanarityplus->Fill(Acoplanarity,weight);
    }
    if(haspi_minus && haspi_plus){
      TLorentzVector tauporig=taup;
      TLorentzVector taumorig=taum;
      
      // now boost to Higgs frame
      pi_plus.Boost(-1*Zboson.BoostVector());
      pi_minus.Boost(-1*Zboson.BoostVector());
      
      taup.Boost(-1*Zboson.BoostVector());
      taum.Boost(-1*Zboson.BoostVector());
      
      if(abs(boson->pdgId())==PdtPdgMini::Higgs0){
    TauSpinEffectsH_pipiAcollinearity->Fill(acos(pi_plus.Vect().Dot(pi_minus.Vect())/(pi_plus.P()*pi_minus.P())));
    TauSpinEffectsH_pipiAcollinearityzoom->Fill(acos(pi_plus.Vect().Dot(pi_minus.Vect())/(pi_plus.P()*pi_minus.P())));
      }
      
      double proj_m=taum.Vect().Dot(pi_minus.Vect())/(taum.P()*taum.P());
      double proj_p=taup.Vect().Dot(pi_plus.Vect())/(taup.P()*taup.P());
      TVector3 Tau_m=taum.Vect();
      TVector3 Tau_p=taup.Vect();
      Tau_m*=proj_m;
      Tau_p*=proj_p;
      TVector3 Pit_m=pi_minus.Vect()-Tau_m;
      TVector3 Pit_p=pi_plus.Vect()-Tau_p;
      
      double Acoplanarity=acos(Pit_m.Dot(Pit_p)/(Pit_p.Mag()*Pit_m.Mag()));
      TVector3 n=Pit_p.Cross(Pit_m);
      if(n.Dot(Tau_m)/Tau_m.Mag()>0){Acoplanarity*=-1; Acoplanarity+=2*TMath::Pi();}
      // fill histograms
      if(abs(boson->pdgId())==PdtPdgMini::Higgs0) TauSpinEffectsH_pipiAcoplanarity->Fill(Acoplanarity,weight);
      taup=tauporig;
      taum=taumorig;
    }
    if(nSinglePionDecays == 2 && tautau.M()!= 0) {
      for(int i=0;i<zsbins;i++){
    double zslow=((double)i)*(zsmax-zsmin)/((double)zsbins)+zsmin; 
    double zsup=((double)i+1)*(zsmax-zsmin)/((double)zsbins)+zsmin;
    double aup=Zstoa(zsup), alow=Zstoa(zslow);
    if(x2-x1>alow && x2-x1<aup){
      double zs=(zsup+zslow)/2;
      if(abs(boson->pdgId())==PdtPdgMini::Z0)     TauSpinEffectsZ_Zs->Fill(zs,weight);
      if(abs(boson->pdgId())==PdtPdgMini::Higgs0) TauSpinEffectsH_Zs->Fill(zs,weight);
      break;
    }
      }
      if(abs(boson->pdgId())==PdtPdgMini::Z0)     TauSpinEffectsZ_MVis->Fill(pipi.M()/tautau.M(),weight);
      if(abs(boson->pdgId())==PdtPdgMini::Higgs0) TauSpinEffectsH_MVis->Fill(pipi.M()/tautau.M(),weight);
      
      if(x1!=0){
     const std::vector<const reco::GenParticle*> m=GetMothers(boson);
    int q(0),qbar(0);
    TLorentzVector Z(0,0,0,0);
    for(unsigned int i=0;i<m.size();i++){
      if(m.at(i)->pdgId()==PdtPdgMini::d      || m.at(i)->pdgId()==PdtPdgMini::u      ){q++;}
      if(m.at(i)->pdgId()==PdtPdgMini::anti_d || m.at(i)->pdgId()==PdtPdgMini::anti_u ){qbar++;}
    }
    if(q==1 && qbar==1){// assume q has largest E (valence vs see quarks) 
      if(taum.Vect().Dot(Zboson.Vect())/(Zboson.P()*taum.P())>0){
        if(abs(boson->pdgId())==PdtPdgMini::Z0)      TauSpinEffectsZ_Xf->Fill(x1,weight);
        if(abs(boson->pdgId())==PdtPdgMini::Higgs0) TauSpinEffectsH_Xf->Fill(x1,weight);
      }
      else{
        if(abs(boson->pdgId())==PdtPdgMini::Z0)      TauSpinEffectsZ_Xb->Fill(x1,weight);
        if(abs(boson->pdgId())==PdtPdgMini::Higgs0) TauSpinEffectsH_Xb->Fill(x1,weight);
      }
    }
      }
    }
  }
}

int TauValidation::tauDecayChannel ( const reco::GenParticle *  tau, int  jak_id, unsigned int  TauBitMask, double  weight  )

private

Definition at line 354 of file TauValidation.cc.

References a1, countParticles(), electron, MonitorElement::Fill(), K, Kstar, muon, pi, pi1pi0, pinpi0, rho, stable, reco::LeafCandidate::status(), TauDecayChannels, tripi, tripinpi0, and undetermined.

Referenced by analyze().

                                                                                                               {
  int channel = undetermined; 
  if(tau->status() == 1) channel = stable; 
  int allCount   = 0, 
    eCount     = 0, 
    muCount    = 0, 
    pi0Count   = 0, 
    piCount    = 0, 
    rhoCount   = 0, 
    a1Count    = 0, 
    KCount     = 0, 
    KstarCount = 0; 

  countParticles(tau,allCount,eCount,muCount,pi0Count,piCount,rhoCount,a1Count,KCount,KstarCount);

  // resonances   
  if(KCount >= 1)     channel = K; 
  if(KstarCount >= 1) channel = Kstar; 
  if(a1Count >= 1)    channel = a1; 
  if(rhoCount >= 1)   channel = rho; 
  if(channel!=undetermined && weight!=0.0) TauDecayChannels->Fill(channel,weight); 
  
  // final state products 
  if(piCount == 1 && pi0Count == 0) channel = pi; 
  if(piCount == 1 && pi0Count == 1) channel = pi1pi0; 
  if(piCount == 1 && pi0Count > 1)  channel = pinpi0; 
  if(piCount == 3 && pi0Count == 0) channel = tripi; 
  if(piCount == 3 && pi0Count > 0)  channel = tripinpi0; 
  if(eCount == 1)                   channel = electron; 
  if(muCount == 1)                  channel = muon; 
  if(weight!=0.0) TauDecayChannels->Fill(channel,weight); 
  return channel; 
} 

int TauValidation::tauMother ( const reco::GenParticle *  tau, double  weight  )

private

Definition at line 334 of file TauValidation.cc.

References A0, funct::abs(), B, D, MonitorElement::Fill(), findMother(), gamma, H0, Hpm, HSM, diffTwoXMLs::label, other, reco::LeafCandidate::pdgId(), TauMothers, W, and Z.

Referenced by analyze().

                                                                     {
  if(abs(tau->pdgId()) != 15 ) return -3;
  int mother_pid = findMother(tau);
  if(mother_pid == -2) return -2;
  int label = other;
  if(abs(mother_pid) == 24) label = W;
  if(abs(mother_pid) == 23) label = Z;
  if(abs(mother_pid) == 22) label = gamma;
  if(abs(mother_pid) == 25) label = HSM;
  if(abs(mother_pid) == 35) label = H0;
  if(abs(mother_pid) == 36) label = A0;
  if(abs(mother_pid) == 37) label = Hpm;
  int mother_shortpid=(abs(mother_pid)%10000);
  if(mother_shortpid>500 && mother_shortpid<600 )label = B;
  if(mother_shortpid>400 && mother_shortpid<500)label = D;
  TauMothers->Fill(label,weight);
  if(label==B || label == D || label == other) return -1;
  return mother_pid;
}

int TauValidation::tauProngs ( const reco::GenParticle *  , double  weight  )

private

double TauValidation::visibleTauEnergy ( const reco::GenParticle *  tau )

private

Definition at line 699 of file TauValidation.cc.

References funct::abs(), reco::CompositeRefCandidateT< D >::daughter(), reco::LeafCandidate::energy(), i, reco::CompositeRefCandidateT< D >::numberOfDaughters(), p4, reco::LeafCandidate::pdgId(), sysUtil::pid, reco::LeafCandidate::px(), reco::LeafCandidate::py(), and reco::LeafCandidate::pz().

                                                                {
  TLorentzVector p4(tau->px(),tau->py(),tau->pz(),tau->energy());
  for(unsigned int i = 0; i <tau->numberOfDaughters(); i++){
    const reco::GenParticle *dau=static_cast<const reco::GenParticle*>(tau->daughter(i));
    int pid = dau->pdgId();
    if(abs(pid) == 15) return visibleTauEnergy(dau);
    if(abs(pid) == 12 || abs(pid) == 14 || abs(pid) == 16) {
      p4-=TLorentzVector(dau->px(),dau->py(),dau->pz(),dau->energy());
    }
  }
  return p4.E();
}

double TauValidation::Zstoa ( double  zs )

private

Definition at line 669 of file TauValidation.cc.

References a, and mathSSE::sqrt().

Referenced by spinEffectsZH().

                                    {
  double a=1-sqrt(fabs(1.0-2*fabs(zs)));
  if(zs<0){
    a*=-1.0;
  }
  return a;
}

Member Data Documentation

MonitorElement * TauValidation::DecayLength

private

Definition at line 92 of file TauValidation.h.

Referenced by analyze(), and bookHistograms().

edm::ESHandle<HepPDT::ParticleDataTable> TauValidation::fPDGTable

private

PDT table.

Definition at line 89 of file TauValidation.h.

Referenced by dqmBeginRun(), and spinEffectsZH().

edm::InputTag TauValidation::genparticleCollection_

private

Definition at line 86 of file TauValidation.h.

Referenced by TauValidation().

edm::EDGetTokenT<reco::GenParticleCollection> TauValidation::genparticleCollectionToken_

private

Definition at line 111 of file TauValidation.h.

Referenced by analyze(), and TauValidation().

MonitorElement * TauValidation::LifeTime

private

Definition at line 92 of file TauValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TauValidation::MODEID

private

Definition at line 105 of file TauValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<std::vector<MonitorElement *> > TauValidation::MODEInvMass

private

Definition at line 106 of file TauValidation.h.

Referenced by analyze(), and bookHistograms().

unsigned int TauValidation::NMODEID

private

Definition at line 104 of file TauValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * TauValidation::nPrimeTaus

private

Definition at line 91 of file TauValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TauValidation::nTaus

private

Definition at line 91 of file TauValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * TauValidation::TauBremPhotonsN

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and photons().

MonitorElement * TauValidation::TauBremPhotonsPt

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and photons().

MonitorElement * TauValidation::TauBremPhotonsPtSum

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and photons().

MonitorElement * TauValidation::TauDecayChannels

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and tauDecayChannel().

MonitorElement * TauValidation::TauEta

private

Definition at line 92 of file TauValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * TauValidation::TauFSRPhotonsN

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and photons().

MonitorElement * TauValidation::TauFSRPhotonsPt

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and photons().

MonitorElement * TauValidation::TauFSRPhotonsPtSum

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and photons().

MonitorElement * TauValidation::TauMothers

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and tauMother().

MonitorElement * TauValidation::TauPhi

private

Definition at line 92 of file TauValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * TauValidation::TauProngs

private

Definition at line 92 of file TauValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TauValidation::TauPt

private

Definition at line 92 of file TauValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * TauValidation::TauSpinEffectsH_eX

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_muX

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_MVis

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_pipiAcollinearity

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_pipiAcollinearityzoom

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_pipiAcoplanarity

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_rhorhoAcoplanarityminus

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_rhorhoAcoplanarityplus

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_X

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_Xb

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_Xf

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_Zs

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsHpm_eX

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsHpm_muX

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsHpm_UpsilonA1

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsHpm_UpsilonRho

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsHpm_X

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsW_eX

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsW_muX

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsW_UpsilonA1

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsW_UpsilonRho

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsW_X

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsZ_eX

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_muX

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_MVis

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_X

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_X100to120

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_X120UP

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_X50to75

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_X75to88

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_X88to100

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_Xb

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_Xf

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_Zs

private

Definition at line 92 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

int TauValidation::zsbins

private

Definition at line 108 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

double TauValidation::zsmax

private

Definition at line 109 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().

double TauValidation::zsmin

private

Definition at line 109 of file TauValidation.h.

Referenced by bookHistograms(), and spinEffectsZH().