PndPiPiGenerator.h

Go to the documentation of this file.

//----------------------------------------------------//
//            PndPiPiGenerator header file            //
//----------------------------------------------------//

#ifndef PND_PIPIGENERATOR_H
#define PND_PIPIGENERATOR_H

// c/c++
#include <stdlib.h>
#include <iostream>
#include <fstream>
#include <sstream>
#include <cmath>
#include <algorithm>
#include <time.h>

//ranlux
#include "ranlxs.h"
#include "ranlxd.h"

void ranlxs(float r[], int n);
void ranlxd(double r[], int n);

void rlxs_init(int level, int seed);
void rlxd_init(int level, int seed);

//mz (general)
#include "mzparameters.h"
#include "mzfunctions.h"

//mz for pp_to_pipi_vandewi
#include "mzfunctions_pp_to_pipi_vandewi.h"

#include "FairGenerator.h"

using namespace std;

class FairPrimaryGenerator;

class PndPiPiGenerator : public FairGenerator
{
public:
  PndPiPiGenerator();

  virtual ~PndPiPiGenerator();

  void SetBeamMom(Double_t in_P)                 {P = in_P;};                           // set pbar momentum
  void SetSeed(Int_t in_seed)                    {seed = in_seed;};                     // set seed
  void SetCosThetaMin(Double_t in_cos_theta_min) {cos_theta_min = in_cos_theta_min;};   // set cos(theta*) min in the pbarp CM frame
  void SetCosThetaMax(Double_t in_cos_theta_max) {cos_theta_max = in_cos_theta_max;};   // set cos(theta*) max in the pbarp CM frame

  virtual Bool_t ReadEvent(FairPrimaryGenerator* primGen); 

  Bool_t Init();

  ClassDef(PndPiPiGenerator,1);

  Double_t P;               // pbar momentum
  Int_t seed;               // seed
  Double_t cos_theta_min;   // cos(theta*) min in the pbarp CM frame
  Double_t cos_theta_max;   // cos(theta*) max in the pbarp CM frame
  Double_t piplus_p[4];     // lab frame: 4-mom pi+
  Double_t piminus_p[4];    // lab frame: 4-mom pi-
};

#endif