[u/mrichter/AliRoot.git] / FMD / AliFMDCalibFaker.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, 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$ */
/** @file    AliFMDCalibFaker.cxx
    @author  Christian Holm Christensen <cholm@nbi.dk>
    @date    Sun Mar 26 18:29:21 2006
    @brief   Make fake calibration data 
    @ingroup FMD_util
*/
//____________________________________________________________________
//                                                                          
// Forward Multiplicity Detector based on Silicon wafers. 
//
// This task creates fake calibrations. Which calibration, depends on
// the bit mask passed to the constructor, or added by `AddCalib'.
//
// The default is to write all calibration parameters to a local
// storage `local://$ALICE_ROOT/OCDB' which is were the sources live (sigh!
// - why oh why do we need to shit where we eat - it's just not
// healty).
//                                                       
#include "AliFMDDebug.h"		// ALIFMDDEBUG_H ALILOG_H
#include "AliFMDCalibFaker.h"      // ALIFMDCALIBFAKER_H
#include "AliFMDCalibGain.h"       // ALIFMDCALIBGAIN_H
#include "AliFMDCalibPedestal.h"   // ALIFMDCALIBPEDESTAL_H
#include "AliFMDCalibSampleRate.h" // ALIFMDCALIBPEDESTAL_H
#include "AliFMDAltroMapping.h"    // ALIFMDALTROMAPPING_H
#include "AliFMDCalibStripRange.h" // ALIFMDCALIBSTRIPRANGE_H
#include <AliCDBManager.h>         // ALICDBMANAGER_H
#include <AliCDBEntry.h>           // ALICDBMANAGER_H
//#include <Riostream.h>
#include <TSystem.h>
// #include <TMath.h>
#include <TROOT.h>
#include <TRandom.h>
#include <TF1.h>

//====================================================================
ClassImp(AliFMDCalibFaker)
#if 0
  ; // This is here to keep Emacs for indenting the next line
#endif

//____________________________________________________________________
AliFMDCalibFaker::AliFMDCalibFaker(Int_t mask, const char* loc) 
  : TTask("FMDCalibFaker", loc),
    fMask(mask),
    fGain(-1),
    fThresholdFactor(.1),
    fThreshold(-1),
    fPedestalMin(20),
    fPedestalMax(30), 
    fDeadChance(0),
    fRate(4),
    fZeroThreshold(0),
    fRunMin(0),
    fRunMax(10),
    fStripMin(0),
    fStripMax(127)
{
  // Default constructor 
}


#define MAKE_META(meta) \
  do { \
    meta = new AliCDBMetaData; \
    meta->SetResponsible(gSystem->GetUserInfo()->fRealName.Data()); \
    meta->SetAliRootVersion(gROOT->GetVersion()); \
    meta->SetBeamPeriod(1); \
    meta->SetComment("Dummy data for testing"); } while (false);
 

//__________________________________________________________________
void
AliFMDCalibFaker::Exec(Option_t*)
{
  // Make the objects. 
  AliCDBManager*     cdb      = AliCDBManager::Instance();
  AliFMDParameters*  param    = AliFMDParameters::Instance();
  Float_t            maxADC   = 1.F*param->GetAltroChannelSize();
  TObjArray          cleanup;

  if (GetTitle() && GetTitle()[0] != '\0') { 
    AliInfo(Form("Setting default storage to '%s'", GetTitle()));
    cdb->SetDefaultStorage(GetTitle());
  }
  
    
  AliCDBMetaData* meta = 0;
  if (TESTBIT(fMask, kPulseGain)) {
    // Info("Exec","Default gain to %f = %d * %f / %d", 
    //      (param->GetVA1MipRange() * param->GetEdepMip() / maxADC),
    //      param->GetVA1MipRange(), param->GetEdepMip(), Int_t(maxADC));
    if (fGain <= 0) {
      fGain = (param->GetVA1MipRange() * param->GetEdepMip() / maxADC);
    }
    fThreshold = fThresholdFactor * param->GetEdepMip();
    AliFMDCalibGain* gain = MakePulseGain();
    AliCDBId         id(AliFMDParameters::PulseGainPath(), fRunMin, fRunMax);
    MAKE_META(meta);
    meta->SetProperty("key1", gain);
    cdb->Put(gain, id, meta);
    cleanup.Add(gain);
    cleanup.Add(meta);
  }  
  if (TESTBIT(fMask, kPedestal)) {
    fPedestalMin = TMath::Max(TMath::Min(fPedestalMin, maxADC), 0.F);
    fPedestalMax = TMath::Max(TMath::Min(fPedestalMax, maxADC), fPedestalMin);
    AliFMDCalibPedestal* pedestal = MakePedestal();
    AliCDBId             id(AliFMDParameters::PedestalPath(),fRunMin,fRunMax);
    MAKE_META(meta);
    meta->SetProperty("key1", pedestal);
    cdb->Put(pedestal, id, meta);
    cleanup.Add(pedestal);
    cleanup.Add(meta);
  }
  if (TESTBIT(fMask, kDeadMap)) {
    fDeadChance = TMath::Max(TMath::Min(fDeadChance, 1.F), 0.F);
    AliFMDCalibDeadMap* deadMap = MakeDeadMap();
    AliCDBId            id(AliFMDParameters::DeadPath(), fRunMin, fRunMax);
    MAKE_META(meta);
    meta->SetProperty("key1", deadMap);
    cdb->Put(deadMap, id, meta);
    cleanup.Add(deadMap);
    cleanup.Add(meta);
  }
  if (TESTBIT(fMask, kZeroSuppression)) {
    fZeroThreshold = TMath::Min(fZeroThreshold, UShort_t(maxADC));
    AliFMDCalibZeroSuppression* zeroSup = MakeZeroSuppression();
    AliCDBId                    id(AliFMDParameters::ZeroSuppressionPath(), 
				   fRunMin, fRunMax);
    MAKE_META(meta);
    meta->SetProperty("key1", zeroSup);
    cdb->Put(zeroSup, id, meta);
    cleanup.Add(zeroSup);
    cleanup.Add(meta);
  }
  if (TESTBIT(fMask, kSampleRate)) {
    fRate = TMath::Max(TMath::Min(fRate, UShort_t(8)), UShort_t(1));
    AliFMDCalibSampleRate* rate = MakeSampleRate();
    AliCDBId               id(AliFMDParameters::SampleRatePath(),
			      fRunMin,fRunMax);
    MAKE_META(meta);
    meta->SetProperty("key1", rate);
    cdb->Put(rate, id, meta);
    cleanup.Add(rate);
    cleanup.Add(meta);
  }
  if (TESTBIT(fMask, kStripRange)) {
    fRate = TMath::Max(TMath::Min(fRate, UShort_t(8)), UShort_t(1));
    AliFMDCalibStripRange* range = MakeStripRange();
    AliCDBId               id(AliFMDParameters::StripRangePath(),
			      fRunMin,fRunMax);
    MAKE_META(meta);
    meta->SetProperty("key1", range);
    cdb->Put(range, id, meta);
    cleanup.Add(range);
    cleanup.Add(meta);
  }
  if (TESTBIT(fMask, kAltroMap)) {
    AliFMDAltroMapping* altroMap = MakeAltroMap();
    AliCDBId            id(AliFMDParameters::AltroMapPath(), fRunMin, fRunMax);
    MAKE_META(meta);
    meta->SetProperty("key1", altroMap);
    cdb->Put(altroMap, id, meta);
    cleanup.Add(altroMap);
    cleanup.Add(meta);
  }
  cdb->Destroy();
  cleanup.Delete();
}


//__________________________________________________________________
AliFMDCalibGain*
AliFMDCalibFaker::MakePulseGain() const
{
  // Make the actual data
  AliFMDCalibGain*  gain  = new AliFMDCalibGain;
  // Set threshold 
  gain->Set(fThreshold);
  for (UShort_t det = 1; det <= 3; det++) {
    Char_t rings[] = { 'I', (det == 1 ? '\0' : 'O'), '\0' };
    for (Char_t* ring = rings; *ring != '\0'; ring++) {
      UShort_t nSec = ( *ring == 'I' ? 20  :  40 );
      UShort_t nStr = ( *ring == 'I' ? 512 : 256 );
      for (UShort_t sec = 0; sec < nSec; sec++) {
        for (UShort_t str = 0; str < nStr; str++) {
          gain->Set(det, *ring, sec, str,
                    gRandom->Gaus(fGain, .01 * fGain));
        }
      }
    }
  }
  return gain;
}

//__________________________________________________________________
Float_t 
AliFMDCalibFaker::MakeNoise(Char_t ring, UShort_t str) const
{
  const UShort_t innerN    = 512;
  const UShort_t outerN    = 256;
  const UShort_t innerCut  = 350;
  const UShort_t outerCut  = 190;
  const Float_t  innerBase =   1.2;
  const Float_t  outerBase =   2.1;
  const Float_t  innerInc  =   0.5;
  const Float_t  outerInc  =   0.8;
  Float_t cut, base, inc, n;
  switch (ring) {
  case 'I': 
    cut = innerCut; base = innerBase; inc = innerInc; n = innerN; break;
  case 'O': 
    cut = outerCut; base = outerBase; inc = outerInc; n = outerN; break;
  default:
    return -1;
  }
  Float_t bare = base + (str < cut ? 
			 str / cut * inc : 
			 inc  - (str - cut) / (n - cut) * inc);
  return bare + gRandom->Uniform(-.07, .07);
}
  
//__________________________________________________________________
AliFMDCalibPedestal*
AliFMDCalibFaker::MakePedestal() const
{
  // Make the actual data
  AliFMDCalibPedestal*  pedestal  = new AliFMDCalibPedestal;
  
  for (UShort_t det = 1; det <= 3; det++) {
    Char_t rings[] = { 'I', 'O', '\0' };
    for (Char_t* ring = rings; *ring != '\0'; ring++) {
      if (*ring == 'O' && det == 1) continue;
      UShort_t nSec = ( *ring == 'I' ? 20  :  40 );
      UShort_t nStr = ( *ring == 'I' ? 512 : 256 );
      for (UShort_t sec = 0; sec < nSec; sec++) {
        for (UShort_t str = 0; str < nStr; str++) {
	  Float_t noise = MakeNoise(*ring, str);
	  Float_t ped   = gRandom->Uniform(fPedestalMin, fPedestalMax);
          pedestal->Set(det, *ring, sec, str, ped, noise);
	}
      }
    }
  }
  return pedestal;
}

//__________________________________________________________________
AliFMDCalibDeadMap*
AliFMDCalibFaker::MakeDeadMap() const
{
  // Make the actual data
  AliFMDCalibDeadMap*  deadmap  = new AliFMDCalibDeadMap;
  for (UShort_t det = 1; det <= 3; det++) {
    Char_t rings[] = { 'I', (det == 1 ? '\0' : 'O'), '\0' };
    for (Char_t* ring = rings; *ring != '\0'; ring++) {
      UShort_t nSec = ( *ring == 'I' ? 20  :  40 );
      UShort_t nStr = ( *ring == 'I' ? 512 : 256 );
      for (UShort_t sec = 0; sec < nSec; sec++) {
        for (UShort_t str = 0; str < nStr; str++) {
          deadmap->operator()(det, *ring, sec, str) = 
	    gRandom->Uniform(0, 1) < fDeadChance;
        }
      }
    }
  }
  return deadmap;
}

//__________________________________________________________________
AliFMDCalibZeroSuppression*
AliFMDCalibFaker::MakeZeroSuppression() const
{
  // Make the actual data
  AliFMDCalibZeroSuppression*  zs  = new AliFMDCalibZeroSuppression;
  for (UShort_t det = 1; det <= 3; det++) {
    Char_t rings[] = { 'I', (det == 1 ? '\0' : 'O'), '\0' };
    for (Char_t* ring = rings; *ring != '\0'; ring++) {
      UShort_t nSec = ( *ring == 'I' ? 20  :  40 );
      UShort_t nStr = ( *ring == 'I' ? 512 : 256 );
      for (UShort_t sec = 0; sec < nSec; sec++) {
        for (UShort_t str = 0; str < nStr; str++) {
          zs->operator()(det, *ring, sec, str) =  fZeroThreshold;
        }
      }
    }
  }
  return zs;
}

//__________________________________________________________________
AliFMDCalibSampleRate*
AliFMDCalibFaker::MakeSampleRate() const
{
  // Make sample rates 
  AliFMDCalibSampleRate*  sampleRate  = new AliFMDCalibSampleRate;
  for (UShort_t det = 1; det <= 3; det++) {
    Char_t rings[] = { 'I', (det == 1 ? '\0' : 'O'), '\0' };
    for (Char_t* ring = rings; *ring != '\0'; ring++) {
      UShort_t nSec = ( *ring == 'I' ? 20  :  40 );
      for (UShort_t sec = 0; sec < nSec; sec++) {
	sampleRate->Set(det, *ring, sec, 0, fRate);
      }
    }
  }
  return sampleRate;
}

//__________________________________________________________________
AliFMDCalibStripRange*
AliFMDCalibFaker::MakeStripRange() const
{
  // Make strip ranges 
  AliFMDCalibStripRange*  striprange  = new AliFMDCalibStripRange;
  for (UShort_t det = 1; det <= 3; det++) {
    Char_t rings[] = { 'I', (det == 1 ? '\0' : 'O'), '\0' };
    for (Char_t* ring = rings; *ring != '\0'; ring++) {
      UShort_t nSec = ( *ring == 'I' ? 20  :  40 );
      for (UShort_t sec = 0; sec < nSec; sec++) {
	striprange->Set(det, *ring, sec, 0, fStripMin, fStripMax);
      }
    }
  }
  return striprange;
}

//__________________________________________________________________
AliFMDAltroMapping*
AliFMDCalibFaker::MakeAltroMap() const
{
  // Make hardware mapping 
  AliFMDAltroMapping*  m  = new AliFMDAltroMapping;
  return m;
}
  
  
//____________________________________________________________________
//
// EOF
//
Commit	Line	Data
20345ac5	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
20345ac5	15	/* $Id$ */
c2fc1258	16	/** @file AliFMDCalibFaker.cxx
	17	@author Christian Holm Christensen <cholm@nbi.dk>
	18	@date Sun Mar 26 18:29:21 2006
	19	@brief Make fake calibration data
02a27b50	20	@ingroup FMD_util
c2fc1258	21	*/
20345ac5	22	//____________________________________________________________________
	23	//
	24	// Forward Multiplicity Detector based on Silicon wafers.
	25	//
	26	// This task creates fake calibrations. Which calibration, depends on
	27	// the bit mask passed to the constructor, or added by `AddCalib'.
	28	//
	29	// The default is to write all calibration parameters to a local
162637e4	30	// storage `local://$ALICE_ROOT/OCDB' which is were the sources live (sigh!
6169f936	31	// - why oh why do we need to shit where we eat - it's just not
6169f936	32	// healty).
20345ac5	33	//
f95a63c4	34	#include "AliFMDDebug.h" // ALIFMDDEBUG_H ALILOG_H
20345ac5	35	#include "AliFMDCalibFaker.h" // ALIFMDCALIBFAKER_H
	36	#include "AliFMDCalibGain.h" // ALIFMDCALIBGAIN_H
	37	#include "AliFMDCalibPedestal.h" // ALIFMDCALIBPEDESTAL_H
	38	#include "AliFMDCalibSampleRate.h" // ALIFMDCALIBPEDESTAL_H
	39	#include "AliFMDAltroMapping.h" // ALIFMDALTROMAPPING_H
c2fc1258	40	#include "AliFMDCalibStripRange.h" // ALIFMDCALIBSTRIPRANGE_H
20345ac5	41	#include <AliCDBManager.h> // ALICDBMANAGER_H
20345ac5	42	#include <AliCDBEntry.h> // ALICDBMANAGER_H
02a27b50	43	//#include <Riostream.h>
20345ac5	44	#include <TSystem.h>
02a27b50	45	// #include <TMath.h>
1e8f773e	46	#include <TROOT.h>
20345ac5	47	#include <TRandom.h>
a9579262	48	#include <TF1.h>
20345ac5	49
	50	//====================================================================
	51	ClassImp(AliFMDCalibFaker)
	52	#if 0
	53	; // This is here to keep Emacs for indenting the next line
	54	#endif
	55
	56	//____________________________________________________________________
	57	AliFMDCalibFaker::AliFMDCalibFaker(Int_t mask, const char* loc)
	58	: TTask("FMDCalibFaker", loc),
	59	fMask(mask),
	60	fGain(-1),
	61	fThresholdFactor(.1),
	62	fThreshold(-1),
	63	fPedestalMin(20),
	64	fPedestalMax(30),
	65	fDeadChance(0),
2aeec17d	66	fRate(4),
20345ac5	67	fZeroThreshold(0),
20345ac5	68	fRunMin(0),
c2fc1258	69	fRunMax(10),
	70	fStripMin(0),
	71	fStripMax(127)
20345ac5	72	{
	73	// Default constructor
	74	}
	75
	76
	77	#define MAKE_META(meta) \
	78	do { \
	79	meta = new AliCDBMetaData; \
	80	meta->SetResponsible(gSystem->GetUserInfo()->fRealName.Data()); \
	81	meta->SetAliRootVersion(gROOT->GetVersion()); \
	82	meta->SetBeamPeriod(1); \
	83	meta->SetComment("Dummy data for testing"); } while (false);
	84
	85
	86	//__________________________________________________________________
	87	void
	88	AliFMDCalibFaker::Exec(Option_t*)
	89	{
	90	// Make the objects.
	91	AliCDBManager* cdb = AliCDBManager::Instance();
	92	AliFMDParameters* param = AliFMDParameters::Instance();
8ec606c2	93	Float_t maxADC = 1.F*param->GetAltroChannelSize();
20345ac5	94	TObjArray cleanup;
20345ac5	95
c2fc1258	96	if (GetTitle() && GetTitle()[0] != '\0') {
	97	AliInfo(Form("Setting default storage to '%s'", GetTitle()));
	98	cdb->SetDefaultStorage(GetTitle());
	99	}
	100
20345ac5	101
	102	AliCDBMetaData* meta = 0;
	103	if (TESTBIT(fMask, kPulseGain)) {
8ec606c2	104	// Info("Exec","Default gain to %f = %d * %f / %d",
	105	// (param->GetVA1MipRange() * param->GetEdepMip() / maxADC),
	106	// param->GetVA1MipRange(), param->GetEdepMip(), Int_t(maxADC));
	107	if (fGain <= 0) {
	108	fGain = (param->GetVA1MipRange() * param->GetEdepMip() / maxADC);
	109	}
20345ac5	110	fThreshold = fThresholdFactor * param->GetEdepMip();
20345ac5	111	AliFMDCalibGain* gain = MakePulseGain();
02a27b50	112	AliCDBId id(AliFMDParameters::PulseGainPath(), fRunMin, fRunMax);
20345ac5	113	MAKE_META(meta);
	114	meta->SetProperty("key1", gain);
	115	cdb->Put(gain, id, meta);
	116	cleanup.Add(gain);
	117	cleanup.Add(meta);
	118	}
	119	if (TESTBIT(fMask, kPedestal)) {
	120	fPedestalMin = TMath::Max(TMath::Min(fPedestalMin, maxADC), 0.F);
	121	fPedestalMax = TMath::Max(TMath::Min(fPedestalMax, maxADC), fPedestalMin);
	122	AliFMDCalibPedestal* pedestal = MakePedestal();
8ec606c2	123	AliCDBId id(AliFMDParameters::PedestalPath(),fRunMin,fRunMax);
20345ac5	124	MAKE_META(meta);
	125	meta->SetProperty("key1", pedestal);
	126	cdb->Put(pedestal, id, meta);
	127	cleanup.Add(pedestal);
	128	cleanup.Add(meta);
	129	}
	130	if (TESTBIT(fMask, kDeadMap)) {
	131	fDeadChance = TMath::Max(TMath::Min(fDeadChance, 1.F), 0.F);
	132	AliFMDCalibDeadMap* deadMap = MakeDeadMap();
02a27b50	133	AliCDBId id(AliFMDParameters::DeadPath(), fRunMin, fRunMax);
20345ac5	134	MAKE_META(meta);
	135	meta->SetProperty("key1", deadMap);
	136	cdb->Put(deadMap, id, meta);
	137	cleanup.Add(deadMap);
	138	cleanup.Add(meta);
	139	}
	140	if (TESTBIT(fMask, kZeroSuppression)) {
	141	fZeroThreshold = TMath::Min(fZeroThreshold, UShort_t(maxADC));
	142	AliFMDCalibZeroSuppression* zeroSup = MakeZeroSuppression();
02a27b50	143	AliCDBId id(AliFMDParameters::ZeroSuppressionPath(),
20345ac5	144	fRunMin, fRunMax);
	145	MAKE_META(meta);
	146	meta->SetProperty("key1", zeroSup);
	147	cdb->Put(zeroSup, id, meta);
	148	cleanup.Add(zeroSup);
	149	cleanup.Add(meta);
	150	}
	151	if (TESTBIT(fMask, kSampleRate)) {
	152	fRate = TMath::Max(TMath::Min(fRate, UShort_t(8)), UShort_t(1));
	153	AliFMDCalibSampleRate* rate = MakeSampleRate();
8ec606c2	154	AliCDBId id(AliFMDParameters::SampleRatePath(),
8ec606c2	155	fRunMin,fRunMax);
20345ac5	156	MAKE_META(meta);
	157	meta->SetProperty("key1", rate);
	158	cdb->Put(rate, id, meta);
	159	cleanup.Add(rate);
	160	cleanup.Add(meta);
	161	}
c2fc1258	162	if (TESTBIT(fMask, kStripRange)) {
	163	fRate = TMath::Max(TMath::Min(fRate, UShort_t(8)), UShort_t(1));
	164	AliFMDCalibStripRange* range = MakeStripRange();
8ec606c2	165	AliCDBId id(AliFMDParameters::StripRangePath(),
8ec606c2	166	fRunMin,fRunMax);
c2fc1258	167	MAKE_META(meta);
	168	meta->SetProperty("key1", range);
	169	cdb->Put(range, id, meta);
	170	cleanup.Add(range);
	171	cleanup.Add(meta);
	172	}
20345ac5	173	if (TESTBIT(fMask, kAltroMap)) {
20345ac5	174	AliFMDAltroMapping* altroMap = MakeAltroMap();
02a27b50	175	AliCDBId id(AliFMDParameters::AltroMapPath(), fRunMin, fRunMax);
20345ac5	176	MAKE_META(meta);
	177	meta->SetProperty("key1", altroMap);
	178	cdb->Put(altroMap, id, meta);
	179	cleanup.Add(altroMap);
	180	cleanup.Add(meta);
	181	}
	182	cdb->Destroy();
	183	cleanup.Delete();
	184	}
	185
	186
	187	//__________________________________________________________________
	188	AliFMDCalibGain*
02a27b50	189	AliFMDCalibFaker::MakePulseGain() const
20345ac5	190	{
	191	// Make the actual data
	192	AliFMDCalibGain* gain = new AliFMDCalibGain;
	193	// Set threshold
	194	gain->Set(fThreshold);
	195	for (UShort_t det = 1; det <= 3; det++) {
	196	Char_t rings[] = { 'I', (det == 1 ? '\0' : 'O'), '\0' };
	197	for (Char_t* ring = rings; *ring != '\0'; ring++) {
	198	UShort_t nSec = ( *ring == 'I' ? 20 : 40 );
	199	UShort_t nStr = ( *ring == 'I' ? 512 : 256 );
	200	for (UShort_t sec = 0; sec < nSec; sec++) {
	201	for (UShort_t str = 0; str < nStr; str++) {
	202	gain->Set(det, *ring, sec, str,
	203	gRandom->Gaus(fGain, .01 * fGain));
	204	}
	205	}
	206	}
	207	}
	208	return gain;
	209	}
	210
a9579262	211	//__________________________________________________________________
	212	Float_t
	213	AliFMDCalibFaker::MakeNoise(Char_t ring, UShort_t str) const
	214	{
	215	const UShort_t innerN = 512;
	216	const UShort_t outerN = 256;
	217	const UShort_t innerCut = 350;
	218	const UShort_t outerCut = 190;
	219	const Float_t innerBase = 1.2;
	220	const Float_t outerBase = 2.1;
	221	const Float_t innerInc = 0.5;
	222	const Float_t outerInc = 0.8;
	223	Float_t cut, base, inc, n;
	224	switch (ring) {
	225	case 'I':
	226	cut = innerCut; base = innerBase; inc = innerInc; n = innerN; break;
	227	case 'O':
	228	cut = outerCut; base = outerBase; inc = outerInc; n = outerN; break;
	229	default:
	230	return -1;
	231	}
	232	Float_t bare = base + (str < cut ?
	233	str / cut * inc :
	234	inc - (str - cut) / (n - cut) * inc);
	235	return bare + gRandom->Uniform(-.07, .07);
	236	}
	237
20345ac5	238	//__________________________________________________________________
20345ac5	239	AliFMDCalibPedestal*
02a27b50	240	AliFMDCalibFaker::MakePedestal() const
20345ac5	241	{
	242	// Make the actual data
	243	AliFMDCalibPedestal* pedestal = new AliFMDCalibPedestal;
a9579262	244
20345ac5	245	for (UShort_t det = 1; det <= 3; det++) {
a9579262	246	Char_t rings[] = { 'I', 'O', '\0' };
20345ac5	247	for (Char_t* ring = rings; *ring != '\0'; ring++) {
a9579262	248	if (*ring == 'O' && det == 1) continue;
20345ac5	249	UShort_t nSec = ( *ring == 'I' ? 20 : 40 );
	250	UShort_t nStr = ( *ring == 'I' ? 512 : 256 );
	251	for (UShort_t sec = 0; sec < nSec; sec++) {
	252	for (UShort_t str = 0; str < nStr; str++) {
a9579262	253	Float_t noise = MakeNoise(*ring, str);
	254	Float_t ped = gRandom->Uniform(fPedestalMin, fPedestalMax);
	255	pedestal->Set(det, *ring, sec, str, ped, noise);
	256	}
20345ac5	257	}
	258	}
	259	}
	260	return pedestal;
	261	}
	262
	263	//__________________________________________________________________
	264	AliFMDCalibDeadMap*
02a27b50	265	AliFMDCalibFaker::MakeDeadMap() const
20345ac5	266	{
	267	// Make the actual data
	268	AliFMDCalibDeadMap* deadmap = new AliFMDCalibDeadMap;
	269	for (UShort_t det = 1; det <= 3; det++) {
	270	Char_t rings[] = { 'I', (det == 1 ? '\0' : 'O'), '\0' };
	271	for (Char_t* ring = rings; *ring != '\0'; ring++) {
	272	UShort_t nSec = ( *ring == 'I' ? 20 : 40 );
	273	UShort_t nStr = ( *ring == 'I' ? 512 : 256 );
	274	for (UShort_t sec = 0; sec < nSec; sec++) {
	275	for (UShort_t str = 0; str < nStr; str++) {
	276	deadmap->operator()(det, *ring, sec, str) =
	277	gRandom->Uniform(0, 1) < fDeadChance;
	278	}
	279	}
	280	}
	281	}
	282	return deadmap;
	283	}
	284
	285	//__________________________________________________________________
	286	AliFMDCalibZeroSuppression*
02a27b50	287	AliFMDCalibFaker::MakeZeroSuppression() const
20345ac5	288	{
	289	// Make the actual data
	290	AliFMDCalibZeroSuppression* zs = new AliFMDCalibZeroSuppression;
	291	for (UShort_t det = 1; det <= 3; det++) {
	292	Char_t rings[] = { 'I', (det == 1 ? '\0' : 'O'), '\0' };
	293	for (Char_t* ring = rings; *ring != '\0'; ring++) {
	294	UShort_t nSec = ( *ring == 'I' ? 20 : 40 );
	295	UShort_t nStr = ( *ring == 'I' ? 512 : 256 );
	296	for (UShort_t sec = 0; sec < nSec; sec++) {
	297	for (UShort_t str = 0; str < nStr; str++) {
	298	zs->operator()(det, *ring, sec, str) = fZeroThreshold;
	299	}
	300	}
	301	}
	302	}
	303	return zs;
	304	}
	305
	306	//__________________________________________________________________
	307	AliFMDCalibSampleRate*
02a27b50	308	AliFMDCalibFaker::MakeSampleRate() const
20345ac5	309	{
02a27b50	310	// Make sample rates
20345ac5	311	AliFMDCalibSampleRate* sampleRate = new AliFMDCalibSampleRate;
c2fc1258	312	for (UShort_t det = 1; det <= 3; det++) {
	313	Char_t rings[] = { 'I', (det == 1 ? '\0' : 'O'), '\0' };
	314	for (Char_t* ring = rings; *ring != '\0'; ring++) {
	315	UShort_t nSec = ( *ring == 'I' ? 20 : 40 );
	316	for (UShort_t sec = 0; sec < nSec; sec++) {
	317	sampleRate->Set(det, *ring, sec, 0, fRate);
	318	}
	319	}
	320	}
20345ac5	321	return sampleRate;
	322	}
	323
c2fc1258	324	//__________________________________________________________________
c2fc1258	325	AliFMDCalibStripRange*
02a27b50	326	AliFMDCalibFaker::MakeStripRange() const
c2fc1258	327	{
02a27b50	328	// Make strip ranges
c2fc1258	329	AliFMDCalibStripRange* striprange = new AliFMDCalibStripRange;
	330	for (UShort_t det = 1; det <= 3; det++) {
	331	Char_t rings[] = { 'I', (det == 1 ? '\0' : 'O'), '\0' };
	332	for (Char_t* ring = rings; *ring != '\0'; ring++) {
	333	UShort_t nSec = ( *ring == 'I' ? 20 : 40 );
	334	for (UShort_t sec = 0; sec < nSec; sec++) {
	335	striprange->Set(det, *ring, sec, 0, fStripMin, fStripMax);
	336	}
	337	}
	338	}
	339	return striprange;
	340	}
	341
20345ac5	342	//__________________________________________________________________
20345ac5	343	AliFMDAltroMapping*
02a27b50	344	AliFMDCalibFaker::MakeAltroMap() const
20345ac5	345	{
02a27b50	346	// Make hardware mapping
20345ac5	347	AliFMDAltroMapping* m = new AliFMDAltroMapping;
	348	return m;
	349	}
	350
	351
	352
	353	//____________________________________________________________________
	354	//
	355	// EOF
	356	//