#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::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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 DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private 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 edm::one::EDProducerBase
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)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

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

29 {}

Member Function Documentation

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

overridevirtual

Gathering the reco::GenParticleCollection information

Reimplemented from DQMEDAnalyzer.

Definition at line 156 of file TauValidation.cc.

References funct::abs(), TauDecay_GenParticle::AnalyzeTau(), 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 &
	)

overridevirtual

Setting the DQM top directories

Implements DQMEDAnalyzer.

Definition at line 35 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 388 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++;
     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
	)

override

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(), 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 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(), 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 285 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 255 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 264 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 278 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 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(), 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 694 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 712 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 408 of file TauValidation.cc.

References a1, funct::abs(), 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().

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

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation