diff options
Diffstat (limited to 'fileio/root.C')
| -rw-r--r-- | fileio/root.C | 49 | 
1 files changed, 25 insertions, 24 deletions
| diff --git a/fileio/root.C b/fileio/root.C index 8f7d0a9..7052a71 100644 --- a/fileio/root.C +++ b/fileio/root.C @@ -4,10 +4,10 @@  #include <string>  struct Photon { -  double t; -  TVector3 pos; -  TVector3 dir; -  TVector3 pol; +  double time; +  TVector3 position; +  TVector3 direction; +  TVector3 polarization;    double wavelength; // nm    unsigned int history;    int last_hit_triangle; @@ -15,19 +15,20 @@ struct Photon {  struct Track {    std::string particle; -  double t; -  TVector3 pos; -  TVector3 dir; +  double time; +  TVector3 position; +  TVector3 direction;    double start_time;    double total_energy;  };  struct MC {    std::string particle; -  TVector3 gen_pos; -  TVector3 gen_dir; -  double gen_total_e; +  TVector3 gen_position; +  TVector3 gen_direction; +  double gen_total_energy; +  int nphoton;    std::vector<Track> subtrack;    std::vector<Photon> photon_start; @@ -36,10 +37,10 @@ struct MC {  };  struct Channel { -  Channel() : channel_id(-1), t(-9999.0), q(-9999.0) { }; +  Channel() : channel_id(-1), time(-9999.0), charge(-9999.0) { };    int channel_id; -  double t; -  double q; +  double time; +  double charge;    unsigned int mc_history;  }; @@ -65,8 +66,8 @@ struct Event {        if (channel_id < nentries) {  	hit[channel_id] = 1; -	time[channel_id] = channel[i].t; -	charge[channel_id] = channel[i].q; +	time[channel_id] = channel[i].time; +	charge[channel_id] = channel[i].charge;        }      }    } @@ -83,10 +84,10 @@ void fill_photons(Event *ev, bool start,    for (unsigned int i=0; i < nphotons; i++) {      Photon &photon = photons[i]; -    photon.t = t0[i]; -    photon.pos.SetXYZ(pos[3*i], pos[3*i + 1], pos[3*i + 2]); -    photon.dir.SetXYZ(dir[3*i], dir[3*i + 1], dir[3*i + 2]); -    photon.pol.SetXYZ(pol[3*i], pol[3*i + 1], pol[3*i + 2]); +    photon.time = t0[i]; +    photon.position.SetXYZ(pos[3*i], pos[3*i + 1], pos[3*i + 2]); +    photon.direction.SetXYZ(dir[3*i], dir[3*i + 1], dir[3*i + 2]); +    photon.polarization.SetXYZ(pol[3*i], pol[3*i + 1], pol[3*i + 2]);      photon.wavelength = wavelength[i];      if (histories)        photon.history = histories[i]; @@ -101,18 +102,18 @@ void fill_photons(Event *ev, bool start,  } -void fill_hits(Event *ev, unsigned int nchannels, float *t, -	       float *q, unsigned int *history) +void fill_hits(Event *ev, unsigned int nchannels, float *time, +	       float *charge, unsigned int *history)  {    ev->channel.resize(0);    ev->nhit = 0;    Channel ch;    for (unsigned int i=0; i < nchannels; i++) { -    if (t[i] < 1e8) { +    if (time[i] < 1e8) {        ev->nhit++;        ch.channel_id = i; -      ch.t = t[i]; -      ch.q = q[i]; +      ch.time = time[i]; +      ch.charge = charge[i];        ch.mc_history = history[i];        ev->channel.push_back(ch);      } | 
