#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

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
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

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

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) 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	spinEffectsWHpm (const HepMC::GenParticle , int, int, std::vector< HepMC::GenParticle > &part, double weight)

void	spinEffectsZH (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
DQMStore *	dbe
	ME's "container". More...

MonitorElement *	DecayLength

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

edm::InputTag	hepmcCollection_

edm::EDGetTokenT < edm::HepMCProduct >	hepmcCollectionToken_

MonitorElement *	JAKID

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

MonitorElement *	LifeTime

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_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_Xb

MonitorElement *	TauSpinEffectsZ_Xf

MonitorElement *	TauSpinEffectsZ_Zs

WeightManager	wmanager_

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::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 34 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 38 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 52 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 19 of file TauValidation.cc.

References dbe, hepmcCollection_, hepmcCollectionToken_, and cppFunctionSkipper::operator.

                                                         : 
   wmanager_(iPSet,consumesCollector())
   ,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->();
 
   hepmcCollectionToken_=consumes<HepMCProduct>(hepmcCollection_);
 }

TauValidation::~TauValidation ( )

virtual

Definition at line 34 of file TauValidation.cc.

34 {}

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

References funct::abs(), TauDecay_CMSSW::AnalyzeTau(), beta, trackerHits::c, DecayLength, MonitorElement::Fill(), TauDecay_CMSSW::Get_TauDecayProducts(), edm::Event::getByToken(), hepmcCollectionToken_, PdtPdgMini::Higgs0, i, isLastTauinChain(), TauDecay::isTauFinalStateParticle(), TauDecay::JAK_A1_3PI, TauDecay::JAK_KPIK, TauDecay::JAK_KPIPI, JAKID, JAKInvMass, PdtPdgMini::K_plus, LifeTime, NJAKID, nPrimeTaus, TauDecay::nProng(), nTaus, PdtPdgMini::nu_e, PdtPdgMini::nu_mu, PdtPdgMini::nu_tau, photons(), PdtPdgMini::pi_plus, rtau(), spinEffectsWHpm(), spinEffectsZH(), tauDecayChannel(), TauEta, tauMother(), TauPhi, tauProngs(), TauPt, WeightManager::weight(), histoStyle::weight, wmanager_, and PdtPdgMini::Z0.

 { 
   edm::Handle<HepMCProduct> evt;
   iEvent.getByToken(hepmcCollectionToken_, 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())==PdtPdgMini::Z0 || abs((*iter)->pdg_id())==PdtPdgMini::Higgs0){
        spinEffectsZH(*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();
          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)->pdg_id())!=PdtPdgMini::nu_tau && TD.isTauFinalStateParticle(part.at(i)->pdg_id()) && !hasDL){
          PV=TVector3((*iter)->production_vertex()->point3d().x(),(*iter)->production_vertex()->point3d().y(),(*iter)->production_vertex()->point3d().z());
          SV=TVector3(part.at(i)->production_vertex()->point3d().x(),part.at(i)->production_vertex()->point3d().y(),part.at(i)->production_vertex()->point3d().z());
          TVector3 DL=SV-PV;
          DecayLength->Fill(DL.Mag(),weight);
          double c(2.99792458E8),Ltau(DL.Mag()/1000),beta((*iter)->momentum().rho()/(*iter)->momentum().m());
          LifeTime->Fill( Ltau/(c*beta),weight);
            }
 
           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 36 of file TauValidation.cc.

References A0, a1, DQMStore::book1D(), D, dbe, DecayLength, electron, gamma, H0, Hpm, HSM, i, TauDecay::JAK_A1_3PI, TauDecay::JAK_KPIK, TauDecay::JAK_KPIPI, JAKID, JAKInvMass, K, Kstar, LifeTime, muon, NJAKID, nPrimeTaus, nTaus, other, Pi, 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_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_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}");
 
     DecayLength = dbe->book1D("DecayLength","#tau Decay Length", 100 ,0,20);  DecayLength->setAxisTitle("L_{#tau} (mm)");
     LifeTime =  dbe->book1D("LifeTime","#tau LifeTime ", 100 ,0,1000E-15);     LifeTime->setAxisTitle("#tau_{#tau}s)");
 
     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");
 
     TauSpinEffectsH_pipiAcoplanarity = dbe->book1D("TauSpinEffectsH_pipiAcoplanarity","H Acoplanarity for #pi^{-}#pi^{+}", 50 ,0,2*TMath::Pi()); TauSpinEffectsH_pipiAcoplanarity->setAxisTitle("Acoplanarity"); 
 
     TauSpinEffectsH_pipiAcollinearity = dbe->book1D("TauSpinEffectsH_pipiAcollinearity","H Acollinearity for #pi^{-}#pi^{+}", 50 ,0,TMath::Pi()); TauSpinEffectsH_pipiAcollinearity->setAxisTitle("Acollinearity");
     TauSpinEffectsH_pipiAcollinearityzoom = dbe->book1D("TauSpinEffectsH_pipiAcollinearityzoom","H Acollinearity for #pi^{-}#pi^{+}", 50 ,3,TMath::Pi()); TauSpinEffectsH_pipiAcollinearityzoom->setAxisTitle("Acollinearity");
 
     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_X= dbe->book1D("TauSpinEffectsZX","X of #tau^{-}", 25 ,0,1.0);           TauSpinEffectsZ_X->setAxisTitle("X");
     TauSpinEffectsH_X= dbe->book1D("TauSpinEffectsH_X","X of #tau^{-}", 25 ,0,1.0);           TauSpinEffectsH_X->setAxisTitle("X");
 
     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");
 
      TauSpinEffectsH_rhorhoAcoplanarityminus = dbe->book1D("TauSpinEffectsH_rhorhoAcoplanarityminus","#phi^{*-} (acoplanarity) for Higgs #rightarrow #rho-#rho (y_{1}*y_{2}<0)", 32 ,0,2*TMath::Pi());     TauSpinEffectsH_rhorhoAcoplanarityminus->setAxisTitle("#phi^{*-} (Acoplanarity)");
      TauSpinEffectsH_rhorhoAcoplanarityplus = dbe->book1D("TauSpinEffectsH_rhorhoAcoplanarityplus","#phi^{*+} (acoplanarity) for Higgs #rightarrow #rho-#rho (y_{1}*y_{2}>0)", 32 ,0,2*TMath::Pi());     TauSpinEffectsH_rhorhoAcoplanarityplus->setAxisTitle("#phi^{*+} (Acoplanarity)");
 
     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 165 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 161 of file TauValidation.cc.

                           {
   return;
 }

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

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 171 of file TauValidation.cc.

171 {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 350 of file TauValidation.cc.

References funct::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 313 of file TauValidation.cc.

References parents.

Referenced by spinEffectsZH(), 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 374 of file TauValidation.cc.

References funct::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 338 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 288 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 300 of file TauValidation.cc.

References parents.

Referenced by spinEffectsZH().

                                                                                          {
   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 326 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 767 of file TauValidation.cc.

References leadingPionP4().

Referenced by rtau().

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

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

private

Definition at line 771 of file TauValidation.cc.

References funct::abs(), p4, and sysUtil::pid.

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

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

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

Referenced by spinEffectsWHpm().

                                                                  {
   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 828 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 495 of file TauValidation.cc.

References funct::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::spinEffectsWHpm	(	const HepMC::GenParticle *	tau,
		int	mother,
		int	decay,
		std::vector< HepMC::GenParticle * > &	part,
		double	weight
	)

private

Definition at line 511 of file TauValidation.cc.

References a1, funct::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::spinEffectsZH	(	const HepMC::GenParticle *	boson,
		double	weight
	)

private

Definition at line 560 of file TauValidation.cc.

References funct::abs(), PdtPdgMini::anti_d, PdtPdgMini::anti_u, DeDxDiscriminatorTools::charge(), PdtPdgMini::d, electron, MonitorElement::Fill(), findMother(), fPDGTable, GetMothers(), PdtPdgMini::Higgs0, i, leadingPionP4(), m, muon, n, Pi, pi, pi1pi0, sysUtil::pid, lumiQueryAPI::q, rho, tauDecayChannel(), 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_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);
    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->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(tauDecayChannel(*des) == pi){
       if(abs(boson->pdg_id())==PdtPdgMini::Z0) TauSpinEffectsZ_X->Fill(LVpi.P()/LVtau.E(),weight);
       if(abs(boson->pdg_id())==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();}
      }
        if(abs(findMother(*des)) != 15 &&  abs(pid) == 15 && (tauDecayChannel(*des) == rho || tauDecayChannel(*des) == pi1pi0) ){
      if ( (*des)->end_vertex() ) {
        HepMC::GenVertex::particle_iterator tauprod;
        TLorentzVector LVtau((*des)->momentum().px(),(*des)->momentum().py(),(*des)->momentum().pz(),(*des)->momentum().e());
        if(pid == 15)taum=LVtau; 
        if(pid ==-15)taup=LVtau;
        for(tauprod = (*des)->end_vertex()->particles_begin(HepMC::descendants); tauprod!= (*des)->end_vertex()->particles_end(HepMC::descendants);++tauprod ) {
          int pid_d = (*tauprod)->pdg_id();
          if(abs(pid_d)==211 || abs(pid_d)==111){
            TLorentzVector LV((*tauprod)->momentum().px(),(*tauprod)->momentum().py(),(*tauprod)->momentum().pz(),(*tauprod)->momentum().e());
            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(abs(findMother(*des)) != 15 &&  abs(pid) == 15 && (tauDecayChannel(*des) == pi) ){
          if ( (*des)->end_vertex() ) {
        HepMC::GenVertex::particle_iterator tauprod;
            TLorentzVector LVtau((*des)->momentum().px(),(*des)->momentum().py(),(*des)->momentum().pz(),(*des)->momentum().e());
            if(pid == 15)taum=LVtau;
            if(pid ==-15)taup=LVtau;
            for(tauprod = (*des)->end_vertex()->particles_begin(HepMC::descendants); tauprod!= (*des)->end_vertex()->particles_end(HepMC::descendants);++tauprod ) {
              int pid_d = (*tauprod)->pdg_id();
              if(abs(pid_d)==211 || abs(pid_d)==111){
                TLorentzVector LV((*tauprod)->momentum().px(),(*tauprod)->momentum().py(),(*tauprod)->momentum().pz(),(*tauprod)->momentum().e());
                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_rhominus.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->pdg_id())==PdtPdgMini::Higgs0 && y1*y2<0) TauSpinEffectsH_rhorhoAcoplanarityminus->Fill(Acoplanarity,weight);
      if(abs(boson->pdg_id())==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->pdg_id())==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()*taup.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->pdg_id())==PdtPdgMini::Higgs0) TauSpinEffectsH_pipiAcoplanarity->Fill(Acoplanarity,weight);
      taup=tauporig;
      taum=taumorig;
 
    }
 
 
   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 442 of file TauValidation.cc.

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

Referenced by analyze(), and spinEffectsZH().

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

References A0, funct::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 422 of file TauValidation.cc.

References funct::abs(), DeDxDiscriminatorTools::charge(), MonitorElement::Fill(), fPDGTable, sysUtil::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 802 of file TauValidation.cc.

References funct::abs(), p4, and sysUtil::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 758 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

DQMStore* TauValidation::dbe

private

ME's "container".

Definition at line 104 of file TauValidation.h.

Referenced by beginJob(), and TauValidation().

MonitorElement * TauValidation::DecayLength

private

Definition at line 107 of file TauValidation.h.

Referenced by analyze(), and beginJob().

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

private

PDT table.

Definition at line 101 of file TauValidation.h.

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

edm::InputTag TauValidation::hepmcCollection_

private

Definition at line 96 of file TauValidation.h.

Referenced by TauValidation().

edm::EDGetTokenT<edm::HepMCProduct> TauValidation::hepmcCollectionToken_

private

Definition at line 126 of file TauValidation.h.

Referenced by analyze(), and TauValidation().

MonitorElement* TauValidation::JAKID

private

Definition at line 120 of file TauValidation.h.

Referenced by analyze(), and beginJob().

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

private

Definition at line 121 of file TauValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * TauValidation::LifeTime

private

Definition at line 107 of file TauValidation.h.

Referenced by analyze(), and beginJob().

unsigned int TauValidation::NJAKID

private

Definition at line 119 of file TauValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * TauValidation::nPrimeTaus

private

Definition at line 106 of file TauValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement* TauValidation::nTaus

private

Definition at line 106 of file TauValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * TauValidation::TauBremPhotonsN

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and photons().

MonitorElement * TauValidation::TauBremPhotonsPt

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and photons().

MonitorElement * TauValidation::TauBremPhotonsPtSum

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and photons().

MonitorElement * TauValidation::TauDecayChannels

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and tauDecayChannel().

MonitorElement * TauValidation::TauEta

private

Definition at line 107 of file TauValidation.h.

Referenced by analyze(), and beginJob().

double TauValidation::tauEtCut

private

Definition at line 98 of file TauValidation.h.

Referenced by rtau().

MonitorElement * TauValidation::TauFSRPhotonsN

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and photons().

MonitorElement * TauValidation::TauFSRPhotonsPt

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and photons().

MonitorElement * TauValidation::TauFSRPhotonsPtSum

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and photons().

MonitorElement * TauValidation::TauMothers

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and tauMother().

MonitorElement * TauValidation::TauPhi

private

Definition at line 107 of file TauValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * TauValidation::TauProngs

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and tauProngs().

MonitorElement* TauValidation::TauPt

private

Definition at line 107 of file TauValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * TauValidation::TauRtauHpm

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and rtau().

MonitorElement * TauValidation::TauRtauW

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and rtau().

MonitorElement * TauValidation::TauSpinEffectsH_eX

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_muX

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_MVis

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_pipiAcollinearity

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_pipiAcollinearityzoom

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_pipiAcoplanarity

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_rhorhoAcoplanarityminus

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_rhorhoAcoplanarityplus

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_X

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_Xb

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_Xf

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_Zs

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsHpm_eX

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsHpm_muX

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsHpm_UpsilonA1

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsHpm_UpsilonRho

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsHpm_X

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsW_eX

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsW_muX

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsW_UpsilonA1

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsW_UpsilonRho

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsW_X

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsZ_eX

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_muX

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_MVis

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_X

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_Xb

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_Xf

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_Zs

private

Definition at line 107 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

WeightManager TauValidation::wmanager_

private

Definition at line 73 of file TauValidation.h.

Referenced by analyze().

int TauValidation::zsbins

private

Definition at line 123 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

double TauValidation::zsmax

private

Definition at line 124 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

double TauValidation::zsmin

private

Definition at line 124 of file TauValidation.h.

Referenced by beginJob(), 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