The information about a hadronic interaction needed to calculate weights. More...

#include <InteractionData.h>

Public Member Functions
	InteractionData ()
	Default Constructor. More...

	InteractionData (int genid, double incMom[], int incPdg, double prodMom[], int prodPdg, std::string volname, std::string procname, double vtx[])
	Constructor given kinematic of the interaction. More...

virtual	~InteractionData ()

std::ostream &	print (std::ostream &os) const

Public Attributes
int	gen
	generation index More...

int	Inc_pdg
	pdg code of the incident particle More...

int	Prod_pdg
	pdg code of the produced particle More...

double	Inc_P
	Momentum magnitude of the incident particle. More...

double	Inc_P4 [4]
	Momentum 4 vector of the incident particle, E=p[3]. More...

double	Prod_P
	Momentum magnitude of the produced particle. More...

double	Prod_P4 [4]
	Momentum 4 vector of the produced particle, E=p[3]. More...

double	Vtx [3]
	Location of the interaction. More...

double	Inc_Mass
	Mass of the incident particle. More...

double	Prod_Mass
	Mass of the produced particle. More...

double	xF
	Feynmann-x of the produced particle: $x_{F} = 2P_{L}/\sqrt(s)$ . More...

double	Pz
	Longitudinal momentum (GeV/c) of the produced particle. More...

double	Theta
	Angle (rad) of the produced particle. More...

double	Pt
	Transversal momentum (GeV/c) of the produced particle. More...

double	Ecm
	Center of mass energy of the collision indident particle - nuclear proton. More...

double	Betacm
	$\beta_{CM}$ More...

double	Gammacm
	$\gamma_{CM}$ More...

std::string	Vol
	Interaction volume. More...

std::string	Proc
	Interaction process. More...

Private Attributes
TDatabasePDG *	particle

Detailed Description

The information about a hadronic interaction needed to calculate weights.

Information about the kinematics of the interaction, the identity of the projectile and target, the producted particles, etc. This looks like KinProd to me.

Definition at line 19 of file InteractionData.h.

Constructor & Destructor Documentation

NeutrinoFluxReweight::InteractionData::InteractionData ( )

Default Constructor.

Definition at line 8 of file InteractionData.cpp.

                                   {
     
     particle = TDatabasePDG::Instance();
     
     gen     = 0;
     Inc_pdg = 0;
     Prod_pdg= 0;
     
     Inc_P     = -1000.;
     Prod_P    = -1000.;
     Inc_Mass  = -1000.;
     Prod_Mass = -1000.;
 
     for(int i=0; i<4; i++){
       Inc_P4[i]=0;
       Prod_P4[i]=0;
       if(i<3) Vtx[i]=0;
     }
         
     xF      = -1000.;
     Pz      = -1000.;
     Theta   = -1000.;
     Pt      = -1000.;
     
     Ecm     = -1000.;
     Betacm  = -1000.;
     Gammacm = -1000.;
     
     Vol = "NoDefinied";
     
   }

NeutrinoFluxReweight::InteractionData::InteractionData	(	int	genid,
		double	incMom[],
		int	incPdg,
		double	prodMom[],
		int	prodPdg,
		std::string	volname,
		std::string	procname,
		double	vtx[]
	)

Constructor given kinematic of the interaction.

Definition at line 40 of file InteractionData.cpp.

                                                                                                                                                             {
 
     particle = TDatabasePDG::Instance();
     // Z direction along the direction of the incident particle
     // Cos between incMom and prodMom:
     // The units are in GeV
     
     InteractionData::gen      = genid;
     
     InteractionData::Inc_pdg  = incPdg;
     InteractionData::Prod_pdg = prodPdg;
     
     InteractionData::Inc_P  = std::sqrt(incMom[0]*incMom[0] + incMom[1]*incMom[1] +incMom[2]*incMom[2]);
     InteractionData::Prod_P = std::sqrt(prodMom[0]*prodMom[0] + prodMom[1]*prodMom[1] +prodMom[2]*prodMom[2]);
     
     double cos_theta = (incMom[0]*prodMom[0]+incMom[1]*prodMom[1]+incMom[2]*prodMom[2])/(Inc_P*Prod_P);
     double sin_theta = std::sqrt(1.-pow(cos_theta,2.0));
     
     //Theta in rads:  
     InteractionData::Theta = std::acos(cos_theta);
     
     InteractionData::Pt = Prod_P*sin_theta;
     InteractionData::Pz = Prod_P*cos_theta; 
     
     if(Inc_pdg != 1000010020)Inc_Mass = particle->GetParticle(Inc_pdg)->Mass();
     else {Inc_Mass = 1.875;}
     
     if(Prod_pdg != 1000010020)Prod_Mass = particle->GetParticle(Prod_pdg)->Mass();
     else{Prod_Mass = 1.875;}
      
     InteractionData::Inc_Mass  = particle->GetParticle(Inc_pdg)->Mass();
     InteractionData::Prod_Mass = particle->GetParticle(Prod_pdg)->Mass();
 
 
     //Ecm, gamma:
     double inc_E_lab = std::sqrt(Inc_P*Inc_P + pow(Inc_Mass,2));
     InteractionData::Ecm       = std::sqrt(2.*pow(Inc_Mass,2)+2.*inc_E_lab*Inc_Mass); 
     InteractionData::Betacm    = std::sqrt(pow(inc_E_lab,2)-pow(Inc_Mass,2.0))/(inc_E_lab + Inc_Mass); 
     InteractionData::Gammacm   = 1./std::sqrt(1.-pow(Betacm,2.0));
     
     //xF:
     double prod_E_lab  = std::sqrt(Prod_P*Prod_P + pow(Prod_Mass,2));
     double PL          = Gammacm*(Pz-Betacm*prod_E_lab);  // PL is measured in CM frame
     InteractionData::xF  = PL*2./Ecm;
 
     //4 momenta:
     Inc_P4[3]=inc_E_lab;
     Prod_P4[3]=prod_E_lab;
     for(int i=0; i<3; i++) {Inc_P4[i]=incMom[i]; Prod_P4[i]=prodMom[i];}
 
     
     //Volume:
     InteractionData::Vol = volname;
 
     //Process:
     InteractionData::Proc = procname;
 
     //Vertex:
     for(int i=0; i<3; i++) Vtx[i]=vtx[i];
 
   }

NeutrinoFluxReweight::InteractionData::~InteractionData ( )

virtual

Definition at line 103 of file InteractionData.cpp.

                                    {
     
   }

Member Function Documentation

std::ostream & NeutrinoFluxReweight::InteractionData::print ( std::ostream & os ) const

Definition at line 107 of file InteractionData.cpp.

                                                          {
     using namespace std;
     os<<"in:"<<setw(5)<<Inc_pdg
       <<"|p3:";
     for(int i=0; i<3; i++) {
       os<<setiosflags(ios::fixed) << setprecision(2)<<setw(6)<<Inc_P4[i]<<" ";
     }
     os<<"||out:"<<setw(5)<<Prod_pdg
       <<"|p3:"<<setiosflags(ios::fixed) << setprecision(2);
     for(int i=0; i<3; i++) {
       os<<setiosflags(ios::fixed) << setprecision(2)<<setw(6)<<Prod_P4[i]<<" ";
     }
     os <<"|v3:";
     for(int i=0; i<3; i++) {
       os<<setiosflags(ios::fixed) << setprecision(2)<<setw(5)<<Vtx[i]<<" ";
     }
     os<<"xF,pT:"<<xF<<","<<Pt;
     os<<endl;
     return os;
   }

Member Data Documentation

double NeutrinoFluxReweight::InteractionData::Betacm

$\beta_{CM}$

Definition at line 77 of file InteractionData.h.

double NeutrinoFluxReweight::InteractionData::Ecm

Center of mass energy of the collision indident particle - nuclear proton.

Definition at line 74 of file InteractionData.h.

double NeutrinoFluxReweight::InteractionData::Gammacm

$\gamma_{CM}$

Definition at line 80 of file InteractionData.h.

int NeutrinoFluxReweight::InteractionData::gen

generation index

Definition at line 32 of file InteractionData.h.

double NeutrinoFluxReweight::InteractionData::Inc_Mass

Mass of the incident particle.

Definition at line 56 of file InteractionData.h.

double NeutrinoFluxReweight::InteractionData::Inc_P

Momentum magnitude of the incident particle.

Definition at line 41 of file InteractionData.h.

double NeutrinoFluxReweight::InteractionData::Inc_P4[4]

Momentum 4 vector of the incident particle, E=p[3].

Definition at line 44 of file InteractionData.h.

int NeutrinoFluxReweight::InteractionData::Inc_pdg

pdg code of the incident particle

Definition at line 35 of file InteractionData.h.

TDatabasePDG* NeutrinoFluxReweight::InteractionData::particle

private

Definition at line 91 of file InteractionData.h.

std::string NeutrinoFluxReweight::InteractionData::Proc

Interaction process.

Definition at line 86 of file InteractionData.h.

double NeutrinoFluxReweight::InteractionData::Prod_Mass

Mass of the produced particle.

Definition at line 59 of file InteractionData.h.

double NeutrinoFluxReweight::InteractionData::Prod_P

Momentum magnitude of the produced particle.

Definition at line 47 of file InteractionData.h.

double NeutrinoFluxReweight::InteractionData::Prod_P4[4]

Momentum 4 vector of the produced particle, E=p[3].

Definition at line 50 of file InteractionData.h.

int NeutrinoFluxReweight::InteractionData::Prod_pdg

pdg code of the produced particle

Definition at line 38 of file InteractionData.h.

double NeutrinoFluxReweight::InteractionData::Pt

Transversal momentum (GeV/c) of the produced particle.

Definition at line 71 of file InteractionData.h.

double NeutrinoFluxReweight::InteractionData::Pz

Longitudinal momentum (GeV/c) of the produced particle.

Definition at line 65 of file InteractionData.h.

double NeutrinoFluxReweight::InteractionData::Theta

Angle (rad) of the produced particle.

Definition at line 68 of file InteractionData.h.

std::string NeutrinoFluxReweight::InteractionData::Vol

Interaction volume.

Definition at line 83 of file InteractionData.h.

double NeutrinoFluxReweight::InteractionData::Vtx[3]

Location of the interaction.

Definition at line 53 of file InteractionData.h.

double NeutrinoFluxReweight::InteractionData::xF

Feynmann-x of the produced particle: $x_{F} = 2P_{L}/\sqrt(s)$ .

Definition at line 62 of file InteractionData.h.

The documentation for this class was generated from the following files:

ppfx/include/InteractionData.h
ppfx/src/InteractionData.cpp

Public Member Functions

Public Attributes

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation