TauValidation Class Reference

#include <TauValidation.h>

Inheritance diagram for TauValidation:

Public Types
enum	{   undetermined, electron, muon, pi,   rho, a1, K, Kstar,   pi1pi0, pinpi0, tripi, tripinpi0,   stable }
enum	{   other, B, D, gamma,   Z, W, HSM, H0,   A0, Hpm }
Public Types inherited from edm::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 ()

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
WeightManager	_wmanager
DQMStore *	dbe
	ME's "container". More...
MonitorElement *	DecayLength
edm::ESHandle < HepPDT::ParticleDataTable >	fPDGTable
	PDT table. More...
edm::InputTag	hepmcCollection_
MonitorElement *	JAKID
std::vector< std::vector < MonitorElement * > >	JAKInvMass
MonitorElement *	LifeTime
std::vector< int >	n_taupinu
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 *	TauSpinEffectstautau_polvxM
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
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
CurrentProcessingContext const *	currentContext () const

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 19 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)
  ,n_taupinu(200,0)
{    
  dbe = 0;
  dbe = edm::Service<DQMStore>().operator->();
}

TauValidation::~TauValidation ( )

virtual

Definition at line 33 of file TauValidation.cc.

{}

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

References _wmanager, abs, TauDecay_CMSSW::AnalyzeTau(), beta, trackerHits::c, DecayLength, MonitorElement::Fill(), TauDecay_CMSSW::Get_TauDecayProducts(), edm::Event::getByLabel(), hepmcCollection_, 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(), CommonMethods::weight(), and PdtPdgMini::Z0.

{ 
  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())==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 35 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, TauSpinEffectstautau_polvxM, 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)");



     TauSpinEffectstautau_polvxM = dbe->book1D("TauSpinEffectsZ_polvxM","#tau Polarization ", 40,0,200);  TauSpinEffectstautau_polvxM->setAxisTitle("Mass (GeV)");


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

                          {
  return;
}

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

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 175 of file TauValidation.cc.

{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 354 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 317 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 378 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 342 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 292 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 304 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 330 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 793 of file TauValidation.cc.

References leadingPionP4().

Referenced by rtau().

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

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

private

Definition at line 797 of file TauValidation.cc.

References abs, p4, and evf::utils::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 823 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 854 of file TauValidation.cc.

References MonitorElement::Fill(), findFSRandBrem(), FindPhotosFSR(), findTauList(), i, TauBremPhotonsN, TauBremPhotonsPt, TauBremPhotonsPtSum, TauFSRPhotonsN, TauFSRPhotonsPt, TauFSRPhotonsPtSum, and CommonMethods::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 499 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::spinEffectsWHpm ( const HepMC::GenParticle *  tau, int  mother, int  decay, std::vector< HepMC::GenParticle * > &  part, double  weight  )

private

Definition at line 515 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 CommonMethods::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 564 of file TauValidation.cc.

References abs, PdtPdgMini::anti_d, PdtPdgMini::anti_u, DeDxDiscriminatorTools::charge(), PdtPdgMini::d, electron, MonitorElement::Fill(), findMother(), fPDGTable, GetMothers(), MonitorElement::getTH1(), PdtPdgMini::Higgs0, interpolateCardsSimple::histo, i, leadingPionP4(), m, timingPdfMaker::mean, muon, n, n_taupinu, Pi, pi, pi1pi0, evf::utils::pid, lumiQueryAPI::q, rho, mathSSE::sqrt(), 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, TauSpinEffectstautau_polvxM, TauSpinEffectsZ_eX, TauSpinEffectsZ_muX, TauSpinEffectsZ_MVis, TauSpinEffectsZ_X, TauSpinEffectsZ_Xb, TauSpinEffectsZ_Xf, TauSpinEffectsZ_Zs, PdtPdgMini::u, CommonMethods::weight(), vdt::x, 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(tautau.M()>2.0){
        // for both z and h compute tau polarization
        TH1* histo = TauSpinEffectstautau_polvxM->getTH1();
        int n=histo->FindBin(tautau.M());
        double mean=histo->GetBinContent(n);
        double err=histo->GetBinError(n);
        double x=(6.0*(LVpi.P()/LVtau.E()-0.5));
        double M2=err*err*(n_taupinu[n]-1)*n_taupinu[n];
        //use Knuth algorithm: Donald E. Knuth (1998). The Art of Computer Programming, 
        //volume 2: Seminumerical Algorithms, 3rd edn., p. 232. Boston: Addison-Wesley. 
        n_taupinu[n]+=1.0;
        double d=x-mean; 
        mean+=d/n_taupinu[n];
        //std::cout << n << M2 << " " << d << " " << (x-mean) << std::endl;
        M2+=d*(x-mean);
        if(n_taupinu[n]>1){err=sqrt((M2/(n_taupinu[n]-1))/n_taupinu[n]);}
        else{ err=sqrt(M2);}
        //save values
        //std::cout << n << " " << tautau.M() << " " << mean << "+/-" << err << std::endl; 
        histo->SetBinContent(n,mean);
        histo->SetBinError(n,err);
      }
    }
    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 446 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 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 401 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 426 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 828 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 784 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

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

MonitorElement * TauValidation::DecayLength

private

Definition at line 109 of file TauValidation.h.

Referenced by analyze(), and beginJob().

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

private

PDT table.

Definition at line 103 of file TauValidation.h.

Referenced by beginRun(), spinEffectsZH(), 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 122 of file TauValidation.h.

Referenced by analyze(), and beginJob().

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

private

Definition at line 123 of file TauValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * TauValidation::LifeTime

private

Definition at line 109 of file TauValidation.h.

Referenced by analyze(), and beginJob().

std::vector<int> TauValidation::n_taupinu

private

Definition at line 127 of file TauValidation.h.

Referenced by spinEffectsZH().

unsigned int TauValidation::NJAKID

private

Definition at line 121 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 spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_muX

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_MVis

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_pipiAcollinearity

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_pipiAcollinearityzoom

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_pipiAcoplanarity

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_rhorhoAcoplanarityminus

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_rhorhoAcoplanarityplus

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_X

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_Xb

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_Xf

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsH_Zs

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsHpm_eX

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsHpm_muX

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsHpm_UpsilonA1

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsHpm_UpsilonRho

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsHpm_X

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectstautau_polvxM

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsW_eX

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsW_muX

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsW_UpsilonA1

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsW_UpsilonRho

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsW_X

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsWHpm().

MonitorElement * TauValidation::TauSpinEffectsZ_eX

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_muX

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_MVis

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_X

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_Xb

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_Xf

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

MonitorElement * TauValidation::TauSpinEffectsZ_Zs

private

Definition at line 109 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

int TauValidation::zsbins

private

Definition at line 125 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

double TauValidation::zsmax

private

Definition at line 126 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().

double TauValidation::zsmin

private

Definition at line 126 of file TauValidation.h.

Referenced by beginJob(), and spinEffectsZH().