|
|
|
#include <RecoLocalCalo/Castor/src/CastorTowerProducer.cc>
Public Member Functions | |
| CastorTowerProducer (const edm::ParameterSet &) | |
| ~CastorTowerProducer () | |
Private Types | |
| typedef std::vector < reco::CastorCell > | CastorCellCollection |
| typedef edm::RefVector < CastorCellCollection > | CastorCellRefVector |
| 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 | beginJob () |
| virtual void | ComputeTowerVariable (const reco::CastorCellRefVector &usedCells, double &Ehot, double &depth) |
| virtual void | endJob () |
| virtual void | produce (edm::Event &, const edm::EventSetup &) |
Private Attributes | |
| std::string | input_ |
| double | towercut_ |
Description: CastorTower Reconstruction Producer. Produce CastorTowers from CastorCells. Implementation:
Definition at line 44 of file CastorTowerProducer.cc.
typedef std::vector<reco::CastorCell> CastorTowerProducer::CastorCellCollection [private] |
Definition at line 59 of file CastorTowerProducer.cc.
typedef edm::RefVector<CastorCellCollection> CastorTowerProducer::CastorCellRefVector [private] |
Definition at line 61 of file CastorTowerProducer.cc.
typedef std::vector<reco::CastorTower> CastorTowerProducer::CastorTowerCollection [private] |
Definition at line 60 of file CastorTowerProducer.cc.
typedef ROOT::Math::RhoZPhiPoint CastorTowerProducer::CellPoint [private] |
Definition at line 58 of file CastorTowerProducer.cc.
typedef math::XYZPointD CastorTowerProducer::Point [private] |
Definition at line 56 of file CastorTowerProducer.cc.
typedef ROOT::Math::RhoEtaPhiPoint CastorTowerProducer::TowerPoint [private] |
Definition at line 57 of file CastorTowerProducer.cc.
| CastorTowerProducer::CastorTowerProducer | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 80 of file CastorTowerProducer.cc.
: input_(iConfig.getUntrackedParameter<std::string>("inputprocess","CastorCellReco")), towercut_(iConfig.getUntrackedParameter<double>("towercut",1.)) { //register your products produces<CastorTowerCollection>(); //now do what ever other initialization is needed }
| CastorTowerProducer::~CastorTowerProducer | ( | ) |
Definition at line 90 of file CastorTowerProducer.cc.
{
// do anything here that needs to be done at desctruction time
// (e.g. close files, deallocate resources etc.)
}
| void CastorTowerProducer::beginJob | ( | void | ) | [private, virtual] |
Reimplemented from edm::EDProducer.
Definition at line 235 of file CastorTowerProducer.cc.
References LogDebug.
{
LogDebug("CastorTowerProducer")
<<"Starting CastorTowerProducer";
}
| void CastorTowerProducer::ComputeTowerVariable | ( | const reco::CastorCellRefVector & | usedCells, |
| double & | Ehot, | ||
| double & | depth | ||
| ) | [private, virtual] |
Definition at line 246 of file CastorTowerProducer.cc.
References edm::RefVector< C, T, F >::begin(), edm::RefVector< C, T, F >::end(), and dt_offlineAnalysis_common_cff::reco.
Referenced by produce().
{
using namespace reco;
double Etot = 0;
// loop over the cells used in the tower k
for (CastorCell_iterator it = usedCells.begin(); it != usedCells.end(); it++) {
reco::CastorCellRef cell_p = *it;
double Ecell = cell_p->energy();
double zcell = cell_p->z();
if(Ecell > Ehot) Ehot = Ecell;
depth+=Ecell*zcell;
Etot+=Ecell;
}
depth/=Etot;
Ehot/=Etot;
}
| void CastorTowerProducer::endJob | ( | void | ) | [private, virtual] |
Reimplemented from edm::EDProducer.
Definition at line 241 of file CastorTowerProducer.cc.
References LogDebug.
{
LogDebug("CastorTowerProducer")
<<"Ending CastorTowerProducer";
}
| void CastorTowerProducer::produce | ( | edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [private, virtual] |
Implements edm::EDProducer.
Definition at line 102 of file CastorTowerProducer.cc.
References ComputeTowerVariable(), edm::Event::getByLabel(), i, input_, j, gen::k, LogDebug, edm::Event::put(), dt_offlineAnalysis_common_cff::reco, and towercut_.
{
using namespace edm;
using namespace reco;
using namespace TMath;
// Produce CastorTowers from CastorCells
edm::Handle<CastorCellCollection> InputCells;
iEvent.getByLabel(input_,InputCells);
std::auto_ptr<CastorTowerCollection> OutputTowers (new CastorTowerCollection);
// get and check input size
int nCells = InputCells->size();
LogDebug("CastorTowerProducer")
<<"2. entering CastorTowerProducer"<<std::endl;
if (nCells==0)
LogDebug("CastorTowerProducer") <<"Warning: You are trying to run the Tower algorithm with 0 input cells.";
// declare castor array
// (0,x): Energies - (1,x): emEnergies - (2,x): hadEnergies - (3,x): phi position
double poscastortowerarray[4][16];
double negcastortowerarray[4][16];
CastorCellRefVector poscastorusedcells[16];
CastorCellRefVector negcastorusedcells[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.;
}
// loop over cells to build castortowerarray[4][16] and castorusedcells[16]
for (unsigned int i = 0; i < InputCells->size(); i++) {
reco::CastorCellRef cell_p = reco::CastorCellRef(InputCells, i);
double Ecell = cell_p->energy();
double zcell = cell_p->z();
double phicell = cell_p->phi();
// loop over the 16 towers possibilities
for ( int j=0;j<16;j++) {
// phi matching condition
if (TMath::Abs(phicell - poscastortowerarray[3][j]) < 0.0001) {
// condition over cell z value
if (zcell > 0.) {
poscastortowerarray[0][j]+=Ecell;
if (TMath::Abs(zcell) < 14488) poscastortowerarray[1][j]+=Ecell;
else poscastortowerarray[2][j]+=Ecell;
poscastorusedcells[j].push_back(cell_p);
}
else {
negcastortowerarray[0][j]+=Ecell;
if (TMath::Abs(zcell) < 14488) negcastortowerarray[1][j]+=Ecell;
else negcastortowerarray[2][j]+=Ecell;
negcastorusedcells[j].push_back(cell_p);
} // end condition over cell z value
} // end phi matching condition
} // end loop over the 16 towers possibilities
} // end loop over cells to build castortowerarray[4][16] and castorusedcells[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 > Ecut
if (poscastortowerarray[0][k] > towercut_) {
fem = poscastortowerarray[1][k]/poscastortowerarray[0][k];
CastorCellRefVector usedCells = poscastorusedcells[k];
ComputeTowerVariable(usedCells,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,
poscastorusedcells[k]);
OutputTowers->push_back(newtower);
} // end select the positive towers with E > Ecut
// select the negative towers with E > Ecut
if (negcastortowerarray[0][k] > towercut_) {
fem = negcastortowerarray[1][k]/negcastortowerarray[0][k];
CastorCellRefVector usedCells = negcastorusedcells[k];
ComputeTowerVariable(usedCells,Ehot,depth);
LogDebug("CastorTowerProducer")
<<"tower "<<k+1<<": 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,
negcastorusedcells[k]);
OutputTowers->push_back(newtower);
} // end select the negative towers with E > Ecut
} // end loop over the 16 towers possibilities
iEvent.put(OutputTowers);
}
std::string CastorTowerProducer::input_ [private] |
Definition at line 62 of file CastorTowerProducer.cc.
Referenced by produce().
double CastorTowerProducer::towercut_ [private] |
Definition at line 63 of file CastorTowerProducer.cc.
Referenced by produce().
1.7.3