AraSim_old/Event.h 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
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
////////////////////////////////////////////////////////////////////////////////////////////////
//class Event:
////////////////////////////////////////////////////////////////////////////////////////////////

#ifndef EVENT_H
#define EVENT_H

#include <vector>
#include "Vector.h"
#include <string>
#include "Primaries.h"
#include "Report.h"
//#include "Neutrino.h"


using namespace std;

class Settings;
class Spectra;
class Primaries;
class Interaction;
//class Report;

class Event {

  public:

      int Event_type;   // 0 : neutrino event,  1 : blackhole event,  2 : monopole event,... etc

      int inu_thrown; // event number. in case we save triggered events only, this event number could be useful

      int inu_passed; // event number. in case we save triggered events only, this event number could be useful


      double pnu;   // energy of neutrino
      //Vector nnu;   // direction of neutrino
      string nuflavor;  // flavor of neutrino
      int nuflavorint;  // 1 : nue,  2 : numu,  3 : nutau
      int n_interactions;    // number of interaction inside the ice
      int IsCalpulser;

    void Initialize ();

      void Choose_Evt_Type(Settings *settings1);    // choose the event type depend on the settings->EVENT_TYPE value

      std::vector <Interaction> Nu_Interaction;

      //vector <Report> test_report;

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

      ~Event(); // destructor


      void delete_all(); // delete all vectors


//--------------------------------------------------
//   Event(int, Spectra*);
//   //  Event(double);
//   ~Event();
//   int iev; // index for this event
//   int part_type; // incident particle type (neutrino, magnetic monopole)
//   vector<double> n_incident;  // unit vector pointing in direction of incident particle trajectory (like a neutrino)
//   double e_incident;  // energy of particle incident at the interaction point
//
//   // something that describes the radiation that emerges from the shower
//   // this may be a root file or text file that is generated from jaime's program
//   vector<double> eField1m; // E field - 3d vector as a function of frequency at 1 m from the interaction
//
//   vector<double> getEField1m();
//
//
//  Event(double) : pos_sh(),int_type(0),iev(0),part_type(0),n_incident(3), e_incident(1.E18), eField1m(3) {}
//--------------------------------------------------

      ClassDef(Event,1);

};

#endif //EVENT_H


/* class Event { */
/* public: */

/* }; */

/* class NeutrinoEvent : public Event { */
/* public: */

/* }; */

/* class BlackHoleEvent : public Event { */
/* public: */

/* }; */

/* class MonopoleEvent : public Event { */
/* public: */

/* }; */
//class NeutrinoEvent;

/* template<class T, class U> */
/*   class NeutrinoEvent : public Event { */
/*  public: */
/*   NeutrinoEvent(T this_ev,U spectra1): Event(this_ev,spectra1) { */

/*   } */
/* }; */


/*
 * Create all available events and print their prices
 */
/* void event_information( EventFactory::EventType eventtype,int this_ev,Spectra spectra1 ) */
/* { */
/*   Event* event = EventFactory::createEvent(eventtype,this_ev,spectra1); */
/* 	//        std::cout << "Price of " << eventtype << " is " << event->getPrice() << std::endl; */
/*         delete event; */
/* } */