DTTTrigCorrectionFirst Class Reference

#include <DTTTrigCorrectionFirst.h>

Inheritance diagram for DTTTrigCorrectionFirst:

Public Member Functions
virtual void	analyze (const edm::Event &event, const edm::EventSetup &setup)
virtual void	beginJob ()
virtual void	beginRun (const edm::Run &run, const edm::EventSetup &setup)
	DTTTrigCorrectionFirst (const edm::ParameterSet &pset)
	Constructor. More...
virtual void	endJob ()
virtual	~DTTTrigCorrectionFirst ()
	Destructor. More...
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Private Attributes
std::string	dbLabel
bool	debug
edm::ESHandle< DTGeometry >	muonGeom
double	rmsLimit
const DTTtrig *	tTrigMap
double	ttrigMax
double	ttrigMin

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)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 22 of file DTTTrigCorrectionFirst.h.

Constructor & Destructor Documentation

DTTTrigCorrectionFirst::DTTTrigCorrectionFirst

(

const edm::ParameterSet &

pset

)

Constructor.

Definition at line 30 of file DTTTrigCorrectionFirst.cc.

References dtTTrigAnalyzer_cfg::dbLabel, debug, and edm::ParameterSet::getUntrackedParameter().

                                                                       {
 
  debug = pset.getUntrackedParameter<bool>("debug",false);
  ttrigMin = pset.getUntrackedParameter<double>("ttrigMin",0);
  ttrigMax = pset.getUntrackedParameter<double>("ttrigMax",5000);
  rmsLimit = pset.getUntrackedParameter<double>("rmsLimit",5.);

  dbLabel  = pset.getUntrackedParameter<string>("dbLabel", "");
}

DTTTrigCorrectionFirst::~DTTTrigCorrectionFirst ( )

virtual

Destructor.

Definition at line 40 of file DTTTrigCorrectionFirst.cc.

{}

Member Function Documentation

virtual void DTTTrigCorrectionFirst::analyze ( const edm::Event &  event, const edm::EventSetup &  setup  )

inlinevirtual

Implements edm::EDAnalyzer.

Definition at line 34 of file DTTTrigCorrectionFirst.h.

{}

virtual void DTTTrigCorrectionFirst::beginJob ( void  )

inlinevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 32 of file DTTTrigCorrectionFirst.h.

{}

void DTTTrigCorrectionFirst::beginRun ( const edm::Run &  run, const edm::EventSetup &  setup  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 42 of file DTTTrigCorrectionFirst.cc.

References gather_cfg::cout, dtTTrigAnalyzer_cfg::dbLabel, and edm::EventSetup::get().

                                                                                     {
  ESHandle<DTTtrig> tTrig;
  setup.get<DTTtrigRcd>().get(dbLabel,tTrig);
  tTrigMap = &*tTrig;
  cout << "[DTTTrigCorrection]: TTrig version: " << tTrig->version() << endl;

  setup.get<MuonGeometryRecord>().get(muonGeom);
}

void DTTTrigCorrectionFirst::endJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 51 of file DTTTrigCorrectionFirst.cc.

References abs, PDRates::average, gather_cfg::cout, debug, DTTimeUnits::ns, record, plotscripts::rms(), DTTtrig::set(), and mathSSE::sqrt().

                                    {
  // Create the object to be written to DB
  DTTtrig* tTrigNewMap = new DTTtrig();  
  //Get the superlayers list
  vector<DTSuperLayer*> dtSupLylist = muonGeom->superLayers();

  //Loop on all superlayers to compute the mean 
  double average = 0.;
  double average2 = 0.;
  double rms = 0.;
  double averageSigma = 0.;
  double average2Sigma = 0.;
  double rmsSigma = 0.;
  double counter = 0.;
  double averagekfactor = 0;
  float kfactor = 0;

  for (vector<DTSuperLayer*>::const_iterator sl = dtSupLylist.begin();
       sl != dtSupLylist.end(); sl++) {
    float ttrigMean = 0;
    float ttrigSigma = 0;   
    tTrigMap->get((*sl)->id(),
          ttrigMean,
          ttrigSigma,
                  kfactor,
          DTTimeUnits::ns);
    if( ttrigMean < ttrigMax && ttrigMean > ttrigMin ) {
      average += ttrigMean;
      averageSigma += ttrigSigma;
      if(ttrigMean > 0)
    counter +=  1.;
    }
  }//End of loop on superlayers 

  average =  average/ counter;
  averageSigma = averageSigma/counter;
  averagekfactor = kfactor;

  //  cout << " average counter "<< average << " "<< counter <<endl;

  for (vector<DTSuperLayer*>::const_iterator sl = dtSupLylist.begin();
       sl != dtSupLylist.end(); sl++) {
    float ttrigMean = 0;
    float ttrigSigma = 0;
    float kfactor = 0;
    tTrigMap->get((*sl)->id(),
          ttrigMean,
          ttrigSigma,
                  kfactor,
          DTTimeUnits::ns);
    if( ttrigMean < ttrigMax && ttrigMean > ttrigMin ) {
      average2 += (ttrigMean-average)*(ttrigMean-average);
      average2Sigma += (ttrigSigma-averageSigma)*(ttrigSigma-averageSigma);
    }
  }//End of loop on superlayers 

     rms = average2 /(counter-1);
     rmsSigma = average2Sigma /(counter-1);
     rms = sqrt(rms);
     rmsSigma = sqrt(rmsSigma);
  cout << "average averageSigma counter rms "<< average <<" " << averageSigma << " " << counter << " " << rms  <<endl;


  for (vector<DTSuperLayer*>::const_iterator sl = dtSupLylist.begin();
       sl != dtSupLylist.end(); sl++) {
    //Compute new ttrig
    double newTTrigMean =  0;
    double newTTrigSigma =  0;
    double newkfactor =0;
    float tempttrigMean = 0;
    float tempttrigSigma = 0;
    float tempkfactor = 0;
    float ttrigMean = 0;
    float ttrigSigma = 0;
    float kfactor = 0;
    tTrigMap->get((*sl)->id(),
          ttrigMean,
          ttrigSigma,
                  kfactor,
          DTTimeUnits::ns);
    int chamber = (*sl)->id().chamberId().station();

    //check if ttrigMean is similar to the mean
    if (abs(ttrigMean - average) < rmsLimit*rms ){
      newTTrigMean = ttrigMean;
      newTTrigSigma = ttrigSigma;
      newkfactor = kfactor;
    } else {
      // do not consider if ttrig == 0
      if(ttrigMean > 0) {
    //cout << "ttrig chamber " << ttrigMean <<" "<<chamber<<endl;
    if(((*sl)->id().superlayer()) == 1){
      //cout << " superlayer " << ((*sl)->id().superlayer()) <<endl; 
      DTSuperLayerId slId((*sl)->id().chamberId(),3);
      tTrigMap->get(slId,
            tempttrigMean,
            tempttrigSigma,
                        tempkfactor,
            DTTimeUnits::ns);
      if (abs(tempttrigMean - average) < rmsLimit*rms) {
        newTTrigMean = tempttrigMean;
        newTTrigSigma = tempttrigSigma;
            newkfactor = tempkfactor;
        cout <<"Chamber "<< chamber << " sl " << ((*sl)->id().superlayer()) << "has ttrig "<< ttrigMean  
         << " -> takes value of sl 3 " <<  newTTrigMean <<endl;
      } else if(chamber == 4){
        cout << "No correction possible within same chamber (sl1) "  <<endl; 
        newTTrigMean = average;
        newTTrigSigma = averageSigma;
            newkfactor = averagekfactor;
        cout<<"####### Bad SL: " << (*sl)->id() << " from "<< ttrigMean <<" to "<< newTTrigMean <<endl;
      } else if(chamber != 4){    
        DTSuperLayerId slId((*sl)->id().chamberId(),2);
        tTrigMap->get(slId,
              tempttrigMean,
              tempttrigSigma,
                          tempkfactor,
              DTTimeUnits::ns);
        if (abs(tempttrigMean - average) < rmsLimit*rms) {
          newTTrigMean = tempttrigMean;
          newTTrigSigma = tempttrigSigma;
              newkfactor = tempkfactor;
          cout <<"Chamber "<< chamber << " sl " << ((*sl)->id().superlayer()) << "has ttrig "<< ttrigMean  
           <<" -> takes value of sl 2 " <<  newTTrigMean <<endl;
        } else {
          cout << "No correction possible within same chamber (sl1) "  <<endl;
          newTTrigMean = average;
          newTTrigSigma = averageSigma;
              newkfactor = averagekfactor;
          cout<<"####### Bad SL: " << (*sl)->id() << " from "<< ttrigMean <<" to "<< newTTrigMean <<endl;
        }
      }
    } else if (((*sl)->id().superlayer()) == 2) {
      //cout << " superlayer " << ((*sl)->id().superlayer()) <<endl; 
      DTSuperLayerId slId((*sl)->id().chamberId(),1);
      tTrigMap->get(slId,
            tempttrigMean,
            tempttrigSigma,
                        tempkfactor,
            DTTimeUnits::ns);
      if (abs(tempttrigMean - average) < rmsLimit*rms) {
        newTTrigMean = tempttrigMean;
        newTTrigSigma = tempttrigSigma;
            newkfactor = tempkfactor;
        cout <<"Chamber "<< chamber << " sl " << ((*sl)->id().superlayer()) << "has ttrig "<< ttrigMean  
         << " -> takes value of sl 1 " <<  newTTrigMean <<endl;
      } else {
        DTSuperLayerId slId((*sl)->id().chamberId(),3);
        tTrigMap->get(slId,
              tempttrigMean,
              tempttrigSigma,
                          tempkfactor,
              DTTimeUnits::ns);
        if (abs(tempttrigMean - average) < rmsLimit*rms) {
          newTTrigMean = tempttrigMean;
          newTTrigSigma = tempttrigSigma;
              newkfactor = tempkfactor;
          cout <<"Chamber "<< chamber << " sl " << ((*sl)->id().superlayer()) << "has ttrig "<< ttrigMean  
           << "-> takes value of sl 3 " <<  newTTrigMean <<endl;
          } else {
          cout << "No correction possible within same chamber (sl2)  "  <<endl;
          newTTrigMean = average;
          newTTrigSigma = averageSigma;
              newkfactor = averagekfactor;
          cout<<"####### Bad SL: " << (*sl)->id() << " from "<< ttrigMean <<" to "<< newTTrigMean <<endl;
        }
      } 
      
      
    } else if (((*sl)->id().superlayer()) == 3) {
      //cout << " superlayer " << ((*sl)->id().superlayer()) <<endl; 
      DTSuperLayerId slId((*sl)->id().chamberId(),1);
      tTrigMap->get(slId,
            tempttrigMean,
            tempttrigSigma,
                        tempkfactor,
            DTTimeUnits::ns);
      if (abs(tempttrigMean - average) < rmsLimit*rms) {
        newTTrigMean = tempttrigMean;
        newTTrigSigma = tempttrigSigma;
            newkfactor = tempkfactor;
        cout <<"Chamber "<< chamber << " sl " << ((*sl)->id().superlayer()) << "has ttrig "<< ttrigMean  
         << " -> takes value of sl 1 " <<  newTTrigMean <<endl;
      } else if(chamber == 4) {
        cout << "No correction possible within same chamber (sl3)"  <<endl;
        newTTrigMean = average;
        newTTrigSigma = averageSigma;
            newkfactor = averagekfactor;
        cout<<"####### Bad SL: " << (*sl)->id() << " from "<< ttrigMean <<" to "<< newTTrigMean <<endl;
      } else if(chamber != 4){
        DTSuperLayerId slId((*sl)->id().chamberId(),2);
        tTrigMap->get(slId,
              tempttrigMean,
              tempttrigSigma,
                          tempkfactor,
              DTTimeUnits::ns);
        if (abs(tempttrigMean - average) < rmsLimit*rms) {
          newTTrigMean = tempttrigMean;
          newTTrigSigma = tempttrigSigma;
              newkfactor = tempkfactor;
          cout <<"Chamber "<< chamber << " sl " << ((*sl)->id().superlayer()) << "has ttrig "<< ttrigMean  
           << " -> takes value of sl 2 " <<  newTTrigMean <<endl;
        } else {
          cout << "No correction possible within same chamber (sl3) "  <<endl;
          newTTrigMean = average;
          newTTrigSigma = averageSigma;
              newkfactor = averagekfactor;
          cout<<"####### Bad SL: " << (*sl)->id() << " from "<< ttrigMean <<" to "<< newTTrigMean <<endl;
        }
      } 
      
    }
      } else {
    //  cout << "SL not present " << ttrigMean << " " <<((*sl)->id().superlayer()) <<endl;
    newTTrigMean = average;
    newTTrigSigma = averageSigma;
        newkfactor = averagekfactor;
    cout<<"####### NotPresent SL: " << (*sl)->id() << " from "<< ttrigMean <<" to "<< newTTrigMean <<endl;
      }
    }

    //Store new ttrig in the new map
    tTrigNewMap->set((*sl)->id(), newTTrigMean, newTTrigSigma, newkfactor, DTTimeUnits::ns);
    if(debug){
      cout<<"New tTrig: " << (*sl)->id()
          << " from "<< ttrigMean <<" to "<< newTTrigMean <<endl;
      cout<<"New tTrigSigma: " << (*sl)->id()
          << " from "<<ttrigSigma  <<" to "<< newTTrigSigma <<endl;

    }
  }//End of loop on superlayers 

  
  //Write object to DB
  cout << "[DTTTrigCorrection]: Writing ttrig object to DB!" << endl;
  string record = "DTTtrigRcd";
   DTCalibDBUtils::writeToDB<DTTtrig>(record, tTrigNewMap);
} 

Member Data Documentation

std::string DTTTrigCorrectionFirst::dbLabel

private

Definition at line 44 of file DTTTrigCorrectionFirst.h.

bool DTTTrigCorrectionFirst::debug

private

Definition at line 46 of file DTTTrigCorrectionFirst.h.

edm::ESHandle<DTGeometry> DTTTrigCorrectionFirst::muonGeom

private

Definition at line 42 of file DTTTrigCorrectionFirst.h.

double DTTTrigCorrectionFirst::rmsLimit

private

Definition at line 47 of file DTTTrigCorrectionFirst.h.

const DTTtrig* DTTTrigCorrectionFirst::tTrigMap

private

Definition at line 41 of file DTTTrigCorrectionFirst.h.

double DTTTrigCorrectionFirst::ttrigMax

private

Definition at line 47 of file DTTTrigCorrectionFirst.h.

double DTTTrigCorrectionFirst::ttrigMin

private

Definition at line 47 of file DTTTrigCorrectionFirst.h.