#include "FCS_Cell.h"
#include "../ISF_FastCaloSimParametrization/MeanAndRMS.h"
#include "Identifier/Identifier.h"
#include "CaloDetDescr/CaloDetDescrElement.h"
#include "CaloGeometryFromFile.h"
#include "CaloHitAna.h"
#include "../ISF_FastCaloSimParametrization/CaloGeometry.h"
#include "../ISF_FastCaloSimParametrization/FCAL_ChannelMap.h"
#include "TTree.h"
#include <string>
#include <sstream>
#include <iostream>
#include "TSystem.h"
#include "TString.h"
#include "TFile.h"
#include <stdlib.h>
#include "TLorentzVector.h"
#include "TH1.h"
#include "TH2.h"
#include "TH1F.h"
#include "TH2F.h"
#include <vector>
#include "TCanvas.h"
Go to the source code of this file.
◆ TestFCalIdentifiers()
void TestFCalIdentifiers |
( |
TString |
sampling = "Sampling_0" | ) |
|
Definition at line 31 of file TestFCalIdentifiers.cxx.
34 geo->LoadGeometryFromFile(
"/afs/cern.ch/atlas/groups/Simulation/FastCaloSim/ATLAS-GEO-20-00-01.root",
"ATLAS-GEO-20-00-01");
35 geo->LoadFCalGeometryFromFiles(
"FCal1-electrodes.sorted.HV.09Nov2007.dat",
"FCal2-electrodes.sorted.HV.April2011.dat",
"FCal3-electrodes.sorted.HV.09Nov2007.dat");
55 TFile *
inputFile = TFile::Open(
"/eos/atlas/user/s/schaarsc/FCS/user.fladias.428137.FastCalo_pid11_E65536_etam35_35_zv_m100.e4001_s2864_r7736.w0_162706_matched_output.root/user.fladias.8834798._000001.matched_output.root");
64 std::cout <<
"Sampling is " << sampling << std::endl;
67 std::cout <<
"nEvt " << nEvt << std::endl;
71 stringstream ssIdentifier;
73 for (Int_t ientry=0; ientry<nEvt; ientry++){
77 std::cout <<
"Processing event # " << ientry << std::endl;
82 for (UInt_t j=0; j<(*vec).size(); j++){
89 cout << cell_ID << hex <<
" 0x"<< cell_ID << endl;
105 TLorentzVector *hitVec =
new TLorentzVector(
x,
y,
z,
t);
111 if(hitSampling<21)foundCell=
geo->getDDE(hitSampling,hitVec->Eta(),hitVec->Phi());
112 else if(hitSampling<24)foundCell=
geo->getFCalDDE(hitSampling,
x,
y,
z);
114 cout << endl <<
"Warning: Found hit with sampling > 23 !!!!!!!!!!!!!!" << endl << endl;
126 channelMap->
getTileID(hitSampling - 20,
x,
y,ieta,iphi);
127 ssIdentifier.str(
"");
134 << (ieta << 16)+iphi;
138 if(
z > 0) cout <<
"Side A " << hex << foundCell->
identify() <<
" " << cell_ID <<
" " << ((
FCS_matchedcell)(*
vec)[j]).hit[ihit].cell_identifier << dec <<
" " << cell_ID <<
" " << ((
FCS_matchedcell)(*vec)[j]).hit[ihit].cell_identifier <<
" " << ieta <<
" " << iphi <<
" " << (ieta << 16)+iphi << endl;
139 else cout <<
"Side C " << hex << foundCell->
identify() <<
" " << cell_ID <<
" " << ((
FCS_matchedcell)(*
vec)[j]).hit[ihit].cell_identifier << dec<<
" " << cell_ID <<
" " << ((
FCS_matchedcell)(*vec)[j]).hit[ihit].cell_identifier <<
" " << ieta <<
" " << iphi <<
" " << (ieta << 16)+iphi << endl;
144 else cout <<
"Cell not found!" << endl;