Page MenuHomeHEPForge
Files F19244700

DrellYanPT.cc
No OneTemporary
Actions

Size
2 KB
Referenced Files
None
Subscribers
None

DrellYanPT.cc
View Options

// -*- C++ -*-
//
// This is the implementation of the non-inlined, non-templated member
// functions of the DrellYanPT class.
//
#include "DrellYanPT.h"
#include "ThePEG/Interface/ClassDocumentation.h"
#include "ThePEG/Repository/EventGenerator.h"
#include "ThePEG/EventRecord/Particle.h"
#include "ThePEG/EventRecord/Event.h"
#include "ThePEG/PDT/EnumParticles.h"
#include "ThePEG/Persistency/PersistentOStream.h"
#include "ThePEG/Persistency/PersistentIStream.h"
using namespace Herwig;
DrellYanPT::DrellYanPT()
: _Zpt(0.,250.,250), _Wppt(0.,250.,250), _Wmpt(0.,250.,250) {}
void DrellYanPT::dofinish() {
AnalysisHandler::dofinish();
ofstream outZ ("pt_Z.dat");
_Zpt.normaliseToCrossSection();
_Zpt.simpleOutput(outZ,true);
ofstream outWm ("pt_Wm.dat");
_Wmpt.normaliseToCrossSection();
_Wmpt.simpleOutput(outWm,true);
ofstream outWp ("pt_Wp.dat");
_Wppt.normaliseToCrossSection();
_Wppt.simpleOutput(outWp,true);
}
void DrellYanPT::analyze(tEventPtr event, long ieve, int loop, int state) {
AnalysisHandler::analyze(event, ieve, loop, state);
// Rotate to CMS, extract final state particles and call analyze(particles).
StepVector::const_iterator sit =event->primaryCollision()->steps().begin();
StepVector::const_iterator send=event->primaryCollision()->steps().end();
for(;sit!=send;++sit) {
ParticleSet part=(**sit).all();
ParticleSet::const_iterator iter=part.begin();
ParticleSet::const_iterator end =part.end();
for( ;iter!=end;++iter) {
if(((**iter).id()==ParticleID::Z0||(**iter).id()==ParticleID::gamma)
&& (**iter).children().size()==2) {
_Zpt.addWeighted((**iter).momentum().perp()/GeV,event->weight());
} else if ((**iter).id()==ParticleID::Wplus && (**iter).children().size()==2) {
_Wppt.addWeighted((**iter).momentum().perp()/GeV,event->weight());
} else if ((**iter).id()==ParticleID::Wminus && (**iter).children().size()==2) {
_Wmpt.addWeighted((**iter).momentum().perp()/GeV,event->weight());
}
}
}
}
NoPIOClassDescription<DrellYanPT> DrellYanPT::initDrellYanPT;
// Definition of the static class description member.
void DrellYanPT::Init() {
static ClassDocumentation<DrellYanPT> documentation
("Analyses the pt of weak bosons produces in Drell-Yan processes.");
}

File Metadata

Mime Type
text/x-c
Expires
Tue, Sep 30, 4:44 AM (1 d, 7 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
6357167
Default Alt Text
DrellYanPT.cc (2 KB)

Event Timeline

HEPForge.org · .