Fixing the decoding of regional header bits.

[u/mrichter/AliRoot.git] / ITS / AliITSCalibrationSSD.cxx

/**************************************************************************
 * Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
 *                                                                         *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/
/* $Id$  */

#include "AliITSCalibrationSSD.h"
//////////////////////////////////////////////////////
//  Calibration class for set:ITS                   //
//  Specific subdetector implementation             //
//  for silicon strips detectors                    //
//                                                  //
//                                                  //
//////////////////////////////////////////////////////

const Double_t AliITSCalibrationSSD::fgkNoiseNDefault = 625.;
const Double_t AliITSCalibrationSSD::fgkNoisePDefault = 420.;
const Int_t AliITSCalibrationSSD::fgkNParDefault = 6;
const Double_t AliITSCalibrationSSD::fgkSigmaPDefault = 3.;
const Double_t AliITSCalibrationSSD::fgkSigmaNDefault = 2.;

ClassImp(AliITSCalibrationSSD)

//______________________________________________________________________
AliITSCalibrationSSD::AliITSCalibrationSSD():
fNPar(0),
fDetPar(0),
fNoiseP(0),
fNoiseN(0),
fSigmaP(0),
fSigmaN(0),
fNoise(0),
fPedestal(),
fGain(0),
fBadChannels(0),
fIsBad(kFALSE),
fGainP(0),
fGainN(0),
fNoisP(0),
fNoisN(0),
fPedP(),
fPedN(),
fNoisePThreshold(0),
fNoisyPChannelsList(0),
fNoiseNThreshold(0),
fNoisyNChannelsList(0),
fDeadNChannelsList(0),
fDeadPChannelsList(0){
    // Default Constructor

    SetNoiseParam(fgkNoisePDefault,fgkNoiseNDefault);
    for(Int_t i=0;i<fgkChipsPerModule;i++){
      fIsChipBad[i]=kFALSE;
    }
}
//______________________________________________________________________
AliITSCalibrationSSD::AliITSCalibrationSSD(const char *dataType):
fNPar(0),
fDetPar(0),
fNoiseP(0),
fNoiseN(0),
fSigmaP(0),
fSigmaN(0),
fNoise(0),
fPedestal(0),
fGain(0),
fBadChannels(0),
fIsBad(kFALSE),
fGainP(0),
fGainN(0),
fNoisP(0),
fNoisN(0),
fPedP(0),
fPedN(0),
fNoisePThreshold(0),
fNoisyPChannelsList(0),
fNoiseNThreshold(0),
fNoisyNChannelsList(0),
fDeadNChannelsList(0),
fDeadPChannelsList(0){
    // constructor

    SetNoiseParam(fgkNoisePDefault,fgkNoiseNDefault);
    SetDataType(dataType);
    SetSigmaSpread(fgkSigmaPDefault,fgkSigmaNDefault);
    SetNDetParam(fgkNParDefault);   // Sets fNPar=6 by default.
    fDetPar = new Double_t[fNPar];
    if (fNPar==6) {
	fDetPar[0]=10.;
	fDetPar[1]=5.;
	fDetPar[2]=0.02;
	fDetPar[3]=0.02;
	fDetPar[4]=0.02;
	fDetPar[5]=0.03;
    } // end if
    for(Int_t i=0;i<fgkChipsPerModule;i++){
      fIsChipBad[i]=kFALSE;
    }
}
//______________________________________________________________________
AliITSCalibrationSSD::~AliITSCalibrationSSD(){
    // destructor
 
    delete [] fDetPar;
}
//______________________________________________________________________
void AliITSCalibrationSSD::SetDetParam(Double_t  *par){
    // set det param
    Int_t i;

    for (i=0; i<fNPar; i++) {
	fDetPar[i]=par[i];
	//printf("\n CompressPar %d %d \n",i,fCPar[i]);    
    } // end for i
}
//______________________________________________________________________
void AliITSCalibrationSSD::GetDetParam(Double_t  *par) const {
    // get det param
    Int_t i;

    for (i=0; i<fNPar; i++) {
	par[i]=fDetPar[i];
    } // end for i
}

//______________________________________________________________________
void AliITSCalibrationSSD::FillBadChipMap() {

  Int_t mc=0;
  Int_t cc[12];

  // P-side
  for(Int_t i=0; i<6; i++){
    cc[i]=0;
    for(Int_t j=0; j<ChannelsPerChip(); j++) {
      if(IsPChannelBad(i*ChannelsPerChip()+j)) cc[i]++;
    }
    if(cc[i]==ChannelsPerChip()) { SetChipBad(i); mc++; }
  }
  
  // N-side
  for(Int_t i=6; i<11; i++){
    cc[i]=0;
    for(Int_t j=0; j<ChannelsPerChip(); j++) {
      if(IsNChannelBad(1535-i*ChannelsPerChip()-j)) cc[i]++;      
    }
    if(cc[i]==ChannelsPerChip()) { SetChipBad(i); mc++; }
  }
  
  if(mc==ChannelsPerChip()) fIsBad=kTRUE;
}