OltTimePlot.cxx

Go to the documentation of this file.

// $Id: OltTimePlot.cxx,v 1.6 2006/05/22 19:33:13 rhatcher Exp $
//
// OltBlockPlot.cxx
//
// Block Plot
//
// Author: Phil Adamson 5/2002

#include "JobControl/JobCModuleRegistry.h"
#include "MinosObjectMap/MomNavigator.h"
#include "OfflineTrigger/OltTimePlot.h"
#include "MessageService/MsgService.h"
#include "RawData/RawRecord.h"
#include "RawData/RawHeader.h"
#include "RawData/RawDataBlock.h"
#include "RawData/RawDaqHeaderBlock.h"
#include "RawData/RawSnarlHeaderBlock.h"
#include "RawData/RawChannelId.h"
#include "RawData/RawVaDigit.h"
#include "RawData/RawQieDigit.h"
#include "RawData/RawDigitDataBlock.h"
#include "RawData/RawRunCommentBlock.h"
#include "RawData/RawDeadChipBlock.h"

#include "Persistency/PerOutputStreamManager.h"
#include "Persistency/PerOutputStream.h"
#include "Rotorooter/RotoObjectifier.h"
#include "Plex/PlexHandle.h"
#include "Plex/PlexSEIdAltL.h"
#include "Validity/VldContext.h"

#include "TString.h"
#include "TFile.h"
#include "TH1.h"
#include "TH2.h"
#include <string>

#include <algorithm>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <fstream>
#include <map>
#ifndef __CINT__
#include "OnlineUtil/mdBlockDefinitions.h"
#endif
ClassImp(OltTimePlot)

 // ......

//CVSID("$Id: OltTimePlot.cxx,v 1.6 2006/05/22 19:33:13 rhatcher Exp $");
JOBMODULE(OltTimePlot,"OltTimePlot","The Offline Trigger Time Plotner");


OltTimePlot::OltTimePlot()
{
  TDirectory *save = gDirectory;
  fFile = new TFile("/home/pa/timeplot.root","RECREATE");
  save->cd();
  fPlex = 0;
}

OltTimePlot::~OltTimePlot()
{
  TDirectory *save = gDirectory;
  fFile->cd();
  fFile->Close();
  save->cd();
}

JobCResult OltTimePlot::Reco(MomNavigator* mom)
{
  int snarl=0,tf=0;
  TIter iter = mom->FragmentIter();
  while (TObject *obj = iter.Next()) {
    RawRecord *rawrec = dynamic_cast<RawRecord *>(obj);
    if (rawrec) {
      TIter recit = rawrec->GetRawBlockIter();
      while (TObject *obj = recit.Next()) {
        RawDataBlock *rawdata = dynamic_cast<RawDataBlock *>(obj);
        if (rawdata) {
          if (rawdata->InheritsFrom("RawSnarlHeaderBlock")) {
            RawSnarlHeaderBlock *rshb =dynamic_cast<RawSnarlHeaderBlock *>(rawdata);
             tf = rshb->GetTimeFrameNo();
             snarl = rshb->GetSnarl();
             if (!fPlex) fPlex = new PlexHandle(rshb->GetVldContext());
          }
          
          RawDigitDataBlock *rddb = dynamic_cast<RawDigitDataBlock *>(rawdata);
          if (rddb) {
            RawChannelId id(Detector::kCalDet,
                            ElecType::kVA,0,0x1602);
            id.SetModeBits(true,true,true);
            map<int,int> times;
            int t0 = 0;
            TIter it = rddb->GetDatumIter();
            while (TObject *obj = it.Next()) {
              RawVaDigit *rvd = dynamic_cast<RawVaDigit *>(obj);
              if (rvd) {
                RawChannelId rcid = rvd->GetChannel();
                PlexSEIdAltL altl= fPlex->GetSEIdAltL(rcid);

                times.insert(make_pair(rvd->GetTDC(),altl.GetPlane()));
//cout<<rvd->GetChannel()<<" "<<id<<endl;
                if (rcid == id) t0 = rvd->GetTDC();
              }
            }
            TDirectory *save = gDirectory;
            char hname[100];
            sprintf(hname,"snarl%d",snarl);
            TH1F h(hname,hname,200,-1000,1000);
            sprintf(hname,"psnarl%d",snarl);
            TH2F hh(hname,hname,50,-1000,1000,30,0,60);
            for (map<int,int>::iterator it =
times.begin();it!=times.end();it++) {
              h.Fill((*it).first-t0);
              hh.Fill((*it).first-t0,(*it).second);
//cout<<(*it)<<"  "<<t0<<endl;
            }
            fFile->cd();
            h.Write();
            hh.Write();
            save->cd();
          }
        }
      }
    }
  }
  return JobCResult::kAOK;
}
 
const Registry& OltTimePlot::DefaultConfig() const
{

  static Registry r; // Default configuration for module


  std::string name = this->JobCModule::GetName();
  name += ".config.default";
  r.SetName(name.c_str());

  // Set values in configuration
  r.UnLockValues();

  r.LockValues();

  return r;
}

//......................................................................

void OltTimePlot::Config(const Registry& /* r */)
{
  //int   tmpb;

}

---