#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 Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override

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

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

	TauValidation (const edm::ParameterSet &)

	~TauValidation () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

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

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

override

Definition at line 28 of file TauValidation.cc.

28 {}

Member Function Documentation

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

override

Gathering the reco::GenParticleCollection information

Definition at line 368 of file TauValidation.cc.

References funct::abs(), TauDecay_GenParticle::AnalyzeTau(), pfBoostedDoubleSVAK8TagInfos_cfi::beta, EnergyCorrector::c, DecayLength, MonitorElement::Fill(), genparticleCollectionToken_, GenHFHadronMatcher_cfi::genParticles, TauDecay_GenParticle::Get_TauDecayProducts(), edm::Event::getByToken(), PdtPdgMini::Higgs0, mps_fire::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, nanoDQM_cfi::PV, spinEffectsWHpm(), spinEffectsZH(), nanoDQM_cfi::SV, tauDecayChannel(), TauEta, tauMother(), TauPhi, TauProngs, TauPt, mps_merge::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 &
	)

override

Setting the DQM top directories

Definition at line 32 of file TauValidation.cc.

References A0, a1, B, DQMHelper::book1dHisto(), D, DecayLength, TauDecay::DecayMode(), electron, gamma, H0, Hpm, HSM, K, Kstar, LifeTime, TauDecay::MODE_3PI, TauDecay::MODE_KPIK, TauDecay::MODE_KPIPI, TauDecay::MODE_PI2PI0, MODEID, MODEInvMass, 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, 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 636 of file TauValidation.cc.

References funct::abs(), reco::CompositeRefCandidateT< D >::daughter(), mps_fire::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++;
     else if (abs(pid) == 13)
       muCount++;
     else if (abs(pid) == 111)
       pi0Count++;
     else if (abs(pid) == 211)
       piCount++;
     else if (abs(pid) == 213)
       rhoCount++;
     else if (abs(pid) == 20213)
       a1Count++;
     else if (abs(pid) == 321)
       KCount++;
     else 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
	)

override

Definition at line 30 of file TauValidation.cc.

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

30 { c.getData(fPDGTable); }

edm::EventSetup::getData

bool getData(T &iHolder) const

Definition: EventSetup.h:111

TauValidation::fPDGTable

edm::ESHandle< HepPDT::ParticleDataTable > fPDGTable

PDT table.

Definition: TauValidation.h:89

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 506 of file TauValidation.cc.

References funct::abs(), reco::CompositeRefCandidateT< D >::daughter(), mps_fire::i, isLastTauinChain(), reco::CompositeRefCandidateT< D >::numberOfDaughters(), common_cff::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 486 of file TauValidation.cc.

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

Referenced by spinEffectsZH(), and tauMother().

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

reco::LeafCandidate::pdgId

int pdgId() const final

PDG identifier.

Definition: LeafCandidate.h:161

TauValidation::GetMother

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

Definition: TauValidation.cc:465

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

private

Definition at line 533 of file TauValidation.cc.

References funct::abs(), reco::CompositeRefCandidateT< D >::daughter(), GetMother(), mps_fire::i, funct::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 496 of file TauValidation.cc.

References mps_fire::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 465 of file TauValidation.cc.

References mps_fire::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 475 of file TauValidation.cc.

References mps_fire::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 488 of file TauValidation.cc.

References reco::CompositeRefCandidateT< D >::daughter(), mps_fire::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 1022 of file TauValidation.cc.

References leadingPionP4().

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

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

private

Definition at line 1026 of file TauValidation.cc.

References funct::abs(), reco::CompositeRefCandidateT< D >::daughter(), reco::LeafCandidate::energy(), mps_fire::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 1041 of file TauValidation.cc.

References reco::LeafCandidate::energy(), GetMother(), funct::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 1060 of file TauValidation.cc.

References MonitorElement::Fill(), findFSRandBrem(), FindPhotosFSR(), findTauList(), mps_fire::i, TauBremPhotonsN, TauBremPhotonsPt, TauBremPhotonsPtSum, TauFSRPhotonsN, TauFSRPhotonsPt, TauFSRPhotonsPtSum, and mps_merge::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.empty()) {
     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 670 of file TauValidation.cc.

References a1, funct::abs(), randomXiThetaGunProducer_cfi::energy, MonitorElement::Fill(), gamma, mps_fire::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 mps_merge::weight.

Referenced by analyze().

                                                                                                                 {
   // polarization only for 1-prong hadronic taus with no neutral pions
   if (decay == TauDecay::MODE_PION || decay == TauDecay::MODE_MUON || decay == TauDecay::MODE_ELECTRON) {
     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);
       else if (abs(mother) == 37)
         TauSpinEffectsHpm_X->Fill(energy, weight);
     }
     else if (decay == TauDecay::MODE_MUON) {
       if (abs(mother) == 24)
         TauSpinEffectsW_muX->Fill(energy, weight);
       else if (abs(mother) == 37)
         TauSpinEffectsHpm_muX->Fill(energy, weight);
     }
     else if (decay == TauDecay::MODE_ELECTRON) {
       if (abs(mother) == 24)
         TauSpinEffectsW_eX->Fill(energy, weight);
       else 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;
       }
       else if (abs(part.at(i)->pdgId()) == PdtPdgMini::pi0) {
         rho += LV;
       }
     }
     if (abs(mother) == 24)
       TauSpinEffectsW_UpsilonRho->Fill(2 * pi.P() / rho.P() - 1, weight);
     else 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++;
       }
       else 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;
     else 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);
     else if (abs(mother) == 37)
       TauSpinEffectsHpm_UpsilonA1->Fill(gamma, weight);
   }
 }

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

private

Definition at line 744 of file TauValidation.cc.

References funct::abs(), TauDecay_GenParticle::AnalyzeTau(), PdtPdgMini::anti_d, PdtPdgMini::anti_u, ALCARECOTkAlJpsiMuMu_cff::charge, PdtPdgMini::d, reco::CompositeRefCandidateT< D >::daughter(), reco::LeafCandidate::energy(), MonitorElement::Fill(), findMother(), fPDGTable, TauDecay_GenParticle::Get_TauDecayProducts(), GetMothers(), PdtPdgMini::Higgs0, mps_fire::i, createfilelist::int, leadingPionP4(), funct::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, mps_merge::weight, globals_cff::x1, globals_cff::x2, 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 595 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 563 of file TauValidation.cc.

References A0, funct::abs(), B, D, MonitorElement::Fill(), findMother(), gamma, H0, Hpm, HSM, 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 1046 of file TauValidation.cc.

References funct::abs(), reco::CompositeRefCandidateT< D >::daughter(), reco::LeafCandidate::energy(), mps_fire::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 1014 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().

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation