| 1 |
// $Id: PamVMCPrimaryGenerator.h,v 1.0 2007/06/03 |
| 2 |
|
| 3 |
#ifndef PAMVMC_PRIMARY_GENERATOR_H |
| 4 |
#define PAMVMC_PRIMARY_GENERATOR_H |
| 5 |
|
| 6 |
#include <iostream> |
| 7 |
#include <TVirtualMCApplication.h> |
| 8 |
#include <TClonesArray.h> |
| 9 |
#include <TVector3.h> |
| 10 |
#include <TVector2.h> |
| 11 |
#include <TMath.h> |
| 12 |
#include <TRandom3.h> |
| 13 |
|
| 14 |
#include "PamRootManager.h" |
| 15 |
|
| 16 |
class TVirtualMCStack; |
| 17 |
|
| 18 |
using std::cout; |
| 19 |
using std::endl; |
| 20 |
using TMath::Sqrt; |
| 21 |
|
| 22 |
class PamVMCPrimary : public TObject |
| 23 |
{ |
| 24 |
|
| 25 |
public: |
| 26 |
|
| 27 |
PamVMCPrimary(): fPDG(0), fX0(0.), fY0(0.), fZ0(0.), fTHETA(0.), fPHI(0.), fP0(0.), fGOOD(kFALSE) { }; |
| 28 |
|
| 29 |
|
| 30 |
Int_t fPDG; |
| 31 |
Double_t fX0, fY0, fZ0; |
| 32 |
Double_t fTHETA, fPHI; // RADIANTS |
| 33 |
Double_t fP0; |
| 34 |
Bool_t fGOOD; |
| 35 |
|
| 36 |
void Clean() { fPDG=0; fX0=fY0=fZ0=fTHETA=fPHI=fP0=0.; fGOOD=kFALSE;} |
| 37 |
|
| 38 |
void Clear(const Option_t * =""){ Clean(); } |
| 39 |
|
| 40 |
void Print(const Option_t * ="") const |
| 41 |
{ cout<<"PRIMARY particle infromation:"<<endl; |
| 42 |
cout<<"Pdg="<<fPDG<<endl; //add name later |
| 43 |
cout<<"Position: "<<"("<<fX0<<","<<fY0<<","<<fZ0<<")"<<endl; |
| 44 |
cout<<"P0, Theta, Phi: "<<fP0<<","<<fTHETA*180./3.1415<<","<<fPHI*180./3.1415<<endl; |
| 45 |
cout<<"GOOD Single Track:"<<fGOOD<<endl; |
| 46 |
} |
| 47 |
|
| 48 |
ClassDef(PamVMCPrimary,1); |
| 49 |
}; |
| 50 |
|
| 51 |
PamVMCPrimary & operator+=(PamVMCPrimary &a, const PamVMCPrimary &b); |
| 52 |
|
| 53 |
class PamVMCPrimaryGenerator : public TObject |
| 54 |
{ |
| 55 |
public: |
| 56 |
PamVMCPrimaryGenerator(TVirtualMCStack* stack, UInt_t seed); |
| 57 |
//PamVMCPrimaryGenerator(UInt_t seed); |
| 58 |
PamVMCPrimaryGenerator(); |
| 59 |
|
| 60 |
virtual ~PamVMCPrimaryGenerator(); |
| 61 |
|
| 62 |
// methods |
| 63 |
|
| 64 |
void GeneratePrimary(); |
| 65 |
// set methods |
| 66 |
void SetParticle(Int_t pdg); |
| 67 |
void SetPosition(Double_t vx, Double_t vy, Double_t vz) |
| 68 |
{fprim.fX0=vx; fprim.fY0=vy; fprim.fZ0=vz;}; |
| 69 |
|
| 70 |
void SetKinEnergy(Double_t kinEnergy) {fprim.fP0=KinEToMomentum(kinEnergy);}; |
| 71 |
void SetRigidity(Double_t rigidity) {fprim.fP0=RigToMomentum(rigidity);}; |
| 72 |
void SetMomentum(Double_t momentum) {fprim.fP0=momentum; }; |
| 73 |
|
| 74 |
void SetDirection(Double_t theta, Double_t phi) {fprim.fTHETA=theta; fprim.fPHI=phi; }; |
| 75 |
void SetMomentum(Double_t px, Double_t py, Double_t pz); |
| 76 |
|
| 77 |
void SetGood(Bool_t isgood) {fprim.fGOOD=isgood; }; |
| 78 |
|
| 79 |
// gen methods |
| 80 |
|
| 81 |
void GenPosition(Double_t xmin, Double_t xmax, Double_t ymin, Double_t ymax, Double_t zmin, Double_t zmax){ |
| 82 |
SetPosition(frnd->Uniform(xmin,xmax),frnd->Uniform(ymin,ymax),frnd->Uniform(zmin,zmax)); } |
| 83 |
//Spherical distribution in radius, thetamin, thetamax, phimin, phimax. |
| 84 |
void GenSphDist(Double_t, Double_t, Double_t, Double_t, Double_t); |
| 85 |
|
| 86 |
void GenDirection(Double_t thetamin, Double_t thetamax, Double_t phimin, Double_t phimax){ |
| 87 |
SetDirection(frnd->Uniform(thetamin,thetamax),frnd->Uniform(phimin,phimax)); } |
| 88 |
|
| 89 |
//flat spectra generator |
| 90 |
void GenSpe(Double_t PEmin, Double_t PEmax, Bool_t isEnergy=kFALSE); |
| 91 |
//power law spectra, gamma - differential spectral index |
| 92 |
void GenSpe(Double_t PEmin, Double_t PEmax, Double_t gamma, Bool_t isEnergy=kFALSE); |
| 93 |
|
| 94 |
//spherical dist. from radius r, allow only those particles going to cross 3 TOFs |
| 95 |
void GenSphDist(Double_t); |
| 96 |
|
| 97 |
|
| 98 |
void GenSphPhiThe(Double_t , Double_t , Double_t , Double_t , |
| 99 |
Double_t , Double_t ); // cecilia |
| 100 |
|
| 101 |
// get methods |
| 102 |
Int_t GetParticle(){ return fprim.fPDG; }; |
| 103 |
void GetPositon(TVector3 & v){ v.SetXYZ(fprim.fX0,fprim.fY0,fprim.fZ0); }; |
| 104 |
void GetPosition(Double_t X0, Double_t Y0,Double_t Z0) |
| 105 |
{X0=fprim.fX0; Y0=fprim.fY0; Z0=fprim.fZ0;}; |
| 106 |
|
| 107 |
Double_t GetKinEnergy() { return MomentumToKinE(fprim.fP0); }; |
| 108 |
Double_t GetRigidity() { return MomentumToRig(fprim.fP0); }; |
| 109 |
Double_t GetMomentum() { return fprim.fP0; }; |
| 110 |
|
| 111 |
void GetDirection(TVector2 & v) {v.Set(fprim.fTHETA,fprim.fPHI);} |
| 112 |
|
| 113 |
Bool_t Getgood(){ return fprim.fGOOD; }; |
| 114 |
|
| 115 |
//work with collection of primaries |
| 116 |
|
| 117 |
void Register(){ |
| 118 |
PamRootManager::Instance()-> |
| 119 |
Register("PRIM","TClonesArray", &fprimColl); |
| 120 |
} |
| 121 |
|
| 122 |
void ClearPrimCol(){ |
| 123 |
fevno=0; |
| 124 |
(*fprimColl).Clear(); |
| 125 |
} |
| 126 |
|
| 127 |
void Print(const Option_t * ="") const {fprim.Print();}; |
| 128 |
|
| 129 |
private: |
| 130 |
// methods |
| 131 |
|
| 132 |
Double_t MomentumToKinE(Double_t P0); |
| 133 |
Double_t KinEToMomentum(Double_t E0); |
| 134 |
Double_t MomentumToRig(Double_t P0) { return P0/fcharge; }; |
| 135 |
Double_t RigToMomentum(Double_t R0){ return R0*fcharge; }; |
| 136 |
|
| 137 |
// data members |
| 138 |
TVirtualMCStack* fStack; |
| 139 |
Int_t fevno; |
| 140 |
PamVMCPrimary fprim; |
| 141 |
Double_t fmass; |
| 142 |
Double_t fcharge; |
| 143 |
TClonesArray* fprimColl; |
| 144 |
|
| 145 |
TRandom3 * frnd; |
| 146 |
|
| 147 |
ClassDef(PamVMCPrimaryGenerator,2) //PamVMCPrimaryGenerator |
| 148 |
}; |
| 149 |
|
| 150 |
// inline functions |
| 151 |
|
| 152 |
|
| 153 |
inline Double_t PamVMCPrimaryGenerator::MomentumToKinE(Double_t P0) |
| 154 |
{ return Sqrt(P0*P0+fmass*fmass)-fmass; } |
| 155 |
|
| 156 |
inline Double_t PamVMCPrimaryGenerator::KinEToMomentum(Double_t E0) |
| 157 |
{ return Sqrt(E0*E0+2.*E0*fmass); } |
| 158 |
|
| 159 |
#endif //PAMVMC_PRIMARY_GENERATOR_H |
| 160 |
|
Parent Directory
| 
Revision Log