/data/refman/pasoursint/CMSSW_5_3_3/src/CalibMuon/CSCCalibration/interface/CSCBadWiresConditions.h

Go to the documentation of this file.

#ifndef _CSCBADWIRESCONDITIONS_H
#define _CSCBADWIRESCONDITIONS_H

#include <memory>
#include <cmath>
#include "FWCore/Framework/interface/SourceFactory.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/ESProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/EventSetupRecordIntervalFinder.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include <DataFormats/MuonDetId/interface/CSCDetId.h>
#include "CondFormats/CSCObjects/interface/CSCBadWires.h"
#include "CondFormats/DataRecord/interface/CSCBadWiresRcd.h"

class CSCBadWiresConditions: public edm::ESProducer, public edm::EventSetupRecordIntervalFinder  {
 public:
  CSCBadWiresConditions(const edm::ParameterSet&);
  ~CSCBadWiresConditions();
  

  inline static CSCBadWires *  prefillBadWires();

  typedef const  CSCBadWires * ReturnType;
  
  ReturnType produceBadWires(const CSCBadWiresRcd&);
  
 private:
  // ----------member data ---------------------------
  void setIntervalFor(const edm::eventsetup::EventSetupRecordKey &, const edm::IOVSyncValue&, edm::ValidityInterval & );
  CSCBadWires *cndbBadWires ;

};

#include<fstream>
#include<vector>
#include<iostream>

// to workaround plugin library
inline CSCBadWires *  CSCBadWiresConditions::prefillBadWires()
{
  CSCBadWires * cndbbadwires = new CSCBadWires();

  int new_index,new_chan;
  int new_layer,new_channel, new_flag1,  new_flag2, new_flag3,new_pointer;
  std::vector<int> new_index_id;
  std::vector<short int> new_layer_id;
  std::vector<short int> new_chan_id;
  std::vector<int> new_badchannels;
  std::vector<short int> new_flag1_id;
  std::vector<short int> new_flag2_id;
  std::vector<short int> new_flag3_id;
  std::vector<int> new_cham_id;
  std::vector<int> new_point;

  //when bad channels exist use the following
   int new_nrlines1=0;
   int new_nrlines2=0;
  //when empty
  //int new_nrlines1=-1;
  //int new_nrlines2=-1;
 
  std::ifstream newdata1;
  std::ifstream newdata2;

  newdata1.open("badwires1.dat",std::ios::in); 
  if(!newdata1) {
    std::cerr <<"Error: badwires1.dat -> no such file!"<< std::endl;
    exit(1);
  }
  newdata2.open("badwires2.dat",std::ios::in); //channelcontainer
  if(!newdata2) {
    std::cerr <<"Error: badwires2.dat -> no such file!"<< std::endl;
    exit(1);
  }
  
  while (!newdata1.eof() ) { 
    newdata1 >> new_index >> new_pointer >> new_chan; 
    new_index_id.push_back(new_index);
    new_point.push_back(new_pointer);
     new_badchannels.push_back(new_chan);
    new_nrlines1++;
  }
  newdata1.close();
  
 while (!newdata2.eof() ) { 
    newdata2 >> new_layer >> new_channel >> new_flag1 >> new_flag2 >> new_flag3; 
    new_layer_id.push_back(new_layer);
    new_chan_id.push_back(new_channel);
    new_flag1_id.push_back(new_flag1);
    new_flag2_id.push_back(new_flag2);
    new_flag3_id.push_back(new_flag3);
    new_nrlines2++;
  }
  newdata2.close();

 
  CSCBadWires::BadChamberContainer & itemvector1 = cndbbadwires->chambers;
  itemvector1.resize(new_nrlines1);
  

  CSCBadWires::BadChannelContainer & itemvector2 = cndbbadwires->channels;
  itemvector2.resize(new_nrlines2);
  
  cndbbadwires->numberOfBadChannels = new_nrlines2;

  for(int i=0; i<new_nrlines1;i++){
    itemvector1[i].chamber_index = new_index_id[i];
    itemvector1[i].pointer = new_point[i];
    itemvector1[i].bad_channels = new_badchannels[i];
    
  }

  for(int j=0; j<new_nrlines2;++j){
       itemvector2[j].layer = new_layer_id[j];
       itemvector2[j].channel = new_chan_id[j];
       itemvector2[j].flag1 = new_flag1_id[j];
       itemvector2[j].flag2 = new_flag2_id[j];
       itemvector2[j].flag3 = new_flag3_id[j];
  }

  return cndbbadwires;
}

#endif