[u/mrichter/AliRoot.git] / MUON / MUONGenerateBusPatch.C

/**************************************************************************
 * 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$

#if !defined(__CINT__) || defined(__MAKECINT__)
// ROOT includes
#include <Riostream.h>
#include <TROOT.h>
#include <TSystem.h>
#include <TString.h>

// MUON includes
#include "AliMpBusPatch.h"

#endif

// Generates buspatch id and DDL id for given detection element id
// station 1 & 2 assuming 24 buspatches per quadrant
// station345, reading from DetElemIdToSlatType.dat file and calculates
// the number of bus patches per slat (and also number of Translator and Bridge Boards).
// Generates an output file DetElemIdToBusPatch.dat.out, preserve from overwriting
// (Ch. Finck, July 05)
// (Nov. 05, added DDL)
// (June 06, correction for St123)
// (Feb. 07, add 1st manu list for St12 (starting on NB !) and new ddl sharing for station 3)
// (June. 07, new numbering of station 345, and buspatch starts at 1)


void MUONGenerateBusPatch()
{

  TString dirPath2 = gSystem->Getenv("ALICE_ROOT");
  dirPath2 += "/MUON/mapping/data/"; 

  TString dirPath1 = dirPath2 + "station345/";

  TString infile =  dirPath1 + "DetElemIdToSlatType.dat";
  TString outfile = dirPath2 + "DetElemIdToBusPatch.dat.out";

  ifstream in(infile, ios::in);
  ofstream out(outfile, ios::out);
  char line[80];
  char nameSlat[82];
  Int_t idDE;
  Int_t i;
  Int_t cursor = 0;
  Int_t begin[900];
  Int_t end[900];
  Int_t nbBusPatch;
  Int_t nbTB = 0;
  Int_t nbBB = 0;
  Int_t idSt12[] = {100, 101, 102, 103, 
		    200, 201, 202, 203, 
		    300, 301, 302, 303,
		    400, 401, 402, 403};

  Char_t manuListSt1[] 
      = " 1,27,53,79,105,131,157,183,201,214,224,232,1025,1051,1077,1103,1129,1155,1181,1207,1225,1238,1249,1257";

  Char_t manuListSt2[] 
      = " 1,27,53,79,105,131,157,183,201,214,226,246,1025,1051,1077,1103,1129,1155,1181,1207,1225,1238,1251,1269";


  Int_t iDDL = 0;
  // station 1 & 2
  nbBusPatch = 24;
  Int_t nbHalfBusPatch =  nbBusPatch/2;
  cout << "#DE BusPatch DDL SlatName" << endl;
  out << "#DE BusPatch DDL  1st manu in buspatch" << endl;

  for (Int_t j = 0; j < 16; j++) {

    idDE = idSt12[j];
    if (idDE % 100 == 0) {
      iDDL++;
      begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1);
      cout << "# Chamber " << idDE/100 << endl;
      out  << "# Chamber " << idDE/100 << endl;
    }
    if (idDE % 100 == 1) { 
      iDDL++;
      begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1);
    }
   
    if (idDE % 100 == 3) {
      iDDL--;
      begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1) + nbBusPatch;
    }

    end[cursor]     = begin[cursor] + nbBusPatch - 1;
    ++cursor; // it seems that vec[++i] does not work in Cint
    begin[cursor] = end[cursor-1] + 1;

    cout << idDE << " " << begin[cursor-1] + nbHalfBusPatch << "-" <<end[cursor-1] << ";" <<
	begin[cursor-1]	<< "-" << end[cursor-1] - nbHalfBusPatch << " " << iDDL-1 << endl;

    if (idDE < 300 )
	out << idDE  << " " << begin[cursor-1] + nbHalfBusPatch << "-" <<end[cursor-1] << ";" <<
	    begin[cursor-1]	<< "-" << end[cursor-1] - nbHalfBusPatch << " " << iDDL-1 <<
	    manuListSt1 << endl;
    else
	out << idDE  << " " << begin[cursor-1] + nbHalfBusPatch << "-" <<end[cursor-1] << ";" <<
	    begin[cursor-1]	<< "-" << end[cursor-1] - nbHalfBusPatch << " " <<  iDDL-1 <<
	    manuListSt2  << endl;
    if (idDE % 100 == 3) iDDL++;

  }
  
  // station 345
  nbBusPatch = 0;
  Int_t idCh5swp1 = 501; // 1/4 chamber for DDL on horizontal
  Int_t idCh5swp2 = 505; // 1/4 chamber for DDL on horizontal
  Int_t idCh5swp3 = 510; // 1/4 chamber for DDL on horizontal
  Int_t idCh5swp4 = 514; // 1/4 chamber for DDL on horizontal

  Int_t idCh6swp1 = 600; // 1/4 chamber for DDL on horizontal
  Int_t idCh6swp2 = 605; // 1/4 chamber for DDL on horizontal
  Int_t idCh6swp3 = 609; // 1/4 chamber for DDL on horizontal
  Int_t idCh6swp4 = 614; // 1/4 chamber for DDL on horizontal

  Int_t idSt45swp1 = 7;  // half chamber for DDL in vertical cutting twice the official numbering
  Int_t idSt45swp2 = 20;

  Int_t offsetBusCh5swp1 = 13; // number of buspatches between DE 501-504
  Int_t offsetBusCh5swp2 = 17; // number of buspatches before DE 510

  Int_t offsetBusCh6swp1 = 17; // number of buspatches before DE 600-604
  Int_t offsetBusCh6swp2 = 13; // number of buspatches between DE 614-617

  Int_t offsetBusSt4 = 25;     // number of buspatches before DE 700-800
  Int_t offsetBusSt4swp = 46;  // number of buspatches between 801-806 (701-706)

  Int_t offsetBusSt5 = 27;     // number of buspatches before DE 900-1000
  Int_t offsetBusSt5swp = 50;  // number of buspatches between 901-906 (1001-1006)

  // reads from file
  while ( in.getline(line,80) ) {

    if ( line[0] == '#' || line[0] == '\0' ) continue;
    i = 0;
    sscanf(line, "%d %s", &idDE, nameSlat);

    // number of PCB's
    while(nameSlat[i] < '4' && nameSlat[i] !='0') i++;
    Int_t len = strlen(nameSlat);

    // number of buspatch 
    if (i == 2 || i == 3)
      nbBusPatch = 2;
    else
      nbBusPatch = 4;

    // calculate the number of TB & BB
    // R1: 2 bridges more per PCB
    // R2 & R3: 1 translator more per PCB
    if (i == 2) nbBB += 2; 
    if (i == 3) nbBB += 4; 
    if (i == 4) nbBB += 4; 
    if (i == 5) nbBB += 6;
    if (i == 6) nbBB += 8;

    nbTB += nbBusPatch;
    if (nameSlat[len-1] == '1') {
      nbTB += 2;
      nbBusPatch+=2;
    }
    if (nameSlat[len-1] == '2' || nameSlat[len-1] == '3') {
      nbTB += 1;
      nbBusPatch++;
    }

    // station 3
    // for buspatch length reasons, one ddl connects one 1/4 of two chambers
    // really messy isn't it ? 
    // much more with the new DDL sharing for station 3

    if (idDE < 700 ) {
  
      if (idDE % 100 == 0) {
	cout << "# Chamber " << idDE/100 << endl;
	out  << "# Chamber " << idDE/100 << endl;
      }

      // chamber 5
      if (idDE  == 500 ) {
	iDDL = (idDE/100)*2+2;
	begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1);
	
	end[cursor]   = begin[cursor] + nbBusPatch - 1;
	++cursor;
	begin[cursor] = end[cursor-1] + 1;

	cout << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	out << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] <<" " << iDDL-1 <<endl;
      }

      if (idDE  >= idCh5swp1 && idDE < idCh5swp2) {
	
	if ( idDE ==  idCh5swp1) {
	  iDDL = (idDE/100)*2 ;
	  begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusCh5swp1, iDDL-1);
	}
	
	end[cursor]   = begin[cursor] - nbBusPatch + 1;
	++cursor;
	begin[cursor] = end[cursor-1] - 1;

	cout << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	out << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] <<" " << iDDL-1 <<endl;
      }

      if (idDE >=  idCh5swp2  && idDE  < idCh5swp3) {

	if (idDE ==  idCh5swp2) {
	  iDDL = (idDE/100)*2-1;
	      begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1);	  
	}

	end[cursor]   = begin[cursor] + nbBusPatch - 1;
	++cursor;
	begin[cursor] = end[cursor-1] + 1;

	cout << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	out << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] <<" " << iDDL-1 <<endl;
      }

      if (idDE >=  idCh5swp3  && idDE  < idCh5swp4) {

	if (idDE ==  idCh5swp3) {
	  iDDL = (idDE/100)*2+1;
	      begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1);	  
	}

	end[cursor]   = begin[cursor] + nbBusPatch - 1;
	++cursor;
	begin[cursor] = end[cursor-1] + 1;

	cout << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	out << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] <<" " << iDDL-1 <<endl;
      }

     if (idDE  >= idCh5swp4 && idDE < 600) {
	
	if ( idDE ==  idCh5swp4) {
	  iDDL = (idDE/100)*2+2 ;
	  begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusCh5swp2, iDDL-1);
	}
	
	end[cursor]   = begin[cursor] - nbBusPatch + 1;
	++cursor;
	begin[cursor] = end[cursor-1] - 1;

	cout << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	out << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] <<" " << iDDL-1 <<endl;
      }

     // chamber 6
      if (idDE  >= idCh6swp1 && idDE < idCh6swp2) {
	
	if ( idDE ==  idCh6swp1) {
	  iDDL = (idDE/100)*2-2 ;
	  begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusCh5swp1+offsetBusCh6swp1, iDDL-1);
	}
	
	end[cursor]   = begin[cursor] - nbBusPatch + 1;
	++cursor;
	begin[cursor] = end[cursor-1] - 1;

	cout << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	out << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] <<" " << iDDL-1 <<endl;
      }

      if (idDE >=  idCh6swp2  && idDE  < idCh6swp3) {

	if (idDE ==  idCh6swp2) {
	  iDDL = (idDE/100)*2-3;
	      begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusCh5swp2+1, iDDL-1);	  
	}

	end[cursor]   = begin[cursor] + nbBusPatch - 1;
	++cursor;
	begin[cursor] = end[cursor-1] + 1;

	cout << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	out << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] <<" " << iDDL-1 <<endl;
      }

      if (idDE >=  idCh6swp3  && idDE  < idCh6swp4) {

	if (idDE ==  idCh6swp3) {
	  iDDL = (idDE/100)*2-1;
	      begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusCh5swp1+1, iDDL-1);	  
	}

	end[cursor]   = begin[cursor] + nbBusPatch - 1;
	++cursor;
	begin[cursor] = end[cursor-1] + 1;

	cout << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	out << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] <<" " << iDDL-1 <<endl;
      }

     if (idDE  >= idCh6swp4 && idDE < 700) {
	
	if ( idDE ==  idCh6swp4) {
	  iDDL = (idDE/100)*2 ;
	  begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusCh5swp2+offsetBusCh6swp2, iDDL-1);
	}
	
	end[cursor]   = begin[cursor] - nbBusPatch + 1;
	++cursor;
	begin[cursor] = end[cursor-1] - 1;

	cout << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	out << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] <<" " << iDDL-1 <<endl;
      }

    }
 
    // station 4 & 5
    // back to normal one ddl connects one 1/2 chamber
    // finally !
 
    if (idDE >= 700) {

      if (idDE % 100 == 0) {
	cout << "# Chamber " << idDE/100 << endl;
	out  << "# Chamber " << idDE/100 << endl;
      }

      if ((idDE % 100) >= idSt45swp1 && (idDE % 100) < idSt45swp2) {
	if ((idDE % 100) == idSt45swp1) {
	  iDDL = (idDE/100)*2-1;
	  begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1);
	}
	end[cursor]   = begin[cursor] + nbBusPatch - 1;
	++cursor;
	begin[cursor] = end[cursor-1] + 1;

	cout << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	out << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] <<" " << iDDL-1 <<endl;
      }

      if (idDE % 100 >= idSt45swp2) {

	if ((idDE % 100) == idSt45swp2 ) {
	  iDDL = (idDE/100)*2 ;
	  if (idDE/100 == 7 || idDE/100 == 8)
	      begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusSt4swp, iDDL-1);
	  else
	      begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusSt5swp, iDDL-1);

	}
	end[cursor]   = begin[cursor] - nbBusPatch + 1;
	++cursor;
	begin[cursor] = end[cursor-1] - 1;

	cout << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	out << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] <<" " << iDDL-1 <<endl;
      }

      if ((idDE % 100) >= 0 && (idDE % 100) < idSt45swp1) {

	if ((idDE % 100) == 0 ) {
	  iDDL = (idDE/100)*2;
	  if (idDE/100 == 7 || idDE/100 == 8)
	      begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusSt4, iDDL-1);
	  else 
	      begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusSt5, iDDL-1);
	  
	}
	end[cursor]   = begin[cursor] - nbBusPatch + 1;
	++cursor;
	begin[cursor] = end[cursor-1] - 1;

	cout << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	out << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] <<" " << iDDL-1 <<endl;
      }
 
    }

  }
  printf(" Slat: number of TB %d and BB %d\n", nbTB, nbBB);
  in.close();
  out.close();  
}
Commit	Line	Data
	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
	18	#if !defined(__CINT__) \|\| defined(__MAKECINT__)
	19	// ROOT includes
	20	#include <Riostream.h>
	21	#include <TROOT.h>
	22	#include <TSystem.h>
	23	#include <TString.h>
	24
	25	// MUON includes
	26	#include "AliMpBusPatch.h"
	27
	28	#endif
	29
	30	// Generates buspatch id and DDL id for given detection element id
	31	// station 1 & 2 assuming 24 buspatches per quadrant
	32	// station345, reading from DetElemIdToSlatType.dat file and calculates
	33	// the number of bus patches per slat (and also number of Translator and Bridge Boards).
	34	// Generates an output file DetElemIdToBusPatch.dat.out, preserve from overwriting
	35	// (Ch. Finck, July 05)
	36	// (Nov. 05, added DDL)
	37	// (June 06, correction for St123)
	38	// (Feb. 07, add 1st manu list for St12 (starting on NB !) and new ddl sharing for station 3)
	39	// (June. 07, new numbering of station 345, and buspatch starts at 1)
	40
	41
	42	void MUONGenerateBusPatch()
	43	{
	44
	45	TString dirPath2 = gSystem->Getenv("ALICE_ROOT");
	46	dirPath2 += "/MUON/mapping/data/";
	47
	48	TString dirPath1 = dirPath2 + "station345/";
	49
	50	TString infile = dirPath1 + "DetElemIdToSlatType.dat";
	51	TString outfile = dirPath2 + "DetElemIdToBusPatch.dat.out";
	52
	53	ifstream in(infile, ios::in);
	54	ofstream out(outfile, ios::out);
	55	char line[80];
	56	char nameSlat[82];
	57	Int_t idDE;
	58	Int_t i;
	59	Int_t cursor = 0;
	60	Int_t begin[900];
	61	Int_t end[900];
	62	Int_t nbBusPatch;
	63	Int_t nbTB = 0;
	64	Int_t nbBB = 0;
	65	Int_t idSt12[] = {100, 101, 102, 103,
	66	200, 201, 202, 203,
	67	300, 301, 302, 303,
	68	400, 401, 402, 403};
	69
	70	Char_t manuListSt1[]
	71	= " 1,27,53,79,105,131,157,183,201,214,224,232,1025,1051,1077,1103,1129,1155,1181,1207,1225,1238,1249,1257";
	72
	73	Char_t manuListSt2[]
	74	= " 1,27,53,79,105,131,157,183,201,214,226,246,1025,1051,1077,1103,1129,1155,1181,1207,1225,1238,1251,1269";
	75
	76
	77
	78	Int_t iDDL = 0;
	79	// station 1 & 2
	80	nbBusPatch = 24;
	81	Int_t nbHalfBusPatch = nbBusPatch/2;
	82	cout << "#DE BusPatch DDL SlatName" << endl;
	83	out << "#DE BusPatch DDL 1st manu in buspatch" << endl;
	84
	85	for (Int_t j = 0; j < 16; j++) {
	86
	87	idDE = idSt12[j];
	88	if (idDE % 100 == 0) {
	89	iDDL++;
	90	begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1);
	91	cout << "# Chamber " << idDE/100 << endl;
	92	out << "# Chamber " << idDE/100 << endl;
	93	}
	94	if (idDE % 100 == 1) {
	95	iDDL++;
	96	begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1);
	97	}
	98
	99	if (idDE % 100 == 3) {
	100	iDDL--;
	101	begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1) + nbBusPatch;
	102	}
	103
	104	end[cursor] = begin[cursor] + nbBusPatch - 1;
	105	++cursor; // it seems that vec[++i] does not work in Cint
	106	begin[cursor] = end[cursor-1] + 1;
	107
	108	cout << idDE << " " << begin[cursor-1] + nbHalfBusPatch << "-" <<end[cursor-1] << ";" <<
	109	begin[cursor-1] << "-" << end[cursor-1] - nbHalfBusPatch << " " << iDDL-1 << endl;
	110
	111	if (idDE < 300 )
	112	out << idDE << " " << begin[cursor-1] + nbHalfBusPatch << "-" <<end[cursor-1] << ";" <<
	113	begin[cursor-1] << "-" << end[cursor-1] - nbHalfBusPatch << " " << iDDL-1 <<
	114	manuListSt1 << endl;
	115	else
	116	out << idDE << " " << begin[cursor-1] + nbHalfBusPatch << "-" <<end[cursor-1] << ";" <<
	117	begin[cursor-1] << "-" << end[cursor-1] - nbHalfBusPatch << " " << iDDL-1 <<
	118	manuListSt2 << endl;
	119	if (idDE % 100 == 3) iDDL++;
	120
	121	}
	122
	123	// station 345
	124	nbBusPatch = 0;
	125	Int_t idCh5swp1 = 501; // 1/4 chamber for DDL on horizontal
	126	Int_t idCh5swp2 = 505; // 1/4 chamber for DDL on horizontal
	127	Int_t idCh5swp3 = 510; // 1/4 chamber for DDL on horizontal
	128	Int_t idCh5swp4 = 514; // 1/4 chamber for DDL on horizontal
	129
	130	Int_t idCh6swp1 = 600; // 1/4 chamber for DDL on horizontal
	131	Int_t idCh6swp2 = 605; // 1/4 chamber for DDL on horizontal
	132	Int_t idCh6swp3 = 609; // 1/4 chamber for DDL on horizontal
	133	Int_t idCh6swp4 = 614; // 1/4 chamber for DDL on horizontal
	134
	135	Int_t idSt45swp1 = 7; // half chamber for DDL in vertical cutting twice the official numbering
	136	Int_t idSt45swp2 = 20;
	137
	138	Int_t offsetBusCh5swp1 = 13; // number of buspatches between DE 501-504
	139	Int_t offsetBusCh5swp2 = 17; // number of buspatches before DE 510
	140
	141	Int_t offsetBusCh6swp1 = 17; // number of buspatches before DE 600-604
	142	Int_t offsetBusCh6swp2 = 13; // number of buspatches between DE 614-617
	143
	144	Int_t offsetBusSt4 = 25; // number of buspatches before DE 700-800
	145	Int_t offsetBusSt4swp = 46; // number of buspatches between 801-806 (701-706)
	146
	147	Int_t offsetBusSt5 = 27; // number of buspatches before DE 900-1000
	148	Int_t offsetBusSt5swp = 50; // number of buspatches between 901-906 (1001-1006)
	149
	150	// reads from file
	151	while ( in.getline(line,80) ) {
	152
	153	if ( line[0] == '#' \|\| line[0] == '\0' ) continue;
	154	i = 0;
	155	sscanf(line, "%d %s", &idDE, nameSlat);
	156
	157	// number of PCB's
	158	while(nameSlat[i] < '4' && nameSlat[i] !='0') i++;
	159	Int_t len = strlen(nameSlat);
	160
	161	// number of buspatch
	162	if (i == 2 \|\| i == 3)
	163	nbBusPatch = 2;
	164	else
	165	nbBusPatch = 4;
	166
	167	// calculate the number of TB & BB
	168	// R1: 2 bridges more per PCB
	169	// R2 & R3: 1 translator more per PCB
	170	if (i == 2) nbBB += 2;
	171	if (i == 3) nbBB += 4;
	172	if (i == 4) nbBB += 4;
	173	if (i == 5) nbBB += 6;
	174	if (i == 6) nbBB += 8;
	175
	176	nbTB += nbBusPatch;
	177	if (nameSlat[len-1] == '1') {
	178	nbTB += 2;
	179	nbBusPatch+=2;
	180	}
	181	if (nameSlat[len-1] == '2' \|\| nameSlat[len-1] == '3') {
	182	nbTB += 1;
	183	nbBusPatch++;
	184	}
	185
	186	// station 3
	187	// for buspatch length reasons, one ddl connects one 1/4 of two chambers
	188	// really messy isn't it ?
	189	// much more with the new DDL sharing for station 3
	190
	191	if (idDE < 700 ) {
	192
	193	if (idDE % 100 == 0) {
	194	cout << "# Chamber " << idDE/100 << endl;
	195	out << "# Chamber " << idDE/100 << endl;
	196	}
	197
	198	// chamber 5
	199	if (idDE == 500 ) {
	200	iDDL = (idDE/100)*2+2;
	201	begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1);
	202
	203	end[cursor] = begin[cursor] + nbBusPatch - 1;
	204	++cursor;
	205	begin[cursor] = end[cursor-1] + 1;
	206
	207	cout << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	208	out << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] <<" " << iDDL-1 <<endl;
	209	}
	210
	211	if (idDE >= idCh5swp1 && idDE < idCh5swp2) {
	212
	213	if ( idDE == idCh5swp1) {
	214	iDDL = (idDE/100)*2 ;
	215	begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusCh5swp1, iDDL-1);
	216	}
	217
	218	end[cursor] = begin[cursor] - nbBusPatch + 1;
	219	++cursor;
	220	begin[cursor] = end[cursor-1] - 1;
	221
	222	cout << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	223	out << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] <<" " << iDDL-1 <<endl;
	224	}
	225
	226	if (idDE >= idCh5swp2 && idDE < idCh5swp3) {
	227
	228	if (idDE == idCh5swp2) {
	229	iDDL = (idDE/100)*2-1;
	230	begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1);
	231	}
	232
	233	end[cursor] = begin[cursor] + nbBusPatch - 1;
	234	++cursor;
	235	begin[cursor] = end[cursor-1] + 1;
	236
	237	cout << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	238	out << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] <<" " << iDDL-1 <<endl;
	239	}
	240
	241	if (idDE >= idCh5swp3 && idDE < idCh5swp4) {
	242
	243	if (idDE == idCh5swp3) {
	244	iDDL = (idDE/100)*2+1;
	245	begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1);
	246	}
	247
	248	end[cursor] = begin[cursor] + nbBusPatch - 1;
	249	++cursor;
	250	begin[cursor] = end[cursor-1] + 1;
	251
	252	cout << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	253	out << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] <<" " << iDDL-1 <<endl;
	254	}
	255
	256	if (idDE >= idCh5swp4 && idDE < 600) {
	257
	258	if ( idDE == idCh5swp4) {
	259	iDDL = (idDE/100)*2+2 ;
	260	begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusCh5swp2, iDDL-1);
	261	}
	262
	263	end[cursor] = begin[cursor] - nbBusPatch + 1;
	264	++cursor;
	265	begin[cursor] = end[cursor-1] - 1;
	266
	267	cout << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	268	out << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] <<" " << iDDL-1 <<endl;
	269	}
	270
	271	// chamber 6
	272	if (idDE >= idCh6swp1 && idDE < idCh6swp2) {
	273
	274	if ( idDE == idCh6swp1) {
	275	iDDL = (idDE/100)*2-2 ;
	276	begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusCh5swp1+offsetBusCh6swp1, iDDL-1);
	277	}
	278
	279	end[cursor] = begin[cursor] - nbBusPatch + 1;
	280	++cursor;
	281	begin[cursor] = end[cursor-1] - 1;
	282
	283	cout << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	284	out << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] <<" " << iDDL-1 <<endl;
	285	}
	286
	287	if (idDE >= idCh6swp2 && idDE < idCh6swp3) {
	288
	289	if (idDE == idCh6swp2) {
	290	iDDL = (idDE/100)*2-3;
	291	begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusCh5swp2+1, iDDL-1);
	292	}
	293
	294	end[cursor] = begin[cursor] + nbBusPatch - 1;
	295	++cursor;
	296	begin[cursor] = end[cursor-1] + 1;
	297
	298	cout << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	299	out << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] <<" " << iDDL-1 <<endl;
	300	}
	301
	302	if (idDE >= idCh6swp3 && idDE < idCh6swp4) {
	303
	304	if (idDE == idCh6swp3) {
	305	iDDL = (idDE/100)*2-1;
	306	begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusCh5swp1+1, iDDL-1);
	307	}
	308
	309	end[cursor] = begin[cursor] + nbBusPatch - 1;
	310	++cursor;
	311	begin[cursor] = end[cursor-1] + 1;
	312
	313	cout << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	314	out << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] <<" " << iDDL-1 <<endl;
	315	}
	316
	317	if (idDE >= idCh6swp4 && idDE < 700) {
	318
	319	if ( idDE == idCh6swp4) {
	320	iDDL = (idDE/100)*2 ;
	321	begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusCh5swp2+offsetBusCh6swp2, iDDL-1);
	322	}
	323
	324	end[cursor] = begin[cursor] - nbBusPatch + 1;
	325	++cursor;
	326	begin[cursor] = end[cursor-1] - 1;
	327
	328	cout << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	329	out << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] <<" " << iDDL-1 <<endl;
	330	}
	331
	332	}
	333
	334	// station 4 & 5
	335	// back to normal one ddl connects one 1/2 chamber
	336	// finally !
	337
	338	if (idDE >= 700) {
	339
	340	if (idDE % 100 == 0) {
	341	cout << "# Chamber " << idDE/100 << endl;
	342	out << "# Chamber " << idDE/100 << endl;
	343	}
	344
	345	if ((idDE % 100) >= idSt45swp1 && (idDE % 100) < idSt45swp2) {
	346	if ((idDE % 100) == idSt45swp1) {
	347	iDDL = (idDE/100)*2-1;
	348	begin[cursor] = AliMpBusPatch::GetGlobalBusID(1, iDDL-1);
	349	}
	350	end[cursor] = begin[cursor] + nbBusPatch - 1;
	351	++cursor;
	352	begin[cursor] = end[cursor-1] + 1;
	353
	354	cout << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	355	out << idDE << " " << begin[cursor-1]<<"-"<<end[cursor-1] <<" " << iDDL-1 <<endl;
	356	}
	357
	358	if (idDE % 100 >= idSt45swp2) {
	359
	360	if ((idDE % 100) == idSt45swp2 ) {
	361	iDDL = (idDE/100)*2 ;
	362	if (idDE/100 == 7 \|\| idDE/100 == 8)
	363	begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusSt4swp, iDDL-1);
	364	else
	365	begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusSt5swp, iDDL-1);
	366
	367	}
	368	end[cursor] = begin[cursor] - nbBusPatch + 1;
	369	++cursor;
	370	begin[cursor] = end[cursor-1] - 1;
	371
	372	cout << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	373	out << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] <<" " << iDDL-1 <<endl;
	374	}
	375
	376	if ((idDE % 100) >= 0 && (idDE % 100) < idSt45swp1) {
	377
	378	if ((idDE % 100) == 0 ) {
	379	iDDL = (idDE/100)*2;
	380	if (idDE/100 == 7 \|\| idDE/100 == 8)
	381	begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusSt4, iDDL-1);
	382	else
	383	begin[cursor] = AliMpBusPatch::GetGlobalBusID(offsetBusSt5, iDDL-1);
	384
	385	}
	386	end[cursor] = begin[cursor] - nbBusPatch + 1;
	387	++cursor;
	388	begin[cursor] = end[cursor-1] - 1;
	389
	390	cout << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] << " " << iDDL-1 << " " <<nameSlat <<endl;
	391	out << idDE << " " << end[cursor-1]<<"-"<<begin[cursor-1] <<" " << iDDL-1 <<endl;
	392	}
	393
	394	}
	395
	396	}
	397	printf(" Slat: number of TB %d and BB %d\n", nbTB, nbBB);
	398	in.close();
	399	out.close();
	400	}