Print

Print


Commit in lcsim-contrib/src/main/java/org/lcsim/contrib/Cassell/recon/analysis on MAIN
SPFinderSingleParticleTest.java+511.1 -> 1.2
Divide by Energy

lcsim-contrib/src/main/java/org/lcsim/contrib/Cassell/recon/analysis
SPFinderSingleParticleTest.java 1.1 -> 1.2
diff -u -r1.1 -r1.2
--- SPFinderSingleParticleTest.java	6 Oct 2010 19:09:48 -0000	1.1
+++ SPFinderSingleParticleTest.java	6 Oct 2010 19:35:03 -0000	1.2
@@ -45,6 +45,7 @@
     int mx = 10000000;
     int nmax = 10000000;
     String[] pre0 = {"Barrel/","Cross/","Endcap/"};
+    double[] Evals = {1.,2.,5.,10.,20.,50.,100.};
     public SPFinderSingleParticleTest()
     {
 //       add(new ReconDriver());
@@ -94,6 +95,9 @@
            ntr++;
            return;
         }
+       double E = thisp.getEnergy();
+       int Est = (int) (E + .5);
+       String pre2 = "E="+Est+"/";
         List<TrackXtrapInfo> txil = new ArrayList<TrackXtrapInfo>(event.get(TrackXtrapInfo.class,"TrackXtrapInfo"));
         TrackXtrapInfo txi = txil.get(0);
          List<XtrapPosition> xpl = txi.getPositions();
@@ -141,6 +145,7 @@
            if(name.contains("ShowerFinder"))
            {
               String ppre = name +"/";
+              String ppre0 = ppre+pre2;
               if(cl.size() != 1)
               {
                  aida.cloud1D(ppre+"No mip/Endpoint layer",nmax).fill(lay);
@@ -155,12 +160,29 @@
                  aida.cloud1D(ppre+pre0+pre1+"No mip/Endpoint step",nmax).fill(lastind);
                  aida.cloud1D(ppre+pre0+pre1+"No mip/MC energy",nmax).fill(thisp.getEnergy());
                  aida.cloud1D(ppre+pre0+pre1+"No mip/MC cos theta",nmax).fill(ct);
+
+                 aida.cloud1D(ppre0+"No mip/Endpoint layer",nmax).fill(lay);
+                 aida.cloud1D(ppre0+"No mip/Endpoint step",nmax).fill(lastind);
+                 aida.cloud1D(ppre0+"No mip/MC energy",nmax).fill(thisp.getEnergy());
+                 aida.cloud1D(ppre0+"No mip/MC cos theta",nmax).fill(ct);
+                 aida.cloud1D(ppre0+pre0+"No mip/Endpoint layer",nmax).fill(lay);
+                 aida.cloud1D(ppre0+pre0+"No mip/Endpoint step",nmax).fill(lastind);
+                 aida.cloud1D(ppre0+pre0+"No mip/MC energy",nmax).fill(thisp.getEnergy());
+                 aida.cloud1D(ppre0+pre0+"No mip/MC cos theta",nmax).fill(ct);
+                 aida.cloud1D(ppre0+pre0+pre1+"No mip/Endpoint layer",nmax).fill(lay);
+                 aida.cloud1D(ppre0+pre0+pre1+"No mip/Endpoint step",nmax).fill(lastind);
+                 aida.cloud1D(ppre0+pre0+pre1+"No mip/MC energy",nmax).fill(thisp.getEnergy());
+                 aida.cloud1D(ppre0+pre0+pre1+"No mip/MC cos theta",nmax).fill(ct);
                  continue;
               }
               int nhim = cl.get(0).getCalorimeterHits().size();
               aida.cloud1D(ppre+"# hits in mip",nmax).fill(nhim);
               aida.cloud1D(ppre+pre0+"# hits in mip",nmax).fill(nhim);
               aida.cloud1D(ppre+pre0+pre1+"# hits in mip",nmax).fill(nhim);
+
+              aida.cloud1D(ppre0+"# hits in mip",nmax).fill(nhim);
+              aida.cloud1D(ppre0+pre0+"# hits in mip",nmax).fill(nhim);
+              aida.cloud1D(ppre0+pre0+pre1+"# hits in mip",nmax).fill(nhim);
               double maxR = 0.;
               CalorimeterHit lh = null;
               Hep3Vector lp = new BasicHep3Vector();
@@ -195,6 +217,10 @@
               aida.cloud1D(ppre+"# Recon particles",nmax).fill(nrp);
               aida.cloud1D(ppre+pre0+"# Recon particles",nmax).fill(nrp);
               aida.cloud1D(ppre+pre0+pre1+"# Recon particles",nmax).fill(nrp);
+
+              aida.cloud1D(ppre0+"# Recon particles",nmax).fill(nrp);
+              aida.cloud1D(ppre0+pre0+"# Recon particles",nmax).fill(nrp);
+              aida.cloud1D(ppre0+pre0+pre1+"# Recon particles",nmax).fill(nrp);
               int dlay = foundlay - lay;
               int dstep = foundstep - lastind;
               aida.cloud1D(ppre+"Endpoint layer",nmax).fill(lay);
@@ -221,6 +247,31 @@
               aida.cloud1D(ppre+pre0+pre1+"Found step",nmax).fill(foundstep);
               aida.cloud1D(ppre+pre0+pre1+"Delta step",nmax).fill(dstep);
               aida.cloud2D(ppre+pre0+pre1+"Dstep vs Estep",nmax).fill(lastind,dstep);
+
+              aida.cloud1D(ppre0+"Endpoint layer",nmax).fill(lay);
+              aida.cloud1D(ppre0+"Found layer",nmax).fill(foundlay);
+              aida.cloud1D(ppre0+"Delta layer",nmax).fill(dlay);
+              aida.cloud2D(ppre0+"Dlay vs Elay",nmax).fill(lay,dlay);
+              aida.cloud1D(ppre0+"Endpoint step",nmax).fill(lastind);
+              aida.cloud1D(ppre0+"Found step",nmax).fill(foundstep);
+              aida.cloud1D(ppre0+"Delta step",nmax).fill(dstep);
+              aida.cloud2D(ppre0+"Dstep vs Estep",nmax).fill(lastind,dstep);
+              aida.cloud1D(ppre0+pre0+"Endpoint layer",nmax).fill(lay);
+              aida.cloud1D(ppre0+pre0+"Found layer",nmax).fill(foundlay);
+              aida.cloud1D(ppre0+pre0+"Delta layer",nmax).fill(dlay);
+              aida.cloud2D(ppre0+pre0+"Dlay vs Elay",nmax).fill(lay,dlay);
+              aida.cloud1D(ppre0+pre0+"Endpoint step",nmax).fill(lastind);
+              aida.cloud1D(ppre0+pre0+"Found step",nmax).fill(foundstep);
+              aida.cloud1D(ppre0+pre0+"Delta step",nmax).fill(dstep);
+              aida.cloud2D(ppre0+pre0+"Dstep vs Estep",nmax).fill(lastind,dstep);
+              aida.cloud1D(ppre0+pre0+pre1+"Endpoint layer",nmax).fill(lay);
+              aida.cloud1D(ppre0+pre0+pre1+"Found layer",nmax).fill(foundlay);
+              aida.cloud1D(ppre0+pre0+pre1+"Delta layer",nmax).fill(dlay);
+              aida.cloud2D(ppre0+pre0+pre1+"Dlay vs Elay",nmax).fill(lay,dlay);
+              aida.cloud1D(ppre0+pre0+pre1+"Endpoint step",nmax).fill(lastind);
+              aida.cloud1D(ppre0+pre0+pre1+"Found step",nmax).fill(foundstep);
+              aida.cloud1D(ppre0+pre0+pre1+"Delta step",nmax).fill(dstep);
+              aida.cloud2D(ppre0+pre0+pre1+"Dstep vs Estep",nmax).fill(lastind,dstep);
            }
         }
    }
CVSspam 0.2.8