diff --git a/macros/CCEScan/landauFit.C b/macros/CCEScan/landauFit.C
index 6a4a17e..8090a7a 100755
--- a/macros/CCEScan/landauFit.C
+++ b/macros/CCEScan/landauFit.C
@@ -116,6 +116,8 @@
   bool notSave  = false;
   bool temSave  = false;
   double radius = -10.0;
+  double radmin = 0.0;
+  double radmax = 2000.0;
   
   for (int i=1; i<argc; i++) {
     if (strcmp(argv[i],"-d")==0 && i+1<argc) {
@@ -184,9 +186,22 @@
   // Set the radius to an irrelevant value
   if (radius<0.0) {
     radius = 2000.0;
-  }else{
+  } else {
     vn_result  +=Form("_r%d",(int)radius);
     fn_graph_s +=Form("_r%d",(int)radius);
+
+    // Elena - implement circular crowns
+    if (radius == 45) {
+        radmax = 45.0;
+    } else if (radius == 75) {
+        radmax = 75.0;
+        radmin = 45.0;
+    } else if (radius == 40) {
+        // Placeholder for "rest of the sector"
+        // r < 40 mm has too little statistics anyways!
+        radmin = 75.0;
+        radmax = 2000.0;
+    }
   }
   
   fn_graph_s += ".pdf";
@@ -264,7 +279,7 @@
   TCut cu_sel = STTool::get_track_selection(det, fill);
   
   // Get full selection for signal tracks/hits including radius
-  TCut cu_sig = cu_bkg && cu_sel && TCut(Form("(TMath::Sqrt(xHit*xHit+yHit*yHit)<%f)",radius));
+  TCut cu_sig = cu_bkg && cu_sel && TCut(Form("(TMath::Sqrt(xHit*xHit+yHit*yHit)>%f && (TMath::Sqrt(xHit*xHit+yHit*yHit)<%f)", radmin, radmax));
 
   
   // Fill Histograms