| 1 |
ceballos |
1.5 |
// $Id: DiTauSystem.cc,v 1.4 2008/12/08 15:21:27 loizides Exp $
|
| 2 |
loizides |
1.1 |
|
| 3 |
|
|
#include "MitPhysics/Utils/interface/DiTauSystem.h"
|
| 4 |
loizides |
1.2 |
#include "MitAna/DataTree/interface/CompositeParticle.h"
|
| 5 |
|
|
#include "MitAna/DataTree/interface/Met.h"
|
| 6 |
ceballos |
1.5 |
#include "MitAna/DataTree/interface/Particle.h"
|
| 7 |
loizides |
1.1 |
|
| 8 |
|
|
using namespace mithep;
|
| 9 |
|
|
|
| 10 |
|
|
//--------------------------------------------------------------------------------------------------
|
| 11 |
ceballos |
1.5 |
DiTauSystem::DiTauSystem(Particle *t1, Particle *t2, const Met *met) :
|
| 12 |
loizides |
1.1 |
fT1(t1),
|
| 13 |
|
|
fT2(t2),
|
| 14 |
loizides |
1.2 |
fMet(met),
|
| 15 |
|
|
fRecoMass(0),
|
| 16 |
|
|
fVisMass(0),
|
| 17 |
|
|
fMT(0),
|
| 18 |
|
|
fETll(0),
|
| 19 |
|
|
fETnn(0)
|
| 20 |
loizides |
1.1 |
{
|
| 21 |
loizides |
1.2 |
// Constructor.
|
| 22 |
|
|
|
| 23 |
|
|
Init();
|
| 24 |
loizides |
1.1 |
}
|
| 25 |
|
|
|
| 26 |
|
|
//--------------------------------------------------------------------------------------------------
|
| 27 |
loizides |
1.2 |
void DiTauSystem::Init()
|
| 28 |
loizides |
1.1 |
{
|
| 29 |
loizides |
1.2 |
// Calculate the kinematical variables.
|
| 30 |
|
|
|
| 31 |
|
|
CompositeParticle tt;
|
| 32 |
|
|
tt.AddDaughter(fT1);
|
| 33 |
|
|
tt.AddDaughter(fT2);
|
| 34 |
|
|
|
| 35 |
|
|
CompositeParticle higgs;
|
| 36 |
|
|
higgs.AddDaughter(fT1);
|
| 37 |
|
|
higgs.AddDaughter(fT2);
|
| 38 |
|
|
higgs.AddDaughter(fMet);
|
| 39 |
|
|
|
| 40 |
|
|
Double_t xvar[3];
|
| 41 |
|
|
xvar[0] = higgs.Px()*fT2->Py()-higgs.Py()*fT2->Px();
|
| 42 |
|
|
xvar[1] = higgs.Py()*fT1->Px()-higgs.Px()*fT1->Py();
|
| 43 |
loizides |
1.1 |
xvar[2] = fT1->Px()*fT2->Py()-fT1->Py()*fT2->Px();
|
| 44 |
|
|
|
| 45 |
loizides |
1.2 |
for (Int_t i=0; i<2; ++i)
|
| 46 |
loizides |
1.1 |
xvar[i]==0 ? fXTau[i]=0 : fXTau[i]=xvar[2]/xvar[i];
|
| 47 |
|
|
|
| 48 |
loizides |
1.2 |
fVisMass = tt.Mass();
|
| 49 |
loizides |
1.1 |
if (fXTau[0] > 0 && fXTau[1] > 0)
|
| 50 |
loizides |
1.2 |
fRecoMass = fVisMass / TMath::Sqrt(fXTau[0]*fXTau[1]);
|
| 51 |
loizides |
1.1 |
else
|
| 52 |
|
|
fRecoMass = 0;
|
| 53 |
|
|
|
| 54 |
loizides |
1.2 |
Double_t visMassS = fVisMass*fVisMass;
|
| 55 |
|
|
Double_t ptll = tt.Pt();
|
| 56 |
|
|
Double_t ptmis = fMet->Pt();
|
| 57 |
|
|
if (visMassS > 0) {
|
| 58 |
|
|
fETll = TMath::Sqrt(ptll*ptll + visMassS);
|
| 59 |
|
|
fETnn = TMath::Sqrt(ptmis*ptmis + visMassS);
|
| 60 |
|
|
fMT = (fETll+fETnn)*(fETll+fETnn)-(ptll+ptmis)*(ptll+ptmis);
|
| 61 |
|
|
(fMT > 0) ? fMT=TMath::Sqrt(fMT) : fMT=0;
|
| 62 |
|
|
}
|
| 63 |
loizides |
1.1 |
}
|
