StTpcCalibrationMaker/StTpcFEEGainCalibMaker.cxx

// $Id: StTpcFEEGainCalibMaker.cxx,v 1.8 2007/04/28 17:57:17 perev Exp $
// $Log: StTpcFEEGainCalibMaker.cxx,v $
// Revision 1.8  2007/04/28 17:57:17  perev
// Redundant StChain.h removed
//
// Revision 1.7  2003/09/02 17:59:11  perev
// gcc 3.2 updates + WarnOff
//
// Revision 1.6  2003/07/18 18:31:49  perev
// test for nonexistance of XXXReader added
//
// Revision 1.5  2000/08/04 21:03:53  perev
// Leaks + Clear() cleanup
//
// Revision 1.4  1999/10/11 08:04:17  fretiere
// Fix bugg + add README, LOG and ID
//
//                                                                      //
//                   DEAD CHANNEL FINDER                                //
// For the time being, the finder look for corrupted pads :
// It makes histograms of the number of coorupted time bin as a function
// of the pad.
//                                                                      //
// C/C++
#include <Stiostream.h>
#include "Stiostream.h"
#include <stdio.h>
#include <unistd.h>
#include <errno.h>
#include <sys/types.h>
#include <sys/stat.h>
#include <fcntl.h>
#include <math.h>
#include <string.h>
// Root
#include "TFile.h"
#include "TH1.h"
#include "TH2.h"
#include "TH3.h"
#include "TNtuple.h"
// STAR
#include "St_DataSetIter.h"
//#include "StTpcDb/StTpcDb.h"
#include "StDAQMaker/StDAQReader.h"
#include "StDaqLib/TPC/fee_pin.h"
// Local
#include "StTpcFEEGainCalibMaker.h"
#include "StTpcCalibSetup.h"
#include "StTpcCalibSector.h"
// root again
ClassImp(StTpcFEEGainCalibMaker)
//  
//_____________________________________________________________________
// --- Constructor
//
StTpcFEEGainCalibMaker::StTpcFEEGainCalibMaker(const char *aMakerName,       
                               const StTpcCalibSetup* aSetup)
:
  StMaker(aMakerName),
  mSetup(aSetup)
{

}
//
//_____________________________________________________________________
// --- Destructor
//
StTpcFEEGainCalibMaker::~StTpcFEEGainCalibMaker(){
delete [] mTpcCalibSector;
}
//_____________________________________________________________________
// --- Init
//
Int_t StTpcFEEGainCalibMaker::Init(){

  //  cout << "StTpcFEEGainCalibMaker::Init()" << endl;
  //int tNumberOfSector=gStTpcDb->Dimensions()->numberOfSectors();
  int tNumberOfSector=24;

  static int tNPadAtRow[45]={
  88,96,104,112,118,126,134,142,150,158,166,174,182,
  98,100,102,104,106,106,108,110,112,112,114,116,118,120,122,122,
  124,126,128,128,130,132,134,136,138,138,140,142,144,144,144,144};
  ifstream* tBadFile = new ifstream(mSetup->getBadFileName());
  if(tBadFile->fail()) { 
    cout << "Couldn't open bad file " 
         << mSetup->getBadFileName() << endl
         << "Please run this macro with CalibType=1 before" << endl;
    return kStFatal;
  }
  ifstream* tDeadFile = new ifstream(mSetup->getDeadFileName());
  if(tDeadFile->fail()) { 
    cout << "Couldn't open bad file " 
         << mSetup->getDeadFileName()  << endl
         << "Please run this macro with CalibType=2 before" << endl;
    return kStFatal;
  }

  mTpcCalibSector = new StTpcCalibSector*[tNumberOfSector];
  for(int tiSect=0;
      tiSect<tNumberOfSector;
      tiSect++){
    mTpcCalibSector[tiSect]= new StTpcCalibSector(mSetup,(tiSect+1),
                                                  tNPadAtRow);
      mTpcCalibSector[tiSect]->readBadTable(tBadFile);
      mTpcCalibSector[tiSect]->readDeadTable(tDeadFile);
  }
  tBadFile->close();
  tDeadFile->close();
  delete tBadFile;

  return kStOK;
}
//
//_____________________________________________________________________
// --- Make calculation
//
Int_t StTpcFEEGainCalibMaker::Make(){
  //int tNumberOfSector=gStTpcDb->Dimensions()->numberOfSectors();
  int tNumberOfSector=24;
  //  cout << "StTpcFEEGainCalibMaker::Make()" << endl;
  static St_DataSet * tDataSet=GetDataSet("StDAQReader");
  //cout << "Got the data set " << tDataSet<< endl;
  static StDAQReader *tDAQReader=(StDAQReader*)(tDataSet->GetObject());
  //cout << "Got the daq reader " << tDAQReader<< endl;
  static StTPCReader *tZSupReader=tDAQReader->getTPCReader(); 
  //cout << "Got the zero suppressed reader " << tZSupReader<< endl;
  if(!tZSupReader) {cout << "Coudln't get the reader " << endl;return kStWarn;}
  for(int tiSect=0;
      tiSect< tNumberOfSector;
      tiSect++){
    mTpcCalibSector[tiSect]->updateGain(tZSupReader);
  }

  return kStOK;
}
//_____________________________________________________________________________
void StTpcFEEGainCalibMaker::PrintInfo(){
  printf("**************************************************************\n");
  printf("* $Id: StTpcFEEGainCalibMaker.cxx,v 1.8 2007/04/28 17:57:17 perev Exp $\n");
  printf("**************************************************************\n");
  if (Debug()) StMaker::PrintInfo();
}

void StTpcFEEGainCalibMaker::CalcElectronicConvertor(int** aPadToFeeConvertor, 
                                             int** aPadToRDOConvertor){

  // Eventually this function should be replaced by a database function
  for(int tiFee=0;tiFee<182;tiFee++){
    for(int tiPin=0;tiPin<32;tiPin++){
      if(row_vs_fee[tiFee][tiPin]!=0 && pad_vs_fee[tiFee][tiPin]!=0){
        aPadToFeeConvertor[(row_vs_fee[tiFee][tiPin]-1)]
                          [(pad_vs_fee[tiFee][tiPin]-1)]=(tiFee+1);
        aPadToRDOConvertor[(row_vs_fee[tiFee][tiPin]-1)]
                          [(pad_vs_fee[tiFee][tiPin]-1)]=
          rdo_vs_fee[tiFee][tiPin];
      }
    }
  }
}
//_____________________________________________________________________________
Int_t StTpcFEEGainCalibMaker::Finish(){
  cout << "StTpcFEEGainCalibMaker::Finish()" << endl;

  //int tNumberOfSector=gStTpcDb->Dimensions()->numberOfSectors();
  int tNumberOfSector=24;
  //  mHPercOfDead = new TH1D("HPerDead","Percentage of bad pads",24,0.5,24.5);
  TFile *tHFile = new TFile(mSetup->getRootOutFileName(),"recreate");
  ofstream* tCalibFile = new ofstream(mSetup->getGainCalibFileName());

  for(int tiSect=0;
      tiSect<tNumberOfSector;
      tiSect++){
    mTpcCalibSector[tiSect]->calcGainCoeficient();
    mTpcCalibSector[tiSect]->writeGainHisto();
    mTpcCalibSector[tiSect]->writeCalibCoefTable(tCalibFile);

  }
  //mHPercOfDead->Write();
  tHFile->Close();
  tCalibFile->close();
  return kStOK;
}

---