100pb file

We processed all the 100pb file and here is a glance at how they look like.
You can see the Energy distribution, the number of the electrons in a xtal and the EoP distribution.
Good night!

Sunday's news

my first post here!

committed alcaelectronstest.h and AlCaElectronsTest_filtered.root
I've tried to run with 100 events, every TH1F and TH2f looks to work properly.

we got

OBJ: TH1F pDistro_endcap pDistro_endcap : 0 at: 0x9421028
KEY: TH2F m_barrelGlobalCrystalsMap;1 m_barrelGlobalCrystalsMap
KEY: TH2F m_barrelLocalCrystalsMap;1 m_barrelLocalCrystalsMap
KEY: TH2F m_endcapGlobalCrystalsMap;1 m_endcapGlobalCrystalsMap
KEY: TH2F m_endcapLocalCrystalsMap;1 m_endcapLocalCrystalsMap
KEY: TH2F m_barrelGlobalCrystalsEnergy;1 m_barrelGlobalCrystalsEnergy
KEY: TH2F m_barrelLocalCrystalsEnergy;1 m_barrelLocalCrystalsEnergy
KEY: TH2F m_endcapGlobalCrystalsEnergy;1 m_endcapGlobalCrystalsEnergy
KEY: TH2F m_endcapLocalCrystalsEnergy;1 m_endcapLocalCrystalsEnergy
KEY: TH2F m_barrelGlobalCrystalsEnergyMap;1 m_barrelGlobalCrystalsEnergyMap
KEY: TH2F m_endcapGlobalCrystalsEnergyMap;1 m_endcapGlobalCrystalsEnergyMap
KEY: TProfile2D m_barrelGlobalCrystalsMomentum;1 m_barrelGlobalCrystalsMomentum
KEY: TH1F pDistro_barrel;1 pDistro_barrel
KEY: TH1F Edistro_barrel;1 Edistro_barrel
KEY: TH1F EoP_barrel;1 EoP_barrel
KEY: TProfile2D m_endcapGlobalCrystalsMomentum;1 m_endcapGlobalCrystalsMomentum
KEY: TH1F pDistro_endcap;1 pDistro_endcap
KEY: TH1F Edistro_endcap;1 Edistro_endcap
KEY: TH1F EoP_endcap;1 EoP_endcap

as a example:

Saturday's news

The 100pb have been processed and put where i told before on castor. They won't be on cmsmi3 till this evening.
422Kevents succeded while ~129K still needs to be reprocessed...Hoping in castor.

Data processing

The 10pb data have almost all been processed. Some are still not available on castor but we have 422Kevents. The are in /castor/cern.ch/user/p/presotto/10pb/
They are currently been transferred also in cmsmi3: /data/calib/csa07/alcaReco/10pb/
Each file is about 7MB.
We'll try to recover another ~100Kevents from the processes aborted because castor didn't feed the data fast enough for the LSF.

Now the batch system is processing, actually really fastly, the 100 pb datas.
They will be on /castor/cern.ch/user/p/presotto/100pb/ Then we'll transfer that on /data/calib/csa07/alcaReco/100pb/

As now we have 130Kevents processed, and they have been processed just in 30 minutes.
Tomorrow we'll process the frontier.

Then we'll run the alcaElectron test on all the data to see how does the data look like and well start deep testing the InvMatrixLooper.

Hi!

The draw of a Region

Now we tested the region Definition.
In this plot we show with different color even and odd regions.
The number are the id of the xtals in that specific region

CalibBlock's Test

This morning the EcalCalibBlock class underwent some test.

We filled a map of 25 entries with random number that acted like the "energy deposit"
We summed all the entries and this was the "momentum".
With this values we filled the calibBlocks and we repeated the procedure for many times. (with only one filling the system isn't determined)
That we invoked the solve method and then read back the coefficients. If the class is working properly the results should all be one.
And that's what happened.
Then we multiplied the entries of the map for some random coefficients, after summing this numbers for the "momentum". If the class was working the results should have been the inverse of our input coeffs. And that's what happened.
The program that does this test is in EcalAlCaRecoProducers/bin/BlockTest.cc

Ring definitions

Now we are ready.
Finally we solved all the majors errors in the InvRingCalib code about the first part that is about getting the geometry from the event setup and use it to define what is a ring in the Ecal Endcap.
Here are some graphs, as always we apologize for the horrible colours that root gives us.
In the cfg we specified all the rings (from 0 to 39) with regions 4 ring wide.
In the barrel eta start 3, eta end 159, 4 rings wide.
In the endcap 0(green) is where there aren't xtals, 10 in some regions and -10 in the following ones.