Commit 88df6798 authored by Luca Mastrolorenzo's avatar Luca Mastrolorenzo

Merge branch 'master' of gitlab.doc.ic.ac.uk:vpalladi/HgcTpgSim

parents cd0771c4 133f9edf
......@@ -19,8 +19,8 @@ class SubAnalysisPhotonTrgHit {
public:
enum {
nMCut=10,
nRCut=10
nMCut=9,
nRCut=9
};
SubAnalysisPhotonTrgHit(const std::string &sRoot="",
......@@ -140,9 +140,10 @@ public:
assert(vPhoton.size()==1);
for(unsigned j(0);j<nMCut;j++) {
double mCut(5.0+5.0*j);
double mCut(8.0+8.0*j);
for(unsigned k(0);k<nRCut;k++) {
double rCut(2.0+2.0*k);
double rCutSq(rCut*rCut);
TVectorD vd[2]={TVectorD(Coefficients::kNumberOfCoefficients),
TVectorD(Coefficients::kNumberOfCoefficients)};
......@@ -159,11 +160,18 @@ public:
for(unsigned c(0);c<Geometry::numberOfTriggerCells(l,w);c++) {
Point pCell(Geometry::point(vTrgHit[c]));
Point dp(pCell-proj.position());
if(dp.rho()<rCut) {
if(vTrgHit[c].mips()>=mCut) ml+=vTrgHit[c].mips();
if(vTrgHit[c].transverseMips()>=0.25*mCut) mlt+=vTrgHit[c].transverseMips();
}
//Point dp(pCell-proj.position());
//if(dp.rho()<rCut)
double dRho(pCell.deltaRho(proj.position()));
double dRhoPhi(pCell.rhoDeltaPhi(proj.position()));
//std::cout << "dp.rho() = " << dp.rho() << " =? " << sqrt(dRho*dRho+dRhoPhi*dRhoPhi) << std::endl;
if(3.0*3.0*dRho*dRho+dRhoPhi*dRhoPhi<rCutSq)
{
if(vTrgHit[c].mips()>=mCut) ml+=vTrgHit[c].mips();
if(vTrgHit[c].transverseMips()>=0.25*mCut) mlt+=vTrgHit[c].transverseMips();
}
}
}
}
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment