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

typedef WorkerT< EDAnalyzer >	WorkerType

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

virtual void	beginJob ()

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

virtual void	endJob ()

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

	TauValidation (const edm::ParameterSet &)

virtual	~TauValidation ()

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
void	findFSRandBrem (const HepMC::GenParticle p, bool doBrem, std::vector< const HepMC::GenParticle > &ListofFSR, std::vector< const HepMC::GenParticle * > &ListofBrem)

int	findMother (const HepMC::GenParticle *)

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

void	findTauList (const HepMC::GenParticle tau, std::vector< const HepMC::GenParticle > &TauList)

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

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

bool	isLastTauinChain (const HepMC::GenParticle *tau)

double	leadingPionMomentum (const HepMC::GenParticle *, double weight)

TLorentzVector	leadingPionP4 (const HepMC::GenParticle *)

TLorentzVector	motherP4 (const HepMC::GenParticle *)

void	photons (const HepMC::GenParticle *, double weight)

void	rtau (const HepMC::GenParticle *, int, int, double weight)

void	spinEffects (const HepMC::GenParticle , int, int, std::vector< HepMC::GenParticle > &part, double weight)

void	spinEffectsZ (const HepMC::GenParticle *boson, double weight)

int	tauDecayChannel (const HepMC::GenParticle *, double weight=0.0)

int	tauMother (const HepMC::GenParticle *, double weight)

int	tauProngs (const HepMC::GenParticle *, double weight)

double	visibleTauEnergy (const HepMC::GenParticle *)

double	Zstoa (double zs)

Private Attributes
WeightManager	_wmanager

DQMStore *	dbe
	ME's "container". More...

edm::ESHandle < HepPDT::ParticleDataTable >	fPDGTable
	PDT table. More...

edm::InputTag	hepmcCollection_

MonitorElement *	JAKID

std::vector< std::vector < MonitorElement * > >	JAKInvMass

unsigned int	NJAKID

MonitorElement *	nPrimeTaus

MonitorElement *	nTaus

MonitorElement *	TauBremPhotonsN

MonitorElement *	TauBremPhotonsPt

MonitorElement *	TauBremPhotonsPtSum

MonitorElement *	TauDecayChannels

MonitorElement *	TauEta

double	tauEtCut

MonitorElement *	TauFSRPhotonsN

MonitorElement *	TauFSRPhotonsPt

MonitorElement *	TauFSRPhotonsPtSum

MonitorElement *	TauMothers

MonitorElement *	TauPhi

MonitorElement *	TauProngs

MonitorElement *	TauPt

MonitorElement *	TauRtauHpm

MonitorElement *	TauRtauW

MonitorElement *	TauSpinEffectsH_eX

MonitorElement *	TauSpinEffectsH_muX

MonitorElement *	TauSpinEffectsH_MVis

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_Xb

MonitorElement *	TauSpinEffectsZ_Xf

MonitorElement *	TauSpinEffectsZ_Zs

int	zsbins

double	zsmax

double	zsmin

Additional Inherited Members
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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

CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 36 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 40 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 54 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 20 of file TauValidation.cc.

References dbe, and cppFunctionSkipper::operator.

                                                         : 
   _wmanager(iPSet)
   ,hepmcCollection_(iPSet.getParameter<edm::InputTag>("hepmcCollection"))
   ,tauEtCut(iPSet.getParameter<double>("tauEtCutForRtau"))
   ,NJAKID(22)
   ,zsbins(20)
   ,zsmin(-0.5)
   ,zsmax(0.5)
 {    
   dbe = 0;
   dbe = edm::Service<DQMStore>().operator->();
 }

TauValidation::~TauValidation ( )

virtual

Definition at line 33 of file TauValidation.cc.

33 {}

Member Function Documentation

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

virtual

Gathering the HepMCProduct information

Implements edm::EDAnalyzer.

Definition at line 158 of file TauValidation.cc.

References _wmanager, abs, TauDecay_CMSSW::AnalyzeTau(), MonitorElement::Fill(), TauDecay_CMSSW::Get_TauDecayProducts(), edm::Event::getByLabel(), hepmcCollection_, i, isLastTauinChain(), TauDecay::isTauFinalStateParticle(), TauDecay::JAK_A1_3PI, TauDecay::JAK_KPIK, TauDecay::JAK_KPIPI, JAKID, JAKInvMass, PdtPdgMini::K_plus, NJAKID, nPrimeTaus, TauDecay::nProng(), nTaus, PdtPdgMini::nu_e, PdtPdgMini::nu_mu, PdtPdgMini::nu_tau, photons(), PdtPdgMini::pi_plus, rtau(), spinEffects(), spinEffectsZ(), tauDecayChannel(), TauEta, tauMother(), TauPhi, tauProngs(), TauPt, WeightManager::weight(), and histoStyle::weight.

 { 
   edm::Handle<HepMCProduct> evt;
   iEvent.getByLabel(hepmcCollection_, evt);
 
   //Get EVENT
   HepMC::GenEvent *myGenEvent = new HepMC::GenEvent(*(evt->GetEvent()));
 
   double weight = _wmanager.weight(iEvent);
 
   /*
   edm::Handle<double> WT;
   iEvent.getByLabel(edm::InputTag("TauSpinnerGen","TauSpinnerWT"),WT);
   weight = 1.0;
   if(*(WT.product())>1e-3 && *(WT.product())<=10.0) weight=(*(WT.product()));
   else {weight=1.0;}
   */
 
   // find taus
   for(HepMC::GenEvent::particle_const_iterator iter = myGenEvent->particles_begin(); iter != myGenEvent->particles_end(); iter++) {
     if(abs((*iter)->pdg_id())==23){
       spinEffectsZ(*iter,weight);
     }
     if(abs((*iter)->pdg_id())==15){
       if(isLastTauinChain(*iter)){
     nTaus->Fill(0.5,weight);
     int mother  = tauMother(*iter,weight);
     int decaychannel = tauDecayChannel(*iter,weight);
         tauProngs(*iter, weight);
     if(mother>-1){ // exclude B, D and other non-signal decay modes
       nPrimeTaus->Fill(0.5,weight);
       TauPt->Fill((*iter)->momentum().perp(),weight);
       TauEta->Fill((*iter)->momentum().eta(),weight);
       TauPhi->Fill((*iter)->momentum().phi(),weight);
       rtau(*iter,mother,decaychannel,weight);
       photons(*iter,weight);
     }
     //Adding JAKID and Mass information
     //
         TauDecay_CMSSW TD;
         unsigned int jak_id, TauBitMask;
         if(TD.AnalyzeTau((*iter),jak_id,TauBitMask,false,false)){
       JAKID->Fill(jak_id,weight);
       if(jak_id<=NJAKID){
         int tcharge=(*iter)->pdg_id()/abs((*iter)->pdg_id());
         std::vector<HepMC::GenParticle*> part=TD.Get_TauDecayProducts();
         spinEffects(*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;
         for(unsigned int i=0;i<part.size();i++){
           if(TD.isTauFinalStateParticle(part.at(i)->pdg_id()) &&
          abs(part.at(i)->pdg_id())!=PdtPdgMini::nu_e &&
          abs(part.at(i)->pdg_id())!=PdtPdgMini::nu_mu &&
          abs(part.at(i)->pdg_id())!=PdtPdgMini::nu_tau ){
         TLorentzVector LV(part.at(i)->momentum().px(),part.at(i)->momentum().py(),part.at(i)->momentum().pz(),part.at(i)->momentum().e());
         LVQ+=LV;
         if(jak_id==TauDecay::JAK_A1_3PI ||
            jak_id==TauDecay::JAK_KPIK ||
            jak_id==TauDecay::JAK_KPIPI
            ){
           if((tcharge==part.at(i)->pdg_id()/abs(part.at(i)->pdg_id()) && TD.nProng(TauBitMask)==3) || (jak_id==TauDecay::JAK_A1_3PI && TD.nProng(TauBitMask)==1 && abs(part.at(i)->pdg_id())==PdtPdgMini::pi_plus) ){
             LVS13+=LV;
             LVS23+=LV;
           }
           else{
             LVS12+=LV;
             if(!haspart1 && ((jak_id==TauDecay::JAK_A1_3PI)  || (jak_id!=TauDecay::JAK_A1_3PI && abs(part.at(i)->pdg_id())==PdtPdgMini::K_plus) )){
               LVS13+=LV;
               haspart1=true;
             }
             else{
               LVS23+=LV;
             }
           }
         }
           }
         }
         part.clear();
         JAKInvMass.at(jak_id).at(0)->Fill(LVQ.M(),weight);
         if(jak_id==TauDecay::JAK_A1_3PI ||
            jak_id==TauDecay::JAK_KPIK ||
            jak_id==TauDecay::JAK_KPIPI
            ){
           JAKInvMass.at(jak_id).at(1)->Fill(LVS13.M(),weight);
           JAKInvMass.at(jak_id).at(2)->Fill(LVS23.M(),weight);
           JAKInvMass.at(jak_id).at(3)->Fill(LVS12.M(),weight);
         }
       }
     }
     else{
       JAKID->Fill(jak_id,weight);  
     }
       }
     }
   }
   delete myGenEvent;
 }//analyze

void TauValidation::beginJob ( void )

virtual

Setting the DQM top directories

Reimplemented from edm::EDAnalyzer.

Definition at line 35 of file TauValidation.cc.

References A0, a1, DQMStore::book1D(), D, dbe, electron, gamma, H0, Hpm, HSM, i, TauDecay::JAK_A1_3PI, TauDecay::JAK_KPIK, TauDecay::JAK_KPIPI, JAKID, JAKInvMass, K, Kstar, muon, NJAKID, nPrimeTaus, nTaus, other, pi, pi1pi0, pinpi0, rho, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), stable, TauBremPhotonsN, TauBremPhotonsPt, TauBremPhotonsPtSum, TauDecayChannels, TauEta, TauFSRPhotonsN, TauFSRPhotonsPt, TauFSRPhotonsPtSum, TauMothers, TauPhi, TauProngs, TauPt, TauRtauHpm, TauRtauW, TauSpinEffectsH_eX, TauSpinEffectsH_muX, TauSpinEffectsH_MVis, 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_Xb, TauSpinEffectsZ_Xf, TauSpinEffectsZ_Zs, tmp, tripi, tripinpi0, undetermined, W, Z, zsbins, zsmax, and zsmin.

 {
   if(dbe){
     dbe->setCurrentFolder("Generator/Tau");
     
     // Number of analyzed events
     nTaus = dbe->book1D("nTaus", "n analyzed Taus", 1, 0., 1.);
     nPrimeTaus = dbe->book1D("nPrimeTaus", "n analyzed prime Taus", 1, 0., 1.);
 
     //Kinematics
     TauPt            = dbe->book1D("TauPt","Tau pT", 100 ,0,100);
     TauEta           = dbe->book1D("TauEta","Tau eta", 100 ,-2.5,2.5);
     TauPhi           = dbe->book1D("TauPhi","Tau phi", 100 ,-3.14,3.14);
     TauProngs        = dbe->book1D("TauProngs","Tau n prongs", 7 ,0,7);
     TauDecayChannels = dbe->book1D("TauDecayChannels","Tau decay channels", 13 ,0,13);
     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        = dbe->book1D("TauMothers","Tau mother particles", 10 ,0,10);
     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}");
 
     TauRtauW          = dbe->book1D("TauRtauW","W->Tau p(leading track)/E(visible tau)", 50 ,0,1);     TauRtauW->setAxisTitle("rtau");
     TauRtauHpm        = dbe->book1D("TauRtauHpm","Hpm->Tau p(leading track)/E(visible tau)", 50 ,0,1); TauRtauHpm->setAxisTitle("rtau");
 
     TauSpinEffectsW_X   = dbe->book1D("TauSpinEffectsWX","Pion energy in W rest frame", 50 ,0,1);     TauSpinEffectsW_X->setAxisTitle("X");
     TauSpinEffectsHpm_X = dbe->book1D("TauSpinEffectsHpmX","Pion energy in Hpm rest frame", 50 ,0,1); TauSpinEffectsHpm_X->setAxisTitle("X");
 
     TauSpinEffectsW_eX   = dbe->book1D("TauSpinEffectsWeX","e energy in W rest frame", 50 ,0,1);     TauSpinEffectsW_eX->setAxisTitle("X");
     TauSpinEffectsHpm_eX = dbe->book1D("TauSpinEffectsHpmeX","e energy in Hpm rest frame", 50 ,0,1); TauSpinEffectsHpm_eX->setAxisTitle("X");
 
     TauSpinEffectsW_muX   = dbe->book1D("TauSpinEffectsWmuX","mu energy in W rest frame", 50 ,0,1);     TauSpinEffectsW_muX->setAxisTitle("X");
     TauSpinEffectsHpm_muX = dbe->book1D("TauSpinEffectsHpmmuX","mu energy in Hpm rest frame", 50 ,0,1); TauSpinEffectsHpm_muX->setAxisTitle("X");
 
     TauSpinEffectsW_UpsilonRho   = dbe->book1D("TauSpinEffectsWUpsilonRho","#Upsilon for #rho", 50 ,-1,1);     TauSpinEffectsW_UpsilonRho->setAxisTitle("#Upsilon");
     TauSpinEffectsHpm_UpsilonRho = dbe->book1D("TauSpinEffectsHpmUpsilonRho","#Upsilon for #rho", 50 ,-1,1);   TauSpinEffectsHpm_UpsilonRho->setAxisTitle("#Upsilon");
 
     TauSpinEffectsW_UpsilonA1   = dbe->book1D("TauSpinEffectsWUpsilonA1","#Upsilon for a1", 50 ,-1,1);       TauSpinEffectsW_UpsilonA1->setAxisTitle("#Upsilon");
     TauSpinEffectsHpm_UpsilonA1 = dbe->book1D("TauSpinEffectsHpmUpsilonA1","#Upsilon for a1", 50 ,-1,1);     TauSpinEffectsHpm_UpsilonA1->setAxisTitle("#Upsilon");
 
     TauSpinEffectsZ_MVis   = dbe->book1D("TauSpinEffectsZMVis","Mass of pi+ pi-", 25 ,0,1.1);       TauSpinEffectsZ_MVis->setAxisTitle("M_{#pi^{+}#pi^{-}}");
     TauSpinEffectsH_MVis   = dbe->book1D("TauSpinEffectsHMVis","Mass of pi+ pi-", 25 ,0,1.1);       TauSpinEffectsZ_MVis->setAxisTitle("M_{#pi^{+}#pi^{-}}");
 
     TauSpinEffectsZ_Zs   = dbe->book1D("TauSpinEffectsZZs","Z_{s}", zsbins ,zsmin,zsmax);        TauSpinEffectsZ_Zs->setAxisTitle("Z_{s}");
     TauSpinEffectsH_Zs   = dbe->book1D("TauSpinEffectsHZs","Z_{s}", zsbins ,zsmin,zsmax);        TauSpinEffectsZ_Zs->setAxisTitle("Z_{s}");
 
     TauSpinEffectsZ_Xf   = dbe->book1D("TauSpinEffectsZXf","X of forward emitted #tau^{-}", 25 ,0,1.0);           TauSpinEffectsZ_Xf->setAxisTitle("X_{f}");
     TauSpinEffectsH_Xf   = dbe->book1D("TauSpinEffectsHXf","X of forward emitted #tau^{-}", 25 ,0,1.0);           TauSpinEffectsZ_Xf->setAxisTitle("X_{f}");
 
     TauSpinEffectsZ_Xb   = dbe->book1D("TauSpinEffectsZXb","X of backward emitted #tau^{-}", 25 ,0,1.0);           TauSpinEffectsZ_Xb->setAxisTitle("X_{b}");
     TauSpinEffectsH_Xb   = dbe->book1D("TauSpinEffectsHXb","X of backward emitted #tau^{-}", 25 ,0,1.0);           TauSpinEffectsZ_Xb->setAxisTitle("X_{b}");
 
     TauSpinEffectsZ_eX   = dbe->book1D("TauSpinEffectsZeX","e energy in Z rest frame", 50 ,0,1);     TauSpinEffectsZ_eX->setAxisTitle("X");
     TauSpinEffectsH_eX = dbe->book1D("TauSpinEffectsHeX","e energy in H rest frame", 50 ,0,1); TauSpinEffectsH_eX->setAxisTitle("X");
 
     TauSpinEffectsZ_muX   = dbe->book1D("TauSpinEffectsZmuX","mu energy in z rest frame", 50 ,0,1);     TauSpinEffectsZ_muX->setAxisTitle("X");
     TauSpinEffectsH_muX = dbe->book1D("TauSpinEffectsHmuX","mu energy in H rest frame", 50 ,0,1); TauSpinEffectsH_muX->setAxisTitle("X");
 
 
     TauFSRPhotonsN=dbe->book1D("TauFSRPhotonsN","FSR Photons radiating from/with tau (Gauge Boson)", 5 ,-0.5,4.5);
     TauFSRPhotonsN->setAxisTitle("N FSR Photons radiating from/with tau");
     TauFSRPhotonsPt=dbe->book1D("TauFSRPhotonsPt","Pt of FSR Photons radiating from/with tau (Gauge Boson)", 100 ,0,100);
     TauFSRPhotonsPt->setAxisTitle("P_{t} of FSR Photons radiating from/with tau [per tau]");
     TauFSRPhotonsPtSum=dbe->book1D("TauFSRPhotonsPtSum","Pt of FSR Photons radiating from/with tau (Gauge Boson)", 100 ,0,100);
     TauFSRPhotonsPtSum->setAxisTitle("P_{t} of FSR Photons radiating from/with tau [per tau]");
 
     TauBremPhotonsN=dbe->book1D("TauBremPhotonsN","Brem. Photons radiating in tau decay", 5 ,-0.5,4.5);
     TauBremPhotonsN->setAxisTitle("N FSR Photons radiating from/with tau");
     TauBremPhotonsPt=dbe->book1D("TauBremPhotonsPt","Sum Brem Pt ", 100 ,0,100);
     TauFSRPhotonsPt->setAxisTitle("P_{t} of Brem. Photons radiating in tau decay");    
     TauBremPhotonsPtSum =dbe->book1D("TauBremPhotonsPtSum","Sum of Brem Pt ", 100 ,0,100);
     TauFSRPhotonsPtSum->setAxisTitle("Sum P_{t} of Brem. Photons radiating in tau decay");
 
     JAKID =dbe->book1D("JAKID","JAK ID",NJAKID+1,-0.5,NJAKID+0.5);
     for(unsigned int i=0; i<NJAKID+1;i++){
       JAKInvMass.push_back(std::vector<MonitorElement *>());
       TString tmp="JAKID";
       tmp+=i;
       JAKInvMass.at(i).push_back(dbe->book1D("M"+tmp,"M_{"+tmp+"} (GeV)", 80 ,0,2.0));
       if(i==TauDecay::JAK_A1_3PI ||
      i==TauDecay::JAK_KPIK ||
      i==TauDecay::JAK_KPIPI ){
     JAKInvMass.at(i).push_back(dbe->book1D("M13"+tmp,"M_{13,"+tmp+"} (GeV)", 80 ,0,2.0));
     JAKInvMass.at(i).push_back(dbe->book1D("M23"+tmp,"M_{23,"+tmp+"} (GeV)", 80 ,0,2.0));
     JAKInvMass.at(i).push_back(dbe->book1D("M12"+tmp,"M_{12,"+tmp+"} (GeV)", 80 ,0,2.0));
       }
     }
   }
   
   return;
 }

void TauValidation::beginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

virtual

Get PDT Table

Reimplemented from edm::EDAnalyzer.

Definition at line 151 of file TauValidation.cc.

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

 {
   iSetup.getData( fPDGTable );
   return;
 }

void TauValidation::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 147 of file TauValidation.cc.

                           {
   return;
 }

void TauValidation::endRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 157 of file TauValidation.cc.

157 {return;}

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

private

Definition at line 325 of file TauValidation.cc.

References abs, and isLastTauinChain().

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->pdg_id())==15){
     if(isLastTauinChain(p)){ doBrem=true;}
     else{ doBrem=false;}
   }
     int photo_ID=22;
   if ( p->end_vertex() ) {
     HepMC::GenVertex::particle_iterator dau;
     for (dau = p->end_vertex()->particles_begin(HepMC::children); dau!= p->end_vertex()->particles_end(HepMC::children); dau++ ) {
       //if(doBrem) std::cout << "true " << (*dau)->pdg_id() << " "  << std::endl;
       //else std::cout << "false " << (*dau)->pdg_id() << " " << std::endl;
       if(abs((*dau)->pdg_id()) == abs(photo_ID) && !doBrem){ListofFSR.push_back(*dau);}
       if(abs((*dau)->pdg_id()) == abs(photo_ID) && doBrem){ListofBrem.push_back(*dau);}
       if((*dau)->end_vertex() && (*dau)->end_vertex()->particles_out_size()>0 && abs((*dau)->pdg_id()) != 111 && abs((*dau)->pdg_id()) != 221/* remove pi0 and eta decays*/){
     findFSRandBrem(*dau,doBrem,ListofFSR,ListofBrem);
       }
     }
   }
 }

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

private

Definition at line 288 of file TauValidation.cc.

References parents.

Referenced by spinEffectsZ(), and tauMother().

                                                         {
   int mother_pid = 0;
   if ( tau->production_vertex() ) {
     HepMC::GenVertex::particle_iterator mother;
     for (mother = tau->production_vertex()->particles_begin(HepMC::parents); mother!= tau->production_vertex()->particles_end(HepMC::parents); mother++ ) {
       mother_pid = (*mother)->pdg_id();
       if(mother_pid == tau->pdg_id()) return findMother(*mother); //mother_pid = -1; Make recursive to look for last tau in chain
     }
   }
   return mother_pid;
 }

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

private

Definition at line 349 of file TauValidation.cc.

References abs, configurableAnalysis::GenParticle, GetMother(), and m.

Referenced by photons().

                                                                                                                            {
   BosonScale=0.0;
   const HepMC::GenParticle* m=GetMother(p);
   double mother_pid=m->pdg_id();
   if(m->end_vertex() && mother_pid!=p->pdg_id()){
     HepMC::GenVertex::particle_iterator dau;
     for (dau = m->end_vertex()->particles_begin(HepMC::children); dau!= m->end_vertex()->particles_end(HepMC::children); dau++ ) {
       int dau_pid = (*dau)->pdg_id();
       if(fabs(dau_pid) == 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 HepMC::GenParticle *	tau,
		std::vector< const HepMC::GenParticle * > &	TauList
	)

private

Definition at line 313 of file TauValidation.cc.

References parents, and metsig::tau.

Referenced by photons().

                                                                                                       {
   TauList.insert(TauList.begin(),tau);
   if ( tau->production_vertex() ) {
     HepMC::GenVertex::particle_iterator mother;
     for (mother = tau->production_vertex()->particles_begin(HepMC::parents); mother!= tau->production_vertex()->particles_end(HepMC::parents);mother++) {
       if((*mother)->pdg_id() == tau->pdg_id()){
     findTauList(*mother,TauList);
       }
     }
   }
 }

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

private

Definition at line 263 of file TauValidation.cc.

References parents, and metsig::tau.

Referenced by FindPhotosFSR(), and motherP4().

                                                                            {
   if ( tau->production_vertex() ) {
     HepMC::GenVertex::particle_iterator mother;
     for (mother = tau->production_vertex()->particles_begin(HepMC::parents); mother!= tau->production_vertex()->particles_end(HepMC::parents); mother++ ) {
       if((*mother)->pdg_id() == tau->pdg_id()) return GetMother(*mother);
       return (*mother);
     }
   }
   return tau;
 }

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

private

Definition at line 275 of file TauValidation.cc.

References parents.

Referenced by spinEffectsZ().

                                                                                          {
   std::vector<HepMC::GenParticle*> mothers;
   if ( boson->production_vertex() ) {
     HepMC::GenVertex::particle_iterator mother;
     for (mother = boson->production_vertex()->particles_begin(HepMC::parents); mother!= boson->production_vertex()->particles_end(HepMC::parents); mother++ ) {
       if((*mother)->pdg_id() == boson->pdg_id()) return GetMothers(*mother);
       mothers.push_back(*mother);
     }
   }
   return mothers;
 }

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

private

Definition at line 301 of file TauValidation.cc.

Referenced by analyze(), and findFSRandBrem().

                                                                {
   if ( tau->end_vertex() ) {
     HepMC::GenVertex::particle_iterator dau;
     for (dau = tau->end_vertex()->particles_begin(HepMC::children); dau!= tau->end_vertex()->particles_end(HepMC::children); dau++ ) {
       int dau_pid = (*dau)->pdg_id();
       if(dau_pid == tau->pdg_id()) return false;
     }
   }
   return true;
 }

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

private

Definition at line 618 of file TauValidation.cc.

References leadingPionP4().

Referenced by rtau().

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

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

private

Definition at line 622 of file TauValidation.cc.

References abs, p4, and evf::utils::pid.

Referenced by leadingPionMomentum(), spinEffects(), and spinEffectsZ().

                                                                       {
 
     TLorentzVector p4(0,0,0,0);
 
         if ( tau->end_vertex() ) {
               HepMC::GenVertex::particle_iterator des;
               for(des = tau->end_vertex()->particles_begin(HepMC::descendants);
                   des!= tau->end_vertex()->particles_end(HepMC::descendants);++des ) {
                         int pid = (*des)->pdg_id();
 
                         if(abs(pid) == 15) return leadingPionP4(*des);
 
                         if(!(abs(pid)==211 || abs(pid)==13 || abs(pid)==11)) continue;
  
             if((*des)->momentum().rho() > p4.P()) {
                 p4 = TLorentzVector((*des)->momentum().px(),
                                                     (*des)->momentum().py(),
                                                     (*des)->momentum().pz(),
                                                     (*des)->momentum().e());
             } 
                 }
         }
 
     return p4;
 }

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

private

Definition at line 648 of file TauValidation.cc.

References configurableAnalysis::GenParticle, GetMother(), and m.

Referenced by spinEffects().

                                                                  {
   const HepMC::GenParticle* m=GetMother(tau);
   return TLorentzVector(m->momentum().px(),m->momentum().py(),m->momentum().pz(),m->momentum().e());
 }

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

private

Definition at line 679 of file TauValidation.cc.

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

Referenced by analyze().

                                                                      {
   // Find First tau in chain
   std::vector<const HepMC::GenParticle*> TauList;
   findTauList(tau,TauList);
 
   // Get List of Gauge Boson to tau(s) FSR and Brem
   bool passedW=false;
   std::vector<const HepMC::GenParticle*> ListofFSR;  ListofFSR.clear();
   std::vector<const HepMC::GenParticle*> ListofBrem; ListofBrem.clear();
   std::vector<const HepMC::GenParticle*> FSR_photos; FSR_photos.clear();
   double BosonScale(1);
   if(TauList.size()>0){
     TauValidation::findFSRandBrem(TauList.at(0),passedW,ListofFSR,ListofBrem);
     TauValidation::FindPhotosFSR(TauList.at(0),FSR_photos,BosonScale);
 
     // Add the Tau Brem. information
     TauBremPhotonsN->Fill(ListofBrem.size(),weight);
     double photonPtSum=0;
     for(unsigned int i=0;i<ListofBrem.size();i++){
       photonPtSum+=ListofBrem.at(i)->momentum().perp();
       TauBremPhotonsPt->Fill(ListofBrem.at(i)->momentum().perp(),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)->momentum().perp();
     TauFSRPhotonsPt->Fill(ListofFSR.at(i)->momentum().perp(),weight);
       }
       double FSR_photosSum(0);
       for(unsigned int i=0;i<FSR_photos.size();i++){
     FSR_photosSum+=FSR_photos.at(i)->momentum().perp();
     TauFSRPhotonsPt->Fill(FSR_photos.at(i)->momentum().perp()/BosonScale,weight*BosonScale);
       }
       TauFSRPhotonsPtSum->Fill(photonPtSum+FSR_photosSum/BosonScale,weight);
     }
   }
 }

void TauValidation::rtau	(	const HepMC::GenParticle *	tau,
		int	mother,
		int	decay,
		double	weight
	)

private

Definition at line 470 of file TauValidation.cc.

References abs, MonitorElement::Fill(), leadingPionMomentum(), pi1pi0, tauEtCut, TauRtauHpm, TauRtauW, and visibleTauEnergy().

Referenced by analyze().

                                                                                         {
 
     if(decay != pi1pi0) return; // polarization only for 1-prong hadronic taus with one neutral pion to make a clean case
 
     if(tau->momentum().perp() < tauEtCut) return; // rtau visible only for boosted taus
     
     double rTau = 0;
     double ltrack = leadingPionMomentum(tau, weight);
     double visibleTauE = visibleTauEnergy(tau);
 
     if(visibleTauE != 0) rTau = ltrack/visibleTauE;
 
     if(abs(mother) == 24) TauRtauW->Fill(rTau,weight);
         if(abs(mother) == 37) TauRtauHpm->Fill(rTau,weight); 
 }

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

private

Definition at line 486 of file TauValidation.cc.

References a1, abs, relval_parameters_module::energy, MonitorElement::Fill(), gamma, i, TauDecay::JAK_A1_3PI, TauDecay::JAK_ELECTRON, TauDecay::JAK_MUON, TauDecay::JAK_PION, TauDecay::JAK_RHO_PIPI0, leadingPionP4(), motherP4(), pi, PdtPdgMini::pi0, PdtPdgMini::pi_minus, PdtPdgMini::pi_plus, rho, TauSpinEffectsHpm_eX, TauSpinEffectsHpm_muX, TauSpinEffectsHpm_UpsilonA1, TauSpinEffectsHpm_UpsilonRho, TauSpinEffectsHpm_X, TauSpinEffectsW_eX, TauSpinEffectsW_muX, TauSpinEffectsW_UpsilonA1, TauSpinEffectsW_UpsilonRho, TauSpinEffectsW_X, and histoStyle::weight.

Referenced by analyze().

                                                                                                                                   {
   if(decay == TauDecay::JAK_PION || decay == TauDecay::JAK_MUON || decay == TauDecay::JAK_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::JAK_PION){
       if(abs(mother) == 24) TauSpinEffectsW_X->Fill(energy,weight); 
       if(abs(mother) == 37) TauSpinEffectsHpm_X->Fill(energy,weight);
     }
     if(decay == TauDecay::JAK_MUON){
       if(abs(mother) == 24) TauSpinEffectsW_muX->Fill(energy,weight);
       if(abs(mother) == 37) TauSpinEffectsHpm_muX->Fill(energy,weight);
     }
     if(decay == TauDecay::JAK_ELECTRON){
       if(abs(mother) == 24) TauSpinEffectsW_eX->Fill(energy,weight);
       if(abs(mother) == 37) TauSpinEffectsHpm_eX->Fill(energy,weight);
     }
 
   }
   else if(decay==TauDecay::JAK_RHO_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)->momentum().px(),part.at(i)->momentum().py(),part.at(i)->momentum().pz(),part.at(i)->momentum().e());
       if(abs(part.at(i)->pdg_id())==PdtPdgMini::pi_plus){pi+=LV; rho+=LV;}
       if(abs(part.at(i)->pdg_id())==PdtPdgMini::pi0){rho+=LV;}
     }
     if(abs(mother) == 24) TauSpinEffectsW_UpsilonRho->Fill(2*pi.Pt()/rho.Pt()-1,weight);
     if(abs(mother) == 37) TauSpinEffectsHpm_UpsilonRho->Fill(2*pi.Pt()/rho.Pt()-1,weight);
   }
   else if(decay==TauDecay::JAK_A1_3PI){ // 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)->momentum().px(),part.at(i)->momentum().py(),part.at(i)->momentum().pz(),part.at(i)->momentum().e());
       if(part.at(i)->pdg_id()==PdtPdgMini::pi_plus){ pi_p+=LV; a1+=LV; nplus++;}
       if(part.at(i)->pdg_id()==PdtPdgMini::pi_minus){pi_m+=LV; a1+=LV; nminus++;}
     }
     double gamma=0;
     if(nplus+nminus==3 && nplus==1)  gamma=2*pi_p.Pt()/a1.Pt()-1;
     if(nplus+nminus==3 && nminus==1) gamma=2*pi_m.Pt()/a1.Pt()-1;
     else{
       pi_p+=pi_m; gamma=2*pi_p.Pt()/a1.Pt()-1;
     }
     if(abs(mother) == 24) TauSpinEffectsW_UpsilonA1->Fill(gamma,weight);
     if(abs(mother) == 37) TauSpinEffectsHpm_UpsilonA1->Fill(gamma,weight);
   }
 }

void TauValidation::spinEffectsZ	(	const HepMC::GenParticle *	boson,
		double	weight
	)

private

Definition at line 535 of file TauValidation.cc.

References abs, PdtPdgMini::anti_d, PdtPdgMini::anti_u, DeDxDiscriminatorTools::charge(), PdtPdgMini::d, electron, MonitorElement::Fill(), findMother(), fPDGTable, GetMothers(), PdtPdgMini::Higgs0, i, leadingPionP4(), m, muon, pi, evf::utils::pid, lumiQueryAPI::q, tauDecayChannel(), TauSpinEffectsH_eX, TauSpinEffectsH_muX, TauSpinEffectsH_MVis, TauSpinEffectsH_Xb, TauSpinEffectsH_Xf, TauSpinEffectsH_Zs, TauSpinEffectsZ_eX, TauSpinEffectsZ_muX, TauSpinEffectsZ_MVis, TauSpinEffectsZ_Xb, TauSpinEffectsZ_Xf, TauSpinEffectsZ_Zs, PdtPdgMini::u, histoStyle::weight, Z, PdtPdgMini::Z0, zsbins, zsmax, zsmin, and Zstoa().

Referenced by analyze().

                                                                             {
 
   TLorentzVector tautau(0,0,0,0);
   TLorentzVector pipi(0,0,0,0);
   TLorentzVector taum(0,0,0,0);
   int nSinglePionDecays(0),nSingleMuonDecays(0),nSingleElectronDecays(0);
   double x1(0),x2(0); 
   TLorentzVector Zboson(boson->momentum().px(),boson->momentum().py(),boson->momentum().pz(),boson->momentum().e());
   if ( boson->end_vertex() ) {
     HepMC::GenVertex::particle_iterator des;
     for(des = boson->end_vertex()->particles_begin(HepMC::children);    des!= boson->end_vertex()->particles_end(HepMC::children);++des ) {
       int pid = (*des)->pdg_id();
       if(abs(findMother(*des)) != 15 &&  abs(pid) == 15 && (tauDecayChannel(*des) == pi || tauDecayChannel(*des) == muon || tauDecayChannel(*des) == electron )){
     if(tauDecayChannel(*des) == pi)nSinglePionDecays++;
     if(tauDecayChannel(*des) == muon)nSingleMuonDecays++;
     if(tauDecayChannel(*des) == electron)nSingleElectronDecays++;
     TLorentzVector LVtau((*des)->momentum().px(),(*des)->momentum().py(),(*des)->momentum().pz(),(*des)->momentum().e());
     tautau += LVtau;
     TLorentzVector LVpi=leadingPionP4(*des);
     pipi+=LVpi;
     const HepPDT::ParticleData*  pd = fPDGTable->particle((*des)->pdg_id ());
     int charge = (int) pd->charge();
     LVtau.Boost(-1*Zboson.BoostVector());
     LVpi.Boost(-1*Zboson.BoostVector());
     if(charge<0){x1=LVpi.P()/LVtau.E(); taum=LVtau;}
     else{ x2=LVpi.P()/LVtau.E();}
      }
     }
   }
   if(nSingleMuonDecays==2){
     if(abs(boson->pdg_id())==PdtPdgMini::Z0)     TauSpinEffectsZ_muX->Fill(x1,weight);
     if(abs(boson->pdg_id())==PdtPdgMini::Higgs0) TauSpinEffectsH_muX->Fill(x1,weight);
   }
   if(nSingleElectronDecays==2){
     if(abs(boson->pdg_id())==PdtPdgMini::Z0)     TauSpinEffectsZ_eX->Fill(x1,weight);
     if(abs(boson->pdg_id())==PdtPdgMini::Higgs0) TauSpinEffectsH_eX->Fill(x1,weight);
   }
   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->pdg_id())==PdtPdgMini::Z0)     TauSpinEffectsZ_Zs->Fill(zs,weight);
     if(abs(boson->pdg_id())==PdtPdgMini::Higgs0) TauSpinEffectsH_Zs->Fill(zs,weight);
     break;
       }
     }
     if(abs(boson->pdg_id())==PdtPdgMini::Z0)     TauSpinEffectsZ_MVis->Fill(pipi.M()/tautau.M(),weight);
     if(abs(boson->pdg_id())==PdtPdgMini::Higgs0) TauSpinEffectsH_MVis->Fill(pipi.M()/tautau.M(),weight);
  
     if(x1!=0){
       const std::vector<HepMC::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)->pdg_id()==PdtPdgMini::d      || m.at(i)->pdg_id()==PdtPdgMini::u      ){q++;}
     if(m.at(i)->pdg_id()==PdtPdgMini::anti_d || m.at(i)->pdg_id()==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->pdg_id())==PdtPdgMini::Z0)      TauSpinEffectsZ_Xf->Fill(x1,weight);
       if(abs(boson->pdg_id())==PdtPdgMini::Higgs0) TauSpinEffectsH_Xf->Fill(x1,weight);
     }
     else{
       if(abs(boson->pdg_id())==PdtPdgMini::Z0)      TauSpinEffectsZ_Xb->Fill(x1,weight);
       if(abs(boson->pdg_id())==PdtPdgMini::Higgs0) TauSpinEffectsH_Xb->Fill(x1,weight);
     }
       }
     }
   }
 }

int TauValidation::tauDecayChannel	(	const HepMC::GenParticle *	tau,
		double	weight = `0.0`
	)

private

Definition at line 417 of file TauValidation.cc.

References a1, abs, electron, MonitorElement::Fill(), K, Kstar, muon, pi, pi1pi0, evf::utils::pid, pinpi0, rho, stable, TauDecayChannels, tripi, tripinpi0, and undetermined.

Referenced by analyze(), and spinEffectsZ().

                                                                             {
 
   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;
   
   if ( tau->end_vertex() ) {
     HepMC::GenVertex::particle_iterator des;
     for(des = tau->end_vertex()->particles_begin(HepMC::descendants);
     des!= tau->end_vertex()->particles_end(HepMC::descendants);++des ) {
       int pid = (*des)->pdg_id();
       if(abs(pid) == 15) return tauDecayChannel(*des,weight);
       
       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++;
       
     }
   }
   // 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 HepMC::GenParticle *	tau,
		double	weight
	)

private

Definition at line 372 of file TauValidation.cc.

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

Referenced by analyze().

                                                                       {
 
   if(abs(tau->pdg_id()) != 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 HepMC::GenParticle *	tau,
		double	weight
	)

private

Definition at line 397 of file TauValidation.cc.

References abs, DeDxDiscriminatorTools::charge(), MonitorElement::Fill(), fPDGTable, evf::utils::pid, and TauProngs.

Referenced by analyze().

                                                                       {
   int nProngs = 0;
   if ( tau->end_vertex() ) {
     HepMC::GenVertex::particle_iterator des;
     for(des = tau->end_vertex()->particles_begin(HepMC::descendants);
     des!= tau->end_vertex()->particles_end(HepMC::descendants);++des ) {
       int pid = (*des)->pdg_id();
       if(abs(pid) == 15) return tauProngs(*des, weight);
       if((*des)->status() != 1) continue; // dont count unstable particles
       
       const HepPDT::ParticleData*  pd = fPDGTable->particle((*des)->pdg_id ());
       int charge = (int) pd->charge();
       if(charge == 0) continue;
       nProngs++;
     }
   }
   TauProngs->Fill(nProngs,weight);
   return nProngs;
 }

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

private

Definition at line 653 of file TauValidation.cc.

References abs, p4, and evf::utils::pid.

Referenced by rtau().

                                                                  {
     TLorentzVector p4(tau->momentum().px(),
                           tau->momentum().py(),
                           tau->momentum().pz(),
                           tau->momentum().e());
 
         if ( tau->end_vertex() ) {
               HepMC::GenVertex::particle_iterator des;
               for(des = tau->end_vertex()->particles_begin(HepMC::descendants);
                   des!= tau->end_vertex()->particles_end(HepMC::descendants);++des ) {
                         int pid = (*des)->pdg_id();
 
                         if(abs(pid) == 15) return visibleTauEnergy(*des);
 
                         if(abs(pid) == 12 || abs(pid) == 14 || abs(pid) == 16) {
                 p4 -= TLorentzVector((*des)->momentum().px(),
                                          (*des)->momentum().py(),
                                          (*des)->momentum().pz(),
                                          (*des)->momentum().e());
             }
                 }
         }
 
     return p4.E();
 }

double TauValidation::Zstoa ( double zs )

private

Definition at line 609 of file TauValidation.cc.

References a, and mathSSE::sqrt().

Referenced by spinEffectsZ().

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

Member Data Documentation

WeightManager TauValidation::_wmanager

private

Definition at line 96 of file TauValidation.h.

Referenced by analyze().

DQMStore* TauValidation::dbe

private

ME's "container".

Definition at line 106 of file TauValidation.h.

Referenced by beginJob(), and TauValidation().

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

private

PDT table.

Definition at line 103 of file TauValidation.h.

Referenced by beginRun(), spinEffectsZ(), and tauProngs().

edm::InputTag TauValidation::hepmcCollection_

private

Definition at line 98 of file TauValidation.h.

Referenced by analyze().

MonitorElement* TauValidation::JAKID

private

Definition at line 119 of file TauValidation.h.

Referenced by analyze(), and beginJob().

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

private

Definition at line 120 of file TauValidation.h.

Referenced by analyze(), and beginJob().

unsigned int TauValidation::NJAKID

private

Definition at line 118 of file TauValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * TauValidation::nPrimeTaus

private

Definition at line 108 of file TauValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement* TauValidation::nTaus

private

Definition at line 108 of file TauValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * TauValidation::TauBremPhotonsN

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and photons().

MonitorElement * TauValidation::TauBremPhotonsPt

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and photons().

MonitorElement * TauValidation::TauBremPhotonsPtSum

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and photons().

MonitorElement * TauValidation::TauDecayChannels

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and tauDecayChannel().

MonitorElement * TauValidation::TauEta

private

Definition at line 109 of file TauValidation.h.

Referenced by analyze(), and beginJob().

double TauValidation::tauEtCut

private

Definition at line 100 of file TauValidation.h.

Referenced by rtau().

MonitorElement * TauValidation::TauFSRPhotonsN

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and photons().

MonitorElement * TauValidation::TauFSRPhotonsPt

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and photons().

MonitorElement * TauValidation::TauFSRPhotonsPtSum

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and photons().

MonitorElement * TauValidation::TauMothers

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and tauMother().

MonitorElement * TauValidation::TauPhi

private

Definition at line 109 of file TauValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * TauValidation::TauProngs

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and tauProngs().

MonitorElement* TauValidation::TauPt

private

Definition at line 109 of file TauValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * TauValidation::TauRtauHpm

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and rtau().

MonitorElement * TauValidation::TauRtauW

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and rtau().

MonitorElement * TauValidation::TauSpinEffectsH_eX

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZ().

MonitorElement * TauValidation::TauSpinEffectsH_muX

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZ().

MonitorElement * TauValidation::TauSpinEffectsH_MVis

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZ().

MonitorElement * TauValidation::TauSpinEffectsH_Xb

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZ().

MonitorElement * TauValidation::TauSpinEffectsH_Xf

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZ().

MonitorElement * TauValidation::TauSpinEffectsH_Zs

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZ().

MonitorElement * TauValidation::TauSpinEffectsHpm_eX

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffects().

MonitorElement * TauValidation::TauSpinEffectsHpm_muX

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffects().

MonitorElement * TauValidation::TauSpinEffectsHpm_UpsilonA1

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffects().

MonitorElement * TauValidation::TauSpinEffectsHpm_UpsilonRho

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffects().

MonitorElement * TauValidation::TauSpinEffectsHpm_X

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffects().

MonitorElement * TauValidation::TauSpinEffectsW_eX

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffects().

MonitorElement * TauValidation::TauSpinEffectsW_muX

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffects().

MonitorElement * TauValidation::TauSpinEffectsW_UpsilonA1

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffects().

MonitorElement * TauValidation::TauSpinEffectsW_UpsilonRho

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffects().

MonitorElement * TauValidation::TauSpinEffectsW_X

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffects().

MonitorElement * TauValidation::TauSpinEffectsZ_eX

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZ().

MonitorElement * TauValidation::TauSpinEffectsZ_muX

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZ().

MonitorElement * TauValidation::TauSpinEffectsZ_MVis

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZ().

MonitorElement * TauValidation::TauSpinEffectsZ_Xb

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZ().

MonitorElement * TauValidation::TauSpinEffectsZ_Xf

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZ().

MonitorElement * TauValidation::TauSpinEffectsZ_Zs

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZ().

int TauValidation::zsbins

private

Definition at line 122 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZ().

double TauValidation::zsmax

private

Definition at line 123 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZ().

double TauValidation::zsmin

private

Definition at line 123 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZ().

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