EventData.cc

00001 #include "EventData.hh"
00002 #include "Message.hh"
00003 
00004 void EventData::Print(int verbosity)
00005 {
00006     if (verbosity >= 1)
00007     {
00008         Message m(INFO);
00009         m<<std::endl;
00010         m<<"************************************************************************"<<std::endl;
00011         m<<"************************* EVENT INFORMATION ****************************"<<std::endl;
00012         m<<"************************************************************************"<<std::endl;
00013         m<<"Run ID: "<<run_id<<std::endl;
00014         m<<"Event ID: "<<event_id<<std::endl;
00015         m<<"Status: "<<status<<std::endl;
00016         m<<"Trigger Count: "<<trigger_count<<std::endl;
00017         m<<"Time Stamp: "<<timestamp<<std::endl;
00018         m<<"Dt: "<<dt<<std::endl;
00019         m<<"Event Time: "<<event_time<<std::endl;
00020         m<<"N. Chans: "<<nchans<<std::endl;
00021         m<<"Saturated: "<<saturated<<std::endl;
00022         m<<"S1 Valid: "<<s1_valid<<std::endl;
00023         m<<"S2 Valid: "<<s2_valid<<std::endl;
00024         m<<"S1S2 Valid: "<<s1s2_valid<<std::endl;
00025         m<<"S1 Start Time: "<<s1_start_time<<std::endl;
00026         m<<"S1 End Time: "<<s1_end_time<<std::endl;
00027         m<<"S2 Start Time: "<<s2_start_time<<std::endl;
00028         m<<"S2 End Time: "<<s2_end_time<<std::endl;
00029         m<<"S1 Full: "<<s1_full<<std::endl;
00030         m<<"S1 Fixed: "<<s1_fixed<<std::endl;
00031         m<<"S2 Full: "<<s2_full<<std::endl;
00032         m<<"S2 Fixed: "<<s2_fixed<<std::endl;
00033         m<<"Max S1: "<<max_s1<<std::endl;
00034         m<<"Max S2: "<<max_s2<<std::endl;
00035         m<<"Max S1 Channel: "<<max_s1_chan<<std::endl;
00036         m<<"Max S2 Channel: "<<max_s2_chan<<std::endl;
00037         m<<"F90 Full: "<<f90_full<<std::endl;
00038         m<<"F90 Fixed: "<<f90_fixed<<std::endl;
00039         m<<"Gatti: "<<gatti<<std::endl;
00040         m<<"LL_R: "<<ll_r<<std::endl;
00041         m<<"Drift Time: "<<drift_time<<std::endl;
00042         m<<"Position Valid: "<<position_valid<<std::endl;
00043         m<<"X Position: "<<x<<std::endl; 
00044         m<<"Y Position: "<<y<<std::endl;
00045         m<<"Z Position: "<<z<<std::endl;
00046         m<<"Barycentre Valid: "<<bary_valid<<std::endl;
00047         m<<"X Barycentre: "<<bary_x<<std::endl;
00048         m<<"Y Barycentre: "<<bary_y<<std::endl;
00049         m<<"************************************************************************"<<std::endl;
00050     }
00051     if (verbosity == 2)
00052     {
00053         {
00054             Message m(INFO);
00055             m<<std::endl;
00056             m<<"************************************************************************"<<std::endl;
00057             m<<"************************** CHANNEL SUMMARY *****************************"<<std::endl;
00058             m<<"************************************************************************"<<std::endl;
00059             m<<"\t\tCh. Baseline\t\t\t\t\t"<<std::endl;
00060             m<<"Ch.\tPulses\tFound\tMean\tInt.Min\tS1-Full\tS1-Fix\tS2-Full\tS2-Fix"<<std::endl;
00061             m<<"************************************************************************"<<std::endl;
00062             for (size_t ch = 0; ch < channels.size(); ch++)
00063             {
00064                 m<<std::setw(2)<<ch<<"\t"<<std::setw(2)<<channels[ch].npulses<<"\t"
00065                  <<std::setw(2)<<channels[ch].baseline.found_baseline<<"\t"<<std::setw(6)<< std::setprecision(4)<<channels[ch].baseline.mean<<"\t"      
00066                  <<std::setw(7)<< std::setprecision(5) <<channels[ch].integral_min<<"\t"
00067                  <<std::setw(4)<< std::setprecision(4) <<channels[ch].s1_full<<"\t"<<std::setw(4)<< std::setprecision(4) <<channels[ch].s1_fixed<<"\t"
00068                  <<std::setw(4)<< std::setprecision(4) <<channels[ch].s2_full<<"\t"<<std::setw(4)<< std::setprecision(4) <<channels[ch].s2_fixed
00069                  <<std::endl;
00070             }
00071             m<<"***********************************************************************"<<std::endl;
00072             //Print details for sum of integrals for each pulse
00073             for (size_t p = 0; p < sum_of_int.size(); p++)
00074             {
00075                 sum_of_int.at(p).Print(p);
00076             }
00077             //Print details for sum channel
00078             if (GetChannelByID(ChannelData::CH_SUM))
00079             {
00080                 m<<std::endl;
00081                 m<<"************************************************************************"<<std::endl;
00082                 m<<"*********************** SUM CHANNEL INFORMATION ************************"<<std::endl;
00083                 m<<"************************************************************************"<<std::endl;
00084             }
00085         }
00086         if (GetChannelByID(ChannelData::CH_SUM))
00087             GetChannelByID(ChannelData::CH_SUM)->Print(verbosity);
00088     }
00089     if (verbosity == 3)
00090     {
00091         {
00092             Message m(INFO);
00093             m<<std::endl;
00094             m<<"************************************************************************"<<std::endl;
00095             m<<"************************* CHANNEL INFORMATION **************************"<<std::endl;
00096             m<<"************************************************************************"<<std::endl;
00097         }
00098         for (size_t ch = 0; ch < channels.size(); ch++)
00099         {
00100             channels[ch].Print(verbosity);
00101         }
00102     }
00103 }