|
|
|
#include <TauValidation.h>
Public Types | |
| enum | { undetermined, electron, muon, pi, K, pi1pi0, pinpi0, tripi, tripinpi0, stable } |
| enum | { other, gamma, Z, W, HSM, H0, A0, Hpm } |
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 () |
Private Member Functions | |
| double | leadingPionMomentum (const HepMC::GenParticle *) |
| TLorentzVector | leadingPionP4 (const HepMC::GenParticle *) |
| TLorentzVector | motherP4 (const HepMC::GenParticle *) |
| void | photons (const HepMC::GenParticle *) |
| void | rtau (const HepMC::GenParticle *, int, int) |
| void | spinEffects (const HepMC::GenParticle *, int, int) |
| int | tauDecayChannel (const HepMC::GenParticle *) |
| int | tauMother (const HepMC::GenParticle *) |
| int | tauProngs (const HepMC::GenParticle *) |
| double | visibleTauEnergy (const HepMC::GenParticle *) |
Private Attributes | |
| DQMStore * | dbe |
| ME's "container". | |
| edm::ESHandle < HepPDT::ParticleDataTable > | fPDGTable |
| PDT table. | |
| edm::InputTag | hepmcCollection_ |
| MonitorElement * | nEvt |
| int | nTaus |
| int | nTausWithPhotons |
| double | photonFromTauPtSum |
| MonitorElement * | TauDecayChannels |
| MonitorElement * | TauEta |
| double | tauEtCut |
| MonitorElement * | TauMothers |
| MonitorElement * | TauPhotons |
| MonitorElement * | TauProngs |
| MonitorElement * | TauPt |
| double | tauPtSum |
| MonitorElement * | TauRtauHpm |
| MonitorElement * | TauRtauW |
| MonitorElement * | TauSpinEffectsHpm |
| MonitorElement * | TauSpinEffectsW |
Definition at line 34 of file TauValidation.h.
| anonymous enum |
| anonymous enum |
| TauValidation::TauValidation | ( | const edm::ParameterSet & | iPSet | ) | [explicit] |
Definition at line 20 of file TauValidation.cc.
References dbe, and cmsCodeRules::cppFunctionSkipper::operator.
: hepmcCollection_(iPSet.getParameter<edm::InputTag>("hepmcCollection")), tauEtCut(iPSet.getParameter<double>("tauEtCutForRtau")) { dbe = 0; dbe = edm::Service<DQMStore>().operator->(); }
| TauValidation::~TauValidation | ( | ) | [virtual] |
Definition at line 28 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 98 of file TauValidation.cc.
References abs, MonitorElement::Fill(), edm::Event::getByLabel(), hepmcCollection_, nEvt, photons(), rtau(), spinEffects(), tauDecayChannel(), TauEta, tauMother(), tauProngs(), and TauPt.
{
edm::Handle<HepMCProduct> evt;
iEvent.getByLabel(hepmcCollection_, evt);
//Get EVENT
HepMC::GenEvent *myGenEvent = new HepMC::GenEvent(*(evt->GetEvent()));
nEvt->Fill(0.5);
// find taus
for(HepMC::GenEvent::particle_const_iterator iter = myGenEvent->particles_begin(); iter != myGenEvent->particles_end(); ++iter) {
if((*iter)->status()==3) {
if(abs((*iter)->pdg_id())==15){
TauPt->Fill((*iter)->momentum().perp());
TauEta->Fill((*iter)->momentum().eta());
int mother = tauMother(*iter);
int decaychannel = tauDecayChannel(*iter);
tauProngs(*iter);
rtau(*iter,mother,decaychannel);
spinEffects(*iter,mother,decaychannel);
photons(*iter);
}
}
}
delete myGenEvent;
}//analyze
| void TauValidation::beginJob | ( | void | ) | [virtual] |
Setting the DQM top directories
Reimplemented from edm::EDAnalyzer.
Definition at line 30 of file TauValidation.cc.
References A0, DQMStore::book1D(), dbe, electron, gamma, H0, Hpm, HSM, K, muon, nEvt, nTaus, nTausWithPhotons, other, photonFromTauPtSum, pi, pi1pi0, pinpi0, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), stable, TauDecayChannels, TauEta, TauMothers, TauPhotons, TauProngs, TauPt, tauPtSum, TauRtauHpm, TauRtauW, TauSpinEffectsHpm, TauSpinEffectsW, tripi, tripinpi0, undetermined, W, and Z.
{
if(dbe){
dbe->setCurrentFolder("Generator/Tau");
// Number of analyzed events
nEvt = dbe->book1D("nEvt", "n analyzed Events", 1, 0., 1.);
//Kinematics
TauPt = dbe->book1D("TauPt","Tau pT", 100 ,0,100);
TauEta = dbe->book1D("TauEta","Tau eta", 100 ,-2.5,2.5);
TauProngs = dbe->book1D("TauProngs","Tau n prongs", 7 ,0,7);
TauDecayChannels = dbe->book1D("TauDecayChannels","Tau decay channels", 10 ,0,10);
TauDecayChannels->setBinLabel(1+undetermined,"?");
TauDecayChannels->setBinLabel(1+electron,"e");
TauDecayChannels->setBinLabel(1+muon,"mu");
TauDecayChannels->setBinLabel(1+pi,"#pi^{#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+stable,"Stable");
TauMothers = dbe->book1D("TauMothers","Tau mother particles", 10 ,0,10);
TauMothers->setBinLabel(1+other,"?");
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 = dbe->book1D("TauSpinEffectsW","Pion energy in W rest frame", 50 ,0,1);
TauSpinEffectsW->setAxisTitle("Energy");
TauSpinEffectsHpm = dbe->book1D("TauSpinEffectsHpm","Pion energy in Hpm rest frame", 50 ,0,1);
TauSpinEffectsHpm->setAxisTitle("Energy");
TauPhotons = dbe->book1D("TauPhoton","Photons radiating from tau", 2 ,0,2);
TauPhotons->setBinLabel(1,"Fraction of taus radiating photons");
TauPhotons->setBinLabel(2,"Fraction of tau pt radiated by photons");
}
nTaus = 0;
nTausWithPhotons = 0;
tauPtSum = 0;
photonFromTauPtSum = 0;
return;
}
| void TauValidation::beginRun | ( | const edm::Run & | iRun, |
| const edm::EventSetup & | iSetup | ||
| ) | [virtual] |
Get PDT Table
Reimplemented from edm::EDAnalyzer.
Definition at line 91 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] |
| double TauValidation::leadingPionMomentum | ( | const HepMC::GenParticle * | tau | ) | [private] |
Definition at line 255 of file TauValidation.cc.
References leadingPionP4().
Referenced by rtau().
{
return leadingPionP4(tau).P();
}
| TLorentzVector TauValidation::leadingPionP4 | ( | const HepMC::GenParticle * | tau | ) | [private] |
Definition at line 259 of file TauValidation.cc.
References abs, p4, and evf::utils::pid.
Referenced by leadingPionMomentum(), and spinEffects().
{
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) 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 285 of file TauValidation.cc.
Referenced by spinEffects().
{
TLorentzVector p4(0,0,0,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();
//std::cout << " parent " << mother_pid << std::endl;
p4 = TLorentzVector((*mother)->momentum().px(),
(*mother)->momentum().py(),
(*mother)->momentum().pz(),
(*mother)->momentum().e());
}
}
return p4;
}
| void TauValidation::photons | ( | const HepMC::GenParticle * | tau | ) | [private] |
Definition at line 331 of file TauValidation.cc.
References nTaus, nTausWithPhotons, photonFromTauPtSum, evf::utils::pid, MonitorElement::setBinContent(), TauPhotons, and tauPtSum.
Referenced by analyze().
{
if ( tau->end_vertex() ) {
bool photonFromTau = false;
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(pid == 22) {
photonFromTauPtSum += (*des)->momentum().perp();
photonFromTau = true;
}
}
nTaus++;
if(photonFromTau) {
tauPtSum += tau->momentum().perp();
nTausWithPhotons++;
}
double nFrac = double(nTausWithPhotons)/nTaus;
double ptFrac = 0;
if(tauPtSum > 0) ptFrac = photonFromTauPtSum/tauPtSum;
TauPhotons->setBinContent(1,nFrac);
TauPhotons->setBinContent(2,ptFrac);
}
}
| void TauValidation::rtau | ( | const HepMC::GenParticle * | tau, |
| int | mother, | ||
| int | decay | ||
| ) | [private] |
Definition at line 224 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);
double visibleTauE = visibleTauEnergy(tau);
if(visibleTauE != 0) rTau = ltrack/visibleTauE;
if(abs(mother) == 23) TauRtauW->Fill(rTau);
if(abs(mother) == 37) TauRtauHpm->Fill(rTau);
}
| void TauValidation::spinEffects | ( | const HepMC::GenParticle * | tau, |
| int | mother, | ||
| int | decay | ||
| ) | [private] |
Definition at line 240 of file TauValidation.cc.
References abs, relval_parameters_module::energy, MonitorElement::Fill(), leadingPionP4(), motherP4(), pi, TauSpinEffectsHpm, and TauSpinEffectsW.
Referenced by analyze().
{
if(decay != pi) return; // 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(abs(mother) == 23) TauSpinEffectsW->Fill(energy);
if(abs(mother) == 37) TauSpinEffectsHpm->Fill(energy);
}
| int TauValidation::tauDecayChannel | ( | const HepMC::GenParticle * | tau | ) | [private] |
Definition at line 179 of file TauValidation.cc.
References abs, electron, MonitorElement::Fill(), muon, pi, pi1pi0, evf::utils::pid, pinpi0, stable, TauDecayChannels, tripi, tripinpi0, and undetermined.
Referenced by analyze().
{
int channel = undetermined;
if(tau->status() == 1) channel = stable;
int eCount = 0,
muCount = 0,
pi0Count = 0,
piCount = 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 ) {
//if(abs(tauMother(*des)) != 15) continue;
int pid = (*des)->pdg_id();
//std::cout << " barcode=" << (*des)->barcode() << " pid="
// << pid << " mom=" << tauMother(*des) << " status="
// << (*des)->status() << std::endl;
if(abs(pid) == 15) return tauDecayChannel(*des);
if(abs(pid) == 11) eCount++;
if(abs(pid) == 13) muCount++;
if(abs(pid) == 111) pi0Count++;
if(abs(pid) == 211) piCount++;
}
}
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;
TauDecayChannels->Fill(channel);
return channel;
}
| int TauValidation::tauMother | ( | const HepMC::GenParticle * | tau | ) | [private] |
Definition at line 129 of file TauValidation.cc.
References A0, abs, MonitorElement::Fill(), gamma, H0, Hpm, HSM, label, other, parents, TauMothers, W, and Z.
Referenced by analyze().
{
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();
//std::cout << " parent " << mother_pid << std::endl;
}
}
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;
TauMothers->Fill(label);
return mother_pid;
}
| int TauValidation::tauProngs | ( | const HepMC::GenParticle * | tau | ) | [private] |
Definition at line 155 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);
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;
//std::cout << "TauValidation::tauProngs barcode=" << (*des)->barcode() << " pid="
// << pid << " mom=" << tauMother(*des) << " status="
// << (*des)->status() << " charge=" << charge << std::endl;
nProngs++;
}
}
TauProngs->Fill(nProngs);
return nProngs;
}
| double TauValidation::visibleTauEnergy | ( | const HepMC::GenParticle * | tau | ) | [private] |
Definition at line 305 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();
}
DQMStore* TauValidation::dbe [private] |
ME's "container".
Definition at line 90 of file TauValidation.h.
Referenced by beginJob(), and TauValidation().
PDT table.
Definition at line 87 of file TauValidation.h.
Referenced by beginRun(), and tauProngs().
edm::InputTag TauValidation::hepmcCollection_ [private] |
Definition at line 79 of file TauValidation.h.
Referenced by analyze().
MonitorElement* TauValidation::nEvt [private] |
Definition at line 92 of file TauValidation.h.
Referenced by analyze(), and beginJob().
int TauValidation::nTaus [private] |
Definition at line 84 of file TauValidation.h.
Referenced by beginJob(), and photons().
int TauValidation::nTausWithPhotons [private] |
Definition at line 84 of file TauValidation.h.
Referenced by beginJob(), and photons().
double TauValidation::photonFromTauPtSum [private] |
Definition at line 83 of file TauValidation.h.
Referenced by beginJob(), and photons().
MonitorElement * TauValidation::TauDecayChannels [private] |
Definition at line 93 of file TauValidation.h.
Referenced by beginJob(), and tauDecayChannel().
MonitorElement * TauValidation::TauEta [private] |
Definition at line 93 of file TauValidation.h.
Referenced by analyze(), and beginJob().
double TauValidation::tauEtCut [private] |
Definition at line 81 of file TauValidation.h.
Referenced by rtau().
MonitorElement * TauValidation::TauMothers [private] |
Definition at line 93 of file TauValidation.h.
Referenced by beginJob(), and tauMother().
MonitorElement * TauValidation::TauPhotons [private] |
Definition at line 93 of file TauValidation.h.
Referenced by beginJob(), and photons().
MonitorElement * TauValidation::TauProngs [private] |
Definition at line 93 of file TauValidation.h.
Referenced by beginJob(), and tauProngs().
MonitorElement* TauValidation::TauPt [private] |
Definition at line 93 of file TauValidation.h.
Referenced by analyze(), and beginJob().
double TauValidation::tauPtSum [private] |
Definition at line 83 of file TauValidation.h.
Referenced by beginJob(), and photons().
MonitorElement * TauValidation::TauRtauHpm [private] |
Definition at line 93 of file TauValidation.h.
Referenced by beginJob(), and rtau().
MonitorElement * TauValidation::TauRtauW [private] |
Definition at line 93 of file TauValidation.h.
Referenced by beginJob(), and rtau().
MonitorElement * TauValidation::TauSpinEffectsHpm [private] |
Definition at line 93 of file TauValidation.h.
Referenced by beginJob(), and spinEffects().
MonitorElement * TauValidation::TauSpinEffectsW [private] |
Definition at line 93 of file TauValidation.h.
Referenced by beginJob(), and spinEffects().
1.7.3