PulserCollectorModule.cxx

Go to the documentation of this file.

// $Id: PulserCollectorModule.cxx,v 1.7 2003/02/06 20:02:58 rjn Exp $
//
// PulserCollectorModule.cxx
//
// Collects summaries, writes to database
//
// Author: Phil Adamson 5/2002
//
// $Log: PulserCollectorModule.cxx,v $
// Revision 1.7  2003/02/06 20:02:58  rjn
// Modified PulserSummary.h and added classes to make PulserSummary objects from Raw Data.
//
// Revision 1.6  2002/10/17 05:51:47  rhatcher
// comment out unused variables
//
// Revision 1.5  2002/10/07 10:07:43  phil.adamson
// Try to fix the occasional nan...
//
// Revision 1.4  2002/09/25 08:54:07  phil.adamson
// Add Ryan's bits to make PulserSummary objects from Raw Data
//
// Revision 1.3  2002/08/27 15:16:16  phil.adamson
// Change the PIN diode DB to use the component parts rather than the PIN Id.
//
// Revision 1.2  2002/08/07 16:16:51  phil.adamson
// Minor cosmetic details, plus provide ZeroCorr() methods in PulserSummary.
// Phil.
//
//

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

ClassImp(PulserCollectorModule)

  // .........

CVSID("$Id: PulserCollectorModule.cxx,v 1.7 2003/02/06 20:02:58 rjn Exp $");
JOBMODULE(PulserCollectorModule,"PulserCollectorModule","Collect summaries, write to DB");

PulserCollectorModule::PulserCollectorModule()
{
  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
}

PulserCollectorModule::~PulserCollectorModule()
{

}

JobCResult PulserCollectorModule::ReferencePointReco(MomNavigator *mom) 
{
   //Get PulserSummary...
   PulserSummaryList *psl = dynamic_cast<PulserSummaryList *>
      (mom->GetFragment("PulserSummaryList"));
   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;
   //unused: 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 PulserCollectorModule::Reco(MomNavigator *mom)
{
   if(fDoReferencePoint) {
      return ReferencePointReco(mom);
   }
   //Get PulserSummary...
   PulserSummaryList *psl = dynamic_cast<PulserSummaryList *>
      (mom->GetFragment("PulserSummaryList"));
   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()){
         // Got curve - now write DB
         WriteRawCurveData(psl);
      }
      
   }
   
   return JobCResult::kFailed;
}



void PulserCollectorModule::WriteRawCurveData(PulserSummaryList *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, "PulserCollectorModule");
   VldTimeStamp create = start; // now!
   create.Add(VldTimeStamp(fCodeVersion,0));
   int aggNo = psl->GetAggregateNo();
   DbiWriter<PulserRawGain> writer(vr,aggNo,0,create);
   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,create);
   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,create);
   w << fHighPinCurve;
   w << fLowPinCurve;
   w.Close();
   fNearMapCurve.clear();
   fFarMapCurve.clear();
   fIdent=0;
}

void PulserCollectorModule::WriteReferenceDriftData(PulserSummaryList *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, "PulserCollectorModule");
   VldTimeStamp create = start; 
   create.Add(VldTimeStamp(fCodeVersion,0));
   int aggNo = psl->GetAggregateNo();
   DbiWriter<PulserRefDrift> writer(vr,aggNo,0,create);
   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,create);
   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,create);
   w << fHighPinRef;
   w << fLowPinRef;
   w.Close();
   fNearMapRef.clear();
   fFarMapRef.clear();
}


void PulserCollectorModule::WriteRawDriftData(PulserSummaryList *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, "PulserCollectorModule");
   VldTimeStamp create = start; 
   create.Add(VldTimeStamp(fCodeVersion,0));
   int aggNo = psl->GetAggregateNo();
   DbiWriter<PulserRawDrift> writer(vr,aggNo,0,create);
   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,create);
   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,create);
   w << fHighPinDrift;
   w << fLowPinDrift;
   w.Close();
   fNearMapDrift.clear();
   fFarMapDrift.clear();
}


void PulserCollectorModule::HandleCommand(JobCommand *command) 
{

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

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

}

---