cms::Analyzer_minbias Class Reference

#include <Analyzer_minbias.h>

Inheritance diagram for cms::Analyzer_minbias:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
	Analyzer_minbias (const edm::ParameterSet &)
virtual void	beginJob ()
virtual void	beginRun (const edm::Run &r, const edm::EventSetup &iSetup)
virtual void	endJob ()
virtual void	endRun (const edm::Run &r, const edm::EventSetup &iSetup)
	~Analyzer_minbias ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Private Attributes
int	depth
float	eta
std::string	fOutputFileName
edm::InputTag	hbherecoMB
edm::InputTag	hbherecoNoise
std::string	hcalfile_
TH1F *	hCalo1 [73][43]
TH1F *	hCalo1mom2 [73][43]
TH1F *	hCalo2 [73][43]
TH1F *	hCalo2mom2 [73][43]
edm::InputTag	hfrecoMB
edm::InputTag	hfrecoNoise
TH2F *	hHBHEsize_vs_run
TH2F *	hHFsize_vs_run
edm::InputTag	horecoMB
edm::InputTag	horecoNoise
TFile *	hOutputFile
int	ieta
int	iphi
double	meannoise_min [73][43]
double	meannoise_pl [73][43]
float	mom0_Diff
float	mom0_MB
float	mom0_Noise
float	mom1_Diff
float	mom1_MB
float	mom1_Noise
float	mom2_Diff
float	mom2_MB
float	mom2_Noise
float	mom3_Diff
float	mom3_MB
float	mom3_Noise
float	mom4_Diff
float	mom4_MB
float	mom4_Noise
int	mydet
std::ofstream *	myout_hcal
int	mysubd
TTree *	myTree
double	nevent
int	nevent_run
double	noise_min [73][43]
double	noise_pl [73][43]
float	occup
float	phi
double	theDFFillDetMapMin0 [5][5][73][43]
double	theDFFillDetMapMin1 [5][5][73][43]
double	theDFFillDetMapMin2 [5][5][73][43]
double	theDFFillDetMapPl0 [5][5][73][43]
double	theDFFillDetMapPl1 [5][5][73][43]
double	theDFFillDetMapPl2 [5][5][73][43]
double	theMBFillDetMapMin0 [5][5][73][43]
double	theMBFillDetMapMin1 [5][5][73][43]
double	theMBFillDetMapMin2 [5][5][73][43]
double	theMBFillDetMapMin4 [5][5][73][43]
double	theMBFillDetMapPl0 [5][5][73][43]
double	theMBFillDetMapPl1 [5][5][73][43]
double	theMBFillDetMapPl2 [5][5][73][43]
double	theMBFillDetMapPl4 [5][5][73][43]
double	theNSFillDetMapMin0 [5][5][73][43]
double	theNSFillDetMapMin1 [5][5][73][43]
double	theNSFillDetMapMin2 [5][5][73][43]
double	theNSFillDetMapMin4 [5][5][73][43]
double	theNSFillDetMapPl0 [5][5][73][43]
double	theNSFillDetMapPl1 [5][5][73][43]
double	theNSFillDetMapPl2 [5][5][73][43]
double	theNSFillDetMapPl4 [5][5][73][43]
bool	theRecalib

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 51 of file Analyzer_minbias.h.

Constructor & Destructor Documentation

cms::Analyzer_minbias::Analyzer_minbias ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 33 of file Analyzer_minbias.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hcalLocalRecoNZS_cff::hbherecoMB, ALCARECOHcalCalMinBias_cff::hbherecoNoise, hcalLocalRecoNZS_cff::hfrecoMB, ALCARECOHcalCalMinBias_cff::hfrecoNoise, hcalLocalRecoNZS_cff::horecoMB, ALCARECOHcalCalMinBias_cff::horecoNoise, i, and j.

{
  // get name of output file with histogramms
  fOutputFileName = iConfig.getUntrackedParameter<string>("HistOutFile"); 
  // get names of modules, producing object collections
  
  hbherecoMB = iConfig.getParameter<edm::InputTag>("hbheInputMB");
  horecoMB   = iConfig.getParameter<edm::InputTag>("hoInputMB");
  hfrecoMB   = iConfig.getParameter<edm::InputTag>("hfInputMB");
  
  hbherecoNoise = iConfig.getParameter<edm::InputTag>("hbheInputNoise");
  horecoNoise   = iConfig.getParameter<edm::InputTag>("hoInputNoise");
  hfrecoNoise   = iConfig.getParameter<edm::InputTag>("hfInputNoise");
  
  theRecalib = iConfig.getParameter<bool>("Recalib"); 
   
//
//
  for(int i=0; i<73; i++)
  {
     for(int j=0; j<43; j++)
     {
        noise_min[i][j] = 0.;
    noise_pl[i][j] = 0.;
     } 
  }
//
//
     
}

cms::Analyzer_minbias::~Analyzer_minbias

(

)

Definition at line 64 of file Analyzer_minbias.cc.

{
 
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)

}

Member Function Documentation

void cms::Analyzer_minbias::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

virtual

Implements edm::EDAnalyzer.

Definition at line 329 of file Analyzer_minbias.cc.

References abs, gather_cfg::cout, HcalDetId::depth(), CaloRecHit::energy(), HcalObjRepresent::Fill(), edm::EventSetup::get(), edm::Event::getByLabel(), hcalLocalRecoNZS_cff::hbherecoMB, ALCARECOHcalCalMinBias_cff::hbherecoNoise, hcalLocalRecoNZS_cff::hfrecoMB, ALCARECOHcalCalMinBias_cff::hfrecoNoise, i, HcalDetId::ieta(), HcalDetId::iphi(), j, gen::k, prof2calltree::l, LogDebug, nevent, funct::pow(), edm::ESHandle< class >::product(), DetId::rawId(), mathSSE::return(), and edm::Event::run().

{

  std::cout<<" Start Analyzer_minbias::analyze "<<nevent<<std::endl;
  nevent++; 
  nevent_run++;
  using namespace edm;

  float rnnum = (float)iEvent.run(); 
/*
  std::vector<Provenance const*> theProvenance;
  iEvent.getAllProvenance(theProvenance);

  for( std::vector<Provenance const*>::const_iterator ip = theProvenance.begin();
                                                      ip != theProvenance.end(); ip++)
  {
     cout<<" Print all module/label names "<<(**ip).moduleName()<<" "<<(**ip).moduleLabel()<<
     " "<<(**ip).productInstanceName()<<endl;
  }
*/

   edm::Handle<L1GlobalTriggerReadoutRecord> gtRecord;
   iEvent.getByLabel("gtDigis", gtRecord);
   
   if (!gtRecord.isValid()) {

//     LogDebug("L1GlobalTriggerRecordProducer")
//       << "\n\n Error: no L1GlobalTriggerReadoutRecord found with input tag "
//       << m_l1GtReadoutRecord
//       << "\n Returning empty L1GlobalTriggerRecord.\n\n"
//       << std::endl;
      cout<<" No L1 trigger record "<<endl;
//     return;
   }


  const HcalRespCorrs* myRecalib=0;
  if( theRecalib ) {
// Radek:   
  edm::ESHandle <HcalRespCorrs> recalibCorrs;
  iSetup.get<HcalRespCorrsRcd>().get("recalibrate",recalibCorrs);
  myRecalib = recalibCorrs.product();
// end
  } // theRecalib

// Noise part for HB HE

     double tmpNSFillDetMapPl1[5][5][73][43]; 
     double tmpNSFillDetMapPl2[5][5][73][43];
     double tmpNSFillDetMapMin1[5][5][73][43]; 
     double tmpNSFillDetMapMin2[5][5][73][43];

   for (int i=0; i<5;i++)
   {
    for (int j=0; j<5;j++)
    {
     for (int k=0; k<73;k++)
     {
       for (int l=0; l<43;l++)
       {
        tmpNSFillDetMapPl1[i][j][k][l] = 0.;
        tmpNSFillDetMapPl2[i][j][k][l] = 0.;
        tmpNSFillDetMapMin1[i][j][k][l] = 0.;
        tmpNSFillDetMapMin2[i][j][k][l] = 0.;
       }
     }  
    }
   }    
   edm::Handle<HBHERecHitCollection> hbheNormal;
   iEvent.getByLabel("hbhereco", hbheNormal);
   if(!hbheNormal.isValid()){  
     cout<<" hbheNormal failed "<<endl;
   } else {
       cout<<" The size of the normal collection "<<hbheNormal->size()<<endl;
   } 


   edm::Handle<HBHERecHitCollection> hbheNS;
   iEvent.getByLabel(hbherecoNoise, hbheNS);


   if(!hbheNS.isValid()){
     LogDebug("") << "HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
     cout<<" No HBHE MS "<<endl;
     return ;
   }

  
   const HBHERecHitCollection HithbheNS = *(hbheNS.product());
   cout<<" HBHE NS size of collection "<<HithbheNS.size()<<endl;
   hHBHEsize_vs_run->Fill(rnnum,(float)HithbheNS.size());

   if(HithbheNS.size()!= 5184) {
          cout<<" HBHE problem "<<rnnum<<" "<<HithbheNS.size()<<endl;
          return;
   }
   edm::Handle<HBHERecHitCollection> hbheMB;
   iEvent.getByLabel(hbherecoMB, hbheMB);

   if(!hbheMB.isValid()){
     LogDebug("") << "HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
     cout<<" No HBHE MB"<<endl;
     return ;
   }

  const HBHERecHitCollection HithbheMB = *(hbheMB.product());
  cout<<" HBHE MB size of collection "<<HithbheMB.size()<<endl;
  if(HithbheMB.size()!= 5184) {
     cout<<" HBHE problem "<<rnnum<<" "<<HithbheMB.size()<<endl;
     return;
  }

   edm::Handle<HFRecHitCollection> hfNS;
   iEvent.getByLabel(hfrecoNoise, hfNS);

   if(!hfNS.isValid()){
     LogDebug("") << "HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
     cout<<" No HF NS "<<endl;
     return ;
   }

   const HFRecHitCollection HithfNS = *(hfNS.product());
   cout<<" HFE NS size of collection "<<HithfNS.size()<<endl;
   hHFsize_vs_run->Fill(rnnum,(float)HithfNS.size());
   if(HithfNS.size()!= 1728) {
          cout<<" HF problem "<<rnnum<<" "<<HithfNS.size()<<endl;
          return;
   }

   edm::Handle<HFRecHitCollection> hfMB;
   iEvent.getByLabel(hfrecoMB, hfMB);

   if(!hfMB.isValid()){
     LogDebug("") << "HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
     cout<<" No HBHE MB"<<endl;
     return ;
   }

  const HFRecHitCollection HithfMB = *(hfMB.product());
  cout<<" HF MB size of collection "<<HithfMB.size()<<endl;
   if(HithfMB.size()!= 1728) {
          cout<<" HF problem "<<rnnum<<" "<<HithfMB.size()<<endl;
          return;
   }



  for(HBHERecHitCollection::const_iterator hbheItr=HithbheNS.begin(); hbheItr!=HithbheNS.end(); hbheItr++)
        {
// Recalibration of energy
         float icalconst=1.;     
         DetId mydetid = hbheItr->id().rawId();
     if( theRecalib ) icalconst=myRecalib->getValues(mydetid)->getValue();
    
     HBHERecHit aHit(hbheItr->id(),hbheItr->energy()*icalconst,hbheItr->time());
     
         double energyhit = aHit.energy();
     
     DetId id = (*hbheItr).detid(); 
     HcalDetId hid=HcalDetId(id);
 
         int mysu = ((hid).rawId()>>25)&0x7;
         if( hid.ieta() > 0 ) {
     theNSFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theNSFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()]+ 1.;
         theNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()]+energyhit;
     theNSFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theNSFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()]+pow(energyhit,2);
         theNSFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theNSFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()]+pow(energyhit,4);
     
         tmpNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = energyhit;
     tmpNSFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] = pow(energyhit,2);
     
     
         } else {
     theNSFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] = theNSFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())]+ 1.;
         theNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] = theNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())]+energyhit;
     theNSFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] = theNSFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())]+pow(energyhit,2);
         theNSFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] = theNSFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())]+pow(energyhit,4);


         tmpNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = energyhit;
     tmpNSFillDetMapMin2[mysu][hid.depth()][hid.iphi()][hid.ieta()] = pow(energyhit,2);
     
     }  
     
         if(hid.depth() == 1) {
         if( hid.ieta() > 0 ) {
          hCalo1[hid.iphi()][hid.ieta()]->Fill(energyhit-noise_pl[hid.iphi()][hid.ieta()]);
          hCalo1mom2[hid.iphi()][hid.ieta()]->Fill(pow(energyhit,2));
         } else {
          hCalo2[hid.iphi()][abs(hid.ieta())]->Fill(energyhit-noise_min[hid.iphi()][abs(hid.ieta())]);
          hCalo2mom2[hid.iphi()][abs(hid.ieta())]->Fill(pow(energyhit,2));
         } // eta><0
     } // depth=1
        } // HBHE_NS


// Signal part for HB HE
     

  for(HBHERecHitCollection::const_iterator hbheItr=HithbheMB.begin(); hbheItr!=HithbheMB.end(); hbheItr++)
        {
// Recalibration of energy
         float icalconst=1.;     
         DetId mydetid = hbheItr->id().rawId();
     if( theRecalib ) icalconst=myRecalib->getValues(mydetid)->getValue();
    
     HBHERecHit aHit(hbheItr->id(),hbheItr->energy()*icalconst,hbheItr->time());
     
         double energyhit = aHit.energy();
     
     DetId id = (*hbheItr).detid(); 
     HcalDetId hid=HcalDetId(id);
 
         int mysu = ((hid).rawId()>>25)&0x7;
         if( hid.ieta() > 0 ) {
     theMBFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theMBFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()]+ 1.;
         theMBFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theMBFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()]+energyhit;
     theMBFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theMBFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()]+pow(energyhit,2);
         theMBFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theMBFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()]+pow(energyhit,4);
     float mydiff = energyhit - tmpNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()];
     
     
         theDFFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theDFFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()]+mydiff;
     theDFFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theDFFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()]+pow(mydiff,2);
         } else {
     theMBFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] = theMBFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())]+ 1.;
         theMBFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] = theMBFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())]+energyhit;
     theMBFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] = theMBFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())]+pow(energyhit,2);
         theMBFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] = theMBFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())]+pow(energyhit,4);


     float mydiff = energyhit - tmpNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()];
         theDFFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theDFFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()]+mydiff;
     theDFFillDetMapMin2[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theDFFillDetMapMin2[mysu][hid.depth()][hid.iphi()][hid.ieta()]+pow(mydiff,2);
     }  

     
         if(hid.depth() == 1) {
         if( hid.ieta() > 0 ) {
          hCalo1[hid.iphi()][hid.ieta()]->Fill(energyhit);
          hCalo1mom2[hid.iphi()][hid.ieta()]->Fill(pow(energyhit,2));
         } else {
          hCalo2[hid.iphi()][abs(hid.ieta())]->Fill(energyhit);
          hCalo2mom2[hid.iphi()][abs(hid.ieta())]->Fill(pow(energyhit,2));
         } // eta><0
     } // depth=1
        } // HBHE_MB
    
// HF
 
  for(HFRecHitCollection::const_iterator hbheItr=HithfNS.begin(); hbheItr!=HithfNS.end(); hbheItr++)
        {
// Recalibration of energy
         float icalconst=1.;     
         DetId mydetid = hbheItr->id().rawId();
     if( theRecalib ) icalconst=myRecalib->getValues(mydetid)->getValue();
    
     HFRecHit aHit(hbheItr->id(),hbheItr->energy()*icalconst,hbheItr->time());
     
         double energyhit = aHit.energy();
     
     DetId id = (*hbheItr).detid(); 
     HcalDetId hid=HcalDetId(id);
 
         int mysu = ((hid).rawId()>>25)&0x7;
         if( hid.ieta() > 0 ) {
     theNSFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theNSFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()]+ 1.;
         theNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()]+energyhit;
     theNSFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theNSFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()]+pow(energyhit,2);
         theNSFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theNSFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()]+pow(energyhit,4);
     
         tmpNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = energyhit;
     tmpNSFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] = pow(energyhit,2);
     
     
         } else {
     theNSFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] = theNSFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())]+ 1.;
         theNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] = theNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())]+energyhit;
     theNSFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] = theNSFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())]+pow(energyhit,2);
         theNSFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] = theNSFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())]+pow(energyhit,4);

 
         tmpNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = energyhit;
     tmpNSFillDetMapMin2[mysu][hid.depth()][hid.iphi()][hid.ieta()] = pow(energyhit,2);
     
     }  
     
         if(hid.depth() == 1) {
         if( hid.ieta() > 0 ) {
          hCalo1[hid.iphi()][hid.ieta()]->Fill(energyhit-noise_pl[hid.iphi()][hid.ieta()]);
          hCalo1mom2[hid.iphi()][hid.ieta()]->Fill(pow(energyhit,2));
         } else {
          hCalo2[hid.iphi()][abs(hid.ieta())]->Fill(energyhit-noise_min[hid.iphi()][abs(hid.ieta())]);
          hCalo2mom2[hid.iphi()][abs(hid.ieta())]->Fill(pow(energyhit,2));
         } // eta><0
     } // depth=1
        } // HBHE_NS


// Signal part for HB HE

  for(HFRecHitCollection::const_iterator hbheItr=HithfMB.begin(); hbheItr!=HithfMB.end(); hbheItr++)
        {
// Recalibration of energy
         float icalconst=1.;     
         DetId mydetid = hbheItr->id().rawId();
     if( theRecalib ) icalconst=myRecalib->getValues(mydetid)->getValue();
    
     HFRecHit aHit(hbheItr->id(),hbheItr->energy()*icalconst,hbheItr->time());
     
         double energyhit = aHit.energy();
     
     DetId id = (*hbheItr).detid(); 
     HcalDetId hid=HcalDetId(id);
 
         int mysu = ((hid).rawId()>>25)&0x7;
         if( hid.ieta() > 0 ) {
     theMBFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theMBFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()]+ 1.;
         theMBFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theMBFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()]+energyhit;
     theMBFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theMBFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()]+pow(energyhit,2);
         theMBFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theMBFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()]+pow(energyhit,4);


         theDFFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theDFFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()]+energyhit-tmpNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()];
     theDFFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
     theDFFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()]+pow((energyhit-tmpNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()]),2);
         } else {
     theMBFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] = theMBFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())]+ 1.;
         theMBFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] = theMBFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())]+energyhit;
     theMBFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] = theMBFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())]+pow(energyhit,2);
         theMBFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] = theMBFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())]+pow(energyhit,4);
 
         theDFFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = theDFFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()]+energyhit-tmpNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()];
     theDFFillDetMapMin2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
     theDFFillDetMapMin2[mysu][hid.depth()][hid.iphi()][hid.ieta()]+pow((energyhit-tmpNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()]),2);
     }  

     
         if(hid.depth() == 1) {
         if( hid.ieta() > 0 ) {
          hCalo1[hid.iphi()][hid.ieta()]->Fill(energyhit);
          hCalo1mom2[hid.iphi()][hid.ieta()]->Fill(pow(energyhit,2));
         } else {
          hCalo2[hid.iphi()][abs(hid.ieta())]->Fill(energyhit);
          hCalo2mom2[hid.iphi()][abs(hid.ieta())]->Fill(pow(energyhit,2));
         } // eta><0
     } // depth=1
        } // HF_MB

   std::cout<<" Event is finished "<<std::endl;
}

void cms::Analyzer_minbias::beginJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 81 of file Analyzer_minbias.cc.

References gather_cfg::cout, eta(), i, j, gen::k, prof2calltree::l, nevent, phi, and mathSSE::return().

{
   
   hOutputFile   = new TFile( fOutputFileName.c_str(), "RECREATE" ) ;
   
   myTree = new TTree("RecJet","RecJet Tree");
   myTree->Branch("mydet",  &mydet, "mydet/I");
   myTree->Branch("mysubd",  &mysubd, "mysubd/I");
   myTree->Branch("depth",  &depth, "depth/I");
   myTree->Branch("ieta",  &ieta, "ieta/I");
   myTree->Branch("iphi",  &iphi, "iphi/I");
   myTree->Branch("eta",  &eta, "eta/F");
   myTree->Branch("phi",  &phi, "phi/F");
   
   myTree->Branch("mom0_MB",  &mom0_MB, "mom0_MB/F");
   myTree->Branch("mom1_MB",  &mom1_MB, "mom1_MB/F");
   myTree->Branch("mom2_MB",  &mom2_MB, "mom2_MB/F");
   myTree->Branch("mom4_MB",  &mom4_MB, "mom4_MB/F");
   
   myTree->Branch("mom0_Noise",  &mom0_Noise, "mom0_Noise/F");
   myTree->Branch("mom1_Noise",  &mom1_Noise, "mom1_Noise/F");
   myTree->Branch("mom2_Noise",  &mom2_Noise, "mom2_Noise/F");
   myTree->Branch("mom4_Noise",  &mom4_Noise, "mom4_Noise/F");
   
   myTree->Branch("mom0_Diff",  &mom0_Diff, "mom0_Diff/F");
   myTree->Branch("mom1_Diff",  &mom1_Diff, "mom1_Diff/F");
   myTree->Branch("mom2_Diff",  &mom2_Diff, "mom2_Diff/F");

   myTree->Branch("occup",  &occup, "occup/F");
   
   std::cout<<" Before ordering Histos "<<std::endl;
  
   char str0[15];
   char str1[15];

   char str10[15];
   char str11[15];

   int k=0;
   nevent = 0;
// Size of collections

   hHBHEsize_vs_run = new TH2F("hHBHEsize_vs_run","hHBHEsize_vs_run",500,111500.,112000.,6101,-100.5,6000.5);
   hHFsize_vs_run = new TH2F("hHFsize_vs_run","hHFsize_vs_run",500,111500.,112000.,6101,-100.5,6000.5); 

   for(int i=1;i<73;i++){
    for(int j=1;j<43;j++){

       meannoise_pl[i][j] = 0.;
       meannoise_min[i][j] = 0.;

//     for(int l=1;l<5;l++){
        k = i*1000+j;
        sprintf(str0,"mpl%d",k);
        sprintf(str1,"mmin%d",k);

        sprintf(str10,"vpl%d",k);
        sprintf(str11,"vmin%d",k);
//      cout<<" "<<i<<" "<<j<<endl;
    if( j < 30 )
    {
// first order moment
    hCalo1[i][j] = new TH1F(str0, "h0", 320, -10., 10.);
    hCalo2[i][j] = new TH1F(str1, "h1", 320, -10., 10.);

// second order moment
    hCalo1mom2[i][j] = new TH1F(str10, "h10", 320, 0., 20.);
    hCalo2mom2[i][j] = new TH1F(str11, "h11", 320, 0., 20.);
    }
      else
      {
// HF
// first order moment
//   cout<<" "<<i<<" "<<j<<" "<<k<<endl;
   if(j < 40)
   {
    hCalo1[i][j] = new TH1F(str0, "h0", 320, -10., 10.);
    hCalo2[i][j] = new TH1F(str1, "h1", 320, -10., 10.);
//
// second order moment
    hCalo1mom2[i][j] = new TH1F(str10, "h10", 320, 0., 40.);
    hCalo2mom2[i][j] = new TH1F(str11, "h11", 320, 0., 40.);
   }
     else
     {
    hCalo1[i][j] = new TH1F(str0,"h0" , 320, -10., 10.);
    hCalo2[i][j] = new TH1F(str1, "h1", 320, -10., 10.);

// second order moment
    hCalo1mom2[i][j] = new TH1F(str10, "h10", 320, 0., 120.);
    hCalo2mom2[i][j] = new TH1F(str11, "h11", 320, 0., 120.);

     }
    } // HE/HF boundary
//     } // l
    } // j
   } // i


   std::cout<<" After ordering Histos "<<std::endl;

   std::string ccc = "noise_0.dat";

   myout_hcal = new ofstream(ccc.c_str());
   if(!myout_hcal) cout << " Output file not open!!! "<<endl;

//
   for (int i=0; i<5;i++)
   {
    for (int j=0; j<5;j++)
    {
     for (int k=0; k<73;k++)
     {
       for (int l=0; l<43;l++)
       {
        theMBFillDetMapPl0[i][j][k][l] = 0.;
        theMBFillDetMapPl1[i][j][k][l] = 0.;
        theMBFillDetMapPl2[i][j][k][l] = 0.;
        theMBFillDetMapPl4[i][j][k][l] = 0.;

        theMBFillDetMapMin0[i][j][k][l] = 0.;
        theMBFillDetMapMin1[i][j][k][l] = 0.;
        theMBFillDetMapMin2[i][j][k][l] = 0.;
        theMBFillDetMapMin4[i][j][k][l] = 0.;


        theNSFillDetMapPl0[i][j][k][l] = 0.;
        theNSFillDetMapPl1[i][j][k][l] = 0.;
        theNSFillDetMapPl2[i][j][k][l] = 0.;
        theNSFillDetMapPl4[i][j][k][l] = 0.;

        theNSFillDetMapMin0[i][j][k][l] = 0.;
        theNSFillDetMapMin1[i][j][k][l] = 0.;
        theNSFillDetMapMin2[i][j][k][l] = 0.;
        theNSFillDetMapMin4[i][j][k][l] = 0.;

        theDFFillDetMapPl0[i][j][k][l] = 0.;
        theDFFillDetMapPl1[i][j][k][l] = 0.;
        theDFFillDetMapPl2[i][j][k][l] = 0.;
        theDFFillDetMapMin0[i][j][k][l] = 0.;
        theDFFillDetMapMin1[i][j][k][l] = 0.;
        theDFFillDetMapMin2[i][j][k][l] = 0.;
       }
     }  
    }
   }    
    
   return ;
}

void cms::Analyzer_minbias::beginRun ( const edm::Run &  r, const edm::EventSetup &  iSetup  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 72 of file Analyzer_minbias.cc.

{
  nevent_run = 0;
}

void cms::Analyzer_minbias::endJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 233 of file Analyzer_minbias.cc.

References gather_cfg::cout, i, j, gen::k, prof2calltree::l, and mathSSE::return().

{
   int ii=0;
      
   for (int i=1; i<5;i++)
   {
    for (int j=1; j<5;j++)
    {
     for (int k=1; k<73;k++)
     {
       for (int l=1; l<43;l++)
       {
      if(theMBFillDetMapPl0[i][j][k][l] > 0)
      { 
            mom0_MB = theMBFillDetMapPl0[i][j][k][l];
            mom1_MB = theMBFillDetMapPl1[i][j][k][l];
            mom2_MB = theMBFillDetMapPl2[i][j][k][l];
            mom4_MB = theMBFillDetMapPl4[i][j][k][l];
            mom0_Noise = theNSFillDetMapPl0[i][j][k][l];
            mom1_Noise = theNSFillDetMapPl1[i][j][k][l];
            mom2_Noise = theNSFillDetMapPl2[i][j][k][l];
            mom4_Noise = theNSFillDetMapPl4[i][j][k][l];
            mom0_Diff = theDFFillDetMapPl0[i][j][k][l];
            mom1_Diff = theDFFillDetMapPl1[i][j][k][l];
            mom2_Diff = theDFFillDetMapPl2[i][j][k][l];
        
            mysubd = i;
            depth = j;
            ieta = l;
            iphi = k;
            cout<<" Result Plus= "<<mysubd<<" "<<ieta<<" "<<iphi<<" mom0  "<<mom0_MB<<" mom1 "<<mom1_MB<<" mom2 "<<mom2_MB<<endl;
            myTree->Fill();
            ii++;
       } // Pl > 0
    
         
      if(theMBFillDetMapMin0[i][j][k][l] > 0)
      { 
            mom0_MB = theMBFillDetMapMin0[i][j][k][l];
            mom1_MB = theMBFillDetMapMin1[i][j][k][l];
            mom2_MB = theMBFillDetMapMin2[i][j][k][l];
            mom4_MB = theMBFillDetMapMin4[i][j][k][l];
            mom0_Noise = theNSFillDetMapMin0[i][j][k][l];
            mom1_Noise = theNSFillDetMapMin1[i][j][k][l];
            mom2_Noise = theNSFillDetMapMin2[i][j][k][l];
            mom4_Noise = theNSFillDetMapMin4[i][j][k][l];
            mom0_Diff = theDFFillDetMapMin0[i][j][k][l];
            mom1_Diff = theDFFillDetMapMin1[i][j][k][l];
            mom2_Diff = theDFFillDetMapMin2[i][j][k][l];
        
        
            mysubd = i;
            depth = j;
            ieta = -1*l;
            iphi = k;
            cout<<" Result Minus= "<<mysubd<<" "<<ieta<<" "<<iphi<<" mom0  "<<mom0_MB<<" mom1 "<<mom1_MB<<" mom2 "<<mom2_MB<<endl;
            myTree->Fill();
            ii++;
        
      } // Min>0  
      } // ieta
     } // iphi  
    } // depth
   } //subd    
   
   
   
   cout<<" Number of cells "<<ii<<endl; 
      
   hOutputFile->Write() ;   
   hOutputFile->cd();
   hHBHEsize_vs_run->Write() ;
   hHFsize_vs_run->Write() ;
   for(int i=1;i<73;i++){
    for(int j=1;j<43;j++){
    hCalo1[i][j]->Write();
    hCalo2[i][j]->Write();
    hCalo1mom2[i][j]->Write();
    hCalo2mom2[i][j]->Write();
   }
  }
  
   
   myTree->Write();
   hOutputFile->Close() ;
   
   return ;
}

void cms::Analyzer_minbias::endRun ( const edm::Run &  r, const edm::EventSetup &  iSetup  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 76 of file Analyzer_minbias.cc.

References gather_cfg::cout, and edm::RunBase::run().

{
 std::cout<<" Runnumber "<<r.run()<<" Nevents  "<<nevent_run<<std::endl;
}

Member Data Documentation

int cms::Analyzer_minbias::depth

private

Definition at line 80 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::eta

private

Definition at line 81 of file Analyzer_minbias.h.

std::string cms::Analyzer_minbias::fOutputFileName

private

Definition at line 64 of file Analyzer_minbias.h.

edm::InputTag cms::Analyzer_minbias::hbherecoMB

private

Definition at line 121 of file Analyzer_minbias.h.

edm::InputTag cms::Analyzer_minbias::hbherecoNoise

private

Definition at line 125 of file Analyzer_minbias.h.

std::string cms::Analyzer_minbias::hcalfile_

private

Definition at line 65 of file Analyzer_minbias.h.

TH1F* cms::Analyzer_minbias::hCalo1[73][43]

private

Definition at line 72 of file Analyzer_minbias.h.

TH1F* cms::Analyzer_minbias::hCalo1mom2[73][43]

private

Definition at line 74 of file Analyzer_minbias.h.

TH1F* cms::Analyzer_minbias::hCalo2[73][43]

private

Definition at line 73 of file Analyzer_minbias.h.

TH1F* cms::Analyzer_minbias::hCalo2mom2[73][43]

private

Definition at line 75 of file Analyzer_minbias.h.

edm::InputTag cms::Analyzer_minbias::hfrecoMB

private

Definition at line 123 of file Analyzer_minbias.h.

edm::InputTag cms::Analyzer_minbias::hfrecoNoise

private

Definition at line 127 of file Analyzer_minbias.h.

TH2F* cms::Analyzer_minbias::hHBHEsize_vs_run

private

Definition at line 76 of file Analyzer_minbias.h.

TH2F* cms::Analyzer_minbias::hHFsize_vs_run

private

Definition at line 77 of file Analyzer_minbias.h.

edm::InputTag cms::Analyzer_minbias::horecoMB

private

Definition at line 122 of file Analyzer_minbias.h.

edm::InputTag cms::Analyzer_minbias::horecoNoise

private

Definition at line 126 of file Analyzer_minbias.h.

TFile* cms::Analyzer_minbias::hOutputFile

private

Definition at line 70 of file Analyzer_minbias.h.

int cms::Analyzer_minbias::ieta

private

Definition at line 80 of file Analyzer_minbias.h.

int cms::Analyzer_minbias::iphi

private

Definition at line 80 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::meannoise_min[73][43]

private

Definition at line 88 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::meannoise_pl[73][43]

private

Definition at line 88 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom0_Diff

private

Definition at line 84 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom0_MB

private

Definition at line 82 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom0_Noise

private

Definition at line 83 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom1_Diff

private

Definition at line 84 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom1_MB

private

Definition at line 82 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom1_Noise

private

Definition at line 83 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom2_Diff

private

Definition at line 84 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom2_MB

private

Definition at line 82 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom2_Noise

private

Definition at line 83 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom3_Diff

private

Definition at line 84 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom3_MB

private

Definition at line 82 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom3_Noise

private

Definition at line 83 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom4_Diff

private

Definition at line 84 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom4_MB

private

Definition at line 82 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::mom4_Noise

private

Definition at line 83 of file Analyzer_minbias.h.

int cms::Analyzer_minbias::mydet

private

Definition at line 80 of file Analyzer_minbias.h.

std::ofstream* cms::Analyzer_minbias::myout_hcal

private

Definition at line 66 of file Analyzer_minbias.h.

int cms::Analyzer_minbias::mysubd

private

Definition at line 80 of file Analyzer_minbias.h.

TTree* cms::Analyzer_minbias::myTree

private

Definition at line 71 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::nevent

private

Definition at line 93 of file Analyzer_minbias.h.

int cms::Analyzer_minbias::nevent_run

private

Definition at line 79 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::noise_min[73][43]

private

Definition at line 89 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::noise_pl[73][43]

private

Definition at line 89 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::occup

private

Definition at line 82 of file Analyzer_minbias.h.

float cms::Analyzer_minbias::phi

private

Definition at line 81 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theDFFillDetMapMin0[5][5][73][43]

private

Definition at line 117 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theDFFillDetMapMin1[5][5][73][43]

private

Definition at line 118 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theDFFillDetMapMin2[5][5][73][43]

private

Definition at line 119 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theDFFillDetMapPl0[5][5][73][43]

private

Definition at line 114 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theDFFillDetMapPl1[5][5][73][43]

private

Definition at line 115 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theDFFillDetMapPl2[5][5][73][43]

private

Definition at line 116 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theMBFillDetMapMin0[5][5][73][43]

private

Definition at line 99 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theMBFillDetMapMin1[5][5][73][43]

private

Definition at line 100 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theMBFillDetMapMin2[5][5][73][43]

private

Definition at line 101 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theMBFillDetMapMin4[5][5][73][43]

private

Definition at line 102 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theMBFillDetMapPl0[5][5][73][43]

private

Definition at line 94 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theMBFillDetMapPl1[5][5][73][43]

private

Definition at line 95 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theMBFillDetMapPl2[5][5][73][43]

private

Definition at line 96 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theMBFillDetMapPl4[5][5][73][43]

private

Definition at line 97 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theNSFillDetMapMin0[5][5][73][43]

private

Definition at line 109 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theNSFillDetMapMin1[5][5][73][43]

private

Definition at line 110 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theNSFillDetMapMin2[5][5][73][43]

private

Definition at line 111 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theNSFillDetMapMin4[5][5][73][43]

private

Definition at line 112 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theNSFillDetMapPl0[5][5][73][43]

private

Definition at line 104 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theNSFillDetMapPl1[5][5][73][43]

private

Definition at line 105 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theNSFillDetMapPl2[5][5][73][43]

private

Definition at line 106 of file Analyzer_minbias.h.

double cms::Analyzer_minbias::theNSFillDetMapPl4[5][5][73][43]

private

Definition at line 107 of file Analyzer_minbias.h.

bool cms::Analyzer_minbias::theRecalib

private

Definition at line 130 of file Analyzer_minbias.h.