|
|
|
#include <TauValidation.h>
Public Types | |
| enum | { undetermined, electron, muon, pi, rho, a1, K, Kstar, pi1pi0, pinpi0, tripi, stable } |
| enum | { other, B, D, 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 | |
| int | findMother (const HepMC::GenParticle *) |
| int | findTauDecayChannel (const HepMC::GenParticle *) |
| double | leadingPionMomentum (const HepMC::GenParticle *, double weight) |
| TLorentzVector | leadingPionP4 (const HepMC::GenParticle *) |
| TLorentzVector | motherP4 (const HepMC::GenParticle *) |
| void | photons (const HepMC::GenParticle *, double weight) |
| void | rtau (const HepMC::GenParticle *, int, int, double weight) |
| void | spinEffects (const HepMC::GenParticle *, int, int, double weight) |
| void | spinEffectsZ (const HepMC::GenParticle *, double weight) |
| int | tauDecayChannel (const HepMC::GenParticle *, double weight) |
| int | tauMother (const HepMC::GenParticle *, double weight) |
| int | tauProngs (const HepMC::GenParticle *, double weight) |
| double | visibleTauEnergy (const HepMC::GenParticle *) |
Private Attributes | |
| WeightManager | _wmanager |
| DQMStore * | dbe |
| ME's "container". | |
| edm::ESHandle < HepPDT::ParticleDataTable > | fPDGTable |
| PDT table. | |
| edm::InputTag | hepmcCollection_ |
| MonitorElement * | JAKID |
| std::vector< std::vector < MonitorElement * > > | JAKInvMass |
| MonitorElement * | nEvt |
| unsigned int | NJAKID |
| MonitorElement * | TauDecayChannels |
| MonitorElement * | TauEta |
| double | tauEtCut |
| MonitorElement * | TauMothers |
| MonitorElement * | TauPhi |
| MonitorElement * | TauPhotonsN |
| MonitorElement * | TauPhotonsPt |
| MonitorElement * | TauProngs |
| MonitorElement * | TauPt |
| MonitorElement * | TauRtauHpm |
| MonitorElement * | TauRtauW |
| MonitorElement * | TauSpinEffectsHpm |
| MonitorElement * | TauSpinEffectsW |
| MonitorElement * | TauSpinEffectsZ |
Definition at line 36 of file TauValidation.h.
| anonymous enum |
| anonymous enum |
| TauValidation::TauValidation | ( | const edm::ParameterSet & | iPSet | ) | [explicit] |
Definition at line 22 of file TauValidation.cc.
References dbe, and cppFunctionSkipper::operator.
: _wmanager(iPSet) ,hepmcCollection_(iPSet.getParameter<edm::InputTag>("hepmcCollection")) ,tauEtCut(iPSet.getParameter<double>("tauEtCutForRtau")) ,NJAKID(22) { 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 124 of file TauValidation.cc.
References _wmanager, abs, TauDecay_CMSSW::AnalyzeTau(), MonitorElement::Fill(), TauDecay_CMSSW::Get_TauDecayProducts(), edm::Event::getByLabel(), hepmcCollection_, i, TauDecay::isTauFinalStateParticle(), TauDecay::JAK_A1_3PI, TauDecay::JAK_KPIK, TauDecay::JAK_KPIPI, JAKID, JAKInvMass, PdtPdgMini::K_plus, nEvt, NJAKID, TauDecay::nProng(), 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);
nEvt->Fill(0.5,weight);
// find taus
for(HepMC::GenEvent::particle_const_iterator iter = myGenEvent->particles_begin(); iter != myGenEvent->particles_end(); ++iter) {
if(abs((*iter)->pdg_id())==15){
if(tauMother(*iter,weight)!=-1){ // exclude B, D and other non-signal decay modes
TauPt->Fill((*iter)->momentum().perp(),weight);
TauEta->Fill((*iter)->momentum().eta(),weight);
TauPhi->Fill((*iter)->momentum().phi(),weight);
int mother = tauMother(*iter,weight);
int decaychannel = tauDecayChannel(*iter,weight);
tauProngs(*iter, weight);
rtau(*iter,mother,decaychannel,weight);
spinEffects(*iter,mother,decaychannel,weight);
photons(*iter,weight);
}
if(abs((*iter)->pdg_id())==23){
spinEffectsZ(*iter,weight);
}
//Adding JAKID and Mass information
//
TauDecay_CMSSW TD;
unsigned int jak_id, TauBitMask;
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();
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);
}
}
}
}
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, nEvt, NJAKID, other, pi, pi1pi0, pinpi0, rho, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), stable, TauDecayChannels, TauEta, TauMothers, TauPhi, TauPhotonsN, TauPhotonsPt, TauProngs, TauPt, TauRtauHpm, TauRtauW, TauSpinEffectsHpm, TauSpinEffectsW, TauSpinEffectsZ, tmp, tripi, 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);
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", 12 ,0,12);
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 = 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");
TauSpinEffectsZ = dbe->book1D("TauSpinEffectsZ","Mass of pi+ pi-", 22 ,0,1.1);
TauSpinEffectsZ->setAxisTitle("M_{#pi^{+}#pi^{-}}");
TauPhotonsN = dbe->book1D("TauPhotonsN","Photons radiating from tau", 2 ,0,2);
TauPhotonsN->setBinLabel(1,"Number of taus");
TauPhotonsN->setBinLabel(2,"Number of taus radiating photons");
TauPhotonsPt = dbe->book1D("TauPhotonsPt","Photon pt radiating from tau", 2 ,0,2);
TauPhotonsPt->setBinLabel(1,"Sum of tau pt");
TauPhotonsPt->setBinLabel(2,"Sum of tau pt radiated by photons");
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 117 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] |
| int TauValidation::findMother | ( | const HepMC::GenParticle * | tau | ) | [private] |
Definition at line 215 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()) mother_pid = -1;//findMother(*mother);
//std::cout << " parent " << mother_pid << std::endl;
}
}
return mother_pid;
}
| int TauValidation::findTauDecayChannel | ( | const HepMC::GenParticle * | tau | ) | [private] |
Definition at line 279 of file TauValidation.cc.
References a1, abs, electron, K, Kstar, muon, pi, pi1pi0, evf::utils::pid, pinpi0, rho, stable, tripi, and undetermined.
Referenced by spinEffectsZ(), and tauDecayChannel().
{
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 ) {
//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 findTauDecayChannel(*des);
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++;
}
}
if(KCount == 1 && allCount == 2) channel = K;
if(KstarCount == 1 && allCount == 2) channel = Kstar;
if(a1Count == 1 && allCount == 2) channel = a1;
if(rhoCount == 1 && allCount == 2) channel = rho;
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;
return channel;
}
| double TauValidation::leadingPionMomentum | ( | const HepMC::GenParticle * | tau, |
| double | weight | ||
| ) | [private] |
Definition at line 409 of file TauValidation.cc.
References leadingPionP4().
Referenced by rtau().
{
return leadingPionP4(tau).P();
}
| TLorentzVector TauValidation::leadingPionP4 | ( | const HepMC::GenParticle * | tau | ) | [private] |
Definition at line 413 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) 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 439 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, |
| double | weight | ||
| ) | [private] |
Definition at line 485 of file TauValidation.cc.
References MonitorElement::Fill(), evf::utils::pid, TauPhotonsN, and TauPhotonsPt.
Referenced by analyze().
{
if ( tau->end_vertex() ) {
double photonFromTauPtSum = 0;
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;
}
}
TauPhotonsN->Fill(0.5,weight);
//doesn't seems like it makes sense to use a weight below
TauPhotonsPt->Fill(0.5,tau->momentum().perp());
if(photonFromTau) {
TauPhotonsN->Fill(1.5,weight);
//doesn't seems like it makes sense to use a weight below
TauPhotonsPt->Fill(1.5,photonFromTauPtSum);
}
}
}
| void TauValidation::rtau | ( | const HepMC::GenParticle * | tau, |
| int | mother, | ||
| int | decay, | ||
| double | weight | ||
| ) | [private] |
Definition at line 344 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, | ||
| double | weight | ||
| ) | [private] |
Definition at line 360 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) == 24) TauSpinEffectsW->Fill(energy,weight);
if(abs(mother) == 37) TauSpinEffectsHpm->Fill(energy,weight);
}
| void TauValidation::spinEffectsZ | ( | const HepMC::GenParticle * | boson, |
| double | weight | ||
| ) | [private] |
Definition at line 375 of file TauValidation.cc.
References abs, MonitorElement::Fill(), findMother(), findTauDecayChannel(), leadingPionP4(), pi, evf::utils::pid, TauSpinEffectsZ, and CommonMethods::weight().
Referenced by analyze().
{
TLorentzVector tautau(0,0,0,0);
TLorentzVector pipi(0,0,0,0);
int nSinglePionDecays = 0;
if ( boson->end_vertex() ) {
HepMC::GenVertex::particle_iterator des;
for(des = boson->end_vertex()->particles_begin(HepMC::descendants);
des!= boson->end_vertex()->particles_end(HepMC::descendants);++des ) {
int pid = (*des)->pdg_id();
/*std::cout << " barcode=" << (*des)->barcode() << " pid="
<< pid << " mom=" << findMother(*des) << " status="
<< (*des)->status() << " px="
<< (*des)->momentum().px() << " decay="
<< findTauDecayChannel(*des) << std::endl;
*/
if(abs(findMother(*des)) != 15 &&
abs(pid) == 15 && findTauDecayChannel(*des) == pi){
nSinglePionDecays++;
tautau += TLorentzVector((*des)->momentum().px(),
(*des)->momentum().py(),
(*des)->momentum().pz(),
(*des)->momentum().e());
pipi += leadingPionP4(*des);
}
}
}
if(nSinglePionDecays == 2 && tautau.M() != 0) {
TauSpinEffectsZ->Fill(pipi.M()/tautau.M(),weight);
}
}
| int TauValidation::tauDecayChannel | ( | const HepMC::GenParticle * | tau, |
| double | weight | ||
| ) | [private] |
Definition at line 338 of file TauValidation.cc.
References MonitorElement::Fill(), findTauDecayChannel(), and TauDecayChannels.
Referenced by analyze().
{
int channel = findTauDecayChannel(tau);
TauDecayChannels->Fill(channel,weight);
return channel;
}
| int TauValidation::tauMother | ( | const HepMC::GenParticle * | tau, |
| double | weight | ||
| ) | [private] |
Definition at line 230 of file TauValidation.cc.
References A0, abs, B, D, MonitorElement::Fill(), findMother(), gamma, H0, Hpm, HSM, label, other, TauMothers, W, and Z.
Referenced by analyze().
{
if(abs(tau->pdg_id()) != 15 ) return -1;
int mother_pid = findMother(tau);
if(mother_pid == -1) return -1;
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 255 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;
//std::cout << "TauValidation::tauProngs barcode=" << (*des)->barcode() << " pid="
// << pid << " mom=" << tauMother(*des) << " status="
// << (*des)->status() << " charge=" << charge << std::endl;
nProngs++;
}
}
TauProngs->Fill(nProngs,weight);
return nProngs;
}
| double TauValidation::visibleTauEnergy | ( | const HepMC::GenParticle * | tau | ) | [private] |
Definition at line 459 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();
}
WeightManager TauValidation::_wmanager [private] |
Definition at line 89 of file TauValidation.h.
Referenced by analyze().
DQMStore* TauValidation::dbe [private] |
ME's "container".
Definition at line 99 of file TauValidation.h.
Referenced by beginJob(), and TauValidation().
PDT table.
Definition at line 96 of file TauValidation.h.
Referenced by beginRun(), and tauProngs().
edm::InputTag TauValidation::hepmcCollection_ [private] |
Definition at line 91 of file TauValidation.h.
Referenced by analyze().
MonitorElement* TauValidation::JAKID [private] |
Definition at line 107 of file TauValidation.h.
Referenced by analyze(), and beginJob().
std::vector<std::vector<MonitorElement *> > TauValidation::JAKInvMass [private] |
Definition at line 108 of file TauValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement* TauValidation::nEvt [private] |
Definition at line 101 of file TauValidation.h.
Referenced by analyze(), and beginJob().
unsigned int TauValidation::NJAKID [private] |
Definition at line 106 of file TauValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement * TauValidation::TauDecayChannels [private] |
Definition at line 102 of file TauValidation.h.
Referenced by beginJob(), and tauDecayChannel().
MonitorElement * TauValidation::TauEta [private] |
Definition at line 102 of file TauValidation.h.
Referenced by analyze(), and beginJob().
double TauValidation::tauEtCut [private] |
Definition at line 93 of file TauValidation.h.
Referenced by rtau().
MonitorElement * TauValidation::TauMothers [private] |
Definition at line 102 of file TauValidation.h.
Referenced by beginJob(), and tauMother().
MonitorElement * TauValidation::TauPhi [private] |
Definition at line 102 of file TauValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement * TauValidation::TauPhotonsN [private] |
Definition at line 102 of file TauValidation.h.
Referenced by beginJob(), and photons().
MonitorElement * TauValidation::TauPhotonsPt [private] |
Definition at line 102 of file TauValidation.h.
Referenced by beginJob(), and photons().
MonitorElement * TauValidation::TauProngs [private] |
Definition at line 102 of file TauValidation.h.
Referenced by beginJob(), and tauProngs().
MonitorElement* TauValidation::TauPt [private] |
Definition at line 102 of file TauValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement * TauValidation::TauRtauHpm [private] |
Definition at line 102 of file TauValidation.h.
Referenced by beginJob(), and rtau().
MonitorElement * TauValidation::TauRtauW [private] |
Definition at line 102 of file TauValidation.h.
Referenced by beginJob(), and rtau().
MonitorElement * TauValidation::TauSpinEffectsHpm [private] |
Definition at line 102 of file TauValidation.h.
Referenced by beginJob(), and spinEffects().
MonitorElement * TauValidation::TauSpinEffectsW [private] |
Definition at line 102 of file TauValidation.h.
Referenced by beginJob(), and spinEffects().
MonitorElement * TauValidation::TauSpinEffectsZ [private] |
Definition at line 102 of file TauValidation.h.
Referenced by beginJob(), and spinEffectsZ().
1.7.3