OltEnds.cxx

Go to the documentation of this file.

// $Id: OltEnds.cxx,v 1.2 2003/01/20 10:13:58 west Exp $
//
// OltEnds.cxx
//
// Block Scan
//
// Author: Phil Adamson 5/2002

#include "JobControl/JobCModuleRegistry.h"
#include "MinosObjectMap/MomNavigator.h"
#include "OfflineTrigger/OltEnds.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 "Persistency/PerOutputStreamManager.h"
#include "Persistency/PerOutputStream.h"
#include "Rotorooter/RotoObjectifier.h"
#include "Plex/PlexHandle.h"
#include "Plex/PlexSEIdAltL.h"
#include "Plex/PlexStripEndId.h"
#include "Validity/VldContext.h"

#include "TString.h"
#include "TFile.h"
#include "TTree.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(OltEnds)

  // ......

CVSID("$Id: OltEnds.cxx,v 1.2 2003/01/20 10:13:58 west Exp $");
JOBMODULE(OltEnds,"OltEnds","A Hack");


OltEnds::OltEnds()
{
  fRecord = 0;
  TDirectory *save = gDirectory;
  fFile = new TFile("/home/pa/ends.root","RECREATE");
  fBothEnds = new TH2F("bothends","ADC: Both Ends Hit",500,0,10000,500,0,10000);
  save->cd();
}

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

JobCResult OltEnds::Reco(MomNavigator* mom)
{
  TIter iter = mom->FragmentIter();
  while (TObject *obj = iter.Next()) {
    RawRecord *rawrec = dynamic_cast<RawRecord *>(obj);
    if (rawrec) {
      ++fRecord;
      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);
            fTreeTF = rshb->GetTimeFrameNo();
            fVld = rshb->GetVldContext();
            MSG("Olt",Msg::kInfo)<<"Timeframe "<<fTreeTF<<endl;
          } 
          fTreeIndex = 0;
          RawDigitDataBlock *rddb = dynamic_cast<RawDigitDataBlock *>(rawdata);
          if (rddb) {
            multimap<int,RawDigit *> hitmap[2]; // map from (plane*25 + strip)
            map<int,int> adc[2];
            PlexHandle plex(fVld);
            TIter it = rddb->GetDatumIter();
            while (TObject *obj = it.Next()) {
              RawVaDigit *rvd = dynamic_cast<RawVaDigit *>(obj);
              if (rvd) {
                //cout<<*rvd<<endl;
                PlexSEIdAltL altl = plex.GetSEIdAltL(rvd->GetChannel());
                if (altl.size()>0) {
                  PlexStripEndId seid = altl.GetBestSEId();
                  //cout<<seid<<endl;
                  int plane = seid.GetPlane();
                  int strip = seid.GetStrip();
                  int end = int(seid.GetEnd());
                  if (end>0 && end<3) {
                    cout<<seid.GetEnd()<<endl;
                    cout<<"V "<<end<<" "<<plane<<" "<<strip<<endl;
                    hitmap[end-1].insert(pair<int,RawDigit *>(plane*25 + strip,rvd));
                    adc[end-1][plane*25 + strip] +=rvd->GetADC();
                  }
                }
              }
              RawQieDigit *rqd = dynamic_cast<RawQieDigit *>(obj);
              if (rqd) {
                //cout<<*rqd<<endl;
                PlexSEIdAltL altl = plex.GetSEIdAltL(rqd->GetChannel());
                if (altl.size()>0) {
                  PlexStripEndId seid = altl.GetBestSEId();
                  //cout<<seid<<endl;
                  int plane = seid.GetPlane();
                  int strip = seid.GetStrip();
                  int end = int(seid.GetEnd());
                  if (end>0 && end<3) {
                    cout<<seid.GetEnd()<<endl;
                    cout<<"Q "<<end<<" "<<plane<<" "<<strip<<endl;
                    hitmap[end-1].insert(pair<int,RawDigit *>(plane*25 + strip,rqd));           
                    adc[end-1][plane*25 + strip] +=(rqd->GetADC()-50);
                  }
                }
              }
            }
            // Now match up hits...
            map<int,int>::iterator mit[2];
            for (int i=0;i<2;i++) mit[i] = adc[i].begin();
            for (;;) {
             if (mit[0]==adc[0].end() && mit[1]==adc[1].end()) break;
             if (mit[0]==adc[0].end()) {
                // Side 1 only left:
                ++(mit[1]);
              }
              else if (mit[1] == adc[1].end()) {
                ++(mit[0]);
              }
              else {
                int ind0 = (*(mit[0])).first;
                int ind1 = (*(mit[1])).first;
                
                if (ind0==ind1) {
                  //Matched hits
                  fBothEnds->Fill((*(mit[0])).second,(*(mit[1])).second);
                  cout<<"Matched a hit"<<endl;
                  ++(mit[0]); ++(mit[1]);
                }
                else if (ind0 < ind1 ) {
                  // only on side 0
                  ++(mit[0]);
                }
                else if (ind0 > ind1) {
                  // only on side 1
                  ++(mit[1]);
                }
              }
             
            }
          }
          RawRunCommentBlock *rc = dynamic_cast<RawRunCommentBlock *>(rawdata);
          if (rc) cout<<rc->GetRunComment()<<endl;
        }
      }
    }
  }
  return JobCResult::kAOK;
}


 
const Registry& OltEnds::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.Set("DumpData", true);
  r.LockValues();

  return r;
}

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

void OltEnds::Config(const Registry& r)
{
  int   tmpb;
  cout<<"Teapots\n";
  cout<<r<<endl;
  if (r.Get("DumpData",tmpb)) { cout<<"Hello "<<tmpb<<endl;fDump = bool(tmpb); }
}

---