PndFTSCAParam.cxx File Reference

#include "PndFTSCAParam.h"
#include "PndFTSCAMath.h"
#include <iostream>
#include "FairRunAna.h"
#include "FairField.h"
#include "debug.h"
#include "TVector3.h"
#include "TVectorD.h"
#include "TMatrixD.h"
#include "TH2F.h"
#include "BinaryStoreHelper.h"

Go to the source code of this file.

Functions
std::istream &	operator>> (std::istream &in, PndFTSCAParam &p)

Function Documentation

std::istream& operator>>	(	std::istream &	in,
		PndFTSCAParam &	p
	)

Definition at line 159 of file PndFTSCAParam.cxx.

References PndFTSCAParam::cBz(), FTSCAStation::CellLength, FTSCAStripInfo::cos, cos(), exit(), FTSCAStation::f, PndFTSCAParam::fBz, PndFTSCAParam::fMaxR, PndFTSCAParam::fMaxZ, PndFTSCAParam::fMinR, PndFTSCAParam::fMinZ, PndFTSCAParam::fNStations, PndFTSCAParam::fStations, PndFTSCAParam::fVtxFieldValue, i, log(), FTSCAStation::NDF, FTSCAStripInfo::sin, sin(), FTSCAStation::x0, FTSCAStation::xOverX0, and FTSCAStation::xTimesRho.

{
  // Read settings from the file
  in >> p.fNStations;
  if(p.fStations) delete [] p.fStations;
  p.fStations = new FTSCAStation[p.fNStations];
  in >> p.fBz;

  for(int i=0; i<p.fNStations; i++)
  {
    float inttmp;
    in >> inttmp;
    if( inttmp != i ) {
      cout << "ERROR: Settings.data format is wrong! Station: " << i-1 << endl;
      exit(0);
    }
    in >> p.fStations[i].x0;
#ifdef PANDA_FTS
    in >> p.fStations[i].materialInfo.RadThick;
    p.fStations[i].materialInfo.RadThick *= 10.f; //TODO !!!!!!! put correct material
    p.fStations[i].materialInfo.logRadThick = log( p.fStations[i].materialInfo.RadThick );
    float xTimesRho; 
    in >> xTimesRho;
    p.fStations[i].materialInfo.thick = xTimesRho/1.39; // ~54 mum mular
    p.fStations[i].materialInfo.RL = p.fStations[i].materialInfo.thick/p.fStations[i].materialInfo.RadThick;
#else
    in >> p.fStations[i].xOverX0;
    in >> p.fStations[i].xTimesRho;
#endif
    
    float beta;
    in >> beta;
    p.fStations[i].f.sin = sin(beta);
    p.fStations[i].f.cos = cos(beta);
    //std::cout<<"p.fStations[i].f.sin "<<p.fStations[i].f.sin<<" p.fStations[i].f.cos "<<p.fStations[i].f.cos<<std::endl;
    in >> inttmp;
    p.fStations[i].NDF = inttmp;
    in >> inttmp;
    p.fStations[i].CellLength = inttmp;
  }

#ifdef STAR_HFT
  p.fMaxZ = 35;
  p.fMaxR = 25;
  p.fMinR = 0;
  p.fMinZ = -p.fMaxZ;
#elif ALICE_ITS
  p.fMaxZ = 60;
  p.fMaxR = 60;
  p.fMinR = 0;
  p.fMinZ = -p.fMaxZ;
#elif PANDA_STT
  p.fMaxZ = 75+20;
  p.fMaxR = 41;
  p.fMinR = 0;
  p.fMinZ = -75+20;
#elif PANDA_FTS
  p.fMinX = -200; // more precisely 196
  p.fMaxX =  200;
  p.fMinY = -60;
  p.fMaxY = 60;
  p.fMinZ = 290;
  p.fMaxZ = 660;
#endif

#ifdef PANDA_FTS
  // p.fVtxFieldValue[0].y = p.Bz();
  // p.fVtxFieldValue[1].y = p.Bz();
  // p.fZVtxFieldValue[0] = 0.f;
  // p.fZVtxFieldValue[1] = 1.f;
  // for(int i=0; i<p.fNStations; i++)
  //   p.fStations[i].fieldSlice.cy[0] = p.Bz();
//   if(1){ // read field
//     float tmp;
//     for( int i = 0; i < 2; i++ ) {
//       in >> tmp;
//       p.fZVtxFieldValue[i] = tmp;
//       in >> tmp;
//       p.fVtxFieldValue[i].x = tmp;
//       in >> tmp;
//       p.fVtxFieldValue[i].y = tmp;
//       in >> tmp;
//       p.fVtxFieldValue[i].z = tmp;
//     }
// 
//     in >> tmp >> tmp >> tmp >> tmp; // skip one point
//       
//     for(int iSta = 0; iSta < p.fNStations; iSta++) {
// #if 0
//       int N;
//       in >> N;
//       for( int i = 0; i < N; i++ ) {
//         in >> tmp;
//         p.fStations[i].fieldSlice.cx[i] = tmp;
//       }
//       for( int i = 0; i < N; i++ ) {
//         in >> tmp;
//         p.fStations[i].fieldSlice.cy[i] = tmp;
//       }
//       for( int i = 0; i < N; i++ ) {
//         in >> tmp;
//         p.fStations[i].fieldSlice.cz[i] = tmp;
//       }
// #else // 0
//       L1FieldSlice &sl = p.fStations[iSta].fieldSlice;
//       vector<L1FieldSlice> &vSls = p.fStations[iSta].fieldVirtualSlice;
//       int itmp;
//       in >> itmp; // station
//       if ( itmp != iSta ) { cout << "Wrong field" << endl; exit(0); }
//       ASSERT( itmp == iSta, itmp << " = " << iSta );
//       in >> itmp >> p.fStations[iSta].dZVirtualStation; // n virtual stations & dz
//       vSls.resize(itmp);
//       in >> sl.xMin >> sl.dx >> sl.nXBins;
//       in >> sl.yMin >> sl.dy >> sl.nYBins;
//       const int NBins = sl.nXBins*sl.nYBins;
//       for( unsigned int iV = 0; iV < vSls.size(); iV++ ) {
//         L1FieldSlice &vSl = vSls[iV];
//         vSl.xMin = sl.xMin;
//         vSl.dx = sl.dx;
//         vSl.nXBins = sl.nXBins;
//         vSl.yMin = sl.yMin;
//         vSl.dy = sl.dy;
//         vSl.nYBins = sl.nYBins;
//         vSl.AlocateMemory();
//         for( int iB = 0; iB < NBins; iB++ ) {
//           in >> vSl.bX[iB] >> vSl.bY[iB] >> vSl.bZ[iB];
//         }
//       }
// 
//       sl.AlocateMemory();
//       for( int iB = 0; iB < NBins; iB++ ) {
//         in >> sl.bX[iB] >> sl.bY[iB] >> sl.bZ[iB];
//       }
// #endif // 0
//     }
//   }
#else //  PANDA_FTS
  p.fVtxFieldValue[0] = p.cBz();
#endif
  
  return in;
}