gar::evd::RawDataDrawer Class Reference

Aid in the rendering of RawData objects. More...

#include <RawDataDrawer.h>

Public Member Functions
	RawDataDrawer ()
	~RawDataDrawer ()
void	RawDigit3D (art::Event const &evt, evdb::View3D *view)
	Draws raw digit content in 2D wire plane representation. More...

Private Member Functions
void	FillQHisto (gar::raw::RawDigit const &dig, TH1F *histo)
void	DrawRawDigit3D (art::Event const &evt, evdb::View3D *view)
	Performs the 3D drawing. More...
void	GetRawDigits (art::Event const &evt, std::vector< const raw::RawDigit * > &digits)
	Makes sure raw::RawDigit's are available for the current settings. More...

Private Attributes
TH1F *	fEventTQHist
	Charge vs time for all digits in an event. More...
TH1F *	fDigitTQHist
	Charge vs time for a single channel in an event. More...

Detailed Description

Aid in the rendering of RawData objects.

Definition at line 27 of file RawDataDrawer.h.

Constructor & Destructor Documentation

gar::evd::RawDataDrawer::RawDataDrawer

(

)

Definition at line 70 of file RawDataDrawer.cxx.

      {
        fEventTQHist = (TH1F*) new TH1F("eventTQHist","ADC Sum vs Time, All Channels;time (ticks);ADC Sum",100,0,10000);
        fDigitTQHist = (TH1F*) new TH1F("digitTQHist","ADC Sum vs Time, One Channel;time (ticks);ADC Sum",100,0,10000);
      }

gar::evd::RawDataDrawer::~RawDataDrawer

(

)

Definition at line 77 of file RawDataDrawer.cxx.

      {
      }

Member Function Documentation

void gar::evd::RawDataDrawer::DrawRawDigit3D ( art::Event const &  evt, evdb::View3D *  view  )

private

Performs the 3D drawing.

void gar::evd::RawDataDrawer::FillQHisto ( gar::raw::RawDigit const &  dig, TH1F *  histo  )

private

Definition at line 82 of file RawDataDrawer.cxx.

        {

          gar::raw::ADCvector_t uncompressed;
          uncompressed.resize(dig.Samples());
          short ped = dig.Pedestal();
          gar::raw::Uncompress(dig.ADCs(),uncompressed,ped,dig.Compression());
          for(size_t t = 0; t < uncompressed.size(); ++t)
          histo->Fill(t, 1. * uncompressed[t] - 1. * dig.Pedestal());

          return;
        }//end loop over raw hits

void gar::evd::RawDataDrawer::GetRawDigits ( art::Event const &  evt, std::vector< const raw::RawDigit * > &  digits  )

private

Makes sure raw::RawDigit's are available for the current settings.

Definition at line 206 of file RawDataDrawer.cxx.

            {
              digits.clear();

              art::ServiceHandle<evd::RawDrawingOptions> rawopt;

              std::vector<const raw::RawDigit*> temp;

              try{
                evt.getView(rawopt->fRawDataLabel, temp);
                for(size_t t = 0; t < temp.size(); ++t){
                  digits.push_back(temp[t]);
                }
              }
              catch(cet::exception& e){
                writeErrMsg("GetDigits", e);
              }

              return;
            } // RawDataDrawer::GetRawDigits()

void gar::evd::RawDataDrawer::RawDigit3D	(	art::Event const &	evt,
		evdb::View3D *	view
	)

Draws raw digit content in 2D wire plane representation.

Parameters

evt	source for raw digits
view	target rendered object
plane	number of the plane to be drawn
bZoomToRoI	whether to render only te region of interest

This function performs pre-rendering of the raw digit content into a 3D view of the yzt space The material for rendering is created and sent to view object for actual rendering.

Definition at line 99 of file RawDataDrawer.cxx.

          {
            // Check if we're supposed to draw raw hits at all
            art::ServiceHandle<evd::RawDrawingOptions>   rawopt;
            art::ServiceHandle<evd::ColorDrawingOptions> cdopt;
            art::ServiceHandle<geo::GeometryGAr>            geom;
            auto fTime = gar::providerFrom<detinfo::DetectorClocksServiceGAr>();
            auto detProp = gar::providerFrom<detinfo::DetectorPropertiesService>();

            double driftVelocity = detProp->DriftVelocity(detProp->Efield(),
            detProp->Temperature());

            if(rawopt->fDrawRawOrReco > 0) return;

            std::vector<const raw::RawDigit*> rawhits;
            this->GetRawDigits(evt, rawhits);

            fEventTQHist->Reset();
            fDigitTQHist->Reset();

            //auto const& colorSet = cdopt->RawQ();

            gar::raw::ADCvector_t uncompressed;
            for(auto dig : rawhits){

              this->FillQHisto(*dig, fEventTQHist);
              if(dig->Channel() == rawopt->fChannel)
              this->FillQHisto(*dig, fDigitTQHist);

              // draw the digit as time, y, z with the latter coordinates
              // given by the channel map.
              float xyz[3];

              geom->ChannelToPosition((unsigned int)dig->Channel(),
              xyz);

              //std::cout << "RDGeo: " << dig->Channel() << " " << xyz[0] << " " << xyz[1] << " " << xyz[2] << std::endl;
              double chanposx = xyz[0];

              // draw hit channels on the readout planes in gray.  No timing or x information.
              //int c = kGray;
              //int s=1;
              //int w=1;
              //TPolyLine3D& rdpos = view->AddPolyLine3D(2,c,w,s);
              //rdpos.SetPoint(0,xyz[2],xyz[0],xyz[1]);
              //rdpos.SetPoint(1,xyz[2]+1,xyz[0],xyz[1]);
              //// std::cout << "adding a raw digit: " << xyz[0] << " " << xyz[1] << " " << xyz[2] << std::endl;

              // loop over the ADC values for this digit and draw the value if it
              // is above the threshold
              //TMarker3DBox *box = nullptr;

              uncompressed.resize(dig->Samples());
              short ped = dig->Pedestal();
              gar::raw::Uncompress(dig->ADCs(),uncompressed,ped,dig->Compression());

              for(size_t t = 0; t < dig->Samples(); ++t){
                if(uncompressed[t] < rawopt->fMinSignal) continue;

                double driftdistance = fTime->TPCTick2Time(t) * driftVelocity;

                if (chanposx < 0)
                {
                  xyz[0] = chanposx + driftdistance;
                }
                else
                {
                  xyz[0] = chanposx - driftdistance;
                }

                // nb. coordinates are so Y is up, but ROOT thinks Z is up, so report (z,x,y)
                // somehow these boxes don't work -- draw short lines instead
                //box = &(view->AddMarker3DBox(xyz[2],
                //                            xyz[0],
                //                            xyz[1],
                //                            0.5,    // the extent is 1/2 tick
                //                             0.5 * 0.3, // to fix  geom->ChannelPitch(),
                //                         0.5 * 0.3  // to fix geom->ChannelPitch()
                //                             ));
                //box->SetFillStyle(1001);
                ////auto tmpcolor = colorSet.GetColor(uncompressed[t]);
                //box->SetFillColor(kGreen);
                //box->SetBit(kCannotPick);

                //std::cout << "adding a raw digit: " << xyz[0] << " " << xyz[1] << " " << xyz[2] << std::endl;

                // try drawing with polylines, little green segments


                int c2 = kGreen;
                int s2=1;
                int w2=1;

                // nb. coordinates are so Y is up, but ROOT thinks Z is up, so report (z,x,y)

                TPolyLine3D& rdpos = view->AddPolyLine3D(2,c2,w2,s2);
                rdpos.SetPoint(0,xyz[2],xyz[0]-0.3,xyz[1]);
                rdpos.SetPoint(1,xyz[2],xyz[0]+0.3,xyz[1]);

              } // end loop over ADC values for the digit
            }//end loop over raw digits

            return;
          }

Member Data Documentation

TH1F* gar::evd::RawDataDrawer::fDigitTQHist

private

Charge vs time for a single channel in an event.

Definition at line 67 of file RawDataDrawer.h.

TH1F* gar::evd::RawDataDrawer::fEventTQHist

private

Charge vs time for all digits in an event.

Definition at line 66 of file RawDataDrawer.h.

---

The documentation for this class was generated from the following files:

• garsoft/EventDisplay/EVD/RawDataDrawer.h
• garsoft/EventDisplay/EVD/RawDataDrawer.cxx