@@ -67,7 +67,7 @@ TRestSensitivity::~TRestSensitivity() {}
6767// / \param name The name of the specific metadata. It will be used to find the
6868// / corresponding TRestAxionMagneticField section inside the RML.
6969// /
70- TRestSensitivity::TRestSensitivity (const char * cfgFileName, const std::string & name)
70+ TRestSensitivity::TRestSensitivity (const char * cfgFileName, const std::string& name)
7171 : TRestMetadata(cfgFileName) {
7272 LoadConfigFromFile (fConfigFileName , name);
7373}
@@ -92,7 +92,7 @@ Double_t TRestSensitivity::ApproachByFactor(Double_t g4, Double_t chi0, Double_t
9292 Double_t Chi2 = 0 ;
9393 do {
9494 Chi2 = 0 ;
95- for (const auto & exp : fExperiments ) Chi2 += -2 * UnbinnedLogLikelihood (exp, g4);
95+ for (const auto & exp : fExperiments ) Chi2 += -2 * UnbinnedLogLikelihood (exp, g4);
9696
9797 g4 = factor * g4;
9898 } while (Chi2 - chi0 < target);
@@ -101,7 +101,7 @@ Double_t TRestSensitivity::ApproachByFactor(Double_t g4, Double_t chi0, Double_t
101101 // / Coarse movement to get to Chi2 below target (/2)
102102 do {
103103 Chi2 = 0 ;
104- for (const auto & exp : fExperiments ) Chi2 += -2 * UnbinnedLogLikelihood (exp, g4);
104+ for (const auto & exp : fExperiments ) Chi2 += -2 * UnbinnedLogLikelihood (exp, g4);
105105
106106 g4 = g4 / factor;
107107 } while (Chi2 - chi0 > target);
@@ -114,7 +114,7 @@ Double_t TRestSensitivity::ApproachByFactor(Double_t g4, Double_t chi0, Double_t
114114// /
115115Double_t TRestSensitivity::GetCoupling (Double_t sigma, Double_t precision) {
116116 Double_t Chi2_0 = 0 ;
117- for (const auto & exp : fExperiments ) Chi2_0 += -2 * UnbinnedLogLikelihood (exp, 0 );
117+ for (const auto & exp : fExperiments ) Chi2_0 += -2 * UnbinnedLogLikelihood (exp, 0 );
118118
119119 Double_t target = sigma * sigma;
120120
@@ -131,7 +131,7 @@ Double_t TRestSensitivity::GetCoupling(Double_t sigma, Double_t precision) {
131131// /////////////////////////////////////////////
132132// / \brief It returns the Log(L) for the experiment and coupling given by argument.
133133// /
134- Double_t TRestSensitivity::UnbinnedLogLikelihood (const TRestExperiment * experiment, Double_t g4) {
134+ Double_t TRestSensitivity::UnbinnedLogLikelihood (const TRestExperiment* experiment, Double_t g4) {
135135 Double_t lhood = 0 ;
136136 if (!experiment->IsDataReady ()) {
137137 RESTError << " TRestSensitivity::UnbinnedLogLikelihood. Experiment " << experiment->GetName ()
@@ -146,7 +146,7 @@ Double_t TRestSensitivity::UnbinnedLogLikelihood(const TRestExperiment *experime
146146 if (ROOT::IsImplicitMTEnabled ()) ROOT::DisableImplicitMT ();
147147
148148 std::vector<std::vector<Double_t>> trackingData;
149- for (const auto & var : experiment->GetSignal ()->GetVariables ()) {
149+ for (const auto & var : experiment->GetSignal ()->GetVariables ()) {
150150 auto values = experiment->GetExperimentalDataFrame ().Take <Double_t>(var);
151151 std::vector<Double_t> vDbl = std::move (*values);
152152 trackingData.push_back (vDbl);
@@ -169,11 +169,10 @@ Double_t TRestSensitivity::UnbinnedLogLikelihood(const TRestExperiment *experime
169169// /////////////////////////////////////////////
170170// / \brief
171171// /
172- TH1D *TRestSensitivity::SignalStatisticalTest (Int_t N) {
173-
172+ TH1D* TRestSensitivity::SignalStatisticalTest (Int_t N) {
174173 std::vector<Double_t> couplings;
175174 for (int n = 0 ; n < N; n++) {
176- for (const auto & exp : fExperiments ) exp->GetSignal ()->RegenerateActiveNodeDensity ();
175+ for (const auto & exp : fExperiments ) exp->GetSignal ()->RegenerateActiveNodeDensity ();
177176
178177 Double_t coupling = TMath::Sqrt (TMath::Sqrt (GetCoupling ()));
179178 couplings.push_back (coupling);
@@ -185,7 +184,7 @@ TH1D *TRestSensitivity::SignalStatisticalTest(Int_t N) {
185184
186185 if (fSignalTest ) delete fSignalTest ;
187186 fSignalTest = new TH1D (" SignalTest" , " A signal test" , 100 , 0.9 * min_value, 1.1 * max_value);
188- for (const auto & coup : couplings) fSignalTest ->Fill (coup);
187+ for (const auto & coup : couplings) fSignalTest ->Fill (coup);
189188
190189 return fSignalTest ;
191190}
@@ -197,18 +196,18 @@ void TRestSensitivity::InitFromConfigFile() {
197196 TRestMetadata::InitFromConfigFile ();
198197
199198 int cont = 0 ;
200- TRestMetadata * metadata = (TRestMetadata *)this ->InstantiateChildMetadata (cont, " Experiment" );
199+ TRestMetadata* metadata = (TRestMetadata*)this ->InstantiateChildMetadata (cont, " Experiment" );
201200 while (metadata != nullptr ) {
202201 cont++;
203202 if (metadata->InheritsFrom (" TRestExperimentList" )) {
204- TRestExperimentList * experimentsList = (TRestExperimentList *)metadata;
205- std::vector<TRestExperiment *> exList = experimentsList->GetExperiments ();
203+ TRestExperimentList* experimentsList = (TRestExperimentList*)metadata;
204+ std::vector<TRestExperiment*> exList = experimentsList->GetExperiments ();
206205 fExperiments .insert (fExperiments .end (), exList.begin (), exList.end ());
207206 } else if (metadata->InheritsFrom (" TRestExperiment" )) {
208- fExperiments .push_back ((TRestExperiment *)metadata);
207+ fExperiments .push_back ((TRestExperiment*)metadata);
209208 }
210209
211- metadata = (TRestMetadata *)this ->InstantiateChildMetadata (cont, " Experiment" );
210+ metadata = (TRestMetadata*)this ->InstantiateChildMetadata (cont, " Experiment" );
212211 }
213212
214213 Initialize ();
0 commit comments