CaloGeometry.cxx File Reference

#include "ISF_FastCaloSimParametrization/CaloGeometry.h"
#include <TTree.h>
#include <TVector2.h>
#include <TRandom.h>
#include <TCanvas.h>
#include <TH2D.h>
#include <TGraphErrors.h>
#include <TVector3.h>
#include <TLegend.h>
#include <fstream>
#include <sstream>
#include "CaloDetDescr/CaloDetDescrElement.h"
#include "CaloGeoHelpers/CaloSampling.h"
#include "ISF_FastCaloSimEvent/FastCaloSim_CaloCell_ID.h"

Go to the source code of this file.

Functions
void CaloGeometry::Validate	ATLAS_NOT_THREAD_SAFE (int nrnd)

Variables
const int	CaloGeometry_calocol [24]

Function Documentation

ATLAS_NOT_THREAD_SAFE()

void CaloGeometry::Validate ATLAS_NOT_THREAD_SAFE

(

int

nrnd

)

Definition at line 655 of file CaloGeometry.cxx.

{
  int ntest=0;
  cout<<"start CaloGeometry::Validate()"<<endl;
  for(t_cellmap::iterator ic=m_cells.begin();ic!=m_cells.end();++ic) {
    const CaloDetDescrElement* cell=ic->second;
    int sampling=cell->getSampling();
    //if(sampling>=21) continue;
 
    if(m_debug_identify==cell->identify()) {
      cout<<"CaloGeometry::Validate(), cell "<<ntest<<" id="<<cell->identify()<<endl; 
      m_debug=true;
    }  
    
    for(int irnd=0;irnd<nrnd;++irnd) {
      std::stringstream pos;
      std::stringstream cellpos;
      std::stringstream foundcellpos;
      int steps=0;
      float distance=0;
      bool is_inside;
      bool is_inside_foundcell=false;
      const CaloDetDescrElement* foundcell=nullptr;
      if(sampling<21) {
        float eta=cell->eta()+1.95*(gRandom->Rndm()-0.5)*cell->deta();
        float phi=cell->phi()+1.95*(gRandom->Rndm()-0.5)*cell->dphi();
        foundcell=getDDE(sampling,eta,phi,&distance,&steps);
        
        pos<<"eta="<<eta<<" phi="<<phi;
        cellpos<<"eta="<<cell->eta()<<" eta_raw="<<cell->eta_raw()<<" deta="<<cell->deta()
               <<" ("<<(cell->eta_raw()-cell->eta())/cell->deta()<<") ; "
               <<"phi="<<cell->phi()<<" phi_raw="<<cell->phi_raw()<<" dphi="<<cell->dphi()
               <<" ("<<(cell->phi_raw()-cell->phi())/cell->dphi()<<")";
        if(foundcell) {
          foundcellpos<<"eta="<<foundcell->eta()<<" eta_raw="<<foundcell->eta_raw()<<" deta="<<foundcell->deta()
                      <<" ("<<(foundcell->eta_raw()-foundcell->eta())/foundcell->deta()<<") ; "
                      <<"phi="<<foundcell->phi()<<" phi_raw="<<foundcell->phi_raw()<<" dphi="<<foundcell->dphi()
                      <<" ("<<(foundcell->phi_raw()-foundcell->phi())/cell->dphi()<<")";
          is_inside_foundcell=TMath::Abs( (eta-foundcell->eta())/foundcell->deta() )<0.55 && TMath::Abs( (phi-foundcell->phi())/foundcell->dphi() )<0.55;
        }
        is_inside=TMath::Abs( (eta-cell->eta())/cell->deta() )<0.49 && TMath::Abs( (phi-cell->phi())/cell->dphi() )<0.49;
      } else {
        float x=cell->x()+1.95*(gRandom->Rndm()-0.5)*cell->dx();
        float y=cell->y()+1.95*(gRandom->Rndm()-0.5)*cell->dy();
        float z=cell->z();
        foundcell=getFCalDDE(sampling,x,y,z);
        
        pos<<"x="<<x<<" y="<<y<<" z="<<z;
        cellpos<<"x="<<cell->x()<<" x_raw="<<cell->x_raw()<<" dx="<<cell->dx()
               <<" ("<<(cell->x_raw()-cell->x())/cell->dx()<<") ; "
               <<"y="<<cell->y()<<" y_raw="<<cell->y_raw()<<" dy="<<cell->dy()
               <<" ("<<(cell->y_raw()-cell->y())/cell->dy()<<")";
        if(foundcell) {
          foundcellpos<<"x="<<foundcell->x()<<" x_raw="<<foundcell->x_raw()<<" dx="<<foundcell->dx()
                      <<" ("<<(foundcell->x_raw()-foundcell->x())/foundcell->dx()<<") ; "
                      <<"y="<<foundcell->y()<<" y_raw="<<foundcell->y_raw()<<" dy="<<foundcell->dy()
                      <<" ("<<(foundcell->y_raw()-foundcell->y())/cell->dy()<<")";
          is_inside_foundcell=TMath::Abs( (x-foundcell->x())/foundcell->dx() )<0.75 && TMath::Abs( (y-foundcell->y())/foundcell->dy() )<0.75;
        }       
        is_inside=TMath::Abs( (x-cell->x())/cell->dx() )<0.49 && TMath::Abs( (y-cell->y())/cell->dy() )<0.49;
        //m_debug=true;
      }  
 
      if(m_debug && foundcell) {
        cout<<"CaloGeometry::Validate(), irnd="<<irnd<<": cell id="<<cell->identify()<<", sampling="<<sampling
                                                     <<", foundcell id="<<foundcell->identify()<<", "<<steps<<" steps"<<endl; 
        cout<<"  "<<cellpos.str()<<endl;
      }  
      if(cell==foundcell) {
        if(steps>3 && distance<-0.01) { 
          cout<<"cell id="<<cell->identify()<<", sampling="<<sampling<<" found in "<<steps<<" steps, dist="<<distance<<" "<<pos.str()<<endl;
          cout<<"  "<<cellpos.str()<<endl;
        }
      } else {
        if(!foundcell) {
          cout<<"cell id="<<cell->identify()<<" not found!";
          cout<<" No cell found in "<<steps<<" steps, dist="<<distance<<" "<<pos.str()<<endl;
          cout<<"  input sampling="<<sampling<<" "<<cellpos.str()<<endl;
          //return;
        }
        if( is_inside && foundcell && !is_inside_foundcell) {
          cout<<"cell id="<<cell->identify()<<" not found, but inside cell area!";
          cout<<" Found instead id="<<foundcell->identify()<<" in "<<steps<<" steps, dist="<<distance<<" "<<pos.str()<<endl;
          cout<<"  input sampling="<<sampling<<" "<<cellpos.str()<<endl;
          cout<<"  output sampling="<<foundcell->getSampling()<<" "<<foundcellpos.str()<<endl;
          return;
        }  
      }
    }  
    m_debug=false;
    if(ntest%25000==0) cout<<"Validate cell "<<ntest<<" with "<<nrnd<<" random hits"<<endl;
    ++ntest;
  }
  cout<<"end CaloGeometry::Validate()"<<endl;
}

Variable Documentation

CaloGeometry_calocol

const int CaloGeometry_calocol[24]

Initial value:

={1,2,3,4, 
                                    1,2,3,4, 
                                    1,2,3,4, 
                                    1,2,3,   
                                    -42,-28,6, 
                                    1,2,3,   
                                    1,2,3    
}

Definition at line 27 of file CaloGeometry.cxx.

                                   {1,2,3,4, // LAr barrel
                                    1,2,3,4, // LAr EM endcap
                                    1,2,3,4, // Hadronic end cap cal.
                                    1,2,3,   // Tile barrel
                                    -42,-28,6, // Tile gap (ITC & scint)
                                    1,2,3,   // Tile extended barrel
                                    1,2,3    // Forward EM endcap
};