16 #include <Math/VectorUtil.h> 45 double xfx(
int nset,
double x,
double Q,
int fl);
46 double getXmin(
int nset,
int member);
47 double getXmax(
int nset,
int member);
48 double getQ2min(
int nset,
int member);
49 double getQ2max(
int nset,
int member);
67 <<
" PDF sets requested on input. Using only the first 3 sets and ignoring the rest!!";
74 if (underscore < dot) {
109 int id1 = pdfstuff->
pdf()->
id.first;
110 double x1 = pdfstuff->
pdf()->
x.first;
111 double pdf1 = pdfstuff->
pdf()->
xPDF.first;
113 int id2 = pdfstuff->
pdf()->
id.second;
114 double x2 = pdfstuff->
pdf()->
x.second;
115 double pdf2 = pdfstuff->
pdf()->
xPDF.second;
129 unsigned int gensize = genParticles->size();
131 for (
unsigned int i = 0;
i < gensize; ++
i) {
136 int id = part.
pdgId();
137 if (
id != 23 &&
abs(
id) != 24)
139 mboson = part.
mass();
142 Q =
sqrt(mboson * mboson + Q * Q);
150 std::unique_ptr<std::vector<double> >
weights(
new std::vector<double>);
151 unsigned int nweights = 1;
154 weights->reserve(nweights);
156 for (
unsigned int i = 0;
i < nweights; ++
i) {
160 weights->push_back(newpdf1 / pdf1 * newpdf2 / pdf2);
std::vector< GenParticle > GenParticleCollection
collection of GenParticles
int pdgId() const final
PDG identifier.
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
bool getByToken(EDGetToken token, Handle< PROD > &result) const
std::pair< double, double > x
edm::EDGetTokenT< reco::GenParticleCollection > genToken_
void initPDFSet(int nset, const std::string &filename, int member=0)
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
std::vector< std::string > pdfSetNames_
std::pair< double, double > xPDF
double getXmax(int nset, int member)
double getQ2max(int nset, int member)
#define DEFINE_FWK_MODULE(type)
void extrapolate(bool extrapolate=true)
Abs< T >::type abs(const T &t)
double getXmin(int nset, int member)
PdfWeightProducer(const edm::ParameterSet &)
edm::EDGetTokenT< GenEventInfoProduct > pdfInfoToken_
edm::InputTag pdfInfoTag_
void produce(edm::Event &, const edm::EventSetup &) override
std::vector< std::string > pdfShortNames_
T dot(const Basic3DVector &v) const
Scalar product, or "dot" product, with a vector of same type.
double xfx(int nset, double x, double Q, int fl)
int status() const final
status word
~PdfWeightProducer() override
double getQ2min(int nset, int member)
EDGetTokenT< ProductType > mayConsume(edm::InputTag const &tag)
double mass() const final
mass
void usePDFMember(int nset, int member)