|
|
|
#include <DrellYanValidation.h>
Definition at line 35 of file DrellYanValidation.h.
| DrellYanValidation::DrellYanValidation | ( | const edm::ParameterSet & | iPSet | ) | [explicit] |
Definition at line 20 of file DrellYanValidation.cc.
References dbe, and cmsCodeRules::cppFunctionSkipper::operator.
: _wmanager(iPSet), hepmcCollection_(iPSet.getParameter<edm::InputTag>("hepmcCollection")), _flavor(iPSet.getParameter<int>("decaysTo")), _name(iPSet.getParameter<std::string>("name")) { dbe = 0; dbe = edm::Service<DQMStore>().operator->(); }
| DrellYanValidation::~DrellYanValidation | ( | ) | [virtual] |
Definition at line 30 of file DrellYanValidation.cc.
{}
| void DrellYanValidation::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [virtual] |
Gathering the HepMCProduct information
Implements edm::EDAnalyzer.
Definition at line 77 of file DrellYanValidation.cc.
References _flavor, _wmanager, abs, funct::cos(), cos_theta_gamma_lepton, dilep_mass, dilep_massPeak, dilep_pt, dilep_ptLog, dilep_rap, MonitorElement::Fill(), HepMCValidationHelper::findFSRPhotons(), fPDGTable, gamma_energy, edm::Event::getByLabel(), hepmcCollection_, i, leadeta, leadpt, nEvt, RecoTau_DiTaus_pt_20-420_cfg::ParticleID, perp(), edm::es::products(), seceta, secpt, python::multivaluedict::sort(), HepMCValidationHelper::sortByPt(), lumiQTWidget::t, CommonMethods::weight(), WeightManager::weight(), x, detailsBasic3DVector::y, z, Zdaughters, Zmass, ZmassPeak, Zpt, ZptLog, and Zrap.
{
// we *DO NOT* rely on a Z entry in the particle listings!
edm::Handle<HepMCProduct> evt;
iEvent.getByLabel(hepmcCollection_, evt);
//Get EVENT
const HepMC::GenEvent *myGenEvent = evt->GetEvent();
double weight = _wmanager.weight(iEvent);
//std::cout << "weight: " << weight << std::endl;
nEvt->Fill(0.5,weight);
std::vector<const HepMC::GenParticle*> allproducts;
//requires status 1 for leptons and neutrinos (except tau)
int requiredstatus = (abs(_flavor) == 11 || abs(_flavor) == 12 || abs(_flavor) ==13 || abs(_flavor) ==14 || abs(_flavor) ==16) ? 1 : 3;
bool vetotau = true; //(abs(_flavor) == 11 || abs(_flavor) == 12 || abs(_flavor) ==13 || abs(_flavor) ==14 || abs(_flavor) ==16) ? true : false;
for(HepMC::GenEvent::particle_const_iterator iter = myGenEvent->particles_begin(); iter != myGenEvent->particles_end(); ++iter) {
if (vetotau) {
if ((*iter)->status()==3 && abs((*iter)->pdg_id() == 15) ) return;
}
if((*iter)->status()==requiredstatus) {
if(abs((*iter)->pdg_id())==_flavor)
allproducts.push_back(*iter);
}
}
//nothing to do if we don't have 2 particles
if (allproducts.size() < 2) return;
//sort them in pt
std::sort(allproducts.begin(), allproducts.end(), HepMCValidationHelper::sortByPt);
//get the first element and the first following element with opposite charge
std::vector<const HepMC::GenParticle*> products;
products.push_back(allproducts.front());
const HepPDT::ParticleData* PData1 = fPDGTable->particle(HepPDT::ParticleID(allproducts.front()->pdg_id()));
double charge1 = PData1->charge();
for (unsigned int i = 1; i < allproducts.size(); ++i ){
const HepPDT::ParticleData* PData2 = fPDGTable->particle(HepPDT::ParticleID(allproducts[i]->pdg_id()));
double charge2 = PData2->charge();
if (charge1*charge2 < 0) products.push_back(allproducts[i]);
}
//if we did not find any opposite charge pair there is nothing to do
if (products.size() < 2) return;
assert(products[0]->momentum().perp() >= products[1]->momentum().perp());
//leading lepton with pt > 20.
if (products[0]->momentum().perp() < 20.) return;
//assemble FourMomenta
TLorentzVector lep1(products[0]->momentum().x(), products[0]->momentum().y(), products[0]->momentum().z(), products[0]->momentum().t());
TLorentzVector lep2(products[1]->momentum().x(), products[1]->momentum().y(), products[1]->momentum().z(), products[1]->momentum().t());
TLorentzVector dilepton_mom = lep1 + lep2;
TLorentzVector dilepton_andphoton_mom = dilepton_mom;
//mass > 60.
if (dilepton_mom.M() < 60.) return;
//find possible qed fsr photons
std::vector<const HepMC::GenParticle*> fsrphotons;
HepMCValidationHelper::findFSRPhotons(products, myGenEvent, 0.1, fsrphotons);
Zdaughters->Fill(products[0]->pdg_id(),weight);
Zdaughters->Fill(products[1]->pdg_id(),weight);
std::vector<TLorentzVector> gammasMomenta;
for (unsigned int ipho = 0; ipho < fsrphotons.size(); ++ipho){
TLorentzVector phomom(fsrphotons[ipho]->momentum().x(), fsrphotons[ipho]->momentum().y(), fsrphotons[ipho]->momentum().z(), fsrphotons[ipho]->momentum().t());
dilepton_andphoton_mom += phomom;
Zdaughters->Fill(fsrphotons[ipho]->pdg_id(),weight);
gammasMomenta.push_back(phomom);
}
//Fill Z histograms
Zmass->Fill(dilepton_andphoton_mom.M(),weight);
ZmassPeak->Fill(dilepton_andphoton_mom.M(),weight);
Zpt->Fill(dilepton_andphoton_mom.Pt(),weight);
ZptLog->Fill(log10(dilepton_andphoton_mom.Pt()),weight);
Zrap->Fill(dilepton_andphoton_mom.Rapidity(),weight);
//Fill dilepton histograms
dilep_mass->Fill(dilepton_mom.M(),weight);
dilep_massPeak->Fill(dilepton_mom.M(),weight);
dilep_pt->Fill(dilepton_mom.Pt(),weight);
dilep_ptLog->Fill(log10(dilepton_mom.Pt()),weight);
dilep_rap->Fill(dilepton_mom.Rapidity(),weight);
//Fill lepton histograms
leadpt->Fill(lep1.Pt(),weight);
secpt->Fill(lep2.Pt(),weight);
leadeta->Fill(lep1.Eta(),weight);
seceta->Fill(lep2.Eta(),weight);
//boost everything in the Z frame
TVector3 boost = dilepton_andphoton_mom.BoostVector();
boost*=-1.;
lep1.Boost(boost);
lep2.Boost(boost);
for (unsigned int ipho = 0; ipho < gammasMomenta.size(); ++ipho){
gammasMomenta[ipho].Boost(boost);
}
std::sort(gammasMomenta.begin(), gammasMomenta.end(), HepMCValidationHelper::GreaterByE<TLorentzVector>);
//fill gamma histograms
if (gammasMomenta.size() != 0 && dilepton_andphoton_mom.M() > 50.) {
gamma_energy->Fill(gammasMomenta.front().E(),weight);
double dphi = lep1.DeltaR(gammasMomenta.front()) < lep2.DeltaR(gammasMomenta.front()) ?
lep1.DeltaPhi(gammasMomenta.front()) : lep2.DeltaPhi(gammasMomenta.front());
cos_theta_gamma_lepton->Fill(cos(dphi),weight);
}
}//analyze
| void DrellYanValidation::beginJob | ( | void | ) | [virtual] |
Setting the DQM top directories
Reimplemented from edm::EDAnalyzer.
Definition at line 32 of file DrellYanValidation.cc.
References _name, DQMStore::book1D(), cos_theta_gamma_lepton, dbe, dilep_mass, dilep_massPeak, dilep_pt, dilep_ptLog, dilep_rap, gamma_energy, leadeta, leadpt, nEvt, seceta, secpt, DQMStore::setCurrentFolder(), Zdaughters, Zmass, ZmassPeak, Zpt, ZptLog, and Zrap.
{
if(dbe){
std::string folderName = "Generator/DrellYan";
folderName+=_name;
dbe->setCurrentFolder(folderName.c_str());
// Number of analyzed events
nEvt = dbe->book1D("nEvt", "n analyzed Events", 1, 0., 1.);
//Kinematics
Zmass = dbe->book1D("Zmass","inv. Mass Z", 70 ,0,140);
ZmassPeak = dbe->book1D("ZmassPeak","inv. Mass Z", 80 ,80 ,100);
Zpt = dbe->book1D("Zpt","Z pt",100,0,200);
ZptLog = dbe->book1D("ZptLog","log(Z pt)",100,0.,5.);
Zrap = dbe->book1D("Zrap", "Z y", 100, -5, 5);
Zdaughters = dbe->book1D("Zdaughters", "Z daughters", 60, -30, 30);
dilep_mass = dbe->book1D("dilep_mass","inv. Mass dilepton", 70 ,0,140);
dilep_massPeak = dbe->book1D("dilep_massPeak","inv. Mass dilepton", 80 ,80 ,100);
dilep_pt = dbe->book1D("dilep_pt","dilepton pt",100,0,200);
dilep_ptLog = dbe->book1D("dilep_ptLog","log(dilepton pt)",100,0.,5.);
dilep_rap = dbe->book1D("dilep_rap", "dilepton y", 100, -5, 5);
gamma_energy = dbe->book1D("gamma_energy", "photon energy in Z rest frame", 200, 0., 100.);
cos_theta_gamma_lepton = dbe->book1D("cos_theta_gamma_lepton", "cos_theta_gamma_lepton in Z rest frame", 200, -1, 1);
leadpt = dbe->book1D("leadpt","leading lepton pt", 200, 0., 200.);
secpt = dbe->book1D("secpt","second lepton pt", 200, 0., 200.);
leadeta = dbe->book1D("leadeta","leading lepton eta", 100, -5., 5.);
seceta = dbe->book1D("seceta","second lepton eta", 100, -5., 5.);
}
return;
}
| void DrellYanValidation::beginRun | ( | const edm::Run & | iRun, |
| const edm::EventSetup & | iSetup | ||
| ) | [virtual] |
Get PDT Table
Reimplemented from edm::EDAnalyzer.
Definition at line 70 of file DrellYanValidation.cc.
References fPDGTable, and edm::EventSetup::getData().
| void DrellYanValidation::endJob | ( | void | ) | [virtual] |
| void DrellYanValidation::endRun | ( | const edm::Run & | iRun, |
| const edm::EventSetup & | iSetup | ||
| ) | [virtual] |
int DrellYanValidation::_flavor [private] |
std::string DrellYanValidation::_name [private] |
WeightManager DrellYanValidation::_wmanager [private] |
Definition at line 48 of file DrellYanValidation.h.
Referenced by analyze().
Definition at line 62 of file DrellYanValidation.h.
Referenced by analyze(), and beginJob().
DQMStore* DrellYanValidation::dbe [private] |
ME's "container".
Definition at line 56 of file DrellYanValidation.h.
Referenced by beginJob(), and DrellYanValidation().
MonitorElement* DrellYanValidation::dilep_mass [private] |
Definition at line 60 of file DrellYanValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement * DrellYanValidation::dilep_massPeak [private] |
Definition at line 60 of file DrellYanValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement * DrellYanValidation::dilep_pt [private] |
Definition at line 60 of file DrellYanValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement * DrellYanValidation::dilep_ptLog [private] |
Definition at line 60 of file DrellYanValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement * DrellYanValidation::dilep_rap [private] |
Definition at line 60 of file DrellYanValidation.h.
Referenced by analyze(), and beginJob().
PDT table.
Definition at line 53 of file DrellYanValidation.h.
Referenced by analyze(), and beginRun().
MonitorElement* DrellYanValidation::gamma_energy [private] |
Definition at line 62 of file DrellYanValidation.h.
Referenced by analyze(), and beginJob().
Definition at line 50 of file DrellYanValidation.h.
Referenced by analyze().
MonitorElement * DrellYanValidation::leadeta [private] |
Definition at line 61 of file DrellYanValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement* DrellYanValidation::leadpt [private] |
Definition at line 61 of file DrellYanValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement* DrellYanValidation::nEvt [private] |
Definition at line 58 of file DrellYanValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement * DrellYanValidation::seceta [private] |
Definition at line 61 of file DrellYanValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement * DrellYanValidation::secpt [private] |
Definition at line 61 of file DrellYanValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement * DrellYanValidation::Zdaughters [private] |
Definition at line 59 of file DrellYanValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement* DrellYanValidation::Zmass [private] |
Definition at line 59 of file DrellYanValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement * DrellYanValidation::ZmassPeak [private] |
Definition at line 59 of file DrellYanValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement * DrellYanValidation::Zpt [private] |
Definition at line 59 of file DrellYanValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement * DrellYanValidation::ZptLog [private] |
Definition at line 59 of file DrellYanValidation.h.
Referenced by analyze(), and beginJob().
MonitorElement * DrellYanValidation::Zrap [private] |
Definition at line 59 of file DrellYanValidation.h.
Referenced by analyze(), and beginJob().
1.7.3