Fork me on GitHub

source: git/modules/CscClusterEfficiency.cc

Last change on this file was cc8716b, checked in by GitHub <noreply@…>, 3 years ago

Update to handle CMS endcap muon detector showers for long-lived particles (#103)

Co-authored-by: christinaw97 <christina.wang@…>

  • Property mode set to 100644
File size: 3.5 KB
Line 
1/*
2 * Delphes: a framework for fast simulation of a generic collider experiment
3 * Copyright (C) 2012-2014 Universite catholique de Louvain (UCL), Belgium
4 *
5 * This program is free software: you can redistribute it and/or modify
6 * it under the terms of the GNU General Public License as published by
7 * the Free Software Foundation, either version 3 of the License, or
8 * (at your option) any later version.
9 *
10 * This program is distributed in the hope that it will be useful,
11 * but WITHOUT ANY WARRANTY; without even the implied warranty of
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13 * GNU General Public License for more details.
14 *
15 * You should have received a copy of the GNU General Public License
16 * along with this program. If not, see <http://www.gnu.org/licenses/>.
17 */
18
19 /** \class CscClusterEfficiency
20 *
21 * This module is specific to the CMS paper searching for neutral LLPs in the CMS endcap muon detectors: https://arxiv.org/abs/2107.04838
22 * It is implemented based on the ClusterEfficiency parameterization function provided in the HEPData entry of the paper: https://www.hepdata.net/record/104408
23 *
24 * \author Christina Wang
25 *
26 */
27
28#include "modules/CscClusterEfficiency.h"
29
30#include "classes/DelphesClasses.h"
31#include "classes/DelphesFactory.h"
32#include "classes/DelphesCscClusterFormula.h"
33
34#include "ExRootAnalysis/ExRootClassifier.h"
35#include "ExRootAnalysis/ExRootFilter.h"
36#include "ExRootAnalysis/ExRootResult.h"
37
38#include "TDatabasePDG.h"
39#include "TFormula.h"
40#include "TLorentzVector.h"
41#include "TMath.h"
42#include "TObjArray.h"
43#include "TRandom3.h"
44#include "TString.h"
45
46#include <algorithm>
47#include <iostream>
48#include <sstream>
49#include <stdexcept>
50
51using namespace std;
52
53//------------------------------------------------------------------------------
54
55CscClusterEfficiency::CscClusterEfficiency() :
56 fFormula(0), fItInputArray(0)
57{
58 fFormula = new DelphesCscClusterFormula;
59}
60
61//------------------------------------------------------------------------------
62
63CscClusterEfficiency::~CscClusterEfficiency()
64{
65 if(fFormula) delete fFormula;
66}
67
68//------------------------------------------------------------------------------
69
70void CscClusterEfficiency::Init()
71{
72 // read CscClusterEfficiency formula
73 fFormula->Compile(GetString("EfficiencyFormula", "1.0"));
74
75 // import input array
76
77 fInputArray = ImportArray(GetString("InputArray", "ParticlePropagator/stableParticles"));
78 fItInputArray = fInputArray->MakeIterator();
79
80 // create output array
81
82 fOutputArray = ExportArray(GetString("OutputArray", "stableParticles"));
83}
84
85//------------------------------------------------------------------------------
86
87void CscClusterEfficiency::Finish()
88{
89 if(fItInputArray) delete fItInputArray;
90}
91
92//------------------------------------------------------------------------------
93
94void CscClusterEfficiency::Process()
95{
96 Candidate *candidate;
97 Double_t Ehad, Eem, decayR, decayZ;
98
99 fItInputArray->Reset();
100 while((candidate = static_cast<Candidate *>(fItInputArray->Next())))
101 {
102 const TLorentzVector &candidateDecayPosition = candidate->DecayPosition;
103 decayZ = abs(candidateDecayPosition.Z());
104 decayR = sqrt(pow(candidateDecayPosition.X(),2)+pow(candidateDecayPosition.Y(),2));
105 Ehad = candidate->Ehad;
106 Eem = candidate->Eem;
107 // apply an efficency formula
108 if(gRandom->Uniform() > fFormula->Eval(decayR, decayZ, Ehad, Eem)) continue;
109
110
111 fOutputArray->Add(candidate);
112 }
113}
114
115//------------------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.