SexyPedReco.cxx

Go to the documentation of this file.

#define SexyPedReco_cxx

#include "SexyPedReco.h"

using namespace std;

CVSID("$Id: SexyPedReco.cxx,v 1.8 2007/11/11 05:27:55 rhatcher Exp $");
JOBMODULE(SexyPedReco,"SexyPedReco","");

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

SexyPedReco::SexyPedReco() {

 fFile = 0;

 myPedAna = 0;
 
 fSexyPedTree = 0;
 fPedTree     = 0;
 fSparTree    = 0;

 //fGetTableTree = true;
 //fStoreTree    = true;
 //fOnline       = false;

 fGotCalledPed  = false;
 fGotCalledSpar = false;
}

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

SexyPedReco::~SexyPedReco() {

  return;
}

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

JobCResult SexyPedReco::Ana(const MomNavigator *mom) {

  if(!fFile) {
    char filename[180];
    RawRecord *rr = dynamic_cast<RawRecord *>(mom->GetFragment("RawRecord"));
    if (rr == 0) return JobCResult::kFailed;
      
    const RawDaqHeader* daqHdr = dynamic_cast<const RawDaqHeader*>(rr->GetRawHeader());      

    if (daqHdr) frunNum = daqHdr->GetRun();

    sprintf(filename,"SexyPedestalTree_%d.root",frunNum);


    if(fStoreTree) fFile = new TFile(filename,"RECREATE");

    if(fFile!=0) {
      fSexyPedTree = new TTree("sexyPedTree","SexyPedestal TTree");
      
      fSexyPedTree->Branch("run",  &frunNum,  "run/I");
      fSexyPedTree->Branch("error",&ferror,   "error/I");
      fSexyPedTree->Branch("crate",&fcrate,   "crate/I");
      fSexyPedTree->Branch("varc", &fvarc,    "varc/I");
      fSexyPedTree->Branch("vmm",  &fvmm,     "vmm/I");
      fSexyPedTree->Branch("vfb",  &fvfb,     "vfb/I");
      fSexyPedTree->Branch("va",   &fva,      "va/I");
      fSexyPedTree->Branch("vach", &fvach,    "vach/I");
      fSexyPedTree->Branch("adc",  &fadc,     "adc/I");
      fSexyPedTree->Branch("tdc",  &ftdc,     "tdc/D");
      fSexyPedTree->Branch("plane",&fplane,   "plane/I");
      fSexyPedTree->Branch("tUTC", &ftimeUTC, "tUTC/D");
      
      if(fGetTableTree) {
        
        fPedTree = new TTree("pedTableTree","Online Pedestals Table TTree");
        fPedTree->Branch("crate",&icrate,     "crate/I");
        fPedTree->Branch("varc", &ivarc,      "varc/I");
        fPedTree->Branch("vmm",  &ivmm,       "vmm/I");
        fPedTree->Branch("vfb",  &ivfb,       "vfb/I");
        fPedTree->Branch("va",   &iva,        "va/I");
        fPedTree->Branch("vach", &ivach,      "vach/I");
        fPedTree->Branch("adc",  &fPedTable,  "adc/I");
        fPedTree->Branch("time", &fLoadTime,  "time/D");
        
        fSparTree = new TTree("sparTableTree","Sparsicification Table TTree");
        fSparTree->Branch("crate",&icrate,    "crate/I");
        fSparTree->Branch("varc", &ivarc,     "varc/I");
        fSparTree->Branch("vmm",  &ivmm,      "vmm/I");
        fSparTree->Branch("vfb",  &ivfb,      "vfb/I");
        fSparTree->Branch("va",   &iva,       "va/I");
        fSparTree->Branch("vach", &ivach,     "vach/I");
        fSparTree->Branch("adc",  &fSparTable,"adc/I");
        fSparTree->Branch("time", &fLoadTime, "time/D");
      }
    }
  }

  TObject* tobj;
  TIter    fragiter = mom->FragmentIter();

  while( ( tobj = fragiter.Next() ) ) { //Loop over MOM fragments
    RawRecord* rawrec;
    if( ( rawrec = dynamic_cast<RawRecord*>(tobj) ) ) {  //Loop over RawRecords

      //RAWSNARLHEARDERBLOCK info:
      const RawSnarlHeaderBlock* snarlHdr = 
        dynamic_cast<const RawSnarlHeaderBlock*> (rawrec->FindRawBlock("RawSnarlHeaderBlock"));
      if (snarlHdr) {
        ftimeUTC = snarlHdr->GetTriggerTime().GetSec(); //From 1/1/1970
        MSG("SexyPedReco", Msg::kInfo) << "Timeframe: " << snarlHdr->GetTimeFrameNum() << endl;
      }
      else MSG("SexyPedReco", Msg::kDebug) << "There was not RawSnarlHeaderBlock this time!" << endl;
      
      //Loop over other blocks:
      TIter rdbit = rawrec->GetRawBlockIter();
      TObject *tob;
      while ((tob = rdbit())) {
        
        //Pedestal Table read in:
        RawPedestalTableBlock* rptb;
        if( tob->InheritsFrom("RawPedestalTableBlock") ) {
          rptb = dynamic_cast<RawPedestalTableBlock*>(tob);
          
          this->GetPedestalTable(rptb);
        }
        
        //Sparsifier Table read in:
        RawSparsifierTableBlock* rstb;
        if( tob->InheritsFrom("RawSparsifierTableBlock") ) {
          rstb = dynamic_cast<RawSparsifierTableBlock*>(tob);

          this->GetSparsificationTable(rstb);

        }

        //RawDigits get into SexyPedestal TTree:
        RawDigitDataBlock* rdb;
        if( tob->InheritsFrom("RawDigitDataBlock") ) {
          rdb = dynamic_cast<RawDigitDataBlock*>(tob);

          this->GetSexyPedestals(rdb);
        }
        //return JobCResult::kFailed; 
      }
    }
  }

  return JobCResult::kPassed;
}

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

void SexyPedReco::GetPedestalTable(RawPedestalTableBlock *rptb) {
  
  if(!fGetTableTree || fGotCalledPed) return;

  MSG("SexyPedReco", Msg::kDebug) << endl << "Creating the Pedestal TTree" << endl << endl;
  
  fGotCalledPed = true;
  
  //Load time of the table:
  fLoadTime = rptb->GetLoadTime().GetSec();
  
  MSG("SexyPedReco", Msg::kInfo) << "Pedestal Table Loaded on: " << endl;
  rptb->GetLoadTime().Print();
  cout << endl;
  
  for(icrate = 0; icrate < 16; icrate++) {
    for(ivarc = 0; ivarc < 3; ivarc++) {
      for(ivmm = 0; ivmm < 6; ivmm++) {
        for(ivfb = 0; ivfb < 2; ivfb++) {
          for(iva = 0; iva < 3; iva++) {
            for(ivach = 0; ivach < 22; ivach++) {
              
              const RawChannelId thisRawChannel(Detector::kFar,
                                                ElecType::kVA,
                                                icrate, ivarc, ivmm, ivfb, iva, ivach);
              
              fPedTable = rptb->GetPedValue(thisRawChannel);

              MSG("SexyPedReco", Msg::kDebug) << fPedTable << "\t" 
                                              << icrate    << "\t" 
                                              << ivarc    << "\t" 
                                              << ivmm     << "\t" 
                                              << ivfb     << "\t" 
                                              << iva      << "\t" 
                                              << ivach    << endl;
              
              fPedTree->Fill();
            }
          }
        }
      }
    }
  }
  return;
}

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

void SexyPedReco::GetSparsificationTable(RawSparsifierTableBlock *rstb) {

  if(!fGetTableTree || fGotCalledSpar) return;

  MSG("SexyPedReco", Msg::kDebug) << "Creating the Sparsification TTree" << endl << endl;
  
  fGotCalledSpar = true;
  
  //Load time of the table:           
  fLoadTime = rstb->GetLoadTime().GetSec();
  
  MSG("SexyPedReco", Msg::kInfo) << "Sparsification Table Loaded on: " << endl; 
  rstb->GetLoadTime().Print();
  cout << endl;
    
  for(icrate = 0; icrate < 16; icrate++) {
    for(ivarc = 0; ivarc < 3; ivarc++) { 
      for(ivmm = 0; ivmm < 6; ivmm++) { 
        for(ivfb = 0; ivfb < 2; ivfb++) {
          for(iva = 0; iva < 3; iva++) {
            for(ivach = 0; ivach < 22; ivach++) {
              
              RawChannelId thisRawChannel(Detector::kFar,
                                          ElecType::kVA,
                                          icrate, ivarc, ivmm, ivfb, iva, ivach);
              
              fSparTable = rstb->GetThresholdValue(thisRawChannel);
                                
              MSG("SexyPedReco", Msg::kDebug) << fSparTable << "\t" 
                                              << icrate    << "\t" 
                                              << ivarc    << "\t" 
                                              << ivmm     << "\t" 
                                              << ivfb     << "\t" 
                                              << iva      << "\t" 
                                              << ivach    << endl;

              fSparTree->Fill();
            }
          }
        }
      }
    }
  }

  return;
}

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

void SexyPedReco::GetSexyPedestals(RawDigitDataBlock *rdb) {

  MSG("SexyPedReco", Msg::kDebug) << "Filling the Sexy Pedestal TTree with RawDigits data" << endl;

  TIter rdit = rdb->GetDatumIter();
  RawDigit *rd;
  
  //Get validity context
  while ( (rd = (RawDigit*) rdit()) ) {
    RawChannelId rawid = rd->GetChannel();
    
    fcrate = rawid.GetCrate();
    fvarc  = rawid.GetVarcId();
    fvmm   = rawid.GetVmm();
    fvfb   = rawid.GetVaAdcSel();
    fva    = rawid.GetVaChip();
    fvach  = rawid.GetVaChannel();
    
    fadc   = rd->GetADC();
    ftdc   = rd->GetTDC()*1.5625e-9;
    ferror = rd->GetErrorCode();
    
    if(!fOnline) {
      PlexHandle ph(rdb->GetVldContext());
      fplane = ph.GetSEIdAltL(rawid).GetPlane();
    }
    else fplane = -1;

    //
    //Fill SexyPedestal TTree:
    //
    if(fStoreTree) fSexyPedTree->Fill();

    //
    // Perform Online calculations
    //
    if(fOnline){ //&& ferror&2==2) {

      if(!myPedAna) myPedAna = SexyPedAna::InstantiateMe();

      Int_t key = SexyPedestal::EncriptMeAKey(fcrate,fvarc,fvmm,fvfb,fva,fvach);
      
      //If it does not found (i.e. iter points to end) -> then create it:
      if( myPedAna->mySP.find(key) == myPedAna->mySP.end() ) myPedAna->mySP[key] = new SexyPedestal();

      myPedAna->frun = frunNum;      
      myPedAna->mySP[key]->AddressMe(fcrate,fvarc,fvmm,fvfb,fva,fvach);
      myPedAna->mySP[key]->SetPlane(fplane);
      myPedAna->mySP[key]->SetError(ferror);
      myPedAna->mySP[key]->InputSPEntry(fadc);

      if(ferror>2)
        MSG("SexyPedReco", Msg::kError) << "Found ErrorCode (!=2): " << ferror << endl;

    }
    //if(ferror>2)
    //MSG("SexyPedReco", Msg::kError) << "Found Errors on data: " << ferror << endl;
  }
  return;
}

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

void SexyPedReco::Config(const Registry& r) {

  Int_t tmpi;

  if(r.Get("GetPed-SparsTables", tmpi)) fGetTableTree  = tmpi;
  if(r.Get("StoreSexyTree", tmpi))      fStoreTree     = tmpi;
  if(r.Get("OnlineVersion", tmpi))      fOnline        = tmpi;

  return;
}

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

const Registry& SexyPedReco::DefaultConfig() const { 

  static Registry r;

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

  r.UnLockValues();
  r.Set("GetPed-SparsTables",1);
  r.Set("StoreSexyTree",1);
  r.Set("OnlineVersion",0);
  r.LockValues();

  return r;
}       


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

void SexyPedReco::EndJob() {

  //
  // Store Tree into File
  //
  if(fStoreTree) {

    fFile->Write();
    fFile->Close();
  }

  //
  // Perform Quick Analysis for Online purposes:
  //
  if(fOnline) myPedAna->DoOnlineAnalysis();

  return; 
}

//EOF

---