Author: [log in to unmask]
Date: Tue May 12 03:45:02 2015
New Revision: 2959
Log: (empty)
Modified:
java/trunk/users/src/main/java/org/hps/users/luca/ratetest.java
Modified: java/trunk/users/src/main/java/org/hps/users/luca/ratetest.java
=============================================================================
--- java/trunk/users/src/main/java/org/hps/users/luca/ratetest.java (original)
+++ java/trunk/users/src/main/java/org/hps/users/luca/ratetest.java Tue May 12 03:45:02 2015
@@ -24,10 +24,13 @@
import org.hps.conditions.ecal.EcalChannel.EcalChannelCollection;
import org.hps.conditions.ecal.EcalConditions;
+import org.hps.recon.ecal.triggerbank.AbstractIntData;
+import org.hps.recon.ecal.triggerbank.TIData;
import org.lcsim.event.CalorimeterHit;
import org.lcsim.event.Cluster;
import org.lcsim.event.EventHeader;
+import org.lcsim.event.GenericObject;
import org.lcsim.geometry.Detector;
import org.lcsim.geometry.Subdetector;
import org.lcsim.geometry.subdetector.HPSEcal3;
@@ -46,105 +49,16 @@
private Subdetector ecal;
private EcalChannelCollection channels= null;
AIDA aida = AIDA.defaultInstance();
- ArrayList<IHistogram1D> Top = new ArrayList<IHistogram1D>(6);
- ArrayList<IHistogram1D> Bot = new ArrayList<IHistogram1D>(6);
- ArrayList<IHistogram1D> TopNoCut = new ArrayList<IHistogram1D>(6);
- ArrayList<IHistogram1D> BotNoCut = new ArrayList<IHistogram1D>(6);
- ArrayList<IHistogram2D> TopCry=new ArrayList<IHistogram2D>(6);
- ArrayList<IHistogram2D> BotCry=new ArrayList<IHistogram2D>(6);
- ArrayList<IHistogram2D> CryTop=new ArrayList<IHistogram2D>(5);
- ArrayList<IHistogram2D> CryBot=new ArrayList<IHistogram2D>(5);
- ArrayList<IHistogram1D> TopCryEne = new ArrayList<IHistogram1D>(5);
- ArrayList<IHistogram1D> BotCryEne = new ArrayList<IHistogram1D>(5);
- IHistogram1D clusTot=aida.histogram1D("ALL", 200,0,2.5 );
- IHistogram1D CryId=aida.histogram1D("cry Id",446,1,447);
- double[] timei=new double[6];
- double timef=0;
- double counterTop=0;
- double counterBot=0;
- double counterEcut=0;
+ IHistogram1D c110=aida.histogram1D("110", 200,0,2.5 );
+ IHistogram1D c18=aida.histogram1D("18", 200,0,2.5 );
+ IHistogram1D c60=aida.histogram1D("60", 200,0,2.5 );
+ IHistogram1D c322=aida.histogram1D("322", 200,0,2.5 );
+ IHistogram1D c414=aida.histogram1D("414", 200,0,2.5 );
+ IHistogram1D c364=aida.histogram1D("364", 200,0,2.5 );
- double nevents=0;
-double[] time=new double[6];
-double E0=1.92;
-double[] Ymax={45.91,52.95,60.00,67.01,74.04,81.18};
-double[] Ymin={38.83,45.81,52.95,60.00,67.01,74.04};
-double[] Xmax={-23.59,-23.1,-22.61,-22.11,-21.62,-21.13};
-double[] Xmin={-30.1,-30.59,-31.08,-31.58,-32.07,-32.57};
-
-int[] countertop=new int[6];
-int[] counterbot=new int[6];
-int[]crycountertop=new int[5];
-int[]crycounterbot=new int[5];
-
-double[] vertpos=new double[5];
-
-//e cuts 3393
-
-double[] ecuttopmin={1.20,1.35,1.50,1.50,1.60,1.60};
-double[] ecuttopmax={1.70,1.90,1.90,2.00,2.00,2.00};
-
-double[] ecutbotmin={1.28,1.50,1.50,1.65,1.65,1.50};
-double[] ecutbotmax={1.90,1.90,2.00,2.00,2.00,2.00};
-
-
-double[] cryecuttopmin={1.38,1.38,1.52,1.6,1.6};
-double[] cryecuttopmax={1.8,1.8,1.8,2,2};
-
-double[] cryecutbotmax={2,2,1.9,1.9,2};
-double[] cryecutbotmin= {1.38,1.38,1.5,1.65,1.55};
-
-
-//e cuts 3430
-/*
-double[] ecuttopmin={1.20,1.35,1.50,1.50,1.60,1.60};
-double[] ecuttopmax={1.70,1.80,1.90,2.00,2.00,2.00};
-
-double[] ecutbotmin= {1.20,1.50,1.50,1.55,1.60,1.50};
-double[] ecutbotmax={1.80,1.90,1.90,2.00,2.00,2.00};
-
-
-double[] cryecuttopmin={1.38,1.38,1.52,1.6,1.6};
-double[] cryecuttopmax={1.8,1.8,1.8,2,2};
-
-double[] cryecutbotmax={2,2,1.9,1.9,2};
-double[] cryecutbotmin= {1.38,1.38,1.5,1.65,1.55};
-*/
-//e cuts 3434
-/*
-double[] ecuttopmin={1.25,1.30,1.40,1.50,1.60,1.60};
-double[] ecuttopmax={1.65,1.80,1.90,1.90,2.00,2.00};
-
-double[] ecutbotmin= {1.30,1.40,1.40,1.60,1.60,1.50};
-double[] ecutbotmax={1.70,1.80,1.90,2,2,2};
-
-
-double[] cryecuttopmin={1.38,1.38,1.52,1.6,1.6};
-double[] cryecuttopmax={1.8,1.8,1.8,2,2};
-
-double[] cryecutbotmax={2,2,1.9,1.9,2};
-double[] cryecutbotmin= {1.38,1.38,1.5,1.65,1.55};
-*/
-//e cuts 3435
-/*
-double[] ecuttopmin={1.28,1.30,1.40,1.50,1.50,1.60};
-double[] ecuttopmax={1.7,1.7,1.8,1.8,2,2};
-
-double[] ecutbotmin= {1.3,1.35,1.40,1.50,1.60,1.50};
-double[] ecutbotmax={1.7,1.8,1.9,2,2,2};
-
-
-double[] cryecuttopmin={1.38,1.38,1.52,1.6,1.6};
-double[] cryecuttopmax={1.8,1.8,1.8,2,2};
-
-double[] cryecutbotmax={2,2,1.9,1.9,2};
-double[] cryecutbotmin= {1.38,1.38,1.5,1.65,1.55};
-*/
-
-double vertical;
@@ -158,77 +72,7 @@
@Override
public void startOfData(){
- System.out.println("mo spacco tutto \n");
- //inizializzo il file di uscita
- try{
- //initialize the writers
- writer=new FileWriter(outputFileName);
- writer.write("");
- }
- catch(IOException e ){
- System.err.println("Error initializing output file for event display.");
- }
-
-//inizializzo istogrammi
- for(int t=0;t<5;t++){
- String bin=String.valueOf(t+1);
- String crytopname="FEE in Crystal " + bin;
- String topcryenename="CLu ene in crystal " + bin;
-
- String bin2=String.valueOf(-(t+1));
- String crybotname="FEE in Crystal " + bin2;
- String botcryenename="CLu ene in crystal " + bin2;
- IHistogram2D crytophist=aida.histogram2D(crytopname,46,-24,25,10,-5,5 );
- IHistogram2D crybothist=aida.histogram2D(crybotname,46,-24,25,10,-5,5 );
- IHistogram1D Topcryenehist=aida.histogram1D(topcryenename, 250, 0.0,2.5);
- IHistogram1D botcryenehist=aida.histogram1D(botcryenename,250,0.0,2.5);
- TopCryEne.add(Topcryenehist);
- BotCryEne.add(botcryenehist);
- CryTop.add(crytophist);
- CryBot.add(crybothist);
- crycountertop[t]=0;
- crycounterbot[t]=0;
- }
-
-
- for(int t=0; t<6; t++){
- String bin=String.valueOf(t+1);
- String top="(TOP) FEE in Bin in " + bin;
- String bot="(BOT) FEE in Bin in "+ bin;
- String tope="(TOP) Cluster Energy in Bin " + bin;
- String bote="(BOT) Cluster Energy in Bin "+ bin;
- String topcryname="(TOP) Crystals in Bin " + bin;
- String botcryname="(BOT) Crystals in Bin " + bin;
-
-
-
- IHistogram1D Toppe=aida.histogram1D(top, 250, 0.0,2.5);
- IHistogram1D Botte=aida.histogram1D(bot, 250, 0.0,2.5);
- IHistogram1D TopNoCute=aida.histogram1D(tope, 250, 0.0,2.5);
- IHistogram1D BotNoCute=aida.histogram1D(bote, 250, 0.0,2.5);
- IHistogram2D topcryhist=aida.histogram2D(topcryname,46,-24,25,11,-5,6 );
- IHistogram2D botcryhist=aida.histogram2D(botcryname,46,-24,25,11,-5,6 );
-
-
- Top.add(Toppe);
- Bot.add(Botte);
- TopNoCut.add(TopNoCute);
- BotNoCut.add(BotNoCute);
- TopCry.add(topcryhist);
- BotCry.add(botcryhist);
- }
-
-
-
-for (int i =0;i<6;i++){
- time[i]=0;
- timei[i]=0;
-
- countertop[i]=0;
- counterbot[i]=0;
-
-}
}
@@ -236,144 +80,42 @@
@Override
public void endOfData(){
-double timme=(timef-timei[0])/1000000000;
-
-//for(int i=0;i<6;i++){System.out.println(i + " " + time[i] + "\n");}
-
-try{
-writer.append("Total event in file " + nevents + " Total time = "+ timme +"\n");
-writer.append("TOP \n");
-for(int i=0;i<6;i++){
- int bin=i+1;
- time[i]=timef-timei[i];
-writer.append("Bin " + bin + " Total event " + countertop[i] + " Rate = " + (countertop[i]/time[i])*1000000000 + "\n" );
-}
-for(int i=0;i<6;i++){
- int bin=i+1;
- time[i]=timef-timei[i];
-writer.append("Bin " + bin +" Total event BOT " + counterbot[i] + " Rate = " + (counterbot[i]/time[i])*1000000000 + "\n" );
-}
-
-writer.append("\n \n \n");
-for(int t=0;t<5;t++){
-writer.append("events in row " + t + " = " + crycountertop[t]+ " Rate = " + (crycountertop[t]/timme) + "\n" );
-int menot=-t;
-writer.append("events in row " + menot + " = " + crycounterbot[t]+ " Rate = " + (crycounterbot[t]/timme) + "\n" );
-}
-writer.append("\n verticale = " + vertical +"\n");
-
-}
-
-
-catch(IOException e)
- {
- System.err.println("Non ho scritto sul file");
- }
-
-System.out.println("vertical = " + vertical + "\n");
-/*vertical=vertical-6.5;
-for(int i=0;i<7;i++)
-{double theta=Math.atan2( (vertical+i*13),1397);
-System.out.println("theta = " + theta + "\n");
-}
-*/
-
- try
- {
-//close the file writer.
- writer.close();
- }
- catch(IOException e)
- {
- System.err.println("Error closing utput file for event display.");
- }
+
}
@Override
public void process (EventHeader event){
- timef=event.getTimeStamp();
+ if (event.hasCollection(GenericObject.class, "TriggerBank")) {
+ List<GenericObject> triggerList = event.get(GenericObject.class, "TriggerBank");
+ for (GenericObject data : triggerList)
+ if (AbstractIntData.getTag(data) == TIData.BANK_TAG) {
+ TIData triggerData = new TIData(data);
+ if (!triggerData.isSingle1Trigger())//only process singles0 triggers...
+
+ return;
+ }
+ } else //if (debug)
+ System.out.println(this.getClass().getSimpleName() + ": No trigger bank found...running over all trigger types");
- nevents++;
-
- for(int i=0;i<6;i++){
- if(countertop[i]==0){timei[i]=timef;}}
-
- if(event.hasCollection(Cluster.class,"EcalClusters")){
+
+ if(event.hasCollection(Cluster.class,"EcalClusters")){
List<Cluster> clusters= event.get(Cluster.class,"EcalClusters");
for(Cluster cluster : clusters){
-
+ int id=getDBID(cluster);
+ if(cluster.getEnergy()>0.65 && cluster.getSize()==1){
+ if(id==110){c110.fill(cluster.getEnergy());}
- int ID=getDBID(cluster);
- clusTot.fill(cluster.getEnergy());
- double posY=cluster.getPosition()[1];
- double posX=cluster.getPosition()[0];
- double cryx=cluster.getCalorimeterHits().get(0).getIdentifierFieldValue("ix");
- double cryy=cluster.getCalorimeterHits().get(0).getIdentifierFieldValue("iy");
-
- //if(cryy==1){verpos[0]=posY;}
- //da qui righe e colonne
- if(cryx==-5||cryx==-6){
-
- for(int i=0;i<5;i++){
- //top
- if(cryy==i+1){
- if(cluster.getEnergy()>cryecuttopmin[i]&&cluster.getEnergy()<cryecuttopmax[i]){
- CryTop.get(i).fill(cryx, cryy);
- TopCryEne.get(i).fill(cluster.getEnergy());
- crycountertop[i]++;
- }
- }
-
- //bottom
- if(cryy==-(i+1)){
- if(cluster.getEnergy()>cryecutbotmin[i] && cluster.getEnergy()<cryecutbotmax[i]){
- CryBot.get(i).fill(cryx, cryy);
- BotCryEne.get(i).fill(cluster.getEnergy());
- crycounterbot[i]++;
- }
- }
- }
- }
-
-
- ///da qui bin angolari
-
- for(int i=0;i<6;i++){
- //top
- if(posY>Ymin[i]&&posY<=Ymax[i]){
- if(posX<=Xmax[i] &&posX>=Xmin[i]){
- TopNoCut.get(i).fill(cluster.getEnergy());
- if(cluster.getEnergy()>ecuttopmin[i]&&cluster.getEnergy()<ecuttopmax[i]){
- Top.get(i).fill(cluster.getEnergy());
- countertop[i]++;
- CryId.fill(ID);
- TopCry.get(i).fill(cluster.getCalorimeterHits().get(0).getIdentifierFieldValue("ix"), cluster.getCalorimeterHits().get(0).getIdentifierFieldValue("iy"));
- }//tagli energia
- } //if su X
- }//if pos y
+ }
- //controllo down
- else if(posY<= -Ymin[i]&& posY> -Ymax[i]){
- if(posX<Xmax[i] && posX>Xmin[i]){
- BotNoCut.get(i).fill(cluster.getEnergy());
- if(cluster.getEnergy()>ecutbotmin[i] && cluster.getEnergy()<ecutbotmax[i]){
- Bot.get(i).fill(cluster.getEnergy());
- counterbot[i]++;
- BotCry.get(i).fill(cluster.getCalorimeterHits().get(0).getIdentifierFieldValue("ix"), cluster.getCalorimeterHits().get(0).getIdentifierFieldValue("iy"));
- CryId.fill(ID);
- }//e cut bot
- } //xpos
- }// ypos bot
- }//for sui bin
-
- }//end of for over clusters
- } //and of if has colelction
+ }
+ }//end of if has collection
+
}//end pf process
