Atlfast::EPileupMap Class Reference

Simple std::vector containing EPileupDeposit. More...

#include <EPileupMap.h>

Collaboration diagram for Atlfast::EPileupMap:

[legend]List of all members.

Public Member Functions
	EPileupMap (int)
	~EPileupMap ()
void	fillMap ()
std::vector< EPileupDeposit * > *	getEMMap () const
std::vector< EPileupDeposit * > *	getHadMap () const
Private Member Functions
void	makeCellMap ()
void	setInputFile ()
void	emptyMap ()
Private Attributes
std::vector< EPileupDeposit * > *	m_hadPileupMap
std::vector< EPileupDeposit * > *	m_emPileupMap
std::map< int, double > *	m_phiMap
std::map< int, double > *	m_etaMap
std::string	m_pileupFile
std::string	m_fileDir
int	m_nFiles

---

Detailed Description

Simple std::vector containing EPileupDeposit.

Can randomly select a file to read in and fill vector with EPileupDeposits calculating its eta and phi from the code given.

Author:

Jon Couchman

Definition at line 40 of file EPileupMap.h.

---

Constructor & Destructor Documentation

Atlfast::EPileupMap::EPileupMap

(

int

)

Constructor

Set input epileup file directory and number of files in directory depending on the luminosity

Construnt maps relating the cell ids to etas and phis

Definition at line 13 of file EPileupMap.cxx.

                                      {
    m_emPileupMap = new std::vector<EPileupDeposit*>;
    m_hadPileupMap = new std::vector<EPileupDeposit*>;
    
    if (luminosity == 1){
      m_nFiles = 10401;
      m_fileDir = "allfiles_ll/";
    }else if (luminosity == 2){
      m_nFiles = 10830;
      m_fileDir = "allfiles_hl/";
    }
    this->makeCellMap();
  }

Atlfast::EPileupMap::~EPileupMap

(

)

Destructor

Definition at line 30 of file EPileupMap.cxx.

                         {
    this->emptyMap();
    delete m_emPileupMap;
    delete m_hadPileupMap;
    delete m_phiMap;
    delete m_etaMap;
  }

---

Member Function Documentation

void Atlfast::EPileupMap::fillMap

(

)

Fills vector m_map with new EPileupDeposits from a text file. Deletes any deposits in vector before it fills it with the new ones

Clear map and delete all the EpileupDeposits in it

Set m_pileupFile using random nember generator

Open input file

find whether EM or HAD deposit

Set eta and phi codes from single code read in

Convert eT from MeV to GeV

Put new EPileupDeposit in EM vector

Put new EPileupDeposit in HAD vector

Definition at line 39 of file EPileupMap.cxx.

                          {
    this->emptyMap();
    this->setInputFile();

    std::ifstream input;
    input.open( m_pileupFile.c_str());
    if (input) {
      int code;
      double energy;
      while (input >> code && input >> energy){
        if (code > 0){ 
          std::string type;
          int intType = code/100000;
          int eta = (code - intType*100000)/100;
          int phi =  code - intType*100000 - eta*100;
          double et = energy*GeV;

          if(intType == 1){
            m_emPileupMap->
              push_back(
                        new EPileupDeposit(
                                           (m_etaMap->find(eta))->second,
                                           (m_phiMap->find(phi))->second,
                                           et)
                        );
          }//if(intType == 1)
          
          if(intType == 2){
            m_hadPileupMap->
              push_back(
                        new EPileupDeposit(
                                           (m_etaMap->find(eta))->second,
                                           (m_phiMap->find(phi))->second,
                                           et
                                           )
                        );
          }//if(intType == 2)
          
        }
        
      }
    }else{
      std::cout << "ERROR! COULD NOT OPEN EPILEUP FILE " 
                << m_pileupFile << std::endl;
    } 
    return;
  }

std::vector<EPileupDeposit*>* Atlfast::EPileupMap::getEMMap

(

)

const [inline]

Returns std::vector of EPileupDeposits

Definition at line 50 of file EPileupMap.h.

{return m_emPileupMap;}

std::vector<EPileupDeposit*>* Atlfast::EPileupMap::getHadMap

(

)

const [inline]

Returns std::vector of EPileupDeposits

Definition at line 52 of file EPileupMap.h.

{return m_hadPileupMap;}

void Atlfast::EPileupMap::makeCellMap

(

)

[private]

Makes two maps relating cell ids to eta and phi positions

Maybe these should be setable member variables. As we only have one set of maps at the moment, it would be dangerous for user to fiddle with these

Definition at line 113 of file EPileupMap.cxx.

                              {
    m_phiMap = new std::map<int,double> ;
    m_etaMap = new std::map<int,double> ;
    int EtaCells = 100;
    double EtaRange = 10;
    int PhiCells = 64;
    double PhiRange = 2*3.1415;
    for (int i = 1 ; i <= EtaCells; ++i){
      m_etaMap->
        insert(std::pair<int,double>(i,
                                     (i-(EtaCells/2))/(EtaCells/EtaRange)-
                                      EtaRange/(EtaCells*2.0) 
                                     )
               );
    }
    for (int i = 1 ; i <= PhiCells; ++i){
      m_phiMap->
        insert(std::pair<int,double>(i,
                                     (i-(PhiCells/2))/(PhiCells/PhiRange)-
                                     PhiRange/(PhiCells*2.0)
                                     )
               ) ;
    }
    return;
  }

void Atlfast::EPileupMap::setInputFile

(

)

[private]

Uses a random number to set the value of m_pileupFile

Definition at line 144 of file EPileupMap.cxx.

                               {
    double max = RAND_MAX; 
    int i = (int)(rand()/max*m_nFiles);
    if(i < 10){
      m_pileupFile = m_fileDir + "part0000" + numberToString(i);
     } else if(i < 100){
       m_pileupFile = m_fileDir + "part000" + numberToString(i);
     } else if(i < 1000){
       m_pileupFile = m_fileDir + "part00" + numberToString(i);
     } else if(i < 10000){
       m_pileupFile = m_fileDir + "part0" + numberToString(i);
     } else {
       m_pileupFile = m_fileDir + "part" + numberToString(i);
     }

    m_pileupFile = PathResolver::find_file(m_pileupFile, "DATAPATH");

    return;
  }

void Atlfast::EPileupMap::emptyMap

(

)

[private]

Tidies up m_map

Definition at line 95 of file EPileupMap.cxx.

                           {
    if (m_emPileupMap->size()){
      std::vector<EPileupDeposit*>::iterator itr=m_emPileupMap->begin();
      for (; itr != m_emPileupMap->end();++itr){
        delete (*itr);
      }
    m_emPileupMap->clear();
    }
    if (m_hadPileupMap->size()){
      std::vector<EPileupDeposit*>::iterator itr=m_hadPileupMap->begin();
      for (; itr != m_hadPileupMap->end();++itr){
        delete (*itr);
      }
    m_hadPileupMap->clear();
    }
    return;
  }

---

Member Data Documentation

std::vector<EPileupDeposit*>* Atlfast::EPileupMap::m_hadPileupMap [private]

std::vector containg the EPileupDeposits

Definition at line 62 of file EPileupMap.h.

std::vector<EPileupDeposit*>* Atlfast::EPileupMap::m_emPileupMap [private]

std::vector containg the EPileupDeposits

Definition at line 64 of file EPileupMap.h.

std::map<int,double>* Atlfast::EPileupMap::m_phiMap [private]

std::map relating cell id to phi position

Definition at line 66 of file EPileupMap.h.

std::map<int,double>* Atlfast::EPileupMap::m_etaMap [private]

std::map relating cell id to eta position

Definition at line 68 of file EPileupMap.h.

std::string Atlfast::EPileupMap::m_pileupFile [private]

Name of input text file containg the EPileup values

Definition at line 70 of file EPileupMap.h.

std::string Atlfast::EPileupMap::m_fileDir [private]

Name of directory containg input files

Definition at line 72 of file EPileupMap.h.

int Atlfast::EPileupMap::m_nFiles [private]

Number of files in m_fileDir

Definition at line 74 of file EPileupMap.h.

---

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

• /afs/cern.ch/user/s/sdean/project_builds/AtlasOffline-rel_1/Simulation/Atlfast/AtlfastAlgs/AtlfastAlgs/EPileupMap.h
• /afs/cern.ch/user/s/sdean/project_builds/AtlasOffline-rel_1/Simulation/Atlfast/AtlfastAlgs/src/EPileupMap.cxx

---