ecalInf Class Reference

#include <ecalInf.h>

Inheritance diagram for ecalInf:

Collaboration diagram for ecalInf:

Public Member Functions
	ecalInf ()
Double_t	GetXPos () const
Double_t	GetYPos () const
Double_t	GetZPos () const
Double_t	GetModuleSize () const
Int_t	GetNLayers () const
Double_t	GetLead () const
Double_t	GetScin () const
Double_t	GetTyveec () const
Double_t	GetThicknessLayer () const
Int_t	GetXSize () const
Int_t	GetYSize () const
Double_t	GetContainerXSemiAxiss () const
Double_t	GetContainerYSemiAxiss () const
Double_t	GetElectronCut () const
Double_t	GetHadronCut () const
Double_t	GetEcalSize (Int_t num) const
char	GetType (Int_t x, Int_t y) const
Int_t	GetFastMC () const
void	DumpContainer () const
void	FreeInstance ()
Double_t	GetVariable (const char *key)
Double_t	GetVariableStrict (const char *key)
TString	GetStringVariable (const char *key)
void	AddVariable (const char *key, const char *value)

Static Public Member Functions
static ecalInf *	GetInstance (const char *filename)

Protected Member Functions
	ecalInf (const char *filename)
void	CalculateHoleSize ()
virtual	~ecalInf ()

Static Protected Attributes
static ecalInf *	fInf =NULL
static Int_t	fRefCount =0

Private Member Functions
void	InitVariables ()
Bool_t	ExcludeParameter (TString parname)
	ecalInf (const ecalInf &)
ecalInf &	operator= (const ecalInf &)
	ClassDef (ecalInf, 2)

Private Attributes
TMap *	fVariables
TObjArray	fEcalStr
Double_t	fXPos
Double_t	fYPos
Double_t	fZPos
Int_t	fNLayers
Int_t	fXSize
Int_t	fYSize
Double_t	fModuleSize
Double_t	fLead
Double_t	fScin
Double_t	fTyveec
Double_t	fThicknessLayer
Double_t	fCellSize
Double_t	fEcalSize [3]
Double_t	fECut
Double_t	fHCut
Double_t	fSemiX
Double_t	fSemiY
Int_t	fFastMC
Int_t	fSuccess
TString	fFileName

Detailed Description

Definition at line 19 of file ecalInf.h.

Constructor & Destructor Documentation

ecalInf() [1/3]

ecalInf::ecalInf ( )

inline

This is ROOT requirement

Definition at line 23 of file ecalInf.h.

           : TObject(), fVariables(NULL), fEcalStr(), fXPos(0.), fYPos(0.),
    fZPos(0.), fNLayers(0), fXSize(0), fYSize(0), fModuleSize(0.), fLead(0.),
    fScin(0.), fTyveec(0.), fThicknessLayer(0.), fCellSize(0.), fEcalSize(),
    fECut(0.), fHCut(0.), fSemiX(0.0), fSemiY(0.0), fFastMC(-1), 
    fSuccess(-1), fFileName("") 
    {};

ecalInf() [2/3]

ecalInf::ecalInf ( const char *  filename )

protected

If data from Parameter file differs from ours TODO: should understand a way of storring parameters in SHIP Text file constructor

Constructor reads and parses the ascii file, and fill 
the ECAL geometry container

FairRunAna* ana = FairRunAna::Instance(); if (ana!=NULL) { FairRuntimeDb* rtdb=ana->GetRuntimeDb(); rtdb->getContainer("CbmGeoEcalPar"); }

Definition at line 139 of file ecalInf.cxx.

  : TObject(), 
    fVariables(new TMap(200)), 
    fEcalStr(), 
    fXPos(0.), 
    fYPos(0.),
    fZPos(0.), 
    fNLayers(0), 
    fXSize(0), 
    fYSize(0), 
    fModuleSize(0.), 
    fLead(0.), 
    fScin(0.), 
    fTyveec(0.), 
    fThicknessLayer(0.),
    fCellSize(0.), 
    fEcalSize(), 
    fECut(0.), 
    fHCut(0.), 
    fSemiX(0.),
    fSemiY(0.),
    fFastMC(-1),
    fSuccess(1), 
    fFileName(filename) 
{
  std::ifstream file(filename);
  Int_t linenum;
  Double_t val;
  string buffer;
  string message;
  string variable;
  string value;
  TObjString* str=NULL;
  char** err=NULL;
  char winend[2]={13, 0};
  int ssize=-1;
 
  if (!file) {
    cerr << "ecalInf: Can't open information file " << filename << "!" << endl;
    cerr << "ecalInf: Ecal will not be constructed correctly." << endl;
    fSuccess=0;
    return;
  }
 
  linenum=0;
  while(getline(file,buffer)) {
    linenum++;
    message=buffer.substr(buffer.find_first_not_of("    "));    //Skiping initial spaces
    message=message.substr(0,message.find("#"));    //Removing comments
    // Threat windows end of strings correctly
    message=message.substr(0,message.find(winend));
    if (message.empty()) continue;      //This was just a comment
    variable=message.substr(0,message.find("="));
    if (variable=="structure") {
      while(getline(file,buffer)) {
    linenum++;
    if (buffer.empty()) break;
    message=buffer.substr(buffer.find_first_not_of("    "));    //Skiping initial spaces
    message=message.substr(0,message.find("#"));    //Removing comments
    message=message.substr(0,message.find_last_not_of("     ")+1);  //Skiping ending spaces
 
        // Threat windows end of strings correctly
        message=message.substr(0,message.find(winend));
 
    if (!message.empty()) {
      if (-1==ssize)
        ssize=message.size();
      else
        if (ssize!=(Int_t)message.size()) {
          cerr << "Error in ECAL structure at line " << linenum;
          cerr << "." << endl;
          cerr << "Line length differs from previous one" << endl;
              fSuccess=0;
          file.close();
          return;
          
        }
      
      str=new TObjString(message.c_str());
      fEcalStr.Add(str);
    }
      }
      break;
    }
    if (variable==message) {
      cerr << "Syntax error: File " << filename << ".Line " << linenum << "." << endl;
      fSuccess=0;
      file.close();
      return;
    }
    variable=variable.substr(0,variable.find_first_of("     "));
    value=message.substr(message.find("=")+1);
    value=value.substr(value.find_first_not_of("    "));    //Skiping initial spaces
    value=value.substr(0,value.find_first_of("  "));
/*
    value=value.substr(0,value.find_first_not_of("1234567890-+e."));
    val=strtod(value.c_str(),err);
    if (err) {
      cerr << "Syntax error after =: File " << filename << ".Line " << linenum << "." << endl;
      fSuccess=0;
      file.close();
      return;
    }
*/
    AddVariable(variable.c_str(), value.c_str());
  }
  file.close();
  InitVariables();
}

~ecalInf()

ecalInf::~ecalInf ( )

protectedvirtual

Definition at line 28 of file ecalInf.cxx.

{
    for(Int_t i=0;i<fEcalStr.GetSize();i++)
        delete (TObjString*)(fEcalStr.At(i));
    fEcalStr.Clear();
}

ecalInf() [3/3]

ecalInf::ecalInf ( const ecalInf &  )

private

Member Function Documentation

AddVariable()

void ecalInf::AddVariable	(	const char *	key,
		const char *	value
	)

Definition at line 128 of file ecalInf.cxx.

{
  TObjString* skey=(TObjString*)fVariables->FindObject(key);
  //Value for this key already exists!!!
  if (skey!=NULL) return;
  skey=new TObjString(key);
  skey->String().ToLower();
  TObjString* svalue=new TObjString(value);
  fVariables->Add(skey, svalue);
}

CalculateHoleSize()

void ecalInf::CalculateHoleSize ( )

protected

ClassDef()

ecalInf::ClassDef ( ecalInf  , 2    )

private

DumpContainer()

void ecalInf::DumpContainer

(

)

const

Definition at line 361 of file ecalInf.cxx.

{
  // Print out the ECAL geometry container
 
  if (fVariables)
  {
    TObjString* key;
    TIterator* iter=fVariables->MakeIterator();
    while((key=(TObjString*)iter->Next())!=NULL)
    {
      TObjString* str=(TObjString*)fVariables->GetValue(key);
      cout << key->String() << "=" << str->String() << endl;
    }
  }
  TObjString* key;
  TIterator* iter=fEcalStr.MakeIterator();
 
  Int_t modules=0;
  Int_t channels=0;
  Int_t i;
  Int_t j;
  Int_t m[10];
  char stri[2]={0, 0};
  TString st;
  for(i=0;i<10;i++) m[i]=0;
 
  while((key=(TObjString*)iter->Next())!=NULL)
  {
    st=key->String();
    cout << key->String() << endl;
    for(i=0;i<st.Length();i++)
    {
      stri[0]=st[i];
      j=atoi(stri);
      m[j]++;
      if (j) modules++;
      channels+=j*j;
    }
  }
  cout << "Total modules:  " << modules << endl;
  cout << "Total channels: " << channels << endl;
  for(i=1;i<10;i++)
  {
    if (m[i]==0) continue;
    cout << "   Type " << i << " : modules=" << m[i] << ", channels=" << m[i]*i*i << endl;
  }
 
}

ExcludeParameter()

Bool_t ecalInf::ExcludeParameter ( TString  parname )

private

Ignore a parameter during comparision

Definition at line 260 of file ecalInf.cxx.

{
  if (parname.CompareTo("ecalversion")==0) return kTRUE;
  return kFALSE;
}

FreeInstance()

void ecalInf::FreeInstance ( )

inline

Definition at line 165 of file ecalInf.h.

{
    fRefCount--;
    if (fRefCount==0)
    {
        delete this;
        fInf=NULL;
    }
}

GetContainerXSemiAxiss()

Double_t ecalInf::GetContainerXSemiAxiss ( ) const

inline

Definition at line 47 of file ecalInf.h.

{return fSemiX;}

GetContainerYSemiAxiss()

Double_t ecalInf::GetContainerYSemiAxiss ( ) const

inline

Definition at line 48 of file ecalInf.h.

{return fSemiY;}

GetEcalSize()

Double_t ecalInf::GetEcalSize ( Int_t  num ) const

inline

Definition at line 54 of file ecalInf.h.

                                               {
    if (num>-1&&num<3) return fEcalSize[num];
    return -1;
  }

GetElectronCut()

Double_t ecalInf::GetElectronCut ( ) const

inline

Geant cuts information

Definition at line 51 of file ecalInf.h.

{return fECut;}

GetFastMC()

Int_t ecalInf::GetFastMC ( ) const

inline

Definition at line 59 of file ecalInf.h.

{return fFastMC;}

GetHadronCut()

Double_t ecalInf::GetHadronCut ( ) const

inline

Definition at line 52 of file ecalInf.h.

{return fHCut;}

GetInstance()

ecalInf * ecalInf::GetInstance ( const char *  filename )

static

Definition at line 35 of file ecalInf.cxx.

{
    if (filename==NULL)
    {
        if (fInf!=NULL)
            fRefCount++;
        return fInf;
    }
    TString newname=gSystem->Getenv("VMCWORKDIR");
    newname+="/geometry/";
    newname+=filename;
    if (fInf!=NULL) {
        if (fInf->fFileName==newname) {
            fRefCount++;
            return fInf;
        } else {
          cerr << "ecalInf: Trying create ";
          cerr << "instance of ecalInf with";
          cerr << " name " << filename;
          cerr << ", which is different from ";
          cerr << fInf->fFileName << "." << endl;
          return NULL;
        }
    }
    fInf=new ecalInf(newname);
    //Is something wrong?
    if (fInf->fSuccess==0)
    {
        delete fInf;
        return NULL;
    }
    fRefCount++;
    return fInf;
}

GetLead()

Double_t ecalInf::GetLead ( ) const

inline

Definition at line 39 of file ecalInf.h.

{return fLead;}

GetModuleSize()

Double_t ecalInf::GetModuleSize ( ) const

inline

Definition at line 36 of file ecalInf.h.

{return fModuleSize;}

GetNLayers()

Int_t ecalInf::GetNLayers ( ) const

inline

Definition at line 38 of file ecalInf.h.

{return fNLayers;}

GetScin()

Double_t ecalInf::GetScin ( ) const

inline

Definition at line 40 of file ecalInf.h.

{return fScin;}

GetStringVariable()

TString ecalInf::GetStringVariable

(

const char *

key

)

Definition at line 103 of file ecalInf.cxx.

{
  TObjString* value=(TObjString*)fVariables->GetValue(key);
  if (value==NULL)
  {
    Fatal("GetStringVariable","Can't find variable named %s.", key);
  }
  return value->GetString();
}

GetThicknessLayer()

Double_t ecalInf::GetThicknessLayer ( ) const

inline

Definition at line 42 of file ecalInf.h.

{return fThicknessLayer;}

GetType()

char ecalInf::GetType ( Int_t  x, Int_t  y  ) const

inline

Returns the type of the module with position (x,y)

Definition at line 155 of file ecalInf.h.

{
  if (x<0||y<0||y>fEcalStr.GetLast()) return 0;
  TObjString* str=(TObjString*)fEcalStr.At(y);
  if (str->GetString().Length()<x) return 0;
  char stri[2]={str->GetString()[x],0};
  return atoi(stri);
}

GetTyveec()

Double_t ecalInf::GetTyveec ( ) const

inline

Definition at line 41 of file ecalInf.h.

{return fTyveec;}

GetVariable()

Double_t ecalInf::GetVariable

(

const char *

key

)

key must be lower case. For example, if have in geo file AaaA=90, then you should call GetVariableStrict("aaaa"). If variable not found, will return -1111

Definition at line 115 of file ecalInf.cxx.

{
  TObjString* value=(TObjString*)fVariables->GetValue(key);
  if (value==NULL)
    return -1111;
  Double_t val;
  char* err=NULL;
  val=strtod(value->GetString(),&err);
  if (err[0]!='\0')
    return -1111;
  return val;
}

GetVariableStrict()

Double_t ecalInf::GetVariableStrict

(

const char *

key

)

key must be lower case. For example, if have in geo file AaaA=90, then you should call GetVariableStrict("aaaa"). If variable not found, will generate Fatal

Definition at line 82 of file ecalInf.cxx.

{
  TObjString* value=(TObjString*)fVariables->GetValue(key);
  if (value==NULL)
  {
    cerr << "Can't find variable named \"" << key << "\"";
    Fatal("GetVariableStrict","Exiting...");
  }
  Double_t val;
  char* err=NULL;
  val=strtod(value->GetString(),&err);
  if (err[0]!='\0')
  {
    cerr << "Can't convert variable named \"" << key ;
    cerr << "\" to floating point. Value is \"";
    cerr << value->GetString() << "\"." << endl;
    Fatal("GetVariableStrict","Exiting...");
  }
  return val;
}

GetXPos()

Double_t ecalInf::GetXPos ( ) const

inline

Getters

Definition at line 32 of file ecalInf.h.

{return fXPos;} 

GetXSize()

Int_t ecalInf::GetXSize ( ) const

inline

Size of ECAL in super modules

Definition at line 45 of file ecalInf.h.

{return fXSize;}

GetYPos()

Double_t ecalInf::GetYPos ( ) const

inline

Definition at line 33 of file ecalInf.h.

{return fYPos;} 

GetYSize()

Int_t ecalInf::GetYSize ( ) const

inline

Definition at line 46 of file ecalInf.h.

{return fYSize;}

GetZPos()

Double_t ecalInf::GetZPos ( ) const

inline

Definition at line 34 of file ecalInf.h.

{return fZPos;}

InitVariables()

void ecalInf::InitVariables ( )

private

Init all other variables from fVariables and fEcalStr

Definition at line 324 of file ecalInf.cxx.

{
  TString stri;
  TObjString* str=(TObjString*)fEcalStr.At(0);
 
  fXPos=GetVariableStrict("xpos");
  fYPos=GetVariableStrict("ypos");
  fZPos=GetVariableStrict("zpos");
  fNLayers=(Int_t)GetVariableStrict("nlayers");
  fModuleSize=GetVariableStrict("modulesize");
  fLead=GetVariableStrict("lead");
  fScin=GetVariableStrict("scin");
  fTyveec=GetVariableStrict("tyveec");
  fEcalSize[2]=GetVariableStrict("ecalzsize");
  fECut=GetVariableStrict("ecut");
  fHCut=GetVariableStrict("hcut");
  fFastMC=(Int_t)GetVariableStrict("fastmc");
  fSemiX=GetVariableStrict("xsemiaxis");
  fSemiY=GetVariableStrict("ysemiaxis");
 
  stri=""; stri+=str->GetString().Length();
  AddVariable("xsize", stri);
  stri=""; stri+=fEcalStr.GetLast()+1;
  AddVariable("ysize", stri);
 
  fXSize=(Int_t)GetVariableStrict("xsize");
  fYSize=(Int_t)GetVariableStrict("ysize");
  fThicknessLayer = 2.0*GetTyveec()+GetScin()+GetLead();
  fEcalSize[0]    = GetXSize()*GetModuleSize();
  fEcalSize[1]    = GetYSize()*GetModuleSize();
  stri=""; stri+=fEcalSize[0];
  AddVariable("xecalsize",stri);
  stri=""; stri+=fEcalSize[1];
  AddVariable("yecalsize",stri);
}

operator=()

ecalInf & ecalInf::operator= ( const ecalInf &  )

private

Member Data Documentation

fCellSize

Double_t ecalInf::fCellSize

private

transverse size of ECAL cell for simulation [cm]

Definition at line 130 of file ecalInf.h.

fEcalSize

Double_t ecalInf::fEcalSize[3]

private

Size of ECAL container [cm]

Definition at line 132 of file ecalInf.h.

fEcalStr

TObjArray ecalInf::fEcalStr

private

Structure of ECAL as array of strings This variable should be saved in parameter file

Definition at line 102 of file ecalInf.h.

fECut

Double_t ecalInf::fECut

private

Electron cut for ECAL

Definition at line 134 of file ecalInf.h.

fFastMC

Int_t ecalInf::fFastMC

private

Flag to run Fast (1) or Full (0) MC code

Definition at line 142 of file ecalInf.h.

fFileName

TString ecalInf::fFileName

private

Definition at line 147 of file ecalInf.h.

fHCut

Double_t ecalInf::fHCut

private

Hadron cut for ECAL

Definition at line 136 of file ecalInf.h.

fInf

ecalInf * ecalInf::fInf =NULL

staticprotected

Definition at line 87 of file ecalInf.h.

fLead

Double_t ecalInf::fLead

private

Thickness of lead in one layer [cm]

Definition at line 122 of file ecalInf.h.

fModuleSize

Double_t ecalInf::fModuleSize

private

transverse supermodule size in cm

Definition at line 119 of file ecalInf.h.

fNLayers

Int_t ecalInf::fNLayers

private

Number of ECAL layers

Definition at line 111 of file ecalInf.h.

fRefCount

Int_t ecalInf::fRefCount =0

staticprotected

Definition at line 88 of file ecalInf.h.

fScin

Double_t ecalInf::fScin

private

Thickness of scintillator in one layer [cm]

Definition at line 124 of file ecalInf.h.

fSemiX

Double_t ecalInf::fSemiX

private

Semiaxises keeping volume

Definition at line 138 of file ecalInf.h.

fSemiY

Double_t ecalInf::fSemiY

private

Definition at line 139 of file ecalInf.h.

fSuccess

Int_t ecalInf::fSuccess

private

1 if evething Ok

Definition at line 145 of file ecalInf.h.

fThicknessLayer

Double_t ecalInf::fThicknessLayer

private

Total thickness of layer [cm]

Definition at line 128 of file ecalInf.h.

fTyveec

Double_t ecalInf::fTyveec

private

Thickness of tyvec in one layer [cm]

Definition at line 126 of file ecalInf.h.

fVariables

TMap* ecalInf::fVariables

private

A map containing all variables This variable should be saved in parameter file

Definition at line 99 of file ecalInf.h.

fXPos

Double_t ecalInf::fXPos

private

x-position of ECAL center [cm]

Definition at line 104 of file ecalInf.h.

fXSize

Int_t ecalInf::fXSize

private

x-size of ecal in supermodules

Definition at line 114 of file ecalInf.h.

fYPos

Double_t ecalInf::fYPos

private

y-position of ECAL center [cm]

Definition at line 106 of file ecalInf.h.

fYSize

Int_t ecalInf::fYSize

private

y-size of ecal in supermodules

Definition at line 116 of file ecalInf.h.

fZPos

Double_t ecalInf::fZPos

private

z-position of ECAL front [cm]

Definition at line 108 of file ecalInf.h.

---

The documentation for this class was generated from the following files:

• ecal/ecalInf.h
• ecal/ecalInf.cxx