LIRawNt.cxx

Go to the documentation of this file.

#include <cstdio>
#include <fstream>
#include "LIRawNt.h"
// ROOT includes
#include "TFile.h"
// MINOS includes
#include "MessageService/MsgService.h"     // MSG text output
#include "MinosObjectMap/MomNavigator.h"   // Data access
#include "RawData/RawDigit.h"
#include "JobControl/JobCommand.h"         // JobCommand handling
#include "JobControl/JobCModuleRegistry.h" // JOBMODULE registration macro
#include "RawData/RawChannelId.h"
#include "RawData/RawDaqSnarlHeader.h"
#include "RawData/RawRecord.h"
#include "RawData/RawDaqHeaderBlock.h"
#include "RawData/RawSnarlHeaderBlock.h"
#include "RawData/RawDigitDataBlock.h"
#include "RawData/RawDigit.h"


CVSID("$Id: LIRawNt.cxx,v 1.3 2008/11/04 11:45:35 dja25 Exp $");
JOBMODULE(LIRawNt,"LIRawNt","Create Ntuple for list of raw LI hits\n");

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

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

LIRawNt::LIRawNt() 
  : fPlex(0)
{
  outFile=0;
  string sBase = "liraw";
  string sFileName=sBase+".root";
  ifstream test(sFileName.c_str());
  if(!test){
  outFile = new TFile(sFileName.c_str(),"RECREATE");
  }
  else{
    Int_t count=1;
    while(test){
      test.close();
      string sAppen=Form("%d",count);
      sFileName=sBase+"_"+sAppen+".root";
      test.open(sFileName.c_str());
      ++count;
    }
    outFile=new TFile(sFileName.c_str(),"NEW");
  }

  fLIRawNt = new TTree("liraw", "liraw");
  fLIRawNt -> Branch("run", &nt_run, "run/I");
  fLIRawNt -> Branch("snarl", &nt_snarl, "snarl/I");
  fLIRawNt -> Branch("timesec", &nt_timesec, "timesec/I");
  fLIRawNt -> Branch("timenanosec", &nt_timenanosec, "timenanosec/I");
  fLIRawNt -> Branch("errcode", &nt_errcode, "errcode/I");
  fLIRawNt -> Branch("crate", &nt_crate, "crate/I1");
  fLIRawNt -> Branch("varc", &nt_varc, "varc/I1");
  fLIRawNt -> Branch("vmm", &nt_vmm, "vmm/I1");
  fLIRawNt -> Branch("vaadc", &nt_vaadc, "vaadc/I1");
  fLIRawNt -> Branch("vachip", &nt_vachip, "vachip/I1");
  fLIRawNt -> Branch("vachannel", &nt_vachannel, "vachannel/I1");
  fLIRawNt -> Branch("pixel", &nt_pixel, "pixel/I1");
  fLIRawNt -> Branch("adc", &nt_adc, "adc/I2");
  fLIRawNt -> Branch("tdc", &nt_tdc, "tdc/I");
  fLIRawNt -> Branch("readouttype", &nt_readouttype, "readouttype/I");
  fLIRawNt -> Branch("pbox", &nt_pbox, "pbox/S");
  fLIRawNt -> Branch("pboxopp", &nt_pboxopp, "pboxopp/S");
  fLIRawNt -> Branch("led", nt_led, "led[8]/S");
  fLIRawNt -> Branch("ledopp", nt_ledopp, "led[8]/S");

}
JobCResult LIRawNt::Ana(const MomNavigator *mom) {

  const RawRecord *rawrec = 
    dynamic_cast<const RawRecord *>(mom -> GetFragment("RawRecord"));
  if (!rawrec) { 
    MSG("LIRawNt", Msg::kWarning) << "No RawRecord in MOM." << endl; 
    return JobCResult::kFailed;
  }
 static Int_t count =0;
  const VldContext *vldC = rawrec -> GetVldContext();
  if (!(vldC -> IsValid())) {
    MSG("LIRawNt", Msg::kWarning) << "Invalid context for RawRecord." << endl;
    return JobCResult::kFailed;
  }
 // if(count<10000){
        ++count;
  const Detector::Detector_t detector = vldC -> GetDetector();
  if (fPlex && !(fPlex -> GetVldRange().IsCompatible(vldC))) {
    delete fPlex;
  }
  if (!fPlex) {
    fPlex = new PlexHandle(*vldC);
    MSG("LIRawNt", Msg::kVerbose)
      << "New plex handle; validity"
      << " from " << fPlex -> GetVldRange().GetTimeStart().AsString()
      << " to " << fPlex -> GetVldRange().GetTimeEnd().AsString()
      << endl; 
  }

  const RawDaqSnarlHeader *snarlHdr = 
     dynamic_cast<const RawDaqSnarlHeader*>(rawrec -> GetRawHeader());
  if (!snarlHdr) {
    return JobCResult::kFailed;
  }

  nt_run   = snarlHdr -> GetRun();
  nt_snarl = snarlHdr -> GetSnarl();

  const RawSnarlHeaderBlock *rawsnarlHdrBlock = 
     dynamic_cast<const RawSnarlHeaderBlock *>
     (rawrec -> FindRawBlock("RawSnarlHeaderBlock"));

  nt_timesec     = -1;
  nt_timenanosec = -1;
  if (rawsnarlHdrBlock) {
    nt_timesec     = rawsnarlHdrBlock -> GetTriggerTime().GetSec();
    nt_timenanosec = rawsnarlHdrBlock -> GetTriggerTime().GetNanoSec();
  }

  // Remove this?
  const RawRecord *rr = 
    dynamic_cast<RawRecord *>(mom -> GetFragment("RawRecord"));
  if (!rr) {
    MSG("LIRawNt", Msg::kWarning) << "No RawRecord in MOM." << endl;
    return JobCResult::kFailed;
  }

  // Can this be pulled out from the plex?
  const Int_t vach2pixel[18] = 
    {0, 0, 15, 1, 16, 2, 11, 5, 12, 6, 7, 9, 8, 10, 3, 14, 4, 13};

  TIter rdbit = rr -> GetRawBlockIter();
  TObject *tob;
  while ((tob = rdbit())) {
    //          ++count;//
//      cout<<count<<endl;
        if(count>10)break;//
    if (tob -> InheritsFrom("RawDigitDataBlock")) {
      const RawDigitDataBlock *rdb = 
        dynamic_cast<const RawDigitDataBlock *>(tob);
      TIter rdit = rdb -> GetDatumIter();
      const RawDigit *rd = dynamic_cast<const RawDigit *>(rdit());
      while (rd) {
        nt_pbox = -1;
        nt_pboxopp = -1;
        for (Int_t i = 0; i < nStripPerPixel; ++i) {
          nt_led[i] = -1; nt_ledopp[i] = -1;
        }
        const RawChannelId rawid = rd -> GetChannel();
        if (rawid.GetElecType() == ElecType::kVA/*&&count%10==0*/) {
          nt_crate       = rawid.GetCrate();
          nt_varc        = rawid.GetVarcId();
          nt_vmm         = rawid.GetVmm();
          nt_vaadc       = rawid.GetVaAdcSel();
          nt_vachip      = rawid.GetVaChip();
          nt_vachannel   = rawid.GetVaChannel();
          nt_pixel       = -1;
          if ((nt_vachannel >= 2) && (nt_vachannel <= 17)) {
            nt_pixel     = vach2pixel[nt_vachannel];
          }
          nt_errcode     = rd -> GetErrorCode();
          nt_adc         = rd -> GetADC();
          nt_tdc         = rd -> GetTDC();
          nt_readouttype = static_cast<Int_t>(fPlex -> GetReadoutType(rawid));
          if (nt_readouttype == ReadoutType::kScintStrip) {
            PlexSEIdAltL stripEndList = fPlex -> GetSEIdAltL(rawid);
            Int_t plane = stripEndList.GetPlane();
            if (plane > 485) {} // Ignore shield "planes"
            else if (detector == Detector::kCalDet &&
                     plane > 59) {} // Ignore CalDet muon counter
            else {
              Int_t i = 0;
              for (PlexSEIdAltL::iterator it = stripEndList.begin(); 
                   it != stripEndList.end();
                   it++) {
                PlexStripEndId stripEndId    = (*it).GetSEId();
                PlexStripEndId stripEndIdOpp = stripEndId.GetOppositeSEId();
                PlexLedId ledId    = fPlex -> GetLedId(stripEndId);
                PlexLedId ledIdOpp = fPlex -> GetLedId(stripEndIdOpp);
                if (it == stripEndList.begin()) {
                  nt_pbox = ledId.GetPulserBox();
                  nt_pboxopp = ledIdOpp.GetPulserBox();
                }
                if (i < nStripPerPixel) {
                  nt_led[i] = ledId.GetLedInBox();
                  nt_ledopp[i] = ledIdOpp.GetLedInBox();
                }
                i++;
              }
            }
          }
          else {}
         //       if(count<100000){
          //        ++count;
            fLIRawNt -> Fill();
          //  }
          //      else{return false;}
        }
        rd = dynamic_cast<const RawDigit *>(rdit());
//      }//
      }
    }
  }

  return JobCResult::kPassed;
}


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

void LIRawNt::HandleCommand(JobCommand* cmd)
{
//======================================================================
// Handle a command send from the job controller
// 
// Inputs: cmd - The parsed job command
//======================================================================
  if (cmd->HaveCmd()) {                      // If we have a command...
    string sc = cmd->PopCmd();               // Get the command
    MSG("LIRawNt",Msg::kWarning) << "Don't understand '" << sc.c_str() << "'\n";
  }
}


void LIRawNt::EndJob() 
{
//======================================================================
// At the end of the job print some stuff out and save the histogram
// to a file
//======================================================================
  // fFile->Write();
  outFile->Close();
}

---