[u/mrichter/AliRoot.git] / MUON / AliMUONTrackerIO.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$

#include <cstdlib>
#include "AliMUONTrackerIO.h"

/// \class AliMUONTrackerIO
///
/// Reader class for ASCII calibration files for MUON tracker : 
/// converts those ASCII files into AliMUONVStore (suitable to e.g. feed
/// the OCDB).
///
/// \author Laurent Aphecetche, Subatech

/// \cond CLASSIMP
ClassImp(AliMUONTrackerIO)
/// \endcond

#include "AliLog.h"
#include "AliMUONCalibParamND.h"
#include "AliMUONCalibParamNF.h"
#include "AliMUONVStore.h"
#include "AliMpConstants.h"
#include "AliMpDDLStore.h"
#include "AliMpDetElement.h"
#include <Riostream.h>
#include <TClass.h>
#include <TObjString.h>
#include <TSystem.h>
#include <sstream>

//_____________________________________________________________________________
AliMUONTrackerIO::AliMUONTrackerIO()
{
  /// ctor
}

//_____________________________________________________________________________
AliMUONTrackerIO::~AliMUONTrackerIO()
{
  /// dtor
}

//_____________________________________________________________________________
Int_t 
AliMUONTrackerIO::ReadOccupancy(const char* filename,AliMUONVStore& occupancyMap)
{
  /// Read occupancy file created by online DA
  /// and append values to the occupancyMap store.
  /// Expected format of the file is :
  /// busPatchId manuId sumofn nevt
  
  TString sFilename(gSystem->ExpandPathName(filename));
  
  std::ifstream in(sFilename.Data());
  if (!in.good()) 
  {
    return kCannotOpenFile;
  }
  
  TString datastring;
  ostringstream stream;
  char line[1024];
  while ( in.getline(line,1024) )
  	stream << line << "\n";
  datastring = TString(stream.str().c_str());
  
  in.close();
  
  return DecodeOccupancy(datastring,occupancyMap);
  
}

//_____________________________________________________________________________
Int_t 
AliMUONTrackerIO::DecodeOccupancy(TString data, AliMUONVStore& occupancyMap)
{
  /// Decode occupancy string created append values to the occupancyMap store.
  /// Expected format of the file is :
  /// busPatchId manuId sumofn nevt
 
  if ( ! AliMpDDLStore::Instance(kFALSE) )
  {
    AliErrorClass("Mapping not loaded. Cannot work");
    return 0;
  }
  
  char line[1024];
  istringstream in(data.Data());
  
  Int_t n(0);
  
  while ( in.getline(line,1024) )
  {
    AliDebugClass(3,Form("line=%s",line));
    if ( line[0] == '/' && line[1] == '/' ) continue;
    std::istringstream sin(line);
    
    Int_t busPatchId, manuId;
    Int_t numberOfEvents;
    Double_t sumn;

    sin >> busPatchId >> manuId >> sumn >> numberOfEvents;
    
    Int_t detElemId = AliMpDDLStore::Instance()->GetDEfromBus(busPatchId);

    AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);

    Int_t numberOfChannelsInManu = -1;
    
    if (de) numberOfChannelsInManu = de->NofChannelsInManu(manuId);

    if ( numberOfChannelsInManu <= 0 ) 
    {
      AliErrorClass(Form("BP %5d DE %5d MANU %5d nchannels=%d",busPatchId,detElemId,manuId,numberOfChannelsInManu));
      continue;      
    }
    
    AliMUONVCalibParam* occupancy = 
    static_cast<AliMUONVCalibParam*>(occupancyMap.FindObject(detElemId,manuId));
    if (occupancy) 
    {
      AliErrorClass(Form("DE %5d MANU %5d is already there ?!",detElemId,manuId));
      continue;
    }
        
    occupancy = new AliMUONCalibParamND(5,1,detElemId,manuId,0);

    occupancyMap.Add(occupancy);
    
    occupancy->SetValueAsDouble(0,0,sumn);
    occupancy->SetValueAsDouble(0,1,sumn); // with only 0 and 1s, sumw = sumw2 = sumn
    occupancy->SetValueAsDouble(0,2,sumn);
    occupancy->SetValueAsInt(0,3,numberOfChannelsInManu);
    occupancy->SetValueAsInt(0,4,numberOfEvents);
    ++n;
  }
  
  return n;
}

//_____________________________________________________________________________
Int_t 
AliMUONTrackerIO::ReadPedestals(const char* filename, AliMUONVStore& pedStore)
{
  /// Read pedestal file (produced by the MUONTRKda.exe program for instance)
  /// and append the read values into the given VStore
  /// To be used when the input is a file (for instance when reading data 
  /// from the OCDB).
  
  TString sFilename(gSystem->ExpandPathName(filename));
  
  std::ifstream in(sFilename.Data());
  if (!in.good()) 
  {
    return kCannotOpenFile;
  }
  
  TString datastring;
  ostringstream stream;
  char line[1024];
  while ( in.getline(line,1024) )
  	stream << line << "\n";
  datastring = TString(stream.str().c_str());
  
  in.close();

  return DecodePedestals(datastring,pedStore);
  
}

//_____________________________________________________________________________
Int_t 
AliMUONTrackerIO::DecodePedestals(TString data, AliMUONVStore& pedStore)
{
  /// Read pedestal Data (produced by the MUONTRKda.exe program for instance)
  /// and append the read values into the given VStore
  /// To be used when the input is a TString (for instance when getting data 
  /// from AMORE DB).
  
  char line[1024];
  Int_t busPatchID, manuID, manuChannel;
  Float_t pedMean, pedSigma;
  Int_t n(0);
  istringstream in(data.Data());
  
  while ( in.getline(line,1024) )
  {
    AliDebugClass(3,Form("line=%s",line));
    if ( line[0] == '/' && line[1] == '/' ) continue;
    std::istringstream sin(line);
    sin >> busPatchID >> manuID >> manuChannel >> pedMean >> pedSigma;
    Int_t detElemID = AliMpDDLStore::Instance()->GetDEfromBus(busPatchID);
    AliDebugClass(3,Form("BUSPATCH %3d DETELEMID %4d MANU %3d CH %3d MEAN %7.2f SIGMA %7.2f",
                    busPatchID,detElemID,manuID,manuChannel,pedMean,pedSigma));
		    
    AliMUONVCalibParam* ped = 
      static_cast<AliMUONVCalibParam*>(pedStore.FindObject(detElemID,manuID));
    if (!ped) 
    {
      ped = new AliMUONCalibParamNF(2,AliMpConstants::ManuNofChannels(),
                                    detElemID,manuID,
                                    AliMUONVCalibParam::InvalidFloatValue());  
      pedStore.Add(ped);
    }
    ped->SetValueAsFloat(manuChannel,0,pedMean);
    ped->SetValueAsFloat(manuChannel,1,pedSigma);
    ++n;
  }

  return n;
}

//_____________________________________________________________________________
Int_t 
AliMUONTrackerIO::ReadGains(const char* filename, AliMUONVStore& gainStore,
                            TString& comment)
{
  /// Read gain file (produced by the MUONTRKda.exe program for instance)
  /// and append the read values into the given VStore
  /// To be used when the input is a file (for instance when reading data 
  /// from the OCDB).
  
  comment = "";
  
  TString sFilename(gSystem->ExpandPathName(filename));
  
  std::ifstream in(sFilename.Data());
  if (!in.good()) 
  {
    return kCannotOpenFile;
  }
  
  TString datastring;
  ostringstream stream;
  char line[1024];
  while ( in.getline(line,1024) )
  	stream << line << "\n";
  datastring = TString(stream.str().c_str());
  
  in.close();
  
  return DecodeGains(datastring,gainStore,comment);

}

//_____________________________________________________________________________
Int_t 
AliMUONTrackerIO::DecodeGains(TString data, AliMUONVStore& gainStore,
                            TString& comment)
{
  /// Read gain file (produced by the MUONTRKda.exe program for instance)
  /// and append the read values into the given VStore
  /// To be used when the input is a TString (for instance when getting data 
  /// from AMORE DB).
  
  char line[1024];
  istringstream in(data.Data());
  Int_t busPatchID, manuID, manuChannel;
  Float_t a0, a1;
  Int_t thres;
  UInt_t qual;
  const Int_t kSaturation(3000); // FIXME: how to get this number ?
  Int_t n(0);
  Int_t runNumber(-1);
  Int_t* runs(0x0);
  Int_t* dac(0x0);
  Int_t nDAC(0);
  Int_t iDAC(0);
  
  while ( in.getline(line,1024) )
  {
    if ( strlen(line) < 10 ) continue;
    if ( line[0] == '/' && line[1] == '/' ) 
    {
      TString sline(line);
      if ( sline.Contains("DUMMY") )
      {
        AliDebugClass(1,"Got a dummy file here");
        return kDummyFile;
      }
      if ( sline.Contains("* Run") )
      {
        TObjArray* a = sline.Tokenize(":");
        if ( a->GetLast() >= 1 ) 
        {
          TString s = static_cast<TObjString*>(a->At(1))->String();
          runNumber = s.Atoi();
          AliDebugClass(1,Form("runNumber is %d",runNumber));
        }            
      }
      if ( sline.Contains("DAC values") )
      {
        nDAC = TString(sline(2,sline.Length()-2)).Atoi();
        AliDebugClass(1,Form("# of DAC values = %d",nDAC));
        if ( nDAC > 0 )
        {
          if ( nDAC < 100 ) 
          {
            runs = new Int_t[nDAC];
            dac = new Int_t[nDAC];
            // skip two lines
            in.getline(line,1024);
            in.getline(line,1024);
            // then get run and dac values
            for ( Int_t i = 0; i < nDAC; ++i ) 
            {
              in.getline(line,1024);
              Int_t a,b;
              sscanf(line,"// %d %d",&a,&b);
              runs[iDAC] = a;
              dac[iDAC] = b;
              AliDebugClass(1,Form("RUN %d is DAC %d",runs[iDAC],dac[iDAC]));
              ++iDAC;
            }
          }
          else
          {
            AliErrorClass(Form("Something went wrong, as I get too big nDAC = %d",nDAC));
            nDAC = 0;
            return kFormatError;
          }
        }
      }
      continue;
    }
    
    sscanf(line,"%d %d %d %f %f %d %x",&busPatchID,&manuID,&manuChannel,
           &a0,&a1,&thres,&qual); 
    AliDebugClass(3,Form("line=%s",line));
    Int_t detElemID = AliMpDDLStore::Instance()->GetDEfromBus(busPatchID);
    AliDebugClass(3,Form("BUSPATCH %3d DETELEMID %4d MANU %3d CH %3d A0 %7.2f "
                    "A1 %e THRES %5d QUAL %x",
                    busPatchID,detElemID,manuID,manuChannel,a0,a1,thres,qual));
    if ( qual == 0 ) continue;
    
    AliMUONVCalibParam* gain = 
      static_cast<AliMUONVCalibParam*>(gainStore.FindObject(detElemID,manuID));
    
   if (!gain) 
    {
      gain = new AliMUONCalibParamNF(5,AliMpConstants::ManuNofChannels(),detElemID,manuID,0);
      gainStore.Add(gain);
    }
    gain->SetValueAsFloat(manuChannel,0,a0);
    gain->SetValueAsFloat(manuChannel,1,a1);
    gain->SetValueAsInt(manuChannel,2,thres);
    gain->SetValueAsInt(manuChannel,3,qual);
    gain->SetValueAsInt(manuChannel,4,kSaturation);
    ++n;
  }

  comment = "";
  
  if ( runNumber > 0 )
  {
    comment = Form("RUN %d",runNumber);
  }
  
  for ( Int_t i = 0; i < nDAC; ++i )
  {
    comment += Form(";(RUN %d = DAC %d)",runs[i],dac[i]);
  }
  
  delete[] runs;
  delete[] dac;
  
  return n;
}

//_____________________________________________________________________________
Int_t
AliMUONTrackerIO::ReadCapacitances(const char* file, AliMUONVStore& capaStore)
{
  /// Read capacitance file
  /// and append the read values into the given VStore
  
  ifstream in(gSystem->ExpandPathName(file));
  if (in.bad()) return kCannotOpenFile;
  
  Int_t ngenerated(0);
  
  char line[1024];
  Int_t serialNumber(-1);
  AliMUONVCalibParam* param(0x0);
  
  while ( in.getline(line,1024,'\n') )
  {
    if ( isdigit(line[0]) ) 
    {
      serialNumber = atoi(line);
      param = static_cast<AliMUONVCalibParam*>(capaStore.FindObject(serialNumber));
      if (param)
      {
        AliErrorClass(Form("serialNumber %d appears several times !",serialNumber));
        continue;
      }
      param = new AliMUONCalibParamNF(2,AliMpConstants::ManuNofChannels(),serialNumber,0,1.0);
      Bool_t ok = capaStore.Add(param);
      if (!ok)
      {
        AliErrorClass(Form("Could not set serialNumber=%d",serialNumber));
        continue;
      }      
      continue;
    }
    Int_t channel;
    Float_t capaValue;
    Float_t injectionGain;
    sscanf(line,"%d %f %f",&channel,&capaValue,&injectionGain);
    AliDebugClass(1,Form("SerialNumber %10d Channel %3d Capa %f injectionGain %f",
                    serialNumber,channel,capaValue,injectionGain));
    param->SetValueAsFloat(channel,0,capaValue);
    param->SetValueAsFloat(channel,1,injectionGain);
    ++ngenerated;
  }
  
  in.close();
  
  return ngenerated;
}
Commit	Line	Data
81028269	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	**************************************************************************/
	15
	16	// $Id$
	17
b09247a2	18	#include <cstdlib>
81028269	19	#include "AliMUONTrackerIO.h"
	20
	21	/// \class AliMUONTrackerIO
	22	///
	23	/// Reader class for ASCII calibration files for MUON tracker :
	24	/// converts those ASCII files into AliMUONVStore (suitable to e.g. feed
	25	/// the OCDB).
	26	///
	27	/// \author Laurent Aphecetche, Subatech
	28
	29	/// \cond CLASSIMP
	30	ClassImp(AliMUONTrackerIO)
	31	/// \endcond
	32
	33	#include "AliLog.h"
7eafe398	34	#include "AliMUONCalibParamND.h"
81028269	35	#include "AliMUONCalibParamNF.h"
	36	#include "AliMUONVStore.h"
	37	#include "AliMpConstants.h"
	38	#include "AliMpDDLStore.h"
7eafe398	39	#include "AliMpDetElement.h"
81028269	40	#include <Riostream.h>
	41	#include <TClass.h>
	42	#include <TObjString.h>
	43	#include <TSystem.h>
	44	#include <sstream>
	45
	46	//_____________________________________________________________________________
	47	AliMUONTrackerIO::AliMUONTrackerIO()
	48	{
	49	/// ctor
	50	}
	51
	52	//_____________________________________________________________________________
	53	AliMUONTrackerIO::~AliMUONTrackerIO()
	54	{
	55	/// dtor
	56	}
	57
7eafe398	58	//_____________________________________________________________________________
	59	Int_t
	60	AliMUONTrackerIO::ReadOccupancy(const char* filename,AliMUONVStore& occupancyMap)
	61	{
	62	/// Read occupancy file created by online DA
	63	/// and append values to the occupancyMap store.
	64	/// Expected format of the file is :
	65	/// busPatchId manuId sumofn nevt
	66
	67	TString sFilename(gSystem->ExpandPathName(filename));
	68
	69	std::ifstream in(sFilename.Data());
	70	if (!in.good())
	71	{
	72	return kCannotOpenFile;
	73	}
	74
	75	TString datastring;
	76	ostringstream stream;
	77	char line[1024];
	78	while ( in.getline(line,1024) )
	79	stream << line << "\n";
	80	datastring = TString(stream.str().c_str());
	81
	82	in.close();
	83
	84	return DecodeOccupancy(datastring,occupancyMap);
	85
	86	}
	87
	88	//_____________________________________________________________________________
	89	Int_t
	90	AliMUONTrackerIO::DecodeOccupancy(TString data, AliMUONVStore& occupancyMap)
	91	{
	92	/// Decode occupancy string created append values to the occupancyMap store.
	93	/// Expected format of the file is :
	94	/// busPatchId manuId sumofn nevt
	95
	96	if ( ! AliMpDDLStore::Instance(kFALSE) )
	97	{
	98	AliErrorClass("Mapping not loaded. Cannot work");
	99	return 0;
	100	}
	101
	102	char line[1024];
	103	istringstream in(data.Data());
	104
	105	Int_t n(0);
	106
	107	while ( in.getline(line,1024) )
	108	{
	109	AliDebugClass(3,Form("line=%s",line));
	110	if ( line[0] == '/' && line[1] == '/' ) continue;
	111	std::istringstream sin(line);
	112
	113	Int_t busPatchId, manuId;
	114	Int_t numberOfEvents;
	115	Double_t sumn;
	116
	117	sin >> busPatchId >> manuId >> sumn >> numberOfEvents;
	118
	119	Int_t detElemId = AliMpDDLStore::Instance()->GetDEfromBus(busPatchId);
	120
	121	AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
122
123	Int_t numberOfChannelsInManu = -1;
124
125	if (de) numberOfChannelsInManu = de->NofChannelsInManu(manuId);
126
127	if ( numberOfChannelsInManu <= 0 )
128	{
129	AliErrorClass(Form("BP %5d DE %5d MANU %5d nchannels=%d",busPatchId,detElemId,manuId,numberOfChannelsInManu));
130	continue;
131	}
132
133	AliMUONVCalibParam* occupancy =
134	static_cast<AliMUONVCalibParam*>(occupancyMap.FindObject(detElemId,manuId));
135	if (occupancy)
136	{
137	AliErrorClass(Form("DE %5d MANU %5d is already there ?!",detElemId,manuId));
138	continue;
139	}
140
141	occupancy = new AliMUONCalibParamND(5,1,detElemId,manuId,0);
142
143	occupancyMap.Add(occupancy);
144
145	occupancy->SetValueAsDouble(0,0,sumn);
146	occupancy->SetValueAsDouble(0,1,sumn); // with only 0 and 1s, sumw = sumw2 = sumn
147	occupancy->SetValueAsDouble(0,2,sumn);
148	occupancy->SetValueAsInt(0,3,numberOfChannelsInManu);
149	occupancy->SetValueAsInt(0,4,numberOfEvents);
150	++n;
151	}
152
153	return n;
154	}
155
81028269	156	//_____________________________________________________________________________
	157	Int_t
	158	AliMUONTrackerIO::ReadPedestals(const char* filename, AliMUONVStore& pedStore)
	159	{
	160	/// Read pedestal file (produced by the MUONTRKda.exe program for instance)
	161	/// and append the read values into the given VStore
cba13f7c	162	/// To be used when the input is a file (for instance when reading data
cba13f7c	163	/// from the OCDB).
81028269	164
	165	TString sFilename(gSystem->ExpandPathName(filename));
	166
	167	std::ifstream in(sFilename.Data());
	168	if (!in.good())
	169	{
	170	return kCannotOpenFile;
	171	}
	172
cba13f7c	173	TString datastring;
06a51226	174	ostringstream stream;
	175	char line[1024];
	176	while ( in.getline(line,1024) )
	177	stream << line << "\n";
	178	datastring = TString(stream.str().c_str());
	179
cba13f7c	180	in.close();
	181
	182	return DecodePedestals(datastring,pedStore);
	183
	184	}
	185
	186	//_____________________________________________________________________________
	187	Int_t
	188	AliMUONTrackerIO::DecodePedestals(TString data, AliMUONVStore& pedStore)
	189	{
	190	/// Read pedestal Data (produced by the MUONTRKda.exe program for instance)
	191	/// and append the read values into the given VStore
	192	/// To be used when the input is a TString (for instance when getting data
	193	/// from AMORE DB).
	194
81028269	195	char line[1024];
	196	Int_t busPatchID, manuID, manuChannel;
	197	Float_t pedMean, pedSigma;
	198	Int_t n(0);
cba13f7c	199	istringstream in(data.Data());
81028269	200
	201	while ( in.getline(line,1024) )
	202	{
	203	AliDebugClass(3,Form("line=%s",line));
	204	if ( line[0] == '/' && line[1] == '/' ) continue;
	205	std::istringstream sin(line);
	206	sin >> busPatchID >> manuID >> manuChannel >> pedMean >> pedSigma;
	207	Int_t detElemID = AliMpDDLStore::Instance()->GetDEfromBus(busPatchID);
	208	AliDebugClass(3,Form("BUSPATCH %3d DETELEMID %4d MANU %3d CH %3d MEAN %7.2f SIGMA %7.2f",
	209	busPatchID,detElemID,manuID,manuChannel,pedMean,pedSigma));
cba13f7c	210
81028269	211	AliMUONVCalibParam* ped =
81028269	212	static_cast<AliMUONVCalibParam*>(pedStore.FindObject(detElemID,manuID));
81028269	213	if (!ped)
	214	{
	215	ped = new AliMUONCalibParamNF(2,AliMpConstants::ManuNofChannels(),
	216	detElemID,manuID,
	217	AliMUONVCalibParam::InvalidFloatValue());
	218	pedStore.Add(ped);
	219	}
	220	ped->SetValueAsFloat(manuChannel,0,pedMean);
	221	ped->SetValueAsFloat(manuChannel,1,pedSigma);
	222	++n;
	223	}
81028269	224
	225	return n;
	226	}
	227
	228	//_____________________________________________________________________________
	229	Int_t
	230	AliMUONTrackerIO::ReadGains(const char* filename, AliMUONVStore& gainStore,
	231	TString& comment)
	232	{
	233	/// Read gain file (produced by the MUONTRKda.exe program for instance)
	234	/// and append the read values into the given VStore
cba13f7c	235	/// To be used when the input is a file (for instance when reading data
cba13f7c	236	/// from the OCDB).
81028269	237
	238	comment = "";
	239
	240	TString sFilename(gSystem->ExpandPathName(filename));
	241
	242	std::ifstream in(sFilename.Data());
	243	if (!in.good())
	244	{
	245	return kCannotOpenFile;
	246	}
	247
cba13f7c	248	TString datastring;
	249	ostringstream stream;
	250	char line[1024];
	251	while ( in.getline(line,1024) )
	252	stream << line << "\n";
	253	datastring = TString(stream.str().c_str());
	254
	255	in.close();
	256
	257	return DecodeGains(datastring,gainStore,comment);
	258
	259	}
	260
	261	//_____________________________________________________________________________
	262	Int_t
	263	AliMUONTrackerIO::DecodeGains(TString data, AliMUONVStore& gainStore,
	264	TString& comment)
	265	{
	266	/// Read gain file (produced by the MUONTRKda.exe program for instance)
	267	/// and append the read values into the given VStore
	268	/// To be used when the input is a TString (for instance when getting data
	269	/// from AMORE DB).
	270
81028269	271	char line[1024];
cba13f7c	272	istringstream in(data.Data());
81028269	273	Int_t busPatchID, manuID, manuChannel;
	274	Float_t a0, a1;
	275	Int_t thres;
	276	UInt_t qual;
	277	const Int_t kSaturation(3000); // FIXME: how to get this number ?
	278	Int_t n(0);
	279	Int_t runNumber(-1);
	280	Int_t* runs(0x0);
	281	Int_t* dac(0x0);
	282	Int_t nDAC(0);
	283	Int_t iDAC(0);
	284
	285	while ( in.getline(line,1024) )
	286	{
	287	if ( strlen(line) < 10 ) continue;
	288	if ( line[0] == '/' && line[1] == '/' )
	289	{
	290	TString sline(line);
	291	if ( sline.Contains("DUMMY") )
	292	{
	293	AliDebugClass(1,"Got a dummy file here");
	294	return kDummyFile;
	295	}
	296	if ( sline.Contains("* Run") )
	297	{
	298	TObjArray* a = sline.Tokenize(":");
	299	if ( a->GetLast() >= 1 )
	300	{
	301	TString s = static_cast<TObjString*>(a->At(1))->String();
	302	runNumber = s.Atoi();
	303	AliDebugClass(1,Form("runNumber is %d",runNumber));
	304	}
	305	}
	306	if ( sline.Contains("DAC values") )
	307	{
	308	nDAC = TString(sline(2,sline.Length()-2)).Atoi();
	309	AliDebugClass(1,Form("# of DAC values = %d",nDAC));
	310	if ( nDAC > 0 )
	311	{
	312	if ( nDAC < 100 )
	313	{
	314	runs = new Int_t[nDAC];
	315	dac = new Int_t[nDAC];
	316	// skip two lines
	317	in.getline(line,1024);
	318	in.getline(line,1024);
	319	// then get run and dac values
	320	for ( Int_t i = 0; i < nDAC; ++i )
	321	{
	322	in.getline(line,1024);
	323	Int_t a,b;
	324	sscanf(line,"// %d %d",&a,&b);
	325	runs[iDAC] = a;
	326	dac[iDAC] = b;
	327	AliDebugClass(1,Form("RUN %d is DAC %d",runs[iDAC],dac[iDAC]));
	328	++iDAC;
	329	}
	330	}
	331	else
	332	{
	333	AliErrorClass(Form("Something went wrong, as I get too big nDAC = %d",nDAC));
	334	nDAC = 0;
	335	return kFormatError;
	336	}
337	}
338	}
339	continue;
340	}
341
342	sscanf(line,"%d %d %d %f %f %d %x",&busPatchID,&manuID,&manuChannel,
343	&a0,&a1,&thres,&qual);
344	AliDebugClass(3,Form("line=%s",line));
345	Int_t detElemID = AliMpDDLStore::Instance()->GetDEfromBus(busPatchID);
346	AliDebugClass(3,Form("BUSPATCH %3d DETELEMID %4d MANU %3d CH %3d A0 %7.2f "
347	"A1 %e THRES %5d QUAL %x",
348	busPatchID,detElemID,manuID,manuChannel,a0,a1,thres,qual));
349	if ( qual == 0 ) continue;
350
351	AliMUONVCalibParam* gain =
352	static_cast<AliMUONVCalibParam*>(gainStore.FindObject(detElemID,manuID));
353
cba13f7c	354	if (!gain)
81028269	355	{
	356	gain = new AliMUONCalibParamNF(5,AliMpConstants::ManuNofChannels(),detElemID,manuID,0);
	357	gainStore.Add(gain);
	358	}
	359	gain->SetValueAsFloat(manuChannel,0,a0);
	360	gain->SetValueAsFloat(manuChannel,1,a1);
	361	gain->SetValueAsInt(manuChannel,2,thres);
	362	gain->SetValueAsInt(manuChannel,3,qual);
	363	gain->SetValueAsInt(manuChannel,4,kSaturation);
	364	++n;
	365	}
cba13f7c	366
81028269	367	comment = "";
	368
	369	if ( runNumber > 0 )
	370	{
	371	comment = Form("RUN %d",runNumber);
	372	}
	373
	374	for ( Int_t i = 0; i < nDAC; ++i )
	375	{
	376	comment += Form(";(RUN %d = DAC %d)",runs[i],dac[i]);
	377	}
	378
	379	delete[] runs;
	380	delete[] dac;
	381
	382	return n;
	383	}
	384
	385	//_____________________________________________________________________________
	386	Int_t
	387	AliMUONTrackerIO::ReadCapacitances(const char* file, AliMUONVStore& capaStore)
	388	{
	389	/// Read capacitance file
	390	/// and append the read values into the given VStore
	391
	392	ifstream in(gSystem->ExpandPathName(file));
	393	if (in.bad()) return kCannotOpenFile;
	394
	395	Int_t ngenerated(0);
	396
	397	char line[1024];
	398	Int_t serialNumber(-1);
	399	AliMUONVCalibParam* param(0x0);
	400
	401	while ( in.getline(line,1024,'\n') )
	402	{
	403	if ( isdigit(line[0]) )
	404	{
	405	serialNumber = atoi(line);
	406	param = static_cast<AliMUONVCalibParam*>(capaStore.FindObject(serialNumber));
	407	if (param)
	408	{
	409	AliErrorClass(Form("serialNumber %d appears several times !",serialNumber));
630711ed	410	continue;
81028269	411	}
	412	param = new AliMUONCalibParamNF(2,AliMpConstants::ManuNofChannels(),serialNumber,0,1.0);
	413	Bool_t ok = capaStore.Add(param);
	414	if (!ok)
	415	{
	416	AliErrorClass(Form("Could not set serialNumber=%d",serialNumber));
	417	continue;
	418	}
	419	continue;
	420	}
	421	Int_t channel;
	422	Float_t capaValue;
	423	Float_t injectionGain;
	424	sscanf(line,"%d %f %f",&channel,&capaValue,&injectionGain);
	425	AliDebugClass(1,Form("SerialNumber %10d Channel %3d Capa %f injectionGain %f",
	426	serialNumber,channel,capaValue,injectionGain));
	427	param->SetValueAsFloat(channel,0,capaValue);
	428	param->SetValueAsFloat(channel,1,injectionGain);
	429	++ngenerated;
	430	}
	431
	432	in.close();
	433
	434	return ngenerated;
	435	}