LoadCalDetFilterChannels.cxx

Go to the documentation of this file.

#include "CalDetDST/LoadCalDetFilterChannels.h"
#include "CalDetDST/ChannelsToFilter.h"

#include <iostream>

#include "Validity/VldContext.h"
#include "Validity/VldRange.h"
#include "Plex/PlexHandle.h"
#include "Plex/PlexStripEndId.h"
#include "DatabaseInterface/DbiCascader.h"
#include "DatabaseInterface/DbiTableProxyRegistry.h"
#include "DatabaseInterface/DbiWriter.h"
#include "DatabaseInterface/DbiResultPtr.h"
#include "DatabaseInterface/DbiValidityRec.h"
#include "Conventions/StripEnd.h"

void MakeFilterChannelList(const VldTimeStamp& vts, const VldTimeStamp& vte,
                           const vector<RawChannelId>& v, 
                           bool createtemp,  
                           Detector::Detector_t det,
                           int sim)
{
   if(createtemp){
//  Ask the singleton DbiTableProxyRegistry for the DbiCascader.
      DbiCascader& cascader 
         = DbiTableProxyRegistry::Instance().GetCascader();
      
//  Define the table.
      string tableDescr = "(SeqNo int, Encoded BIGINT UNSIGNED)";
      
//  Ask the cascader to find a database that will accept it.
      Int_t dbNoTemp = cascader.CreateTemporaryTable("CHANNELSTOFILTER",
                                                     tableDescr);
      
      if ( dbNoTemp < 0 ) {
         std::cout << "No database to will accept temporary tables." 
                   << std::endl;
         return;
      }
   }
   
   VldRange vr(det,sim,vts,vte,"MakeFilterChannelList");
   VldTimeStamp now;
   Dbi::Task task=0;
   int aggno=-1; 
   DbiWriter<ChannelsToFilter> writer(vr,aggno,task,now); 
   for(unsigned int i=0; i<v.size(); i++){
      ChannelsToFilter chan(v[i]);
      writer<<chan;
   }
   if(!writer.Close()){
      std::cout<<"\nMakeFilterChannelList(): "
          <<"Output error writing data!"<<std::endl;
   }

}


void FilterChannelsT1102(bool createtemp){

   // construct, by hand, the list of bad strip ends
   // use a plexhandle to convert to rawchannelid
   // fill a vector with the rcids and call makefilterchannellist

   const Detector::Detector_t det = Detector::kCalDet;
//   const int sim = (SimFlag::kData|SimFlag::kMC);
   const VldTimeStamp start(2002, 9, 16, 0, 0, 0);
   const VldTimeStamp end(2002, 10, 4, 0, 0, 0);


   const VldTimeStamp midrun(2002, 9, 24, 0,0,0);
   const VldContext vc(Detector::kCalDet,SimFlag::kData, midrun);
   PlexHandle ph(vc);
   
   vector<RawChannelId> v;
   
   // a complete loss 
   v.push_back(
      ph.GetRawChannelId(PlexStripEndId(det,38,10,StripEnd::kPositive)) );

   // the pedestal problem
   v.push_back(
      ph.GetRawChannelId(PlexStripEndId(det,20,4,StripEnd::kPositive)) );

   v.push_back(
      ph.GetRawChannelId(PlexStripEndId(det,40,4,StripEnd::kPositive)) );

   v.push_back(
      ph.GetRawChannelId(PlexStripEndId(det,21,19,StripEnd::kNegative)) );

   v.push_back(
      ph.GetRawChannelId(PlexStripEndId(det,41,19,StripEnd::kNegative)) );

   MakeFilterChannelList(start,end,v,createtemp);
   

}

void FilterChannelsT702(bool createtemp){

   // construct, by hand, the list of bad strip ends
   // use a plexhandle to convert to rawchannelid
   // fill a vector with the rcids and call makefilterchannellist

   const Detector::Detector_t det = Detector::kCalDet;
//   const int sim = (SimFlag::kData|SimFlag::kMC);
   const VldTimeStamp start(2002, 10, 4, 0, 0, 1);
   const VldTimeStamp end(2002, 10, 17, 0, 0, 1);


   const VldTimeStamp midrun(2002, 10, 10, 0,0,0);
   const VldContext vc(Detector::kCalDet,SimFlag::kData, midrun);
   PlexHandle ph(vc);
   
   vector<RawChannelId> v;
   
   // a complete loss 
   v.push_back(
      ph.GetRawChannelId(PlexStripEndId(det,38,10,StripEnd::kPositive)) );

   // the pedestal problem
   v.push_back(
      ph.GetRawChannelId(PlexStripEndId(det,20,4,StripEnd::kPositive)) );

   v.push_back(
      ph.GetRawChannelId(PlexStripEndId(det,40,4,StripEnd::kPositive)) );

   v.push_back(
      ph.GetRawChannelId(PlexStripEndId(det,21,19,StripEnd::kNegative)) );

   v.push_back(
      ph.GetRawChannelId(PlexStripEndId(det,41,19,StripEnd::kNegative)) );

   MakeFilterChannelList(start,end,v,createtemp);

}

void DumpFilterChannels(const VldContext& vc){
   PlexHandle ph(vc);
   DbiResultPtr<ChannelsToFilter> rptr(vc);   
   std::cout<<"List of Channels: "<<std::endl;
   for(unsigned int i=0; i<rptr.GetNumRows(); i++){
      const ChannelsToFilter* ctof = rptr.GetRow(i);
      std::cout<<(*ctof)<<std::endl;
      std::cout<<ph.GetSEIdAltL(ctof->GetRawChannelId())<<std::endl;
   }
   
}

---