#include <MRCCAna.h>
Public Member Functions | |
| MRCCAna () | |
| ~MRCCAna () | |
| void | ana (ParticleObjectHolder *poh, MRCC *e) |
|
|
Definition at line 6 of file MRCCAna.cxx. 00007 {}
|
|
|
Definition at line 9 of file MRCCAna.cxx. 00010 {}
|
|
||||||||||||
|
Definition at line 13 of file MRCCAna.cxx. References Particle3D::avg_rms_t, MRCC::avg_rms_t, Particle3D::calibrated_energy, MRCC::calibrated_energy, Particle3D::end_u, MRCC::end_u, Particle3D::end_v, MRCC::end_v, Particle3D::end_z, MRCC::end_z, Particle3D::entries, MRCC::entries, MRCC::hasMRCC, ParticleObjectHolder::mrcc, Particle3D::muonfrac, MRCC::muonfrac, MRCC::particle_s, Particle3D::phi, MRCC::phi, MRCCInfo::removedmuon, MRCCInfo::stage, MRCC::stage, Particle3D::start_u, MRCC::start_u, Particle3D::start_v, MRCC::start_v, Particle3D::start_z, MRCC::start_z, Particle3D::sum_e, MRCC::sum_e, Particle3D::theta, MRCC::theta, Particle3D::u, Particle3D::v, and Particle3D::z. Referenced by PRecordAna::ana(). 00014 {
00015 MRCCInfo * mi = poh->mrcc;
00016 if(!mi)
00017 {
00018 e->hasMRCC=0;
00019 return;
00020 }
00021
00022 Particle3D *p=&mi->removedmuon;
00023 if(!p)
00024 {
00025 e->hasMRCC=0;
00026 return;
00027 }
00028
00029 e->hasMRCC=1;
00030
00031 e->stage = mi->stage;
00032 e->start_u=p->start_u;
00033 e->end_u=p->end_u;
00034 e->start_v=p->start_v;
00035 e->end_v=p->end_v;
00036 e->start_z=p->start_z;
00037 e->end_z=p->end_z;
00038
00039 e->sum_e=p->sum_e;
00040
00041 e->muonfrac=p->muonfrac;
00042 e->entries=p->entries;
00043
00044 e->theta=p->theta;
00045 e->phi=p->phi;
00046
00047 e->calibrated_energy=p->calibrated_energy;
00048 e->avg_rms_t=p->avg_rms_t;
00049
00050
00051
00052 double particle_s=0;
00053 for(int j=0;j<p->entries;j++)
00054 {
00055 //double u=p->u[j]-poh->event.vtx_u;
00056 //double v=p->v[j]-poh->event.vtx_v;
00057
00058
00059 if(j>0)
00060 {
00061 double du = p->u[j]-p->u[j-1];
00062 double dv = p->v[j]-p->v[j-1];
00063 double dz = p->z[j]-p->z[j-1];
00064
00065 double dd=sqrt(du*du+dv*dv+dz*dz);
00066 if(dd<100)particle_s+=dd; //in case we have a single bogus value, we will skip it
00067 }
00068
00069 }
00070
00071 e->particle_s=particle_s;
00072
00073 }
|
1.3.9.1