// Smear
// the mu+
Float_t thetas1, phis1, ps1, thetas2, phis2, ps2, pts1, pts2, etas1, etas2;
- Float_t trigeffpL, trigeffpH, trigeffnL, trigeffnH;
+ Float_t trigeffpL, trigeffpH, trigeffnL, trigeffnH, trigeffnA;
res->SetCharge(1);
eff->SetCharge(1);
Float_t effp = eff->Evaluate(pt1, thetad1, phid1);
Float_t accp = acc->Evaluate(pt1, thetad1, phid1);
- trigeff->Evaluate(1, pt1, thetad1, phid1, trigeffpL, trigeffpH);
+ trigeff->Evaluate(1, pt1, thetad1, phid1, trigeffpL, trigeffpH, trigeffnA);
thetas1 *= TMath::Pi()/180.;
phis1 *= TMath::Pi()/180.;
Float_t etas;
res->Evaluate(p2, thetad2, phid2, ps2, thetas2, phis2);
Float_t effn = eff->Evaluate(pt2, thetad2, phid2);
Float_t accn = acc->Evaluate(pt2, thetad2, phid2);
- trigeff->Evaluate(-1, pt2, thetad2, phid2, trigeffnL, trigeffnH);
+ trigeff->Evaluate(-1, pt2, thetad2, phid2, trigeffnL, trigeffnH, trigeffnA);
thetas2 *= TMath::Pi()/180.;
phis2 *= TMath::Pi()/180.;
if (TMath::Tan(thetas2/2.) > 0) {
rl->TreeE()->Fill();
header->SetNprimary(stack->GetNprimary());
header->SetNtrack(stack->GetNtrack());
+ rl->WriteKinematics("OVERWRITE");
} // event loop
TCanvas *c1 = new TCanvas("c1","Canvas 1",400,10,600,700);
massHU->Draw();
gener->SetMomentumRange(0,999);
gener->SetPtRange(0,100.);
- gener->SetPhiRange(-180, 180);
+ gener->SetPhiRange(0., 360.);
gener->SetYRange(2.5,4);
gener->SetCutOnChild(1);
gener->SetChildThetaRange(2.0,9.);