java/trunk/users/src/main/java/org/hps/users/luca
--- java/trunk/users/src/main/java/org/hps/users/luca/CalibTest.java 2014-04-29 11:22:34 UTC (rev 524)
+++ java/trunk/users/src/main/java/org/hps/users/luca/CalibTest.java 2014-04-29 12:36:27 UTC (rev 525)
@@ -3,6 +3,7 @@
//import hep.aida.ITupleColumn.String;
import java.io.IOException;
+import java.io.*;
import java.util.ArrayList;
import java.util.Collection;
import java.util.HashMap;
@@ -14,6 +15,7 @@
import org.hps.recon.ecal.HPSEcalCluster;
import org.lcsim.event.CalorimeterHit;
import org.lcsim.event.EventHeader;
+import org.lcsim.event.ParticleID;
import org.lcsim.geometry.Detector;
import org.lcsim.geometry.subdetector.HPSEcal3;
import org.lcsim.geometry.subdetector.HPSEcal3.NeighborMap;
@@ -31,72 +33,58 @@
// the class has to be derived from the driver class
public class CalibTest extends Driver {
+int counter=0;
// constructor
public CalibTest() {
}
- String ecalName;
- String ecalCollectionName;
- String clusterCollectionName;
- String calorhit;
+
// overwrite the process method
@Override
protected void process(EventHeader event) {
// Get the list of mc particles from the event
- //List<MCParticle> mcParticles = event.getMCParticles();
+ List<MCParticle> mcParticles = event.getMCParticles();
// Print out the number of mc particles
//System.out.println("Event " + event.getEventNumber() + " contains " + mcParticles.size() + " mc particles.");
- // get the list of the lists listed in the event list
- //Set<List> myfirstlist = event.getLists();
- // System.out.println("this is the number of lists: "+ event.getLists().size());
- if(event.hasCollection(SimCalorimeterHit.class,calorhit))
- { List<SimCalorimeterHit> myCalorHit =event.getSimCalorimeterHits(calorhit);
- System.out.println("This file has " + myCalorHit.size() + "SimClaorimeterHits object");}
-
-
-
-
-
-
- /*for(List myhit : myCalorHit)
- { for(Object hit : myhit)
- {System.out.println("this event has energy: "+ hit.getRawEnergy() + "\n");} }*/
-
-
- /*for(List lista : myfirstlist )
- {
- System.out.println("This event contains this list: " + lista + "\n");}*/
- /* List<Track> myTrack = event.getTracks();
- System.out.println("this is what we get printing the tracks: "+ event.getTracks().get(event.getEventNumber()));*/
-
- //List<Cluster> myCluster = event.getClusters();
- /* System.out.println("this is what we get printing the Clusters: "+ myCluster.get(event.getEventNumber()));*/
-
- /* if (event.hasCollection(SimCalorimeterHit.class,ecalName)){
- List<SimCalorimeterHit> mySimCalHits = event.getSimCalorimeterHits(ecalName);
- System.out.println("this is what we get printing the calorimeter hits: "+ mySimCalHits.size());
- }
-
-
- /* List<SimTrackerHit> mySimTrackerHits = event.getSimTrackerHits(String string);
- System.out.println("this is what we get printing the tracker hits: "+ mySimTrackerHits.size());*/
- /* if (event.hasCollection(CalorimeterHit.class, ecalCollectionName)) {
- // Get the list of raw ECal hits.
- List<CalorimeterHit> hits = event.get(CalorimeterHit.class, ecalCollectionName);
-
- // Make a hit map for quick lookup by ID.
- Map<Long, CalorimeterHit> hitMap = new HashMap<Long, CalorimeterHit>();
+ for (MCParticle particle : mcParticles)
+ {
+ if(particle.getPDGID()==11)
+ {if(particle.getEnergy()> 2.1){
+ counter++;
+ }
+ }
+
+ //System.out.println(particle.getPDGID());
+ }
+System.out.println("ho contato" + counter + "elettroni. \n");
+ /* try
+ {
+ FileOutputStream prova = new FileOutputStream("prova.txt");
+ PrintStream scrivi = new PrintStream(prova);*/
+ //if (event.hasCollection(CalorimeterHit.class, "EcalHits")) {
+ // List<CalorimeterHit> hits = event.get(CalorimeterHit.class, "EcalHits");
+
+ // for (CalorimeterHit hit : hits) {
+ // energia=hit.getRawEnergy();
+ //scrivi.print(energia );
- for (CalorimeterHit hit : hits) {
- hitMap.put(hit.getCellID(), hit);
- }
- System.out.println("Number of ECal hits: "+hitMap.size());
-
+ // }
- }*/
+
+
+
+ // }
+ // }
+ /* catch (IOException e)
+ {
+ System.out.println("Errore: " + e);
+ System.exit(1);
+ }*/
+
}
+
}
\ No newline at end of file