[u/mrichter/AliRoot.git] / TRD / AliTRDsimpleGen.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/
 
/*
$Log$                                                          
*/
 
///////////////////////////////////////////////////////////////////////////////
//                                                                           //
//  Particle generator for the simplescopic TRD simulator                    //
//                                                                           //
///////////////////////////////////////////////////////////////////////////////
 
#include <stdlib.h>
 
#include <TRandom.h>

#include "AliTRDsimpleGen.h"
#include "AliTRDsimpleMC.h"
 
ClassImp(AliTRDsimpleGen)
 
//_____________________________________________________________________________
AliTRDsimpleGen::AliTRDsimpleGen():TObject()
{                       
  //
  // AliTRDsimpleGen default constructor
  //

  fPdg    = 211;
  fMomMin = 1.0;
  fMomMax = 1.0;             
                                            
}                                                                               
 
//_____________________________________________________________________________
AliTRDsimpleGen::AliTRDsimpleGen(const AliTRDsimpleGen &g)
{
  //
  // AliTRDsimpleGen copy constructor
  //
 
  ((AliTRDsimpleGen &) g).Copy(*this);
 
}
 
//_____________________________________________________________________________
AliTRDsimpleGen::~AliTRDsimpleGen()
{
  //
  // AliTRDsimpleGen destructor
  //
 
}                                                                               
 
//_____________________________________________________________________________
AliTRDsimpleGen &AliTRDsimpleGen::operator=(const AliTRDsimpleGen &g)
{
  //
  // Assignment operator
  //
 
  if (this != &g) ((AliTRDsimpleGen &) g).Copy(*this);
  return *this;
 
}
 
//_____________________________________________________________________________
void AliTRDsimpleGen::Copy(TObject &g)
{
  //
  // Copy function
  //                             

  ((AliTRDsimpleGen &) g).fPdg     = fPdg;                                                 
  ((AliTRDsimpleGen &) g).fMomMin  = fMomMin;                                                 
  ((AliTRDsimpleGen &) g).fMomMax  = fMomMax;                                                 

}

//_____________________________________________________________________________
void AliTRDsimpleGen::NewParticle(Int_t ievent)
{
  //
  // Generate a new particle and initialize the MC object
  // 

  if (ievent == 0) {
    printf("\n");
    printf("<AliTRDsimpleGen> Generate particles with PDG code %d\n",fPdg);
    if (fMomMax > fMomMin) {
      printf("<AliTRDsimpleGen> Momentum range = %4.2f - %4.2f GeV/c\n"
            ,fMomMin,fMomMax);
    }
    else {
      printf("<AliTRDsimpleGen> Fixed momentum = %4.2f GeV/c\n"
            ,fMomMax);
    }
    printf("\n");
  }

  Double_t p = fMomMax;
  if (fMomMax > fMomMin) {
    p = (fMomMax - fMomMin) * gRandom->Rndm() + fMomMin;
  }

  Double_t px = p;
  Double_t py = 0.0;      
  Double_t pz = 0.0;      

  ((AliTRDsimpleMC *) gMC)->NewTrack(ievent,fPdg,px,py,pz);

}
Commit	Line	Data
16bf9884	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
	16	/*
	17	$Log$
	18	*/
	19
	20	///////////////////////////////////////////////////////////////////////////////
	21	// //
	22	// Particle generator for the simplescopic TRD simulator //
	23	// //
	24	///////////////////////////////////////////////////////////////////////////////
	25
	26	#include <stdlib.h>
	27
	28	#include <TRandom.h>
	29
	30	#include "AliTRDsimpleGen.h"
	31	#include "AliTRDsimpleMC.h"
	32
	33	ClassImp(AliTRDsimpleGen)
	34
	35	//_____________________________________________________________________________
	36	AliTRDsimpleGen::AliTRDsimpleGen():TObject()
	37	{
	38	//
	39	// AliTRDsimpleGen default constructor
	40	//
	41
	42	fPdg = 211;
	43	fMomMin = 1.0;
	44	fMomMax = 1.0;
	45
	46	}
	47
	48	//_____________________________________________________________________________
	49	AliTRDsimpleGen::AliTRDsimpleGen(const AliTRDsimpleGen &g)
	50	{
	51	//
	52	// AliTRDsimpleGen copy constructor
	53	//
	54
	55	((AliTRDsimpleGen &) g).Copy(*this);
	56
	57	}
	58
	59	//_____________________________________________________________________________
	60	AliTRDsimpleGen::~AliTRDsimpleGen()
	61	{
	62	//
	63	// AliTRDsimpleGen destructor
	64	//
65
66	}
67
68	//_____________________________________________________________________________
69	AliTRDsimpleGen &AliTRDsimpleGen::operator=(const AliTRDsimpleGen &g)
70	{
71	//
72	// Assignment operator
73	//
74
75	if (this != &g) ((AliTRDsimpleGen &) g).Copy(*this);
76	return *this;
77
78	}
79
80	//_____________________________________________________________________________
81	void AliTRDsimpleGen::Copy(TObject &g)
82	{
83	//
84	// Copy function
85	//
86
87	((AliTRDsimpleGen &) g).fPdg = fPdg;
88	((AliTRDsimpleGen &) g).fMomMin = fMomMin;
89	((AliTRDsimpleGen &) g).fMomMax = fMomMax;
90
91	}
92
93	//_____________________________________________________________________________
94	void AliTRDsimpleGen::NewParticle(Int_t ievent)
95	{
96	//
97	// Generate a new particle and initialize the MC object
98	//
99
100	if (ievent == 0) {
101	printf("\n");
102	printf("<AliTRDsimpleGen> Generate particles with PDG code %d\n",fPdg);
103	if (fMomMax > fMomMin) {
104	printf("<AliTRDsimpleGen> Momentum range = %4.2f - %4.2f GeV/c\n"
105	,fMomMin,fMomMax);
106	}
107	else {
108	printf("<AliTRDsimpleGen> Fixed momentum = %4.2f GeV/c\n"
109	,fMomMax);
110	}
111	printf("\n");
112	}
113
114	Double_t p = fMomMax;
115	if (fMomMax > fMomMin) {
116	p = (fMomMax - fMomMin) * gRandom->Rndm() + fMomMin;
117	}
118
119	Double_t px = p;
120	Double_t py = 0.0;
121	Double_t pz = 0.0;
122
123	((AliTRDsimpleMC *) gMC)->NewTrack(ievent,fPdg,px,py,pz);
124
125	}