CosmicsGenerator Class Reference

#include <CosmicsGenerator.h>

Inheritance diagram for CosmicsGenerator:

Collaboration diagram for CosmicsGenerator:

Public Member Functions
	CosmicsGenerator ()
virtual	~CosmicsGenerator ()
Bool_t	ReadEvent (FairPrimaryGenerator *)
virtual Bool_t	Init (Bool_t largeMom)

Public Attributes
double	z0
double	yBox
double	xBox
double	zBox
double	xdist
double	zdist
double	minE
int	n_EVENTS

Protected Member Functions
void	GenerateDynamics ()
Bool_t	DetectorBox ()
	ClassDef (CosmicsGenerator, 4)

Protected Attributes
double	P
double	px
double	py
double	pz
double	x
double	y
double	z
double	weighttest
double	weight
double	mass
double	FluxIntegral
double	theta
int	PID
int	nInside
int	nTest
Bool_t	high

Private Attributes
Co3Rng *	fRandomEngine

Detailed Description

Definition at line 38 of file CosmicsGenerator.h.

Constructor & Destructor Documentation

CosmicsGenerator()

CosmicsGenerator::CosmicsGenerator ( )

inline

constructor,destructor

Definition at line 41 of file CosmicsGenerator.h.

{};  

~CosmicsGenerator()

virtual CosmicsGenerator::~CosmicsGenerator ( )

inlinevirtual

Definition at line 42 of file CosmicsGenerator.h.

                               {
        delete fRandomEngine; 
        cout<<nInside<<" events have been generated."<<endl;
        cout<<"There is a total of "<<nInside<<"/"<<nTest<<" muons that passed close enough to the detector."<<endl;
        cout<<"Including the given weight this corresponds to ";
        cout<<weighttest/FluxIntegral/xdist/zdist*10000<<" spills (1 spill = "<<xdist*zdist*FluxIntegral/10000;
        cout<<" real cosmic muons = "<<n_EVENTS<<" simulated events)."<<endl;
    };

Member Function Documentation

ClassDef()

CosmicsGenerator::ClassDef ( CosmicsGenerator  , 4    )

protected

DetectorBox()

Bool_t CosmicsGenerator::DetectorBox ( )

protected

Definition at line 45 of file CosmicsGenerator.cxx.

                                    {
    // check, if a given staring setup x,y,z,px,py,pz will lead to a
    // close enough to the detector
 
    if ((TMath::Abs(x-(y+yBox)*px/py) < xBox && TMath::Abs(z-z0-(y+yBox)*pz/py) < zBox) ||
       (TMath::Abs(x-(y-yBox)*px/py) < xBox && TMath::Abs(z-z0-(y-yBox)*pz/py) < zBox) ||
       (TMath::Abs(y-(x+xBox)*py/px) < yBox && TMath::Abs(z-z0-(x+xBox)*pz/px) < zBox) ||
       (TMath::Abs(y-(x-xBox)*py/px) < yBox && TMath::Abs(z-z0-(x-xBox)*pz/px) < zBox)) {
        return true;
    }
    return false;
}

GenerateDynamics()

void CosmicsGenerator::GenerateDynamics ( )

protected

Definition at line 58 of file CosmicsGenerator.cxx.

                                       {
    // generate starting conditions for CM until the DetectorBox is hit
    do{
        weighttest += weight;   nTest++; //book keeping
        //momentum components
        double phi = fRandomEngine->Uniform(0,2*TMath::Pi());
        theta = fRandomEngine->fTheta->GetRandom();
        px = TMath::Sin(phi)*TMath::Sin(theta);
        pz = TMath::Cos(phi)*TMath::Sin(theta);
        py = -TMath::Cos(theta);
        //staring location
        x = fRandomEngine->Uniform(-xdist/2,xdist/2);
        z = fRandomEngine->Uniform(z0 - zdist/2, z0 + zdist/2);
    }while(!DetectorBox());
    nInside++;
}

Init()

Bool_t CosmicsGenerator::Init ( Bool_t  largeMom )

virtual

Definition at line 75 of file CosmicsGenerator.cxx.

                                            {
    //general
    fRandomEngine = new Co3Rng();
    TDatabasePDG* pdgBase = TDatabasePDG::Instance();
    mass = pdgBase->GetParticle(13)->Mass(); // muons!
    cout<<"----------------------------------------------------------------------"<<endl;
    cout<<"configuration for the CMBG as defined in $FAIRSHIP/python/CMBG_conf.py: "<<endl;
    cout<<"x_dist:  "<<xdist<<endl;
    cout<<"z_dist:  "<<zdist<<endl;
    cout<<"x_box:   "<<xBox<<endl;
    cout<<"y_box:   "<<yBox<<endl;
    cout<<"z_box:   "<<zBox<<endl;
    cout<<"n_EVENTS:"<<n_EVENTS<<endl;
    cout<<"minE:    "<<minE<<endl<<endl;
    if (xdist*zdist*n_EVENTS == 0){cout<<"check the configuration for unphysical behavior."<<endl<<"We stop the execution."<<endl<<endl; return kFALSE;}
 
    high = largeMom;
    if (high) cout<<"Simulation for high momentum"<<endl;
    else cout<<"Simulation for low momentum"<<endl;
 
    // calculating weights for this run
    // weight_flux: expected #muons per spill/ #simulated events per spill: FluxIntegral*xdist*zdist/EVENTS;
    //              the respective integrals are calculated from the fluxes
    // weight_DetectorBox: only consider CM hitting the DetectorBox
    //                     this is gained from a MC test of 10xEVENTS events
    double weight_flux, weight_DetectorBox;
    FluxIntegral = 0;
    if (!high) { // momentum range 1 GeV - 100 GeV
        if (minE > 100) {cout<<"choose minE < 100 !"<<endl; return kFALSE;}
        double dt = TMath::Pi()/2/100;
        for (double t= dt/2; t< TMath::Pi()/2; t += dt){
            FluxIntegral += fRandomEngine->fSpectrumL(t, minE, 0);
        }
        FluxIntegral = 2*TMath::Pi()/3*FluxIntegral*dt*10000;
        cout<< "LowE CM flux with P < minE = "<<minE<<" : "<<FluxIntegral<< "m-2s-1"<<endl;
    }
    else { // momentum range 100 GeV - 1000 GeV
        FluxIntegral = 2*TMath::Pi()/3*fRandomEngine->fSpectrumH->Integral(100,1000);
        cout<< "HighE CM flux: "<<FluxIntegral<< "m-2s-1"<<endl;
    }
    weight_flux = FluxIntegral*xdist*zdist/n_EVENTS/10000;
    nInside = 0;  nTest = 0; weighttest = 0; // book keeping
    y = 1900; //all muons start 19m over beam axis
    for (; nInside < 10*n_EVENTS;){
        GenerateDynamics();
    }
    weight_DetectorBox = 0.10 * nTest/n_EVENTS;
    weight = weight_flux / weight_DetectorBox;
    cout<<"weight_DetectorBox: "<< weight_DetectorBox<<", weight: "<< weight<<endl;
    cout<<"----------------------------------------------------------------------"<<endl<<endl;
    nInside = 0;  nTest = 0; weighttest = 0; // book keeping
 
    return kTRUE;
}

ReadEvent()

Bool_t CosmicsGenerator::ReadEvent

(

FairPrimaryGenerator *

cpg

)

public method ReadEvent

Definition at line 130 of file CosmicsGenerator.cxx.

                                                           {
    // muon or anti-muon
    PID = 26*(fRandomEngine->Uniform(0,1) < 1.0/2.278) - 13;
    // starting conditions
    GenerateDynamics();
    //momentum in the two regions, < or > 100 GeV
    if (!high) P = fRandomEngine->fSpectrumL(theta, minE);
    else P = fRandomEngine->fSpectrumH->GetRandom();
    px = px*P;
    py = py*P;
    pz = pz*P;
    // transfer to Geant4
    cpg->AddTrack(PID,px,py,pz,x,y,z,-1,true,TMath::Sqrt(P*P+mass*mass),0,weight);  // -1 = Mother ID, true = tracking, SQRT(x) = Energy, 0 = t
    if (!nInside%10000){cout<<nInside/10000<<"10k events have been simulated"<<endl;}
    return kTRUE;
}

Member Data Documentation

FluxIntegral

double CosmicsGenerator::FluxIntegral

protected

Definition at line 63 of file CosmicsGenerator.h.

fRandomEngine

Co3Rng* CosmicsGenerator::fRandomEngine

private

Definition at line 60 of file CosmicsGenerator.h.

high

Bool_t CosmicsGenerator::high

protected

Definition at line 65 of file CosmicsGenerator.h.

mass

double CosmicsGenerator::mass

protected

Definition at line 63 of file CosmicsGenerator.h.

minE

double CosmicsGenerator::minE

Definition at line 56 of file CosmicsGenerator.h.

n_EVENTS

int CosmicsGenerator::n_EVENTS

Definition at line 57 of file CosmicsGenerator.h.

nInside

int CosmicsGenerator::nInside

protected

Definition at line 64 of file CosmicsGenerator.h.

nTest

int CosmicsGenerator::nTest

protected

Definition at line 64 of file CosmicsGenerator.h.

P

double CosmicsGenerator::P

protected

Definition at line 63 of file CosmicsGenerator.h.

PID

int CosmicsGenerator::PID

protected

Definition at line 64 of file CosmicsGenerator.h.

px

double CosmicsGenerator::px

protected

Definition at line 63 of file CosmicsGenerator.h.

py

double CosmicsGenerator::py

protected

Definition at line 63 of file CosmicsGenerator.h.

pz

double CosmicsGenerator::pz

protected

Definition at line 63 of file CosmicsGenerator.h.

theta

double CosmicsGenerator::theta

protected

Definition at line 63 of file CosmicsGenerator.h.

weight

double CosmicsGenerator::weight

protected

Definition at line 63 of file CosmicsGenerator.h.

weighttest

double CosmicsGenerator::weighttest

protected

Definition at line 63 of file CosmicsGenerator.h.

x

double CosmicsGenerator::x

protected

Definition at line 63 of file CosmicsGenerator.h.

xBox

double CosmicsGenerator::xBox

Definition at line 56 of file CosmicsGenerator.h.

xdist

double CosmicsGenerator::xdist

Definition at line 56 of file CosmicsGenerator.h.

y

double CosmicsGenerator::y

protected

Definition at line 63 of file CosmicsGenerator.h.

yBox

double CosmicsGenerator::yBox

Definition at line 56 of file CosmicsGenerator.h.

z

double CosmicsGenerator::z

protected

Definition at line 63 of file CosmicsGenerator.h.

z0

double CosmicsGenerator::z0

Definition at line 56 of file CosmicsGenerator.h.

zBox

double CosmicsGenerator::zBox

Definition at line 56 of file CosmicsGenerator.h.

zdist

double CosmicsGenerator::zdist

Definition at line 56 of file CosmicsGenerator.h.

---

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

• shipgen/CosmicsGenerator.h
• shipgen/CosmicsGenerator.cxx