Page MenuHomeHEPForge
Files F8309676

CrossSection.h
No OneTemporary
Actions

CrossSection.h
View Options

//==============================================================================
// CrossSection.h
//
// Copyright (C) 2010-2017 Tobias Toll and Thomas Ullrich
//
// This file is part of Sartre.
//
// This program is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation.
// This program is distributed in the hope that it will be useful,
// but without any warranty; without even the implied warranty of
// merchantability or fitness for a particular purpose. See the
// GNU General Public License for more details.
// You should have received a copy of the GNU General Public License
// along with this program. If not, see <http://www.gnu.org/licenses/>.
//
// Author: Thomas Ullrich
// Last update:
// $Date: 2018-08-15 13:19:49 +0100 (Wed, 15 Aug 2018) $
// $Author: ttoll $
//==============================================================================
//
// Functor class.
//
// operator() returns d^3sig/(dt dQ2 dW2) in nb/GeV^6.
//
//===============================================================================
#ifndef CrossSection_h
#define CrossSection_h
#include "Enumerations.h"
#include "PhotonFlux.h"
class TRandom3;
class EventGeneratorSettings;
class TableCollection;
class CrossSection {
public:
CrossSection(TableCollection* = 0, TableCollection* = 0);
~CrossSection();
double operator()(double t, double Q2, double W2);
double operator()(double t, double xpom); // UPC version
double operator()(const double*); // array of t, Q2, W2 or t, xpom for UPC
double unuranPDF(const double*); // for UNU.RAN using log(Q2) (or log(xpom) for UPC)
// and returning log of cross-section
void setTableCollection(TableCollection*);
void setProtonTableCollection(TableCollection*);
GammaPolarization polarizationOfLastCall() const;
DiffractiveMode diffractiveModeOfLastCall() const;
void setCheckKinematics(bool);
protected:
friend class Sartre; // mostly for debugging and QA
double dsigdtdQ2dW2_total_checked(double t, double Q2, double W2);
double dsigdtdxp_total_checked(double t, double xpom);
double dsigdt_total(double t, double Q2, double W2, GammaPolarization) const; // modified
double dsigdt_coherent(double t, double Q2, double W2, GammaPolarization) const;
double dsigdt_incoherent(double t, double Q2, double W2, GammaPolarization) const; // new
double dsigdtdQ2dW2_total(double t, double Q2, double W2, GammaPolarization) const;
double dsigdtdQ2dW2_coherent(double t, double Q2, double W2, GammaPolarization) const;
double dsigdt_total(double t, double xpom) const; // UPC version
double dsigdt_coherent(double t, double xpom) const; // UPC version
double dsigdt_incoherent(double t, double xpom) const; // UPC version
double dsigdtdxp_total(double t, double xpom) const; // UPC
double dsigdtdxp_coherent(double t, double xpom) const; // UPC
double logDerivateOfAmplitude(double t, double Q2, double W2, GammaPolarization) const;
double logDerivateOfGluonDensity(double t, double Q2, double W2, GammaPolarization) const;
double realAmplitudeCorrection(double t, double Q2, double W2, GammaPolarization pol) const;
double skewednessCorrection(double t, double Q2, double W2, GammaPolarization pol) const;
double logDerivateOfAmplitude(double t, double xpom) const; // UPC version
double logDerivateOfGluonDensity(double t, double xpom) const; // UPC version
double realAmplitudeCorrection(double t, double xpom) const; // UPC version
double skewednessCorrection(double t, double xpom) const; // UPC version
double skewednessCorrection(double lambda) const;
double realAmplitudeCorrection(double lambda) const;
double UPCPhotonFlux(double t, double xpom) const;
private:
TRandom3* mRandom;
GammaPolarization mPolarization;
DiffractiveMode mDiffractiveMode;
PhotonFlux mPhotonFlux;
TableCollection* mTableCollection;
TableCollection* mProtonTableCollection;
EventGeneratorSettings* mSettings;
double mS;
double mVmMass;
bool mCheckKinematics;
};
#endif

File Metadata

Mime Type
text/x-c++
Expires
Sat, Dec 21, 4:10 PM (21 h, 30 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4023394
Default Alt Text
CrossSection.h (4 KB)

Event Timeline