PndPiPiGenerator.h

Go to the documentation of this file.

//****************************************************************************
//*                   This file is part of PandaRoot.                        *
//*                                                                          *
//*            PandaRoot is distributed under the terms of the               *
//*              GNU General Public License (GPL) version 3,                 *
//*                 copied verbatim in the file "LICENSE".                   *
//*                                                                          *
//*  Copyright (C) 2006 - 2024 FAIR GmbH and copyright holders of PandaRoot  *
//*     The copyright holders are listed in the file "COPYRIGHTHOLDERS".     *
//*               The authors are listed in the file "AUTHORS".              *
//****************************************************************************

//----------------------------------------------------//
//            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