---

testXANAGeometry.cpp

Go to the documentation of this file.

#include <XANADOO/domain/interface/XANAEsdEvent.h>
#include <XANADOO/XANAClusters/interface/XANACluster.h>
#include <XANADOO/XANAClusters/interface/XANACaloRecHit.h>
#include <XANADOO/XANAClusters/interface/XANAEmCluster.h>
#include <XANADOO/XANAClusters/interface/XANASuperCluster.h>
#include <XANADOO/XANAClusters/interface/XANAHadCluster.h>
#include <XANADOO/XANATracks/interface/XANATrack.h>
#include <XANADOO/XANAMcInfo/interface/XANAGeneratorEvent.h>
#include <XANADOO/XANAMcInfo/interface/XANAGeantEvent.h>
#include <XANADOO/XANAMcInfo/interface/XANAGeneratorParticle.h>
#include <XANADOO/XANAElectronCandidate/interface/XANAElectronCandidate.h>
#include <XANADOO/XANAElectronCandidate/interface/XANAElectronTrack.h>
#include <XANADOO/XANAElectronCandidate/interface/XANAElectronSeed.h>
#include <XANADOO/XANAMuons/interface/XANAMuonCandidate.h>

#include <XANADOO/XANAMcInfo/interface/XANAGeantTrack.h>
#include <XANADOO/XANAMcInfo/interface/XANAGeantVertex.h>

#include <XANADOO/XANAGeometry/interface/XANAEcalGeometry.h>

#include <string>
#include <vector>

#include <TH1.h>
#include <TH2.h>

using namespace std;

int main(int argc, char** argv)
{
  std::cout<<"[OVAL]essai dif"<<endl;
  TFile *f2 = new TFile("testXANAGeometry.root", "recreate");
  f2->cd();

  TH1F * hIetaORx = new TH1F("hIetaORx","IetaORx",170, 0., 170.);
  TH2F * hIetaORx2 = new TH2F("hIetaORx2","IetaORx",170, 0., 170., 1, 0., 170.);

  TH1F * hIphiORy = new TH1F("hIphiORy","IphiORy",360, 0., 360.);
  TH2F * hIphiORy2 = new TH2F("hIphiORy2","IphiORy",360, 0., 360., 1, 0., 360.);

  TH2F * testEnd=new TH2F("testEnd","ENDCAP",100, 0., 100., 100, 0., 100.);
  TH2F * mapEndcap= new TH2F("mapEndcap", "map", 100, 0., 100., 100, 0., 100);
  Float_t position[3];
  XANAEcalGeometry * geom = new XANAEcalGeometry();
  Int_t ordre=1;

  for(int xtal=0; xtal<10000; xtal++)
    {
      if (geom->isXtal(xtal,-1)==0)
        {
          Int_t Ieta=xtal%100;
          Int_t Iphi=int(xtal/100);
          mapEndcap->Fill(Ieta, Iphi, 100.);
        }
      
    }

  position[0]=-1.55589;
  position[1]=33.8957;
  position[2]=-316.538;
  //This is the position of a XANACellID in the Endcap
  

  HepPoint3D pos(position[0],position[1],position[2]);
  
  XANACellID cell = XANACellID(pos);
  Int_t Ieta=geom->getEtaORxIndex(cell);
  Int_t Iphi=geom->getPhiORyIndex(cell);
  
  testEnd->Fill(Ieta, Iphi, ordre);
  ordre++;
    
  XANACellID east=geom->getEast(cell);
  
  while(!(east==cell))
    {
      cell=east;
      east=geom->getEast(cell);
    }   
  
  Ieta=geom->getEtaORxIndex(cell);
  hIetaORx->Fill(Ieta);
  hIetaORx2->Fill(Ieta,1);

  Iphi=geom->getPhiORyIndex(cell);
  hIphiORy->Fill(Iphi);
  hIphiORy2->Fill(Iphi,1);
  
  testEnd->Fill(Ieta, Iphi, ordre);
  ordre++;
  
  XANACellID north=geom->getNorth(cell);
  
  while(!(north==cell))
    {
      cell=north;
    }
  
  Iphi=geom->getPhiORyIndex(cell);
  Ieta=geom->getEtaORxIndex(cell);
  testEnd->Fill(Ieta, Iphi, ordre);
  ordre++;
  
  cell=geom->getNorth(cell);
  Iphi=geom->getPhiORyIndex(cell);
  Ieta=geom->getEtaORxIndex(cell);
  testEnd->Fill(Ieta, Iphi, ordre);
  ordre++;
  
  XANACellID west=geom->getWest(cell);
  
  while(!(west==cell))
    {
      cell=west;
      west=geom->getWest(cell);
    }   
  Iphi=geom->getPhiORyIndex(cell);
  Ieta=geom->getEtaORxIndex(cell);
  testEnd->Fill(Ieta, Iphi, ordre);
  ordre++;
  
  cell=geom->getSouth(cell);
  Iphi=geom->getPhiORyIndex(cell);
  Ieta=geom->getEtaORxIndex(cell);
  testEnd->Fill(Ieta, Iphi, ordre);
  ordre++; std::cout<<"[OVAL]  Ieta="<<Ieta<<"  Iphi="<<Iphi<<std::endl;
  
  
  
  ordre=1;
  cell = XANACellID(pos);
  
  for(int dec=1; dec<25; dec++)
    {
      cell=geom->getWest(cell);
    }
  
  Iphi=geom->getPhiORyIndex(cell);
  Ieta=geom->getEtaORxIndex(cell);
  testEnd->Fill(Ieta, Iphi, ordre);
  ordre++; std::cout<<"[OVAL]  Ieta="<<Ieta<<"  Iphi="<<Iphi<<std::endl;
  
  cell=geom->getNorthWest(cell);
  Iphi=geom->getPhiORyIndex(cell);
  Ieta=geom->getEtaORxIndex(cell);
  testEnd->Fill(Ieta, Iphi, ordre);
  ordre++; std::cout<<"[OVAL]  Ieta="<<Ieta<<"  Iphi="<<Iphi<<std::endl;
  cell=geom->getWest(cell);
  Iphi=geom->getPhiORyIndex(cell);
  Ieta=geom->getEtaORxIndex(cell);
  testEnd->Fill(Ieta, Iphi, ordre);
  ordre++; std::cout<<"[OVAL]  Ieta="<<Ieta<<"  Iphi="<<Iphi<<std::endl;
  cell=geom->getSouthWest(cell);
  Iphi=geom->getPhiORyIndex(cell);
  Ieta=geom->getEtaORxIndex(cell);
  testEnd->Fill(Ieta, Iphi, ordre);
  ordre++; std::cout<<"[OVAL]  Ieta="<<Ieta<<"  Iphi="<<Iphi<<std::endl;
  cell=geom->getSouth(cell);
  Iphi=geom->getPhiORyIndex(cell);
  Ieta=geom->getEtaORxIndex(cell);
  testEnd->Fill(Ieta, Iphi, ordre);
  ordre++; std::cout<<"[OVAL]  Ieta="<<Ieta<<"  Iphi="<<Iphi<<std::endl;
  cell=geom->getSouthEast(cell);
  Iphi=geom->getPhiORyIndex(cell);
  Ieta=geom->getEtaORxIndex(cell);
  testEnd->Fill(Ieta, Iphi, ordre);
  ordre++; std::cout<<"[OVAL]  Ieta="<<Ieta<<"  Iphi="<<Iphi<<std::endl;
  cell=geom->getEast(cell);
  Iphi=geom->getPhiORyIndex(cell);
  Ieta=geom->getEtaORxIndex(cell);
  testEnd->Fill(Ieta, Iphi, ordre);
  ordre++; std::cout<<"[OVAL]  Ieta="<<Ieta<<"  Iphi="<<Iphi<<std::endl;
  cell=geom->getSouthEast(cell);
  Iphi=geom->getPhiORyIndex(cell);
  Ieta=geom->getEtaORxIndex(cell);
  testEnd->Fill(Ieta, Iphi, ordre);
  ordre++; std::cout<<"[OVAL]  Ieta="<<Ieta<<"  Iphi="<<Iphi<<std::endl;
  cell=geom->getSouth(cell);
  Iphi=geom->getPhiORyIndex(cell);
  Ieta=geom->getEtaORxIndex(cell);
  testEnd->Fill(Ieta, Iphi, ordre);
  ordre++; std::cout<<"[OVAL]  Ieta="<<Ieta<<"  Iphi="<<Iphi<<std::endl;
  cell=geom->getSouth(cell);
  Iphi=geom->getPhiORyIndex(cell);
  Ieta=geom->getEtaORxIndex(cell);
  testEnd->Fill(Ieta, Iphi, ordre);
  ordre++; std::cout<<"[OVAL]  Ieta="<<Ieta<<"  Iphi="<<Iphi<<std::endl;
  cell=geom->getSouthWest(cell);
  Iphi=geom->getPhiORyIndex(cell);
  Ieta=geom->getEtaORxIndex(cell);
  testEnd->Fill(Ieta, Iphi, ordre);
  ordre++; std::cout<<"[OVAL]  Ieta="<<Ieta<<"  Iphi="<<Iphi<<std::endl;
  cell=geom->getWest(cell);
  Iphi=geom->getPhiORyIndex(cell);
  Ieta=geom->getEtaORxIndex(cell);
  testEnd->Fill(Ieta, Iphi, ordre);
  ordre++; std::cout<<"[OVAL]  Ieta="<<Ieta<<"  Iphi="<<Iphi<<std::endl;
  cell=geom->getNorthWest(cell);
  Iphi=geom->getPhiORyIndex(cell);
  Ieta=geom->getEtaORxIndex(cell);
  testEnd->Fill(Ieta, Iphi, ordre);
  ordre++; std::cout<<"[OVAL]  Ieta="<<Ieta<<"  Iphi="<<Iphi<<std::endl;
  cell=geom->getNorth(cell);
  Iphi=geom->getPhiORyIndex(cell);
  Ieta=geom->getEtaORxIndex(cell);
  testEnd->Fill(Ieta, Iphi, ordre);
  ordre++; std::cout<<"[OVAL]  Ieta="<<Ieta<<"  Iphi="<<Iphi<<std::endl;
  
  
  f2->Write();
  f2->Close();
  
}

---