---

ElectronEcConfigurer.cxx

Go to the documentation of this file.

#include "FastShowerUtils/ElectronEcConfigurer.h"


#include "FastShowerUtils/Samplers/S0EmEcalEc1.h"
#include "FastShowerUtils/Samplers/PeakMeanS0EmEcalEc1.h"
#include "FastShowerUtils/Samplers/PeakSigmaS0EmEcalEc1.h"
#include "FastShowerUtils/Samplers/TailFracS0EmEcalEc1.h"
#include "FastShowerUtils/Samplers/TailAuxS0EmEcalEc1.h"
//Cell0
#include "FastShowerUtils/Samplers/C0EmEcalEc1.h"
#include "FastShowerUtils/Samplers/PeakMeanC0EmEcalEc1.h"
#include "FastShowerUtils/Samplers/PeakSigmaC0EmEcalEc1.h"
#include "FastShowerUtils/Samplers/TailFracC0EmEcalEc1.h"
#include "FastShowerUtils/Samplers/TailAuxC0EmEcalEc1.h"
//SN1
//#include "FastShowerUtils/SamplerRouterC0Sample.h"
#include "FastShowerUtils/SNSamplerRouterC0Sample.h"
//
#include "FastShowerUtils/Samplers/SNEmEcalEc1.h"
#include "FastShowerUtils/Samplers/PeakMeanSNEmEcalEc1.h"
#include "FastShowerUtils/Samplers/PeakSigmaSNEmEcalEc1.h"
#include "FastShowerUtils/Samplers/TailFracSNEmEcalEc1.h"
//SN2
#include "FastShowerUtils/Samplers/SNEmEcalEc2.h"
#include "FastShowerUtils/Samplers/TailFracSNEmEcalEc2.h"
//SN3 - as SN1
#include "FastShowerUtils/Samplers/SNEmEcalEc3.h"
//AN1
//#include "FastShowerUtils/SamplerRouterSNSample.h"
#include "FastShowerUtils/ANSamplerRouterSNSample.h"
#
#include "FastShowerUtils/Samplers/ANEmEcalEc1.h"
#include "FastShowerUtils/Samplers/PeakMeanANEmEcalEc1.h"
#include "FastShowerUtils/Samplers/PeakSigmaANEmEcalEc1.h"
//AN2
#include "FastShowerUtils/Samplers/ANEmEcalEc2.h"
#include "FastShowerUtils/Samplers/PeakMeanANEmEcalEc2.h"
#include "FastShowerUtils/Samplers/PeakSigmaANEmEcalEc2.h"
//
//Normalisers
#include "FastShowerUtils/Normalisers/EmEtaResponse.h"
#include "FastShowerUtils/Normalisers/EmHcalLeakageEc.h"
//
#include "FastShowerUtils/PtrConverter.h"
#include "FastShowerUtils/ISampler.h"
#include "FastShowerUtils/SamplerPtr.h"
#include "FastShowerUtils/INormaliser.h"
#include "FastShowerUtils/NormaliserPtr.h"
#include "FastShowerUtils/IDepositor.h"
#include "FastShowerUtils/DepositorPtr.h"
#include "FastShowerUtils/CoreDepositor.h"
#include "FastShowerUtils/HaloDepositor.h"
#include "FastShowerUtils/PencilDepositor.h"
//
#include "FastShowerUtils/PtrAdapter.h"
//
#include "FastShowerUtils/Normalisers/EmEcalNormaliser.h"
#include "FastShowerUtils/Normalisers/EmHcalNormaliser.h"
//
#include <string>
#include <algorithm>

namespace FastShower{

  //
  ElectronEcConfigurer::ElectronEcConfigurer(Moni m):ConfigurerBase(m){
    this->initialise();
  }
  ElectronEcConfigurer::ElectronEcConfigurer(NoMoni nm):ConfigurerBase(nm){
    this->initialise();
  }
  void ElectronEcConfigurer::initialise(){
    //Slice0
    this->addFn("S0PeakMean",        new  PeakMeanS0EmEcalEc1);
    this->addFn("S0PeakSigma",       new  PeakSigmaS0EmEcalEc1);
    this->addFn("S0TailFrac",        new  TailFracS0EmEcalEc1);
    this->addFn("S0TailLowerBound",  new  TailAuxS0EmEcalEc1);
    //Cell0
    this->addFn("C0PeakMean",        new  PeakMeanC0EmEcalEc1);
    this->addFn("C0PeakSigma",       new  PeakSigmaC0EmEcalEc1);
    this->addFn("C0TailFrac",        new  TailFracC0EmEcalEc1);
    this->addFn("C0TailLowerBound",  new  TailAuxC0EmEcalEc1);
    //SN1
    this->addFn("SN1PeakMean",       new  PeakMeanSNEmEcalEc1);
    this->addFn("SN1PeakSigma",      new  PeakSigmaSNEmEcalEc1);
    this->addFn("SN1TailFrac",       new  TailFracSNEmEcalEc1);
    //SN2 - uses SN1 functions for mean and sigma
    this->addFn("SN2PeakMean",       new  PeakMeanSNEmEcalEc1);
    this->addFn("SN2PeakSigma",      new  PeakSigmaSNEmEcalEc1);
    this->addFn("SN2TailFrac",       new  TailFracSNEmEcalEc2);
    //SN3 - uses SN1 functions for mean and sigma
    this->addFn("SN3PeakMean",       new  PeakMeanSNEmEcalEc1);
    this->addFn("SN3PeakSigma",      new  PeakSigmaSNEmEcalEc1);
    //AN1 
    this->addFn("AN1PeakMean",       new  PeakMeanANEmEcalEc1);
    this->addFn("AN1PeakSigma",      new  PeakSigmaANEmEcalEc1);
    //AN2 
    this->addFn("AN2PeakMean",       new  PeakMeanANEmEcalEc2);
    this->addFn("AN2PeakSigma",      new  PeakSigmaANEmEcalEc2);
    //
    //*******************************************************
    //*              Functions for Normalisers              *
    //*******************************************************
    this->addFnPP("EtaResponse",    new  EmEtaResponse);
    //*******************************************************
    //*              Functions2 for Normalisers             *
    //*******************************************************
    this->addFnPP2("E/HcalFractions", new  EmHcalLeakageEc);
    //
    //*******************************************************
    //*       IUpdating Gaussians (Peaks for Samplers)      *
    //*******************************************************
    this->addIUG("AN1Peak", this->makeUpdatingGaussian0LP("AN1") );
    this->addIUG("AN2Peak", this->makeUpdatingGaussian0LP("AN2") );
    this->addIUG("C0Peak",  this->makeUpdatingGaussian0LP("C0")  );
    this->addIUG("SN1Peak", this->makeUpdatingGaussian0LP("SN1") );
    this->addIUG("SN2Peak", this->makeUpdatingGaussian0LP("SN2") );
    this->addIUG("SN3Peak", this->makeUpdatingGaussian0LP("SN3") );
    this->addIUG("S0Peak",  this->makeUpdatingGaussian0LP("S0")  );
    //
    //*******************************************************
    //*       Samplers                                      *
    //*******************************************************
    this->addSampler("S0",  new S0EmEcalEc1(this, "")  );
    this->addSampler("C0",  new C0EmEcalEc1(this, "")  );
    this->addSampler("SN1", new SNEmEcalEc1(this, "") );
    this->addSampler("SN2", new SNEmEcalEc2(this, "") );
    this->addSampler("SN3", new SNEmEcalEc3(this, "") );
    this->addSampler("AN1", new ANEmEcalEc1(this, "") );
    this->addSampler("AN2", new ANEmEcalEc2(this, "") );
    this->addSampler("SN",  this->makeSN() );
    this->addSampler("AN",  this->makeAN() );
    //
    //*******************************************************
    //*       Depositors                                    *
    //*******************************************************
    this->addDepositorsE();
    this->addDepositorsH();
    //*******************************************************
    //*       Normalisers                                   *
    //*******************************************************
    this->normalisers();
   }
  //
  ISampler* ElectronEcConfigurer::makeSN() const {
    
    std::vector< std::pair<double, ISampler*> > v;
    v.push_back(std::pair<double, ISampler*>(0.50, this->findSampler("SN1"))); 
    v.push_back(std::pair<double, ISampler*>(0.95, this->findSampler("SN2"))); 
    v.push_back(std::pair<double, ISampler*>(1.00, this->findSampler("SN3"))); 
    //return new SamplerRouterC0Sample(v, 0.0, "SNRouter");
    return new SNSamplerRouterC0Sample(v, 0.0, "SNRouter");
    
  }
  ISampler* ElectronEcConfigurer::makeAN() const {

    std::vector< std::pair<double, ISampler*> > v;
    v.push_back(std::pair<double, ISampler*>(0.05, this->findSampler("AN1"))); 
    v.push_back(std::pair<double, ISampler*>(1.00, this->findSampler("AN2"))); 
    //return new SamplerRouterSNSample(v, 0.0, "ANRouter");
    return new ANSamplerRouterSNSample(v, 0.0, "ANRouter");

  }
  void ElectronEcConfigurer::addDepositorsE(){
    
    this->addEdepositor( "EcalCore",
                         new CoreDepositor(
                                           this->findSampler("S0"),
                                           this->findSampler("C0"),
                                           this->findSampler("SN"),
                                           this->findSampler("AN")
                                           )
                         );

  }
  void ElectronEcConfigurer::addDepositorsH(){

    this->addHdepositor( "HCalPencil", new PencilDepositor);

  }
  void ElectronEcConfigurer::normalisers(){
    this->addNormaliser("EcalNormaliser",   
                        new EmEcalNormaliser(this,"ElectronEcEcalNormaliser"));
    this->addNormaliser("HcalNormaliser",   
                        new EmHcalNormaliser(this,"ElectronEcHcalNormaliser"));
  }
}//namespace

---