|
|
|
#include <TauValidation.h>
Definition at line 36 of file TauValidation.h.
| anonymous enum |
| anonymous enum |
| 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) { dbe = 0; dbe = edm::Service<DQMStore>().operator->(); }
| TauValidation::~TauValidation | ( | ) | [virtual] |
Definition at line 32 of file TauValidation.cc.
{}
| void TauValidation::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [virtual] |
Gathering the HepMCProduct information
Implements edm::EDAnalyzer.
Definition at line 160 of file TauValidation.cc.
References _wmanager, abs, TauDecay_CMSSW::AnalyzeTau(), MonitorElement::Fill(), TauDecay_CMSSW::Get_TauDecayProducts(), edm::Event::getByLabel(), hepmcCollection_, i, isLastTauinChain(), TauDecay::isTauFinalStateParticle(), TauDecay::JAK_A1_3PI, TauDecay::JAK_KPIK, TauDecay::JAK_KPIPI, JAKID, JAKInvMass, PdtPdgMini::K_plus, NJAKID, nPrimeTaus, TauDecay::nProng(), nTaus, PdtPdgMini::nu_e, PdtPdgMini::nu_mu, PdtPdgMini::nu_tau, photons(), PdtPdgMini::pi_plus, rtau(), spinEffects(), spinEffectsZ(), tauDecayChannel(), TauEta, tauMother(), TauPhi, tauProngs(), TauPt, CommonMethods::weight(), and WeightManager::weight().
{
edm::Handle<HepMCProduct> evt;
iEvent.getByLabel(hepmcCollection_, evt);
//Get EVENT
HepMC::GenEvent *myGenEvent = new HepMC::GenEvent(*(evt->GetEvent()));
double weight = _wmanager.weight(iEvent);
/*
edm::Handle<double> WT;
iEvent.getByLabel(edm::InputTag("TauSpinnerGen","TauSpinnerWT"),WT);
weight = 1.0;
if(*(WT.product())>1e-3 && *(WT.product())<=10.0) weight=(*(WT.product()));
else {weight=1.0;}
*/
// find taus
for(HepMC::GenEvent::particle_const_iterator iter = myGenEvent->particles_begin(); iter != myGenEvent->particles_end(); iter++) {
if(abs((*iter)->pdg_id())==23){
spinEffectsZ(*iter,weight);
}
if(abs((*iter)->pdg_id())==15){
if(isLastTauinChain(*iter)){
nTaus->Fill(0.5,weight);
int mother = tauMother(*iter,weight);
int decaychannel = tauDecayChannel(*iter,weight);
tauProngs(*iter, weight);
if(mother>-1){ // exclude B, D and other non-signal decay modes
nPrimeTaus->Fill(0.5,weight);
TauPt->Fill((*iter)->momentum().perp(),weight);
TauEta->Fill((*iter)->momentum().eta(),weight);
TauPhi->Fill((*iter)->momentum().phi(),weight);
rtau(*iter,mother,decaychannel,weight);
photons(*iter,weight);
}
//Adding JAKID and Mass information
//
TauDecay_CMSSW TD;
unsigned int jak_id, TauBitMask;
if(TD.AnalyzeTau((*iter),jak_id,TauBitMask,false,false)){
JAKID->Fill(jak_id,weight);
if(jak_id<=NJAKID){
int tcharge=(*iter)->pdg_id()/abs((*iter)->pdg_id());
std::vector<HepMC::GenParticle*> part=TD.Get_TauDecayProducts();
spinEffects(*iter,mother,jak_id,part,weight);
TLorentzVector LVQ(0,0,0,0);
TLorentzVector LVS12(0,0,0,0);
TLorentzVector LVS13(0,0,0,0);
TLorentzVector LVS23(0,0,0,0);
bool haspart1=false;
for(unsigned int i=0;i<part.size();i++){
if(TD.isTauFinalStateParticle(part.at(i)->pdg_id()) &&
abs(part.at(i)->pdg_id())!=PdtPdgMini::nu_e &&
abs(part.at(i)->pdg_id())!=PdtPdgMini::nu_mu &&
abs(part.at(i)->pdg_id())!=PdtPdgMini::nu_tau ){
TLorentzVector LV(part.at(i)->momentum().px(),part.at(i)->momentum().py(),part.at(i)->momentum().pz(),part.at(i)->momentum().e());
LVQ+=LV;
if(jak_id==TauDecay::JAK_A1_3PI ||
jak_id==TauDecay::JAK_KPIK ||
jak_id==TauDecay::JAK_KPIPI
){
if((tcharge==part.at(i)->pdg_id()/abs(part.at(i)->pdg_id()) && TD.nProng(TauBitMask)==3) || (jak_id==TauDecay::JAK_A1_3PI && TD.nProng(TauBitMask)==1 && abs(part.at(i)->pdg_id())==PdtPdgMini::pi_plus) ){
LVS13+=LV;
LVS23+=LV;
}
else{
LVS12+=LV;
if(!haspart1 && ((jak_id==TauDecay::JAK_A1_3PI) || (jak_id!=TauDecay::JAK_A1_3PI && abs(part.at(i)->pdg_id())==PdtPdgMini::K_plus) )){
LVS13+=LV;
haspart1=true;
}
else{
LVS23+=LV;
}
}
}
}
}
part.clear();
JAKInvMass.at(jak_id).at(0)->Fill(LVQ.M(),weight);
if(jak_id==TauDecay::JAK_A1_3PI ||
jak_id==TauDecay::JAK_KPIK ||
jak_id==TauDecay::JAK_KPIPI
){
JAKInvMass.at(jak_id).at(1)->Fill(LVS13.M(),weight);
JAKInvMass.at(jak_id).at(2)->Fill(LVS23.M(),weight);
JAKInvMass.at(jak_id).at(3)->Fill(LVS12.M(),weight);
}
}
}
else{
JAKID->Fill(jak_id,weight);
}
}
}
}
delete myGenEvent;
}//analyze
| void TauValidation::beginJob | ( | void | ) | [virtual] |
Setting the DQM top directories
Reimplemented from edm::EDAnalyzer.
Definition at line 34 of file TauValidation.cc.
References A0, a1, DQMStore::book1D(), D, dbe, electron, gamma, H0, Hpm, HSM, i, TauDecay::JAK_A1_3PI, TauDecay::JAK_KPIK, TauDecay::JAK_KPIPI, JAKID, JAKInvMass, K, Kstar, muon, NJAKID, nPrimeTaus, nTaus, other, pi, pi1pi0, pinpi0, rho, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), stable, TauBremPhotonsN, TauBremPhotonsPt, TauBremPhotonsPtSum, TauDecayChannels, TauEta, TauFSRPhotonsN, TauFSRPhotonsPt, TauFSRPhotonsPtSum, TauMothers, TauPhi, TauProngs, TauPt, TauRtauHpm, TauRtauW, TauSpinEffectsH_eX, TauSpinEffectsH_muX, TauSpinEffectsH_MVis, TauSpinEffectsH_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}");
TauRtauW = dbe->book1D("TauRtauW","W->Tau p(leading track)/E(visible tau)", 50 ,0,1); TauRtauW->setAxisTitle("rtau");
TauRtauHpm = dbe->book1D("TauRtauHpm","Hpm->Tau p(leading track)/E(visible tau)", 50 ,0,1); TauRtauHpm->setAxisTitle("rtau");
TauSpinEffectsW_X = dbe->book1D("TauSpinEffectsWX","Pion energy in W rest frame", 50 ,0,1); TauSpinEffectsW_X->setAxisTitle("X");
TauSpinEffectsHpm_X = dbe->book1D("TauSpinEffectsHpmX","Pion energy in Hpm rest frame", 50 ,0,1); TauSpinEffectsHpm_X->setAxisTitle("X");
TauSpinEffectsW_eX = dbe->book1D("TauSpinEffectsWeX","e energy in W rest frame", 50 ,0,1); TauSpinEffectsW_eX->setAxisTitle("X");
TauSpinEffectsHpm_eX = dbe->book1D("TauSpinEffectsHpmeX","e energy in Hpm rest frame", 50 ,0,1); TauSpinEffectsHpm_eX->setAxisTitle("X");
TauSpinEffectsW_muX = dbe->book1D("TauSpinEffectsWmuX","mu energy in W rest frame", 50 ,0,1); TauSpinEffectsW_muX->setAxisTitle("X");
TauSpinEffectsHpm_muX = dbe->book1D("TauSpinEffectsHpmmuX","mu energy in Hpm rest frame", 50 ,0,1); TauSpinEffectsHpm_muX->setAxisTitle("X");
TauSpinEffectsW_UpsilonRho = dbe->book1D("TauSpinEffectsWUpsilonRho","#Upsilon for #rho", 50 ,-1,1); TauSpinEffectsW_UpsilonRho->setAxisTitle("#Upsilon");
TauSpinEffectsHpm_UpsilonRho = dbe->book1D("TauSpinEffectsHpmUpsilonRho","#Upsilon for #rho", 50 ,-1,1); TauSpinEffectsHpm_UpsilonRho->setAxisTitle("#Upsilon");
TauSpinEffectsW_UpsilonA1 = dbe->book1D("TauSpinEffectsWUpsilonA1","#Upsilon for a1", 50 ,-1,1); TauSpinEffectsW_UpsilonA1->setAxisTitle("#Upsilon");
TauSpinEffectsHpm_UpsilonA1 = dbe->book1D("TauSpinEffectsHpmUpsilonA1","#Upsilon for a1", 50 ,-1,1); TauSpinEffectsHpm_UpsilonA1->setAxisTitle("#Upsilon");
TauSpinEffectsZ_MVis = dbe->book1D("TauSpinEffectsZMVis","Mass of pi+ pi-", 25 ,0,1.1); TauSpinEffectsZ_MVis->setAxisTitle("M_{#pi^{+}#pi^{-}}");
TauSpinEffectsH_MVis = dbe->book1D("TauSpinEffectsHMVis","Mass of pi+ pi-", 25 ,0,1.1); TauSpinEffectsZ_MVis->setAxisTitle("M_{#pi^{+}#pi^{-}}");
TauSpinEffectsZ_Zs = dbe->book1D("TauSpinEffectsZZs","Z_{s}", zsbins ,zsmin,zsmax); TauSpinEffectsZ_Zs->setAxisTitle("Z_{s}");
TauSpinEffectsH_Zs = dbe->book1D("TauSpinEffectsHZs","Z_{s}", zsbins ,zsmin,zsmax); TauSpinEffectsZ_Zs->setAxisTitle("Z_{s}");
TauSpinEffectsZ_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");
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 153 of file TauValidation.cc.
References fPDGTable, and edm::EventSetup::getData().
| void TauValidation::endJob | ( | void | ) | [virtual] |
| void TauValidation::endRun | ( | const edm::Run & | iRun, |
| const edm::EventSetup & | iSetup | ||
| ) | [virtual] |
| 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 327 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 290 of file TauValidation.cc.
References parents.
Referenced by spinEffectsZ(), and tauMother().
{
int mother_pid = 0;
if ( tau->production_vertex() ) {
HepMC::GenVertex::particle_iterator mother;
for (mother = tau->production_vertex()->particles_begin(HepMC::parents); mother!= tau->production_vertex()->particles_end(HepMC::parents); mother++ ) {
mother_pid = (*mother)->pdg_id();
if(mother_pid == tau->pdg_id()) return findMother(*mother); //mother_pid = -1; Make recursive to look for last tau in chain
}
}
return mother_pid;
}
| void TauValidation::FindPhotosFSR | ( | const HepMC::GenParticle * | p, |
| std::vector< const HepMC::GenParticle * > & | ListofFSR, | ||
| double & | BosonScale | ||
| ) | [private] |
Definition at line 351 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 315 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 265 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 277 of file TauValidation.cc.
References parents.
Referenced by spinEffectsZ().
{
std::vector<HepMC::GenParticle*> mothers;
if ( boson->production_vertex() ) {
HepMC::GenVertex::particle_iterator mother;
for (mother = boson->production_vertex()->particles_begin(HepMC::parents); mother!= boson->production_vertex()->particles_end(HepMC::parents); mother++ ) {
if((*mother)->pdg_id() == boson->pdg_id()) return GetMothers(*mother);
mothers.push_back(*mother);
}
}
return mothers;
}
| bool TauValidation::isLastTauinChain | ( | const HepMC::GenParticle * | tau | ) | [private] |
Definition at line 303 of file TauValidation.cc.
Referenced by analyze(), and findFSRandBrem().
| double TauValidation::leadingPionMomentum | ( | const HepMC::GenParticle * | tau, |
| double | weight | ||
| ) | [private] |
Definition at line 624 of file TauValidation.cc.
References leadingPionP4().
Referenced by rtau().
{
return leadingPionP4(tau).P();
}
| TLorentzVector TauValidation::leadingPionP4 | ( | const HepMC::GenParticle * | tau | ) | [private] |
Definition at line 628 of file TauValidation.cc.
References abs, p4, and evf::utils::pid.
Referenced by leadingPionMomentum(), spinEffects(), and spinEffectsZ().
{
TLorentzVector p4(0,0,0,0);
if ( tau->end_vertex() ) {
HepMC::GenVertex::particle_iterator des;
for(des = tau->end_vertex()->particles_begin(HepMC::descendants);
des!= tau->end_vertex()->particles_end(HepMC::descendants);++des ) {
int pid = (*des)->pdg_id();
if(abs(pid) == 15) return leadingPionP4(*des);
if(!(abs(pid)==211 || abs(pid)==13 || abs(pid)==11)) continue;
if((*des)->momentum().rho() > p4.P()) {
p4 = TLorentzVector((*des)->momentum().px(),
(*des)->momentum().py(),
(*des)->momentum().pz(),
(*des)->momentum().e());
}
}
}
return p4;
}
| TLorentzVector TauValidation::motherP4 | ( | const HepMC::GenParticle * | tau | ) | [private] |
Definition at line 654 of file TauValidation.cc.
References configurableAnalysis::GenParticle, GetMother(), and m.
Referenced by spinEffects().
{
const HepMC::GenParticle* m=GetMother(tau);
return TLorentzVector(m->momentum().px(),m->momentum().py(),m->momentum().pz(),m->momentum().e());
}
| void TauValidation::photons | ( | const HepMC::GenParticle * | tau, |
| double | weight | ||
| ) | [private] |
Definition at line 685 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 472 of file TauValidation.cc.
References abs, MonitorElement::Fill(), leadingPionMomentum(), pi1pi0, tauEtCut, TauRtauHpm, TauRtauW, and visibleTauEnergy().
Referenced by analyze().
{
if(decay != pi1pi0) return; // polarization only for 1-prong hadronic taus with one neutral pion to make a clean case
if(tau->momentum().perp() < tauEtCut) return; // rtau visible only for boosted taus
double rTau = 0;
double ltrack = leadingPionMomentum(tau, weight);
double visibleTauE = visibleTauEnergy(tau);
if(visibleTauE != 0) rTau = ltrack/visibleTauE;
if(abs(mother) == 24) TauRtauW->Fill(rTau,weight);
if(abs(mother) == 37) TauRtauHpm->Fill(rTau,weight);
}
| void TauValidation::spinEffects | ( | const HepMC::GenParticle * | tau, |
| int | mother, | ||
| int | decay, | ||
| std::vector< HepMC::GenParticle * > & | part, | ||
| double | weight | ||
| ) | [private] |
Definition at line 488 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::spinEffectsZ | ( | const HepMC::GenParticle * | boson, |
| double | weight | ||
| ) | [private] |
Definition at line 537 of file TauValidation.cc.
References abs, PdtPdgMini::anti_d, PdtPdgMini::anti_u, DeDxDiscriminatorTools::charge(), PdtPdgMini::d, electron, MonitorElement::Fill(), findMother(), fPDGTable, GetMothers(), PdtPdgMini::Higgs0, i, leadingPionP4(), m, muon, pi, evf::utils::pid, lumiQueryAPI::q, tauDecayChannel(), TauSpinEffectsH_eX, TauSpinEffectsH_muX, TauSpinEffectsH_MVis, TauSpinEffectsH_X, TauSpinEffectsH_Xb, TauSpinEffectsH_Xf, TauSpinEffectsH_Zs, TauSpinEffectsZ_eX, TauSpinEffectsZ_muX, TauSpinEffectsZ_MVis, TauSpinEffectsZ_X, TauSpinEffectsZ_Xb, TauSpinEffectsZ_Xf, TauSpinEffectsZ_Zs, PdtPdgMini::u, CommonMethods::weight(), Z, PdtPdgMini::Z0, zsbins, zsmax, zsmin, and Zstoa().
Referenced by analyze().
{
TLorentzVector tautau(0,0,0,0);
TLorentzVector pipi(0,0,0,0);
TLorentzVector taum(0,0,0,0);
int nSinglePionDecays(0),nSingleMuonDecays(0),nSingleElectronDecays(0);
double x1(0),x2(0);
TLorentzVector Zboson(boson->momentum().px(),boson->momentum().py(),boson->momentum().pz(),boson->momentum().e());
if ( boson->end_vertex() ) {
HepMC::GenVertex::particle_iterator des;
for(des = boson->end_vertex()->particles_begin(HepMC::children); des!= boson->end_vertex()->particles_end(HepMC::children);++des ) {
int pid = (*des)->pdg_id();
if(abs(findMother(*des)) != 15 && abs(pid) == 15 && (tauDecayChannel(*des) == pi || tauDecayChannel(*des) == muon || tauDecayChannel(*des) == electron )){
if(tauDecayChannel(*des) == pi)nSinglePionDecays++;
if(tauDecayChannel(*des) == muon)nSingleMuonDecays++;
if(tauDecayChannel(*des) == electron)nSingleElectronDecays++;
TLorentzVector LVtau((*des)->momentum().px(),(*des)->momentum().py(),(*des)->momentum().pz(),(*des)->momentum().e());
tautau += LVtau;
TLorentzVector LVpi=leadingPionP4(*des);
pipi+=LVpi;
const HepPDT::ParticleData* pd = fPDGTable->particle((*des)->pdg_id ());
int charge = (int) pd->charge();
LVtau.Boost(-1*Zboson.BoostVector());
LVpi.Boost(-1*Zboson.BoostVector());
if(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(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 419 of file TauValidation.cc.
References a1, abs, electron, MonitorElement::Fill(), K, Kstar, muon, pi, pi1pi0, evf::utils::pid, pinpi0, rho, stable, TauDecayChannels, tripi, tripinpi0, and undetermined.
Referenced by analyze(), and spinEffectsZ().
{
int channel = undetermined;
if(tau->status() == 1) channel = stable;
int allCount = 0,
eCount = 0,
muCount = 0,
pi0Count = 0,
piCount = 0,
rhoCount = 0,
a1Count = 0,
KCount = 0,
KstarCount = 0;
if ( tau->end_vertex() ) {
HepMC::GenVertex::particle_iterator des;
for(des = tau->end_vertex()->particles_begin(HepMC::descendants);
des!= tau->end_vertex()->particles_end(HepMC::descendants);++des ) {
int pid = (*des)->pdg_id();
if(abs(pid) == 15) return tauDecayChannel(*des,weight);
allCount++;
if(abs(pid) == 11) eCount++;
if(abs(pid) == 13) muCount++;
if(abs(pid) == 111) pi0Count++;
if(abs(pid) == 211) piCount++;
if(abs(pid) == 213) rhoCount++;
if(abs(pid) == 20213) a1Count++;
if(abs(pid) == 321) KCount++;
if(abs(pid) == 323) KstarCount++;
}
}
// resonances
if(KCount >= 1) channel = K;
if(KstarCount >= 1) channel = Kstar;
if(a1Count >= 1) channel = a1;
if(rhoCount >= 1) channel = rho;
if(channel!=undetermined && weight!=0.0) TauDecayChannels->Fill(channel,weight);
// final state products
if(piCount == 1 && pi0Count == 0) channel = pi;
if(piCount == 1 && pi0Count == 1) channel = pi1pi0;
if(piCount == 1 && pi0Count > 1) channel = pinpi0;
if(piCount == 3 && pi0Count == 0) channel = tripi;
if(piCount == 3 && pi0Count > 0) channel = tripinpi0;
if(eCount == 1) channel = electron;
if(muCount == 1) channel = muon;
if(weight!=0.0) TauDecayChannels->Fill(channel,weight);
return channel;
}
| int TauValidation::tauMother | ( | const HepMC::GenParticle * | tau, |
| double | weight | ||
| ) | [private] |
Definition at line 374 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 399 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 659 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 615 of file TauValidation.cc.
References a, and mathSSE::sqrt().
Referenced by spinEffectsZ().
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().
PDT table.
Definition at line 103 of file TauValidation.h.
Referenced by beginRun(), spinEffectsZ(), and tauProngs().
edm::InputTag TauValidation::hepmcCollection_ [private] |
Definition at line 98 of file TauValidation.h.
Referenced by analyze().
MonitorElement* TauValidation::JAKID [private] |
Definition at line 119 of file TauValidation.h.
Referenced by analyze(), and beginJob().
std::vector<std::vector<MonitorElement *> > TauValidation::JAKInvMass [private] |
Definition at line 120 of file TauValidation.h.
Referenced by analyze(), and beginJob().
unsigned int TauValidation::NJAKID [private] |
Definition at line 118 of file TauValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement * TauValidation::nPrimeTaus [private] |
Definition at line 108 of file TauValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement* TauValidation::nTaus [private] |
Definition at line 108 of file TauValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement * TauValidation::TauBremPhotonsN [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and photons().
MonitorElement * TauValidation::TauBremPhotonsPt [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and photons().
MonitorElement * TauValidation::TauBremPhotonsPtSum [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and photons().
MonitorElement * TauValidation::TauDecayChannels [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and tauDecayChannel().
MonitorElement * TauValidation::TauEta [private] |
Definition at line 109 of file TauValidation.h.
Referenced by analyze(), and beginJob().
double TauValidation::tauEtCut [private] |
Definition at line 100 of file TauValidation.h.
Referenced by rtau().
MonitorElement * TauValidation::TauFSRPhotonsN [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and photons().
MonitorElement * TauValidation::TauFSRPhotonsPt [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and photons().
MonitorElement * TauValidation::TauFSRPhotonsPtSum [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and photons().
MonitorElement * TauValidation::TauMothers [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and tauMother().
MonitorElement * TauValidation::TauPhi [private] |
Definition at line 109 of file TauValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement * TauValidation::TauProngs [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and tauProngs().
MonitorElement* TauValidation::TauPt [private] |
Definition at line 109 of file TauValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement * TauValidation::TauRtauHpm [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and rtau().
MonitorElement * TauValidation::TauRtauW [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and rtau().
MonitorElement * TauValidation::TauSpinEffectsH_eX [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffectsZ().
MonitorElement * TauValidation::TauSpinEffectsH_muX [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffectsZ().
MonitorElement * TauValidation::TauSpinEffectsH_MVis [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffectsZ().
MonitorElement * TauValidation::TauSpinEffectsH_X [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffectsZ().
MonitorElement * TauValidation::TauSpinEffectsH_Xb [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffectsZ().
MonitorElement * TauValidation::TauSpinEffectsH_Xf [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffectsZ().
MonitorElement * TauValidation::TauSpinEffectsH_Zs [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffectsZ().
MonitorElement * TauValidation::TauSpinEffectsHpm_eX [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffects().
MonitorElement * TauValidation::TauSpinEffectsHpm_muX [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffects().
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffects().
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffects().
MonitorElement * TauValidation::TauSpinEffectsHpm_X [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffects().
MonitorElement * TauValidation::TauSpinEffectsW_eX [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffects().
MonitorElement * TauValidation::TauSpinEffectsW_muX [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffects().
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffects().
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffects().
MonitorElement * TauValidation::TauSpinEffectsW_X [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffects().
MonitorElement * TauValidation::TauSpinEffectsZ_eX [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffectsZ().
MonitorElement * TauValidation::TauSpinEffectsZ_muX [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffectsZ().
MonitorElement * TauValidation::TauSpinEffectsZ_MVis [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffectsZ().
MonitorElement * TauValidation::TauSpinEffectsZ_X [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffectsZ().
MonitorElement * TauValidation::TauSpinEffectsZ_Xb [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffectsZ().
MonitorElement * TauValidation::TauSpinEffectsZ_Xf [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffectsZ().
MonitorElement * TauValidation::TauSpinEffectsZ_Zs [private] |
Definition at line 109 of file TauValidation.h.
Referenced by beginJob(), and spinEffectsZ().
int TauValidation::zsbins [private] |
Definition at line 122 of file TauValidation.h.
Referenced by beginJob(), and spinEffectsZ().
double TauValidation::zsmax [private] |
Definition at line 123 of file TauValidation.h.
Referenced by beginJob(), and spinEffectsZ().
double TauValidation::zsmin [private] |
Definition at line 123 of file TauValidation.h.
Referenced by beginJob(), and spinEffectsZ().
1.7.3