AraSim_old/Event.cc at master · ara-software/AraSim_old · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
#include "Event.h"
#include <string>
#include "Settings.h"
#include "Spectra.h"
#include "Primaries.h"
#include "TMath.h"
#include "Report.h"

ClassImp(Event);

using namespace std;


Event::Event () {
    //default constructor
    Initialize ();
}

void Event::Initialize () {

    Nu_Interaction.clear();
    //test_report.clear();
    n_interactions = 1;  // total number of interaction (including secondaries)

}


//Event::Event (Settings *settings1, Spectra *spectra1, Primaries *primary1, IceModel *icemodel, Detector *detector, Signal *signal, Secondaries *sec1 ) {
Event::Event (Settings *settings1, Spectra *spectra1, Primaries *primary1, IceModel *icemodel, Detector *detector, Signal *signal, Secondaries *sec1, int event_num ) {

    Initialize ();

    inu_thrown = event_num;

    Choose_Evt_Type (settings1);

    if (Event_type == 0) { // if only neutrino events exist


        pnu = spectra1->GetNuEnergy();
//        cout << pnu << endl;
/*
        double hereTheta = 10.;
        Vector output;
        output.SetX(-1.*TMath::Sin(hereTheta*3.1415926535/180.));
        output.SetY(0.);
        output.SetZ(TMath::Cos(hereTheta*3.1415926535/180.));
        nnu = output;
*/
        nuflavor = primary1->GetNuFlavor();

        /*
        if (settings1->NNU_THIS_THETA==1) {    // set specific theta angle for nnu
            nnu = primary1->GetThatDirection(settings1->NNU_THETA, settings1->NNU_D_THETA);
        }
        else { // nnu angle random
            nnu = primary1->GetAnyDirection();
        }
        */

        if (nuflavor=="nue")
            nuflavorint=1;
        else if (nuflavor=="numu")
            nuflavorint=2;
        else if (nuflavor=="nutau")
            nuflavorint=3;

        Interaction *Nu_temp;
        //Report *report_tmp;

        Nu_temp = new Interaction (pnu, nuflavor, n_interactions, icemodel, detector, settings1, primary1, signal, sec1 );
        //report_tmp = new Report(detector ,settings1);

        Nu_Interaction.push_back(*Nu_temp);  // for the first interaction
        //test_report.push_back(*report_tmp);

        delete Nu_temp;

        // for multiple interactions...
/*
        while (interaction_count < n_interactions) {    // not sure if this will work???

            Nu_tmp = new Interaction (...., n_interactions );

            Nu_Interaction.push_back( Nu_tmp );
        }
*/


    }

    IsCalpulser = primary1->IsCalpulser;


}

Event::~Event() {
    Nu_Interaction.clear();
}


void Event::delete_all() {

    Nu_Interaction.clear();
}


void Event::Choose_Evt_Type (Settings *settings1) {

    if (settings1->EVENT_TYPE==0) {
        //cout<<"Only Neutrino Evnets!"<<endl;
        Event_type = 0;
    }
    else {
        cout<<"Currently, only neutrino events possible!"<<endl;
        cout<<"Change Evt_type from "<<settings1->EVENT_TYPE<<" to 0"<<endl;
        Event_type = 0;
    }

}