CastorTowerProducer Class Reference

#include <RecoLocalCalo/Castor/src/CastorTowerProducer.cc>

Inheritance diagram for CastorTowerProducer:

Public Member Functions
	CastorTowerProducer (const edm::ParameterSet &)
	~CastorTowerProducer () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

Private Types
typedef edm::SortedCollection< CastorRecHit >	CastorRecHitCollection
typedef edm::RefVector< CastorRecHitCollection >	CastorRecHitRefVector
typedef std::vector< reco::CastorTower >	CastorTowerCollection
typedef ROOT::Math::RhoZPhiPoint	CellPoint
typedef math::XYZPointD	Point
typedef ROOT::Math::RhoEtaPhiPoint	TowerPoint

Private Member Functions
virtual void	ComputeTowerVariable (const edm::RefVector< edm::SortedCollection< CastorRecHit > > &usedRecHits, double &Ehot, double &depth)
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
double	maxtime_
double	mintime_
edm::EDGetTokenT< CastorRecHitCollection >	tok_input_
double	towercut_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Description: CastorTower Reconstruction Producer. Produce CastorTowers from CastorCells. Implementation:

Definition at line 50 of file CastorTowerProducer.cc.

Member Typedef Documentation

typedef edm::SortedCollection<CastorRecHit> CastorTowerProducer::CastorRecHitCollection

private

Definition at line 65 of file CastorTowerProducer.cc.

typedef edm::RefVector<CastorRecHitCollection> CastorTowerProducer::CastorRecHitRefVector

private

Definition at line 67 of file CastorTowerProducer.cc.

typedef std::vector<reco::CastorTower> CastorTowerProducer::CastorTowerCollection

private

Definition at line 66 of file CastorTowerProducer.cc.

typedef ROOT::Math::RhoZPhiPoint CastorTowerProducer::CellPoint

private

Definition at line 64 of file CastorTowerProducer.cc.

typedef math::XYZPointD CastorTowerProducer::Point

private

Definition at line 62 of file CastorTowerProducer.cc.

typedef ROOT::Math::RhoEtaPhiPoint CastorTowerProducer::TowerPoint

private

Definition at line 63 of file CastorTowerProducer.cc.

Constructor & Destructor Documentation

CastorTowerProducer::CastorTowerProducer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 86 of file CastorTowerProducer.cc.

References edm::ParameterSet::getParameter(), AlCaHLTBitMon_QueryRunRegistry::string, and tok_input_.

    : towercut_(iConfig.getParameter<double>("towercut")),
      mintime_(iConfig.getParameter<double>("mintime")),
      maxtime_(iConfig.getParameter<double>("maxtime")) {
  tok_input_ = consumes<CastorRecHitCollection>(iConfig.getParameter<std::string>("inputprocess"));
  //register your products
  produces<CastorTowerCollection>();
  //now do what ever other initialization is needed
}

CastorTowerProducer::~CastorTowerProducer ( )

override

Definition at line 96 of file CastorTowerProducer.cc.

                                          {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
}

Member Function Documentation

void CastorTowerProducer::ComputeTowerVariable ( const edm::RefVector< edm::SortedCollection< CastorRecHit > > &  usedRecHits, double &  Ehot, double &  depth  )

privatevirtual

Definition at line 281 of file CastorTowerProducer.cc.

References DEFINE_FWK_MODULE, edm::Ref< C, T, F >::id(), and module::module().

Referenced by produce().

                                                              {
  using namespace reco;

  double Etot = 0;

  // loop over the cells used in the tower k
  for (CastorRecHitRefVector::iterator it = usedRecHits.begin(); it != usedRecHits.end(); it++) {
    edm::Ref<CastorRecHitCollection> rechit_p = *it;

    double Erechit = rechit_p->energy();
    HcalCastorDetId id = rechit_p->id();
    int module = id.module();
    double zrechit = 0;
    if (module < 3)
      zrechit = -14390 - 24.75 - 49.5 * (module - 1);
    if (module > 2)
      zrechit = -14390 - 99 - 49.5 - 99 * (module - 3);

    if (Erechit > Ehot)
      Ehot = Erechit;
    depth += Erechit * zrechit;
    Etot += Erechit;
  }

  depth /= Etot;
  Ehot /= Etot;
}

void CastorTowerProducer::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 106 of file CastorTowerProducer.cc.

References Abs(), ewkTauDQM_cfi::channels, ComputeTowerVariable(), LEDCalibrationChannels::depth, edm::EventSetup::get(), CastorCondObjectContainer< Item >::getAllChannels(), edm::Event::getByToken(), mps_fire::i, edm::Ref< C, T, F >::id(), dqmiolumiharvest::j, dqmdumpme::k, LogDebug, maxtime_, mintime_, module::module(), eostools::move(), AlCaHLTBitMon_ParallelJobs::p, edm::RefVector< C, T, F >::push_back(), edm::Event::put(), DetId::rawId(), edm::SortedCollection< T, SORT >::size(), findQualityFiles::size, mathSSE::sqrt(), tok_input_, and towercut_.

                                                                               {
  using namespace edm;
  using namespace reco;
  using namespace TMath;

  // Produce CastorTowers from CastorCells

  edm::Handle<CastorRecHitCollection> InputRecHits;
  iEvent.getByToken(tok_input_, InputRecHits);

  auto OutputTowers = std::make_unique<CastorTowerCollection>();

  // get and check input size
  int nRecHits = InputRecHits->size();

  LogDebug("CastorTowerProducer") << "2. entering CastorTowerProducer" << std::endl;

  if (nRecHits == 0)
    LogDebug("CastorTowerProducer") << "Warning: You are trying to run the Tower algorithm with 0 input rechits.";

  // declare castor array
  // (0,x): Energies - (1,x): emEnergies - (2,x): hadEnergies - (3,x): phi position

  double poscastortowerarray[4][16];
  double negcastortowerarray[4][16];

  CastorRecHitRefVector poscastorusedrechits[16];
  CastorRecHitRefVector negcastorusedrechits[16];

  // set phi values and everything else to zero
  for (int j = 0; j < 16; j++) {
    poscastortowerarray[3][j] = -2.94524 + j * 0.3927;
    poscastortowerarray[0][j] = 0.;
    poscastortowerarray[1][j] = 0.;
    poscastortowerarray[2][j] = 0.;

    negcastortowerarray[3][j] = -2.94524 + j * 0.3927;
    negcastortowerarray[0][j] = 0.;
    negcastortowerarray[1][j] = 0.;
    negcastortowerarray[2][j] = 0.;
  }

  // retrieve the channel quality lists from database
  edm::ESHandle<CastorChannelQuality> p;
  iSetup.get<CastorChannelQualityRcd>().get(p);
  std::vector<DetId> channels = p->getAllChannels();

  // loop over rechits to build castortowerarray[4][16] and castorusedrechits[16]
  for (unsigned int i = 0; i < InputRecHits->size(); i++) {
    edm::Ref<CastorRecHitCollection> rechit_p = edm::Ref<CastorRecHitCollection>(InputRecHits, i);

    HcalCastorDetId id = rechit_p->id();
    DetId genericID = (DetId)id;

    // first check if the rechit is in the BAD channel list
    bool bad = false;
    for (std::vector<DetId>::iterator channel = channels.begin(); channel != channels.end(); channel++) {
      if (channel->rawId() == genericID.rawId()) {
        // if the rechit is found in the list, set it bad
        bad = true;
        break;
      }
    }
    // if bad, continue the loop to the next rechit
    if (bad)
      continue;

    double Erechit = rechit_p->energy();
    int module = id.module();
    int sector = id.sector();
    double zrechit = 0;
    if (module < 3)
      zrechit = -14390 - 24.75 - 49.5 * (module - 1);
    if (module > 2)
      zrechit = -14390 - 99 - 49.5 - 99 * (module - 3);
    double phirechit = -100;
    if (sector < 9)
      phirechit = 0.19635 + (sector - 1) * 0.3927;
    if (sector > 8)
      phirechit = -2.94524 + (sector - 9) * 0.3927;

    // add time conditions for the rechit
    if (rechit_p->time() > mintime_ && rechit_p->time() < maxtime_) {
      // loop over the 16 towers possibilities
      for (int j = 0; j < 16; j++) {
        // phi matching condition
        if (TMath::Abs(phirechit - poscastortowerarray[3][j]) < 0.0001) {
          // condition over rechit z value
          if (zrechit > 0.) {
            poscastortowerarray[0][j] += Erechit;
            if (module < 3) {
              poscastortowerarray[1][j] += Erechit;
            } else {
              poscastortowerarray[2][j] += Erechit;
            }
            poscastorusedrechits[j].push_back(rechit_p);
          } else {
            negcastortowerarray[0][j] += Erechit;
            if (module < 3) {
              negcastortowerarray[1][j] += Erechit;
            } else {
              negcastortowerarray[2][j] += Erechit;
            }
            negcastorusedrechits[j].push_back(rechit_p);
          }  // end condition over rechit z value
        }    // end phi matching condition
      }      // end loop over the 16 towers possibilities
    }        // end time conditions

  }  // end loop over rechits to build castortowerarray[4][16] and castorusedrechits[16]

  // make towers of the arrays

  double fem, Ehot, depth;
  double rhoTower = 88.5;

  // loop over the 16 towers possibilities
  for (int k = 0; k < 16; k++) {
    fem = 0;
    Ehot = 0;
    depth = 0;

    // select the positive towers with E > sqrt(Nusedrechits)*Ecut
    if (poscastortowerarray[0][k] > sqrt(poscastorusedrechits[k].size()) * towercut_) {
      fem = poscastortowerarray[1][k] / poscastortowerarray[0][k];
      CastorRecHitRefVector usedRecHits = poscastorusedrechits[k];
      ComputeTowerVariable(usedRecHits, Ehot, depth);

      LogDebug("CastorTowerProducer") << "tower " << k + 1 << ": fem = " << fem << " ,depth = " << depth
                                      << " ,Ehot = " << Ehot << std::endl;

      TowerPoint temptowerposition(rhoTower, 5.9, poscastortowerarray[3][k]);
      Point towerposition(temptowerposition);

      CastorTower newtower(poscastortowerarray[0][k],
                           towerposition,
                           poscastortowerarray[1][k],
                           poscastortowerarray[2][k],
                           fem,
                           depth,
                           Ehot,
                           poscastorusedrechits[k]);
      OutputTowers->push_back(newtower);
    }  // end select the positive towers with E > Ecut

    // select the negative towers with E > sqrt(Nusedrechits)*Ecut
    if (negcastortowerarray[0][k] > sqrt(negcastorusedrechits[k].size()) * towercut_) {
      fem = negcastortowerarray[1][k] / negcastortowerarray[0][k];
      CastorRecHitRefVector usedRecHits = negcastorusedrechits[k];
      ComputeTowerVariable(usedRecHits, Ehot, depth);

      LogDebug("CastorTowerProducer") << "tower " << k + 1 << " energy = " << negcastortowerarray[0][k]
                                      << "EM = " << negcastortowerarray[1][k] << "HAD = " << negcastortowerarray[2][k]
                                      << "phi = " << negcastortowerarray[3][k] << ": fem = " << fem
                                      << " ,depth = " << depth << " ,Ehot = " << Ehot << std::endl;

      TowerPoint temptowerposition(rhoTower, -5.9, negcastortowerarray[3][k]);
      Point towerposition(temptowerposition);

      CastorTower newtower(negcastortowerarray[0][k],
                           towerposition,
                           negcastortowerarray[1][k],
                           negcastortowerarray[2][k],
                           fem,
                           depth,
                           Ehot,
                           negcastorusedrechits[k]);
      OutputTowers->push_back(newtower);
    }  // end select the negative towers with E > Ecut

  }  // end loop over the 16 towers possibilities

  iEvent.put(std::move(OutputTowers));
}

Member Data Documentation

double CastorTowerProducer::maxtime_

private

Definition at line 71 of file CastorTowerProducer.cc.

Referenced by produce().

double CastorTowerProducer::mintime_

private

Definition at line 70 of file CastorTowerProducer.cc.

Referenced by produce().

edm::EDGetTokenT<CastorRecHitCollection> CastorTowerProducer::tok_input_

private

Definition at line 68 of file CastorTowerProducer.cc.

Referenced by CastorTowerProducer(), and produce().

double CastorTowerProducer::towercut_

private

Definition at line 69 of file CastorTowerProducer.cc.

Referenced by produce().