---

EtMissTrivialTest.cpp

Go to the documentation of this file.

#include <Utilities/Configuration/interface/Architecture.h>

#include "CARF/G3Event/interface/G3EventProxy.h"
#include "Utilities/Notification/interface/Observer.h"

// Calorimetry stuff
#include "Calorimetry/CaloRecHit/interface/CaloRecHit.h"
#include "Calorimetry/EcalPlusHcalTower/interface/EcalPlusHcalTower.h"
#include "Calorimetry/EcalPlusHcalTower/src/EcalPlusHcalTowerBuilder.h"
#include "Calorimetry/CaloCommon/interface/CaloItr.h"
#include "Calorimetry/CaloDetector/interface/CellID.h"

// reconstructor stuff
#include "CARF/Reco/interface/RecItr.h"
#include "CARF/Reco/interface/Reconstructor.h"

#include "Utilities/Notification/interface/PackageInitializer.h"

#include "CLHEP/Geometry/Point3D.h"
#include "CLHEP/Units/PhysicalConstants.h"

// EtMiss class
#include <XANADOO/XANAEtMiss/interface/XANAEtMiss.h>

#include <iostream>

class EtMissTrivialTest : private Observer<G3EventProxy *>
{
public:
  EtMissTrivialTest(){    Observer<G3EventProxy *>::init() ;}
  ~EtMissTrivialTest(){ cout<<"End of EtMissTrivialTest"<<endl;}
  void FillingTest(G3EventProxy *ev);  
private:
  void upDate(G3EventProxy *ev){if (ev!=0) FillingTest(ev);}
  void FillEtMiss(G3EventProxy *ev, EtMiss *a);
  void Print(EtMiss a);
};

void EtMissTrivialTest::FillingTest(G3EventProxy *ev)
{
  EtMiss* myMet = new EtMiss();
  FillEtMiss (ev,myMet);
  Print(*myMet);
  delete myMet;
}

void EtMissTrivialTest::Print(EtMiss a)
{
  cout<<"\n ********* ECAL ******************"<<endl;
  for (int i=0;i<5;i++){cout<<a.getEECal(i)<<" ";}
  cout<<endl;
  for (int i=0;i<5;i++){cout<<a.getEScalarECal(i)<<" ";}

  cout<<"\n ********* HCAL ******************"<<endl;
  for (int i=0;i<5;i++){cout<<a.getEHCal(i)<<" ";}
  cout<<endl;
  for (int i=0;i<5;i++){cout<<a.getEScalarHCal(i)<<" ";}

  cout<<"\n ********* Total ******************"<<endl;
  for (int i=0;i<5;i++){cout<<a.getE(i)<<" ";}
  cout<<endl;
  for (int i=0;i<5;i++){cout<<a.getEScalar(i)<<" ";}

  cout<<"\n ********* MC ******************"<<endl;
  for (int i=0;i<5;i++){cout<<a.getETrue(i)<<" ";}
  cout<<endl<<endl;

}

  void EtMissTrivialTest::FillEtMiss(G3EventProxy *ev, EtMiss* a)
{
  a->SetEECal(1,2.5,4.5);
  a->SetEHCal(10,20.5,40.5);

  a->SetEScalECal(11,12,13,55,56);
  a->SetEScalECal(110,120,130,70,80);
  a->SetETotal(50,60,70);
  a->SetEScalTotal(0.9,0.4,0.3,10.3,20.3);

}

static PKBuilder<EtMissTrivialTest> pippo("EtMissTrivialTest") ;

---