PulserFromRawCollectorModule.cxx

Go to the documentation of this file.

// $Id: PulserFromRawCollectorModule.cxx,v 1.1 2003/02/06 21:38:33 rjn Exp $
//
// PulserFromRawCollectorModule.cxx
//
// Collects summaries, writes to database
//
// Author: Ryan Nichol
//

#include "JobControl/JobCModuleRegistry.h"
#include "JobControl/JobCommand.h"
#include "PulserCalibration/PulserFromRawCollectorModule.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "PulserCalibration/PulserFromRawSummaryList.h"
#include "PulserCalibration/PulserSummary.h"
#include "Validity/VldContext.h"
#include "Validity/VldRange.h"
#include "Validity/VldTimeStamp.h"
#include "DatabaseInterface/DbiWriter.h"

ClassImp(PulserFromRawCollectorModule)

  // .........

CVSID("$Id: PulserFromRawCollectorModule.cxx,v 1.1 2003/02/06 21:38:33 rjn Exp $");
JOBMODULE(PulserFromRawCollectorModule,"PulserFromRawCollectorModule","Collect summaries, write to DB");

PulserFromRawCollectorModule::PulserFromRawCollectorModule()
{
  fDoReferencePoint=0;
  fIdent = 0; 
  fGCValid = 3600*24*62; // 2 months validity for Gaincurve
  fDPValid = 3600*24*1;  // Driftpoint valid up to a day?
  fLongTime = 3600*24*365*20; // 20 yrs
}

PulserFromRawCollectorModule::~PulserFromRawCollectorModule()
{

}

JobCResult PulserFromRawCollectorModule::ReferencePointReco(MomNavigator *mom) 
{
   //Get PulserSummary...
   PulserFromRawSummaryList *psl = dynamic_cast<PulserFromRawSummaryList *>
      (mom->GetFragment("PulserFromRawSummaryList"));
   if (psl==0) {
      MSG("Pulser",Msg::kVerbose) <<"No pulser summary list here\n";
      return JobCResult::kFailed;
   }
   MSG("Pulser",Msg::kDebug) <<"Got pulser summary "<<psl->GetPointIndex()
                             <<" of "<<psl->GetExpectedPoints()<<endl;
   // int point = psl->GetPointIndex();
   switch(psl->GetExpectedPoints()) {
   case 0:
      // Noise - no pulsing
      break;
  case 1:
     // Doug Point
     MSG("Pulser",Msg::kDebug) <<"Doug Point"<<endl;
     ;
     fHighPinRef.New(psl->GetAggregateNo(),
                     psl->GetHighPin().GetPinDiodeId(),
                     psl->GetHighPin());
     fLowPinRef.New(psl->GetAggregateNo(),
                    psl->GetLowPin().GetPinDiodeId(),
                    psl->GetLowPin());
     // Loop over near end data
     for (std::map<Int_t,PulserSummary>::const_iterator it=psl->GetNearBegin();
          it!=psl->GetNearEnd(); it++) {
        fNearMapRef[(*it).first].New(psl->GetAggregateNo(),
                                     (*it).second.GetStripEndId().BuildPlnStripEndKey(),
                                     (*it).second);
     }
     for (std::map<Int_t,PulserSummary>::const_iterator it=psl->GetFarBegin();
          it!=psl->GetFarEnd(); it++) {
        fFarMapRef[(*it).first].New(psl->GetAggregateNo(),
                                    (*it).second.GetStripEndId().BuildPlnStripEndKey(),
                                    (*it).second);
     }
     WriteReferenceDriftData(psl);
     break;
   default:
      return JobCResult::kFailed;
   }
   return JobCResult::kPassed;
}



JobCResult PulserFromRawCollectorModule::Reco(MomNavigator *mom)
{
   if(fDoReferencePoint) {
      return ReferencePointReco(mom);
   }
   //Get PulserSummary...
   PulserFromRawSummaryList *psl = dynamic_cast<PulserFromRawSummaryList *>
      (mom->GetFragment("PulserFromRawSummaryList"));
   if (psl==0) {
      MSG("Pulser",Msg::kVerbose) <<"No pulser summary list here\n";
      return JobCResult::kFailed;
   }
   MSG("Pulser",Msg::kDebug) <<"Got pulser summary "<<psl->GetPointIndex()
                             <<" of "<<psl->GetExpectedPoints()<<endl;
   int point = psl->GetPointIndex();
   switch(psl->GetExpectedPoints()) {
   case 0:
      // Noise - no pulsing
      break;
   case 1:
      // Doug Point
      MSG("Pulser",Msg::kDebug) <<"Doug Point"<<endl;
      fHighPinDrift.New(psl->GetAggregateNo(),
                        psl->GetHighPin().GetPinDiodeId(),
                        psl->GetHighPin());
      fLowPinDrift.New(psl->GetAggregateNo(),
                       psl->GetLowPin().GetPinDiodeId(),
                       psl->GetLowPin());
      // Loop over near end data
      for (std::map<Int_t,PulserSummary>::const_iterator it=psl->GetNearBegin();
           it!=psl->GetNearEnd(); it++) {
         fNearMapDrift[(*it).first].New(psl->GetAggregateNo(),
                                        (*it).second.GetStripEndId().BuildPlnStripEndKey(),
                                        (*it).second);
      }
      for (std::map<Int_t,PulserSummary>::const_iterator it=psl->GetFarBegin();
           it!=psl->GetFarEnd(); it++) {
         fFarMapDrift[(*it).first].New(psl->GetAggregateNo(),
                                       (*it).second.GetStripEndId().BuildPlnStripEndKey(),
                                       (*it).second);
      }
      WriteRawDriftData(psl);
      break;
   default:
      // Gain Curve
      // Check we're on the right curve...
      if (fIdent!=0 && fIdent!=psl->GetIdent()) {
         MSG("Pulser",Msg::kWarning)<< "Conflicting gaincurves in same place!\n";
      }
      if (fIdent==0) {
         fHighPinCurve.New(psl->GetAggregateNo(),psl->GetExpectedPoints(),
                           psl->GetHighPin().GetPinDiodeId());
         fLowPinCurve.New(psl->GetAggregateNo(),psl->GetExpectedPoints(),
                          psl->GetLowPin().GetPinDiodeId());
         fIdent = psl->GetIdent();
      }
      MSG("Pulser",Msg::kDebug) <<"Gain Curve"<<endl;
      // Loop over near end data
      for (std::map<Int_t,PulserSummary>::const_iterator it=psl->GetNearBegin();
           it!=psl->GetNearEnd(); it++) {
         if (!fNearMapCurve.count((*it).first)) {
            fNearMapCurve[(*it).first].New(psl->GetAggregateNo(),
                                           psl->GetExpectedPoints(),
                                           (*it).first);
         }
         fNearMapCurve[(*it).first].AddPoint(point,(*it).second);
    }
      for (std::map<Int_t,PulserSummary>::const_iterator it=psl->GetFarBegin();
           it!=psl->GetFarEnd(); it++) {
         if (!fFarMapCurve.count((*it).first))    {
            fFarMapCurve[(*it).first].New(psl->GetAggregateNo(),
                                          psl->GetExpectedPoints(),
                                          (*it).first);
         }
         fFarMapCurve[(*it).first].AddPoint(point,(*it).second);
      }
      fHighPinCurve.AddPoint(point,psl->GetHighPin());
      fLowPinCurve.AddPoint(point,psl->GetLowPin());
      if (point == psl->GetExpectedPoints()||point==37){
         // Got curve - now write DB
         WriteRawCurveData(psl);
      }
      
   }
   return JobCResult::kPassed;

}



void PulserFromRawCollectorModule::WriteRawCurveData(PulserFromRawSummaryList *psl)
{
   VldContext vld = psl->GetVldContext();
   VldTimeStamp start = vld.GetTimeStamp();
   VldTimeStamp end = start;
   end.Add(VldTimeStamp(fGCValid,0));
   VldRange vr(vld.GetDetector(),vld.GetSimFlag(),start,
               end, "PulserFromRawCollectorModule");
   VldTimeStamp create; // now!
   int aggNo = psl->GetAggregateNo();
   DbiWriter<PulserRawGain> writer(vr,aggNo,0,start);
   MSG("Pulser",Msg::kDebug)<<"Writing database "<<psl->GetVldContext()<<endl;
   for (std::map <Int_t,PulserRawGain>::const_iterator 
           it = fNearMapCurve.begin(); it!=fNearMapCurve.end(); it++) {
      writer << (*it).second;
   }
   writer.Close();
   DbiWriter<PulserRawGainF> fwriter(vr,aggNo,0,start);
   for (std::map <Int_t,PulserRawGainF>::const_iterator 
           it = fFarMapCurve.begin(); it!=fFarMapCurve.end(); it++) {
      fwriter << (*it).second;
   }
   fwriter.Close();
   DbiWriter<PulserRawGainPin> w(vr,aggNo,0,start);
   w << fHighPinCurve;
   w << fLowPinCurve;
   w.Close();
   fNearMapCurve.clear();
   fFarMapCurve.clear();
   fIdent=0;
}

void PulserFromRawCollectorModule::WriteReferenceDriftData(PulserFromRawSummaryList *psl)
{
   VldContext vld = psl->GetVldContext();
   //   VldTimeStamp start = vld.GetTimeStamp();
   //VldTimeStamp start = vld.GetTimeStamp();
   //VldTimeStamp end = start;
   //end.Add(VldTimeStamp(fDPValid,0));
   VldTimeStamp start = VldTimeStamp(2002,9,4,15,0,1); //Valid For all time
   VldTimeStamp end = VldTimeStamp(2002,10,2,0,0,0);
   //end.Add(VldTimeStamp(fLongTime,0));
   VldRange vr(vld.GetDetector(),vld.GetSimFlag(),start,
               end, "PulserFromRawCollectorModule");
   VldTimeStamp create; // now!
   int aggNo = psl->GetAggregateNo();
   DbiWriter<PulserRefDrift> writer(vr,aggNo,0,start);
   MSG("Pulser",Msg::kDebug)<<"Writing database "<<psl->GetVldContext()<<endl;
   for (std::map <Int_t,PulserRefDrift>::const_iterator 
           it = fNearMapRef.begin(); it!=fNearMapRef.end(); it++) {
      //MSG("Pulser",Msg::kWarning)<<"RefDrift: " << (*it).first << endl;
      writer << (*it).second;
   }
   writer.Close();
   DbiWriter<PulserRefDriftF> fwriter(vr,aggNo,0,start);
   for (std::map <Int_t,PulserRefDriftF>::const_iterator 
           it = fFarMapRef.begin(); it!=fFarMapRef.end(); it++) {
      fwriter << (*it).second;
   }
   fwriter.Close();
   DbiWriter<PulserRefDriftPin> w(vr,aggNo,0,start);
   w << fHighPinRef;
   w << fLowPinRef;
   w.Close();
   fNearMapRef.clear();
   fFarMapRef.clear();
}


void PulserFromRawCollectorModule::WriteRawDriftData(PulserFromRawSummaryList *psl)
{
   VldContext vld = psl->GetVldContext();
   VldTimeStamp start = vld.GetTimeStamp();
   VldTimeStamp end = start;
   end.Add(VldTimeStamp(fDPValid,0));
   VldRange vr(vld.GetDetector(),vld.GetSimFlag(),start,
               end, "PulserFromRawCollectorModule");
   VldTimeStamp create; // now!
   int aggNo = psl->GetAggregateNo();
   DbiWriter<PulserRawDrift> writer(vr,aggNo,0,start);
   MSG("Pulser",Msg::kDebug)<<"Writing database "<<psl->GetVldContext()<<endl;
   for (std::map <Int_t,PulserRawDrift>::const_iterator 
           it = fNearMapDrift.begin(); it!=fNearMapDrift.end(); it++) {
      writer << (*it).second;
   }
   writer.Close();
   DbiWriter<PulserRawDriftF> fwriter(vr,aggNo,0,start);
   for (std::map <Int_t,PulserRawDriftF>::const_iterator 
           it = fFarMapDrift.begin(); it!=fFarMapDrift.end(); it++) {
      fwriter << (*it).second;
   }
   fwriter.Close();
   DbiWriter<PulserRawDriftPin> w(vr,aggNo,0,start);
   w << fHighPinDrift;
   w << fLowPinDrift;
   w.Close();
   fNearMapDrift.clear();
   fFarMapDrift.clear();
}


void PulserFromRawCollectorModule::HandleCommand(JobCommand *command) 
{

// Process configuration commands
  TString cmd = command->PopCmd();

  MSG("Pulser", Msg::kDebug)
                   << "PulserFromRawCollectorModule:HandleCommand: " << cmd << endl;
  if (cmd == "Reference") {
     fDoReferencePoint=1;
  }
  else {
     MSG("Pulser", Msg::kWarning)
                          << "Unrecognized Command: " << cmd << endl;
  }

}

---