lcsim-contrib/src/main/java/org/lcsim/contrib/Cassell/recon/analysis
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);
}
}
}