totcal = 0 if totmap = 0 to avoid artifact in cal level

MartinaFeijoo · jose-luis-rs · commit 795d681c4ce9 · 2025-06-10T09:19:32.000+02:00
diff --git a/califa/calibration/R3BCalifaMapped2CrystalCal.cxx b/califa/calibration/R3BCalifaMapped2CrystalCal.cxx
@@ -259,7 +259,8 @@ void R3BCalifaMapped2CrystalCal::Exec(Option_t* /*option*/)
         {
             double a0 = params_tot.at(fNumTotParams * (crystalId - 1));
             double a1 = params_tot.at(fNumTotParams * (crystalId - 1) + 1);
-            TotCal = a0 * TMath::Exp(Tot / a1);
+            // TotCal = a0 * TMath::Exp(Tot / a1);
+            TotCal = (Tot == 0) ? 0 : a0 * TMath::Exp(Tot / a1);
         }
         AddCalData(crystalId, cal[en], cal[Nf], cal[Ns], wrts, TotCal);
     }

Original file line number	Diff line number	Diff line change
`@@ -259,7 +259,8 @@ void R3BCalifaMapped2CrystalCal::Exec(Option_t* /option/)`
`259`	`259`	`{`
`260`	`260`	`double a0 = params_tot.at(fNumTotParams * (crystalId - 1));`
`261`	`261`	`double a1 = params_tot.at(fNumTotParams * (crystalId - 1) + 1);`
`262`		`- TotCal = a0 * TMath::Exp(Tot / a1);`
	`262`	`+ // TotCal = a0 * TMath::Exp(Tot / a1);`
	`263`	`+ TotCal = (Tot == 0) ? 0 : a0 * TMath::Exp(Tot / a1);`
`263`	`264`	`}`
`264`	`265`	`AddCalData(crystalId, cal[en], cal[Nf], cal[Ns], wrts, TotCal);`
`265`	`266`	`}`